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

УДК 528.516

А. М. Алешечкин

ПЕРЕБОРНЫЙ МЕТОД РАЗРЕШЕНИЯ НЕОДНОЗНАЧНОСТИ В МНОГОШКАЛЬНОЙ ФАЗОВОЙ РАДИОНАВИГАЦИОННОЙ СИСТЕМЕ*

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

Ключевые слова: радионавигационные системы, разрешение неоднозначности, вероятность.

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

Примером РНС наземного базирования, использующей фазовые методы определения РНП, является РНС «Крабик», разработка и серийный выпуск которой с 2001 г. осуществляется предприятием ФГУП «НІ III «Радиосвязь» (Красноярск) совместно с Сибирским федеральным университетом.

В состав РНС «Крабик» входят 3.. .6 опорных станций (ОС), размещенных на берегу в точках с известными ко -ординатами. РНС обеспечивает определение координат бортовых станций (БС) в дальномерном, разностно-даль-номерном и комбинированном режимах определения места [1].

В настоящее время частотный план РНС «Крабик» задан следующими значениями несущих частот: /0 = 421 МГ ц -основная несущая частота; /1 = 421,01 МГц, /2 = 421,1 МГц, /3 = 422 МГц, /4 = 426 МГц, /5 = 431 МГц -вспомогательные несущие частоты. Значения фазовых сдвигов (ФС) измеряются ортогональным фазометром с ограниченными сигналами [2], выполняющим оценку ФС принимаемых сигналов по отношению к их копиям, сформированным синтезатором частоты принимающей станции. Измеренные значения ФС оказываются пропорциональными измеряемым значениям радионавигационных параметров: дальностей БС-ОС - в дальномерном режиме работы, разностям дальностей - в разностно-дально-мерном режиме, разностям и суммам дальностей - в комбинированном режиме.

С целью уменьшения систематических погрешностей измерений отсчет значений ФС осуществляется на разностных (метрических) частотах РНС, получаемых путем вычитания основной несущей частоты /0 из значений вспомогательных частот /х.../5. В результате этого сетка метрических частот образует следующий ряд значений:

= /\ - /0 = 10 кГц,

= /2 - /о = 100 кГц,

^ = /з - /о = 1 МГц,

= /4 - /о = 5 МГц,

= /5 - /0 = 10 МГц.

Определение РНП осуществляется на самой высокой

метрической частоте

значение которой составляет

10 МГц. Диапазон однозначных измерений РНП для сигнала частоты Рть соответствует длине волны данного сигнала, которая составляет около з0 м. Процедура нахождения целого числа длин волн, укладывающихся в от-счетном значении РНП, называется разрешением неоднозначности (PH) фазовых измерений.

Для Р Н может быть применен известный метод пересчета измерений (МПИ), описанный в [3]. Определение значения полного фазового сдвига на самой высокой

метрической частоте

осуществляется путем после-

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

Грубые ступени РН характеризуются большими значениями длины волны 1 и, следовательно, большими значениями периода однозначных измерений. Согласно МПИ, переход от грубых к точным ступеням РН осуществляется с использованием соотношений

ф, = м, + ф,

N.. =

.ф,-ф, + 0,5

где і = 1,..., т - номер текущей шкалы; т - общее число шкал, для РНС «Крабик» равное 5; [.] означает операцию выделения целой части числа; 1, - длина волны сигнала на текущей шкале; ф, - значение измеренного ФС на ,-й шкале, выраженное в фазовых циклах; , -1 - номер предыдущей шкалы с длиной волны 1,; Ф,-1 - значение полного ФС, полученного на предыдущей шкале; Мі - значение целочисленной неоднозначности для ,-й шкалы; Ф, - значение полного ФС на текущей шкале.

*Работа выполнена при финансовой поддержке Российского фонда фундаментальных исследований (код проекта 08-08-00849).

Длины волн в (3) определяются по значениям метри-

ческих частот Рщ...Рть как

1, =

¥_

-т

N =

Зт,

-т

• ф,-1 -Ф, + 0,5

ненты, называемого квадратом суммарной невязки. При независимости [10] и равноточности измерений ФС суммарная невязка Q(х, у) может быть определена по формуле

где , = 1, ...,5 - номер метрической частоты РНС «Кра-бик»; V - значение скорости распространения радиосигналов.

