СПОСОБ АВТОМАТИЧЕСКОЙ СТАБИЛИЗАЦИИ ДИАГРАММЫ НАПРАВЛЕННОСТИ АНТЕННЫ ПРИ АДАПТИВНОЙ КОМПЕНСАЦИИ ПОМЕХ Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

https://doi.org/10.38013/2542-0542-2022-4-42-51 УДК 621.396.969

Способ автоматической стабилизации диаграммы направленности антенны при адаптивной компенсации помех

П. Е. Шимкин1,2

1 Акционерное общество «Всероссийский научно-исследовательский институт радиотехники», Москва, Российская Федерация

2 Национальный исследовательский университет «Московский энергетический институт», Москва, Российская Федерация

Предложено применимое на практике решение по реализации адаптивного компенсатора (АКП) активных шумовых помех (АШП) с автоматической стабилизацией формы диаграммы направленности антенны (ДНА), а также сниженным числом обусловленности промежуточных матричных вычислений. Приведено сравнительное исследование влияния дискретной и адаптивной регуляризации помеховой ковариационной матрицы или исходной выборки на параметры ДНА при работе радиолокационной станции с АКП АШП. Путем применения методов имитационного статистического моделирования и программирования подтверждена эффективность предложенного решения.

Ключевые слова: адаптивный компенсатор активных шумовых помех, регуляризация помеховой ковариационной матрицы, стабилизация формы диаграммы направленной антенны, снижение числа обусловленности

Для цитирования: Шимкин П. Е. Способ автоматической стабилизации диаграммы направленности антенны при адаптивной компенсации помех // Вестник Концерна ВКО «Алмаз - Антей». 2022. № 4. С. 42-51. https://doi.org/10.38013/2542-0542-2022-4-42-51

For citation: Shimkin P. E. A method of automatic stabilization of an antenna radiation pattern with adaptive jammer cancellation // Vestnik Koncerna VKO "Almaz - Antey". 2022. No. 4. P. 42-51. https://doi.org/10.38013/2542-0542-2022-4-42-51

Поступила 27.11.2022 Отрецензирована 14.12.2022 Одобрена 23.12.2022 Опубликована 17.01.2023

см см о см

< I

