Перов А.И., Харисов В.Н. ГЛОНАСС. Принципы построения и функционирования (4-е издание, 2010) (1151865), страница 83
Текст из файла (страница 83)
Поэтому, подставив опорный сигнал, заранее «развернутый» во времени в (13.20), можно рассчитать все требуемые корреляционные суммы. Таким образом, из сказанного и соотношения (13.19) получается следующий алгоритм обработки: рассчитывается спектральная плотность опорного сигнала, развернутого во времени; рассчитывается спектральная плотность входного сигнала; выполняется поэлементное перемножение массивов спектральных плотностей входного и развернутого во времени опорного сигнала; выполняется обратное преобразование Фурье спектральной плотности, полученной на предыдущем шаге. Отметим, что «развернуть» во времени опорный сигнал удобно в спектральной области, выполнив комплексное сопряжение исходной спектральной плотности. Структурная схема алгоритма изображена на рис. 13.10.
Рис 13.10. Схема блока поиска с использованием дискретной свертки во временной области Следует отметить, что в описанном алгоритме поиска вычисляется циклическая свертка, то есть сигнал сдвигается циклически. Такая процедура хорошо 502 Аппаратура потребителей подходит для поиска периодических сигналов, которыми являются СТ-сигналы СРНС ГЛОНАСС. Если в сигнале используется непериодическая кодовая последовательность, например сигналы для специальных потребителей, то использование циклической свертки неправомерно.
Для таких сигналов можно использовать метод добавления нулей (Лего радйпд), заключающийся в том, что один из сигналов заполняет полностью интервал, анализируемый БФП, а другой сигнал заполняет только половину этого интервала, а его вторая половина заполняется нулями. В этом случае половина полученного в результате преобразования сигнала будет содержать требуемые корреляционные суммы. Поэтому, при использовании данного метода для непериодических требуется удвоенный объем БПФ. Описанный алгоритм имеет еще одну особенность. Существует теорема о том, что спектральная плотность произведения сигнала на комплексный гармонический сигнал равна частотному сдвигу спектральной плотности сигнала.
Это позволяет упростить поиск сигнала по частоте. Перед предварительным накоплением сигнал умножается на гармоническое колебание, частота которого соответствует центральной частоте ячейки поиска по частоте. Далее, после преобразования Фурье, полученная спектральная плотность циклически сдвигает по частоте, тем самым перебирая каналы поиска по частоте с шагом 1! Т.
После этого сформированная спектральная плотность умножается на комплексно сопряженный спектр опорного сигнала и выполняется обратное преобразование Фурье. В данном алгоритме прямое преобразование Фурье входного сигнала выполняется один раз, а обратное преобразование Фурье выполняется столько раз, сколько есть частотных каналов. Это дает выигрыш по вычислительным затратам. При программной организации поиска описанный алгоритм наиболее логичен ввиду простоты реализации БПФ. Однако сложность заключается в том, что наиболее эффективная его реализация требует хранения большого объема промежуточных данных, что может наталкиваться на ограниченностью аппаратных ресурсов НАП. Поиск с использованием БПФ ио час>поте Существует метод быстрого поиска, использующий последовательное сканирование по частоте и параллельное — по задержке.
Для параллельного поиска по частоте используется БПФ, реализуемое, как правило, программно. Схема блока поиска, основанного на данном методе, приведена на рис. 13.11. Входной сигнал умножается на опорный сигнал, представляющий собой дальномерной код для заданной ячейки разрешения по задержке. Далее используется коррелятор с квадратурными каналами, предназначенный для накопления на интервале То = 1/Л~,„р — — тоТ, длительность которого определяет- 503 Глава 13 сЯ зоной поиска по частоте 1з1д,р . ПаРаметРы опоРного сигнала, постУпаюЩего на коррелятор, соответствуют значениям задержки и частоты анализируемой ячейки разрешения.
Накопленные на требуемом количестве У1. интервалов длительностью То (равном %1 — — ф„ /6/,, где а ~'„— шаг анализа по частоте) отсчеты поступают на блок ВПФ, на выходе которого формируются У отсчетов, соответствующих одному значению задержки (одной ячейке поиска по задержке) и У1 значениям частоты (ячейкам поиска по частоте). У(ч, Текущее значение задержки Рис 13.11.
Схема блока поиска с использованием БПФ по частоте Достоинством данного метода является относительная простота его реализации в аппаратуре — нужно ввести специальный режим работы коррелятора с накоплением на интервалах времени То и организовать передачу накопленного массива в процессор, а остальное можно реализовать программно.
Преимущество использования данного метода по сравнению с методом последовательного поиска на корреляторах состоит в использовании эффективного алгоритма БПФ, имеющего вычислительную сложность М.1од2(М) вместо Х~ для прямого вычисления корреляционных сумм. Типичный объем БПФ в таком методе М = 64...256. Однако при практической реализации с учетом разрядности вычисления данный метод может проиграть методу, указанному выше методу. Перечисленные выше методы и алгоритмы иллюстрируют лишь некоторые из возможных подходов к построению блока быстрого поиска. Данные алгоритмы можно комбинировать друг с другом. В качестве примера на рис. 13.12 приведена схема алгоритма, представляющего собой комбинацию согласованного фильтра и БПФ по частоте.
504 Аппаратура потребителей Текущая анализируемая частота Й Рис 13.12. Схема блока поиска с использованием комбинированного алгоритма Данный алгоритм содержит несколько коротких согласованных фильтров, процессы с выхода которых подвергаются БПФ и на каждом такте работы на выходе БПФ формируются все корреляционные суммы по частоте для одного значения задержки. Таким образом, возможны различные варианты реализации алгоритма быстрого поиска. Следует отметить, что основная проблема проектирования таких алгоритмов заключается не в подборе наиболее подходящей структуры, а возможностями реализации того или иного алгоритма с учетом аппаратных ограничений.
Основные вопросы при этом: выбор разрядности представления чисел на каждом из шагов обработки, выбор темпа обработки для отдельных блоков и оптимизация требований к памяти. Как результат, при правильной реализации большинство из рассмотренных алгоритмов дают примерно одинаковый результат по скорости поиска. 13.7.2. Алгоритм работы и схема слежения за фазой сигнала В когерентном режиме работы НАП осуществляется слежение за фазой сигнала НС.
Подробное описание принципов построения следящих систем, в том числе и за фазой сигнала (ССФ), приведено в п. 6.3.2, из которого следует, что ССФ включает фазовый дискриминатор, сглаживающий фильтр и генератор опорного сигнала. В п. 6.3.4.1 показано, что в зависимости от условий работы ССФ оптимальными (с точки зрения обеспечения наилучшей точности слежения) являются различные типы ФД. В табл. 13.2 приведена итоговая сводка различных типов ФД и их основные характеристики.
505 Глава 13 "Габлица 13.2. Фазовые дискриминаторы и их характеристики Общие свойства Зависимость ДХ от фазовой ошибки Алгоритм работы Оптимальный ФД при произвольном отношении сигнал/шум. Крутизна ДХ пропорциональна ам- плит де А. Большие вычислительные зат аты. -1Ь(/р) Ор Близок к оптимальному при большом отношении сигнал/шум. Крутизна ДХ пропорциональна амплиде А . Минимальные вычислительные за аты.
яп(с>ф) — ядп(1р)Др Близок к оптимальному при малом отношении сигнал/шум. Крутизна ДХ пропорциональна А . Уме- 2 енные вычислительные зат аты. яп (2аф) — 1 Ц, Субоптимален, но имеет хорошие характеристики при большом и малом отношении сигнал/шум. Крутизна ДХ не зависит от амплитуды А. Большие вычислительные затраты и необходима проверка операции деления при ошибке +90' (деление на ноль . 1фЮФ) Главный угол арктангенса. Оптимален в смысле оценок максимума функции правдоподобия [5.1, 5.21 при произвольном отношении сигнал/шум.
Крутизна ДХ не зависит от амплитуды А. Наибольшие вычислительные зат аты. — Агап Яр/1р ) Рис. 13.13. Дискриминационные характеристики ФД 506 Дискриминационные характеристики ФД приведены на рис. 13.13 (по оси ординат характеристики нормированы к размерности «градус»). Аппаратура потребитепей (1е — 1~ ) 1р + © Д~ ) Я, — близок к оптимальному при когерентном приеме при малом отношении сигнал/шум; хорошие характеристики (кривая 4 на рис. 13.15) при модуле ошибки слежения менее 0,5г,; минимальные вычислительные затраты. Рис.
13.15. Дискриминационные характеристики дискриминаторов задержки сигнала Типичная схема ССЗ с одним из типов дискриминаторов задержки приведена на рис. 13.16. Коррелятор Рис. 13.16. Схема ССЗ В схеме рис. 13.16 управление положением опорного сигнала осуществляется приращением оценки задержки сигнала на соседних тактах работы ССЗ.