Международный информационно-аналитический журнал «Crede Experto: транспорт, общество, образование, язык». № 2 (13). Июнь 2017 (http://ce.if-mstuca.ru)
УДК 621.396.98 ББК 32.95 К438
В. В. Кирюшкин Воронеж, Россия Е. А. Дьяконов Воронеж, Россия
ИССЛЕДОВАНИЕ ТОЧНОСТИ ОПРЕДЕЛЕНИЯ КООРДИНАТ ВОЗДУШНОЙ ЦЕЛИ В СПУТНИКОВОЙ-ПСЕВДОСПУТНИКОВОЙ МНОГОПОЗИЦИОННОЙ СИСТЕМЕ НАБЛЮДЕНИЯ
В настоящей работе проведено исследование точности определения координат воздушной цели в многопозиционной системе наблюдения, использующей для подсвета воздушных целей сигналы навигационных спутников глобальной навигационной спутниковой системы и сигналы псевдоспутников. Показано, что использование сигналов псевдоспутников позволяет в 8-10 раз уменьшить пиковую и в 2-3 раза среднюю погрешность измерения координат воздушной цели в спутниковой многопозиционной системе наблюдения.
Ключевые слова: глобальная навигационная спутниковая система, многопозиционная система наблюдения, погрешность измерения координат, зона действия.
V. V. Kiryushkin Voronezh, Russia E. A. D'yakonov Voronezh, Russia
A STUDY OF THE ACCURACY OF THE DEFINITION OF COORDINATES
© Кирюшкин В. В., Дьяконов Е. А., 2017
OF AIR TARGETS IN SATELLITE-PSEUDOLITE MULTIPOSI-TIONAL SURVEILLANCE SYSTEM
The article describes the study of the accuracy of determination of coordinates of air targets in multipositional surveillance system using the signals of navigation satellites of global navigation satellite systems and pseudolite signals to illuminate air targets. It is shown that the use of pseudolite signals allows to reduce the peak error 810 times and decrease the average error of measurement of coordinates of air targets 2-3 times in a multipositional satellite surveillance system.
Key words: global navigation satellite system, multipositional surveillance system, the error of measurement of coordinates, operating zone.
В многопозиционной системе наблюдения {рис. 1), использующей для подсвета воздушных целей сигналы навигационных спутников глобальной навигационной спутниковой системы (ГНСС) Cj (j=1, 2, ..N, где N > 4) с координатами
(x y ■ z■)
j j j' [Kiryushkin, 2016] и сигналы псевдоспутников ПС/ (¿=1, 2, ..., M) с
(xf, yfc, zfc ) u
координатами v г г ' в наземном приемнике П с известными координатами (%п'Уп'Zn), наряду с навигационными сигналами прямого распространения принимаются сигналы, рассеянные воздушной целью Ц с неизвестными координатами (x'y'z), находящейся в зоне действия многопозиционной системы наблюдения.
Рис. 1. Многопозиционная система наблюдения «навигационные спутники/псевдоспутники - воздушная цель - наземный приемник» Навигационные сигналы, рассеянные воздушной целью, выделяют на фоне навигационных сигналов прямого распространения одним из оценочно-корреляционно-компенсационных методов [Кирюшкин, 2016 а]. По рассеянным сигналам осуществляется измерение дальностей 3 вдоль пути распространения «/-и
навигационный спутник - воздушная цель - наземный приемник»
«
'• = Г/ + Г0 + г, (1)
г,
где 3
п1/2
"(х, - х)2 + (у. - у)2 + (г, - г)21
3 3 3 -I - дальность пути «/-и навигационныи
спутник - воздушная цель»,
го =
_(хп - х)2 + (уп - у)2 + (гп - г)21
дальность
пути «воздушная цель - наземный приемник»; г - случайная погрешность из- гР8
мерения дальности, и измерение дальностей г вдоль пути распространения «г-й псевдоспутник - воздушная цель - наземный приемник»
г РЯ = г ПС + + г г =' + г0 + г, (2)
Г.пс
(хГ - X)2 + (у?с - у)2 + (г?с - г)21 _
г г -1 - дальность пути «г-и псевдо-
спутник - воздушная цель».
В канале вторичной обработки приемника на основе измеренных дальностей / /
3 и 3 с использованием разностно-дальномерного метода и метода наименьших квадратов определяются неизвестные координаты воздушной цели. Для этого одно из измерений, например, измерение для 1-го спутника, выбирается в качестве опорного и вычитается из измерений сигналов всех остальных спутников и псевдоспутников.
Тогда для N спутников, для которых выполняется условие просветной радиолокации (значение бистатического угла лежит в пределах -130 -180 [ Черняк,
1993], см. рис. 1), получают систему из N-1 уравнений
АГ5 = г5 - г5 + г 'з1 Г3 г +6, (3)
Аг р
где } = 2, 3, ..N - номер спутника; п - разность измеренных дальностей от цели
до у-го и до 1-го навигационных спутников.
Соответственно, для М псевдоспутников, для которых выполняется условие просветной радиолокации, получают систему из М уравнений
Аг^ _ Р5 _ 1 ,р
Аг 1 _г г\ +£, (4)
Аг
где I = 1, 2, ..., М - номер псевдоспутника; 11 - разность измеренных дальностей от цели до г-го псевдоспутника и до 1-го навигационного спутника.
Для решения полученной системы нелинейных уравнений воспользуемся итерационным методом наименьших квадратов [Шабшеевич, 1982], для чего запишем системы нелинейных уравнений (3) и (4) в обобщенном виде
К _К(я,05), (5)
Кря _ КРЯ (Ч, 0РЯ ) , (6)
та
где 5 - (^1)-мерный вектор измерений разностей дальностей до спутников;
та
ж - М-мерный вектор измерений разностей дальностей до псевдоспутников;
а _\х у I"Т - Ос
^ J - вектор истинных координат воздушной цели; - матрица ко-
ординат навигационных спутников; 0р5 - матрица координат псевдоспутников.
Результирующий (#+М-1)-мерный вектор К измерений разностей дальностей определяется как
К _\КякРЗ ] , (7)
а результирующее уравнение измерений в матричном виде запишется
к _ К(а, О), (8)
где 0 _ \О50р ] - общая матрица координат навигационных маяков (спутников
и псевдоспутников).
Решение уравнения (8) представляет собой процесс многократной (итерационной) обработки результатов навигационных измерений [Шабшеевич, 1982]:
я к _ я к _1+( С^С _1 )_1 С^К, _1, (9)
где
Як_1 _ ( Хк_1, ук_1, к_1 )
оценка вектора координат цели на к-1 итерации вы-
_ та _ та „
числений; к~1 к~1 - вектор разностей (невязок) измеренных К и рас-
К С
считанных на к-1 итерации к_1 разностей дальностей; Ск_1 - матрица наблюдений; Р - матрица, обратная корреляционной матрице погрешностей измерений дальностей.
Для бистатических разностно-дальномерных измерений элементы вектора
ста
невязок к-1 будут вычисляться по формулам (10) - для '-го спутника и (11) -для г-го псевдоспутника:
Я(Агкк_1) _ А' _ г5к_1+ Г1_1, (10)
s(АгpS )_Агр5 _ Гр + Г5
где
г5
1,к _1
ГРЯ _
' г,к_1
( _хк_1) +(У' _ ук_1) + (¿' _ ¿к_)
(х1 _ хк _ )2 + (У1 _ ук _ )2 + (¿1 _ 1к _1 )2
(хГ _ Хк_ )2+(У?с _ Ук_ )2+(_ ¿к_)
1/2
(12)
(13)
(14)
С
Матрица наблюдений к_ имеет размерность (Ы+М- 1)х3 и определяется следующим образом
Ск _1 _ С21,к _1 •••С '1,к _1 •••С т,к _1 •••С( N+г )1,к _1 "•С( N+М )1,к _1
С '1,к _1
(ео8 а5 _ ео8 а?) (ео8 (( _ ео8 (() (ео8 у? _ ео8 у?)
5 _ У' _ ук_1
Я х' хк _1 о
ео8а; - со8(
Г5
Г5
С08 у.
Я _ ¿' ¿к_1
Г5
',к _1
со8а? _со8(( _^ ™у _
'5 Г',к _1
'5 Г',к _1
С08 у1
'Я Г' ,к _1
СР5
С( N+г )1,к _1
(С08 арЯ _ С08 а? ) (С08 ( _ С08 (( ) (С08 урЯ _ С08 у?)
(15)
(16)
(17)
(18) (19)
5
ПС _ ~ ПС _ ~ ПС _ ~
cosa?S = x .psXk-1 cos= -1 cosy?s = z .psZk-1
ri ,k -1 ; ri,k -1 ; -1 (20)
Точность оценки вектора q координат воздушной цели определяется корреляционной матрицей погрешностей [Шабшеевич, 1982]:
K9 = (cLPCk-1) . (21)
Описанный алгоритм реализован в среде программирования MATLAB. Исследование проводилось для истинной траектории воздушной цели, сформированной в среде авиасимулятора FlightGear и показанной на рис. 2 сплошной кривой. Жирным треугольником на рис. 2 обозначен наземный приемник, а звездочками - псевдоспутники. Среднеквадратическое отклонение (СКО) ошибки раз-ностно-дальномерных измерений было выбрано 5 м и считалось постоянным на протяжении всего сеанса измерений.
Рис. 2. Геометрия численного эксперимента Результаты проведенного численного эксперимента показаны на рис. 3 в виде временной зависимости абсолютной радиальной погрешности формируемой оценки
АЯ = Г Ах2 + Ау2 + Аг21Ш ^ . /
координат I- и СКО формируемой оценки Я (жирная ли-
ния), где Ах = х - х, АУ = у - У, Аг = £ - г, ^ + *У + Г, ^у2^ - элементы формируемой корреляционной матрицы К4 погрешностей.
Рис. 3. Временная зависимость радиальной погрешности оценки местоположения воздушной цели с использованием псевдоспутников
Аналогичные зависимости, полученные без использования псевдоспутников [Кирюшкин, 2016 б], показаны на рис. 4.
Анализ рис. 3-4 показывает, что использование сигналов псевдоспутников дополнительно с сигналами навигационных спутников позволяет в 8-10 раз уменьшить пиковую и в 2-3 раза среднюю погрешность измерения координат воздушной цели.
Однако и в первом и во втором случае точность формируемых оценок варьируется в очень широких пределах, от единиц до сотен метров, что связано с изменением геометрических условий наблюдения.
В соответствии с [Кирюшкин, 2016 в] решающим условием в увеличении значения геометрического фактора в полученной системе наблюдения является количество и взаимное расположение спутников, попадающих в узкий диапазон телесных углов, образуемых конусом с вершиной в точке расположения приемника и центром основания, совпадающим с местом положения воздушной цели. При этом значение телесного угла конуса определяется бистатическими услови-
Я = 130°
ями (^тш 130 ), а положение конуса будет изменяться в процессе полета воздушной цели. Такой конус можно рассматривать как «луч диаграммы направленности антенны» полученной системы наблюдения.
Рис. 4. Временная зависимость радиальной погрешности оценки местоположения воздушной цели без использования псевдоспутников Для оценки зоны действия полученной системы наблюдения было зафиксировано время эксперимента = 2000 с и высота полета цели к = 2000 м. Изменяя азимут Ф и угол места 9 наблюдения цели из места положения наземного приемника и определяя истинные координаты цели (х'у'г) как точку пересечения «луча диаграммы направленности антенны» с плоскостью на высоте к, для каждого значения азимута Ф было получено значение минимального угла места 9, при котором выполняются бистатические условия наблюдения и возможна оценка координат цели предложенным способом. Полученные зоны действия показаны на рис. 5.
Анализ рис. 5 показывает, что использование спутников только одной ГНСС ГЛОНАСС позволяет успешно оценивать координаты воздушной цели только в узком диапазоне азимутов и углов места (для рассмотренного момента времени соответственно 80°-160° и 45°-90°, см. рис. 5, а).
Увеличение орбитальной группировки навигационных спутников позволяет существенно расширить зону действия системы наблюдения. Так, при использовании
двух систем ГЛОНАСС и GPS обеспечивается успешное определение координат воздушной цели, наблюдаемой на всех азимутах, однако в диапазоне азимутов 190°-360° зона действия ограничена минимальными значениями угла места 30°-45° (рис. 5, б).
При использовании полной перспективной группировки навигационных спутников ГНСС ГЛОНАСС, GPS, GALILEO и COMPASS минимальные значения угла места уменьшаются до 15°-30° (рис. 5, <?).
Использование дополнительно к полной группировке спутников восьми псевдоспутников позволяет расширить зону действия многопозиционной системы наблюдения практически до полусферы (рис. 5, г).
270
180
б)
180
180
В)
а) спутников ГЛОНАСС; б) спутников ГЛОНАСС и GPS; в) спутников ГЛО-НАСС, GPS, GALILEO и COMPASS; г) спутников ГЛОНАСС, GPS, GALILEO
Библиографический список
1. Кирюшкин В. В. Способ обнаружения и оценки радионавигационных параметров сигнала космической системы навигации, рассеянного воздушной целью, и устройство его реализации / В. В. Кирюшкин,. Д. А. Черепанов, А. А. Дисенов, С. С. Ткаченко // Заявка №2014101847 (002722) от 21.01.2014 г. Патент №2591052 от 17.06.2016 а.
2. Кирюшкин, В. В. Определение координат воздушной цели в многопозиционной системе наблюдения «навигационные спутники - воздушная цель - наземный приемник» / В. В. Кирюшкин, Д. А. Черепанов, Е. А Дьяконов // Теория и техника радиосвязи. № 2. 2016 в. С. 29-35.
3. Кирюшкин В. В. Влияние геометрических условий наблюдения на возможность и точность оценки координат воздушной цели в многопозиционной системе «навигационные спутники - воздушная цель - наземный приемник» / В. В. Кирюшкин, Д. А. Черепанов, В. И. Ко-стылев // Радиолокация, навигация и связь: XXII международная научно-техническая конференция. Воронеж. 2016 б. Т. 2. С. 1395-1408.
4. Черняк В. С. Многопозиционная радиолокация / В. С. Черняк. М.: Радио и связь, 1993. 416
с.
5. Шабшеевич В. С. Сетевые спутниковые радионавигационные системы / Под ред. П. П. Дмитриева, В. С. Шабшеевича. М.: Радио и связь, 1982. 272 с.
и COMPASS и 8 псевдоспутников
6. Kiryushkin V. V. Coordinates estimation of the air target in the multiitem observation system «navigation satellites - the air target - the ground receiver» / Kiryushkin V.V., Cherepanov D.A. // J. Sib. Fed. Univ. Eng. technol. 2016. 9(8). P. 1172-1182. DOI: 10.17516/1999-494X-2016-9-8-1172-1182. References
1. Kiryushkin, V. V. (2016). The method of detection and evaluation of navigation parameters of the signal space navigation systems, ambient due to the air target, and device for its implementation / V. V. Kiryushkin,. D. A. Cherepanov, A. A. Disenov, S. S. Tkachenko // Application No. 2014101847 (002722) from 21.01.2014 G. Patent No. 2591052 from 17.06.2016. (In Russian).
2. Kiryushkin, V. V. (2016). Determination of coordinates of air targets in the multipositional surveillance system "navigation satellites - aerial target - ground-receiver" / V. V. Kiryushkin, D. A. Cherepanov, Ye.A. D'yakonov // Theory and technique of radio communication. No. 2. 2016 V. P. 29-35. (In Russian).
3. Kiryushkin, V. V. (2016). Influence of geometrical conditions of observation and the accuracy of the estimation of coordinates of air targets in a multi-positional system "navigation satellites - aerial target - ground-receiver" / V. V. Kiryushkin, D. A. Cherepanov, V. I. Kostylev // Radar, navigation and communication: XXII international scientific-technical conference. Voronezh. 2016 b. Vol. 2. P. 1395-1408. (In Russian).
4. Chernyak, V. S. (1993). Multipositional radiolocation / V. S. Chernyak. Moscow: Radio and communication, 1993. 416 p. (In Russian).
5. Shabsheyevich V. S. (1982). Network satellite navigation system / Ed. by P. P. Dmitriev, V. S. Shabsheyevich. Moscow: Radio and communication, 1982. 272 p. (In Russian).
6. Kiryushkin V. V. Coordinates estimation of the air target in the multiitem observation system «navigation satellites - the air target - the ground receiver» / Kiryushkin V.V., Cherepanov D.A. //J. Sib. Fed. Univ. Eng. technol. 2016. 9(8). P. 1172-1182. DOI: 10.17516/1999-494X-2016-9-8-1172-1182.