CC BY
111УДК 621.396.626
ЦИФРОВОЕ ФОРМИРОВАНИЕ НУЛЕЙ ДИАГРАММЫ НАПРАВЛЕННОСТИ ФАР ПРИ ПЕРЕМЕЩЕНИИ ПОСТАНОВЩИКА АКТИВНЫХ ПОМЕХ
Штрунова Екатерина Сергеевна
кандидат технических наук, доцент ФГБОУ ВПО «Рязанский государственный радиотехнический университет». E-mail: shtrunova.e.s@rsreu.ru. Адрес: 390005, Рязань, ул. Гагарина, 59/1.
Аннотация: В условиях радиоэлектронного противодействия угловое перемещение постановщика активных помех за интервал времени между этапами адаптации приводит к изменению направления прихода активной шумовой помехи. В связи с этим помеха выходит из зоны ранее сформированных нулей диаграммы направленности приемной фазированной антенной решетки радиолокационной системы. В статье оценено снижение коэффициента улучшения отношения сигнал - (помеха+шум) оптимального алгоритма непосредственного обращения выборочной межканальной корреляционной матрицы помехи при приращении направляющих косинусов активной помехи за время между этапами адаптации по сравнению с неподвижной помехой. При угловой скорости перемещения постановщика активных помех на эффективность алгоритмов адаптивной обработки наибольшее влияние оказывает активная помеха по ближнему лепестку диаграммы направленности фазированной антенной решетки при приближении ее углового направления к области главного лепестка. Для уменьшения влияния угловой скорости постановщика целесообразно минимизировать время между интервалами адаптации таким образом, чтобы можно было пренебречь изменением углового положения помехи либо же подстраивать положение нуля диаграммы направленности в зависимости от угловой скорости постановщика с учетом максимизации отношения сигнал - (помеха+шум).
Ключевые слова: активные шумовые помехи, адаптация фазированной антенной решетки, угловая скорость постановщика помех.
Введение
Задача адаптивной антенной решетки - обеспечение эффективного функционирования систем радиолокации, радионавигации в условиях радиоэлектронного противодействия, в частности, на фоне многоточечных по пространству активных шумовых помех (АТТТП) Применение фазированных антенных решеток (ФАР) и устройств цифрового запоминания пространственных частот фЯБМ^) в станциях постановки помех позволяет с высокой скоростью изменять направление прихода, количество АТТТП, их мощность; поддерживая определенный закон переключения АТТТП, формировать мерцающие активные помехи [1, 2].
Поэтому к адаптивным алгоритмам формирования диаграммы направленности (ДН) антенны РЛС предъявляют жесткие требования по времени адаптации, частоте настройки вектора весовых коэффициентов (ВВК), характеризующих амплитудно-фазовое распределение на раскрыве приемной ФАР. Дополнительным
фактором, влияющим на эффективность алгоритмов помехозащиты на основе пространственной обработки сигналов, является движение постановщика активных помех (ПАП), так как происходит смещение направления прихода АТТТП за интервал времени между этапами адаптации относительно предварительно сформированного нуля ДН ФАР.
Одним из методов борьбы с данным недостатком является применение алгоритмов формирования ДН с расширенными нулями в направлении действия источников АТТТП Основой этих алгоритмов является введение дополнительных условий получения ВВК.
Применительно к линейной адаптивной ФАР в литературе рассмотрены следующие подходы по формированию расширенных нулей ДН:
1) в [3] вносится предположение о наличии дополнительных дискретных равномощных ПАП, расположенных на одной линии, а в [4] о непрерывном распределении источников;
2) весовое оконное сглаживание направляющих векторов [5];
3) условие равенства нулю производной ДН [6];
4) усреднение пространственной межканальной корреляционной матрицы (МКМ) по угловым координатам сигналов за интервал времени [7].
К недостаткам алгоритмов формирования ДН с расширенными нулями относится снижение глубины формируемого нуля при увеличении ширины зоны подавления АТТТП, а при действии АТТТП по ближним боковым лепесткам может произойти подавление полезного сигнала. В [8] перемещение ПАП учитывалось путем изменения МКМ помехи во времени, но точного закона изменения модели МКМ при движении не предложено.
Цель работы - исследование влияния угловой скорости постановщика помех на эффективность адаптации цифровой антенной решетки.
Пространственно-временное представление сигнала и АШП
В общем случае АТТТП представляет собой широкополосный в пространственно-временном смысле процесс. При условии нормальности и стационарности принимаемого воздействия распределение сигнала и АТТТП полностью определяется функцией пространственно-временной корреляции (ФПВК) R(tb t2, d1, d2), которая связана обратным четырехмерным преобразованием Фурье с частотно-волновым спектром, характеризующим распределение мощности по временной f и пространственным q частотам [1]:
R(i,p) = J ШGf,q) exp(i2rc(f i -qp))dfdqxdqydqz
-'Xi
= J JJJG(f,q) exp(i2xfi-kp))dfdkxdkydkz,
-да
где q - вектор линейных пространственных частот с проекциями qx = ux / X, qy = uy / X, qz = uz / X , ux, uy - направляющие векторы (проекции единичного вектора на оси x и y), ux= cos ux, uy= cos uy, ux,y= x,y/X - угловая пространственная частота линейного раскрыва по
координате х и у соответственно; к - волновой вектор, к = 2п q, к = (кх, ку, к); р - разность координат двух точек наблюдения, р = (рх, ру, р2), р = d2 - d1, dг■ -радиус-вектор, определяющий положение /-го антенного элемента на апертуре ФАР.
Для обработки узкополосных пространственно-временных сигналов входной реализации возможна факторизация [9, 10], поэтому функцию пространственно-временной корреляции представляют отдельно как временную и пространственную: R(т, р) = R(т) R(р). Функция пространственной корреляции с учетом обобщения теоремы Винера - Хинчина на пространственный случайный процесс [1]
R(р) = {{{ G(q) ехр^2^р))^Зр,
—ад
где G(q) — угловой спектр мощности принимаемого процесса, характеризующий распределение мощности по пространственным частотам.
В адаптивных алгоритмах пространственной обработки сигналов, основанных на прямом обращении выборочной МКМ и рекуррентном оценивании обратной МКМ [11], ре-гуляризованных, проекционных [12, 13] алгоритмах определение оценки ВВК XV основано на предварительном нахождении выборочной МКМ или обратной МКМ.
Для адаптации Л-элементной ФАР Л-мерная оцениваемая МКМ аддитивной смеси коррелированной по пространству АШП и собственного шума приемных каналов представляется через вектор характеризующий пространственное положение помехи, как [11]
Кщл = I + у К Ь (1)
где у - отношение АТТТП - шум по мощности; I- Л-мерная корреляционная матрица собственного шума приемных каналов ФАР (единичная); ^ = (ехр(^под),..., ехр^пкД , ехр^^щ^]}* - направляющий вектор помехи для Л-элементной плоской ФАР прямоугольного раскрыва при условии совпадения начала координат с расположением нулевого элемента; vпkj - фаза помехи, пришедшей на к]-й элемент ФАР, Vпкj = 2л/поА^пк,;-; /по - частота излучения АШП; А^ - время запаздывания фрон-
та волны АТТТП, принимаемой к/-м элементом относительно начала координат для угловых координат ап, 0п прихода АТТТП, Л^пАу = [к dxux+j dy иу]/с, dx, dy - межэлементное расстояние по координате х и у соответственно, их, иу - направляющие косинусы (проекции единичного вектора на оси х и у), их = sm9пcosап, иу = sin9п sinап; ап, 9п
- азимут и угол места АТТТП соответственно, i
- мнимая единица.
На рис. 1 представлено изменение координат ПАП при его перемещении.
HZ
Рис. 1.
С учетом влияния угловой скорости ПАП оценка МКМ помехи примет вид:
Й пш= Киш+ЛЯ, (2)
АЯ - добавка к МКМ из-за движения ПАП на расстоянии D от РЛС.
Если ПАП имеет скорость перемещения такую, что приращение направляющих косинусов и расстояния от фазового центра ФАР до ПАП на интервале наблюдения постоянны, то координаты ПАП будут в (/+1)-й момент времени равны:
(+1):
(0^du*
: ±
dt
t, и
(i+1) :
= uy® ±.
du
dt
-t,
Eft+V = ± dDt.
dt
dD
где--приращение расстояния от фазового
dt
du du у
центра ФАР до ПАП; х у
dt dt
приращения
направляющих косинусов по координате х, у соответственно за время между этапами адаптации, величина которых зависит от тангенциальных составляющих вектора скорости по азимуту и углу места соответственно [9]:
Va = D
ve = D
du* dt
du*
-cosa --
du , dt
du ,
-sina
/cose,
-sina —
cosa
/ sine.
dt dt
\ у
Тогда временена задержки помехи в (/+1)-й
момент времени
A^(í+1) = [kdx(u(í) ±du^t) + jdy(ut)]/c.
doy dt
Цифровое формирование ДН ФАР
В отечественной и зарубежной литературе особое внимание уделяется адаптивному формированию нулей ДН антенны в направлении на ПАП, реализуемому в цифровом виде [14 - 17]. Принцип цифрового формирования ДН ФАР основан на дискретном преобразовании Фурье (ДПФ) амплитудно-фазового распределения ФАР. С учетом того, что пространственная избирательность ФАР определяется множителем решетки [18], то для плоской ФАР с прямоугольным раскрывом его можно представить в виде двумерного ДПФ амплитудно-фазового распределения:
N-1M-1
2imd, 2nmdy Л
-x ux +--u,
X x X y ,
МЧх, и у) КУ" - е
п=0 т=0
где Жпт, ' пт - амплитуда и фаза соответственно оценки весового коэффициента для элемента ФАР п-й строки и т-го столбца, N М - количество строк и столбцов в ФАР соответственно.
Пространственную выборку входной реализации можно записать:
В(их, иу, /) = {¿(/ + тг )е12л /д t + р(1 )}х
N-1М-1 л , 12ят[| 12л
хЦ Кт^'пт е 1 ' ^ .
п=0 т=0
где + хг) - комплексная огибающая отраженного от цели сигнала, содержащая
и
x
информацию о дальности объекта; тг - время задержки отраженного сигнала относительно переданного, определяемое наклонной дальностью до цели; / - частота Доплера, р(0 - АШП.
Для оптимального алгоритма непосредственного обращения матрицы оценка ВВК определяется как собственный вектор, соответствующий максимальному
собственному значению матричного
произведения XV <-К—ш Яс, где обратная
МКМ (1) в соответствии с тождеством
уТ Г
/и п
Вудбери: Rnm =1 -
1 + Yfn fn
Rc - МКМ
сигнала.
В случае движущейся АШП оценка МКМ при адаптации не совпадает с МКМ выражения (2). Тогда мощность помехи на выходе системы обработки при адаптации ФАР без учета угловой скорости ПАП:
( ; ( ; ;н ^
Pn+m = Rc
yf fc
< П П
1 + yfn fn
R„
I
yf fc
¡ П П
1 +yfn fn
Rc
п п / V 1 п п /
Для оценки влияния угловой скорости на качество подавления активной помехи применялось снижение Дц коэффициента улучшения отношения (сигнал - помеха+шум) (ОСПШ) оптимального алгоритма
непосредственного обращения МКМ в случае неподвижной помехой по сравнению с угловым смещением АШП.
На рис. 2 и 3 представлены зависимости снижения Дц коэффициента улучшения ОСПШ алгоритма формирования ДН от приращения направляющего косинуса за время между этапами для линейной эквидистантной ФАР с N = 15 при действии помехи по ближнему боковому лепестку ихп=0,1 и по лепесткам с ихп=0,2; 0,3; 0,4 соответственно.
Из рис. 2, 3 можно сделать вывод, что при использовании оптимального алгоритма непосредственного обращения МКМ для линейной ФАР N=15 наибольшее внимание стоит уделять случаям действия АШП по ближнему лепестку ДН ФАР при смещении помехи в область главного лепестка, так как уменьшение коэффициента улучшения происходит на 1-5 дБ.
Дц>
ДБ 1.5
0.5
0
г__=0.4 - /
---------------- г^хп = 0-2
_ __ __ = о-з^
-0.01 -5x10"3 0 Рис. 3.
5x10"3 dux/dt
Заключение
Для увеличения эффективности адаптивной обработки целесообразно выбирать время между интервалами адаптации таким образом, чтобы можно было пренебречь смещением помехи, либо же подстраивать положение нуля ДН в зависимости от угловой скорости ПАП с учетом максимизации отношения сигнал -(АШП+шум).
Литература
1. Добыкин В.Д., Куприянов А.И., Пономарев В.Г., Шустов Л.Н. Радиоэлектронная борьба. Цифровое запоминание и воспроизведение радиосигналов и электромагнитных волн / под ред. А. И. Куприянова. - М: Вузовская книга, 2009. -360 с:
2. Радиоэлектронное подавление информационных каналов систем управления
I
оружием / под ред. Ю.М. Перунова. -М.: Радиотехника, 2008. - 416 с.
3. Mailloux R.J. Covariance matrix augmentation to produce adaptive array pattern troughs // Electron. Lett. - 1995. - Vol. 31. - P. 771-772.
4. Zatman M. Production of adaptive array troughs by dispersion synthesis // Electron. Lett. - 1995. - Vol. 31, No. 25. - P. 2141-2142.
5. Нечаев Ю.Б., Борисов Д.Н., Пешков И.В. Алгоритм формирования диаграммы направленности кольцевой антенной решетки, устойчивый к многолучевому распространению и нестационарным источникам помех // Известия вузов. Радиоэлектроника. - 2011. - Т. 54, № 11. -С. 26-34.
6. Gershman A.B., Nickel U., Buhme J.F. Adaptive beamforming algorithms with robustness against jammer motion // IEEE Trans. Signal Processing. -1997. - Vol. 4. - P. 1878-1885.
7. Riba J., Goldberg J., Vazquez G. Robust beamforming for interference rejection in mobile communications // IEEE Trans. Signal Processing. -1997. - Vol. 45. - P. 271-275.
8. Ширман Я.Д., Манжос В.Н. Теория и техника обработки радиолокационной информации на фоне помех.- М.: Радио и связь, 1981. - 416 с.
9. Коростелев А.А. Пространственно-временная теория радиосистем. - Радио и связь, 1987. - 320 с.
10.0бработка сигналов в радиотехнических системах / под ред. А.П. Лукошкина. - Л.: изд-во Ленингр. ун-та, 1987. - 400 с.
11.Монзинго Р.А., Миллер Т.У. Адаптивные антенные решетки. Введение в теорию / пер. с англ.
Поступила 20 мая 2015 г.
под ред. В.А. Лексаченко. - М.: Радио и связь, 1986.- 446 с.
12.Черемисин О.П., Ратынский М.В., Комов А.А., Пушин А.Е. Эффективный проекционный алгоритм адаптивной
пространственной фильтрации // Радиотехника и электроника. - 1994. - Т. 39, №2. - С. 259-263.
13.Ратынский М.В. Адаптация и сверхразрешение в адаптивных антенных решетках.
- М.: Радио и связь, 2003. - 200 с.
14.Григорьев Л.Н. Цифровое формирование диаграммы направленности в фазированных антенных решетках. - М.: Радиотехника, 2010. -144 с.
15.Van Trees Harry L. Optimum Array Processing. Part IV of Detection, Estimation, and Modulation Theory. - New York: JOHN WILEY & SONS. - 2002.
- 1470 pp.
16.Алексеев С.Д. Цифровое диаграммообразование в фазированной антенной решетке с использованием сигнальных процессоров // Цифровая обработка сигналов. - 2002. - № 2. -С. 54-55.
17.Шмачилин П.А. Характеристики направленности антенных решеток при цифровой обработке сигналов на несущей // Антенны. - 2011.
- №3. - С. 47-56.
18. Журавлев А.К., Лукошкин А.П., Поддубный С.С. Обработка сигналов в адаптивных антенных решетках. - Л.: Изд-во Ленингр. ун-та, 1983. - 240 с.
English
Phased antenna array nulling in case of active jammer re-location
Ekaterina Sergeyevna Shtrunova - Candidate of Engineering, Associate Professor "Ryazan State Radio Engineering University". E-mail: shtrunova.e.s@rsreu.ru. Address: 390005, Ryazan, Gagarin St., 59/1.
Abstract: Angular motion of active jammer at the time interval between adaptation stages leads to rerouting of direction of active jamming in the environment of radio-electronic counteraction. In regard to this the jamming leaves area of earlier nulled beam pattern of receiving phased antenna array of the radar-location system. The article gives evaluation of reduction of signal-to-noise ratio improvement factor -(jamming+noise) optimum algorithm of direct inversion of selective interchannel correlation jamming matrix with direction cosines incrementation of active jamming between adaptation stages in comparison with passive jamming. Active jamming in near pattern lobe of phased antenna array when its angular direction approaches area of the main lobe has the greatest impact on efficiency of algorithms of adaptive processing at angular motion rate of active jammer. To reduce the impact of jammer angular rate it is expedient to minimize time between adaptation intervals so that it could be possible to neglect change of jamming angular position or to adjust pattern null position depending on jammer angular rate in view of maximizing signal-to-jamming ratio- (jamming+noise).
Key words: active noise jamming, adaptation of phased antenna array, jammer angular rate.
References
1. Dobykin V.D., Kupriyanov A.I., Ponomarev V. G., Shustov L.N. Radio-electronic warfare. Digital storing and reproduction of radio signals and electromagnetic waves. - Ed.by A. I. Kupriyanov. - M: Vuzovskaya kniga, 2009. - 360 p:
2. Radio-electronic suppression of data channels of weapon control systems. - Ed.by Yu.M. Perunov. - M.: Radiotekhnika, 2008. - 416 p.
3. Mailloux R.J. Covariance matrix augmentation to produce adaptive array pattern troughs. - Electron. Lett. -1995. - Vol. 31. - P. 771-772.
4. Zatman M. Production of adaptive array troughs by dispersion synthesis. - Electron. Lett. - 1995. - Vol. 31, No. 25. - P. 2141-2142.
5. Nechayev Yu.B., Borisov D. N., Peshkov I.V. Algorithm of pattern beam forming of ring antenna array, resistant to multibeam transmission and to jamming non-stationary sources. - Izvestiya vuzov. Radioelektronika. -2011. - T. 54, No. 11. - P. 26-34.
6. Gershman A.B., Nickel U., Buhme J.F. Adaptive beamforming algorithms with robustness against jammer motion. - IEEE Trans. Signal Processing. - 1997. - Vol. 4. - P. 1878-1885.
7. Riba J., Goldberg J., Vazquez G. Robust beamforming for interference rejection in mobile communications. -IEEE Trans. Signal Processing. - 1997. - Vol. 45. - P. 271-275.
8. Shirman Ya.D., Manzhos V. N. Theory and technology of radar data processing amid jamming. - M.:Radio i svyaz, 1981. - 416 p.
9. Kupriyanov A.I., Sakharov A.V. Basic theory of radio-electronic warfare. - M.: Vuzovskaya kniga, 2007. -356 p.
10. Korostelev A.A. Spatially-temporal theory of radio systems. - Radio i svyaz, 1987. - 320 p.
11. Processing of signals in radio technical systems. - Ed.by A.P. Lukoshkin. - L.: izdatelstvovo Leningradskogo universiteta, 1987. - 400 p.
12. Monzingo R. A., Miller T.U. Adaptive antenna arrays. Introduction to theory. - Transl. from Engl. Ed.by V.A. Leksachenko. - M.: Radio i svyaz, 1986. - 446 p.
13. Cheremisin O. P., Ratynsky M. V., Komov A.A., Pushin A.E. Effective projective algorithm of adaptive spatial filtration. - Radiotekhnika i elektronika. - 1994. - T. 39, No. 2. - P. 259-263.
14. Ratynsky M. V. Adaptation and superresolution in adaptive antenna arrays. - M.: Radio i svyaz, 2003. - 200 p.
15. Grigoriev L.N. Digital pattern beam forming in phased antenna arrays. - M.: Radiotekhnika, 2010. - 144 p.
16. Van Trees Harry L. Optimum Array Processing. Part IV of Detection, Estimation, and Modulation Theory. -New York: JOHN WILEY & SONS. - 2002. - 1470 pp.
17. Alekseev S. D. A Digital pattern beam forming in phased antenna array with using signal processors. -Tsifrovaya obrabotka signalov. - 2002. - No. 2. - P. 54-55.
18. Shmachilin P. A. Characteristics of antenna arrays' directivity with digital processing of carrier's signals. -Antenny. - 2011. - No. 3. - P. 47-56.
19. Zhuravlev A.K., Lukoshkin A.P., Poddubny S. S. Signals processing in adaptive antenna arrays. - L.: izdatelstvovo Leningradskogo universiteta, 1983. - 240 p.
