Радиолокация и радионавигация
УДК 621.396.96
Р. С. Фадеев, А. В. Мякиньков, В. Н. Буров, А. Г. Огурцов Нижегородский государственный технический университет им. Р. Е. Алексеева
Возможности обнаружения и определения координат малозаметных целей в многопозиционных радиолокационных системах с размещением позиций на борту беспилотных летательных аппаратов1
Рассчитаны просветные эффективные площади рассеяния некоторых широко распространенных моделей беспилотных летательных аппаратов (БЛА), а также зона обнаружения просветной радиолокационной станции (РЛС) при работе по целям указанного типа. Предложена конфигурация просветной многопозиционной РЛС, повышающая эффективность обнаружения малозаметных целей. Построена общая структурная схема системы обработки сигналов. Получены оценки потенциальной точности определения координат целей в многопозиционных РЛС с размещением позиций на борту БЛА.
Многопозиционная просветная радиолокационная система, подвижные позиции, беспилотный летательный аппарат, зона обнаружения, компенсация пассивных помех, многоканальный когерентный накопитель, метод максимального правдоподобия, информационная матрица Фишера
В последнее время все более широкое применение для решения задач воздушной разведки и нанесения ударов с воздуха находят беспилотные летательные аппараты (БЛА). БЛА, применяемые для разведки, имеют малые габариты; БЛА, используемые для нанесения высокоточных ракетных ударов, обладают большими размерами (порядка 10 м), но изготавливаются как правило, по технологии "Стелс". Для изготовления корпуса БЛА часто используются композитные материалы с высоким коэффициентом поглощения. Из литературных источников известно, что эффективная площадь рассеяния (ЭПР) современных БЛА со-
2
ставляет порядка 0.001 м [1], вследствие чего их обнаружение наземными моностатическими радиолокационными станциями (РЛС) затруднено. С другой стороны, использование для их обнаружения крупногабаритных РЛС с мощными передатчиками, обеспечивающими приемлемую зону обнаружения, может оказаться затруднительным при необходимости оперативного изменения зоны действия.
Эффективность обнаружения малозаметных целей может быть обеспечена при использовании
1 Работа выполнена при поддержке Министерства образования и (Постановление Правительства РФ от 9 апреля 2010 г. № 218).
© Фадеев Р. С., Мякиньков А. В., Буров В. Н., Огурцов А. Г.,
принципов просветной радиолокации за счет эффекта резкого увеличения эффективной площади рассеяния (ЭПР) при бистатическом угле (угле между линиями, соединяющими цель с передающей и приемной позициями), близком к 180°. При такой конфигурации просветная ЭПР цели не зависит от наличия радиопоглощающего покрытия, а определяется только теневым силуэтом цели [2], поэтому РЛС с обнаружением "на просвет" позволяют эффективно обнаруживать и измерять координаты малоразмерных целей, в том числе имеющих указанное покрытие [3]. Расширение зоны обнаружения просветной РЛС может быть достигнуто использованием многопозиционной конфигурации [2], при этом также возможно повышение разрешающей способности и точности измерения координат целей [2].
Одной из актуальных задач при проектировании является обеспечение быстрого развертывания радиолокационного комплекса в том месте, где его применение необходимо в данный момент. Одним из вариантов решения указанной задачи является размещение приемных и передающих позиций на
науки РФ в рамках договора № 02.G25.31.0061 от 12 февраля 2013 г.
2014
29
борту БЛА, которые, передвигаясь в заданной области воздушного пространства, обеспечивали бы зону обнаружения необходимых размера и конфигурации. При этом вследствие малой ЭПР БЛА достигается хорошая скрытность самой системы. Кроме того, такое построение многопозиционной РЛС позволяет наилучшим образом использовать энергетический потенциал системы за счет выбора оптимального угла облучения целей.
Эффективные площади рассеяния беспилотных летательных аппаратов. На рис. 1 приведены контуры некоторых известных моделей БЛА: "MQ-1 Predator" (а), "RQ-2 Pioneer" (б), "Boeing X48B" (в) и "ZALA 42116" (г), а на рис. 2 -результаты расчета просветной ЭПР этих целей для пары ортогональных ракурсов (углов обзора). Контуры 1 (вид "сбоку") соответствуют просвет-ной РЛС с расположением позиций на земле, контуры 2 (вид "сверху") - при размещении передающей позиции в воздухе, а приемной - на земле. Зависимости просветных ЭПР от угловой координаты (азимута) ß для этих же БЛА при движении цели поперек линии базы представлены на рис. 2. Расчеты проведены для длины волны Х = 0.13 м методом аппроксимации, описанным в [3]-[5].
Кривые 1, 3 на рис. 2 соответствуют контуру 1 на рис 1, кривые 2, 4 - контуру 2. Кривые 1, 2 получены для минимального расстояния между траекторией цели и линией базы Дгтщ = 0, кривые 3, 4 - для Armin = 1 км.
При Агт^д = 0 (траектория цели пересекает линию базы) просветная ЭПР цели при обзоре "сверху" значительно больше, чем при обзоре "сбоку". Кроме того, ширина основного лепестка больше для контура 2 (см. рис. 1). Это объясняется тем, что размер крыльев цели в направлении координаты х (см. рис. 1) намного меньше, чем в направлении координаты у, а ширина главного лепестка определяется в основном размером по координате х. Таким образом, при обзоре "сверху" наблюдается сумма доминирующего широкого основного лепестка просветной ЭПР крыльев и узких основного и боковых лепестков просветной ЭПР фюзеляжа. В результате формируется достаточно широкий главный лепесток с колебаниями, которые определяются лепестками диаграммы рассеяния контура фюзеляжа.
При смещении относительно линии базы усредненное ЭПР уменьшается, причем оно становится практически одинаковым для обоих направлений обзора, т. е. уменьшение ЭПР гораздо сильнее для обзора "сверху". Это объясняется узостью основного лепестка диаграммы прямого рассеяния цели при облучении "сверху", которая обусловлена большей площадью контура 2 (см. рис. 1).
Учитывая полученные результаты, можно предложить геометрию многопозиционной просветной РЛС (рис. 3). В этой конфигурации воздушные приемные позиции Пр1, Пр2 используются сов-
У, м
12-
4-
4.1
а
а
Д'
дб/м~ 40 25 10 -5
У, м 3.9 2.6 1.3
.3- 1__h_
0 2 3 г
а
2.3
б
У, м 4.8
0 1 2 х, м в
Рис. 1
а
ц'
дб/м~ 35 25 15 5
д
Ч 2 \
1 'Ч
-7.0 -3.5 0 3.5 ß, ..
а
-7.0 -3.5 0 3.5 ß, -7.0 -3.5 0 3.5 ß,
б в Рис. 2
У, м
а
ц
дБ/м 22 12
2
0.6
г
2
....
J_I_L
-7.0 -3.5 0 3.5 ß, .
г
8
2
1
0
х, м
х, м
4
Л—
Пд
/3 '.........
/
^ /4 /3
Пр1
Поверхность Земли
Рис. 3
местно с наземной приемной позицией Пр3. Подсвет цели Ц производится воздушной передающей позицией Пд. На рис. 3 указаны также расстояния Ь - Ьз между Пд и Пр1-Пр3. Сплошными линиями показано распространение электромагнитной энергии от Пд до Ц и далее распространение отраженной энергии к Пр1-Пр3. Высота полета воздушных позиций должна быть не меньше высоты полета целей. Предложенная конфигурация обладает рядом преимуществ перед конфигурацией с наземным размещением позиций. Во-первых, в ней возможно облучение цели при выгодном ракурсе ("сверху"), во-вторых, поскольку лучи диаграмм направленности наземных приемных антенн ориентированы на воздушные позиции, влияние отражений от поверхности Земли минимально.
При размещении позиций на подвижных носителях спектр пассивной помехи (ПП) от местных предметов и подстилающей поверхности сдвигается за счет взаимного перемещения приемной и передающей антенн, причем для разных участков отражающей поверхности этот сдвиг оказывается различным. Вследствие ограниченности геометрических размеров БЛА на Пр1 и Пр2 предполагается использование слабонаправленных антенн, поэтому измерение угловых координат целей (азимута и угла места) не производится, а измеряемыми параметрами являются временные задержки сигналов. Каждая пара, включающая передающую и одну из приемных позиций, характеризуется пространственными элементами разреше-
ния в виде эллипсоидов равной суммарной дальности. На поверхности Земли эти эллипсоиды отображаются эллипсами равной дальности (см. рис. 3, 2), пересечение которых определяют элементы разрешения системы 1. Лучи /1, /{ — /4, /4 отображают формирование мешающего сигнала отражения от подстилающей поверхности. Антенна наземной позиции Пр3 направлена в сторону передающей позиции, в связи с чем диаграмма направленности имеет положительный угол места. При этом влияние поверхности будет существенно меньше, чем для случая воздушных приемных позиций, и влиянием подстилающей позиции на наземную приемную позицию можно пренебречь.
Важно отметить, что различным участкам поверхности соответствует различная задержка отраженного сигнала. Таким образом, для компенсации ПП необходимо, чтобы РЛС обладала разрешающей способностью по задержке, для чего следует использовать сигналы с широким спектром. Анализ показывает, что ширина полосы сигнала для обеспечения необходимого разрешения по задержке при длине линии базы порядка 10 км и длине волны около 0.5 м должна составлять более 100 МГц, т. е. должны применяться широкополосные или сверхширокополосные зондирующие сигналы. При этом размер элемента разрешения подстилающей поверхности будет составлять, приблизительно, 100*100 м.
Принципы построения системы обработки сигналов. Рассмотрим общие принципы постро-
Ь
2
4
Рис. 4
ения системы обработки широкополосных импульсных сигналов на приемной позиции биста-тической (просветной) РЛС. Общая структурная схема системы обработки изображена на рис. 4.
Сигнал промежуточной частоты /Пч с выхода супергетеродинного приемника поступает на первый преобразователь с понижением несущей частоты (ППНЧ1). В ППНЧ1 спектр сигнала переносится на нулевую частоту. Однако вследствие доплеровского смещения частоты, обусловленного взаимным движением позиций, и нестабильности частоты гетеродина /г, к тому же не точно совпадающей с центральной частотой сигнала /Пч , прямой сигнал на выходе ППНЧ1 имеет центральную частоту А/, отличную от нуля.
Необходимо оценить частоту А/ и дополнительно перенести сигнал на нулевую частоту. Возможны два варианта: использование ППНЧ2 (ход сигнала показан штрихпунктирными линиями) или обратной связи для подстройки /г в ППНЧ1 (система частотной автоподстройки частоты) (ход сигнала показан штриховыми линиями).
используется режекторный фильтр РФг- подавления ПП в элементе разрешения, в каждый момент времени настроенный таким образом, что центральная частота зоны режекции /Пи соответствует текущим параметрам движения позиций, а ширина зоны режекции А/Пп определяется разбросом доплеровских частот ПП в элементе разрешения (рис. 6). В рассматриваемом диапазоне длин волн при скоростях движения позиций до 200 м/с ширина спектра помехи в элементе разрешения может составлять до 300 Гц.
После РФ включены когерентный многоканальный (по частоте Доплера рассеянного целью сигнала) накопитель пачки и схема обнаружения цели. Сигналы от каналов обнаружения поступают в схему вторичной обработки. По совокупности обнаружений и измерений параметров сигналов производится вторичная обработка, в результате которой строятся траектории движения целей.
Пример расчета зоны обнаружения. Методика расчета зоны обнаружения подробно изложена в [3]. В качестве примера определим срезы
Вторичная обработка
Рис. 5
Далее сигнал поступает в оптимальный накопитель-обнаружитель (рис. 5), состоящий из согласованного фильтра (СФ) для одиночного сигнала и каналов обнаружения для каждого элемента разрешения. Результаты первичных измерений поступают на вход модуля вторичной обработки.
В бистатической системе доплеровский сдвиг частоты обусловливается как движением позиций, так и взаимными скоростями позиций и цели. Поэтому в каждом г-м канале обнаружения
Частотная характееистикаРФ
Рис. 6
областей обнаружения цели БЛА "Boeing X48B" для одной пары "передатчик-приемник", которую можно рассматривать как отдельную бистатиче-скую ячейку многопозиционной системы. При расчетах длина линии базы полагалась равной b = 5 км, угловые положения цели -15 < ац < 15°.
Рабочая длина волны Х = 0.13 м, импульсная мощность излучения передатчика /щ = 50 Вт.
Расчет проведен для двух приемных позиций: приемной позиции на БЛА и наземной приемной позиции. В первом случае вследствие ограниченности геометрических размеров БЛА предусматривалось использование слабонаправленной антенны с коэффициентом направленного действия = 16 дБ, эквивалентными размерами рас-
крыва в азимутальной плоскости Аэ = 0.52 м, в угломестной плоскости Вэ = 0.13 м. При расчете параметров системы относительно наземной приемной позиции на последней предусматривалось использование антенны с коэффициентом направленного действия £Пр = 33 дБ и размерами рас-
крыва Аэ = 2.5 м, Вэ = 1.3 м.
На рис. 7, a представлен срез области обнаружения для случая размещения приемной позиции на БЛА при характерном для этого случая обзоре "сбоку". На рис. 7, б показан аналогичный срез наземной приемной позиции и характерным для этого случая обзором "сверху". Ширина зоны обнаружения цели А одной парой позиций при размещении приемной позиции на БЛА составляет 400 м, при наземной приемной позиции - 1500 м.
Используя несколько пар передающих и приемных позиций, входящих в состав МПРЛС, можно обеспечить формирование общей зоны обнаружения системы.
Потенциальная точность определения координат. Выясним потенциальную точность определения координат целей для МПРЛС с конфигура-
цией, представленной на рис. 3. На практике число передающих позиций может быть и больше. В этом случае необходимо раздельно принимать сигналы, излученные разными передатчиками, на каждой из приемных позиций. Для этого передающие позиции должны излучать ортогональные (или квазиортогональные) сигналы.
Рассмотрим случай, когда первичными параметрами являются суммарные дальности "передатчик-цель-приемник" I (к = 1, 3 - номер
приемной позиций; г - момент выполнения /-го измерения первичных параметров), которые определяются по измерениям временных задержек рассеянных целью сигналов относительно прямого сигнала передатчика. Будем полагать, что цель движется по прямолинейной траектории. Связь между первичными параметрами и параметрами траектории в декартовой системе координат имеет следующий вид:
1к, г (хп ) =
= -хПрк г )2 + ( -Шрк, )2 + ( -^Прк г )2 +
где хп =[хп, уп ,хп ]т - вектор координат цели при п-м последовательном измерении; х/, у/, г. -координаты цели в г-й момент; хпрк ; , УПрк> ^, ¿Прк ■ - координаты приемных позиций в г-й момент; хпд., Упд, ¿Пд - координаты передаю. „ "т" щей позиции в г-й момент, причем - символ
транспонирования.
Положим, что выполнено N последовательных измерений суммарных дальностей "передатчик-цель-приемник" /1 /, ¡2 /, ¡з /, образующих вектор измерений:
у, км 0.12 0
-0.12 -0.24
у, км 0. 5 0
-0.5 -1.0
Рис. 7
а
^ = [11Ь Ьъ l31, , 12N, ^] =
= 1N (хN ) + А1, (2)
где 1 п - нелинейная векторная функция, компоненты которой описываются (1); А1п - вектор ошибок измерений.
Будем считать, что ошибки измерения первичных параметров (суммарных дальностей) независимы для каждого момента времени и распределены по гауссовскому закону с известными
дисперсиями ст2 'к
полагаемыми далее одинако-
выми для всех моментов времени. При этом считаем, что координаты приемников и передатчика в каждый момент времени точно известны.
При использовании метода максимального правдоподобия потенциальная точность определения координат цели определяется информационной матрицей Фишера. Для описанной модели наблюдения элементы матрицы Фишера [2], [6] запишем в виде
' -2 1
■шк (хп) = "М
&
1п Р ( хп )1
50шп0кп 1 1 V Э/1,- (хп )э/1,- (хп)
2 до СТ/ -=1 00шп
эе
кп
1 V Э/2,- (хп ) Э/2,- (хп )
Ст/2 -=1
эешп &екп
_1_ Ъ Э'з,- (хп ) Э'з,- (хп ) + ст2 Ъ
СТ13 -=1
эеш
эекп
(3)
где М {•} - знак статистического усреднения; Р (п | хп ) - условная плотность вероятности вектора измерений 1п; 0ш, Оь, ш, к = 1, 3 - элементы вектора хп : е1п = хп; е2п = Уп; е3п = Ч.
2
СТ У, N
Используя (3), можно найти дисперсии стх N,
2
стг N ошибок определения параметров
траектории как диагональные элементы матрицы, обратной матрице Фишера.
Расчеты выполнены для следующих значений параметров системы: Ь = 5020 м, = 5030 м, ¿3 = = 5025 м, скоростях цели Уц = 80 м/с, передающей позиции Уцд = 10 м/с, приемных позиций Уир1 = Уир2 = 10 м/с, Уир3 = 0. Цель и приемные
позиции Пр1 и Пр2 движутся вдоль оси 0у, начальная высота полета цели 500 м, ее высота полета в конце наблюдения 776 м, полет Пр1 и Пр2 происходит на постоянной высоте. Начальное положение цели задавалось координатами влета в зону обнаружения: Хц0 = 2500 м, уц0 =
= 200 м. Начальное положение приемных позиций пр1: хПр10 = УПр10 = 500 м гПр10 = 500 м;
пр2: хир20 =0, Уир20 =15 00 м, ^пр20 =500 м;
пр3: хПр30 = 0 УПр30 = 1000 м гПр30 = 0, начальное положение передающей позиции: хцд0 = = 5000 м, Упд0 = 950 м, ^цд0 = 1000 м. Средне-
квадратическое отклонение (СКО) ошибок измерения суммарных дальностей положено ст/ = 0.5 м.
На рис. 8 представлены зависимости средне-квадратических ошибок (СКО) определения координат цели: х (а), У (б) и г (в) от координаты У местоположения цели. При расчете кривых 1 в каждой точке траектории движения цели учитывалось только одно измерение набора первичных параметров (N = 1), т. е. в (3) суммы заменялась единственными элементами. На зависимостях 1 наблюдается резкий рост ошибок определения координат при пересечения целью линий баз системы. Это обусловлено тем, что эллипс равной суммарной дальности в момент пересечения линии базы вырожда-
48
32
16 0
200 600
1000
а
±
2
1400 у, м 200 600
=1=
1000 б
Рис. 8
6
3=.
0
1400 У, м 200 600
1000 в
1400 У, м
ах, м
СТУ, м
СТ, М
ется в прямую, что приводит к обращению в ноль производной суммарной дальности по координате x.
Кривыми 2 на рис. 8 показаны зависимости СКО ошибок определения координат от координаты y местоположения цели при учете каждой точке траектории движения вектора из N измерений. При расчете интервал первичных измерений принимался равным T = 0.1 с, интервал времени наблюдения за целью составлял 20 с. Таким образом, в конце наблюдения вектор первичных измерений состоял из 200 элементов. Как следует из полученных зависимостей, использование всех накопленных к текущему моменту измерений первичных параметров при наличии априорной информации о характере траектории движения цели позволяет существенно повысить точность определения координат.
На основании проведенных исследований следует сделать вывод, что цели с малой ЭПР (такие как БЛА) могут эффективно обнаруживаться с помощью многопозиционной просветной РЛС. Для обеспечения требуемой зоны обнаружения передающие позиции должны располагаться на борту воздушных позиций. В качестве таких позиций могут
выступать БЛА или другие летательные аппараты, способные в течение продолжительного времени оставаться на высотах, превышающих высоту движения потенциальных целей. Одним из основных преимуществ этой конфигурации является выгодный угол обзора целей. Обзор "сверху" значительно увеличивает ЭПР малозаметных целей типа БЛА по сравнению с обзором "сбоку", характерном для про-светных РЛС с расположением позиций на земле.
Использование нескольких воздушных приемопередающих позиций позволит перекрыть требуемую зону обнаружения, при этом на борту каждой позиции достаточно разместить передатчик небольших размеров с мощностью в пределах десятков ватт.
Потенциальная точность определения координат в такой системе (на примере координаты х местоположения цели и при расстоянии между позициями порядка 5 км) составляет десятки метров (при использовании одного измерения набора первичных параметров) и единицы метров (при использовании вектора измерения набора первичных параметров).
СПИСОК ЛИТЕРАТУРЫ
1. Moses Al., Rutherford M. J., Valavanis K. P. Radar-based detection and identification for miniature air vehicles, // Proc. of int. conf. on control applications (CCA). Denver, CO, USA, 28-30 Sept. 2011. IEEE: Los Alamitos, CA, USA, 2011. P. 933-940.
2. Черняк В. С. Многопозиционная радиолокации. М.: Радио и связь, 1993. 416 с.
3. Bistatic radar: principles and practice / А. B. Blyakh-man, A. V. Myakinkov, V. I. Kostylev et al. / ed. by M. Cher-niakov. Chichester, England: John Wiley & Sons, Ltd, 2007. 504 p.
4. Blyakhman A. B. Forward scattering bistatic radar // PIERS Workshop on advances in radar methods, Baveno, Italy, 20-22 July1998. Luxemburg: Europ. com. joint research centre. P. 107-113.
5. Blyakhman A. B., Runova I. A. Forward scattering radiolocation bistatic RCS and target detection // Proc. of the 1999 IEEE radar conference. Waltham, USA, April 1999. IEEE: Hoboken, USA, 1999. P. 203-208.
6. Тихонов В. И. Статистическая радиотехника. М.: Радио и связь, 1982. 624 с.
R. S. Fadeyev, A. V. Myakinkov, V. N. Burov, A. G. Ogurtsov Nizhny Novgorod state technical university n.a. R. E. Alekseev
The capabilities of detection and coordinate measurement of low-observable target in multi-static radar with the placement of positions onboard of the UAV
Forward scatter radar cross section of the widely spread models of unmanned aerial vehicles are calculated. The coverage of forward scatter radar when detecting targets of this type is estimated. The configuration of multi-static forward scatter radar enabling the improvement of low-observable target detection efficiency is proposed. The general block diagram of digital signal processing system is designed. The estimations of potential coordinate measurement accuracy in multi-static radar with the placement of positions on the UAV are obtained.
Multi-static forward scattering radar, mobile positions, unmanned aerial vehicle, radar coverage, clutter cancellation, multi-channel coherent integrator, maximum-likelihood method, Fisher Information matrix
Статья поступила в редакцию 22 ноября 2014 г.