способ посадки летательных аппаратов с использованием спутниковой навигационной системы и система посадки на его основе

область техники

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

известны способы посадки JIA, позволяющие определять координаты JIA при посадке с использованием спутниковой навигационной системы, в которых прием сигналов спутников осуществляется на земле локальной контрольно-корректирующей станцией (JIKKC) с известными координатами, вычисляющей дифференциальные поправки к определяемым на борту JIA по сигналам спутников значениям псевдодальности и, соответственно, значениям координат JIA [1,2].

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

в то же время о необходимости учета ошибок многолучевости и на JIKKC, и на борту JIA указывается в ряде международных документов [3,4], однако в них не предлагается конкретного способа определения ошибки многолучевости на JIKKC, а учет ошибок многолучевости на борту предлагается проводить по способу, в котором ошибка многолучевости заведомо значительно выше реальной. в этих документах также указывается на необходимость оценки целостности системы посадки - достоверности (малой вероятности потери необходимой точности) определения координат JIA. для контроля достоверности необходимо иметь

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

известен ряд источников информации [5-7], в которых предлагаются различные способы определения ошибки многолучевости, фиксируемой наземным приемником лккс, где вычисляется зеркальная компонента наземной ошибки многолучевости, возникающая вследствие зеркального отражения сигналов спутников от рельефа и окружающих объектов.

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

в способе определения ошибки многолучевости [5] используется корреляция эффектов многолучевости в измерениях псевдодальности и отношения сигнал/шум наземного приемника, при этом ошибки многолучевости оцениваются по текущим измерениям. однако применить этот способ в системе посадки JIA невозможно в связи с тем, что он предназначен для определения составляющих ошибки многолучевости с интервалом корреляции единицы-десятки секунд, а в соответствии с требованиями авиационных стандартов для аппаратуры посадки JIA [3] необходимо определение составляющих ошибки многолучевости с интервалом корреляции более 100 секунд.

наиболее результативным для этой цели является способ [6,7], в котором ошибка многолучевости также формируется по текущим измерениям псевдодальности и отношения сигнал/шум наземного приемника. кроме того, в [6] предлагается методика оценки достоверности определения координат JIA. однако этот способ не дает точных результатов и может приводить к заниженным значениям ошибки многолучевости, так как формируемые в этом способе с помощью полосовых фильтров разности кодовых и фазовых измерений псевдодальностей включают в себя, помимо ошибки многолучевости, удвоенную ионосферную задержку радиосигнала, которая меняется во времени и неэффективно фильтруется, искажая оценку ошибки многолучевости.

необходимо отметить, что все рассмотренные способы [5-7] неприменимы для оценки ошибок многолучевости на борту из-за быстрого изменения угловой ориентации спутников при маневрировании JlA

в качестве ближайшего аналога заявляемого способа выбран способ посадки, описанный в патенте [8], который является наиболее распространенным и реализуется в международной системе посадки [3].

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

в качестве ближайшего аналога заявляемой системы выбрана система посадки, описанная в патенте [8], с помощью которой реализуется описанный в указанном патенте [8] способ посадки.

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

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

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

раскрытие изобретения

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

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

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

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

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

в частном случае выполнения изобретения используют M наземных приемников, для каждого из которых формируют диаграмму объемного

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

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

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

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

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

осуществляют прием сигналов спутников бортовым приемником и определение псевдодальностей на борту летательного аппарата, осуществляют прием сигналов спутников наземным приемником на локальной контрольно-корректирующей станции с известным

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

при этом согласно заявляемому способу производят формирование диаграммы объемного распределения прогнозируемой дисперсии наземной ошибки многолучевости (HOM) путем предварительного измерения отраженных сигналов спутников и введение ее в наземную базу данных, по которой непрерывно определяют значение прогнозируемое дисперсии HOM для каждого спутника в соответствии с его угловыми координатами. далее в процессе текущего приема сигналов спутников определяют значение зеркальной компоненты многолучевости и текущей дисперсии HOM одним из известных способов, например [6,7]. затем производят сравнение прогнозируемой и текущей дисперсий HOM и выбирают максимальное из них . в процессе текущего приема сигналов спутников также определяют дисперсию шумовой ошибки наземного приемника известным способом, например [9]. далее вычисляют значение дисперсии наземной ошибки псевдодальности (ноп), используя дисперсию шумовой ошибки и максимальное значение из сравниваемых дисперсий текущей и прогнозируемой ошибок многолучевости. вычисляют для каждого спутника среднеквадратическое значение дисперсии ноп и транслируют его через лпд на борт ла. одновременно выполняется коррекция в наземной базе данных диаграммы прогнозируемой дисперсии HOM по текущему значению дисперсии HOM .

