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

Изобретение относится к области авиационно-космического приборостроения и может найти применение в системах спутниковой навигации и геодезии. Технический результат - расширение функциональных возможностей. Для этого предложена соответствующая система (фиг.3) и последовательность операций для работы по двум орбитальным группировкам (фиг.4). Кроме того, введен специализированный вычислитель, в котором в системе навигационных уравнений учитывается пятый временной параметр, отвечающий за расхождения часов навигационного приемника от системного времени второй орбитальной группировки. При этом выполняется переход к определенному линейному векторному пространству, в котором определяется решение системы навигационных уравнений в параметрическом виде относительно двух неизвестных, с последующим определением данных неизвестных путем вычисления корней полиномиального уравнения. 2 н. и 1 з.п. ф-лы, 6 ил., 4 табл.

 

Настоящее изобретение относится к области спутниковой навигации и геодезии.

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

Геометрическое расстояние между объектом и j-ым спутником равно:

D j = ( x j x ) 2 + ( y j y ) 2 + ( z j z ) 2 ( 1 )

где {х, y, z} - искомые координаты объекта в момент приема сигнала,

{xj, yj, zj] - координаты j-го спутника в момент излучения сигнала.

Приемник на объекте измеряет время прохождения сигнала с момента его излучения до момента приема, которое, умноженное на скорость света c, равняется так называемой «псевдодальности» ρj. Ввиду того, что сигнал распространяется в атмосфере, а не в вакууме, то измеряемое время включает в себя атмосферные задержки, погрешности, связанные с переотражением сигнала, и др. Кроме того, метка времени, транслируемая спутником в излучаемом сигнале, оказывается смещенной относительно системного времени ГНСС из-за нестабильности частоты бортового генератора. Еще большее отличие псевдодальности ρj от величины Dj связано с уходом часов навигационного приемника относительно системного времени орбитальной группировки.

Обозначим момент излучения (emission) сигнала j-тым спутником t j e , а время приема (reception) приемником - tr. Тогда выражение для результата измерения псевдодальности, который не был обработан или скорректирован (далее просто псевдодальность) можно представить в виде:

ρ ˜ j = c ( t r t j e ) + c ( δ t r δ t j e ) + c T а т м + ε j r ,

или с учетом (1),

ρ ˜ j = ( x j x ) 2 + ( y j y ) 2 + ( z j z ) 2 + c ( δ t r δ t j e ) + c T j а т м + ε j r , ( 2 )

где δtr - отклонение временной шкалы приемника от шкалы системного времени ГНСС в момент приема сигнала приемником,

δ t j e - отклонение временной шкалы j-го спутника от шкалы системного времени ГНСС в момент излучения сигнала,

T j а т м - задержка, возникающая при распространении сигнала j-го спутника в атмосфере,

ε j r - инструментальная погрешность оценки псевдодальности до j-го спутника в приемнике в сумме с некомпенсированными ошибками в момент приема сигнала.

Используя эфемеридную информацию и соответствующие математические модели, можно частично скомпенсировать и учесть в выражении (2) отклонение δ t j e , задержку T j а т м и записать уравнение для скорректированной псевдодальностиу j-го спутника:

ρ j = ( x j x ) 2 + ( y j y ) 2 + ( z j z ) 2 + a + ε j r , ( 3 )

где a = c δ t r ; ( 4 )

a - временной параметр, характеризующий расхождение временной шкалы приемника и системной шкалы времени спутниковой группировки.

Выражение (3) называется навигационным уравнением. В нем имеется четыре неизвестных (x, y, z, δtr).

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

Таким образом, для определения местоположения объекта необходимо, как минимум, четыре уравнения, т.е. четыре одновременно видимых спутника. Для повышения надежности и точности определения координат объекта современные навигационные приемники делаются двух системными, т.е. одновременно принимающими сигналы двух (ГЛОНАСС и GPS) орбитальных группировок.

Обычно данную систему линеаризуют, раскладывая в ряд Тейлора, в окрестности первого приближения (x0, y0, z0, δ t 0 r ), используя только линейные члены, то есть:

ρ j ( x , y , z , δ t r ) = ρ j ( x 0 , y 0 , z 0 , δ t 0 r ) + ρ j ' ( x , y , z , δ t r ) [ Δ x Δ y Δ z Δ δ t r ] , ( 5 )

где [ Δ x Δ y Δ z Δ δ t r ] = [ x y z δ t r ] [ x 0 y 0 z 0 δ t 0 r ]

Известны различные аналитические (безитерационные) способы вычисления координат объекта при использовании 4 спутников: способы, предложенные Краузе [1] в 1987 г, Клеусбергом [2] в 1994 г., Графаредом и Шаном [3] в 1996 г. Недостатками данных способов является невозможность использования в них избыточного количества спутников.

Наиболее близким к заявленному способу является способ, предложенный в 1985 г Банкрофтом [4]. Он предназначен для решения избыточной системы навигационных уравнений в случае одной орбитальной группировки, но не применим при одновременном использовании сигналов от двух группировок. Это связано с тем, что уход часов приемника от системного времени орбитальных группировок ГЛОНАСС и GPS различен. Кроме того, в навигационных приемниках пользователя обычно частоты гетеродинов в каналах ГЛОНАСС и GPS имеют независимые уходы частоты. Поэтому неизвестными становятся пять величин: координаты объекта x, y, z и времена δt1 и δt2 ухода часов приемника (или a 1, a 2 определяемых по формуле 4) в каналах ГЛОНАСС и GPS.

Актуальной является задача определения местонахождения объекта, когда часть навигационных спутников принадлежит орбитальной группировке ГЛОНАСС, а другая - GPS (например, имеется три спутника ГЛОНАСС и два GPS).