(0 те

s

о со

о.

ф

о

о

V

со

см ■ci-io

9 см ■ci-

10 см

Введение

Теоретические основы наиболее известных [16] методов и алгоритмов адаптивной компенсации (АКП) активных шумовых помех (АШП) в импульсно-доплеровских радиолокационных станциях (РЛС) с фазированной антенной решеткой (ФАР) оперируют информацией, содержащейся в помеховой ковариационной матрице (ПКМ) и предполагают точное ее знание. В практических случаях ПКМ априорно неизвестна и в качестве нее используется максимально правдоподобная апостериорная оценка, рассчитываемая по конечной выборке данных с вспомогательных приемных каналов (ВПК).

- © Шимкин П. Е., 2022

Точность оценки ПКМ, связанная с ограничением размера выборки, оказывает влияние на параметры адаптивной ДНА основного приемного канала (ОПК) при дальнейшем ее применении в алгоритмах АКП АШП.

Показателями качества АКП АШП при работе РЛС, помимо значений коэффициента подавления в направлении АШП и отношения сигнал/(помеха+шум) (ОСПШ) на выходе системы АКП АШП, являются также значения параметров адаптивной ДНА ОПК, динамически меняющиеся при обновлении оценки ПКМ и дальнейшем ее применении в алгоритме АКП АШП, а именно: повышение уровня боковых лепестков (УБЛ) и смещение положения максимума относительно исходной ДНА ОПК.

Физической основой дестабилизации вышеуказанных параметров при работе РЛС с включенной системой АКП АШП является разброс шумовых собственных значений ПКМ. Для его уменьшения в [7] предложен метод дискретной регуляризации ПКМ, заключающийся в аддитивной добавке к диагональным элементам оцененной ПКМ некоторого константного значения, сравнимого с уровнем собственного шума в пространственных каналах РЛС. В [3, 4] в качестве частных рекомендаций по выбору вышеуказанной константы предложены ее значения, находящиеся в диапазоне от 1 до 10 от мощности собственного шума в ВПК и ОПК. Работы [8-10] посвящены исследованию влияния значения дискретного коэффициента регуляризации на ОСПШ на выходе системы АКП АШП, а также выводу математических соотношений, теоретически позволяющих определить ограничения величины данного параметра.

В настоящее время интерес представляет исследование и предложение способов и алгоритмов расчета адаптивного значения коэффициента регуляризации ПКМ. Так, в работах [11-14] предложены способы обратного восстановления ПКМ, а в [15-17] исследованы оптимальные алгоритмы, использующие

адаптивную регуляризацию ПКМ на основании вычисления шумовых собственных значений ПКМ.

Несмотря на то что указанные исследования решают указанную научную проблему, их практическая реализация отличается высокой алгоритмической сложностью, например требует точного вычисления собственных чисел ПКМ, и все еще не позволяет в полной мере реализовать адаптивную регуляризацию ПКМ в реальном масштабе времени.

Настоящая работа посвящена исследованию алгоритмически эффективного с практической точки зрения способа вычисления адаптивного коэффициента регуляризации ПКМ, впервые предложенного в работе [18], для его применения к проблеме автоматической стабилизации ДНА ОПК.

Постановка задачи

Обобщенная структурная схема АКП АШП с применением линейной М-элементной ФАР для формирования ОПК и ^-элементной ФАР для формирования N ВПК изображена на рисунке 1.

На каждый элемент ФАР в к-й момент времени принимается аддитивная смесь

VV - V

VV - V

▼ 7акп

Рис. 1. Структурная схема АКП АШП с ОПК и N ВПК

та

X Ф

ч

та 0-

та

О

О.

£

V

ц

(Ч

см см о см

< I

со та

г

о со

о.

<и

о

о <и со

см ■ч-ю

с?

см ■ч-ю см

(П (П

сигнала от цели с направления ац, помехи от Ь источников АШП с направлений ац и шума (1):

XI

(к) = Хц(к) + Хц(к) + ХШ(к), ь

Хц(к) = s(k) а(ац ); Хц(к) = X Ь (к) а(ац);

м (1)

ХШ(к) = п(к),

а(а) =

;2я—вша

\,е х

, е

вша

•УоПК(к) = ЬОПК х(к); хВЦК(к) = СВЦК x(k),

(2)

_УАКП(к) = _УОПК(к) - _УВПК(к); _УВПК(к) = ^ЕШК хВПК(кХ

(3)

где _уАКП(к) - данные в к-й момент времени с выхода АКП АШП; ■№ВПК = [^...^]Т е С^1 -вектор весовых коэффициентов N ВПК.

Оптимальные весовые коэффициенты N ВПК в (3) определяются путем максимизации ОСПШ в ОПК на выходе АКП АШП. Для этого одним из эффективных и распространенных решений является применение метода с непосредственным обращением КМ (НОКМ) (4):

где х(к) = [х^к)... хЛ7(к)]т - выборка данных в к-й момент времени с N элементов ФАР; Хц(к), хц(к), хШ(к) - выборка данных в к-й момент времени сигнала от цели, помехи и шума; s(k), /(к), п(к) - огибающие в к-й момент времени формы сигнала, помехи и шума; а(а) - управляющий вектор ФАР в направлении а; (•)Т -оператор транспонирования.

Далее (1) в процессе формирования ОПК и N ВПК преобразуется в (2):

w

НОКМ

= та-1 гн°км (4)

ВПК ЛВПК жОВПК 5 ги

где Я.ВЦК е CNхN - ковариационная матрица (КМ) N ВПК; КВЦК е С^ - обратная матрица к КМ; г™ = [r1...rN]T е С^х1 - вектор взаимной корреляции ОПК и N ВПК.

В предположении некоррелированности сигнала от цели, помех и шума в (1) и с учетом требований к (2) КМ, фигурирующая в (4), представляется суммой отдельных составляющих и содержит информацию только о мощности и направлениях на Ь АШП и мощности шума (5), т.е. характеризует ПКМ:

где _УоПК(к) - данные в к-й момент времени с выхода схемы формирования ОПК; ЬОЦК = [Ъ1...ЪМ]Т - вектор для формирования ОПК;

хВПК(к) = [хВПК1(к). .. ХВПКдг(к)]Т - выборка данных в к-й момент времени с выхода схемы формирования N ВПК; СВЦК = [с1...сЛ?] - матрица для формирования N ВПК; сВЦК = [с1...сА\Т -вектор для формирования одного из ВПК; (•)н -оператор сопряжения по Эрмиту.

Как правило, матрица для формирования N ВПК в (2) рассчитывается с учетом следующих требований: 1) максимального перекрытия боковых лепестков ДНА ОПК в целях сохранения эффективности АКП с множественных направлений АШП; 2) минимально возможного провала ДНА N ВПК в надире для минимизации отклика N ВПК в направлении на цель в случае присутствия эхосигнала от цели в выборке данных для оценки ПКМ.

В общем случае решением задачи АКП АШП является взвешенное суммирование N ВПК с дальнейшей их компенсацией из ОПК (3):

ЯВЦК СШКЯСВЦК ;

^(Нц + ЯШ)СВЦК;

Я

Яц + Яп + Я

Ш>

4 Е [хц (¿X (к)] = ФК К К);

Кп А (к)4 (*)] - Ь *КК К) =(5)

А

=Хо*яК)яЯК);

