11.12.2011, Mythbusters
Пора разрушить парочку мифов!
Содержание |
Разрушение мифов
Давно хотел поднять эту тему. В нашем деле есть множество положений, которые воспринимаются, как данность. Они либо следуют из общих соображений, либо кажутся очевидными, либо являются результатом недоведённых до конца исследований. Либо это просто заблуждения. Короче говоря, природа их различна.
Некоторые из этих положений вызывают сомнения, но пока не проведено соответствующее исследование, сомнения остаются смутными. Может быть, они и неверны. Но точку поставить надо.
Заниматься исследованиями проблем, которые актуальными в текущий момент не являются, никто не будет. А когда появляется время, о них никто не помнит. Поэтому такие положения годами остаются "повисшими и воздухе".
Предлагаю вести список таких "подозрительных" положений и, при случае, разбираться в них. Если о проблеме никто не помнит, никто её не решит. Беда в том, что о проблеме всегда вспоминают, когда полно других дел.
Разрешение данных проблем можно доверить студентам в качестве работ - польза будет для всех. И такой пример уже есть - Иван Липа сейчас исследует двухэтапную процедуру поиска, подробности ниже.
Прошу присоединяться, "мифов для разрушения" вокруг существует большое количество.
Список "для затравки"
Дискриминаторы: I*Q или sign(I)*Q
Общепринято, что при большом отношении сигнал шум можно использовать дискриминатор sign(I)*Q, а при малом - I*Q.
Корни этого положения растут из оптимального синтеза. Оптимальный синтез алгоритма слежения за фазой даёт дискриминатор th(I)*Q, где th(X) - гипертангенс. Форма функции гипертангенса такова, что если аргумент очень большой (X >> 1), то он похож на знаковую функцию: . Если аргумент очень мал (много меньше единицы), то функция гипертангенса имеет единичную производную, тогда .
Так как систематическая составляющая на выходе коррелятора пропорциональна отношению сигнал/шум, и возникло положение, что при малом отношении следует использовать дискриминатор I*Q, а при большом отношении sign(I)*Q.
Вопрос в том, насколько мало или велико отношение сигнал/шум, о котором идёт речь. Систематическая составляющая на выходе коррелятора упрощённо имеет вид . Если сопоставить данную величину с единицей, то получим: . При времени накопления 1 мс q будет равно 26 дБГц, что ещё является реальной областью работы ФАП. Однако при столь низком отношении обычно используется более длительное время накопления сигнала в корреляторе. Если взять 20 мс, то величина q будет равна 14 дБГц, что уже совсем мало для ФАП. ФАП нет смысла использовать при таком отношении, т.к. выделить данные всё равно будет невозможно.
Таким образом, сомнительной является ситуация "", в которой "оптимальным" является дискриминатор I*Q. Стоит ли вообще когда-либо использовать данный дискримитатор?
Двухэтапная процедура поиска
Поиск сигнала малого уровня требует большой длительности накопления в каждой ячейке анализа. Снижение времени накопления ухудшает вероятностные характеристики, но, зато, сокращает время всего поиска.
Бытует мнение, что для ускорения процедуры поиска её можно разделить на два этапа. На первом этапе осуществлять накопление сигнала в каждой из ячейки поиска на "небольшом" времени накопления, с не очень хорошими характеристиками. В этом случае превышение порога поиска произойдёт в большом количестве ячеек, в том числе, в правильной. Выбирается некоторое количество максимумов, и уже для этих ячеек осуществляется накопление на втором этапе, на более длительном интервале накопления с хорошими вероятностными характеристиками. Т.к. количество ячеек на втором этапе меньше, общее время накопления может быть сокращено.
Этот вопрос сейчас иследует Иван Липа в своей магистерской работе. По предварительным расчётам этот подход эффективен только тогда, когда объём поиска на втором этапе относительно велик и составляет порядка 10% от полного поля поиска. Тогда можно достигнуть сокращения общего времени поиска до 2 раз.
Однако обычно подобные алгоритмы используют, когда объём поиска на втором этапе мал и составляет несколько десятков. В этом случае эффекта от данного подхода нет практически никакого.
Второй недостаток - "выигравшие" на первом этапе позиции идут не "по порядку" в поле поиска, а случайно. Из-за этого затруднительно на втором этапе использовать алгоритмы быстрого поиска, которые, как правило, расчитаны на анализ блока соседних позиций.
Подробности здесь: 24.10.2011_Двухэтапный_поиск_сигналов
Дальномерные коды с "хорошими" свойствами
Это идёт из связи. Для дальномерных кодов стараются подбирать псевдослучайные последовательности с низким уровнем боковых лепестков взаимно-корреляционной фукции.
Если рассмотреть набор просто случайных кодов, среднеквадратическое значение лепестков корреляционной функции будет равно , где N - количество символов кода на интервале накопления. Есть определённые коды, например, коды Касами, свойства которых намного лучше. Уровень боковых лепестков для кодов Касами равен .
По этому поводу было сломано много копий в рамках обсуждений перспективных сигналов. На самом деле это так, то пользы от этого никакой. Все эти "хорошие" свойства остаются в силе только до тех пор, пока два сигнала совпадают по частоте. Как только появляется разность частот сигналов, свойства кодов разрушаются и становяться такими же, как и у просто наборов случайных чисел . Причём разность частот должна быть совсем небольшой. Если рассматривать открытые сигналы на интервале 1 мс, то при разности частот уже 50 Гц свойства кодов полностью разрушаются. А разность частот - она есть всегда. Поэтому на деле пытаться сделать дальномерные коды с "хорошими" свойствами нет смысла.
[ Хронологический вид ]Комментарии
По поводу кодов. Построить бы график перехода 1/N в 1/sqrt(N) как функции ошибки по доплеру. Ведь 50Гц при хорошей энергетике - это очень много. Что будет при 1-5 Гц?
50 Гц - это не ошибка слежения за сигналом. Это разность частот принимаемого сигнала и другого сигнала из той же системы, который создаёт внутрисистемную помеху. Частоты этих сигналов могут отличаться независимо от энергетики. Это определяется взаимным движением спутников и положением потребителя. Графики такие я строил, у меня сейчас просто под рукой его нет.
Понял)) Что-то я чушь сказал, какие 1/N в режиме слежения)
Войдите, чтобы комментировать.