УДК 621.396.96
Метод определения координат сканирующего источника радиоизлучения в многопозиционной радиолокационной системе с некооперируемым подсветом
Лешко Н.А.*, Петроченков Д.М.**, Мозгонов М.Ю.***, Федотов А.В.****
Ярославское высшее военное училище противовоздушной обороны, Московский
проспект, 28, Ярославль, 150001, Россия *e-mail: nikolai_zru@mail. ru * *e-mail: pdm [email protected] ***e-mail: [email protected] ****e-mail: [email protected]
Статья поступила 01.10.2019
Аннотация
В радиолокации широко используются многопозиционные пассивные методы определения координат источников радиоизлучения. Одни из методов требуют точного знания временной и фазовой структуры сигналов, другие используют динамические свойства источника излучения и дополнительные пространственные измерения угловых координат. В обоих случаях следствием является сложность аппаратуры, способной реализовать необходимые измерения с заданным качеством. Поэтому разработка методов измерения координат источников радиоизлучения, позволяющих сократить число измерений является научно-технической проблемой, требующей решения. В статье предложен метод вычисления координат сканирующего источника радиоизлучения, основанный на измерении временных интервалов между последовательным облучением приемных пунктов многопозиционной системы. Основное достоинство предложенного метода
определения координат сканирующего источника - это минимальное количество
измерений, отсутствие измеряемых угловых величин, независимость от
пространственной ориентации пунктов приема относительно друг друга. Для
определения координат источника излучения необходимы измерения временных
задержек последовательного облучения пунктов приема и дальности между ними.
Ключевые слова: многопозиционная радиолокационная система, некооперируемый источник подсвета, метод определения координат источника радиоизлучения.
Введение
Анализ тенденций развития радиолокационной техники показывает, что одним из перспективных направлений является создание систем, объединяющих средства активной и пассивной локации, а также радиотехнической разведки. На основе такой интеграции разрабатываются полуактивные многопозиционные радиолокационные системы (МПРЛС) с кооперируемым (КИП) или некооперируемым источником подсвета (НИП) [1, 2, 3, 4, 5, 6, 7, 8, 9].
При использовании некооперируемого источника подсвета воздушных объектов процедуру обнаружения цели можно условно разделить на два этапа. Первый этап заключается в определении координат и параметров источника подсвета. Второй этап - обнаружение отраженных от целей сигналов источника подсвета на фоне его прямого излучения. Причем, эффективность работы МПРЛС с НИП непосредственно зависит от качества решения задачи на первом этапе.
В [10, 11, 12, 13, 14, 15, 16, 17] предложены различные методы определения
координат источников радиоизлучения, которые на сегодняшний момент уже стали
традиционными. Эти методы основываются на когерентной обработке сигналов,
принятых различными позициями. В [18] рассматривается применение
искусственных нейронных сетей (ИНС) для решения задачи формирования оценки
направления прихода сигнала в пассивных системах радиолокации, в [19]
приводятся оценки погрешности измерения положения источника излучения
относительно измерителей.
В МПРЛС с НИП условия пространственной и временной когерентности
принятых сигналов в различных приемных пунктах, как правило, не выполняются
[20, 21]. Поэтому необходимо искать методы определения координат, не зависящие
от таких условий. Одним из динамических факторов, который может быть
использован для определения местоположения источника подсвета, является режим
сканирования диаграммы направленности (ДНА) его передающей антенны.
В [20, 21, 22] описаны варианты временного способа местоопределения
сканирующих источников радиоизлучения (ИРИ) при наблюдении из нескольких
приемных пунктов. Так, в [20] предложен способ определения дальности до ИРИ
при известной конфигурации размещения приемных пунктов по измеренным
пеленгам и интервалам времени между последовательным облучением приемных
пунктов. В [21] предложен способ определения дальности до ИРИ только по
измеренным интервалам времени между последовательным облучением приемных
пунктов. Однако, указанные выше методы представляют собой решение задачи
определения координат на плоскости и не позволяют определить координаты ИРИ в
пространстве.
Целью работы является разработка метода определения координат ИРИ со сканирующей ДНА на основе измеренных интервалов времени между последовательным облучением приемных пунктов.
Если в системе [21, 22], состоящей из трех приемных пунктов и сканирующего ИРИ, последовательно облучающего их (рис. 1), известны расстояния между приемными пунктами (^31, И12, И32) и углы между направлениями от ИРИ на приемные пункты (у31, у32, у12), которые вычисляются по времени задержки между последовательным облучением приемных пунктов при известной угловой скорости сканирования ИРИ, то местоположение сканирующего ИРИ можно определить как точку пересечения окружностей, описывающих возможное местоположение ИРИ относительно пары приемных пунктов при заданном угле между направлениями от ИРИ на приемные пункты (у31, у32, у12), с радиусами, рассчитываемыми по выражениям (1) - (3):
^31 -
^32 -
К12 -
к31
(/ ) , (1)
^32
2 (/32 ) , (2)
И12
(/12 ) . (3)
Рис. 1 - Определение координат сканирующего источника МПРЛС При известных координатах приемных пунктов ПП1 (х1, у1), ПП2 (х2, у2), ПП3 (х3, у3), координаты ИРИ определяются следующим образом:
1) Рассчитываются координаты центров окружностей {х12, у12}, {х32, у32} {х31, у31} путем решения системы уравнений
(ХППп ХОтп ) Уппп ^Отп ) ^тп ,
[(ХППт ~ ХОтп ) + (УППт ~ ^Отп ) = ^тп ,
(4)
где ХОтп, УОтп - искомые координаты центра соответствующей окружности
О
тп?
хППп, хППт, уППп, уППт- координаты соответствующих приемных пунктов; Ятп - радиус окружности.
2) Определяются координаты ИРИ решением системы уравнений
(* - % )2 +(у -У31 )2 = к,
(* - *12 )2 +(у - У12 )2 = ^
(* - *32 )2 +(У - У 32 )2 = я
2
32,
где х, у- искомые координаты ИРИ;
{х31, у31}, {х12, у12}, {х32, у32} - координаты центров окружностей.
Система уравнений (5) обладает избыточностью. Решение двух любых уравнений из системы (5) будет давать два решения, то есть две точки пересечения двух окружностей из трех. При этом первая точка определяет координаты ИРИ, а вторая точка пересечения будет являться точкой местоположения одного из приемных пунктов системы (рис. 1). По вычисленным координатам определяются наклонные дальности до ИРИ.
Вышеприведенный порядок определения координат ИРИ работоспособен лишь на плоскости. Для определения координат сканирующего ИРИ в пространстве требуется увеличение количества приемных пунктов. Рассмотрим один из вариантов (рис. 2) построения пространственно-распределенной системы, состоящей из пяти приемных пунктов и ИРИ воздушного базирования со сканирующей ДНА.
\
Рис. 2 - Вариант размещения пяти приемных позиций и источника радиоизлучения
в пространстве
Представленную систему можно условно разделить на две подсистемы. Первая включает в себя ПП1, 11112, ПП3, ИРИ, а вторая соответственно ПП1, 11114, ПП5, ИРИ. В каждой подсистеме возможно определить наклонные дальности способом, указанным выше [21], однако углы наклона плоскостей к горизонтальной поверхности можно определить лишь в системе в целом.
По аналогии с решением плоскостной задачи рассмотрим кривые, характеризующие возможное местоположение ИРИ в пространстве при заданных наклонных дальностях 11, 12, 14, для подсистемы ПП1 - ПП2 - ПП3 - ИРИ и 11, 13, 15 для подсистемы ПП1 - ПП4 - ПП5 - ИРИ (рис. 3). Они представляют собой полуокружности, проекции которых на горизонтальную плоскость дают прямые, перпендикулярные прямым, проходящим через точки местоположения приемных пунктов ПП1 - ПП2 - ПП3 и ПП1 - ПП4 - ПП5, соответственно. Точка пересечения полуокружностей определяет местоположение ИРИ в пространстве (рис. 3). Обозначим точки пересечения кривых возможного положения источника излучения с горизонтальной плоскостью как ИРИ1 (для подсистемы ПП1 - ПП4 - ПП5) и ИРИ2 (для подсистемы ПП1 - ПП2 - ПП3) (рис. 3). Точки местоположения ИРИ с максимально возможной высотой над горизонтальной поверхностью обозначим, как ИРИ4 и ИРИ3, соответственно. Проекцию точки пересечения кривых на горизонтальную плоскость точкой ИРИ5. Пересечения проекций кривых с прямыми, проходящими через точки ПП1 и ПП4 обозначим как Ц1, а ПП1 и ПП3 как Ц2.
Рис. 3 - Определение горизонтальной дальности ИРИ до приемного пункта Так как при известных наклонных дальностях в системе из пяти приемных пунктов появляется избыточность, то далее рассмотрим систему, образованную частями подсистем, указанных выше, и состоящую из ПП1, ПП3 от первой подсистемы и ПП1, ПП4 от второй подсистемы.
Вариант расположения такой системы, представленной в виде проекции на горизонтальную плоскость, изображен на рис. 4.
ИРИ1
(Х4; У4)
Ц1 Ж А
(Х7; У7) х ) а1 ПП4
(Х5; У5)
У!
Х Ц2 ПП1
(Х6; У6) (Х1; У1)
Рис. 4 - Геометрия размещения ПП и ИРИ, когда а2 < 90 < а1, р2 < 90 < р1
(а3 =180 - аь Р3 = 180 - Р1) Определим координаты точки ИРИ5 и проекции линий наклонных дальностей
1ГЬ ^^ 1Г3.
На первом этапе требуется определить координаты точек ИРИ1 (х4, у4) и ИРИ2 (х3, у3). Для этого составим две системы уравнений для треугольников ПП4-ПП1-ИРИ1 и ПП 1-ПП3-ИРИ2:
>/(Хз - Х1 )2 +(Уз - У1 )2 = 11 , >/(Х3 - Х2 )2 +(Уз - У2 )2 = 12 ;
(6)
х4 - Х1 )2 +(У4 - У1 )2 = 11 ,
у/(х4 - Х5 )2 +(У4 - У5 )2 = 13 •
(7)
Решим первую систему уравнений относительно неизвестных координат (х3
у3) точки ИРИ2.
<
-M2 +JM22 - 4M,M3 У3 =-2—-—2-, (8)
3 2M W
ГДе Mj = 4 [(x - X )2 + (y - y2 )2
M2 = -4 [У - У2 ) (X2 + У12 - Х' - ^^^ + l22 - к2 ) - 2X (X - X ) (У - У2 ) + 2У2 (X - X ^
M3 = 4(X - X2 )2 (X22 + У22 - ¡22 ) - 4X2 (X - X )(X2 + У12 - *22 - У22 + l22 - О +
+(Xj2 + У12 - X22 - У22 +122 - lj2)2.
Найденное значение y3 может принимать два значения, положительное и отрицательное. С учетом известного направления прихода сигнала от ИРИ отрицательное значение y3 далее не учитывается.
Значение x3 определяется как
_ 2x2 ± ^4x¡ - 4( x22 + У32 - 2 Уз y2 + y22 -122) X - 2 , (9)
найденное значение x3 может принимать два значения, выбор того или иного выполняется аналогично y3.
Таким же способом определяются координаты точки ИРИ1 (x4, y4).
На втором этапе требуется определить координаты точек Ц2 (x6, y6) и Ц1(х7, y7). Из рис. 4 видно, что x6= x3, аy6=y^y2.
Для определения координат точки Ц1 составим систему уравнений
>/(X- X )2 + (У7- У4 )2 - sin (а3 ) l3 >
V(X - X )2 + (У7 - У5 )2 - C0s (а3 ) l3 .
(10)
Решив данную систему, получим координаты х7, у7. Найденное значение у7
может принимать два значения. Для выбора одного из двух решений используем
условия: у7>у5, х7<х5.
На третьем этапе определяются координаты точки ИРИ5.
Составим систему уравнений
| Ах8 + + С1 = О,
I А X + В2 У8 + С2 = 0, ( )
где А1 = У4 - Уб , В1 = Хб - Х4 , С1 = Х4У6 - Х6У4 ; А = У3 - У7 , В2 = Х7 - Х3 , С2 = - ^3 • Решив систему уравнений, получим:
= А^1 - ^ (12)
А1В2 - А2 В1, ' '
X = -КУ - С1 . (13) А1
На четвертом этапе определяется горизонтальная дальность до любого из приемных пунктов, например:
1Г3 =>/(Х5 - Х8 )2 +(У5 - У8 )2 . (14)
При этом высота
Н = ^/32 - /г32 . (15) Разработанный метод позволяет определить координаты ИРИ со сканирующей ДНА по каждому одиночному циклу обзора.
Библиографический список
1. Кобан А.Я., Самотонин Д.К Шучно-технические проблемы развития федеральной системы разведки и контроля воздушного пространства Российской Федерации и пути их решения // Военная мысль. 2017. № 4. С. 14 - 1S.
2. Астапенко Ю., Гуменюк А., Менячихин А. РТВ надо иметь дополнительные средства разведки // Воздушно-космический рубеж. 201S. № 1(3). С. 3S - 42.
3. Ильин Е.М., Климов А.Э., Пащин H.C, Полубехин А.И., Черевко А.Г., Шумский В.К Пассивные локационные системы. Перспективы и решения // Вестник СибГУТИ. 2015. № 2. C. 7 - 20.
4. Канащенков А.И., Меркулов В.И., Самарин О.Ф. Живучесть радиолокационных систем // Радиотехника. 2006. № 11. С. 5 - 10.
5. Ашурков И.С., Какаев В.В., Лешко HA. Оптимизация пространственной структуры многопозиционной радиолокационной системы // Информационно-управляющие системы. 2015. № 6 (79). С. 81 - S5.
6. Ашурков И.С., Лешко HA. Расчет рабочей зоны многопозиционной радиолокационной системы по стороннему источнику подсвета // Труды МАИ. 2015. № 83. URL: http://trudymai.ru/published.php?ID=62297
7. Bergin J., Guerci J.R. MIMO Radar: Theory and Application. Boston, London, Artech House, 201S, 2S0 p.
S. Wang R., Deng Y. Bistatic SAR System and Signal Processing Technology, Springer, 201S, 2S6 p.
9. Griffiths H.D., Baker C.J. An Introduction to Passive Radar, New York, Artech House, 2017, 110 p.
10. Aubain A., Lobert B. Patent 6 417 799 US. Method of locating an interfering transmitter for a satellite telecommunications system, Int. Cl. H 04 D 7/19, 9.07.2002.
11. Webber J.C., Knight C.A. Patent 5 594 452 US. Method and System for Locating an Unknown Transmitter Using Calibrated Oscillator Phases, Int. Cl. H 04 B 7/185, 14.01.1997.
12. Dillman D.J. Patent 6 839 017 US. Method and apparatus for using anti-jam technology to determine the location of an electromagnetic radiation source, Int. Cl. G 01 S 7/36, 4.01.2005.
13. Aubain A. and Lobert B. Patent no. 6417799 US. Method of locating an interfering transmitter for a satellite telecommunications system, 2002.
14. N.J. Willis. Bistatic Radar, 2nd Ed., Silver Spring, MD: Technology Service Corp., 1995. Corrected and republished by Raleigh, NC: SciTech Publishing, Inc., 2005, 337 p.
15. Борисов Е.Г. Определение местоположения воздушных объектов в полистатической радиолокационной системе, паразитирующей на излучении телекоммуникационных систем // Научный вестник московского государственного технического университета гражданской авиации. 2018. Т. 21. № 5. С. 105 - 116.
16. Кирюшкин В.В., Волков Н.С. Межпозиционное отождествление результатов измерений и определение координат воздушных целей в многопозиционной радиолокационной системе на беспилотных летательных аппаратах // Теория и техника радиосвязи. 2019. № 1. С. 107 - 116.
17. Меркулов В.И., Садовский П.А. Оценивание дальности и ее производных в
двухпозиционной пассивной радиолокационной системе // Труды СПИИРАН. 2018. № 1 (56). С. 122 - 143.
18. Ефимов Е.Н., Шевгунов Т.Я. Формирование оценки направления прихода сигнала с использованием искусственных нейронных сетей // Труды МАИ. 2015. № 82. URL: http://trudymai.ru/published.php?ID=58786
19. Анненкова И.Ю., Грибанов А.С. Ошибки линий положения для различных методов определения направления на источник излучения // Труды МАИ. 2011. № 42. URL: http://trudymai.ru/published.php?ID=24260
20. Мельников Ю.П., Попов С.В. Радиотехническая разведка. Методы оценки местоопределения источников излучения. - М.: Радиотехника, 2008. 432 с.
21. Петроченков Д.М., Федотов А.В. Временной метод определения дальности до сканирующего источника радиоизлучения в многопозиционной радиолокационной системе с некооперируемым подсветом // Вестник Ярославского высшего военного училища противовоздушной обороны. 2018. № 2 (3). С. 78 - 84.
22. Петроченков Д.М., Федотов А.В. Организация обзора пространства в многопозиционной радиолокационной системе с некооперируемым источником подсвета на основе комбинаторного принципа // Журнал Сибирского федерального университета. Серия: техника и технологии. 2018. Т. 11. № 7. С. 831 - 841.