м

где Яц, Яц, е САхА - КМ сигнала от цели, АШП и шума; Е[] - оператор математического ожидания; оц - мощность сигнала от цели; оШ - мощность собственного шума; оЦ - мощность у'-го АШЦ; I е (САхА - единичная матрица.

В реальных приложениях ПКМ априорно неизвестна и в качестве нее используется максимально правдоподобная апостериорная оценка по К отсчетам выборки данных с N ВПК (6):

К

ВПК

1 К 'впк]—1>

Л А=1

ВПК

(к)хвпк(к)' (6)

где КВПк е €NхN - оценка ПКМ.

Аналогично (6) для расчета второго множителя в (4) используется (7):

Предлагаемое решение

Связь ПКМ и ее оценки (6) можно представить в виде (10):

¿ВПК = »ВПК + рв, (10)

где в - масштабирующая константа; D е -матрица случайных значений с нулевым средним и единичной дисперсией.

В (10) второе слагаемое является аддитивной ошибкой оценки ПКМ и представляет собой разницу между априорной и апостериорной ПКМ. Для нейтрализации влияния этой ошибки на значения математического ожидания вышеуказанных параметров ДНА ОПК в настоящее время широко применяется метод регуляризации, который для АКП АШП с НОКМ математически представим в виде (11):

ЬЦ^Ь^Ж + 51 = »впк + PD + 51, (11)

где 5 - коэффициент регуляризации; I е -единичная матрица.

Одним из способов вычисления оптимального значения параметра 5 в (11) при условии высоких значений ОПШ или при большом значении К, т.е. ||ЯВПК + 51|| >> Р||Б|| и 5 << КВПК(/,/), I е 1..Лявляется (12) [18-20]:

(м)]<5<<*(квпк)/лт, (12)

^овпк -•^[^пк.Уопк] _ ~ 1>ВПК (^).Уопк (7)

А к=\

Наиболее критичной является оценка (6), т.к. именно она участвует в операции обращения ПКМ, приводит к квазиоптимальному решению задачи по (4) и оказывает влияние на параметры адаптивной ДНА ОПК.