на борту выполняют предварительное формирование диаграммы объемного распределения прогнозируемой дисперсии бортовой ошибки многолучевости (бом) для заданного типа ла и вводят ее в бортовую базу данных. определяют прогнозируемую дисперсию бом с учетом угловой

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

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

прогнозируемое значение дисперсии & mp- _b з \^) для I -того спутника

( Z ⁼ l...iV ). для повышения надежности посадки JlA диаграмма объемного распределения дисперсии ошибки формируется при установке

лккс.

вычисление мгновенного значения зеркальной компоненты HOM и текущей дисперсии HOM необходимо в случае быстрых изменений помеховой обстановки (например, при пролете других JIA или перемещении крупных объектов в районе лккс). методы измерения

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

^л 2 г -\ сравнение и выбор ^mp_з у ' наибольшего из значений прогнозируемой и текущей дисперсий HOM в соответствии с уравнением

^l _P _ ₃ 0) = max{σ^_ ₃ (0,σ^(/)} ₍₁₎ позволяет обеспечить наибольшую безопасность при посадке JIA.

расчет дисперсии шумовой ошибки выполняется с использованием значения сигнал/шум, измеряемого наземным приемником. документ [3] регламентирует эталонный алгоритм фильтрации псев до дал ьности : y(k) = y(k-l) + δ _у2 (k) + α-(y,(k) - y(k-l) - δy ₂ (k)), y(0), (2) где y(k) - результат фильтрации; к - дискретное время; yi(k) - отсчеты псевдо дальности в кодовом канале; δy ₂ (k) = yг(k) - yг(k-l) — приращение фазовых измерений y ₂ (k); α = δt/T - весовой коэффициент; δt = 0,5 c - интервал временной дискретизации, соответствующий периоду выдачи дифференциальных данных; T = IOO c - постоянная времени; y(0) - начальные условия.

алгоритм фильтрации (2) основан на принципе инвариантности к динамике псевдодальности и, по отношению к кодовым измерениям псевдо дальности yi(k), является фильтром нижних частот с постоянной времени T = 100 с.

оценка дисперсии шумовой ошибки борта <$ б_ «ш( V'. >*) / определяется как оценка дисперсии флюктуационной составляющей ошибки на выходе данного фильтра [9]: о)

где ^i \ ^J ^J и ^2 \^ j - оценки дисперсий шумовых ошибок в отсчетах псевдо дальности в кодовом и фазовом каналах; α = 5-1 (г ³ - весовой коэффициент фильтра, описываемого выражением (2); k - дискретное время.

величины σ ^f 2( Vα ' ) _/ , _> входящие в (3), определяются по известной методике [10]:

где χi и χг - коэффициенты энергетических потерь в кодовом и фазовом каналах приемника спутниковых сигналов соответственно; с - скорость света; то - длительность селекторного импульса; TD - период дискретизации; Bc и Bp - эффективные полосы систем слежения за задержкой кода и фазой несущей, соответственно; fц - несущая частота Ll; - отношение сигнал/шум в полосе частот 1 гц; k - дискретное время.

^л 2 / -\

после определения дисперсии HOM G _mp _ ₃ \!) и дисперсии

шумовой ошибки вычисляют дисперсию наземной ошибки псевдодальности : σlr_ ₃ ⁽ о = ^ά lp_з ⁽ о + ^σ i_ ₃ ⁽ о . (6)

среднеквадратическое значение дисперсии ноп pr_з \ ) транслируют на борт ла (по методике [3, 4]).

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

компоненты ном необходима из-за изменения отражающих свойств наземных предметов при перемене погоды и времен года. при коррекции используется известный алгоритм а -фильтра [5]: , (7) где а - весовой коэффициент, характеризующий желаемую скорость обновления наземной диаграммы; I - номер спутника; к - дискретное время обновления базы данных.

формирование объемной диаграммы распределения дисперсии бортовой ошибки многолучевости (бом) выполняется путем теоретического или экспериментального исследования диаграммы рассеяния корпуса выбранного типа JIA. методы экспериментального исследования отражений сигналов от элементов конструкции JIA известны и описаны, например, в [11]. в предложенном способе результаты данных исследований отражений сигналов впервые используются в навигационном процессе для формирования объемной диаграммы прогнозируемого значения дисперсии бом.

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

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

^ mp_б V^» / для / -того спутника в момент времени k.

определение дисперсии шумовой ошибки ^ w б V ' _/ бортового приемника подобно определению шумовой ошибки наземного приемника .

определение дисперсии боп I -того спутника выполняется путем сложения дисперсий σ _б ² (i ₉ k) = σ _m ² _{p б} (i,k) + σl_ _б (i,k) _{. (8)}

затем выполняется определение дисперсии суммарной ошибки системы

^fоf v/ путем сложения дисперсий измерения JIKKC G _{pr з} w , борта

о б V* ^: » O и дополнительных тропосферных ^тр и ионосферных

погрешностей ^ ^" _м , рекомендуемых стандартом [3],

<,(i) = <__(0 + σ* ² (U) + σl + σl _p ⁽⁹⁾

данные о погрешностях измерений используются для вычисления уточненных координат и достоверности их определения . надежность посадки JTA достигается благодаря контролю достоверности (целостности) навигационной системы, который реализуется на борту JIA. результат контроля используется для предупреждения экипажа в случае, если ошибки измерения координат превышают размеры «тoннeля безопасности)), в котором должен находиться JIA при посадке. для выработки сигнала предупреждения на борту согласно [3] выполняется расчет дисперсий ошибок измерения координат и затем защитных уровней по высоте V и боковому отклонению L от заданной глиссады. предупреждение вырабатывается, если защитные уровни превышают пороги Vp и Lp, определяемые нормами безопасности. в бортовом вычислительном устройстве для решения навигационной задачи и задачи контроля целостности используется прямоугольная правосторонняя система координат, начало о которой совмещено с началом взлетно-посадочной полосы; ось Ox направлена вдоль взлетно-посадочной

полосы, касательно к поверхности референсного эллипсоида; ось Oz перпендикулярна указанной поверхности и направлена вне эллипсоида; ось Oy дополняет систему координат до правосторонней.

линеаризованная модель наблюдения, используемая для решения навигационной задачи, имеет следующий вид [3]:

Ap = H - δX + e, где δр - вектор отклонений измеренных псевдодальностей относительно расчетных псевдодальностей, определенных для заданной опорной точки X ₀ ; H - матрица направляющих косинусов в заданной системе координат; δх - искомый вектор приращения координат относительно опорной точки X ₀ ; е - вектор невязок.

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

где S=[H ¹ W ^-1 H] ^-1 H ¹ W ^"1 - проекционная матрица;

W — ^Q шg ^ ⁽ 7 _t0t (Ij,..., ⁽ T _ш (/V )J _ диагональная весовая матрица; S _jjп -

элемент матрицы S, расположенный в i-й строке и п-м столбце; v/ - угол

глиссады; "с ) - коэффициент, обеспечивающий требуемое значение риска потери достоверности Pc;

1 °°

Q(x) =- _f = fexp( -f72 ) λ

благодаря предложенному способу уточняются значения дисперсий

общей ошибки по формуле (9), которая используется для формирования весовой и проекционной матриц при определении

уточненных координат и вычислении достоверности их определения в соответствии с защитными уровнями по формулам (10) и (1 1).

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

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

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

для повышения надежности работы JlKKC целесообразно использовать несколько приемных каналов, разнесенных в пространстве [3]. дифференциальные поправки JlKKC в этом случае формируются как среднеарифметическое значение поправок, выработанных каждым приемником.

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

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

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

равна G ^" mp_з \?> ) , а после вычисления дисперсии шумовой ошибки наземного приемника получают дисперсию шумовой составляющей

ошибки измерения псевдодальности ^ _шj V) ^) . далее вычисляют дисперсию ноп для к -го приемника:

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