Типовая реализация системы определения местоположения объекта, наиболее распространенная в настоящее время, описана в патенте US 4672382 «Система определения местоположения» от 04.12.1985, выбранная в качестве аналога. Блок-схема системы приведена на фиг.1, где:

1. Орбитальная группировка навигационных спутников;

2. Приемная антенна навигационных сигналов;

3. Блок формирования измерений;

4. Блок условия;

5. Блок грубой начальной оценки;

6. Блок итерационных вычислений координат;

7. Блок оценки хранения координат в памяти;

8. Блок управления маршрутом.

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

Система навигационных уравнений (3), связывающая координаты объекта x, y, z и величину δt с текущими значениями координат спутников, является нелинейной. Она решается итерационным методом с помощью линеаризации уравнений в окрестности априорно известных координат объекта. Если количество спутников более четырех, число уравнений в системе становится избыточным. Применением метода наименьших квадратов оно сводится к четырем. Как правило, критерием сходимости итерационного алгоритма решения навигационных уравнений является минимизация некоего функционала.

Последовательность операций для аналога (фиг.1) показана на фиг.2, где:

9. Прием навигационного сигнала;

10. Первичная цифровая обработка навигационных сигналов (с выделением эфемерид и служебной информации) и формирование псевдодальностей;

11. Присвоение весов псевдодальностям, селекция спутников, коррекция псевдодальностей;

12. Прогноз координат спутников по актуальным эфемеридам на моменты времени излучения сигналов;

13. Условия существования первого приближения x, y, z, δtr;

14. Вычисление первого приближения x, y, z, δtr;

15. Решение линеаризованной системы навигационных уравнений, итерация вычисления очередных x, y, z, δtr;

16. Условие выполнимости Δδtr<ε, где ε - допустимая величина невязки по времени между шкалой приемника и системного времени ГНСС;

17. Норма невязки исходной системы должна быть меньше заданной величины;

18. Сглаживание результатов фильтром Калмана;

19. Запись вычислений и вывод x, y, z, δtr для пользователя.

Недостатками итерационных способов решения системы навигационных уравнений являются:

- необходимость задания начального приближения;

- вероятность несходимости вычислений для заданной точности;

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

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

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

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

- принимают сигналы от первой GPS и второй ГЛОНАСС орбитальных группировок навигационных спутников;

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

- рассчитывают скорректированные псевдодальности;

- прогнозируют координаты спутников по актуальным эфемеридам на моменты времени излучения навигационных сигналов;

- определяют координаты x, y, z, временной параметр a 1, соответствующий расхождению часов навигационного приемника от системного времени первой орбитальной группировки и временной параметр a 2, соответствующий расхождению часов навигационного приемника от системного времени второй орбитальной группировки методом наименьших квадратов, используя информацию от всех спутников с одинаковым весом первой и второй орбитальных группировок, путем аналитического решения системы нелинейных уравнений с минимизацией невязок;

- проверяют выполнимость условия a 1<ε, где ε - допустимая величина невязки по времени между шкалой навигационного приемника и системного времени глобальной навигационной спутниковой системы (ГНСС);

- подводят часы навигационного приемника до системного времени ГНСС;

- сглаживают результаты определения фильтром Калмана;

- записывают результаты определения x, y, z, a 1, a 2 и выводят их для пользователя.

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

- первую орбитальную группировку навигационных спутников GPS;

- вторую орбитальную группировку ГЛОНАСС;

- приемную антенну навигационных спутниковых сигналов от первой орбитальную группировки навигационных спутников GPS и второй орбитальной группировки ГЛОНАСС;

- блок формирования измерений;

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

- блок записи и вывода результатов вычислений на пользователя;

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

Блок вычислений координат и временных параметров содержит арифметико-логическое устройство с возможностью операций сложения, умножения и вычитания, первый вход которого соединен с выходом блока формирования измерений, блок расчета комплексно-сопряженной матрицы, вход которого соединен с выходом блока формирования измерений, блок расчета скалярных построчных произведений вектора с матрицей и матрицы с матрицей, вход которого соединен с выходом блока формирования измерений, второй и третий входы арифметико-логического устройства соединены с выходами комплексно-сопряженной матрицы и блока расчета скалярных построчных произведений вектора с матрицей и матрицы с матрицей соответственно, блок расчета корней полинома, вход которого соединен с выходом арифметико-логического устройства с возможностью операций сложения, умножения и вычитания, блок вычисления вектора состояний {x, y, z, a 1, a 2}, вход которого соединен с выходом блока расчета корней полинома, выход блока вычисления вектора состояний {x, y, z, a 1, a 2} является выходом блока вычислений координат и временных параметров.

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

20. Первая орбитальная группировка навигационных спутников GPS;

21. Вторая орбитальная группировка навигационных спутников ГЛОНАСС;

22. Приемная антенна навигационных сигналов;

23. Блок формирования измерений;

24. Блок вычислений координат и временных параметров;

25. Блок записи и вывода результатов вычислений на пользователя.

Блок-схема последовательности определения координат для заявленной системы, работающей по двум ГНСС, показана на фиг.4., где:

26. Прием навигационного сигнала;

27. Первичная цифровая обработка навигационных сигналов (с выделением эфемерид и служебной информации) и формирование псевдодальностей. Выполнение коррекции псевдодальностей;

28. Прогноз координат спутников по актуальным эфемеридам на моменты времени излучения сигналов;

29. Определение x, y, z, a 1, a 2 путем аналитического решения системы нелинейных уравнений с минимизацией невязок;

30. Условие выполнимости a 1<ε, где ε - допустимая величина невязки по времени между шкалой приемника и системного времени ГНСС;

