автореферат диссертации по радиотехнике и связи, 05.12.14, диссертация на тему:Нейтрализация многолучевых помех в РНС космического базирования

кандидата технических наук
Соколов, Андрей Андреевич
город
Санкт-Петербург
год
2012
специальность ВАК РФ
05.12.14
цена
450 рублей
Диссертация по радиотехнике и связи на тему «Нейтрализация многолучевых помех в РНС космического базирования»

Автореферат диссертации по теме "Нейтрализация многолучевых помех в РНС космического базирования"

На правах рукописи

СОКОЛОВ Андрей Андреевич

НЕЙТРАЛИЗАЦИЯ МНОГОЛУЧЕВЫХ ПОМЕХ В РНС КОСМИЧЕСКОГО

БАЗИРОВАНИЯ

Специальность: 05.12.14 - Радиолокация и радионавигация

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата технических наук

^ и п! | | ¿и;

Санкг-Пеггербург - 2012

005018308

005018308

Работа выполнена в Санкт-Петербургском государственном электротехническом университете «ЛЭТИ» им. В.И.Ульянова (Ленина) на кафедре радиотехнических

Научный руководитель:

доктор технических наук, профессор Ипатов Валерий Павлович

Официальные оппоненты:

доктор технических наук Пономаренко Борис Викторович, начальник сектора ОАО «Всероссийский научно-исследовательский институт радиоаппаратуры (ВНИИРА)» (г. Санкт-Петербург)

кандидат технических наук Гугин Виталий Семенович, научный сотрудник СПбГЭТУ «ЛЭТИ», кафедра радиоэлектронных средств

Ведущая организация:

ОАО «Российский институт радионавигации и времени (РИРВ)»

Защита состоится «_18_» апреля. 2012 г. в Л_ часов на заседании совета по защите докторских и кандидатских диссертаций Д 212.238.03 Санкт-Петербургского государственного электротехнического университета «ЛЭТИ» им. В И Ульянова (Ленина) по адресу: 197376, Санкт-Петербург, ул. Проф. Попова, 5

С диссертацией можно ознакомиться в библиотеке СПбГЭТУ.

Автореферат разослан " 16 " марта 2012 г.

Ученый секретарь

Совета по защите докторских и кандидатских диссертаций Д 212.238.03

Д.тл.

Баруздин С А

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальпость темы. В настоящее время спутниковая навигация приобрела огромную популярность. Если раньше использование космических аппаратов в целях определения координат объекта считалось исключительно прерогативой военных, то в последние годы глобальные навигационные спутниковые системы (ГНСС) прочно вошли в повседневную жизнь сотен тысяч обычных людей по всему миру. На современном этапе гражданскому потребителю доступны две равноценные и взаимодополняющие ГНСС: GPS и ГЛОНАСС. На стадии пробных испытаний находится ГНСС Galileo, продвигаемая Европейским сообществом. В ближайшем будущем Китай планирует развернуть ГНСС двойного применения Compass. Кроме того, Япония близка к старту практического развертывания региональной навигационной спутниковой системы QZSS (Quasi-Zenith Satellite System). Подобный рост количества систем спутниковой навигации в пределах одного частотного диапазона L ведет к спектральным коллизиям как между ними самими, так и с системами в смежных спектральных диапазонах. В этом плане традиционные для ГНСС навигационные сигналы (НС) с бинарной фазовой манипуляцией, обладающие широким спектром, оказываются не лучшими с точки зрения электромагнитной совместимости и альтернативой им могут явиться спектрально-эффективные сигналы, в частности, сигналы с минимальной частотной модуляцией (МЧМ). При этом правомочен вопрос о тактических характеристиках ГНСС, использующих такие сигналы, и прежде всего о точностных показателях. Таким образом фундаментально задача сводится к оценке потенциальной точности измерения радионавигационных параметров (РНП), в частности, задержки навигационного сигнала.

Проблема обеспечения качественных навигационных измерений в условиях многолучевости является одной из ключевых в вопросах функционирования ГНСС. Поэтому сформулированную выше задачу целесообразно решать именно для сценария интерференциошюй картины на приемной стороне. При этом с практической точки зрения важным также является вопрос об устойчивости спектрально-эффективных сигналов к многолучевости в сценариях приема с помощью следящих систем.

Повсеместной практикой является моделирование многолучевости одной отраженной репликой прямого сигнала. Такое допущение редко соответствует действительности, что побуждает к усложнению модели в сторону увеличения числа лучей для корректного рассмотрения вопросов потенциальной точности оценки РНП. В то же время вопрос применимости однолучевой модели в разработках навигационных приложений остается открытым, поскольку использование сложных моделей отражений в реальных устройствах не всегда возможно, а адекватность упрощенных версий должна быть в достаточной степени аргументирована. В этом плане полезной может оказаться разработка алгоритмической процедуры нейтрализации многолучевости, поскольку именно алгоритмические методы наиболее чувствительны к модели сигналов, заложенных в их основе.

Названные факторы обуславливают актуальность исследования потенциальной точности оценки РНП спектрально-эффективных сигналов в условиях многолучевого распространения, анализа устойчивости таких сигналов к многолучевости и разработки алгоритмических методов ее нейтрализации.

Цель работы состоит в исследовании потенциальных характеристик измерения времени в ГНСС в том числе для сигналов спектрально-эффективного формата, а также

разработке методов нейтрализации многолучевых ошибок на основе множественных измерений. Для достижения поставленной цели в диссертационной работе предстояло решить следующие задачи:

1. На основе литературных источников рассмотреть методы борьбы с многолучевостью.

2. Принять для дальнейшего анализа конкретные формат спектрально-эффективной модуляции.

3. Аналитически рассчитать точности измерения РНП при оптимальном приеме и ограниченном спектре сигнальных чипов для различных сценариев присутствия многолучевосга, моделируемой одним и двумя лучами.

4. Разработать алгоритмическую процедуру нейтрализации негативных эффектов многолучевосга и детально исследовать ее характеристики путем компьютерного моделирования.

5. Провести экспериментальные исследования с использованием записей реальных сигналов двухчастных навигационных приемоизмерителей ГНСС.

Методы исследования. Для решения поставленных в диссертационной работе задач были использованы методы теории сигналов, теории вероятностей и математической статистики, гармонического анализа, численные методы и методы математического моделирования. Экспериментальные исследования выложены методом полунатурного моделирования с использованием записей реальных сигналов навигационных приемоизмерителей в ГНСС.

Научная новизна. Новыми являются следующие результаты:

1. Обобщены и систематизированы результаты расчетов потенциальных границ точности оценки задержки НС для однолучевой модели отражений.

2. Получены аналитические выражения для расчета дисперсии оценки РНП при оптимальном приеме НС для двухлучевой модели отражений. Указанные равенства получены для ограниченного в частотно-селективном тракте приемника чипа стандартных и спектрально-эффективных сигналов.

3. Установлено влияние второго отраженного луча на потенциальную точность оценки РНП как для традиционных, так и для спектрально-эффективных сигналов.

4. С целью оценки устойчивости к многолучевости рассчитаны огибающие многолучевой ошибки для спектрально-эффективных сигналов в сценариях измерения РНП с помощью классического ранне-позднего дискриминатора и оптимального дискриминатора слежения за характерной точкой. Даны рекомендации по выбору параметров дискриминатора.