С учетом (4) значения неоднозначностей (3) могут быть определены следующим образом:

Использование формулы (3) или (5) совместно с (2) позволяет реализовать процедуру РН с целью нахождения полного ФС Ф5 для самой высокой метрической частоты РНС «Крабик» Рщ. Значение полного ФС Ф5 впоследствии может использоваться для вычисления измеренных значений РНП.

Одной из серьезных проблем при разрешении фазовой неоднозначности в РНС «Крабик» является наличие многозначных отсчетов на самой низкой метрической частоте , равной 10 кГц, длина волны которой составляет 30 км. В результате этого отсчеты РНП могут быть получены в пределах 0.. .30 км, в то время как для рабочей зоны РНС, равной 150 км, значение двойной дальности может достигать 300 км. Для получения однозначного отсчета требуется задание априорных значений РНП с погрешностью ±15 км. Если же в качестве априорных данных вводятся координаты БС, то для определения места РНС в дальномерном режиме предельная погрешность не должна превышать ±7,5 км. Однако выполнение указанных требований не всегда осуществимо на практике.

Для получения однозначного отсчета РНП без необходимости задания априорных координат БС и изменения частотного плана РНС может быть использован переборный метод определения целого числа циклов первой метрической частоты 10 кГц для измеренного значения РНП.

Переборные методы находят применение при определении угловой ориентации по сигналам спутниковых радионавигационных систем [4-6], однако использование данных методов для решения задач определения координат объектов по сигналам многошкальных РНС в литературе освещено недостаточно.

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

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

При нормальном законе распределения погрешностей измерения ФС [8; 9] поиск максимума функции правдоподобия сводится к минимизации показателя ее экспо-

где , = 1,...,п - текущий номер принимающей ОС; п -общее число ОС; !1 - длина волны первой метрической частоты 0,01 МГц; к - перебираемые неоднозначности ФС на первой метрической частоте для ,-й ОС; Ф1, - значения ФС на первой метрической частоте; х1, у, - известные координаты ,-й ОС; х, у - координаты БС.

При использовании результатов измерений ФС на всех метрических частотах значение суммарной невязки Q(х, у) определяется следующим образом:

-7

Q( х, у)=а|Е £ [—Ф

где ] = 1,...,5 - номер фазовой дорожки, заданной метрической частотой Рт^; к,- целочисленная неоднозначность ФС ,-й станции на]-й метрической частоте; 1 -длина волны сигнала]-й метрической частоты; ф, - измеренное значение ФС для ,-й ОС на]-й метрической частоте.

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

При заданных значениях координат БС, приведенных в нижней строке таблицы, дальности БС-ОС1 (Кх) и БС-ОС2 (Я2) составляют соответственно 50 000,56 и 54 000,81 м. Значения суммарной невязки Q(х, у) вычислялись по выражению (7) в диапазоне дальностей Яг и К2 от 2 до 102 км с шагом 1 км, т. е. значения суммарной невязки Q(х, у) были вычислены в 104 точках. При проведении расчетов дальности отбирались и по условиям их реализуемости, а именно: сумма полученных дальностей Яг + К2 должна быть не меньше расстояния между ОС1 и ОС2. По полученным значениям дальностей определялись координаты БС для проверяемой точки и значения дальностей для ОС3-ОС6.

Полученная расстановка ОС, истинное положение БС и положения БС для расчета суммарной невязки приведены на рис. 1, а.

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

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

временных затрат (число перебираемых неоднозначностей составляет 10 для частоты 10 МГц в рабочей зоне

0.. .150 км). В связи с этим целочисленные неоднозначности кр на метрической частоте 10 кГц могут быть определены путем перебора, а на остальных метрических частотах - по МПИ независимо для каждого из РНП. В этом случае для получения отсчета потребуется анализ 121 решения, что существенно сокращает время вычисления.

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

При независимости измеренных значений радионавигационных параметров вероятность правильного определения места БС находится в соответствии с выражением

РЕ = р-Р2-.., Рп,

где р, Р2,..., Р п - вероятности правильного РН при определении значений радионавигационных параметров; п - общее число РНП, использованных для определения места БС.

С целью упрощения дальнейшего анализа можно положить измерения фазовых сдвигов по сигналам разных ОС равноточными. В этом случае выражение (8) может быть представлено в виде

РЕ = Рп,

поскольку Р = Р2 = ... = Р п = Р.