результирующее значение ноп транслируют на борт JlA. заявляемая система посадки (см. фигуру), реализующая заявленный способ, содержит N спутников 1 с известными координатами (Ci,... C _N ), наземную аппаратуру (а), включающую локальную контрольно- корректирующую станцию с по крайней мере одним приемником 2, входом сопряженного с наземной антенной спутниковых сигналов 3, выход которого соединен с первым входом наземного вычислителя дифференциальных поправок 4, на второй вход которого поступает информация о координатах C _A места локальной контрольно-

корректирующей станции, выход которого соединен с первым входом передатчика линии передачи данных 5, второй вход которого соединен с выходом наземной базы данных 6, бортовую аппаратуру (в), включающую бортовой приемник 14, входом сопряженный с бортовой антенной спутниковых сигналов 15, выход которого соединен с первым входом бортового вычислителя текущих координат 16, второй вход которого соединен с первым выходом приемника линии передачи данных 17, второй выход которого соединен с первым входом вычислителя отклонения координат от заданной глиссады 18, второй вход которого соединен с выходом бортового вычислителя текущих координат 16, а выход которого является информационным выходом 8 отклонения от заданной глиссады.

в наземную аппаратуру (а) введены вычислитель дисперсии текущей наземной ошибки многолучевости 7, вычислитель дисперсии шумовой ошибки наземного приемника 8, формирователь диаграммы объемного распределения прогнозируемой дисперсии наземной ошибки многолучевости 9, база данных прогнозируемой дисперсии наземной ошибки многолучевости 10, компаратор 11, определитель максимального значения дисперсии наземной ошибки многолучевости 12, вычислитель среднеквадратического значения дисперсии наземной ошибки псевдо дальности 13, а в бортовую аппаратуру (в) введены база данных объемного распределения прогнозируемой дисперсии бортовой ошибки многолучевости 19, определитель ориентации летательного аппарата 20, вычислитель дисперсии бортовой ошибки многолучевости 21, вычислитель дисперсии шумовой ошибки бортового приемника 22, вычислитель дисперсии бортовой ошибки псевдодальности 23, вычислитель дисперсии суммарной ошибоки псевдодальности системы посадки 24, бортовая база допустимой погрешности определения координат 25 и вычислитель достоверности определения координат 26. при этом реализованы новые структурные связи в системе. первый и второй дополнительные выходы наземного приемника 2 соединены соответственно с входом вычислителя дисперсии текущей наземной ошибки многолучевости 7 и входом вычислителя дисперсии