31. Подвод часов приемника до системного времени ГНСС;

32. Сглаживание результатов фильтром Калмана;

33. Завершение вычисления и вывод x, y, z, a 1, a 2.

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

34. Арифметико-логическое устройство (сложение, умножение, вычитание);

35. Блок расчета комплексно-сопряженной матрицы;

36. Блок расчета скалярного произведения для пар векторов, вектор-матрица (построчно), матрица-матрица (построчно);

37. Блок расчета корней полинома;

38. Блок вычисления вектора состояний {x, y, z, a 1, a 2}.

В соответствии с предлагаемой последовательностью (фиг.4), на первых двух шагах происходит прием навигационной спутниковой антенной 22 и первичная цифровая обработка навигационных сигналов (с выделением эфемерид и служебной информации) и формирование псевдодальностей блоком формирования измерений 23. В блоке 23 так же проводится коррекция измеренной псевдодальности ρj на величину атмосферной задержки по общепринятым моделям, учет ухода часов спутников о системного времени. Окончательным этапом обработки данных в блоке 23 является прогноз координат спутников по актуальным эфемеридам на моменты времени излучения сигналов.

На следующем шаге выполняется определение {x, y, z, a 1, a 2} путем аналитического решения в блоке вычислений координат и временных параметров 24.

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

На последних этапах обработки при необходимости происходит сглаживание информации фильтром Калмана и вывод полученных значений {x, y, z, a 1, a 2}в блок записи и передачи результатов вычислений, выполненного с возможностью записи данных и передачи их пользователю.

Ниже описан алгоритм обработки данных в вычислителе.

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

Для орбитальной группировки под номером 21 (фиг.3.) (ГЛОНАСС), уравнение для j-го спутника в соответствии с (3) имеет вид:

( x j x ) 2 + ( y j y ) 2 + ( z j z ) 2 + a 1 = ρ j , ( 6 )

где [x, y, z] - искомые координаты приемника,

[xj, yj, zj) - координаты j-го спутника,

ρj - скорректированная псевдодальность до j-го спутника,

a 1 - временной параметр, характеризующий расхождение временной шкалы приемника от системного времени первой орбитальной группировки, умноженное на скорость света.

Для орбитальной группировки под номером 20 (фиг.2) (GPS), уравнение для g-го спутника в соответствии с (6) имеет вид:

( x g x ) 2 + ( y g y ) 2 + ( z g z ) 2 + a 2 = ρ g , ( 7 )

где {xg, yg, zg}- координаты g-го спутника,

ρj - скорректированная псевдодальность до g-го спутника,

a 2 - временной параметр, характеризующий расхождение временной шкалы приемника от системного времени второй орбитальной группировки, умноженное на скорость света.

После переноса временного параметра в правую часть уравнений (6), (7) и возведения в квадрат получим:

Далее необходимо учесть, что из каждой ОГ может наблюдаться разное кол-во спутников, пусть в ОГ ГЛОНАСС наблюдается количество спутников =n1, а в GPS - n2, общее кол-во спутников n=n1+n2. Отклонение шкалы приемника от ГЛОНАСС и GPS будут a 1 и а 2 соответственно. В результате получим систему из n уравнений вида (8).

Для решения данной системы необходимо перейти в пятимерное линейное пространство векторов L(R), такое, что для двух векторов r = { r 1 , r 2, r 3 , r 4 , r 5 } и t = { t 1 , t 2, t 3 , t 4 , t 5 } скалярным произведением r , t будет называться произведение вида:

r , t = r 1 t 1 + r 2 t 2 + r 3 t 3 r 4 t 4 r 5 t 5 r 4 t 5 r 5 t 4 ( 9 )

Определим следующий искомый 5-мерный вектор:

Решение системы уравнений вида (8).

Определим следующие 5-мерные векторы:

Где индекс k обозначает размерность единичного вектора e k .

Тогда система будет решаться относительно искомого вектора u .

Для выбранных векторов, в соответствии с (8) имеем:

при любом j

s j , s j = x j 2 + y j 2 + z j 2 p j 2 , при любом j

В этом случае систему навигационных уравнений для ГЛОНАСС и GPS можно записать в следующем виде:

Систему уравнений (14) можно переписать в матричном виде. Для этого необходимо перейти к операциям линейной алгебры, заменив скалярное произведение и арифметические операции с векторами в каждом отдельном уравнении системы на матрично-векторные операции.

Введем обозначения:

Λ = u , u 2 ; B = [ x 1 y 1 z 1 p 1 0 x n y n z n 0 p n ] ; B * = [ x 1 x n y 1 y n z 1 z n p 1 0 0 p n ] ;

P = [ p 1 0 0 0 0 0 0 p n ] ; D = d i a g ( B B * ) = [ x 1 2 + y 1 2 + z 1 2 p 1 2 x n 2 + y n 2 + z n 2 p n 2 ]

Λ - скалярная величина,

B и B* - матрица системы и комплексно сопряженная ей матрица,

D - вектор, составленный из значащих элементов диагональной матрицы, с элементами главной диагонали произведения BB*.

P - диагональная матрица, с элементами pj на главной диагонали.

Введем также обозначение составного вектора, следующего вида:

[ e k 1 e k 2 ] , например, [ e n 1 0 e n 2 ] = [ 1 1 0 0 ] } } n 1 э л е м е н т о в n 2 э л е м е н т о в .

Используя введенные обозначения, получим следующий вид системы навигационных уравнений (14):

,

после выражения обобщенного произведения Лоренца матричными операциями получим:

Данная система уравнений решается однозначно относительно неизвестного вектора u при количестве спутников, равном количеству неизвестных (пять), или по методу наименьших квадратов при избыточном количестве спутников. Примечательно, что если количество орбитальных группировок равно единице, то a 2=n2=0 и система уравнений сводится к системе уравнений по методу Банкрофта, как к частному случаю решения.

