Journal of Siberian Federal University. Engineering & Technologies, 2018, 11(3), 325-335
yflK 621.396.96
Method of Air Objects Status Flag Forming from Aggregate Radar Range Portraits with a Boardhand Sensing
Igor V. Sisigin,
Dmitriy A. Ravdin and Konstantin O. Kolesnikov*
Yaroslavl Higher Military Air Defence College 28 Moskovskiy, Yaroslavl, 150001, Russia
Received 12.12.2017, received in revised form 04.01.2018, accepted 27.02.2018
The article reveals the relevance of the assessment of air object in wideband sensing, suggestion the method of classification of the condition of object on the dynamics of their radar range portraits, described the method of forming aggregate radar range portraits from collapsing and maneuvering air targets with a boardband sensing, presented the results of mathematical simulation of collapsing and maneuvering aircrafts, proposed the structural scheme of the system of assessment of aerial objects on a correlation basis.
Keywords: modeling, detection of a condition of object, correlation coefficient, radar range portrait, wavelet transform.
Citation: Sisigin I.V., Ravdin D.A., Kolesnikov K.O. The method of forming aggregate radar range portraits from collapsing and maneuvering air targets with a boardband sensing, J. Sib. Fed. Univ. Eng. technol., 2018, 11(3), 325-335. DOI: 10.17516/1999-494X-0045.
© Siberian Federal University. All rights reserved
Corresponding author E-mail address: [email protected], [email protected]
Методика формирования признака состояния воздушных объектов по совокупности радиолокационных дальностных портретов при широкополосном зондировании
И.В. Сисигин, Д.А. Равдин, К.О. Колесников
Ярославское высшее военное училище противовоздушной обороны Россия, 150001, Ярославль, пр. Московский, 28
В статье раскрывается актуальность оценки состояния воздушных объектов при широкополосном зондировании, предложен вариант классификации состояния воздушных объектов по динамике их радиолокационных дальностных портретов, описана методика формирования совокупности радиолокационных дальностных портретов маневрирующих и разрушающихся воздушных объектов при широкополосном зондировании, представлены результаты математического моделирования маневрирующих и разрушающихся воздушных объектов.
Ключевые слова: моделирование, распознавание состояния воздушного объекта, коэффициент корреляции, радиолокационный дальностный портрет, вейвлет-преобразование.
В [1, 2] говорится, что наряду с традиционными (обнаружение и разрешение целей, измерение координат и параметров их движения) в число основных задач радиолокационного наблюдения выдвинулось распознавание классов, типов и функционального состояния воздушных объектов (ВО).
В [3-5] достаточно широко представлены методы и алгоритмы распознавания классов и типов ВО. В [6] предложен метод распознавания маневра ВО, но не рассматривается процесс их разрушения. В [7, 8] автор предлагает способ определения момента фрагментации объекта, но рассматривает только процесс кинетического взаимодействия космических аппаратов с астероидами и другим объектами. Однако оценка функционального состояния аэродинамических ВО радиолокационными методами остается нерешенной, поэтому актуальность рассматриваемой темы достаточно очевидна.
Под функциональным состоянием воздушного объекта понимается его положение в пространстве по отношению к радиолокационной станции (РЛС) (прямолинейное движение или совершение какого-либо маневра) и геометрическое расположение его составных частей по отношению друг к другу (разрушение или разделение), которое может быть определено по отраженному сигналу.
Как известно, отраженный от цели сигнал при высоком разрешении имеет достаточно сложную временную структуру, которая определяется геометрическими и электродинамическими особенностями. Такой сигнал после согласованной обработки и амплитудного детектирования дает радиолокационный дальностный портрет (РЛДП). Реализация РЛДП зависит не только от геометрической формы воздушного объекта (ВО), но и от ориентации его элементов
относительно лкниивизирнваниаРЛС.Поаоомунамивыдвинута гипотеза о том, что РЛДП ВО как при маневра твк и при вморушэнир нрутерревает сущэственныз изментиия,одиако при маневрт они моэуо тосиверегрлярныРхарскрерве, приразрушенит-нииричный. Анавиз динамики РРЩПпозвонит ккиссифицнроунсь ВОно срвозянию: пионолинейноздувжениеЗманевр/ разрушение (рис. 1).
Для проверки выдвинутой гипотезы были рассмотрены три возможных состояния ВО: прямолинейное движение, маневр и разрушение.
Для анализа был выбран метод математического моделирования, основанный на преобразовании SD-моделей ВО в дифференциальную эффективную поверхность рассеяния (ЭПР), свертки еес зондирующим сигналомипоследующейсогласованной обработки с вычислением модуля выходного аигеала (рис. 2), которыа подробно описан в [9].
Динамика РЛДП оценивалась по признаку в виде коэффициента корреляции (КК) между текущим и предыдущим РЛДП в соответствии с выражением
г _ Е0рс-*)°вс-у) (1)
ху V2X*fc-*)2Z(yfc-;y)2'
где хк -РЛДПВО в мемснторемсои t; X - среднее значение x по выборке X; yk - РЛДП ВО в молент времени 2 У - среднеезначение y повыборке Y.
Рис. 1. Классификация состояния воздушныхобъектов иописание признаков на основеРЛДП
Fig. 1. Classification of the condition of the objects and a description of the characteristics based on radar longrange portraits(RLRP)
Создание трехмерного
объекта исследования, экспорт модели в формат *.ase
Анализ и считывание данных из *.ase файла в Matlab
Расположение
модели в пространстве, удаление невидимътх примитивов
Рис. 2. Этапы математического моделирования РЛДП Fig. 2. Stages of RLRP mathematical modeling
Преобразование фацетной модели в виде элементарных излучателей
Для проведения эксперимента были разработаны модели процессов разрушения и маневра типового воздушного объекта.
Процесс разрушения ВО в результате кинетического воздействия либо технической неисправности является случайным, поэтому для его анализа были использованы видеоматериалы реальных событий с ВО. В качестве прототипа варианта разрушения был взят видеофайл, отдельные кадры которого представлены на рис. 3.
Время анализа ограничено одной секундой, так как за этот промежуток происходят существенные изменения геометрической формы ВО. На основе имеющихся видеоматериалов в программе 3Dmax разработаны покадровые модели ВО, соответствующие по характеру изменения геометрической формы объекту из видеофайла (рис. 4).
Формирование покадровой последовательности из 25 РЛДП проводилось описанным выше методом математического моделирования [9].
На рис. 5 и 6 представлены совокупности РЛДП разрушающегося ВО №1 при сигналах с шириной спектра (ШС) 30 и 300 МГц соответственно.
Для представленных на рис. 5, 6 совокупностей РЛДП по формуле (1) рассчитаны значения КК. Минимальное значение КК при сигнале с ШС 30 МГц составило 0,95, при сигнале с
I II II -I
0 кадр (0,0 с) 5 кадр (0,2 с) 10 кадр (0,4 с)
\ А
15 кадр (0,6 с) 20 кадр (0,8 с) 24 кадр (1,0 с)
Рис.3. Отдельные кадры процесса разрушенияВО Fig.3. Individualframes oftheprocessofdestructionofaircraft object (AO)
и .и. я
) шлр 0 кадр 10 кадр
.9. ч —. x \V 0 1
15 кадр 20 кадр 24 кадр
Рис. 4. 3D-модели разрушения воздушного объекта №1 Fig. 4. 3D-models of destruction of air object №1
Рис. 5. Совокупность РЛДП разрушающегося 130 № 1при сигнале с ШС 30 МГц
Fig. 5. Set of RLRP of the collapsing air object №1 at the sign al withaspectrumwidth of30 M Hz
10 15 20
Протяженность РЛДП, м
Рис. 6. Совокупность РЛДП разрушающегосяВОаЧ°1при сигнале с ШС 3 00 МГц
Fig. 6. Set of RLRP оf the collapsing tir object №1 at the signal itvitli я spectrum width of 300 MHz
ШС 300 МГц - 0,72. Для получения д остаточной статистики былореализовано несколько вариантов разрушения четы рех типов В О с различных раку рсов . Во всех экспе риментах рузульта-ты имели достаточную сходимост ь: ср еднее значение КК при сигнале с ШС 30 МГц составило не менее 0,9, а для сигнала с ШС300 МГц - не менее 0,7.
Маневр воздушного объекта - это процесс, близкий к детерминированному. Поэтому его можно представить в виде движения поокружности радиусом R (рис. 7).
В случае пилотируемого полета предельное значение длительной перегрузки не может превышать значения 10 g, что соответствует пределу физиологических возможностей человека.
При этом связь между параметрами движенияирадиусомкривизны траектории R определяется выражением
ач = Т= ^ (2)
где Уц - тангенциальная составляющая скорости объекта в точке начала маневра; - угловая скорость движения объекта по окружности; R — радиус кривизны траектории движения объекта.
Начало манёвра
Равномерное прямолинейное движение
Рис. 7. Модель выполнения процесс а маневра ВО Fig. 7. The execution model of the AO maneuver process
4У» 4
R О
Таблица 1 Table 1
Уц, м/с 100 200 300 400 500 600 700 800 900 1000
R,m 102 408 918 1632 2551 3673 5000 6530 8265 10204
Сучееом пнрегрнзки n - отношения абсолютнойвеличины линейного ускорения цели ац вызванного нелравитацивннвшл силами, к ускорению свободнрго падение g, имеем:
ац
П=—,
пд =~t = R-
(3)
(4)
Таким образом, угловая скорость вращения относительно РЛС, которая приводит к изменению РЛДП, зависит от тангенциальной составляющей скорости цели и допустимых значений перегрузки. Рассчитанные радиусы кривизны траектории маневра для скоростей ВО от 100 м/с до 1000 м/с представлены в таблице 1.
Дальнейшие результаты представлены для модели маневра при скорости движения ВО в 300 м/с. Исходя из этого скорость ю будет равна 0,33 рад/с, что соответствует повороту ЛА относительно РЛС со скоростью 19°/с.
На рис. 8 и 9 показаны совокупности РЛДП маневрирующего ВО №1 при углах курса в=0-19° с шагом 1° и ширине спектра зондирующего сигнала 30 и 300 МГц соответственно.
Для представленных на рис. 8 и 9 совокупностей РЛДП по формуле (2) рассчитаны значения коэффициентов корреляции между текущим и предыдущим реализациями РЛДП. Все значения коэффициентов корреляции лежат в пределах значений 0,9-0,99.
На рис. 10 показаны графики изменения коэффициентов корреляции ВО №1 для сигналов с ШС 30 и 300 МГц для маневра и разрушения.
На рис. 10 видно, что КК РЛДП не опускается ниже уровня 0,7 при сигнале с ШС 300 МГц, а при сигнале с ШС 30 МГц - не ниже 0,95.
Отношение коэффициентов корреляции при маневре и разрушении составили 1,04 и 1,36 для 30 и 300 МГц соответственно, что свидетельствует о малой контрастности признака.
Рис. 8. Совокупность РЛДП маневрирующего ВО №1 при углах курса в=0°-19° для сигнала с ШС ширине спектра 30 МГц
Fig. 8. Set of RLDP for maneuvering AO №1 at angles of в=0°-19° for signal spectrum width spectrum width of 30 MHz
ffl
3 2,5
[ 2
:
ч и
de 1
0,5 0
iMiPL
-------------
10 15 20
Протяженность РЛДП, м
25
30
Рис. 9. Совокупность РЛДП маневрирующего ВО №1 при углах курса в=0-19° для сигнала с ШС ширине спектра 300 МГц
Fig.B.SetofRLDP for maneuvaoing АО №1 at angles of /-52°-19° for signal spectrum width spectrum width of 300 MHz
Таким образом, гипотеза о существенном изменении РЛДП при маневре и разрушении подтвердилась, однако значения выбранного признака оказались невелики, поэтому для повышения контрастности выбранного признака (КК) целесообразно провести дополнительную обра ботку.
Для решения этой зодачи былая исоооьзооана двухэтапная процедура вейвлет-преобразования, предложенная в [5] для формирования контрастных признаков распознавания рэдиосякэбибнных целей.
На первом этапе для каждого РЛДП вычисляются вейвлет-коэффициенты Щ,(а, Ь) в соответствии с выражением
И£(а,Ь) = -^Г^Ж^Ж (4)
1.05 1
0.95 0.9 0.85 0.8 0.75 0.7
1 1 1
V / А N /\ ^' / к / / 1 | i
/ / I
\ /
\ / \ / \ /
\/ " V I ! I
РЛДП Су-34 маневр 30 МГц
РЛДП Су-34 маневр 300 МГц
РЛДП Су-34 разрушение 30 МГц
РЛДП Су-34 разрушение 300 МГц
Рис. 10. Графикизменения КК совокупности РЛДПВО №1при сиг нале сШС 30 и 300 МГц для маневра и разрушения
Fig. 10. Graph of changes in the correlation coefficients ofthe set of RLRP of the air object №1 at a signal with a spectrum width of 30 MHz and 300 MHz for maneuvering and destruction
где7(и-вх°дн°й РЛДП; р - параметр масштаба; Ь - параметр временного переноса; ^ -
базисная функция Хаара.
На втором этапе из Ш2(а, Ь) формируются LLE (а, Ь) - линии локальных экстремумов (ЛЛЭ) всоответствиисвыражением
Л | ^(М) | ^
ИЕ(а,Ь) = Г Ь) , (5)
V 1 йЬ 1
где ДШ-заданноепороговоезначениескорости изменения функции Wz(a, Ь).
На рис. 11 представлены совокупности ЛЛЭ ВО №1 при сигнале с ШС 30 МГц для маневра и разрушения.
Характерной особенностью ЛЛЭ является их компактность и малая дисперсия при регулярном изменении РЛДП (рис. 11а). При разрушении цели ЛЛЭ изменяются значительно больше (рис. 116).
Для оценки динамики коэффициента корреляции ЛЛЭ при маневре и разрушении был проведен статистический эксперимент, результаты которого изображены на рис. 12 и 13.
Из графиков на рис. 12, 13 видно, что среднее значение КК^Е при маневре более чем в два раза превышает значение КК11Епри разрушении для сигналов с ШС 30 и 300 МГц, а среднеквадратичное отклонение ККЬЬЕ при маневре в 1,5 раза меньше, чем при разрушении для сигнала с ШС 30 МГц и в 3 раза меньше для сигнала с ШС 300 МГц (0,13 / 0,19; 0,06 / 0,2 соответственно).
Таким образом, предложенная на рис. 14 методика формирования признака маневра/разрушения, основанная на вычислении линий локальных экстремумов от вейвлет-коэффициентов РЛДП, обеспечила его высокую контрастность, что может повысить качество распознавания в задаче оценки функционального состояния объектов.
Рис. 11. Совокупность ЛЛЭ ВО№1: а - приманевре (б=0-19°);б - при разрушении
Fig. 11. Setof lines the local extrema (LLE):a - whenmaneuvering(yS=0-19°);b)with thedestruction
A = 0
0.656 0.6 r
J L„v
A= 1
0.128 0
/ / _____ ---- \ / ______r
—N
/ 1 ч„ /
V
i \ / . i i X_____ /
\ / \ \
\
РЛДП BO №1 маневр, 30 МГц
РЛДП ВО № I
разрушение, 30 МГц
ЛЛЭ ВО №1 маневр, 30 МГц
ЛЛЭ ВО №1 разрушение, 30 МГц
12 к
20
22
24"
Рис. 12. График измене нияКК совокупностиРЛДПиЛЛЭ д^я^^ангад)^ и разрушенияприси гнале с ШС 30 МГц
Fig. 12. Graph of changes incorrelation coefficientsofthesetof RLRPandLLEformaneuvering end destruction at a signal wihda tpectiom wddthat th MHz
Рис. 13. График изменения KKcовокупности РЛДП и ЛЛЭ ВО №1 для маневра и р азрушения при сигнале с ШС 300 МРц
Fig. 13. Graph of time correlation coefficients of the set of- RLRP and LLE to maneuver anddestruction with a spectrum widih of 30 MHz
Рис. 14. Структурная схема методики формированияпризнакаманевра/разрушения Fig. 14. Structural scheme of the formation method of maneuver/destruction sign
Список литературы
[1] Слюсарь Н.М. Эффект вторичной модуляции радиолокационных сигналов: физические основы и практическое применение, Вестник Военной академии Республики Беларусь, 2003, 1, 61-77. [Slyusar N.M. Effect of secondary modulation of radar signals: physical basis and practical application, Bulletin of the military Academy oftheRepublicof Belarus, 2003. 1. 61-77. (in Belarus)]
[2] Слюсарь Н.М. Радиолокационное распознавание, контроль состояний и действий целей в радиотехнических системах вооружения на основе вторичной модуляции отраженных сигналов. Минск: ВА РБ, 2003. 144. [Slyusar N.M. Radar detection, control of States and actions goals, in electronic weapons systems on the basis of secondary modulation of the reflected signals. Minsk: VA RB,2003. 144.(in Belarus)]
[3] Богданов А.В., Бондарев В.Н., Васильев О.В., Гарин Е.Н. и др. Синтез оптимальных алгоритмов распознавания групповых воздушных целей в радиолокационных системах воздушного базирования. Журнал Сибирского федерального ун-та. Техника и технологии 2017, 10(2). 155-168. [Bogdanov A.V., Bondarev V.N., Vasiliev O.V., Garin E.N. and other. Synthesis of optimal algorithms for recognition of group air targets in air-based radar systems. Journal Siberian Federal University , Engineering and Technology, 2017, 10(2). 155-168. (inRussia)]
[4] Ищук И.Н., Степанов Е.А., Бебенин А.А. и др. Способ классификации объектов оптико-электронными системами разведки на основе обработки многоспектрального кубоида изображений. Журнал Сибирского федерального ун-та. Техника и технологии, 2017, 10(2), 183-190. [Ischuk I.N., Stepanov E.A., Bebenin A.A. etc. Method of classification of objects by electro-optical
intelligence systems based on the processing of multispectral cuboid images. Journal Siberian Federal University, Engineering and Technology, 2017, 10(2), 183-190. (in Russia)]
[5] Колесников К.О., Сисигин И.В., Равдин Д.А. и др. Методика формирования признаков распознавания на основе вейвлет-преобразования радиолокационных дальностных портретов. Журнал Сибирского федерального ун-та. Техника и технологии, 2017, 10(2), 277-285. [Kolesnikov K.O., Sisigin I.V., Ravdin D.A. and others. Methods of forming the signs of recognition based on wavelet transform radar dolnestic portraits. Journal Siberian Federal University, Engineering and Technology, 2017, 10(2), 277-285. (in Russia)]
[6] Ковтунов А.Л. Лещенко С.П., Батуринський М.П., Польшина Л.В. Метод обнаружения маневра радиолокационной цели в обзорных РЛС с использованием сверхширокополосных сигналов. Збiрник наукових праць Харювського унiверситету Повтряных Сил, 2013, 3(36). 5559. [Kovtunov A.L. Leshchenko S.P., Baturinsky M.P., Pollina L.V. Method of maneuver detection in radar target surveillance radar using ultra-wideband signals. Collection of scientific papers of Kharkiv University Air Forces, 2013, 3(36). 55-59. (in Ukraine)]
[7] Антропова С.М. Моделирование одного из подходов к решению задачи оценки факта поражения. Сборник докладов НТК «ВВТВКО» в МГТУ имени Н.Э. Баумана, 2014. [Antropova S.M. Modeling of one of the approaches to the problem of assessment of the fact of defeat. Collection of reports of STC "BBT EKR" in MSTU named After N. Uh... Bauman, 2014 (in Russia)]
[8] Антропова С.М. Исследования фрагментации объекта для решения задач контроля факта поражения. Сборник докладов 5-й НТК «Актуальные вопросы развития систем и средств ВКО», 2015. [Antropova S.M. Study of the fragmentation of the object for the decision of tasks of control of the fact of defeat. Collection ofpapers 5-th NTK "Topical issues of development of systems and tools, ASD", 2015. (in Russia)]
[9] Колесников К.О. Сисигин И.В. Равдин Д.А., Сушин Ю.В. Методика моделирования отраженных сигналов от сложных объектов на основе приближения Гюйгенса-Френеля. Вестник воздушно-космической обороны, 2014. [Kolesnikov K.O., Sisigin I.V., Ravdin D.A., Sushin Yu.V. Modeling reflected signals from complex objects on the basis of approximation of the Huygens-Fresnel. Journal of Aerospace Defense, 2014 (in Russia)]