шумовой ошибки наземного приемника 8, выход формирователя диаграммы объемного распределения прогнозируемой дисперсии наземной ошибки многолучевости 9 соединен с входом базы данных диаграммы прогнозируемой наземной ошибки многолучевости 10, выход которой соединен с первым входом компаратора 11 , второй вход которого соединен с выходом вычислителя дисперсии текущей наземной ошибки многолучевости 7, а выход компаратора 11 соединен с входом определителя максимального значения наземной ошибки многолучевости 12, выход которого соединен с первым входом вычислителя среднеквадратического значения наземной ошибки псевдо дальности 13, второй вход которого соединен с выходом вычислителя дисперсии шумовой ошибки наземного приемника 8, а выход которого соединен с третьим входом передатчика линии передачи данных 5. выход вычислителя дисперсии текущей наземной ошибки многолучевости 7 дополнительно соединен с корректирующим входом базы данных прогнозируемой дисперсии наземной ошибки многолучевости 10.

выход базы данных объемного распределения прогнозируемой дисперсии бортовой ошибки многолучевости 19 соединен с первым входом вычислителя дисперсии бортовой ошибки многолучевости 21, второй и третий входы которого соединены соответственно с выходом определителя угловой ориентации летательного аппарата 20 и первым дополнительным выходом бортового приемника 14, второй дополнительный выход которого соединен с входом вычислителя шумовой ошибки бортового приемника 22, выход которого соединен с первым входом вычислителя дисперсии бортовой ошибки псевдодальности 23, второй вход которого соединен с выходом вычислителя дисперсии бортовой ошибки многолучевости 21, а выход которого соединен с первым входом вычислителя дисперсии суммарной ошибки псевдодальности 24, второй вход которого соединен с вторым дополнительным выходом приемника линии передачи данных 17, первый выход 24 соединен с дополнительным входом вычислителя текущих координат 16, второй выход сумматора ошибок псевдодальности 24 соединен с первым входом вычислителя достоверности определения