Решение системы будет иметь следующий вид:

u = B + ( Λ e n + 1 2 [ a 2 2 e n 1 a 1 2 e n 2 ] + a 1 a 2 e n + D 2 ) , ( 16 )

где B+=(BTB)-1BT

a 1; a 2 - параметры при поиске решения, а Λ найдем исходя из его определения Λ = 1 2 u , u . Можно заметить, что В+ в данном случае имеет размерность 5×n, а слагаемые в скобках являются векторами размерности n.

Произведение В+ на любое слагаемое в скобках дает вектор размерности 5. Таким образом (16) является системой из пяти уравнений с шестью неизвестными u = [ x y z a 1 a 2 ] и Λ. Шестое уравнение, позволяющее выразить Λ через параметры a 1, а 2 имеет вид:

Λ = 1 2 u , u ;

Λ = 1 2 B + ( Λ e n + 1 2 [ a 2 2 e n 1 a 1 2 e n 2 ] + a 1 a 2 e n + D 2 ) , B + ( Λ e n + 1 2 [ a 2 2 e n 1 a 1 2 e n 2 ] + a 1 a 2 e n + D 2 ) . ( 17 )

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

Λ = 1 2 B + ( Λ e n + 1 2 [ a 2 2 e n 1 a 1 2 e n 2 ] + a 1 a 2 e n + D 2 ) , B + Λ e n + 1 2 B + ( Λ e n + 1 2 [ a 2 2 e n 1 a 1 2 e n 2 ] + a 1 a 2 e n + D 2 ) , 1 2 B + [ a 2 2 e n 1 a 1 2 e n 2 ] + 1 2 B + ( Λ e n + 1 2 [ a 2 2 e n 1 a 1 2 e n 2 ] + a 1 a 2 e n + D 2 ) , a 1 a 2 B + e n + 1 2 B + ( Λ e n + 1 2 [ a 2 2 e n 1 a 1 2 e n 2 ] + a 1 a 2 e n + D 2 ) , 1 2 B + D ; ( 18 )

Путем алгебраических преобразований получим выражение вида:

Q 2 Λ 2 + Q 1 Λ + Q 0 = 0, ( 20 )

где Q2, Q1, Q0 имеют следующий вид:

Q 2 = 1 2 B + e n , B + e n ;

Решения квадратного уравнения с параметрами (20) имеют следующий вид:

Λ 1,2 = Q 1 ± Q 1 2 4 Q 2 Q 0 2 Q 2

Q 1 2 = 1 4 B + [ a 2 2 e n 1 a 1 2 e n 2 ] , B + e n 2 + a 1 a 2 B + [ a 2 2 e n 1 a 1 2 e n 2 ] , B + e n B + e n , B + e n + 1 2 B + [ a 2 2 e n 1 a 1 2 e n 2 ] , B + e n B + D , B + e n B + [ a 2 2 e n 1 a 1 2 e n 2 ] , B + e n + a 1 2 a 2 2 B + e n , B + e n 2 + a 1 a 2 B + e n , B + e n B + D , B + e n 2 a 1 a 2 B + e n , B + e n + 1 4 B + D , B + e n 2 B + D , B + e n + 1 ;

4 Q 2 Q 0 = 1 4 B + e n , B + e n B + [ a 2 2 e n 1 a 1 2 e n 2 ] , B + [ a 2 2 e n 1 a 1 2 e n 2 ] + a 1 a 2 B + e n , B + e n B + e n , B + [ a 2 2 e n 1 a 1 2 e n 2 ] + 1 2 B + e n , B + e n B + D , B + [ a 2 2 e n 1 a 1 2 e n 2 ] a 1 2 a 2 2 B + e n , B + e n 2 + a 1 a 2 B + e n , B + e n B + D , B + e n + 1 4 B + e n , B + e n B + D , B + D ;

Q 1 2 4 Q 2 Q 0

= 1 4 B + [ a 2 2 e n 1 a 1 2 e n 2 ] , B + e n 2 + 1 2 B + [ a 2 2 e n 1 a 1 2 e n 2 ] , B + e n B + D , B + e n B + [ a 2 2 e n 1 a 1 2 e n 2 ] , B + e n 2 a 1 a 2 B + e n , B + e n + 1 4 B + D , B + e n 2 B + D , B + e n + 1 1 4 B + e n , B + e n B + [ a 2 2 e n 1 a 1 2 e n 2 ] , B + [ a 2 2 e n 1 a 1 2 e n 2 ] 1 2 B + e n , B + e n B + D , B + [ a 2 2 e n 1 a 1 2 e n 2 ] 1 4 B + e n , B + e n B + D , B + D .

Полученные результаты подставляются в (16):

.

Вообще говоря, необходимо решать только последние два уравнения, так как первые три (x, y, z) уже выражены через параметры (a 1, a 2).

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

[ x y z a 1 a 2 ] = Q 1 ± Q 1 2 4 Q 2 Q 0 2 Q 2 B + e n + 1 2 B + [ a 2 2 e n 1 a 1 2 e n 2 ] + a 1 a 2 B + e n + 1 2 B + D

При общем количестве спутников n размерность матрицы В+ будет 5×n. Для выполнения дальнейших выкладок введем следующее обозначение:

B + = [ b 11 b n 1 b 51 b 5 n ] .

Последние два последних уравнения системы (16) в этом случае имеют вид:

Обозначим K = i = 1 n b 5 i / i = 1 n b 4 i и помножим на него уравнение (21), после чего складываем с (22). При этом иррациональность исключается и остается квадратное уравнение содержащее a 1 и а 2. Необходимо решить его, например, относительно a1, выразив его через а 2 и подставив в (22).

