Электронный журнал «Труды МАИ». Выпуск № 71
www.mai.ru/science/trudy/
УДК 621.396.96
Разработка алгоритма получения и обработки спектральных признаков воздушного объекта в интересах идентификации и
селекции движущихся целей
1* 1** 2*** Майоров Д. А. , Перехожев В. А. , Шемяков А.О. .
1 Военная академия войсковой противовоздушной обороны им. A.M. Василевского, ул. Котовского, 2, Смоленск, 214027, Россия Московский Авиационным Институт (национальным исследовательский университет), МАИ, Волоколамское шоссе, 4, Москва, А-80, ГСП-3,
125993, Россия *e-mail:d-mayorov@mail. ru * *e-mail:pvapvo @ mail. ru ***e-mail: [email protected]
Аннотация
Разработан алгоритм получения и обработки спектральных признаков движущегося воздушного объекта в интересах идентификации и селекции движущихся целей при поимпульсной перестройке несущей частоты. Выработаны предложения по реализации способа селекции движущихся целей для перспективных РЛС с перестройкой несущей частоты.
Ключевые слова идентификация, дальностный портрет, селекция движущихся целей, радиальная скорость, сигналы с перестройкой частоты
Введение
Вопросам обработки отраженных от воздушных объектов (ВО) радиолокационных сигналов в интересах их идентификации (радиовидения) в последние десятилетия посвящено множество научных трудов [1, 2]. Интерес к данной теме вызван возможностью определения геометрических особенностей конструкций объектов локации (самолетов, беспилотных летательных аппаратов, ракет) по параметрам отраженных сигналов. Режимы радиовидения всепогодны, работоспособны при больших дальностях ВО, кроме того, они способны существенно расширить информационные возможности радиолокационных станций (РЛС). В результате специальной обработки отраженных сигналов возможно по-
лучение спектральных признаков идентификации ВО, к числу которых относятся радиолокационные изображения (РЛИ), представляющие собой распределения отражательных свойств объектов в радиолокационном диапазоне волн по одной или нескольким координатам [2]. В зависимости от числа используемых координат РЛИ принято делить на одномерные (портреты) и многомерные [3]. Практическую значимость на современном уровне развития техники играют двумерные РЛИ и портреты (дальностные и доплеровские).
Кроме обеспечения идентификации ВО дальностные портреты (ДП) могут быть использованы в перспективных РЛС с поимпульсной перестройкой частоты в интересах обнаружения движущихся на фоне пассивных помех воздушных объектов [4, 5]. Способ, предложенный в [4, 5], предполагает специальную обработку комплексных частотных характеристик (ЧХ) и дальностных портретов ВО для обеспечения селекции движущихся целей (СДЦ). В [4, 5] теоретически обоснована и доказана возможность обеспечения СДЦ в РЛС с поимпульсной перестройкой частоты, но отсутствуют предложения по технической реализации предложенного способа.
Целью работы является выработка предложений по реализации представленного в [4, 5] способа в перспективных РЛС с поимпульсной перестройкой частоты и разработка алгоритма формирования и обработки дальностных портретов ВО в интересах обеспечения идентификации и СДЦ.
Предложенный в [4, 5] способ селекции движущихся целей предназначен для использования в перспективных станциях обнаружения целей с поимпульсной перестройкой несущей частоты. Для его реализации в радиолокационных станциях кругового обзора пространства предлагается разбить зону обнаружения на В=2п/®р азимутальных каналов, где ©в - ширина диаграммы направленности антенны по азимуту. Способ предполагает отсутствие амплитудной модуляции сигналов из состава одной пары пачек СПЧ характеристикой направленности антенны при ее вращении. На современном этапе развития радиолокационной техники указанное требование может быть обеспечено с помощью неподвижных цифровых фазированных антенных решеток прямоугольной (секторного обзора) или цилиндрической формы (кругового обзора).
Для реализации предложенного способа СДЦ в РЛС, использующих вращающуюся антенну, скорость ее вращения ша и длительность пачек сигналов с перестройкой частоты (СПЧ) Тр необходимо выбирать из условий, позволяющих пренебречь указанной амплитудной модуляцией. Из [6, 7] известно, что характеристику направленности антенны в азимутальной плоскости можно описать функцией арочного синуса f (©)=8тс[п Ь 8т(©)/Л,], где © - угол между направлением на расположенную в дальней зоне точку наблюдения и нормалью антенны, Ь - линейный размер антенны. Исследование зависимости модуляции
отраженных сигналов от скорости вращения антенны целесообразно проводить в пределах ширины луча по половинной мощности (ширины диаграммы направленности антенны). Очевидно, что при таком подходе максимальные изменения функции f (О) будут наблюдаться при значениях угла О, близких к ±О^/2. На рис. 1 изображен основной лепесток диаграммы направленности антенны и его фрагмент при значениях О, близких к 0^/2.
-ЛИ -0^/2 о 0/?/2 ЛИ © 0^/2-0^/10 0^/2-0^/20 0^/2-0^/40 0^/2 © Рис. 1. Диаграмма направленности антенны
Из рис. 1 видно, что при изменении переменной О на величину 0^/40 (повороте антенны на угол 0^/40) значение функции f (О) изменяется не более чем на 0,026, при повороте на угол Ов/20 - не более чем на 0,05, при повороте на угол О^/Ю - не более чем на 0,1. Для возможности пренебрежения изменениями функции ее модуляция не должна превышать единицы процентов (до 5). На этой основе выработано требование к излучению одной пары пачек СПЧ за время поворота антенны на угол Ов/20:
<
(О в /20) = 2Т„ 40 Тр
(1)
Следует отметить, что выбор в два раза меньшей скорости вращения антенны
О в
(Ш а =-) уменьшает амплитудную модуляцию сигналов при обзоре пространства,
80 Тр
но при этом увеличивает цикл обзора в 2 раза. В этом случае скорость вращения антенны при круговом обзоре пространства необходимо выбирать с учетом задачи, стоящей перед РЛС обнаружения, но не выше указанной в (1).
Согласованная фильтрация СПЧ, излученных в квазинепрерывном режиме, возможна только при их многоканальной обработке в приемнике РЛС обнаружения. При этом для правильной расфильтровки сигналов длительность импульсов должна выбираться из условия т >1/2А/, где А/ - шаг перестройки частоты. При таком выборе параметров зондирующего сигнала обеспечивается прохождение основного лепестка спектра отраженного импульса п-й частоты в п-й частотный канал. В интересах обеспечения однозначного измерения дальности период повторения Тр импульсов одной частоты (пачек) должен определяться наклонной дальностью до дальней границы зоны обнаружения Яр:
Тр > 2Яп /с. (2)
Обязательным условием работоспособности предложенного в [4, 5] способа СДЦ является излучение пар пачек СПЧ с одинаковым законом перестройки частоты внутри пары. В этом случае легко обеспечить выполнение условия (2), так как в пределах пары период повторения импульсов одной частоты постоянен. Однако при излучении следующей пары пачек СПЧ с другим законом перестройки частоты выполнение условия (2) становится проблематичным, так как период повторения импульсов п-й частоты может уменьшиться. Для обеспечения однозначного измерения дальности необходимо каждую следующую пару пачек СПЧ излучать через время задержки 13=Тр после излучения предыдущей пары (рис. 2). Следовательно, период повторения пар пачек СПЧ составит 3Тр.
О 0^/20 30^/20 60^/20 90^/20 120^/20 150^/20 180^/20 ©^ р
Ов
Рис. 2. Вариант последовательности зондирующих сигналов при ША =-
40 Тр
При формировании такой последовательности зондирующих сигналов с учетом (1) за время поворота антенны на угол О в в направлении ВО будет излучено как минимум 13 пар пачек СПЧ (рис. 2). Указанное число получено для максимально возможной скорости антенны ша = Ов /(40Тр) в результате решения следующей пропорции:
1 пара пачек - 3Тр; х пар пачек - 40 Тр.
При уменьшении скорости вращения антенны количество излученных в направлении ВО пар пачек СПЧ будет пропорционально увеличиваться. Например, при скорости антенны ша = Ов /(60Тр) число излученных в направлении ВО пар пачек СПЧ будет равняться 20.
Согласно предложенному в [4, 5] способу СДЦ решение о наличии ВО принимается после обработки трех пар пачек. Этим объясняется объединение пар пачек в тройки. Принадлежность ВО к одному из азимутальных каналов определяется положением нормали антенны в момент излучения второй пары пачек СПЧ из состава тройки. Обработка 13 пар пачек СПЧ позволит получить не менее 4 решений (назовем их промежуточными)
о наличии ВО в Ъ-м азимутальном канале, где Ъ = 0, В — 1 - номер азимутального канала. Окончательное решение о нахождении ВО в Ъ-м азимутальном канале принимается при выработке хотя бы одного промежуточного решения из четырех возможных. Как уже отмечалось, в интересах исключения явлений «слепых» дальностей и скоростей, присущих режиму квазинепрерывного излучения, период повторения импульсов каждой второй тройки пар пачек рекомендовано выбирать отличающимся от периода повторения остальных троек на т/. При излучении второй тройки пар пачек СПЧ с периодом повторения импульсов Т 2=(21+1) Т', где Q1 - скважность сигналов из состава первой тройки, количество зон «слепых» дальностей и скоростей уменьшается в Q2=(Ql+1) раз, что повышает вероятность обнаружения ВО.
Следует отметить особенности принятия решения о наличии ВО в каком-либо из азимутальных каналов. На рис. 3 показаны различные варианты возможного положения ВО относительно границ азимутальных каналов.
Азимутальный Азимутальный
РЛС
Рис. 3. Варианты положения ВО относительно границ азимутального канала
Правильное решение о принадлежности ВО Ь-му азимутальному каналу будет принято в случае нахождения ВО в центре азимутального канала Ь (рис. 3, ВО №1). В другом предельном случае, когда ВО находится на стыке двух азимутальных каналов (рис. 3, ВО №2), также будет принято правильное решение о нахождении ВО в Ь-м и (Ь+1)-м азимутальных каналах. В остальных случаях будет принято решение о наличии ВО в двух азимутальных каналах, хотя реально ВО находится только в одном из них (рис. 3, ВО №3). Это объясняется тем, что при нахождении нормали антенны на границе одного из азимутальных каналов на вход приемного устройства поступают сигналы, отраженные от ВО, находящейся в соседнем азимутальном канале. Такая ситуация является обычной для РЛС обнаружения. Точность определения азимута ВО составляет 20^.
Таким образом, для реализации предложенного в [4, 5] способа СДЦ необходимо излучать пары пачек сигналов со случайным законом перестройки частоты и периодом следования пар пачек 3Тр. Число импульсов в пачке выбирается Ы=2к, где к - целое число, принимающее значение от 5 до 9, в прямой зависимости от сложности помеховой обстановки. Пары пачек объединяются в тройки (рис. 2), период повторения импульсов каждой четной тройки пар пачек Тг'2 отличается от периода повторения нечетных троек Тц на ти При этом длительность каждой пачки не должна превышать половины интервала угловой корреляции Тук [8]. Длительность пачек из состава четных (Тр2) и нечетных (Тр1) троек, скорость вращения антенны, длительность и период повторения импульсов внутри пачек выбираются из условий:
< т < Т1; Т/1 =Тр1 N =QlТi; с р1 2
т/ >1/2 А/;
Тр2 =Ы Т/2=Ы (Т/1+т^^1+1)т/ Q2 Т/;
Ту1 ^ О в
Тр2 ; Тр2 <-—-; <—-в-
Р 2 2 р2 4У_ (/ ) а 40 Тр
Для излучения сигналов со случайным законом изменения частоты в цифровой вычислительной системе необходимо сформировать последовательность номеров частот излучения, используемых в пачке СПЧ, и распределить их по случайному закону. При этом будет сформирован вектор Г, п-й элемент которого определяет порядковый номер п излучения импульса на п-й частоте. В интересах последующей обработки данные вектора Г сохраняются.
После приема отраженных от ВО сигналов на разных частотах /0+пА/+Гй п, где Гй п - доплеровская частотная добавка отраженного сигнала на п-й частоте, их частоту необходимо понижать до промежуточной /пр+пА/+Г п , где /пр - величина промежуточной частоты. После понижения частоты становится возможным усиление сигналов с помощью широкополосных усилителей промежуточной частоты и их разделение на частотные каналы с помощью набора полосовых фильтров.
Далее предлагается выделять с помощью квадратурных фазовых детекторов квадратурные составляющие ОС в каждом частотном канале, преобразовывать их в цифровую форму с помощью аналого-цифровых преобразователей. Затем необходимо проводить цифровую согласованную фильтрацию принятых сигналов и преобразование откликов согласованных фильтров в комплексную форму.
Затем в каждом дальностном канале формируется вектор 01 (ЧХ) из N элементов, при этом в п-й элемент этого вектора записывается комплексное значение отклика согласованного приемника на п-й частоте в точке опорной дальности й-го канала дальности (или п-й член ЧХ). Данная операция поясняется рисунком 4. Она аналогична операции создания двумерной матрицы из D*N элементов.
Дальностные каналы
«Слепая» дальность
Частотные каналы
Частотные характеристики (векторы С\) Рис. 4. Пояснение к процессу формирования ЧХ в разных каналах дальности
В интересах облегчения толкования алгоритма СДЦ ограничимся рассмотрением процесса обработки информации в ^м канале дальности. Частотная характеристика, полученная в результате обработки первой пачки (вектор 01Д сохраняется в оперативном запоминающем устройстве для ее последующего вычитания из перефазированной частотной характеристики, полученной в результате обработки второй пачки. Для осуществления перефазировки каждый элемент вектора 02^ умножается на комплексный фазовый множитель
Xn, vrПП = exp(-j4пfnVrППТр / c) ,
где vrПП = - УгПП max, - УгППmax+^УгПП ? — ? УгППmax-
В результате будет сформирована двумерная матрица G2d данных из У1=(2УгПП max/АУгПП+1) строк и N столбцов. Далее необходимо производить вычитание n-х элементов вектора G1d из строк матрицы G2d, в результате которого будет сформирована двумерная матрица Sd данных из V1 строк и N столбцов. Перечисленные операции схематично поясняются рисунком 5.
2-я пачка СПЧ Вектор С/2,/
и 0 й1 и2 и N-1
Матрица С2и
*0,0
*0,1
-4,0
П,1
2,0
*2,1
и (2) • Х
ип Ли, УгГ
-1, 0
^-1, 1
-1,У 1-1
Матрица
1-я пачка СПЧ
и0.у 1-1 и1,у 1-1 и2.У1-1
«0,0 «1,0 М2,0 и N-1, 0
и0,1 и1,1 и2,1 и N-1, 1
и0,У1-1 и1,у 1-1 и2,У1-1 и N-1,У 1-1
Рис. 5. Формирование разностных ЧХ в (-м канале дальности
Следующим этапом является перефазировка разностной ЧХ путем умножения каждого п-го элемента на фазовый множитель (&п г = ехр[- ]4п/пУг0 + пТ1)/с]. Значение п берется из сохраненного вектора К В связи с тем, что значение радиальной скорости неизвестно, производится перебор возможных значений скорости уг в диапазоне ±УГтах с шагом ДУ. В результате будет сформирована трехмерная матрица Жданных из Ы*У2 строк и У1*У2 столбцов, где У2=2УГ тах/ДУ+1.
Далее предлагается путем проведения обратного быстрого преобразования Фурье с комплексными векторами данных каждой строки трехмерной матрицы Wd формировать трехмерную матрицу W1( , равную по размерам матрице Wd. После этого необходимо находить максимальное значение модуля комплексного сигнала в матрице W1( и делить комплексные величины всех элементов этой матрицы на найденное значение, т. е. нормировать элементы матрицы W1(, а затем - рассчитывать величину энтропии данных Н для каждой строки матрицы W1(. В результате будет сформирована двумерная матрица Н( из
У1*У2 элементов. Указанные операции поясняются рисунком 6. Первые индексы элемен-
тов матриц, показанных на рис. 6, обозначают номер частоты, вторые - подбираемое значение радиальной скорости ПП, третьи - подбираемое значение радиальной скорости ВО.
Матрица
Рис. 6. Формирование скоростной развертки в d-м канале дальности
На конечном этапе находится номер столбцов v r min и v r пп mn, соответствующих наименьшему значению энтропии Н( v r min, v r пп min), с помощью которых определяется оценка радиальной скорости ВО и ПП.
/V
Vr = -Vr max+ ( V r min - 1)А V ;
Vr ff = -Vr пп max+ ( V r ПП min - 1)А Vr ПП-
Для исключения возможности принятия решения о наличии ВО по «ложному» минимуму [9] производится сравнение значений Н (V г тт, V г пп тп) и Н( V г т[п , V г пп тт ±25/АУг пп), при этом знак вычитания выбирается при
/V /V
Уг > 0, а знак сложения - при Уг < 0. Если разница между указанными значениями не превысила 1, принимается решение об отсутствии ВО, в матрице Нd создается окно запрета Н (V г т'п ±10/А Уг , V г пп т'п) и поиск минимального значения в матрице На повторяется еще
один раз. Полученная оценка радиальной скорости V^р), где р - порядковый номер обрабатываемой пары пачек из состава тройки, сохраняется в специальной двумерной матрице размером 3*Д обновляемой после обработки каждой тройки пар пачек СПЧ.
Для исключения ложной тревоги и принятия правильного решения о наличии ВО в обрабатываемом дальностном канале необходимо излучение и обработка как минимум еще двух пар пачек СПЧ с другими законами изменения частоты. Обработка каждой пары пачек СПЧ предложенным способом позволяет получить еще две оценки радиальной ско-
рости ВО У/и У/ р+2)
. При выполнении условий
У(2) - У/
г г
(1)
<ЛУ и
У/(3) - У/
г г
(2)
<ЛУ и
У/(3) - У/
г г
(1)
<ЛУ
где ДУ - порог, определяемый точностью оценки (измерения) скорости, принимается решение о наличии в ( -м дальностном канале ВО, движущегося с радиальной скоростью
А I .Лч / \ Ль / л \ А , 1 /
V, =Ур) + уг( р+1) + у; р+2) )/ з . Решение о принадлежности обнаруженного ВО к какому-либо из азимутальных каналов принимается с помощью специального кода, несущего в себе информацию о положении нормали антенны в момент излучения второй пары пачек СПЧ из состава тройки.
Следует отметить, что получаемые с помощью такой обработки спектральные признаки ВО (в виде дальностных портретов) могут быть использованы в интересах их идентификации.
Алгоритм, реализующий предложенный в [4, 5] способ в РЛС с поимпульсной перестройкой частоты, представлен на рис. 7.
Начало
Выбор длительности пачек, скорости вращения антенны, периода повторения и количества импульсов в пачке
г 2 ■
Формирование вектора ^ , определяющего закон перестройки частоты для пары пачек
Измерение углового положения нормали антенны в целях определения номера азимутального канала Г 4--
Излучение 2-х пачек сигналов с одинаковым законом перестройки частоты
5-1-
Прием отраженных сигналов,
их преобразование на промежуточную частоту, усиление, разделение на частотные каналы
7 ■
Формирование В векторов G2 из оцифрованных откликов согласованных фильтров
■10
Формирование В двумерных матриц G2 путем домножения на фазовые компоненты, соответствующие подбираемой ¥г пп -1 -
г 12
Вычитание п-х элементов векторов G1 из п-х элементов строк матриц G2 (формирование В двумерных матриц £)
г 14
Нормировка элементов матриц W1, расчет энтропии данных в строках матриц W1 (формирование В матриц Н)
Формирование В векторов G1 из оцифрованных откликов согласованных фильтров
г 6
Выделение квадратурных составляющих отраженных сигналов, их преобразование в цифровую форму
г 9
Согласованная фильтрация принятых сигналов, преобразование откликов согласованных фильтров в комплексную форму
г 13-
Формирование В трехмерных
матриц Ж из двумерных путем дом-ножения на фазовые компоненты, соответствующие подбираемой Уг
г 15
Проведение обратного БПФ с комплексными векторами данных каждой строки матрицы W (формирование В матриц W1)
Поиск минимального значения в матрицах Hd
Vr =Hd (vr min, vr ПП min)
г 18
Создание окна запрета в матрицах Н/, имеющих «ложный» минимум, повторный поиск экстремума в матрицах
Заполнение двумерной матрицы, состоящей из 3 *В элементов
Дальностный канал /, азимутальный канал Ь: «ЦЕЛИ НЕТ»
Формирование вектора V, состоящего из В элементов У/( р ) = У
г г
Дальностный канал /, азимутальный канал Ь: «ЦЕЛЬ ЕСТЬ»; вывод Vr
Уг = (\/г( р ) + У( р+1 + У( р+2) )/зу
Рис. 7. Алгоритм получения и обработки ДП ВО в интересах идентификации и СДЦ
Таким образом, в ходе исследований разработан алгоритм получения и обработки спектральных признаков движущегося воздушного объекта в интересах идентификации и селекции движущихся целей в перспективных РЛС с перестройкой несущей частоты. Для его реализации в РЛС кругового обзора зона обнаружения разбивается на азимутальные и дальностные каналы. Зондирующие сигналы представляют собой последовательность пар пачек СПЧ с периодом повторения между парами, равным трем длительностям пачки. Закон изменения частоты для каждой пары пачек различен, но внутри пары является одинаковым, период повторения импульсов изменяется от тройки пар пачек к тройке. Разделение сигналов на частотные каналы предлагается проводить с помощью набора полосовых фильтров. Решение о наличии ВО принимается на основе результатов обработки каждых трех пар пачек СПЧ в случае, если разница между тремя полученными оценками скорости не превышает установленный порог. Решение о принадлежности ВО Ъ-му азимутальному каналу вырабатывается на основании информации о положении нормали антенны в момент излучения второй пары обрабатываемых пачек СПЧ из состава тройки.
Исследования проводились в рамках проекта №13-07-97505 при финансовой поддержке администрации Смоленской области и Российского фонда фундаментальных исследований.
Библиографический список
1.Митрофанов Д. Г. Нейросетевая идентификация воздушных объектов по радиолокационным изображениям. Москва. Измерительная техника. № 2. 2006. с. 53-59.
2.Митрофанов Д. Г. Способ построения двумерного радиолокационного изображения воздушной цели с произвольной нестабильной траекторией полета. Москва. Измерительная техника. № 8. 2002. с. 7-11.
3.Митрофанов Д. Г. Комплексный адаптивный метод построения радиолокационных изображений в системах управления двойного назначения. Москва. Теория и системы управления. Известия РАН. № 1. 2006. с. 101-118.
4.Патент РФ 2419107, МКИ' 0018. Способ
селекции движущихся целей в режиме поим-пульсной перестройки несущей частоты/Митрофанов Д. Г., Майоров Д. А. и др.; опубл. 20.05.11.
5.Майоров Д. А. Применение сигналов с перестройкой частоты в интересах селекции движущихся целей. М., МАИ, Электронный журнал «Труды МАИ». Выпуск 52.
6.Фрадин А. З. Антенно-фидерные устройства. - М.: Связь, 1977.
7. Антенны УКВ / Под ред. Айзенберга Г. З. - М.: Связь, 1977.
8.Бартон Д. К., Вард Г. Р. Справочник по радиолокационным измерениям. - М.: Сов. радио, 1976.
9.Лифанов Б. А. Математическое моделирование процесса селекции движущихся целей при перестройке несущей частоты. Труды Российского научно-технического общества радиотехники, электроники и связи им. А. С. Попова. Серия: цифровая обработка сигналов и ее применение.2012. Вып. Х1У-1. С. 260-263.