Вероятность правильного РН при оценке одного РНП определяется в соответствии с выражением [10]:

ЛС 0,5 0,5 0,5 0,5

Р =

(2р)2

-0,5 - 0,5 - 0,5 - 0,5

-0,5-

С-Д)йДхйД2йДзйД4,

*я г/ ** г/ *Иг/

4# & ЫН*

по*

Х,м

‘б( У), м

ОС6 \ -

ОС!

ОС4

ОС5 ОС3

: яи1 -

И01-

Я\, м

:> и/ и Г  »:> !'.ь и' и ||‘ <и' *: V

КИР

ИО1-

1б( X, У),

К\, м

в г

Рис. 1. Результаты расчета суммарной невязки в зависимости от координат БС: а - расстановка БС и ОС при проведении расчетов; б - суммарная невязка, полученная по сигналам двух ОС; в - суммарная невязка, полученная по сигналам трех ОС; г - суммарная невязка, полученная по сигналам шести ОС

Координаты объектов, использованные при моделировании

У. м

б

а

м

Станция Координата X, м Координата 7, м

ОС1 6 780 852 409 077

ОС2 6 843 114 387 945

ОС3 6 750 000 370 000

ОС4 6 760 000 280 000

ОС5 6 740 000 310 000

ОС6 6 800 000 270 000

БС 6 796 062 361 446

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

При взаимной независимости и равноточности измерений ФС на разных метрических частотах по сигналам разных опорных станций матрица Кд имеет вид [4]:

Кд =

ф к+1)

-к3 -оф 0 0

-к3 -0ф оф-

к4 -0ф 0

0

к4 -оф

оф •

оф

-к.

к3 = Рт1

0

0

к5 - °ф ^ф •

к4 =

- коэффициенты сопряжения соседних

где

к5 = Рт:

шкал РНС; ^т2, Ет^, ¥т^ - значения метрических

частот РНС, заданных в соответствии с выражением (1); оф - среднеквадратическое отклонение погрешности измерения ФС, принятое одинаковым для всех метрических частот.

Результаты оценки вероятности правильного разрешения неоднозначности при определении места БС по

сигналам РНС «Крабик-БМ» в зависимости от среднеквадратического отклонения погрешности измерений ФС оф для случаев работы по сигналам трех-шести ОС приведены ниже (рис. 2).

На графиках рис. 2 цифрами 1.3 обозначены следующие варианты РН:

- 1 - при определении места БС при независимой оценке РНП по МПИ (при этом неоднозначность по дорожке 10 кГц считается устраненной). В этом случае вероятность определения места БС рассчитывается в соответствии с выражением (9);

- 2 - при определении места БС во всей рабочей зоне путем перебора неоднозначностей для дорожки 10 кГц, РН для остальных метрических частот осуществляется по МПИ, выбор оптимального решения - по критерию (6);

- 3 - то же, что и график 2, только выбор оптимального решения осуществляется по критерию (7).

Значения оценок вероятности правильного разрешения неоднозначности для графиков 2 и 3 получены методом статистического моделирования по результатам обработки 1 000 отсчетов как статистических частот Р*, определяемых по формуле

Р* = т / п,

дР

Л к Ж Ж *'- , 2

—

1

3\*

Р

0.8

0.6

0.4

0.2

0 1.8 3.6 5.4 7.2 9 10.8 12.6 14.4 16.2 18

0 1.8 3.6 5.4 7.2 9 10.8 12.6 14.4 16.2 18

в г

Рис. 2. Вероятность правильного разрешения неоднозначности при работе по вариантам 1-3 при работе по сигналам:

а - трех ОС; б - четырех ОС; в - пяти ОС; г - шести ОС

1

1

где п = 1 000 - общее число испытаний; т - число испытаний, в которых разрешение неоднозначности было выполнено правильно.