После подстановки получим уравнение:

Получаем квадратное уравнение, из которого найдем a 1 как функцию от а 2. Введем следующие обозначения

W 2 = ( i = n 1 + 1 n b 5 i + K i = n 1 + 1 n b 4 i ) ;

W1=-2K

Решение имеет следующий вид:

a 1 = w 1 ± w 1 2 4 w 2 w 0 2 w 2 .

Найденное значение подставляем в выражение (22):

Q 1 2 4 Q 2 Q 0 i = 1 n b 5 i = 2 Q 2 a 2 Q 1 i = 1 n b 5 i + Q 2 a 2 2 i = 1 n 1 b 5 i + Q 2 ( W 1 ± W 1 2 4 W 2 W 0 2 W 2 ) 2 i = n 1 + 1 n b 5 i + 2 Q 2 a 2 W 1 ± W 1 2 4 W 2 W 0 2 W 2 i = 1 n b 5 i + Q 2 i = 1 n ( x i 2 + y i 2 + z i 2 p i 2 ) b 5 i . ( 2 3 )

Подставляя в (23) в символьном виде значения Q2 и Q 1 2 4 Q 2 Q 0 , а также зависимости W0, W1, W2 и возводя обе части уравнения в квадрат получим выражение вида:

S 8 a 2 8 + S 7 a 2 7 + S 6 a 2 6 + S 5 a 2 5 + S 4 a 2 4 + S 3 a 2 3 + S 2 a 2 2 + S 1 a 2 + S 0 = 0 . ( 2 4 )

Коэффициенты полинома (24) не содержат неизвестных и вычисляются с использованием входных данных. В общем виде коэффициенты вычисляемы, но ввиду их громоздкости не приводятся. Используя среды мат. моделирования MATLAB/Maple коэффициенты последнего уравнения можно получить в символьном виде с целью дальнейшей реализации алгоритма на языках низкого уровня (C/C++).

Для нахождения корней полинома восьмой степени целесообразно использовать метод сопровождающей матрицы(companion matrix) [5].

Далее, используя полученные действительные значения а 2 вычислить a 1. А используя полученные параметры, получить решение x, y, z.

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

39. Формирование матрицы B и B+ для двух ОГ;

40. Расчет вспомогательных коэффициентов;

41. Расчет коэффициентов восьмого порядка;

42. Решение уравнения методом сопровождающей матрицы, отсев корней а 2 высокого порядка;

43. Вычисление времени a 1, а так же Q0, Q1, Q2;

44. Вычисление координат приемника x, y, z зная a 1, а 2, Q0, Q1, Q2;

45. Отсев лишних групп решений {x, y, z, a 1, а 2}исходя из физического смысла;

46. Оперативный контроль качества решения по невязкам исходной системы.

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

- частичная аппаратная реализация, где используется блок с оптимизированными векторными вычислениями (например, графический процессор) и стандартное арифметико-логическое устройство;

- полностью программная модель устройства.

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

Рассмотрим качество работы описанного выше способа и сравним с методом Банкрофта в двух случаях:

1) кол-во спутников равно количеству неизвестных (пять);

2) избыточное созвездие.

В качестве предварительного этапа подготовки исходных данных необходимо привести координатные сетки и временные шкалы в одно координатно-временное пространство, выбрав, к примеру, за основу шкалу времени GPS и координаты WGS-84. Рассчитать на текущий момент времени (момент измерения псевдодальностей) координаты спутников GPS и ГЛОНАСС в WGS-84. Большой проблемой в данном случае является матрица перехода из ПЗ-90 в WGS-84, так как нет единого стандарта. Но для метода решения системы нелинейных уравнений данная особенность значения не имеет. Так как с математической точки зрения система решается и решение должно обращать исходную систему в тождество при пяти спутниках и иметь минимум ошибки по МНК в случае избыточности. Далее необходимо составить матрицу В по спутникам GPS и ГЛОНАСС и отдельно по каждой группировке. Возьмем «тяжелый» случай, смешанную систему, в которой имеется всего 3 спутника GPS и 2 спутника ГЛОНАСС. Классический метод Банкрофта в этом случае не сможет дать результат отдельно по каждой из группировок, поэтому будут получены решения на созвездии из 4 спутников GPS (включающем указанные 3) и на созвездии из 4 спутников ГЛОНАСС (включающем указанные 2), далее составлено решение с осредненными координатами x, y, z по системам GPS и ГЛОНАСС, а так же добавлены временные поправки a 1 по системе GPS и a 2 по системе ГЛОНАСС. Невязки при подстановке указанного вектора-решения будут сравниваться с невязками, полученными по обобщенному методу для двух орбитальных группировок.

Исходные матрицы для решения классическим методом Банкрофта имеют следующий вид (x, y, z и скорректированная псевдодальность до спутника):