5. Разработан алгоритмический метод нейтрализации многолучевости, основанный на множественных наблюдениях НС в четырех каналах навигатора - стандартной и высокой точности в двух частотных диапазонах - и однолучевой модели отражений.

6. Экспериментально подтверждена работоспособность упомянутого алгоритма, путем обработки реальпых данных, полученных от приемоизмерителя ГНСС.

Положения, выносимые на защиту.

1. В присутствие многолучевости для узкополосного частотно-селективного тракта приемника (полоса ограничена 1-2 лепестками) при классической минимальной

частотной модуляции выигрыш в точности достигает 1,8 раз в сравнении с используемым в ГНСС форматом модуляции ВРвК..

2. Спектрально-эффектавные сигналы с составной структурой чипа при широкополосном тракте приемника обеспечивают выигрыш в дисперсии оценки задержки НС до 11 дБ в сравнении со стандартными сигналами ГНСС.

3. В рамках традиционного слежения в петле АПВ спектрально-эффективные сигналы обеспечивают меньшую ошибку многолучевое™ при использовании как классического ранне-позднего дискриминатора, так и оптимального дискриминатора слежения за характерной точкой НС.

4. Усложнение модели отражений за счет увеличения числа лучей заметно увеличивает дисперсию оценки задержки сигнала при оптимальном приеме Этот эффект наиболее выражен в случае стандартных сигналов ВРБК и менее заметен в случае спектрально-эффективных сигналов.

5. Навигационные измерения в двух частотных диапазонах в каналах СТ и ВТ позволяют устранить многолучевую ошибку с помощью алгоритмической процедуры, основанной на однолучевой модели отражений.

6. Экспериментально доказано, что однолучевая модель' отражений обладает определенной степенью адекватности и может быть использовала при разработке практических приложений ГНСС.

Практическая ценность работы Основным практическим выходом работы является подтверждение конкурентоспособности спектрально-эффективных сигналов в плане потенциальной точности и устойчивости к многолучевости, аргументация в пользу адекватности однолучевой модели в практических навигационных приложениях и алгоритмический метод нейтрализации многолучевости, не требующий аппаратных изменений стандартного навигатора.

Апробация работы. Результаты работы докладывались и обсуждались на 64-й - 66-й Научной сессии, посвященной Дню Радио (СПб, 2009 - 2011); на VI и VII международной молодежной научно-технической конференции «Современные проблемы радиотехники и телекоммуникаций» (Севастополь, 2010, 2011); на 62-й - 64-й конференции профессорско-преподавательского состава СПбГЭГУ «ЛЭТИ» (СПб, 2009 - 2011); на IX международном

и выставке по электромагнитной совместимости и элекхромапшпюй экологии (СПб, 2011); на Научной сессии «Школа радиэлекгроники» (СПб, 2010).

Публикации. Всего по теме диссертации опубликовано 10 работ. Из них 3 работы опубликованы в центральных рецеширусмых научных журналах, рекомепдованных перечнем ВАК, 7 работ содержатся в сборниках материалов международных, всероссийских и межвузовских научных конференций.

Сгруктура и объем диссертации. Диссертация состоит из введения, четырех глав заключения, библиографического списка. Она изложена на 167 страницах машинописного текста, содержит 81 рисунок, 14 таблиц, библиографический список включает 100 наименований.

ОСНОВНОЕ СОДЕРЖАНИЕ РАБОТЫ

Во введении обоснована актуальность темы, сформулированы цели, направления исследовании и основные научные положения, показана научная новизна и практическая

значимость полученных результатов, содержится краткий обзор литературы по теме диссертации, сформулирована цель работы.

В первой главе проведен обзор современного состояния ГНСС, приведена краткая историческая справка Далее сжато излагаются необходимые для дальнейшего сведения о процедурах измерения времени прихода сигналов в приемниках ГНСС, а также описывается явление многолучевости и приводится обзор известных методов борьбы с ней.

В бюджете погрешностей позиционирования по сигналам ГНСС значительный удельный вес принадлежит ошибкам, обязанным эффектам многолучевого распространения радиоволн на трассе спутник - потребитель.

Известные методы борьбы с многолучевостью основаны на аппаратном усложнении навигационного оборудования, прежде всего антенно-фидерной системы и следящих контуров навигатора. Периодичность повторения многолучевых помех позволяет использовать для их устранения фильтры с большими постоянными времени, а наращивание числа корреляторов в одном канале приемника позволяет приблизится к потенциальным границам точности оценки РНП. В то же время такой подход ведет к увеличению стоимости аппаратуры и снижает интерес со стороны массового потребителя к навигационному оборудованию. В описанной ситуации могут иметь перспективу способы алгоритмической коррекции навигационных измерений, искаженных многолучевой помехой, основанные на обработке измерений стандартного приемника.

В качестве итога можно сделать следующие выводы:

1. При совершенствовании и модернизации ГНСС важнейшим аспектом исследований является проблема устойчивости к многолучевости.

2. На сегодняшний день не существует универсального метода борьбы с многолучевостью.

3. Большинство методов ориентировано на значительное усложнение аппаратной части НАЛ.

4. На данном этапе эффективные чисто алгоритмические методы борьбы с многолучевостью находятся на стадии перспективных разработок.

Во второй главе описывается и подробно рассматривается проблема спектральной совместимости ГНСС ГЛОНАСС с сигналами радиоастрономических наблюдений. В настоящее время просачивание спектральных компонент навигационных сигналов в радиоастрономическое окно строго регламентируется международным союзом электросвязи (МСЭ). Для системы ГЛОНАСС максимально допустимая плотность потока мощности у поверхности Земли установлена на уровне -194 дБ Вт/м2 в полосе ^ = 20 кГц на один навигационный космический аппарат (НКА). В рамках существующих форматов модуляции столь жесткие требования можно соблюсти лишь на основе грубой режекции спектральных компонент НС в радиоастрономическом окне. Кардинальным решением этой проблемы является применение спектрально-эффективных сигналов, в частности, сигналов с минимальной частотной модуляцией МЧМ.

Сопоставление спектров МЧМ сигналов и стандартных БФМ сигналов ГНСС, показывает, что 99% полоса МЧМ сигнала на порядок меньше аналогичной у БФМ. Проблема комплексирования сигналов открытого и санкционированного доступа (ОД и СД) на одной несущей в рамках МЧМ может быть решена путем представления чипов в каждом канале в виде последовательности более коротких чипов («микрочипов»), представляющих собой полуволну косинуса. При этом узкополосный чип сигнала

открытого доступа формируется из микрочипов одной полярности, а более короткий широкополосный чип санкционированного сигнала представляет собой отрезок гармоники (рис. 1). При такой структуре чипов сложение двух потоков микрочипов, предварительно сдвинутых на половину длительности микрочипа, даст МЧМ сигнал.

Составной чип

<----Л ■>,

• ГУУЛ А .

Рис. 1. К построению составных чипов

В заключительной части главы описана методика подбора параметров составного чипа исходя из требований, предъявляемых к спектру сигнала с учетом совместимости с сигналами радиоастрономии. На основе этой методики составлена таблица возможных форматов модуляции (табл. 1), отличающихся длительностью микрочипа.

Таблица 1. Возможные опции МЧМ-форматов сигналов ОД и СД в диапазоне 1Л

№ опции ..2 - 3 4 5

Сигнал ОД MSKi.li ; :М8К(5) М5К(6) МБК(8) мэкоо)

