CC BY
32
Доклады БГУИР
2016
№ 4 (98)
УДК 621.396.969.1
СИНТЕЗ ОПТИМАЛЬНОГО ДВУХКАНАЛЬНОГО УГЛОМЕРА С ЛИНЕЙНЫМ СКАНИРОВАНИЕМ ДИАГРАММЫ НАПРАВЛЕННОСТИ АНТЕННЫ
Представлены результаты синтеза двухканального измерителя угловых координат цели и описание алгоритма его работы в случае использования линейного сканирования диаграммы направленности антенны. Данный метод позволяет снизить влияние флуктуаций отраженного сигнала на ошибки измерения угловых координат, что по точностным характеристикам приближает его к методам мгновенного сравнения.
Ключевые слова: пеленгация, линейное сканирование, оптимальная обработка, корреляция.
Основным недостатком метода линейного сканирования диаграммы направленности антенны является влияние флуктуации пачки сигнала на точность измерения угловой координаты, что особенно характерно для крупногабаритной цели. Наиболее распространенным и, пожалуй, единственно используемым в настоящее время, способом борьбы с флуктуациями является использование методов мгновенного сравнения, которые требуют для своей реализации два канала с разнесенными по угловым координатам диаграммами направленности для метода амплитудного мгновенного сравнения или с разнесенными фазовыми центрами для метода фазового мгновенного сравнения. В настоящей статье предлагается метод, использующий для снижения флуктуационной ошибки второй канал с несканирующей или слабонаправленной антенной. Канал приема с такой антенной присутствует в ряде радиолокационных станций, а значит данный метод достаточно актуален и легкореализуем без дополнительных модернизаций оборудования.
Известна общая процедура обработки флуктуирующей последовательности импульсов с использованием отношения правдоподобия [1 с. 166-168, 2, 3-5, 6 с. 258]. Выражение для функции правдоподобия имеет вид
ММ. БУЛАШ, СВ. ШАЛЯПИН
ОАО «АЛЕВКУРП»
П. Бровки, 6, Минск, 220013, Беларусь
Поступила в редакцию 18 февраля 2016
Введение
Оптимальная обработка пачки импульсов
Т
где ^(9ц) = ,•••, ,•••, ^м } - отсчеты принятого сигнала, N - количество отсчетов за один цикл обработки; Бе1;(А) - определитель матрицы А; Яф- корреляционная матрица радиолокационного фона (смесь помехи и шума); Кс+ф- корреляционная матрица принятого сигнала (смесь отраженного сигнала, помехи и шума); Офи- к,п-й член обратной
корреляционной матрицы
п - к,п-й член обратной корреляционной
матрицы принятого сигнала Qс+ф = (Я.с+ф) 1; - эрмитово сопряжение; 0ц- истинное
угловое положение цели; 0 ц - оцениваемое угловое положение цели.
С учетом замен Qc0 (0ц) = Q* - (0ц) и ас0 = ln
Det(R$)
ч Det(Rc+* (0 ц)) упростится до следующего вида:
выражение (1)
/
Z (0 ц, 0Ц ) = ^ (0ц )Q(;0 (0 ц )^(0Ц ) + «со, (2)
со
где Q - матрица обработки принятого сигнала; асо- смещение.
Максимально правдоподобная оценка угловой координаты цели соответствует уравнению Z (0 ц, 0ц )| = max. (3)
ц ц
Таким образом, оптимальная обработка пачки импульсов состоит из двух основных этапов:
1) вычисление функции правдоподобия Z (0ц, 0ц) в соответствии с (2);
2) нахождение максимума функции правдоподобия Z(0ц, 0ц д = max и,
10 =0 ц ц
соответственно, определение углового положения цели.
Синтез структуры двухканального измерителя угловых координат
Описанная выше методика является общей для определения угловой координаты цели. Для случая использования двух каналов приема (сканирующего и несканирующего) предлагается использовать отсчеты, принятые обоими каналами, тогда вектор отсчетов принятого сигнала за один обзор имеет следующий вид:
4(6ц) = ЙО' £2Ы-\)Т = ' •' • • 50°Д' • • • • (4)
где - отсчет принятого сигнала, промодулированного ДНА сканирующей антенны
(верхний индекс 0) в «-ом периоде обзора; 'П^'1- отсчет принятого сигнала, промодулированного ДНА несканирующей антенны (верхний индекс 1) в «-ом периоде обзора.
Корреляционные матрицы, используемые для обработки принятого сигнала (4), будут иметь размер 2Ы х 2Ы. Это является главным недостатком и платой за повышение точности измерения координаты цели. Так с ростом размерности матрицы значительно увеличивается количество выполняемых операций при обработке.
Для определения матрицы обработки Qс0 (I0ц) необходимо знание корреляционной
матицы фона Я.ф и корреляционной матрицы смеси отраженного сигнала и фона Я.с+ф.
Корреляционная матрица смеси сигнала и фона Я.с+ф может быть представлена следующим образом:
кс+ф = кс (0ц ) + кф = кс (0ц ) + кш + каш + кп, (5)
где Яс (0ц) - корреляционная матрица отраженного от цели сигнала; Я.ш - корреляционная
матрица внутреннего шума на выходе приемника; Я.аш - корреляционная матрица активной
шумовой помехи на выходе приемника; Я.п - корреляционная матрица пассивной помехи на выходе приемника. Структурная схема измерителя представлена на рисунке.
Считается, что корреляционные матрицы, необходимые для обработки сигнала, известны. При этом вычисляются они по выражениям, представленным ниже.
Структурная схема измерителя Корреляционная матрица сигнала определяется следующим образом [7]:
1*°(01Т) = , } = ст? и? (9п)и? (0П)Ь? , гс I, (6)
4 ц/ ( КА^оЛ) ( А"о="о К,"* ц кл ц К к "о л)
где Я£ ^ п - элемент корреляционной матрицы отраженного от цели сигнала, равный коэффициенту корреляции между щ и щ временными отсчетами сигнала в каналах к0 и к;
ко = 0,1, к\ = 0,1 - номера каналов: 0 - канал сканирующий, 1 - канал несканирующий;
щ = 0, N—1, «1 = 0, N -1 - номера отсчетов сигнала в пределах одного канала обработки (для сканирующего канала или несканирующего канала); акп - СКО п-го отсчета сигнала в к-м
канале обработки без учета влияния ДНА; п (9ц, 9дШп) - множитель для п-го периода повторения, определяемый ДНА к-то канала (сканирующего и несканирующего) в зависимости от текущего положения ДНА 9дна п и углового положения цели 0Ц: п, . = И? , ехр(/Аф^ ) -
нормированный коэффициент корреляции между сигналами в соседних каналах обработки, относящимися к одному и тому же временному интервалу щ = «1; , - модуль коэффициента
корреляции И? , ; Л(р/. /. - фаза коэффициента корреляции п, , ; гс = гс ехр(/А\|/„ „ ) -
0 ' 1 ®' ^ 0 э 1 П0 П0 0' т
нормированный коэффициент корреляции между сигналами в одном канале обработки (кд = к[),
относящимися к разным временным отсчетам; г - модуль коэффициента корреляции г
п0 п
гс ;
И„,И, '
А\|/„ п - фаза коэффициента корреляции гс . При этом номер строки и номер столбца, в которых расположены элементы матрицы, определяются соответственно Пстр = ко • N+по, пстб = к • N+щ •
(7)
Корреляционная матрица внутренних шумов на выходе приемника представляет собой диагональную матрицу размером 2N х 2N и имеет вид
Кш = {Якш,к1,п0,п1} = {2< 5п0,п!}, (8)
где аш (аф) - СКО внутреннего шума в сканирующем (несканирующем) канале, независимость
дисперсии шума от номера периода повторения указывает на стационарность внутреннего шума; 5п п - символ Кронекера, который указывает на некоррелированность междупериодных отсчетов шума;
5к к - символ Кронекера, который указывает на некоррелированность межканальных отсчетов шума. Корреляционная матрица активной шумовой помехи на выходе приемника имеет вид
К3ш (9пап) = {ЯГъ } = ЫГ а?1
V пап/ ^ какщ,п1) ( кащ к,
„аш т таш , г\ \ т таш , г\ \ ;„аш ?
(9)
где Я
к0 ,к1 ,п0 ,П1
элемент корреляционной матрицы активной шумовой помехи (АТТТП); а
аш к
СКО п-го отсчета помехи в к-м канале обработки без учета влияния ДНА, зависимость помехи
от индекса «и» указывает на зависимость дисперсии помехи от времени, это означает, что активная помеха нестационарна. Нестационарность такой помехи может проявляться, например, при постановке уводящей по угловой координате помехи угломерам с линейным сканированием
ДНА; И™(0пап,0дши) - множитель для «-го периода повторения, определяемый ДНА к-го канала (сканирующего или несканирующего) в зависимости от текущего положения ДНА 0дна п и углового положения источника АШП 0 (при измерении угловых координат цели этот параметр является малоинформативным, а при измерении угловых координат постановщиков активных помех (ПАП) - информативным); /т^"". ехр(/Аф^ ¡. ) - нормированный
коэффициент корреляции между сигналами в соседних каналах обработки, относящимися к одному и тому же временному интервалу Пд = Щ; - модуль коэффициента корреляции
к^к ' ^п п ~ символ Кронекера, который указывает на некоррелированность междупериодных
отсчетов помехи, при том, что, в отличие от внутреннего шума, существует сильная межканальная корреляция помехи.
Корреляционная матрица пассивной помехи на выходе приемника имеет вид
кп = } = {^А}, (10)
где члены, входящие в выражение (10), имеют тот же смысл, что и в выражении (6). В отличие
от (6) в (10) предполагается, что и^п (0п, 0дна,п) = 1, такое решение принято в связи с тем, что
точно неизвестно угловое положение сосредоточенной пассивной помехи (МП), кроме того работа СНР может осуществляться при наличии распределенной по угловой координате
пассивной помехи (МО); коэффициент межканальной корреляции равен нулю /ц' ^ = 5/( ^ ;
коэффициент междупериодной корреляции зависит не только от номера периода повторения, но и от номера канала, поскольку междупериодная корреляция пассивной распределенной по угловой координате помехи ниже, чем соответствующий коэффициент корреляции в
несканирующем канале приема, т.е. г,",.,, .^ * г-1^ ^ , а при / Ф/ ^ = 0 .
Измерение углового положения цели
Как было сказано ранее, для определения угловой координаты цели необходимо выполнение равенства (3). При цифровой реализации алгоритма обработки сигнала непрерывное изменение
параметра 0ц заменяют набором его дискретных значений 0^ = 0ц + А0ц/, l = 0, N -1 . Тогда члены выражения (2) преобразуются следующим образом:
Q0 (0ц) ^ Q0 (0ц) = Q/c0;
ас0 (0ц ) ^ ас0 (0ц ) = aC 0; ^0 (0ц ) ^ Rc0 (0ц ) = RCо;
Z (0ц) ^ Z (0 ц) = Z1,
а само выражение (2) принимает вид Z1 = ^Q^ + 0^0.
Для определения максимума Zl необходимо в секторе сканирования определить такое l = lmax, при котором Zlmax = max. Математически это соответствует задаче решения уравнения
dZi
d0 ц
= 0, что оказывается крайне сложной задачей. Поэтому целесообразно
0 =0ц ц lma
воспользоваться квазиоптимальным алгоритмом. Следует заметить, что приведенный ниже алгоритм стремится к оптимальному, когда А0ц ^ 0 или когда А0ц меньше СКО
флуктуационной ошибки измерения угла, обусловленной шумом приемника и фоном.
Полагая, что угловое положение максимума главного лепестка ДН сканирующей антенны 0ц близко к истинному значению координаты цели 0ц, необходимо знать лишь
ближайшие несколько значений вблизи от 0ц. При этом, учитывая, что 0ц является
непрерывной величиной, а 0ц дискретной, то необходим некоторый алгоритм интерполяции для более точного определения истинного значения координаты.
Принимается решение Л = /тах - номер отсчета параметра 0ц, при котором величина
т-Л « о т-Л т-Л-1 т-Л+1
Z является максимальной. Затем по отсчету Z и двум соседним отсчетам £ и £
производится интерполяция (оценка) искомого угла. Обозначим: Хо =0ц + Л0ц (Л-1);
V - 7 Л"1-у0 - 7 ;
VI - 7Л; и У2 - 7Л+1
xi -0Ц + А0ЦЛ; x2 -0Ц +А0ц (Л +1).
0 +А0ц (
По трем дискретам (уо, У1, Хо, Х1, Х2) проводится аппроксимация верхушки функции
правдоподобия полиномом второго порядка (аппроксимация подобна интерполяции при помощи
2 2 2 полинома Ньютона): <Эц -0,5Х0(у1 -у2) + Х1 (у2 -у0) + Х2(у0 -у1) .
Х0 (у1 - у2 ) + Х1( у2 - у0 ) + Х2 (у0 - у1)
Заключение
Подробно описан алгоритм работы двухканального угломера для линейного сканирования диаграммы направленности. В данном алгоритме производится оптимальная обработка пачки импульсов по максимуму функции правдоподобия. При этом для поиска самого максимума предлагается квазиоптимальный алгоритм, который при правильном выборе угла дискретизации А0ц стремится к оптимальному. Первые положительные результаты данного метода по сравнению с классическим одноканальным угломером представлены в [8].
SYNTHESIS OF THE OPTIMUM TWO-CHANNEL DIRECTION FINDER WITH LINEAR SCANNING OF THE ANTENNA PATTERN
MM. BULASH, S.V. SHALIAPIN
Abstract
Target angular coordinates two-channel measurer synthesis and the description of operational algorithm in case of use of the antenna pattern linear scanning are described. This algorithm allows to decrease fluctuations of a reflected signal pack that approximates it to instantaneous comparison methods by its accuracy characteristics.
Keywords: direction finding, linear scanning, optimum processing, correlation.
Список литературы
1. Охрименко А.Е. Основы радиолокации и РЭБ. Ч.1. Основы радиолокации. М., 1983.
2. Kurlovich V.I., Geyster S.R. // Electromagnetic Waves & Electronic Systems. 1999. Vol. 4, № 1. P. 11-16.
3. Kurlovich V.I., Shaliapin S.V., Bashkevich V.Ya. // Electromagnetic Waves & Electronic Systems. 1999. Vol. 4, № 1. P. 17-27.
4. Адхам Салим Аль-Муаз, Шаляпин С.В. // Докл. БГУИР. 2007. № 3 (19) С. 43-49.
5. Ярмолик С.Н., Шаляпин С.В. // Докл. БГУИР. 2003. Т. 1, № 3. С. 24-27.
6. Курлович В.И., Шаляпин С.В. Основы теории радиосистем: учеб. пособие. Минск, 2000.
7. Адхам Салим Аль-Муаз, Шаляпин С.В. // Докл. БГУИР. 2007. №1 (17) С. 61-66.
8. Булаш М.М., Малевич И.Ю., Шаляпин С.В. // Тез. докл. 51-й науч. конф. аспирантов, магистрантов и студентов БГУИР. Минск, 7 мая 2015. С. 13.