С точки зрения ОСПШ на выходе системы АКП АШП известна [1-6] зависимость, позволяющая определить математическое ожидание потерь в связи с оценкой ПКМ (8):

= (8)

Так, для ограничения потерь в ОСПШ ОПК после АКП АШП на уровне 0,5 (-3 дБ) необходимо и достаточно провести оценку ПКМ (5) по К = 2N - 3 отсчетам выборки с N ВПК. В практических случаях К выбирается в 10 и более раз превышающем количество ВПК, тем самым ограничивая потери в ОСПШ ОПК после АКП на уровне -0,45 дБ и менее.

По показателю УБЛ ДНА ОПК при работе АКП АШП известна [1-6] зависимость, обратно пропорционально зависящая от К (9):

¿[УБЛ^]^. (9)

При современных требованиях к УБЛ ДНА ОПК на уровне -40 дБ и менее согласно (9) потребуется выборка данных с N ВПК для оценки ПКМ, превосходящая ~104 отсчетов. При работе РЛС в режиме среднего и быстрого темпа обзора пространства с расчетом весовых коэффициентов АКП АШП в каждом периоде зондирования РЛС выполнение данного требования практически нереализуемо, что задает противоречие между параметрами адаптивной ДНА ОПК при работе АКП АШП и заданной кинематикой обзора РЛС.

По показателю смещения максимума ДНА ОПК при работе АКП АШП в настоящее время отсутствует анализ и точное либо приближенное выражение для расчета, в связи с чем в данной работе будут приведены результаты имитационного моделирования с указанием относительных значений данного показателя.

где о(-) - среднеквадратическое отклонение (СКО); RВпк(/,0, i е 1...N- диагональные элементы матрицы RBnK; tr^BnK) - след матрицы

»впк.

Следует также отметить, что при реализации метода НОКМ в аппаратуре с ограничением разрядности вычислений, например в формате плавающей точки одинарной точности (IEEE 754 32 бита), необходимо перейти к вы- _ числению (4) согласно модифицированному | методу обращения КМ (МОКМ) [3, 4] (13), (14): к

МОКМ _-р-1 МОКМ. 5

wBm _ говпк ' (13) ||

те

МОКМ _f\H /1/14 ^

овпк -VlJW* (14) 2

где Q е CKxK, R е CKxN - соответствующие | компоненты результата QR-разложения вы- ё борки x^; R1 = R(1:N,:) е CNxN; Qi = Q (:,1N) J

Хвпк 1 К дкм а, о РЕГ ХрЕГ РАЗ 01 ОТМ И"1 ВЕС ^впк КОМ

1

3, Сщ

АРЕГ

Уопк

АКП

АКП с МОКМ

Хвпк

Хвпк

РКВ

Говпк

Уопк

к

Уопк, АКП

Уакп

Рис. 2. Структурная схема АКП АШП с МОКМ и регуляризацией

см см о см

< I

со те

г

I

о ^

со

о.

<и

о

о <и со

см ■ч-ю

с?

см ■ч-ю см