Сигнал СД вос- ШКСЩ вос- МЖ(2,5;2,5) ВОС- М8К(3,3) ВОС-МЭК(4,4) ВОС-М5К(5,5)

пт ■ 2. . 3 4 5

с 2 '2,5 ;. 3 4 5

птс - 2пт ■4; 5 6 8 10

"р = "„Ч 1 1 1 1 1

Несущая /, МГц 1 562= = 1597,926 1563,5/ = = 1599,460 1561/, = = 1596,903 1562,5/.= = 1598,437 1559,5/ = = 1595,368

Занимаемый диапазон, МГц 1593,6 -1602,3 1594 — -1605 1590,1-1603,7 1589,4-1607,4 1584,1-1606,5

Полоса сигнала СД, МГц 8,7 13,5 18 22,3

Требуемая глубина режекции, дБ 2-3 - 5 4 8 8

Резюмируя содержание главы можно сделать следующие выводы: 1. Показано, что спектрально-эффективная модуляция, в частности МЧМ, может с успехом использоваться в системах спутниковой навигации для предотвращения

проникновения спектра сигналов в радиоастрономическое окно, при этом требуемая глубина режекции компонент НС в нем не составит более 6 8 дБ.

2. Задача комплексирования сигналов открытого и санкционированного доступа может быть решена путем представления их чипов в виде последовательности микрочипов.

3. Варьирование формы составного чипа дает возможность гибко управлять спектрами навигационных сигналов, в частности смещением сигналов на оси частот с заданными характеристиками ширины спектра сигналов ОД и СД.

4. Квадратурное уплотнение сигналов ОД и СД легко осуществить простым суммированием, предварительно сдвинув их на половину длительности микрочипа. Полученный таким образом НС будет свободен от амплитудной модуляции и обладать на порядок более узкой 99-процентой полосой, чем стандартный ФМ сигнал.

В третьей главе приведены результаты сравнительного анализа спектрально-эффективных и стандартных сигналов ГНСС с точки зрения потенциальной точности оценки РНП, а также устойчивости к многолучевой помехе при измерении времени прихода сигнала с помощью традиционных следящих измерителей. Расчеты потенциальной точности выполнены как для однолучевой, так и для двухлучевой модели отражений. Проведен сравнительный анализ подходов, различающихся количеством лучей в модели отражений.

Для суждения о применимости тех или иных сигналов необходимо прежде всего оценить потенциально достижимую точность навигационных измерений с учетом многолучевого характера приема. При расчете потенциальных границ точности, т. е. минимально достижимой дисперсии измерений РНП, в качестве эталона принимались границы точности для стандартной БФМ, иными словами для сигнала с прямоугольной формой чипа, прошедшего частотно-селективный тракт.

При двухлучевой модели отражения многолучевой помехе соответствует четыре параметра - а, в и 6,9- соответственно амплитуда и задержка первого и второго лучей. Для количественной оценки влияния многолучевости на точность оценки задержки НС т были проанализированы три варианта приема, отличающиеся полнотой априорных сведений о параметрах многолучевой помехи:

1. многолучевая помеха полностью детерминирована;

2. оценке подлежат задержки прямого и отраженного сигналов, при известной амплитуде последних;

3. оценке подлежат все параметры совокупного сигнала.

Исходя из определения границы Крамера-Рао, дисперсию оценки задержки прямого сигнала х можно получить путем обращения матрицы Фишера с удержанием её первого диагонального элемента.

Из определения обратной матрицы следует, что

... дяХМ„

уаг,{т}«-—, д»1, (1)

