Перов А.И., Харисов В.Н. ГЛОНАСС. Принципы построения и функционирования (4-е издание, 2010) (1151865), страница 86
Текст из файла (страница 86)
Определение вектора состояния НС в орбитальной прямоугольной системе координат ОХ1Х2, лежащей в орбитальной плоскости с началом в центре Земли, ось Х1 направлена вдоль фокальной оси к перигею, ось Х2 — по нормали к фокальной оси: Соотношения для составляющих скорости НС в орбитальной системе координат получаются дифференцированием координат ~х,'т,хте) по времени Сс учетом (3.17)): ии'Б-е соа~д +' ) х„'Р— 1 — есор Е~ ' ) паял Е1 ") «оР 1 — есор~ Е1 ' ) ) 5. Пересчет ортов е,'",е'„р)1 орбитальной прямоугольной системы координат ОХ~Х~ в инерциальную систему координат ПЗ-90, осуществляется путем трех последовательных поворотов орбитальной системы координат на углы е„*, 1, Я* (см. (3.7)).
Соотношения для проекций орта в,'Р на оси инерциальной системы координат могут быть получены также из (3.18) при г =1, 3 = 0 и замене долготы восходящего узла Й в инерциальной системе на соответствующую долготу восходящего узла Л* в системе координат ПЗ-90: в„=сокау„созЯ вЂ” япв„япА соя~', о о е, =соим„яви +яви„совЛ сов~, ед, =свау„яп~', е,~ — — -япв„созЛ вЂ” совы„япЯ соз1, в,2 = совы„яп1. о е,, = — япв„япЛ +соьв„сояЯ сои, 522 а для проекций орта вД аналогичные соотношения получаются из (3.18) при 9=к/2: Аппаратура потребителей 6. Преобразование вектора состояния НС из орбитальной системы в систему координат ПЗ-90. Введем Х =)Х у д), е1е = ее, ее1 е„(, ет е— - еет еет е~т) Тогда вектор координат НС в системе координат ПЗ-90 определяется соотношением Х=х1 е, +х2 е2.
ор 0 ор 0 Преобразование вектора скорости НС из орбитальной системы координат во вращающуюся систему координат ПЗ-90 проводится в два этапа. Сначала вектор скорости преобразуется в неподвижную систему координат, ось которой смещена относительно оси ХО инерциальной системы координат ОХОУОД на долготу восходящего узла, равную Я, в соответствии с формулой ХУ = ХоРЕО + ХоРЕО 1 1 2 2" т На втором шаге вектор скорости Х' = Х' У' 2' из неподвижной систе- мы координат пересчитывается во вращающуюся систему координат ПЗ-90 в соответствии с общим правилом (3.21) по формулам Х = Х' + в~У, У = У' — ауз Х, У = У' 13.8.3.
Расчет вектора состояния НС на основе оперативной информации Расчеты в НАП проводятся с целью уточнения параметров движения НС на момент времени г, по данным эфемерид (которые содержатся в разделе оперативной информации навигационного сообщения), даются на моменты времени ~„и обновляются через каждые 15 мин, поэтому для времени ~, выполняется условие ~~~ — ~,.~ <15 мин. Процедура пересчета проводится численным интегрированием дифференциальных уравнений орбитального движения (3.24).
Начальными условиями для интегрирования системы уравнений являются данные эфемерид. Интегрирование проводится классическим одношаговым методом Рунге — Кутта четвертого порядка ~4.8), когда при разложении нелинейной функции в ряд удерживаются все члены порядка малости до Л, где Л вЂ” шаг численного интег- 4 рирования. Для этого метода локальная ошибка имеет порядок Л .
5 При использовании этого метода для решения дифференциального уравнения первого порядка вида а1у/а1х = ~(х,у) при начальном условии у(хО) = уО расчеты проводятся по правилу Уо+1 = Уп + КО + 2К1 + 2К2 + К, /6, 523 Глава 13 1 1 где Ко = /(х у.) Л' К~ = /' х, + — Л.У,+ Кв~ Л 2 2 / К =Х~х + — /~,у + — К,~ /~; К, =Х~х + — Л,у + — К ~ Л. 13.8.4. Пересчет координат потребителя из земной в геодезическую систему координат 3) анализируется значение.0: если Р=О,тоВ=г~/(2~4), А=О, Н=~яп — а 1 — е яп В, 2 2 если Р>0, то А„=агсяп(у/Р),приэтом: то Л=г„, то Е=л — А„, если х>0, у>0, если х<О,У>0, если х<О,У<0, если х>0, у<0, г= +л„, то г =2 — Л„; 4) анализируется значение ~: если: =О, то В =О, Н =,0 — а, в других случаях находятся вспомогательные величины г, с, р: =,(к~ <-у~.~~~, с =агск~п(г!г), р = ее а/2~', реализуется итеративный процесс: 524 Алгоритм оценки навигационных параметров (13.32) формирует оценки вектора потребителя в геоцентрической системе координат ОХИ, связанной с Землей.
Потребителя во многих случаях интересуют свои координаты в геодезической системе координат (см. п. 3.1). Поэтому в АП необходимо осуществить пересчет координат из геоцентрической системы в геодезическую. Общие формулы связи двух систем координат имеют вид (3.3), (3.4), а пересчет по этим выражениям может быть реализован следующим вычислительным алгоритмом [13.31: 1) определяются исходные данные — координаты х, у, ~ в геоцентрической системе координат ОХИ, выходные данные — координаты В, Е, Н в геодезической системе координат; 2) вычисляется вспомогательная величина .0= х +у Аппаратура потребителей к,=О; Ь=сек,, к =ассака~ркка(2Ь2кк~1 — сакка Ь), если модуль разности ~з2 — з, ~ < е, где е — требуемая точность вычислений, то В=Ь, В=22сскВскккпВ-сБ — е ккп В, в противном случае полагается з, = з, и вычисления повторяются, начиная с расчета Ь. Во всех случаях выбирается е = 0,0001". После нахождения В,Е,Н становится известной матрица перехода из геоцентрической системы координат в прямоугольную горизонтальную, что позволяет осуществить преобразование составляющих скорости движения подвижного объекта.
13.9. Использование оптимальных алгоритмов обработки сигналов и информации в аппаратуре потребителей Как пример использования оптимальных алгоритмов обработки сигналов и информации в современной НАП СРНС, описанных в гл. 6, рассмотрим образцы приемной аппаратуры, выпускаемые Зачад Мач1дапоп Кукиш Со. (1.ехоп, Магог и др.) и доступные на российском рынке. Эти приемники воплощают в себе наиболее передовые технологии в области приема и обработки спутниковых сигналов ГЛОНАСС и ОРИ.
Перечислим основные особенности реализованных технических решений. Использование ССФ 3-го порядка для слежения за фазой принимаемых сигналов, что обеспечивает устойчивое сопровождение спутников при установке аппаратуры на высокоманевренных объектах. Использование поддержки схемы слежения за задержкой сигнала (ССЗ) от схемы ССФ [см. п.6.3.6.5, 6.3.6.6], что позволяет уменьшить ширину шумовой полосы ССЗ до единиц и даже долей герца и, таким образом, уменьшить флуктуационные ошибки оценки псевдо дальности до уровня единиц-десятков сантиметров.
Соттоп (Со-Ор) Тгасйпд [О]. Представляет собой несколько упрощенный вариант одноэтапного слежения за параметрами сигнала [см. п. 6.5, 17.2.21. Суть алгоритма заключается в разделении динамики фазы сигнала на две составляющие: быстро- (обусловленную динамикой объекта и нестабильностью опорного генератора приемника) и медленноменяющуюся (обусловленную движением спутника и нестабильностью его опорного генератора). Такое разделение позволяет использовать для слежения за медленноменяющейся составляющей ССФ с шириной шумовой полосы порядка 2 Гц, а для слежения за быстроменяющейся составляющей — с полосой порядка 20 Гц. Количество узкополосных ССФ равно количеству сопровождаемых спутников, широкополос- 525 Глава 13 ная ССФ вЂ” одна для всех спутников, и использует для оценки параметров суммарную мощность принимаемых сигналов.
Такой подход позволяет повысить чувствительность приемника примерно на 10 дБ по сравнению с традиционно используемым подходом, основанном на использовании для слежения за каждым из сигналов НС автономных ССФ. Подавление многолучевоети осуществляется на основе двух подходов, которые могут использоваться как совместно, так и по отдельности. Первый подход состоит в использовании специально сконструированной антенны (13.6], обеспечивающей подавление сигналов, отраженных от земной поверхности и объектов, расположенных ниже антенны. Однако, данный подход не обеспечивает подавление сигналов, отраженных от объектов, расположенных над антенной (высокие здания, деревья и т.п.), поэтому в приемнике осуществляется дополнительная сигнальная обработка, позволяющая в существенной степени снизить ошибки оценки параметров сигнала, обусловленных эффектом многолучевости.
Дифференциальный режим навигаиии. Приемник в реальном времени обеспечивает декодирование корректирующих поправок, передаваемых в соответствии с форматом КТСМ 104 (см. гл. 12), и оценивание координат объекта с точностью до десятков сантиметров. Режим КТК (Яеа1 Тле Ктета1к) ~13.71 является дальнейшим расширением режима дифференциальной навигации и позволяет на основе приема данных об измерениях псевдо дальности и псевдо фазы, выполняемых в базовом приемнике, в реальном времени оценивать текущие координаты приемника с точностью до единиц миллиметров в статическом режиме и до единиц сантиметров в динамическом режиме. 1.1 й 1.2.
Приемник может осуществлять слежение за сигналами спутников в двух частотных диапазонах П и Е2, модулированными кодами стандартной точности С1А, кодами высокой точности Р (режим АК выключен) и криптозащищенными кодами высокой точности У (режим АЯ включен). В последнем случае используется специальная технология Кеппсоде1езв 1гасЫпд. Использование измерений в двух частотных диапазонах позволяет выполнить компенсацию погрешностей, обусловленных задержками сигналов в ионосфере, и тем самым повысить точность оценки координат в автономном режиме навигации.