(П (П

е - матрицы сокращенной размерности от Р и Я соответственно; » - треугольная

МОКМ /Г1№1

матрица; г0впк е (ъ - результат произведения и .Уопк.

Как будет показано далее по результатам моделирования, такой переход связан с возрастающим числом обусловленности матриц, участвующих в промежуточных вычислениях (4) в исходном виде при высоких значениях ОПШ на входе АКП АШП, влияющим на значение ОСПШ на выходе АКП АШП.

Согласно (13-14) основным отличием МОКМ от НОКМ является наличие операции ортогонализации исходной выборки данных с N ВПК с помощью QR-разложения, позволяющей снизить число обусловленности матриц, участвующих в дальнейших вычислениях.

При этом для регуляризации исходной выборки данных с N ВПК при реализации АКП

10 0 -10 -20 -30 -40 -50 -60 -70 -80

дБ

-80 -60 -40 -20 0 20 40 60 80

а, °

Рис. 3. ДНА ОПК и ВПК ■ - ОПК;--ВПК 1;--ВПК 2;--ВПК 3

АШП с МОКМ аналогичное (11) будет определяться (15) [4]:

хвт=[хвпк; 8^1], (15)

где [;] - операция конкатенации; I е -единичная матрица.

Тогда схему устройства АКП АШП с МОКМ и регуляризацией можно представить структурой, изображенной на рисунке 2, на котором приняты следующие условные обозначения: ДКМ - устройство вычисления диагональных элементов КМ согласно (6); РЕГ - устройство регуляризации исходной выборки данных с N ВПК согласно (15); РАЗ - устройство вычисления QR-разложения исходной выборки данных с N ВПК; ОТМ - устройство вычисления обратной треугольной матрицы в (13); РКВ - устройство вычисления корреляционного вектора согласно (14); ВЕС - устройство вычисления весовых коэффициентов согласно (8); КОМ - устройство компенсации согласно (3); АРЕГ - признак адаптивной регуляризации согласно (12) (АРЕГ = 1/0 - включена/выключена адаптивная регуляризация; АРЕГ = 0 - включена дискретная регуляризация с фиксированным коэффициентом регуляризации, при этом если 5 = 0, то регуляризация отсутствует); АКП - признак применения рассчитанных весовых коэффициентов (АКП = 1/0 компенсация помех включена/выключена).

Результаты моделирования

В качестве частных параметров моделирования используются следующие значения: М=А = 14; N = 6; бТк = 0,8. дНа ОПК и 6 ВПК (3 попарно одинаковых) показаны на рисунке 3.

Решается одномерная задача АКП в азимутальном направлении нормально распределенной АШП с использованием фиксированной выборки по K = 128 отсчетам данных с 6 ВПК, полученных в конце рабочей дистанции каждого периода зондирования РЛС; количество испытаний Монте-Карло T = 103.

Качественные и количественные результаты моделирования АКП АШП с НОКМ и МОКМ при вычислениях в арифметике IEEE 754 и варьируемом значении ОПШ на входе АКП АШП в направлении 1 АШП показаны на рисунке 4.

Как следует из рисунка 4, при ОПШ на входе АКП АШП более 50 дБ наблюдается деградация ОСПШ на выходе АКП АШП, качественно проявляющаяся в недостаточном

провале ДНА в направлении на АТТТП Этот факт связан с тем, что при вычислении по методу НОКМ число обусловленности ПКМ при относительно высоком ОПШ или высоком отношении максимального к минимальному собственному значению ПКМ превосходит ~104, рисунок 5 а. При прочих равных условиях, в связи с отсутствием вычисления ПКМ в прямой форме, вычисление по методу МОКМ оперирует матрицами с числом обусловленности ~102, т.е. в корень раз меньше, что с используемым форматом арифметики вычислений является достаточным для сохранения линейности зависимости, показанной на рисунке 5б.

Исследуемыми показателями качества АКП АШП, как упоминалось ранее

10 0 -10 -20 -30 -40 -50 -60 -70

дБ

10 0 -10 -20 -30 -40 -50 -60 -70 -80

дБ

20 40 60 80

-80 -60 -40

-60 -40 -20 0

а, ° а, °

аб Рис. 4. ДНА ОПК до и после применения АКП АШП с НОКМ и МОКМ при ОПШ ~55 дБ:

а) исходный масштаб; б) увеличенный масштаб в направлении АШП — до АКП;--после АКП с НОКМ;--после АКП с МОКМ;----направление АШП

10 20 30 40 50 ОПШ до АКП, дБ

а

60 70 80

£

<

80 70 60 50

5 40 о

5 30 ^ 20

10

0 -10

Лг

10 20 30 40 50 ОПШ до АКП, дБ б

60 70 80

Рис. 5. Зависимости от ОПШ на входе АКП: а) значения числа обусловленности; б) значения ОСПШ на выходе