где г = 1, 2, 3 - номер сценария, (ЫФ - определитель матрицы Фишера, А/„ - матрица, полученная из Ф удалением первых строки и столбца (т.е. с!е1 Мп - алгебраическое

дополнение элемента Ф„ в (ЫФ), а д2 - отношение сигнал-шум. Потери в дисперсии из-за многолучевости удобно оценить с помощью коэффициента

Т, ="

уаго М

где уаг0 {т} - дисперсия оценки х в отсутствие многолучевости, т. е. когда а = 0 и 6 = 0 Опираясь на (1) можно получить следующие выражения для трех сценариев: уаг,{т} _

У, =

У 2 =

1

1 + 0о2 +2а0?Р± + 2Ь0?Ш + 2аЛ

Р (0) Р"(0) 0 0 р"(0)

?2р"(0)(ФвеФ9

узг2{Т} ^__

уаг0{т} ~ Ф„ФтФэз -Ф„Ф^Э -Ф,20Ф83 + 2ФяФ98Фд3 - Ф^9Фв(

Уз =

уаг,{т}

(2)

(3)

(5)

(6)

уаг„{т}'

где р(т) -автокорреляционная функция сигнала.

Для однолучевой модели соответствующие равенства для всех трех исследуемых сценариев могут быть получены непосредственно из выражений для двухлучевой модели путем подстановки Ь0 = 0:

1 1

Т) :

_1_

1 + а1 + 2а.

р'(е.)

Р'(0)

ь:

р"(60)

р"(0)

Уз '

1 + [р'(е0)]2 ГРЧ6„)|

Р'(0) 1. р"(0) ]

(7)

На рис. 2 приведены кривые, построенные согласно (7) для третьего варианта приема сигнала с прямоугольной формой чипа при однолучевой модели отражений для БМФ и МЧМ. Видно, что с ростом полосы частотно-селективного тракта IV = к / А. где к - количество лепестков, а Д - длительность чипа, потери в дисперсии для БФМ сигнала быстро уменьшаются, тогда как для МЧМ это падение заметно медленнее. При этом при IV = 1 / А МЧМ обеспечит выигрыш в точности в районе малых, наиболее опасных задержек отражения.

Б МФ 1- 2— т-г 3— ч¥д=ю

1

\

1. 3

р ЛЧМ 1-УУД=1 2—т=2 3— \л/д=1о

?

Дз

...........

0.2 0.4 0.6 0.8 1 1.2 1.4

Рис. 2 - Потери в дисперсии т при отсутствии данных об отражении. Такой эффект легко объясним тем, что с расширением полосы крутизна фронтов чипа БФМ сигнала резко возрастает, тогда как форма МЧМ чипа практически не

меняется. При этом в области малых задержек потери неограниченно растут, поскольку прямой и отраженный лучи практически неразрешимы.

Аналогичные кривые, но для двух лучевой модели представлены на рис. 3. В качестве параметра для построения кривых выступает разность между задержками первого и второго отраженного луча. Видно, что кривые группируются в достаточно плотной области, при этом для определенных вариантов задержек потери в точности заметно превосходят аналогичные, рассчитанные в рамках однолучевой модели

PSK. Variance losses, \№,=1

PS К. Variance losses. Ш=2

Рис. 3

гг а 6

-1 ютери в точности в условиях полной неопределенности.

Двухлучевая модель.

БФМ сигнал

Рис. 4 - Потери в дисперсии оценки задержки при отсутствии данных об отражении. МЧМ сигнал. Двухлучевая модель.

Подобные наблюдения можно провести и в случае МЧМ сигнала (рис. 4). Сравнение с рис. 2 показывает, что с ростом полосы кривые, соответствующие различным задержкам отраженных лучей, плотно группируются, вокруг кривой для однолучевой модели. Таким образом, определенную коррекцию в сторону увеличения дисперсии измерений РНП можно наблюдать в наиболее узкополосном варианте (спектр ограничен одним лепестком).

MSK. Variance losses. WA=1

1 — (W=0 054

2 —(W=025A

3 — (80-е0)-0.354

4 - (s0-e^.st

5 _(90-90)=O.B5i S (8.^0=1.054

MSK. Variance losses Ш=2

В табл. 2 сведены контрольные значения по всем трем сценариям. В качестве контрольных значений б0т!„, выступают такие значения б, при которых потери в точности в энергетическом пересчете составляют -3 дБ.

Таблица 2 - Контрольные значения задержки в сравнении с однолучевой моделью

Оцениваемые параметры w д ФМ МЧМ

1 луч 2 луча 1 луч 2 луча

60mi„/A диапазон 00min / Д диапазон e„min / Д

т 1 0.25 [0.1 0.33] 0.285 [0.07 0.29]

2 0.1 [0.054 0.17] 0.24 [0.07 0.2]

10 0.019 [0.012 0.022] 0.22 [0.05 0.18]

т,е,9 I 0.25 [0.25 0.55] 0.18 [019 0.36]

2 0.11 [0.1 0.2] 0.15 [0.15 0.22]

10 0.02 0.02 0.11 [0.11 0.14]

т да, о, ¿> 1 _0.35 [0.42 0.64] 0.29 [0.25 0.77]

2 0.15 [0.134 0.276] 0.21 [0.19 0.27]

10 0.021 [0.21 0.046] 0.15 [0.15 0.17]

Для сигналов ОД с составной структурой чипов и однолучевой модели кривые для квадрата нормированной дисперсии представлены на рис. 5.

С/А^Дс=5 МЭК УУй =5

Рис. 5 - Зависимости дисперсий оценки т от задержки отражения при IV Ас =5. Однолучевая модель. Сигнал ОД.

Таблица 3 - Отношение дисперсий оценки т по С/А коду GPS и сигналу ОД МЧМ

\&с "X 1 2 3 4 5 6 7 8 9 10

4 0,832 0,923 1,194 4,894 5,941 4,964 4,261 3,994 3,753 3,380

5 0,824 0,878 0,999 1,326 5,848 7,645 6,588 5,765 5,137 4,956

6 0,820 0,855 0,927 1,073 1,452 6,789 9,354 8,247 ,7,334 6,603

В 0,816 0,835 0,870 0,929 1,028 1,215 1,689 8,657 12,776 11,625

10 0,814 0,826 0,847 0,880 0,930 1,005 1,126 1,349 1,915 10,502

Из рисунка видно, что при определенном количестве микрочипов в чипе выигрыш в дисперсии может достигать 11 дБ. В частности, для сигнала с 4 микрочипами на чип

выигрыш в дисперсии появится начиная с IV = 3 / Д, а для 10 микрочипов - при IV = 6 / Д (см. табл. 3). При двухлучевой картине (рис. 6) эти значения несколько изменятся из-за влияния второго отраженного луча. В частности для птс = 4 это значение составит УГ = 5/Д,адая л„=10 -^=10/Д.

MSK(4) W,A=3

MSK(10). W Д =3

Рис. 6 - Зависимости дисперсий оценки х от задержки отражения при различном количестве микрочипов в чипе. Двухлучевая модель. Сигнал ОД.

Сигнал СД необходимо сравнивать с Р-кодом GPS, с длительностью чипа Ър . В этом случае расчет дисперсии показывает, что сопоставимый с Р-кодом GPS по полосе спектрально-эффективный сигнал при пт =10 выигрывает у последнего во всей области задержек отражения в случае, когда фильтр пропускает лишь основной лепесток спектра (1^8^=1) (рис.7). В этой же ситуации выигрыш в заметной части зоны задержек отражения имеет и формат птс= 8. С удвоением же полосы фильтра менее уязвимым к отражению при задержках последнего меньших длительности микрочипа оказывается Р-код GPS.

MSK-BOC. W5 =1 MSK-BOC. WS =2

1.6 1.4 ™ о. 1-2

!' g 0.8

I ! V 1 V2

2 3 - V пт=5

V \

\ \

V 2 Р- код /

i

0.2 0.4 0.6 0.8

Vv а

Рис. 7 - Зависимости дисперсий оценки т. Однолучевая модель. Сигнал СД.

Переход к двухлучевой картине показывает (рис. 8), что введение в модель второго луча практически не влияет на потенциальные характеристики измерений задержки. Лишь в наиболее узкополосном варианте при очень близких задержках двух отраженных лучей

Рис. 8 - Зависимости дисперсий оценки х. Двухлучевая модель. Сигнал СД.

Реализация оптимального устройства, измеряющего время распространения НС в условиях многолучевости, сопряжена с чрезмерными аппаратными затратами. По этой причине производители оборудования отдают предпочтение разнообразным следящим системам, ключевым элементом которых является временной дискриминатор. Вследствие влияния отраженных реплик сигнала характеристика дискриминатора смещается относительно нуля на величину, которая и является ошибкой измерения задержки из-за многолучевости. Зависимость этой ошибки от задержки отражения называют огибающей многолучевости.

В главе рассмотрены два типа дискриминаторов: стандартный ранне-поздний и дискриминатор слежения за характерной точкой сигнала. На рис. 9 изображены дискриминационные характеристики последнего для различных типов сигналов.

БМФ МЧМ

Рис. 9 - Зависимости дисперсий оценки т. Двухлучевая модель. Сигнал СД. Расчет огибающих для дискриминатора слежения за характерной точкой свидетельствует, что вплоть до ограничения спектра сигнала пределами 2W = 10 / 6 = 100 МГц формат сигнала СД не прошрьгеает Р-коду GPS, а в канале ОД в

узкополосном варианте огибающие абсолютно одинаковы. Объяснение этому заключено в том, что в пределах реальной полосы аналоговой части приемника для МЭК чипа

максимальная крутизна фронта остается равной л/2Р% / 5 при \У > 1,2 / 5, а у прямоугольного чипа крутизна прямо пропорциональна IV.

На рис. 10 изображены огибающие для ранне-позднего дискриминатора. Как видно из рис. 10, с расширением полосы частотно-селективного тракта слежение с использованием коррелятора с малым расстоянием между опорами становится более устойчивым к многолучевой помехе в случае МЭК по сравнению с ВРвК. Это обусловлено, прежде всего, тем, что при широкой полосе составной чип передается с меньшими искажениями, и, следовательно, его корреляционная функция (КФ) имеет гораздо более острый пик, чем КФ прямоугольного чипа. Графики рис. 10 свидетельствуют, что, уже начиная с IV = 2/ Ас, максимальный уровень огибающей для М8К в зоне задержек 9 = &с / 2 лежит заметно ниже аналогичного в случае ВР8К.

пхйрй)! етевре тиШра^ етеЬре

Рис. 10 - Огибающие многолучевой ошибки для различной ширины полосы фильтра для различных форматов модуляции.

При полосе IV = 10 / Дс, МЭК существенно выигрывает у стандартной ВРБК. Причем для этой полосы фильтра М8К(8) дает ошибку (по пиковому значению) практически в четыре раза меньшую, чем МЭК(4) и ВРвК. При этом для М8К модуляции с ростом задержки многолучевого сигнала эта ошибка быстро падает (в четыре раза за длительность чипа), чего не наблюдается в случае ВРБК, где ошибка равномерна на всем диапазоне задержек в пределах чипа.

Резюмируя результаты главы можно заключить, что

1. Сигналы спектрально-эффективных форматов не проиг рывают либо выигрывают у стандартных сигналов ГНСС в плане потенциальной точности при реальных полосах частотно-селективного тракта приемника.

2. Для сигналов с составной формой чипов выигрыш в дисперсии измерений может достигать 11 дБ в энергетическом пересчете.

3. Усложнение модели многолучевости за счет ввода дополнительного луча увеличивает дисперсию оценки задержки. Этот эффект наиболее выражен в случае БФМ, и менее заметен в случае МЧМ. Для сигналов с составной структурой чипа это влияние минимально.

4. при охвате фильтром более двух лепестков спектра сигналы со спектрально-эффективной модуляцией заметно выигрывают в плане устойчивости к многолучевой помехе. Расчет огибающей показывает, что при десяти лепестках это преимущество может достигать четырех раз.

Четвертая глава диссертации посвящена разработке и исследованию алгоритмического метода нейтрализации многолучевости, основанного на однолучевой модели отражений и множественных измерениях навигационного сигнала.

При проектировании реального устройства в основу построения соответствующего измерителя придется положить ту или иную модель многолучевой помехи, снабдив принимаемое допущение аргументами в пользу ее адекватности. Естественно, в плане упрощения и коммерциализации навигатора наиболее удобна простейшая - однолучевая модель отражений, так что демонстрация ее приемлемости для достаточно представительного набора реальных сценариев весьма актуальна с точки зрения решений по выбору архитектуры аппаратуры потребителя. С другой стороны, если бы такая модель в точности отвечала истинной ситуации, для устранения многолучевой погрешности можно было бы воспользоваться различием ошибок, вносимых отражением в результаты измерений радионавигационного параметра по грубой и точной дальномерным шкалам, а также по сигналам двух частот системы ГЛОНАСС. Алгоритм подобного тина, не предполагающий вмешательства во внутреннее функционирование стандартного навигатора и оперирующий лишь оценками запаздывания, выдаваемыми следящими петлями последнего, предлагается и детально анализируется в четвертой главе.

Система из четырех уравнений упомянутых выше измерений имеет вид:

У и = * + Ма>0) + и,,

И, + И 1Р(а,в) + п,р

, =у = те + М(а,й,е) + п (8)

где первое слагаемое правых частей - истинная задержка НС, второе - смещение нуля дискриминатора из-за переотражения, третье - шумовая ошибка измерения. Разрешив эту систему, можно найти вектор неизвестных параметров х = (т,а,Р,8)г, а, значит, устранить многолучевую компоненту из результата измерений задержки НС. Стандартным подходом к решению такой системы является применение итерационной процедуры с линеаризацией па каждом шаге. Линеаризованная система (8) имеет вид

У-Уо=С(х0)(х-х0), (9)

где х„ =(т0,а0,р0,90)г -начальное приближение, у0 =тое + М(ао,ро,0о), а С(х0) - 4x4 градиентная матрица, первый столбец которой состоит из единиц, а остальные - из частных производных соответствующих смещений ц по а,р и 0 в точке х = х0. Если

при этом градиентная матрица имеет полный ранг (равный четырем), система (8) имеет решение

=*о + С"'(у-Уо). (10)

используемое далее вместо х0 как начальное приближение на следующей итерации, полностью повторяющей предыдущую и т.д.

Таким образом работа алгоритма характеризуется параметрами сходимости, а именно скоростью и зоной сходимости. Зона сходимости - это область в пространстве а, Р, 9, в которой алгоритм сходится к решению с заданной наперед точностью.

Первый этап исследования заключался в компьютерном моделировании работы алгоритма. Второй - в экспериментальном подтверждении его работоспособности. В результате моделирования алгоритма было выявлено, что:

1. Алгоритм обладает приемлемым быстродействием, сходясь в среднем за пять итераций. В пределах зоны сходимости дли подавляющего большинства комбинаций значений параметров многолучевого сигнала оттслонение требуемого числа итераций от этого среднего не превосходит одного-двух шагов..

2. При оптимальном выборе начальных значений алгоритм успешно сходится в »43%-61% опытов при критерии точности решения Дц-сЮ^Д, При этом

зона, где алгоритм не сходится, ожидаемо соответствует высоким значениям амплитуд на обеих частотах в области малых задержек отражения. С другой стороны типичные значения амплитуды многолучевых компонент |<х|, |р| < -3 дБ, что увеличивает вероятность успешного исхода.

3. При значениях х„ в середине соответствующих диапазонов а, р, 6 вероятность сходимости имеет наибольшее значение.

Далее алгоритм был применен к случаю фильтрованного чипа. Анализ полученных данных показал, что:

1. По-прежнему лучшим выбором начального приближения являются значения из середины рассматриваемого диапазона параметров многолучевости.

2. Алгоритм для случая конечной полосы при а 5 0.75 обладает значительно большей зоной сходимости и вероятностью успешного исхода, при этом последняя при а < 0.45 составляет более 90%, что превосходит аналогичные значения для алгоритма с РГ -> со более, чем на 30%.

3. В отличие от случая с бесконечной полосой зона сходимости практически не зависит от выбранного значения Дц (10~2Д1р или Ю^А.Д

4. По всем параметрам алгоритм инвариантен к измененшо ширины полосы частотно-селективного тракта.

5. При а > 0.85 вероятность и область сходимости резко уменьшается.

6. В зоне сходимости среднее количество затрачиваемых итераций не превосходит 6-7 шагов.

Главный вывод: преобразованный в соответствии с условиями ограничения спектра алгоритм оказался более работоспособен, чем его прототип, разработанный для бесконечной полосы частотно-селективного тракта, практически по всем показателям. Главным недостатком оказалась невозможность коррекции при а >0.9. Однако учитывая, что в реальности отраженные лучи приходят сильно ослабленными по амплитуде, указанный недостаток не относится к критическим.

Далее были исследованы шумовые характеристики алгоритма. Показано, что в типовом диапазоне отношений сигнал-шум алгоритм полностью оправдывает свое применение.

В заключительной части работы приведены результаты полунатурного эксперимента. Алгоритм бы применен к реальным измерениям, записанным с помощью

стандартного приемника с апертурой в один чип. Измерения радионавигационных параметров формировались с темпом 1Гц с 9 ч 46 мин 51 с до 10 ч 28 мин 4 с. Наблюдения проводились за спутником ГНСС ГЛОНАСС №12 (номер литеры «-1») Таким образом в ходе эксперимента было обработано 2474 записи, для каждого из четырех каналов.

По итогам эксперимента можно сформулировать следующие выводы:

1. Алгоритм сходится практически в 70% измерений, при этом скорость сходимости составляет в среднем 10-15 итераций, что ненамного превышает значения, полученные при моделировании (5-6 итераций).

2. Достоверность результатов подтверждается тем, что итерационная процедура сходится с большой точностью к одному и тому же значению в широком диапазоне начальных приближений. Это исключает «случайную» сходимость алгоритма.

3. Результаты эксперимента хорошо соотносятся как с теоретическими границами погрешности многолучевости для данного дискриминатора, так и с результатами компьютерного моделирования

4. В частотном диапазоне L2 алгоритмическая поправка имеет систематическое смещении, которое в пересчете на псевдодальность составляет «7 м.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ РАБОТЫ

1. Усложнение модели многолучевости за счет ввода дополнительного луча увеличивает дисперсию оценки задержки. Этот эффект наиболее выражен в случае БФМ, и менее заметен в случае МЧМ. Для сигналов с составной структурой чипа это влияние минимально.

2. Определены потенциальные границы точности измерения задержки сигнала для стандартных и спектрально-эффективных форматов модуляции для однолучевой и двухлучевой моделей отражений. Показано, что в отсутствии многолучевой помехи для узкополосного радиоприемного тракта (полоса ограничена 1-2 лепестками спектра сигнала) при классической минимально-частотной модуляции выигрыш в точности достигает 1,8 раз.

3. Показано, что применение спектрально-эффективных форматов модуляции с составной структурой чипа даст выигрыш в точности до 11 дБ для широкополосного приемного тракта.

4. При охвате фильтром более двух лепестков спектра сигналы со спектрально-эффективной модуляцией заметно выигрывают в плане устойчивости к многолучевой помехе. Расчет огибающей многолучевости показывает, что при десяти лепестках это преимущество может достигать четырех раз.

5. Разработан алгоритмический метод нейтрализации многолучевой ошибки, основанный на множественных измерениях в четырех следящих контурах. Экспериментально доказана его работоспособность. Результаты эксперимента хорошо соотносятся как с теоретическими границами погрешности многолучевости для данного дискриминатора, так и с результатами компьютерного моделирования.

ПУБЛИКАЦИИ

публикации в изданиях, рекомендованных ВАК

1. Соколов A.A. Алгоритмическое исключение многолучевой погрешности из радионавигационных измерений [Текст] / Ипатов В.П., Соколов A.A., Сошин М.П. И Известия вузов России, Радиоэлектроника, СПб, 2010. - Вып. 6 - С. 42-50.

2. Соколов A.A. Потенциальная точность измерения запаздывания сигнала в присутствии многолучевой помехи [Текст] / Ипатов В.П., Соколов A.A., Шебшаевич Б.В. // Известия вузов России, Радиоэлектроника, СПб, 2011. - Вып. 2 -С. 18-25.

3. Соколов A.A. Потенциальная точность измерения запаздывания сигнала с минимальной частотной модуляцией в присутствии многолучевой помехи [Текст] / Ипатов В.П., Соколов A.A. // Известия вузов России, Радиоэлектроника, СПб, 2011. - Вып. 6 - С. 13-19.

в сборниках трудов международных, всероссийских и межвузовских конференций

4. Соколов A.A. Влияние многолучевой помехи на точность слежения за характерной точкой принимаемого сигнала в дискриминаторах навигационной аппаратуры ГНСС [Текст] / Соколов A.A. // Девятый международный симпозиум по электромагнитной совместимости и электромагнитной экологии , СПб, 2011. - С. 290-293

5. Соколов A.A. Алгоритмические методы борьбы с многолучевостью в навигациоиных системах космического базирования [Текст] / Соколов A.A. // Шестая международная молодежная научно-техническая конференция РТ-2010, Севастополь, 2010. - С. 144

6. Соколов A.A. Потенциальная точность измерения запаздывания сигнала с минимальной частотной модуляцией в присутствии многолучевой помехи [Текст] / Соколов A.A. // Седьмая международная молодежная научно-техническая конференция РТ-20И, Севастополь, 2011. - С. 112.

7. Соколов A.A. Алгоритмические методы нейтрализации многолучевых помех в РНС космического базирования [Текст] / Соколов A.A. // 64-я научно-техническая конференция, посвященная Дню Радио, СПб, 2009. - С.29-30.

8. Соколов A.A. Учет фазовых измерений в компенсации многолучевой погрешности в ГНСС [Текст] / Соколов A.A. // 65-я научно-техническая конференция, посвященная Дню Радио, СПб, 2010. - С. 35-36.

9. Соколов A.A. Иммунитет к многолучевой помехе сигналов со спектрально-эффективной модуляцией [Текст] / Соколов A.A. // 66-я научно-техническая конференция, посвященная Дню Радио, СПб, 2011. - С. 35-36.

10. Соколов A.A. Учет фазовых измерений в устранении многолучевой погрешности при множественных измерениях [Текст] / Соколов A.A. // Сборник тезисов докладов Школы радиоэлектроники, СПб, 2010. - С. 26

Подписано в печать 01.03.12. Формат 60*84 1/16. Бумага офсетная. Печать офсетная. Печ. л. 1,0. Тираж 100 экз. Заказ 14.

Отпечатано с готового оригинал-макета в типографии Издательства СПбГЭТУ "ЛЭТИ"

Издательство СПбГЭТУ "ЛЭТИ" 197376, С.-Петербург, ул. Проф. Попова, 5

Текст работы Соколов, Андрей Андреевич, диссертация по теме Радиолокация и радионавигация

61 12-5/3006

САНКТ-ПЕТЕРБУРГСКИЙ ГОСУДАРСТВЕННЫЙ ЭЛЕКТРОТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ «ЛЭТИ» им. В.И. Ульянова (Ленина)

На правах рукописи

СОКОЛОВ Андрей Андреевич

НЕЙТРАЛИЗАЦИЯ МНОГОЛУЧЕВЫХ ПОМЕХ В РНС КОСМИЧЕСКОГО БАЗИРОВАНИЯ

Специальность 05.12.14 «Радиолокация и радионавигация»

Диссертация на соискание ученой степени кандидата технических наук

Научный руководитель: д.т.н., проф. Ипатов В.П.

Санкт-Петербург - 2012

СОДЕРЖАНИЕ

СПИСОК ИСПОЛЬЗОВАННЫХ СОКРАЩЕНИЙ И ОБОЗНАЧЕНИЙ..........5

ВВЕДЕНИЕ...........................................................................................................11

1. ПРОБЛЕМА БОРЬБЫ С МНОГОЛУЧЕВЫМИ ПОМЕХАМИ В РАДИОНАВИГАЦИИ КОСМИЧЕСКОГО БАЗИРОВАНИЯ.........................20

1.1 Историческая справка....................................................................20

1.2 Основные принципы функционирования и структура спутниковой радионавигационной системы..................................................24

1.3 СРНС ГЛОНАСС............................................................................28

1.4 СРНС НАУ8ТА1ШР8...................................................................29

1.5 Дальномерные сигналы СРНС......................................................30

1.6 Обработка сигналов в аппаратуре потребителей СРНС.............33

1.6.1 Петля автосопровождения по времени (АПВ)......................35

1.6.2 Петля фазовой автоподстройки...............................................38

1.7 Помехи многолучевости в РНС космического базирования......38

1.7.1 Профили систематической ошибки слежения в присутствии

отражения......................................................................................................39

1.8 Современные методы борьбы с эффектами многолучевого распространения сигналов в спутниковых РНС............................................41

1.8.1 Пространственная обработка сигналов..................................42

1.8.2 Узкий коррелятор.....................................................................44

1.8.3 Технологии МЕТ и Ш.С..........................................................45

1.8.4 Технология МЕБЬЬ.................................................................47

1.8.5 Стробовые корреляторы..........................................................47

1.8.6 Алгоритмические методы борьбы с многолучевостью........51

1.9 Выводы по главе.............................................................................56

2. АЛЬТЕРНАТИВНЫЕ ВАРИАНТЫ ФОРМАТА МОДУЛЯЦИИ НАВИГАЦИОННОГО СИГНАЛА.....................................................................58

2.1 Введение..........................................................................................58

2.2 Спектр мощности МЧМ сигнала...................................................59

2.3 Составные чипы на базе МЧМ и их спектры мощности............61

2.4 Конкретизация параметров предлагаемых форматов модуляции.........................................................................................................66

2.5 Выводы по главе.............................................................................71

3. ОЦЕНКА ЗАПАЗДЫВАНИЯ СИГНАЛОВ НОВОГО ФОРМАТА В ПРИСУТСТВИИ МНОГОЛУЧЕВОЙ ПОМЕХИ..............................................73

3.1 Введение..........................................................................................73

3.2 Граница Крамера-Рао и потенциальная точность измерения запаздывания сигнала, искаженного многолучевой помехой......................74

3.2.1 Оценка запаздывания профильтрованного прямоугольного импульса........................................................................................................81

3.2.2 Оценка запаздывания профильтрованного чипа

косинусоидальной формы.............................................................................91

3.2.3 Потенциальная точность оценки запаздывания сигнала ОД в присутствии многолучевости.....................................................................100

3.2.4 Потенциальная точность оценки запаздывания сигнала СД в присутствии многолучевости.....................................................................107

3.3 Огибающие многолучевой погрешности при временном дискриминировании сигналов новых форматов..........................................110

3.3.1 Оптимальный дискриминатор для слежения за характерной точкой сигнала.............................................................................................111

3.3.2 Ранне-поздний дискриминатор.............................................120

3.4 Выводы по главе...........................................................................123

4. МНОЖЕСТВЕННЫЕ ИЗМЕРЕНИЯ И АЛГОРИТМЫ НЕЙТРАЛИЗАЦИИ МНОГОЛУЧЕВОЙ ПОГРЕШНОСТИ ОЦЕНКИ ПСЕВ ДО ДАЛЬНОСТИ......................................................................................125

4.1 Введение........................................................................................125

4.2 Множественные измерения в условиях многолучевости.........126

4.3 Алгоритм устранения влияния многолучевой погрешности измерения задержки НС в когерентной петле AHB....................................127

4.3.1 Влияние значения задержки на сходимость алгоритма......132

4.3.2 Понижение размерности задачи............................................133

4.4 Результаты моделирования..........................................................134

4.4.1 Моделирование в условиях отсутствия шума.....................135

4.4.2 Влияние ширины полосы частотно-селективного тракта на работу алгоритма.........................................................................................138

4.4.3 Шумовая ошибка алгоритма.................................................143

4.4.4 Алгоритм компенсации ошибки многолучевости по измерениям сигналов ОД и СД..................................................................146

4.5 Результаты экспериментального исследования........................148

4.5.1 Цели и задачи..........................................................................148

4.5.2 Исходные условия и методика..............................................148

4.5.3 Результаты эксперимента......................................................149

4.6 Выводы по главе...........................................................................154

ЗАКЛЮЧЕНИЕ...................................................................................................155

БИБЛИОГРАФИЧЕСКИЙ СПИСОК...............................................................157

СПИСОК ИСПОЛЬЗОВАННЫХ СОКРАЩЕНИЙ И ОБОЗНАЧЕНИЙ

АКФ - автокорреляционная функция

АПВ - автоподстройка по времени

АЦП - аналого-цифровой преобразователь

ВТ - высокой точности (канал)

ВЧ - высокочастотный

ГГС - генератор гармонического сигнала

ГДК - генератор дальномерного кода

ГЛОНАСС - глобальная навигационная спутниковая система

ГНСС - глобальная навигационная спутниковая система

ГУН - генератор управляемый напряжением

ИСЗ - искусственный спутник Земли

МВ - метка времени

МЧМ - минимальная частотная модуляция

НАЛ - навигационная аппаратура потребителей

НКА - навигационный космический аппарат

КА - космический аппарат

НП - навигационные параметры

НС - навигационный сигнал

ОД - открытый доступ

ПКУ - подсистема контроля и управления

ПНКА - подсистема навигационных космических аппаратов

ПСП - псевдослучайная последовательность

ПТ - пониженной точности (канал)

ПФ - петлевой фильтр

РНП - радионавигационные параметры

РНТ - радионавигационные точки

СД - санкционированный доступ

ссз - схема слежения за задержкой

ССФ - схема слежения за фазой

ссч - схема слежения за частотой

ст - стандартной точности (канал)

ФАГТЧ - фазовая автоподстройка частоты

ФМ - фазовая модуляция

ци - цифровая информация

вое - binary offset carrier

BPSK - binary phase shift keying

CDMA - code division multiple access

DLL - delay lock loop

FDMA - frequency division multiple access

FIMLA - fast iterative maximum likelihood algorithm

FLL - frequency lock loop

GPS - global positioning system

HRC - high resolution correlator

MEDLL - multipath estimation delay lock loop

MET - multipath elimination technology

MSK - minimum shift keying

PRN - pseudo-random noise

A - эффективная площадь приемной антенны

4 - амплитуда /'-ой компоненты навигационного сигнала

a - амплитуда многолучевой составляющей

- элементы кодовой последовательности

ВССЗ - шумовая полоса схемы слежения за задержкой

с - скорость света

c(t) - дальномерный код

А - расстояние от фазового центра антенны приемника до /-го

спутника

е(в)

¿?£(б,Д,5) Е

Ewl

fc

fxi

Л

/д> fd

К

F

GXf) G(.) h(t)

IE ' IL ' ^P

J(s,A,S)

К Ks

rM

»1,(0

Me(.)

ЛГ

- дискриминационная характеристика

- суммарная дискриминационная характеристика, искаженная многолучевым сигналом

- энергия сигнала

- энергия сигнала после фильтра

- частота грубого канала GPS

- частота /-ой несущей

- допплеровское смещение частоты навигационного сигнала, его оценка

- протяженность зоны неопределенности по оси частот

- контрольная полоса

- спектральная плотность мощности

- градиентная матрица

- импульсная характеристика стационарной линейной системы

- значения корреляционных интегралов в синфазном канале: ранний, поздний, текущий

- интеграл корреляции со стробом

- номер несущей частоты навигационного сигнала

- число видимых космических аппаратов

- коэффициент корреляции

- информационное сообщение z'-ro спутника

- медиана выборки

- общее число опытов при моделировании

- число успешных опытов

- число видимых спутников.

- спектральная плотность мощности белого шума

Nm - количество многолучевых компонент в сигнале

Nc - длина дальномерного кода

n.(t) - шум

п% - эквивалентные флюктуации запаздывания на входе дискри-

минатора

п - количество микрочипов в составном чипе

пт - частота гармонического колебания в составном чипе, отне-

сенная к fc

пр - число периодов гармонического колебания в знакоперемен-

ном чипе

пс - число знакопеременных чипов на длительности чипа С/А

кода GPS

р - пиковая мощность сигнала, вероятность сходимости итера-

ционной процедуры q2 - отношение сигнал-шум по мощности

Qe,Ql,Qp - значения корреляционных интегралов в квадратурном канале: ранний, поздний, текущий

rsd(t0) - коэффициент корреляции

S - крутизна дискриминационной характеристики

s(t) - навигационный сигнал

5.(/) - спектр сигнала

S(t) - комплексная огибающая сигнала

Si(«) - интегральный синус

s0(t) - микрочип

t - время

Т - время накопления отсчетов в корреляторе.

t0 - временное положение характерной точкой слежения за сигналом

««(О

К

K,vy,vz

var,.{•} W

Щ-)

w w

msk ' rr psk

x, x0

X, y,z

y(t)

ЯоУц» У2с>У2р a,P

P^

X

Ac

Дц

p

8s,, 8s_

- сигнал на входе приемника

- скорость потребителя

- компоненты вектора скорости

- дисперсия

- ширина полосы частотно-селективного тракта приемника

- функция правдоподобия

- полосы MSK и PSK сигналов, содержащие 99 % их энергии

- вектор параметров многолучевости и его начальное приближение

- координаты потребителя

- координаты i-го НКА

- реализация наблюдаемого случайного сигнала

- измеренные значения псевдодальности в канале ВТ и ПТ на двух частотах

- амплитуды многолучевости на частоте L1 и L2 соответственно

- доля полной мощности чипа в полуширине полосы

- расхождение шкал времени "потребитель-спутник"

- длительность чипа С/А-кода GPS

- ошибка измерения псевдодальности после алгоритмической коррекции

- среднее геометрическое

- длина чипа Р-кода GPS

- смещение дискриминационной характеристики из-за многолучевости (синфазный и противофазный случай)

- длина микрочипа

- половина длительности строба дискриминатора автопод-

стройки по времени 8 - рассогласование входного и опорного сигналов в корреля-

торе

Ф. - элемент матрицы Фишера

ф(. - начальная фаза /-ой компоненты навигационного сигнала

у. - относительные потери в точности оценки задержки

Х1 - коэффициент Лагранжа

X - длина волны

М- - смещение нуля дискриминационной характеристики

П - плотность потока мощности

0А - задержка к -го многолучевого сигнала по отношению к

прямому

рг(т) - автокорреляционная функция

Ро 1' Ро 9 — десяти- и девяностопроцентный квантиль погрешности алгоритма

су2а - средний квадрат ошибки измерений после коррекции

а2о - средний квадрат суммарной ошибки первичных измерений

с72 - дисперсия шума измерений

т0,а0, 0О - параметры многолучевой помехи

т, х - время распространения навигационного сигнала, его оценка

\|/(») - оператор Тигера-Кайзера

ВВЕДЕНИЕ

Актуальность работы. В настоящее время спутниковая навигация приобрела огромную популярность. Если раньше использование космических аппаратов в целях определения координат объекта считалось исключительно прерогативой военных, то в последние годы глобальные навигационные спутниковые системы (ГНСС) прочно вошли в повседневную жизнь сотен тысяч обычных людей по всему миру. На современном этапе гражданскому потребителю доступны две равноценные и взаимодополняющие ГНСС: GPS и ГЛОНАСС. На стадии пробных испытаний находится ГНСС Galileo, продвигаемая Европейским сообществом. В ближайшем будущем Китай планирует развернуть ГНСС двойного применения Compass. Кроме того, Япония близка к старту практического развертывания региональной навигационной спутниковой системы QZSS (Quasi-Zenith Satellite System). Подобный рост количества систем спутниковой навигации в пределах одного частотного диапазона L ведет к спектральным коллизиям как между ними самими, так и с системами в смежных спектральных диапазонах. В этом плане традиционные для ГНСС навигационные сигналы (НС) с бинарной фазовой манипуляцией, обладающие широким спектром, оказываются не лучшими с точки зрения электромагнитной совместимости и альтернативой им могут явиться спектрально-эффективные сигналы, в частности, сигналы с минимальной частотной модуляцией (МЧМ). При этом правомочен вопрос о тактических характеристиках ГНСС, использующих такие сигналы, и прежде всего о точностных показателях. Говоря конкретно, задача сводится к оценке потенциальной точности измерения радионавигационных параметров (РНП), в частности, задержки навигационного сигнала.

Проблема обеспечения качественных навигационных измерений в условиях многолучевости является одной из ключевых в вопросах функционирования ГНСС. Поэтому сформулированную выше задачу целесообразно решать именно для сценария интерференционной картины на приемной сто-

роне. При этом с практической точки зрения важным также является вопрос об устойчивости спектрально-эффективных сигналов к многолучевости в сценариях приема с помощью следящих систем.

Повсеместной практикой является моделирование многолучевости одной отраженной репликой прямого сигнала. Такое допущение редко соответствует действительности, что побуждает к усложнению модели в сторону увеличения числа лучей для корректного рассмотрения вопросов потенциальной точности оценки РНП. В то же время вопрос применимости однолу-чевой модели в разработках навигационных приложений остается открытым, поскольку использование сложных моделей отражений в реальных устройствах не всегда возможно, а адекватность упрощенных версий должна быть в достаточной степени аргументирована. В этом плане полезной может оказаться разработка алгоритмической процедуры нейтрализации многолучевости, поскольку именно алгоритмические методы наиболее чувствительны к модели сигналов, заложенной в их основе.

Названные факторы обуславливают актуальность исследования потенциальной точности оценки РНП спектрально-эффективных сигналов в условиях многолучевого распространения, анализа устойчивости таких сигналов к многолучевости и разработки алгоритмических методов ее нейтрализации.

Цели и задачи работы. Целью диссертационной работы является исследование потенциальных характеристик измерения времени в ГНСС в том числе для сигналов спектрально-эффективного формата, а также разработка методов нейтрализации многолучевых ошибок на основе множественных измерений. Для достижения поставленной цели в диссертационной работе предстояло решить следующие задачи:

1. На основе литературных источников рассмотреть методы борьбы с многолучевостью.

2. Принять для дальнейшего анализа конкретные форматы спектрально-эффективной модуляции.

3. Аналитически рассчитать точности измерения РНП при оптимальном приеме и ограниченном спектре сигнальных чипов для различных сценариев присутствия многолучевости, моделируемой одним и двумя лучами.

4. Разработать алгоритмическую процедуру нейтрализации негативных эффектов многолучевости и детально исследовать ее характеристики путем компьютерного моделирования.

5. Провести экспериментальные исследования с использованием записей реальных сигналов двухчастотных навигационных приемоизме-рителей ГНСС.

Методы исследования. Для решения поставленных в диссертационной работе задач были использованы методы теории сигналов, теории вероятностей и математической статистики, гармонического анализа, численные методы и методы математического моделирования. Экспериментальные исследования выполнены методом полунатурного моделирования с использованием записей реальных сигналов навигационных приемоизмерителей в ГНСС.

Научная новизна работы. Новыми являются следующие результаты:

1. Обобщены и систематизированы результаты расчетов потенциальных границ точности оценки задержки НС для однолучевой модели отражений.

2. Получены аналитические выражения для расчета дисперсии оценки РНП при оптимальном приеме НС для двухлучевой модели отражений. Указанные равенства получены для ограниченного в частотно-селективном тракте приемника чипа стандартных и спектрально-эффективных сигналов.

3. Установлено влияние второго отраженного луча на потенциальную точность оценки РНП как для стандартных, так и для спектрально-эффективных сигналов.

4. С целью оценки устойчивости к многолучевости, рассчитаны огибающие многолучевой ошибки для спектрально-эффективных сигналов в сценариях измерения РНП с помощью классического ранне-позднего дискриминатора и оптимального дискриминатора слежения за характерной точкой. Даны рекомендации по выбору параметров дискриминатора.

5. Разработан алгоритмический метод нейтрализации многолучевости, основанный на множественных наблюдениях НС в четырех каналах навигатора - стандартной и высокой точности в двух частотных диапазонах - и о дно лучевой модели отражений.

6. Экспериментально подтверждена работоспо