координат 26, второй вход которого соединен с выходом бортовой базы допустимой погрешности определения координат 25, а выход Tp вычислителя достоверности определения координат 26 является информационным выходом о превышении допустимой погрешности определения координат в системе посадки (сигнала "тревога").

на фигуре также показаны отражающие объекты 27 в районе размещения JIKKC и переотражающие элементы конструкции JIA 28. система посадки работает следующим образом. в наземной аппаратуре (а) в вычислителе дифференциальных поправок 4 формируются дифференциальные поправки и через передатчик лпд 5 передаются на борт вместе с координатами, поступающими из наземной базы данных координат заданной глиссады 6. вычисление текущих координат JIA с учетом принятых в приемнике лпд 17 дифференциальных поправок осуществляется вычислителем 16. отклонение л а от выбранной глиссады (по вертикали и горизонтали) определяется вычислителем 18, на входы которого поступают текущие координаты ла из вычислителя 16 и заданные координаты из приемника лпд 17.

объемная диаграмма дисперсий прогнозируемой HOM формируется в формирователе 9 путем предварительного исследования ошибок многолучевости известными методами [6,7], запоминается в базе данных 10 и используется по сигналу текущего времени τ. одновременно с использованием отношения сигнал/шум и выходного сигнала приемника 2 в вычислителе 7 определяется текущее значение дисперсии HOM. сравнение дисперсий выполняется компаратором 11, а выбор наибольшей - определителем максимального значения 12. определение текущей дисперсии HOM в вычислителе 7 позволяет повысить надежность измерения ошибок при отражениях от движущейся наземной техники и пролетающих ла. определение среднеквадратического значения HOM выполняется в вычислителе 13, на входы которого поступают дисперсии шумовой ошибки и наземной ошибки многолучевости. техника измерения дисперсии шумовой ошибки рассмотрена выше.

в бортовой аппаратуре (в) база данных дисперсии прогнозируемой ошибки многолучевости 19 формируется путем исследования отражений для выбранного типа JlA [11]. вычислитель дисперсии бортовой ошибки многолучевости 21 использует данные об угловой ориентации JIA в земной системе координат, поступающие из определителя угловой ориентации JIA 20, и угловых координатах спутников из приемника 14. вычислитель дисперсии шумовой ошибки бортового приемника 22 использует измерения отношения сигнал/шум, получаемые из приемника 14. дисперсия бортовой ошибки псевдодальности определяется в вычислителе 23 путем суммирования дисперсии бом и дисперсии шумовой ошибки бортового приемника. определение дисперсии суммарной ошибки системы посадки выполняется в вычислителе 24, в который поступают значения дисперсии наземной и бортовой ошибок псевдо дальности. дисперсии суммарных ошибок системы поступают из 24 в вычислитель координат 16, где используются при формировании весовой матрицы метода наименьших квадратов [3] для получения уточненных координат и в вычислитель достоверности координат JIA 26, который вырабатывает предупреждение в случае превышения допустимой погрешности определения координат JIA, хранящейся в бортовой базе данных 25. из вычислителя 16 уточненные значения координат поступают в вычислитель 18, где используются для определения уточненных значений отклонений от глиссады.

таким образом, введение формирователей и баз объемных диаграмм прогнозируемых дисперсий наземной и бортовой ошибок многолучевости и их использование в системе посадки в совокупности с блоками, определяющими и обрабатывающими дифференциальные поправки псевдодальностей, и блоками, вычисляющими дисперсии ошибок многолучевого распространения и шумовых ошибок, позволяет существенно уточнить определяемые координаты JIA и достоверность их определения кроме того, как указывалось выше, значительно увеличена надежность системы посадки при отказе текущих измерений HOM благодаря каналу прогнозирования дисперсии и использованию

максимального значения HOM при сравнении прогнозируемой и текущей дисперсий.

предложенная система посадки реализована с использованием приемоизмерителей фирмы «Javad Nаvigаtiоп System»: в качестве бортового приемника - модуль JNSlOO, в качестве наземного приемника - модуль

JGG 12T. также использовались аппаратно-программные модули, построенные на базе широко распространенных стандартных устройств аналоговой и цифровой вычислительной техники фирм «ANALOG