АКП АШП --АКП с НОКМ; - АКП с МОКМ

та

X Ф

ч

та Q.

та

о

о.

£

ф ц

(Ч

0

0

и как можно заметить из рисунка 4а, являются следующие параметры ДНА ОПК при работе РЛС с АКП АПШ: повышение УБЛ и смещение положения максимума ДНА ОПК относительно исходной ДНА ОПК (16):

1 т I

|Апм| = -Хкё^СДНАопк (0]"

Т

-агётах[дНА^(0]|.

делирования, полученные путем применения адаптивной регуляризации исходной выборки данных согласно левой крайней границе в (12), дискретной регуляризации с коэффициентом регуляризации согласно (17) и без регуляризации (5 = 0):

(16)

Качественное влияние регуляризации исходной выборки данных для заданных параметров статистического моделирования изображено на рисунке 6.

Для количественного сравнения на рисунке 7 в виде соответствующих кривых и штриховых линий показаны результаты мо-

Ош < 5 < коп

(17)

где ош - максимальное значение мощности собственного шума среди N ВПК и ОПК; к -масштабирующая константа при дискретной регуляризации.

Кроме того, интерес представляет влияние регуляризации исходной выборки на ОСПШ при работе РЛС в динамике. Так, в РЛС с механическим сканированием по азимуту со скоростью вращения антенны,

см см о см

< I

со та

г

I

о ^

со

.

<и

о

о <и со

см ■ч-ю

с?

см ■ч-ю см

(П (П

10 0 -10 -20 -30 -40 -50 -60 -70

дБ

/ \

\

а 'л

лА' i \

/ А А Г "М1 ш т г

VI V

1

-60 -40 -20 0 20 40 60 80

О

10 0 -10 -20 -30 -40 -50 -60 -70 -80

дБ

-80 -60 -40 -20 0 20 40 60 80

О

аб

Рис. 6. ДНА ОПК до и после включения алгоритма АКП АШП: а) без применения регуляризации;

б) с применением регуляризации --до АКП;--после АКП без РЕГ;--после АКП с РЕГ;----направление АШП

10

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 8/а

10-1

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 8/а

ш

б

Рис. 7. Зависимости параметров ДНА ОПК до и после АКП АШП от значения коэффициента регуляризации 5: а) смещение положения максимума ДНА ОПК; б) увеличение УБЛ ДНА ОПК --АКП с РЕГ по (17);----АКП без РЕГ (5 = 0);----АКП с РЕГ по (12)

а

300

200

й

i 0

<

102

£

101

(-100

-200

О

о 100

и S ft

^ , ¿3 10-1

10-

iziiiiiTizi

V

0,2 0,4

0,6 0,8 1 1,2 t, мс

а

1,1 1,2 1,3 1,4 1,5 1,6 1,7 1,8 1,9 Период обновления коэффициентов, мс

б

Рис. 8. Результаты работы АКП АТТТП «без» и «с» регуляризацией при вращении антенны РЛС: а) временное сечение; б) зависимость потерь в ОСПШ на выходе АКП АТТТП от периода обновления коэффициентов

--до АКП;--после АКП без РЕГ;--после АКП с РЕГ;

--10 об./мин. АКП без РЕГ;--20 об./мин. АКП без РЕГ;--10 об./мин. АКП с РЕГ;

- 20 об./мин. АКП с РЕГ

например от 10 до 20 оборотов в минуту, и дискретной оценкой весовых коэффициентов АКП АТТТП в конце рабочей дистанции РЛС, например один раз в период от 1 до 2 мс, будет наблюдаться снижение ОСПТТТ на выходе АКП АТП, возрастающее к началу рабочей дистанции из-за остаточной «не скомпенсированной» в АКП АТТТП в связи с угловым перемещением постановщика помех по ДНА ОПК. Качественное и количественное проявление данного эффекта во временной области показано на рисунке 8.