B G P S = ( 1 4 2 0 3 7 3 4 . 2 6 2 4 8 6 1 0 1 6 4 1 0 4 0 3 . 4 2 0 0 3 9 8 0 1 ` 5 2 0 3 7 5 7 . 4 5 3 4 2 9 6 0 2 0 7 3 6 0 3 3 . 6 2 6 8 7 3 7 1 4 0 0 3 9 8 6 . 4 9 1 4 0 1 1 0 7 8 6 8 9 3 4 . 5 9 0 0 4 5 8 7 2 1 2 0 4 2 9 7 . 5 9 3 2 5 7 0 2 0 2 7 7 9 1 4 . 1 2 7 3 9 9 6 4 5 0 5 8 3 5 . 5 8 3 9 9 7 1 3 2 3 0 1 2 3 5 3 4 . 0 9 3 8 4 1 3 0 1 2 0 3 4 0 3 0 . 8 0 4 9 2 6 5 0 2 3 0 9 0 6 3 7 . 2 1 2 6 1 4 0 1 1 0 2 0 3 1 2 . 9 4 5 8 0 7 8 1 1 1 3 2 6 9 9 . 3 8 7 2 5 8 1 0 2 1 6 8 5 1 5 2 / 8 0 4 2 8 2 0 0 2 3 2 8 5 2 1 9 / 4 3 0 4 8 0 9 )

Строки матрицы BGPS соответствуют системным номерам спутников GPS 20, 23, 30 и 31 соответственно.

B G L O N A S S = ( 1 0 8 2 7 4 8 5 . 5 1 9 7 6 7 1 5 6 9 6 3 6 5 . 7 5 0 1 6 9 2 1 6 9 6 6 2 0 5 . 6 6 4 9 1 5 9 1 9 5 7 3 0 3 1 . 8 4 2 3 6 7 5 - 5 5 6 4 6 9 4 . 6 4 3 0 5 5 6 1 9 0 5 1 6 4 4 . 9 3 8 6 9 1 7 2 2 3 1 7 3 4 7 4 . 4 5 9 2 6 7 4 2 0 9 5 1 2 6 3 . 2 1 7 9 7 5 6 2 1 9 0 8 2 7 7 . 7 2 5 4 9 1 1 1 2 8 8 2 8 7 7 . 6 1 1 4 2 4 1 - 1 2 8 8 0 0 0 . 6 2 6 6 7 7 2 7 2 2 8 0 7 3 9 2 . 1 2 5 2 6 7 7 1 6 5 8 4 6 5 1 . 9 6 6 4 5 1 9 - 1 5 0 2 2 1 8 6 . 1 3 7 1 0 6 6 1 2 3 1 7 6 1 9 . 1 5 7 9 2 9 6 2 3 1 3 8 0 6 8 . 4 0 2 9 6 2 2 )

Строки матрицы BGL0NASS соответствуют системным номерам спутников ГЛОНАСС 3, 13, 15 и 20 соответственно.

Полученные решения (см. табл.1, все единицы в метрах)

Таблица 1.
GPS ГЛОНАСС
х=2846250.36421181 х=2846135.03277489
y=2201677.85798526 y=2201641.31968233
z=5248681.95352191 z=5248736.41078770
a1=-18.33113694 а2=-129.08714075

Таким образом, можно взять за основу решение, полученное методом Банкрофта вида:

[ x y z a 1 a 2 ] = [ 2 8 4 6 1 9 2 . 6 9 8 4 9 3 3 5 2 2 0 1 6 5 9 . 5 8 8 8 3 3 7 9 5 2 4 8 7 0 9 . 1 8 2 1 5 4 8 0 1 8 . 3 3 1 1 3 6 9 4 1 2 9 . 0 8 7 1 4 0 7 5 ]

Вектор решения по обобщенному методу для двух ОГ, имеющих в составе 3 GPS(23, 30, 31)+2 ГЛОНАСС (15, 20):

[ x y z a 1 a 2 ] = [ 2 8 4 6 2 5 0 . 4 9 4 3 5 8 3 0 2 2 0 1 7 4 4 . 9 5 3 7 6 0 5 6 5 2 4 8 7 6 6 . 3 4 8 5 3 5 9 9 6 6 . 8 9 8 1 0 5 4 5 7 . 3 6 7 2 7 1 1 5 ]

Видно, что решения близкие, а в a 1, а 2 во втором случае вошли ошибки преобразования координат, вектора невязок и их нормы для смешанной исходной системы имеют следующие величины (см. табл.2):

Таблица 2.
Банкрофт Обобщенный метод
Невязки 15.41122860 0.0000000149
по спутникам, -9.89679900 0.0000000112
метры -46.55330177 0.0000000149
-64.55570337 0.0000000149
-19.80951382 0.0000000075
Норма вектора невязки, метры 84.03878972 0.0000000291

Очевидно, что в случае метода Банкрофта невязки, близкие к нулю, получаются только на первых трех уравнениях (GPS), так как координаты и время соответствуют решению только по GPS, а от решения ГЛОНАСС взята только временная поправка.

В случае с обобщенным методом получены близкие к нулю невязки, но не нулевые, так как это обусловлено точностью разрядной сетки. Норма вектора невязки для типовых случаев расположения спутников, как правило, не превышает 10-6 м, при вычислениях с плавающей запятой двойной точности. Данная точность является достаточной для высокоточных приложений. Однако, исходя из значений поправок характеризующих время в векторе решения, следует серьезно проработать значения матриц перехода в единое координатно-временное пространство, так как это будет источником ошибок, уводящим временные поправки от истинных значений.

В случае избыточности ситуация с методом Банкрофта еще сложнее:

Полученные решения (см. табл.3, все единицы в метрах)

Таблица 3.
GPS ГЛОНАСС
х=2 846251.46141769 х=2846072.63300129
y=2 201709.04624132 y=2201646.99218817
z=5 248717.76349039 z=5248658.14254719
a1=18.79691716 а2=-171.02294209

Таким образом, можно взять за основу решение, полученное методом Банкрофта вида:

[ x y z a 1 a 2 ] = [ 2 8 4 6 1 6 2 . 0 4 7 2 0 9 4 9 2 2 0 1 6 7 8 . 0 1 9 2 1 4 7 5 5 2 4 8 6 8 7 . 9 5 3 0 1 8 7 9 1 8 . 7 9 6 9 1 7 1 6 1 7 1 . 0 2 2 9 4 2 0 9 ]

Вектор решения по обобщенному методу для двух ОГ, имеющих в составе избыточную группировку GPS/ГЛОНАСС:

[ x y z a 1 a 2 ] = [ 2 9 4 6 2 2 0 . 9 9 1 9 1 0 2 6 2 2 0 1 7 0 1 . 4 6 6 5 4 8 9 2 5 2 4 8 7 5 8 . 7 5 9 5 8 4 3 9 3 4 . 2 6 2 2 5 1 5 5 6 9 . 5 9 0 6 8 5 2 5 ]

Нормы векторов невязок для смешанной исходной системы имеют следующие величины (см. табл.4):

Таблица 4.
Банкрофт Обобщенный метод
Норма вектора невязки, метры 205.35810742 108.57034441

Как видно, обобщенный метод эффективно работает и может быть рекомендован для вычисления координат точек с минимизацией вектора ошибки по МНК сразу по всей смешанной орбитальной группировке ГЛОНАСС/GPS, чего нельзя добиться используя классический метод Банкрофта.

Список используемой литературы

1. Krause, L.O. (1987): A Direct Solution to GPS-Type Navigation Equations. IEEE Transactions on Aerospace and Electronic Systems, Vol. AES-23, No. 2, pp.225-232

2. Kleusberg, A. (1994): Die direkte Lösung des räumlichen Hyperbelschnitts. Zeitschrift fur Vermessungswesen, Vol. 119, No. 4, pp.188-192

3. Grafarend, E.W. and Chan, J. (1996): A Closed-form Solution of the Nonlinear Pseudo-Ranging Equations (GPS), ARTIFICIAL SATELLITES Planetary Geodesy, No. 28, pp.133-147

4. Bancroft, S. (1985): An Algebraic Solution of the GPS Equations. IEEE Transactions on Aerospace and Electronic Systems, Vol. AES-21, No. 6, pp.56-59

5. Метод сопровождающей матрицы: http://en.wikipedia.org/wiki/Companion_matrix

1. Способ определения местоположения объекта при использовании глобальных навигационных спутниковых систем, заключающийся в том, что:
принимают сигналы от первой GPS и второй ГЛОНАСС орбитальных группировок навигационных спутников;
производят первичную цифровую обработку навигационных сигналов с выделением эфемерид и псевдодальностей;
рассчитывают скорректированные псевдодальности;
прогнозируют координаты спутников по актуальным эфемеридам на моменты времени излучения навигационных сигналов;
определяют координаты x, y, z, временной параметр а1, соответствующий расхождению часов навигационного приемника от системного времени первой орбитальной группировки, и временной параметр а2, соответствующий расхождению часов навигационного приемника от системного времени второй орбитальной группировки методом наименьших квадратов, используя информацию от всех спутников с одинаковым весом первой и второй орбитальных группировок, путем аналитического решения системы нелинейных уравнений с минимизацией невязок;
проверяют выполнимость условия a1<ε, где ε - допустимая величина невязки по времени между шкалой навигационного приемника и системного времени глобальной навигационной спутниковой системы (ГНСС);
подводят часы навигационного приемника до системного времени ГНСС;
сглаживают результаты определения фильтром Калмана;
записывают результаты определения x, y, z, а1, а2 и выводят их для пользователя.

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

3. Система по п.2, отличающаяся тем, что блок вычислений координат и временных параметров содержит арифметико-логическое устройство с возможностью операций сложения, умножения и вычитания, первый вход которого соединен с выходом блока формирования измерений, блок расчета комплексно-сопряженной матрицы, вход которого соединен с выходом блока формирования измерений, блок расчета скалярных построчных произведений вектора с матрицей и матрицы с матрицей, вход которого соединен с выходом блока формирования измерений, второй и третий входы арифметико-логического устройства соединены с выходами комплексно-сопряженной матрицы и блока расчета скалярных построчных произведений вектора с матрицей и матрицы с матрицей соответственно, блок расчета корней полинома, вход которого соединен с выходом арифметико-логического устройства с возможностью операций сложения, умножения и вычитания, блок вычисления вектора состояний {x, y, z, а1, а2}, вход которого соединен с выходом блока расчета корней полинома, выход блока вычисления вектора состояний {x, y, z, а1, а2} является выходом блока вычислений координат и временных параметров.



 

Похожие патенты:

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

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

Изобретение относится к устройству контроля за местонахождением лиц в системах туннелей. .

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

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

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

Изобретение относится к навигации и определению местоположения устройства. .

Изобретение относится к радиолокации и может быть использовано для повышения точности определения местоположения наземных источников радиоизлучений (ИРИ) в пассивных режимах работы радиолокационных станций (РЛС) или станций радиотехнической разведки (СРТР).

Использование: изобретение относится к области звуколокации и радиолокации и может быть использовано для решения научных и прикладных задач, в частности для обнаружения подводных объектов. Сущность: в некоторой точке океана располагается надводный или подводный корабль, который излучает звуковую волну с мощностью Iизл=5*105 Вт/м2, на частоте fзвук=10 кГц. Это излучение распространяется во все стороны и на расстоянии Lдет=30 км от корабля создает звуковое давление порядка p1=17 Вт/м2. Звуковая волна, отражаясь от подводного объекта с коэффициентом отражения котр=10-2, за счет сжимаемости воды создает дифракционную решетку, соответствующую цилиндрической звуковой волне. Высокочастотные генераторы, с мощностью Рген=500 МВт, работающие на частоте fрадио=108 Гц, расположенные на одной группе самолетов, облучают отдельные участки поверхности воды узким лучом радиоволн. Отражение в первом порядке от дифракционной решетки, созданной цилиндрической звуковой волной приводит к появлению отраженных волн. Приемники распространяющегося в узком луче излучения, расположенные на другой группе самолетов, с чувствительностью 3*10-21 Вт, при площади антенн Sант=700 м2, регистрируют мощность принимаемого излучения ~10-19 Вт. Благодаря тому, что рассеяние происходит на бегущей решетке, отраженная от нее электромагнитная волна оказывается Допплеровски сдвинутой на величину δf=100 Гц. По зарегистрированному ифракционному излучению определяют координаты подводного объекта. Технический результат: увеличения дальности обнаружении подводных объектов. 1 ил.

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

Изобретение относится к области радиотехники и касается акустооптического интерферометра. Акустооптический интерферометр состоит из антенной решетки, источника когерентного излучения, коллиматора, акустооптического модулятора с четырьмя пьезопреобразователями, фурье-линзы, матричного фотоприемника и цифрового процессора. Антенная решетка содержит две пары ненаправленных приемных элементов, расположенных в одной плоскости так, что линии, соединяющие приемные элементы каждой пары, перпендикулярны друг другу. Выходы первой пары приемных элементов антенной решетки соединены с первой парой пьезопреобразователей непосредственно, а выходы второй пары приемных элементов антенной решетки соединены со второй парой пьезопреобразователей через фазовращатели на 90°. Технический результат заключается в увеличении сектора однозначно определяемых углов прихода радиоизлучения до 360 градусов. 3 ил.

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

Изобретение предназначено для определения местоположения источников радиоизлучения (ИРИ). Достигаемый технический результат - повышение точности определения координат местоположения ИРИ. Способ основан на использовании измерений на радиоконтрольных постах значений уровней сигналов (УС) на каждой из выбранных частот и обратно пропорциональной зависимости отношений расстояний от поста до ИРИ и соответствующих им УС, на которых усредняют полученные значения и вычисляют текущую дисперсию УС на каждой из выбранных частот до тех пор, пока текущая дисперсия не станет больше предыдущей, затем усредненные значения УС передают на базовый пост, где получают их отношения и составляют три уравнения, каждое из которых описывает окружность с центрами местоположения постов и радиусами обратно пропорциональными УС и является линией положения, а также определитель Кэли-Менгера, по которому и отношениям усредненных значений УС определяют расстояние от ИРИ до постов, а по двум любым парам составленных уравнений определяют текущее среднее значение широты и долготы местоположения ИРИ как координаты точки пересечения радикальных осей окружностей, то есть как координаты радикального центра линий положения. Текущее среднее значение широты и долготы местоположения ИРИ определяют до тех пор, пока разность двух смежных значений текущих сумм дисперсий широты и долготы местоположения ИРИ не изменит свой знак, после чего усредненные значения координат местоположения ИРИ фиксируют как окончательные. 1 ил.

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

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

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

Изобретение относится к способам обработки радиолокационных изображений (РЛИ). Достигаемый технический результат - повышение быстродействия обработки РЛИ. Сущность изобретения состоит в следующем. При зондировании участка земной поверхности с помощью радиолокатора с синтезированной апертурой (РСА), установленного на носителе в виде ЛА, получают отраженный сигнал от земной поверхности, одновременно с получением сигнала определяют с помощью навигационной системы ЛА пространственное положение фазовых центров антенн (ФЦА) и запоминают его. Полученный сигнал на входе РСА представляют в виде суммы радиоизображений объекта, фона и шума наблюдения. При этом фон, на котором расположен объект в совокупности с шумами наблюдения, рассматривают как некоторый эквивалентный шум. Для совместного различения и оценки параметров (координат) используют байесовский метод, предполагающий совместную оптимизацию этих двух операций. В соответствии с байесовским правилом оптимальности необходимо минимизировать апостериорный риск по двум параметрам: оценке дискретного параметра неопределенности i - определить объект, и оценке параметров (координат) объекта, где i - тип объекта. Совместная минимизация риска может быть выполнена в два этапа: сначала по условной оценке параметров (координат) объекта при фиксированном значении i, а затем по всем i. Определение параметров (координат) образа объектов в данном алгоритме предшествует различению самих объектов, однако байесовская оценка формируется после определения i-го объекта. Условную оценку координат объекта получают по методике, приведенной с использованием эталонных моделей объектов, формируемых предварительно. Для различения объектов необходимо выполнить минимизацию апостериорного риска по всем возможным i-типам объектов. При этом алгоритм различения сводится к сравнению усредненных отношений правдоподобия с набором пороговых значений, которые формируются предварительно для всех типов объектов. 1 ил.

Изобретение относится к области маркшейдерско-геодезического мониторинга и может быть использовано для обеспечения безопасности разработки месторождений нефти и газа. Согласно заявленному решению на исследуемой территории проводят геодезические измерения и определяют смещения Ngeod геодезических реперов на север U i n , восток U i e и по вертикали U i v (i=1, 2,…, Ngeod). За тот же интервал времени определяют смещения Nsat устойчиво отражающих площадок в направлении на спутник U j L O S (j=1, 2,…, Nsat) с помощью радарной спутниковой интерферометрии. После чего осуществляют разбивку разрабатываемого месторождения на K элементарных объемов. Рассчитывают смещения в точке j-й устойчиво отражающей площадки в направлении на спутник V j , k L O S , которые возникают в результате увеличения давления на единицу в k-м элементарном объеме, и смещения в точке i-го геодезического репера соответственно на север, восток и по вертикали V i , k n , V i , k e и V i , k v , которые возникают в результате увеличения давления на единицу в k-м элементарном объеме. Определяют в каждом объеме изменения давления ΔPk. После чего определяют три компоненты вектора смещений земной поверхности. Технический результат - повышение точности определения смещений земной поверхности. 4 ил.
Наверх