DEVICES», «AMPRO», «AлEKCAHдEP элEKTPиK», «TEXAS INSTRUMENTS», «MAXIM», «ALTERA».

для разработки программного обеспечения, реализующего необходимые функции упомянутых устройств, использовались стандартные языки программирования («C», «C++»), программно-математическое обеспечений фирм «MICROSOFT», «BORLAND», «QNX» и известные формулы геодезических преобразований.

промышленная применимость

заявляемые способ и система посадки летательных аппаратов являются перспективными для использования на аэродромах, необорудованных инструментальными системами посадки. математическое моделирование, полунатурные и летные испытания системы посадки, которые проводились на самолетах як-42 и ил-76 показали, что использование предложенного способа позволяет повысить точность определения координат на 50%, а достоверность их определения в 2-3 раза.

источники информации:

1. сетевые спутниковые радионавигационные системы [текст] / п.п.дмитриев [и др.]; отв. ред. в.с.шебшаевич. - 2-е изд. - M.: радио и связь, 1993. - 408 с: ил. 2. патент сша JN°5 361 212, кл. G01S5/00, заявл. 11.02.92 г., опубл.

01.11.1994 г.

3. DO-245A, Minimum аviаtiоп Sуstеm реrfогmапсе Stапdагds fог Lосаl аrеа аugmепtаtiоп Sуstеm (LAAS) [еlесtrопiс rеsоurсе] // Rаdiо тесhпiсаl соmmissiоп fог аегопаutiсs, Dесеmbеr 9, 2004. - режим доступа: http://www.rtca.org/doclist.asp. 4. DO-246C, GNSS ваsеd ргесisiоп арргоасh Lосаl аrеа аugmепtаtiоп

Sуstеm (LAAS) - Sigпаl-iп-Sрасе Iпtеrfасе сопtrоl Dосumепt (ICD) [еlесtrопiс rеsоurсе] // Rаdiо тесhпiсаl соmmissiоп fог аегопаutiсs, арril 07, 2005. - режим доступа: http://www.rtca.org/doclist.asp.

5. патент рф JN° 2 237 256, кл. G01S5/00, H04B1/06, заявл. 21.02.2001 г., опубл. 27.09.2004 г.

6. соколов а. а. оценка ошибок дифференциальных поправок JIKKC авиационного назначения. автореферат диссертации на соискание ученой степени кандидата технических наук. с-пб., 2005г.

7. Slееwаеgеп, J. мultiраth мitigаtiоп, вепеfϊts frоm usiпg thе Sigпаl-tо- Nоisе Rаtiо [техt] / J. Slееwаеgеп // юN GPS-1997: ргосееdiпgs оf thе юth

Iпtеrпаtiопаl тесhпiсаl мееtiпg оf thе Sаtеllitе Divisiоп оf thе Iпstitutе оf Nаvigаtiоп. - капsаs сitу, 1997. - P. 531-541.

8. патент японии JNbз 379 958 B2 8503069 а, кл. G01S5/02 заявл. 02.11.1993 г., опубл. 24.02.2003 г. 9. казаринов, ю.м. проектирование устройств фильтрации радиосигналов [текст] / ю.м.казаринов, а.и.соколов, ю.с.юрченко; под общ. ред. ю.м.казаринова. - л.: лэти, 1985. - 160 с: ил.

10. аshjаее, J. аutоmаtеd Dеsigп оf Nаvigаtiоп Rесеivеrs [техt] / J.аshjаее, м.Vоrоbiеv, а.Zhdапоv, м.Zhоdzishskу // юN GPS-1999: ргосееdiпgs оf thе 12th Iпtеrпаtiопаl тесhпiсаl мееtiпg оf thе Sаtеllitе Divisiоп оf thе Iпstitutе оf Nаvigаtiоп. - Nаshvillе, 1999. - P. 821-827. 11. Stеiпgаss, а. а нigh Rеsоlutiоп моdеl fог thе аегопаutiсаl мultiраth Nаvigаtiоп сhаппеl [Text]/A. Stеiпgаss, а.Lеhпеr, F.регеz-Fопtап, е.кubistа, м.маrtiп, в.аrbеssег-Rаstburg // юN GPS-2004: //P. 357-365.