Как следует из рисунков 6-8, применение метода регуляризации с предложенным способом вычисления адаптивного коэффициента регуляризации позволяет решить задачу автоматической стабилизации ДНА ОПК как в статической, так и в динамической работе РЛС с АКП АТП.

Заключение

Для рассмотренного в данной статье примера адаптивная регуляризация при работе РЛС с АКП АТТТП позволит снизить математические ожидания: а) смещения положения максимума ДНА ОПК от ~0,3° до ~0,09°; б) увеличения УБЛ ДНА ОПК с ~6 до ~0,3 дБ.

В то же время дискретная регуляризация при оптимальном, с точки зрения апостериорного минимума смещения положения максимума ДНА ОПК, значении коэффициента

регуляризации ПКМ позволит снизить математические ожидания: а) смещения положения максимума ДНА ОПК от ~ 0,3° до ~ 0,06°; б) увеличения УБЛ ДНА ОПК с ~6 до ~0,6 дБ при работе системы АКП АШП.

Значения параметров смещения положения максимума ДНА ОПК и УБЛ ДНА ОПК при адаптивной регуляризации являются пределами, к которым асимптотически стремятся аналогичные значения при дискретной регуляризации при работе РЛС с АКП АШП.

С практической точки зрения регуляризация с относительно простым с алгоритмической точки зрения способом вычисления адаптивного коэффициента регуляризации может эффективно реализоваться в существующих разработках без изменения аппаратной части устройства АКП.

Список литературы

1. Monzingo R. A., Miller T.W. Introduction to Adaptive Arrays. New York: Wiley, 1980.

2. Johnson D. H., Dudgeon D. E. Array Signal Processing: Concepts and Techniques, Prentice Hall, Englewood Cliffs. New York, 1993.

3. Van Trees H. L. Optimum Array Processing (Detection, Estimation, and Modulation Theory, Part IV). New York: Wiley, 2002.

4. Vinay K. I., Kogon M. S., Manolakis D. Statistical and Adaptive Signal Processing: Spectral Estimation, Signal Modeling,

те

X <D

4 те Q.

те

о

о.

£