Полученные результаты свидетельствуют о том, что использование при вычислении суммарной невязки Q(х, у) измерений ФС на всех частотах обеспечивает повышение вероятности правильного разрешения неоднозначности при малом числе ОС (четырех и особенно трех ОС). При пяти и шести ОС существенной разницы между способами вычисления невязки не наблюдается, что объясняется влиянием ошибок в разрешении неоднозначности при оценке РНП по МПИ при переходе с частоты 10 кГц на частоту 10 МГц. Использование расчета (7) по всем частотам обеспечивает повышение вероятности вследствие использования дополнительных измерений на других метрических частотах. При работе с большим числом ОС избыточность по результатам измерений имеется уже на одной метрической частоте, и при больших значениях среднеквадратической погрешности измерения ФС возникают ошибки разрешения неоднозначности по МПИ, которые влияют на результирующую вероятность правильного разрешения неоднозначности. Вероятность правильного разрешения неоднозначности по МПИ является верхней границей вероятности правильного разрешения неоднозначности для разработанного метода, основанного на переборе метрических частот на самой низкой метрической частоте 10 кГц.

Таким образом, полученные результаты подтверждают возможность правильного РН и однозначной оценки координат объектов в рабочей зоне РНС «Крабик». При этом не нужно задавать априорное место БС с погрешностью, не превышающей ±7,5 км, что требовалось в случае разрешения неоднозначности по МПИ. Необходимым условием реализации алгоритма является наличие избыточных ОС, используемых при определении координат БС. При этом вероятность правильного разрешения неоднозначности при работе по сигналам пятишести ОС приближается к верхней границе, задаваемой вероятностью правильного разрешения неоднозначности при независимой оценке РНП по методу пересчета измерений.

Разработанный переборный метод разрешения неоднозначности может быть использован для определения координат объектов по сигналам морских фазовых РНС с

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

Библиографические ссылки

1. Высокоточная радионавигационная система для морских потребителей / А. М. Алешечкин, В. И. Кокорин, П. Н. Иванов и др. // Г ироскопия и навигация. 2004. № 2. С. 5-12.

2. Пат. 2048676. Рос. Федерация. Устройство для измерения фазы радиосигнала / Гумеров Р. Р. [и др.]. Опубл. 1995, Бюл. № 32.

3. Разрешение неоднозначности в информационноизмерительных многошкальных приборах и системах /

В. А. Пономарев, А. В. Пономарев, Т. М. Пономарева и др. СПб. : Изд-во Воен. инж.-космич. ун-та им. А. Ф. Можайского, 2001.

4. Pat. 4963889 USA. Method and Apparatus for Precision Attitude Determination and Kinematic Positioning / Hatch R. R. Publ. 16.10.90.

5. Pat. 5296861 USA. Method and Apparatus for Maximum Likelihood Estimation Direct Integer Search in Differential Carrier Phase Attitude Determination Systems / Knight D. T. Publ. 22.03.94.

6. Пат. 2379700 Рос. Федерация. Способ угловой ориентации объекта по сигналам спутниковых радионавигационных систем / Алешечкин А. М., Кокорин В. И., Фатеев Ю. Л. Опубл. 20.01.2010, Бюл. № 2.

7. Радионавигационные системы сверхдлинноволно-вого диапазона / С. Б. Болошин, Г. А. Семенов, А. С. Гузман и др. М. : Радио и связь, 1985.

8. Денисов В. П., Дубинин Д. В. Фазовые радиопеленгаторы / Том. ун-т систем упр. и радиоэлектроники. Томск, 2002.

9. Алешечкин А. М. Повышение достоверности оценок радионавигационных параметров в радионавигационных системах с фазовыми датчиками // Датчики и системы. 2009. № 7. С. 25-29.

10. Алешечкин А. М. Вероятность правильного устранения неоднозначности в фазовой радионавигационной системе «Крабик» // Гироскопия и навигация. 2009. № 3. С. 74-82.

A. M. Aleshechkin

SEARCH METHOD OF AMBIGUITY RESOLUTION IN MULTI-SCALE PHASE RADIO-NAVIGATION SYSTEMS

Is examined the algorithm of the single-valued determination of objects coordinates from the results ofphase shifts measurements at five metric frequencies of the radio-navigation system (RNS) “Krabik”, intended for marine users navigation support. In this case single-valued coordinates realizes due to the use of a integer cycles search method of the phase ambiguity resolution and selection of a coordinates, which correspond to the maximum of a likelihood function. It is shown that the solution of the problem presented is achieved at the presence of redundancy in the system of equations, intended for determining the objects coordinates. The calculation the ambiguity resolution success rate, depending on a phase shifts root mean square error for “Krabik” RNS servicing zone is carried out.

Keywords: radio navigation, ambiguity resolution, probability.

©AnemeHKUHA. M., 2010