В последние годы развиваются глобальные спутниковые радионавигационные системы (СРНС) второго поколения, предназначенные в первую очередь для определения координат и скорости подвижных объектов в любой точке планеты и воздушном приземном пространстве. В настоящее время в мире развернуты 2 такие системы: американская система GPS (Global Positioning System — глобальная система местоопределения) и российская спутниковая система ГЛОНАСС (глобальная навигационная спутниковая система) [1,2,3].
Европейскими странами также планируется разработка подобной системы, предназначенной в первую очередь для решения коммерческих задач, в то время как ГЛОНАСС и GPS до недавнего времени были предназначены исключительно для военного использования. Все системы, как действующие, так и планируемые имеют сходные принципы построения, с некоторыми отличиями в структуре пространственного сегмента и навигационных сигналов.
При разработке и начальном использовании этих систем предполагалось, что с их помощью будут решаться задачи определения координат, составляющих вектора скорости и точного времени на борту подвижного объекта.
Между тем существует широкий круг пользователей, которым наряду со знанием координат и составляющих скорости необходима информация о пространственной ориентации. К таким объектам относятся летательные аппараты, в первую очередь вертолеты, для которых пространственная ориентация не всегда совпадает с направлением их движения. Для судов морского и речного флота знание ориентации позволит вычислить угол сноса и своевременно внести коррекцию курса, что позволит с меньшими затратами времени и топлива прийти в заданную маршрутную точку. Кроме того, существуют объекты, для которых знание угловой ориентации является жизненно важным, например спутники, которым необходимо с высокой точностью ориентировать приемопередающие антенны на Землю, а солнечные батареи — на Солнце [4]. Информация об ориентации позволяет повысить точность определения координат исполнительных механизмов некоторых объектов, что является необходимым для ряда областей применения, например при добыче полезных ископаемых [5, 6, 51, 54, 59].
В настоящее время известны два радиотехнических метода измерения направлений: амплитудный и фазовый (интерферометрический).
Амплитудный метод предполагает использование антенной системы с очень узкой, в общем случае веретенообразной диаграммой направленности (ДН), установление слежения по направлению за источником радиосигнала, размещенном на навигационном космическом аппарате (НКА) и измерении углов между осями антенны и осями объекта [7]. Точность определения угловой ориентации зависит от размера апертуры антенны. Как показывают расчеты, для того, чтобы на рабочих частотах СРНС ГЛОНАСС и GPS (1200 -1600 Мгц) обеспечить точность определения угловых направлений порядка 10 угловых минут необходима зеркальная антенна или фазированная антенная решетка с раскрывом 3×3 м [8]. Очевидно, что стоимость таких систем высока и имеются серьезные проблемы с их размещением на подвижных объектах. Однако имеется опыт создания и эксплуатации подобных систем (радиосекстантов, работающих по излучению Солнца) для морских объектов.
Интерферометрический метод определения направлений состоит в том, что несколько разнесенных в пространстве антенн принимают сигнал от одного источника. Расстояние между антеннами (база) примерно равна размеру апертуры антенны, обеспечивающей ту же точность при амплитудном методе. Измерительное устройство оценивает разность хода сигнала до антенн.
Эквивалентом разности хода является разность фаз принятых сигналов. По приведенным в [2, 7] оценкам, при расстоянии между антеннами В-2м, длине волны 1=20 см (длина волны НКА системы ГЛОНАСС) для достижения точности определения ориентации объекта в пределах одной угловой минуты нужно обеспечить погрешность измерения разности фаз в пределах 1°, что является технически достижимым.
Анализ публикаций в технической литературе и рекламных материалов показывает, что вся выпускаемая навигационная аппаратура потребителей (НАЛ) СРНС, выполняющая определение ориентации объектов использует интерферометрический метод измерений. Применение этого метода требует наличия на борту объекта нескольких пространственно разнесенных антенн, принимающих сигналы НКА. На основании измерений фазовых сдвигов сигналов НКА, принятых разнесенными антеннами, пользователь оценивает ориентацию базовой линии, соединяющей приемные антенны НАЛ, а по нейсвою собственную ориентацию.
С целью повышения точности определения ориентации объектов используют НАП определения ориентации с расстоянием между антеннами, достигающим нескольких метров. Существующие фазометры имеют диапазон однозначных измерений в пределах одной длины волны, что порождает проблему оценки количества целых циклов неоднозначности в измеряемых фазовых сдвигах (ФС) для правильной оценки ориентации. Кроме того, вследствие неидентичности фазовых характеристик высокочастотных блоков НАП по каждому из антенных каналов измеренные фазовые сдвиги имеют систематическую погрешность. Также существенное влияние на точность определения ориентации объектов оказывает случайная составляющая погрешности фазовых измерений.
В настоящее время разработаны различные методы разрешения неоднозначности и определения ориентации объектов. Из анализа публикаций были выделены несколько методов определения ориентации объектов при наличии неоднозначностей в измеряемых фазовых сдвигах.
Первый метод не требует прямой оценки вектора неоднозначностей, и основан на переходе от исходной системы уравнений для полных ФС к системе уравнений, составленных по разностям ФС в различные моменты времени. Разности ФС могут накапливаться при неподвижном объекте, за счет движения спутников, либо, если ориентация объекта изменяется, к приращениям фаз за счет движения спутников добавляются приращения, вызванные поворотами объекта. Динамический метод, основанный на изменении ориентации объекта, подробно рассмотрен в [12, 13]. Следует отметить, что для низко динамичных объектов, использующих НАП с большим расстоянием между антеннами, реализация динамического метода затруднена, а подчас и просто невозможна.
Другим методом является применение многоантенных пеленгаторов, с различными значениями расстояния между антеннами. Минимальное расстояние между антеннами пеленгатора составляет около половины длины волны, при этом в измеренном ФС неоднозначность отсутствует. Определение ориентации пеленгатора по отношению к источнику излучения выполняется исходя из принципа максимального правдоподобия [14−16]. Использование данного метода требует наличия дополнительных антенн и каналов обработки их сигналов, что удорожает НАП определения ориентации и усложняет алгоритмы ее работы.
Существуют методы определения ориентации объектов и разрешения фазовой неоднозначности за счет применения дополнительных датчиков пространственной ориентации, например инерциальных. При этом информация, поступающая с датчика, используются для грубой оценки ориентации объекта, а затем, после разрешения фазовой неоднозначности, результат уточняется при помощи НАП определения ориентации. Однако существующие инерциальные датчики ориентации имеют высокую стоимость, значительную потребляемую мощность и габариты. Перечисленные недостатки делают невозможным применение рассмотренного метода на некоторых объектах.
В качестве дополнительных сведений при разрешении неоднозначностей могут быть также использованы измерения ФС сигналов СРНС, полученные на частоте L2, если для определения угловой ориентации объекта используется двухчастотная НАЛ. Располагая измерениями ФС принятых сигналов на частотах Lj и L2, и используя оценку разности фаз принятых сигналов на разностной частоте L1-L2 можно более эффективно решать задачу разрешения неоднозначности фазовых измерений [17]. Данные способы разрешения фазовой неоднозначности находят применение в прецизионных геодезических приемниках GPS, использующих при оценке координат измерения фаз несущих принимаемых сигналов [18]. Недостатком рассмотренного метода является необходимость усложнения НАЛ с целью приема сигналов НКА на частоте Ь2, что приводит к увеличению стоимости аппаратуры.
В зарубежной и отечественной литературе [19, 20, 60] описаны алгоритмы аналитического разрешения неоднозначностей путем перебора всех возможных значений комбинаций целых циклов ФС для нескольких спутников. При этом количество перебираемых решений может исчисляться тысячами, а иногда и миллионами. Из всех возможных комбинаций выбирается только одна, удовлетворяющая некоторым критериям и априорным сведениям о возможных положениях вектора-базы и расстоянии между антеннами. Так в [19] описывается способ определения ориентации, основанный на переборе потенциальных значений решений первоначально по 2-м спутникам. При расстоянии между антеннами, составляющем 1 м и длине волны — 20 см максимальное количество потенциальных решений составляет 188. При привлечении измерений от избыточных спутников число решений может быть уменьшено до одного — верного решения путем нахождения той их комбинации, при которой величина расхождения ориентации объекта, найденной по опорным спутникам и с привлечением дополнительных спутников минимальна. При этом для исключения ложных решений на более ранних этапах поиска используется пороговая величина ошибки оценки ориентации, зависящая от случайной составляющей погрешности измерений ФС принимаемых сигналов. Величина порога должна обеспечивать высокую вероятность отбраковки ложных решений, при низкой вероятности пропуска верных решений.
Фирмой Trimble Navigation в [20] запатентован способ прямого поиска вектора неоднозначностей отличающийся от предыдущего тем, что поиск оптимальной комбинации ведется по всему набору измеренных ФС. Алгоритм оптимизируется с целью исключения ложных комбинаций на возможно более ранних этапах поиска.
Модификации алгоритмов прямого поиска целочисленной неоднозначности исследовались автором и опубликованы в [32, 33]. В результате проведенных исследований выяснилось, что при среднеквадратическом значении погрешности фазовых измерений стф < 3° при расстоянии между антеннами НАЛ равным 1 м вероятность правильного разрешения неоднозначности превышает 95%, при расположении спутников с геометрическим фактором менее 3. Погрешность определения ориентации составляет 5−10 угловых минут в зависимости от геометрического расположения НКА, а время разрешения неоднозначности при использовании опорного созвездия из 3-х НКА (1331 возможных решений) составило около 1с при использовании вычислителя INTEL 486 DX4 — 120 МГц.
Однако для реализации данных алгоритмов перед проведением измерений необходимо произвести калибровку высокочастотных каналов НАЛ с целью исключения аппаратурной составляющей систематической погрешности фазовых измерений, вызванной неидентичностью высокочастотных каналов по сигналам каждой из антенн. Время калибровки при отсутствии априорной информации о взаимном положении антенн, образующих антенную систему и длине фидеров, соединяющих антенны с приемным блоком, может достигать нескольких часов. При этом необходимым условием проведения такой калибровки является неподвижность антенной системы объекта [21]. Кроме того, данные алгоритмы устойчиво работают при расстоянии между антеннами НАЛ не более одного метра. При увеличении расстояния между антеннами алгоритмы прямого поиска неоднозначности дают довольно высокую вероятность ошибки и требуют значительного времени.
Образцы аппаратуры, выполняющей определение ориентации объектов, выпускаются только зарубежными фирмами Sercel (Франция), Trimble Navigatuion, Ashtech (США) и работают по сигналам системы GPS.
В связи с этим важной и актуальной является задача разработки высокоточных алгоритмов определения ориентации, позволяющих выполнять оценку углового положения как по сигналам систем ГЛОНАСС и GPS, так и в совмещенном режиме ГЛОНАСС/GPS при любом расстоянии между антеннами объекта.
Цель работы — разработка высокоточных методов определения ориентации объектов по сигналам спутниковых систем ГЛОНАСС и GPS, некритичных к величине расстояния между антеннами НАП.
Задачи исследования.
В настоящей работе решаются следующие задачи:
1) Разработка алгоритмов определения ориентации по приращениям фазовых сдвигов принимаемых сигналов, вызванных движением НКА, при неподвижной антенной системе объекта.
2) Разработка алгоритмов определения ориентации объектов по измеренным значениям фазовых сдвигов принимаемых сигналов НКА.
3) Исследование погрешностей алгоритмов определения ориентации объектов.
4) Экспериментальные исследования спектрально-корреляционных и статистических характеристик погрешностей измерения фазовых сдвигов сигналов НКА СРНС ГЛОНАСС и GPS.
5) Экспериментальные исследования разработанных алгоритмов определения ориентации при различных расстояниях между антеннами объекта.
Научная новизна.
1) Разработаны алгоритмы определения ориентации неподвижного объекта по приращениям фазовых сдвигов принятых сигналов. Данные алгоритмы обеспечивают уменьшение погрешности определения расстояния между антеннами, азимута и угла места в несколько раз по сравнению с известным способом, основанным на вычислении разностей между первым (опорным) и текущим значениями фазового сдвига. Реализации трех предложенных алгоритмов защищены патентом РФ.
2) Впервые разработаны алгоритмы определения ориентации объекта по измеренным значениям фазовых сдвигов, не требующие разрешения фазовой неоднозначности и обеспечивающие повышение точности определения ориентации в несколько раз. Алгоритмы позволяют определять ориентацию при увеличении расстояния между антеннами до пяти и более метров. Реализация одного из предложенных алгоритмов защищена патентом РФ.
3) Разработаны методы аналитической оценки погрешностей ориентации объекта, позволяющие оценивать точность определения ориентации при неравноточных и взаимно коррелированных измерениях фазовых сдвигов сигналов НКА.
4) Экспериментально исследованы спектрально-корреляционные и статистические характеристики погрешностей измерения фазовых сдвигов сигналов НКА, в результате которых показана корреляция погрешностей измерений по сигналам разных НКА СРНС ГЛОНАСС и GPS, а также коррелированность погрешностей измерений во времени.
5) Выполнено статистическое моделирование и экспериментальные исследования разработанных алгоритмов на образцах НАЛ типа МРК-11 и МРК-15у для различных расстояний между антеннами при работе по сигналам СРНС ГЛОНАСС, GPS и в совмещенном режиме ГЛОНАСС/GPS.
Достоверность результатов. Достоверность результатов подтверждается корректным использованием математического аппарата, результатами моделирования и натурных испытаний на образцах НАЛ типа МРК-11 и МРК-15у.
Защищаемые научные положения:
1) Алгоритмы определения ориентации неподвижного объекта, основанные на измерении приращений всех фазовых сдвигов принимаемых сигналов НКА, которые позволяют повысить точность оценки ориентации в несколько раз.
2) Алгоритмы определения ориентации объекта по измеренным значениям фазовых сдвигов, не требуют разрешения фазовой неоднозначности и позволяют определять ориентацию при увеличении расстояния между антеннами до нескольких метров.
3) Аналитический и статистический методы оценки погрешностей определения ориентации объекта и расстояния между антеннами для разработанных алгоритмов.
4) В результате экспериментальных исследований спектральных, корреляционных и статистических характеристик погрешностей измерений фазовых сдвигов сигналов НКА СРНС ГЛОНАСС и GPS доказана коррелированность погрешностей измерений сигналов разных НКА, а также временная коррелированность погрешностей.
5) Результаты экспериментальных исследований разработанных алгоритмов подтверждают возможность определения ориентации объектов для расстояния между антеннами до пяти и более метров, при работе по сигналам СРНС ГЛОНАСС, GPS и в совмещенном режиме ГЛОНАСС/GPS.
Практическая ценность. Разработаны алгоритмы, позволяющие определять ориентацию объекта с погрешностью 0.1° и менее. Разработаны алгоритмы оценки погрешностей определения ориентации в условиях коррелированных и неравноточных измерений фазовых сдвигов сигналов НКА СРНС. Проведено моделирование и экспериментальные исследования разработанных алгоритмов на образцах НАП МРК-11 и МРК-15у, работающих по сигналам СРНС ГЛОНАСС и GPS. Выполнены экспериментальные исследования статистических, спектральных и корреляционных характеристик погрешностей измерений фазовых сдвигов сигналов НКА СРНС.
Методы исследования. В диссертационной работе использованы методы математического анализа, линейной алгебры, статистической радиотехники, математического и статистического моделирования.
Реализация и внедрение. В диссертацию вошли результаты НИР и ОКР, проводимых в соответствие с решением ГК ВПВ № 2505 от 05.05.91 г., госбюджетных и хоздоговорных работ КГТУ:
Разработка бортовой навигационной аппаратуры потребителя космических радионавигационных систем ГЛОНАСС/ОРЭ для ракет-носителей, разгонных блоков и космических аппаратов;
Исследование методов измерения углового положения объектов на основе глобальных навигационных спутниковых систем по программе «Университеты России» ;
Исследование и разработка аппаратно-программных средств с расширенными функциональными возможностями на основе глобальных навигационных спутниковых систем ГЛОНАСС/ОРЭ;
Разработка бортовой навигационной аппаратуры для космического аппарата «Гонец» ;
Разработка специализированной аппаратуры для ориентации станций связи.
Апробация работы. Результаты работы докладывались и обсуждались на научно-технической конференции «Спутниковые системы связи и навигации» в г. Красноярске (1997г.), на Всероссийской научно-практической конференции студентов, аспирантов и молодых специалистов «Решетневские чтения» в г. Красноярске (1997г.), на региональной научно-технической конференции студентов и молодых ученых Сибири и Дальнего Востока, посвященной 102-й годовщине Дня Радио в г. Красноярске (1997 г.), на Всероссийской с международным участием научно-технической конференции молодых ученых и студентов, посвященной 103-й годовщине Дня Радио в г. Красноярске (1998).
Публикации. Основные результаты диссертации опубликованы в 12 печатных работах, защищены патентом РФ.
Структура и объем диссертации
Диссертация состоит из введения, четырех глав и четырех приложений. Общий объем работы.
Основные результаты, полученные при выполнении работы, можно сформулировать следующим образом:
1. Разработаны алгоритмы определения ориентации неподвижных объектов по приращениям фазовых сдвигов сигналов НКА СРНС ГЛОНАСС и GPS. Рассмотрены варианты определения углового положения по пяти комбинациям разностей фазовых сдвигов при различном их расположении внутри измерительного интервала, обеспечивающие уменьшение погрешности оценки углового положения в несколько раз по сравнению с известным методом.
2. Разработаны алгоритмы определения ориентации объекта по измеренным значениям фазовых сдвигов сигналов НКА, не требующие разрешения фазовой неоднозначности и обеспечивающие повышение точности определения ориентации в несколько раз. Алгоритмы позволяют определять ориентацию при увеличении расстояния между антеннами до пяти и более метров.
3. Разработан метод аналитической оценки погрешностей определения ориентации объекта. Выполнена оценка погрешностей определения ориентации объекта методами аналитической оценки и статистического моделирования.
4. Разработано программное обеспечение для определения углового положения по результатам измерений, выполняемых НАП типа МРК-11 и МРК-15у, проведена отладка программного обеспечения и экспериментальные исследования.
5. Результаты моделирования и экспериментальных исследований подтвердили возможность измерения углового положения неподвижных объектов с погрешностью 0.1° и менее для расстояния между антеннами до пяти и более метров при работе по сигналам НКА СРНС ГЛОНАСС, GPS и в совмещенном режиме ГЛОНАСС/GPS.
Автор выражает благодарность коллективу НИИ радиотехники КГТУ за помощь, оказанную при выполнении данной работы.
ЗАКЛЮЧЕНИЕ
.
В результате проведенных автором исследований, а также научно-исследовательских и опытно-конструкторских работ, проводившихся при его непосредственном участии, были решены задачи по разработке алгоритмов определения ориентации неподвижных объектов, не требующие разрешения фазовой неоднозначности, работающие при любых значениях расстояния между антеннами объекта.