<D (4

1

2

см см о см

< I

со те

s

о

CQ

Adaptive Filtering and Array Processing. Artech, 2005.

5. Melvin W. L., Sheer J. A. Principles of modern radar. Vol. II: Advanced Techniques. SciTech Publishing, 2013.

6. Ратынский М. В. Адаптация и сверхразрешение в антенных решетках. Изд. 2-е, доп. М.: ЛЕНАНД, 2022. 240 с.

7. Carlson B. D. Covariance matrix estimation errors and diagonal loading in adaptive arrays // IEEE Trans. Aerospace and Electronic Systems. 1988. No. 24(4). P. 397-401.

8. Черемисин О. П. Эффективность адаптивного алгоритма с регуляризацией выборочной корреляционной матрицы // Радиотехника и электроника. 1982. Т. 27. № 10. С. 1933-1941.

9. Черемисин О. П. О выборе параметра для регуляризованного метода адаптивной оптимизации фильтров // Радиотехника и электроника. 1985. Т. 30. № 12. С. 2369-2377.

10. Ратынский М. В. Выбор регуляризато-ра в задаче адаптивной пространственной фильтрации // Успехи современной радиоэлектроники. 2016. № 7. С. 53-63.

11. Gu Y., Leshem A. Robust adaptive beam-forming based on interference covariance matrix reconstruction and steering vector estimation // IEEE Transactions on Signal Processing. 2012. Vol. 60, no. 7. P. 3881-3885.

12. Yang X., Yan L., Sun Y., Zeng T. Improved orthogonal projection approach utilising interference covariance matrix reconstruction for adaptive beamforming // IEEE Electronics Letters. 2014. Vol. 50, no. 20. P. 1446-1447.

13. Lu Z., Li Y., Gao M., Zhang Y. Interference covariance matrix reconstruction via steering vectors estimation for robust adaptive beamforming //

IEEE Electronics Letters. 2013. Vol. 49, no. 22. P. 1373-1374.

14. Wu X., Zhu W., Yan J., Toeplitz A. Covariance Matrix Reconstruction Approach for Direction-of-Arrival Estimation // IEEE Trans. Veh. Technol. 2017. Vol. 66. P. 8223-8237.

15. Besson O., Vincent F. Performance analysis of beamformers using generalized loading of the co-variance matrix in the presence of random steering vector errors // IEEE Trans. Signal Process. 2005. Vol. 53. P. 452-459.

16. Gu J., Wolfe P. J. Robust adaptive beamforming using variable loading, in Proceedings of the 4th IEEE Sensor Array and Multichannel Signal Processing Workshop Proceedings. SAM, 2006. P. 1-5.

17. Du L., Li J., Stoica P. Fully automatic computation of diagonal loading levels for robust adaptive beamforming // IEEE Transactions on Aerospace and Electronic Systems. 2010. Vol. 46, no. 1. P. 449-458.

18. Ning Ma, Joo Thiam Goh. Efficient method to determine diagonal loading value // IEEE International Conference on Acoustics, Speech, and Signal Processing. 2003. Proceedings. (ICASSP '03). P. V-341. DOI: 10.1109/ ICASSP.2003.1199948

19. Xiao Yuteng, Yin Jihang, Qi Honggang, Yin Hongsheng, Hua Gang. MVDR Algorithm Based on Estimated Diagonal Loading for Beamforming // Mathematical Problems in Engineering. 2017. P. 1-7.

20. Muhammad M., Li M., Abbasi Q., Goh C., Imran M. A. A Covariance Matrix Reconstruction Approach for Single Snapshot Direction of Arrival Estimation. Sensors. 2022. Vol. 22, no. 8. P. 3096.

Q.

<D

О

о

<D CQ

CM ■Clio 9

CM ■Clio

CM

Об авторе

Шимкин Павел Евгеньевич - кандидат технических наук, заместитель начальника отдела Акционерного общества «Всероссийский научно-исследовательский институт радиотехники», Москва, Российская Федерация; доцент кафедры радиотехнических приборов и антенных систем Национального исследовательского университета «Московский энергетический институт», Москва, Российская Федерация.

Область научных интересов: разработка методов и алгоритмов цифровой обработки сигналов в радиолокационных системах и комплексах.

(П (П

A method of automatic stabilization of an antenna radiation pattern with adaptive jammer cancellation

Shimkin P. E.1'2

1 JSC "All-Russian Scientific and Research Institute of Radio Engineering" Moscow, Russian Federation

2 National Research University "Moscow Power Engineering Institute", Moscow, Russian Federation

The paper describes an applicable solution that allows to implement an adaptive jammer canceller (AJC) with automatic stabilization of the antenna radiation pattern (ARP) shape along with the reduced condition number of intermediate matrix computations. A comparative analysis is conducted to investigate how discrete and adaptive regularization of the interference covariance matrix or the initial sampling data may affect the ARP parameters during operation of a radar with the AJC. The proposed solution proves efficient by applying the statistical simulation methods and programming.

Keywords: adaptive jammer canceller, regularization of interference covariance matrix, stabilization of antenna radiation pattern shape, reducing the condition number

Information about the author

Shimkin Pavel Evgenievich - Candidate of Engineering Sciences, Deputy Head of Department, JSC "All-Russian Scientific and Research Institute of Radio Engineering"; Associate Professor of Department of Radio Devices and Antenna Systems, National Research University "Moscow Power Engineering Institute", Moscow, Russian Federation. Science research interests: development of methods and algorithms of digital signal processing in radar systems.