Научная статья на тему 'Повышение эффективности подавления активных помех в РЛС перехвата воздушных целей в режиме обзора'

Повышение эффективности подавления активных помех в РЛС перехвата воздушных целей в режиме обзора Текст научной статьи по специальности «Электротехника, электронная техника, информационные технологии»

CC BY
1172
205
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
АКТИВНАЯ ШУМОВАЯ ПОМЕХА / АДАПТИВНАЯ КОМПЕНСАЦИЯ ПОМЕХ / КОВАРИАЦИОННАЯ МАТРИЦА ПОМЕХ / ВЕКТОР ВЕСОВЫХ КОЭФФИЦИЕНТОВ / МЕТОД ПОСЛЕДОВАТЕЛЬНОЙ РЕГРЕССИИ

Аннотация научной статьи по электротехнике, электронной технике, информационным технологиям, автор научной работы — Маркевич В.Э.

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

i Надоели баннеры? Вы всегда можете отключить рекламу.

Похожие темы научных работ по электротехнике, электронной технике, информационным технологиям , автор научной работы — Маркевич В.Э.

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Increase of efficiency of suppression of the active noise jamming in radar interception of air targets in the review mode

The modified algorithm of the adaptive compensating of the active noise jamming, based on usage of the recurrence estimation of inverse covariance matrix of jamming in the course of an self-tuning is considered. The matrix estimation is formed by a method of a serial regression (SRM), a consequence application of a lemma about inversing of a matrix to differential to the difference equation for a rounded estimation of direct covariance matrix of jamming (CMJ). The single estimation of straight line CMJ is formed as result dyad products of the instant vectors of noises in handling channels. In partial channels the adaptive transversal filters with the self-tuning on value of a signal of an aggregate error are implemented. In article skeleton diagrams of devices of the compensating implemented by various methods are resulted. Mathematical modeling of operation of the algorithms, implemented by a method of the least squares (criterion of a minimum RMS error) and a method of a serial regression is performed. Results of mathematical modeling of operation of the device of compensating of active noises jamming in the conditions of the active radio-electronic counteraction at multipath propagation of signals both non-identical amplitude and phase characteristics of channels of handling are resulted. Optimal values of technical parameters of devices of compensating of the active noises jamming providing obtaining of demanded values of coefficient of suppression are given. The comparative analysis of indexes of efficiency of classical and offered algorithms of compensating is made. Advantage of algorithm of compensating of the active noises jamming, using the recurrence estimation of a inverse covariance matrix of jamming is shown. The algorithm of functioning on the basis of SRM allows to lower essentially requirements to specialized VLSI (FPGA) at implementation is hardware program complexes for preprocessing of the radar, radio communication and radio navigation information in the conditions of influence of the active noises jamming in the presence of effects of multipath propagation of signals, in the presence of delay of signals on an antenna system aperture, and also at unremovable technological dispersion amplitude and phase characteristics of the microwave of channels of handling.

Текст научной работы на тему «Повышение эффективности подавления активных помех в РЛС перехвата воздушных целей в режиме обзора»

УДК 621.396

В. Э. МАРКЕВИЧ

ПОВЫШЕНИЕ ЭФФЕКТИВНОСТИ ПОДАВЛЕНИЯ АКТИВНЫХ ПОМЕХ В РЛС ПЕРЕХВАТА ВОЗДУШНЫХ ЦЕЛЕЙ

В РЕЖИМЕ ОБЗОРА

оао «алевкурп»

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

В статье приведены структурные схемы устройств компенсации, реализованных различными способами. Выполнено математическое моделирование работы алгоритмов, реализованного методом наименьших квадратов (критерий минимума среднеквадратической ошибки) и методом последовательной регрессии.

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

Алгоритм функционирования на основе МПР позволяет существенно снизить требования к специализиро-ванным СБИС (ПЛИС) при реализации аппаратно-программных комплексов для первичной обработки радиолокационной, радиосвязной и радионавигационной информации в условиях воздействия активных шумовых помех при наличии эффектов многолучевого распространения сигналов, при наличии запаздывания сигналов по раскрыву антенной системы, а также при неустранимом технологическом разбросе амплитудно-фазовых характеристик СВЧтракта каналов обработки.

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

Введение

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

1. Высокий уровень мощности АТТТП в основном и дополнительных каналах обработки

по сравнению с мощностью сигнала и уровнем внутренних шумов приемного устройства, в том числе - при приеме по дальним боковым лепесткам.

2. Расстояние между дополнительными приемными (компенсационными) антеннами удовлетворяет условию разделения обработки на пространственную и временную (межэлементное запаздывание сигналов по апертуре антенны пренебрежимо мало);

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

4. Многолучевое распространение АШП, обусловленное наличием на трассе участков

локального отражения и приводящее к искажению волнового фронта помехи на апертуре антенных систем, отсутствует;

5. Угловые направления помехи находятся за пределами основного лепестка диаграммы направленности антенной системы (ДНА АС) основного канала (ОК), число активных помех, воздействующих с различных угловых направлений, не превышает число антенных элементов.

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

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

В предлагаемой статье рассмотрен алгоритм, обеспечивающий заданную эффективность работы многофункциональной РЛС в режиме обзора при разбросе амплитудно-фазовых характеристик тракта обработки и воздействия нескольких скользящих по частоте хаотических импульсных помех (ХИП) с учетом многолучевого распространения сигналов.

Исходные данные и начальные условия для моделирования

Величина коэффициента подавления определяется следующим образом:

Кп = 101ё

N-1 ^

1 -П

, ]=0 У

> К тр =

(1)

V Р1 У

= 101ё -2 =(20...30) дБ,

ный коэффициент корреляции; Рь Р2 - мощность АТТТП в основном канале обработки на входе и выходе АКМИ соответственно; J - число независимых факторов, определяющих величину коэффициент подавления.

Для интегральной оценки эффективности подавления используется коэффициент подавления АТТТП, усредненный на интервале обзора по угловым координатам, длительностью 62 мс:

( (К-1

Кп = 10

Л

^ £ Р2к

V к=0 У

(к-1 ^Л

- ^ £ г*

V к=0 УУ

(2)

Р1к = хк Хк, Р1к гк,

где Р1к, Р2к - усредненные за период повторения импульсов мощности АТТТП на входе и выходе АКМИ; Хк, £к - вектор отсчетов входного и выходного сигнала АКМИ в основном канале обработки для к-го периода повторения (ПП)

зондирующих импульсов; выходных К - об-

Н Н

щее число периодов повторения; Хк , £к - эр-митово-сопряженные вектора хк, гк.

Из соотношения оценим минимальное значение требуемого коэффициента корреляции г0, обеспечивающего заданный коэффициент подавления АКМИ:

г0 = (1 -10-01Кп )05,

(3)

где Кп, КГр - фактический и требуемый коэффициенты подавления помехи; г]- - парциаль-

Так, для обеспечения коэффициента подавления Кп = 20 дБ результирующий коэффициент межканальной корреляции должен составлять г0 = 0.995, а для Кп = 30 дБ г0 = 0.9995. Видно, что требования, предъявляемые к степени идентичности АТТТП в каналах обработки являются достаточно жесткими и требуют реализации специальных мер по обеспечению заданного коэффициента подавления.

Рассмотрим упрощенную структурную схему устройства АКМИ, представленную на рис. 1, представлена одним основным и четырьмя дополнительными идентичными пространственными каналами приема.

Основной пространственный канал обработки представлен приемо-передающей зеркальной антенной с электромеханическим сканированием и трактом обработки, включающем СВЧ-приемник, УПЧ, АЦП и блок цифровой обработки (БЦО). Компенсационные каналы представлены слабонаправленными антеннами в количестве N = 4, основной лепесток диаграммы направленности (ДНА) которых

Рис. 1. Упрощенная структурная схема устройства АКМИ, реализующего алгоритм минимизации СКО ошибки

оценивания

перекрывает сектор 15...17 градусов с центром в направлении оптической оси зеркальной антенны основного канала приема. Алгоритм адаптивной компенсации помех реализуется в БЦО в виде классической схемы с корреляционной обратной связью и самонастройкой [1-6].

Для адаптации к уровню мощности полезного сигнала используется модифицированный алгоритм самонастройки с дополнительным сигналом смещения, вводимого в контур обратной связи.

Дополнительная цепь отрицательной обратной связи с коэффициентом усиления а используется для исключения возможности подавления всех сигналов (включая помехи), уровень которых не превосходит определенной величины. В [1] показано, что за счет выбора параметра а можно исключить (минимизировать) подавление полезного сигнала и обеспечить одновременно компенсацию помех с высоким уровнем мощности. Поэтому выходное отношение сигнал - (помеха + шум) будет выше, чем в случае использования в контуре корреляционной обратной связи только интегратора.

Максимальная эффективность использования описанной схемы реализуется при выполнении условия:

Рп « Р, « Р}, (4)

где Рп - мощность внутренних шумов в приемных каналах; Р, - мощность полезного сигнала; PJ - мощность АТТТП в приемных каналах.

Алгоритм формирования ВВК при использовании дополнительного сигнала в цепи корреляционной обратной связи, будет рассмотрен ниже.

На структурной используются устройства АКМИ следующие обозначения:

х^) = ( х0 ^) х1 ^) х2 ^) х3 ^)) - вектор комплексных амплитуд сигналов на входах компенсационных каналов, образованный парциальными сигналами;; Уо(0 - сигнал на входе основного (опорного) канала обработки с зеркальной антенной и высокими направленными свойствами; а - скалярный множитель, используемый для регулировки порогового уровня для исключения подавления слабых сигналов и помех; ^ = (2-2..2-7) - постоянная скалярная величина (эквивалент постоянной времени

сглаживания), определяющая начальный размер шага и параметры фильтра для сглаженной оценки мощности АТП; е(?) - выходной сигнал устройства АКМИ после подавления помехи; N е (1, 4) - число дополнительных каналов приема и обработки.

Характеристики антенной системы основного канала:

- ширина основного лепестка ДНА ОК, град - 1.2;

- уровень ближних боковых лепестков, дБ -минус 23;

- угловая скорость обзора, град/с - 278;

Характеристики антенных систем компенсационных каналов:

- ширина основного лепестка ДНА ОК, град -15...17;

- уровень ближних боковых лепестков, дБ -минус 12;

- уровень кросс-поляризационной ДНА, дБ -минус 17;

Характеристики отраженного сигнала:

- длина волны зондирующего сигнала, м -0.032;

- длительность зондирующего импульса, мкс - 0.3;

- период повторения импульсов, мкс - 280;

Характеристики активной шумовой помехи:

- эффективная ширина спектра, МГц - 3.4;

- длительность импульса, мкс - 0.3.30;

- период повторения импульсов, мкс -1.100;

- спектральная плотность мощности АТТТП, кВт/МГц - 2.4;

Характеристики приемного устройства:

- частота дискретизации сигнала на выходе фильтра, МГц - 5.83(3);

- число импульсов в пределах цикла обзора по углам - 128;

- размер изображения на индикаторе дальность-угол, мхград - 42000x10.

Модель функционирования устройства адаптивной компенсации АШП на основе алгоритма минимимума среднеквадратической ошибки

Рассмотрим наиболее важный для практики случай, когда не выполняются условия оптимальной обработки № 3 и № 4:

- существуют взаимные амплитудно-фазовые, частотные и временные неидентичности

в основном и дополнительных парциальных каналах обработки сигнала;

- прием сигнала и помехи происходит в условиях многолучевого распространения из-за переотражения АТТТП от подстилающей поверхности на трассе распространения (например, при воздействии специальной АТТТП типа «антипод»).

В качестве исходного алгоритма расчета ВВК используется модифицированный градиентный алгоритм функционирования адаптивного фильтра пространственной режекции помех (с учетом дополнительного контура отрицательной обратной связи), синтезированного по критерию минимума средней квадратиче-ской ошибки (минимума СКО) [1-6, 10]:

) = М | 8 2(^ -№) -1 awH w

дw

(5)

= 0, 8^) = У0^) - wH ^)Х^)

где М(8 w)) - среднее значение квадрата функции 8(?,w); w - комплексный вектор весовых коэффициентов (ВВК); wH = w*T - эр-митово-сопряженный ВВК.

Работа контуров управления автоматической подстройкой весовых коэффициентов осуществляется в соответствии с дифференциальным уравнением:

^^ (Г )=(У „ (^, w ))-aw )=-2ц(

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

х* (Г) Рх (t)

8^) - aw

Рх (t) = у 0§РХ (t), 8РХ (t) = Рх (t) - Рх (t), Рх (t) =

хН(t)ха), У0 = Ц1 ((5РХ(t) < 0)? 2- : 1), / е (1...4),

(6)

где V№(Е,^, w)) - оценка градиента функции w); Рх (t) - оценка мгновенной мощности помехи в каналах обработки; w ц0 = (2-2...2-8) -постоянная скалярная величина (коэффициент сходимости), определяющая размер шага и регулирующая скорость сходимости и устойчивость процесса самонастройки АКМИ.

Дальнейшим развитием представленного на рис. 1 АКМИ является устройство адаптивной компенсации АТП с использованием КИХ-фильтров в каналах обработки.

Основные преимущества использования адаптивных КИХ - фильтров в устройствах пространственной компенсации АТП:

Рис. 2. Упрощенная структурная схема устройства АКМ, на основе адаптивных фильтров с КИХ в каналах обработки

- идентификация и выравнивание искажений импульсной характеристики трассы распространения сигнала (помехи), обусловленных наличием многолучевых каналов, в реальном масштабе времени;

- выравнивание амплитудно-фазовых характеристик каналов обработки, обусловленных различной электрической длиной трактов СВЧ и неидентичностью аналоговых фильтров СВЧ и УПЧ;

- компенсация задержки фронта сигнала, обусловленной распространением по раскры-ву малоэлементной решетки, образованной дополнительными антенной элементами.

В этом случае ВВК будет представлен в виде составного вектора, образованного набором из N парциальных векторов размерностью L (где N - число каналов, L - размерность импульсной характеристики адаптивного фильтра в каналах обработки):

/ t w = I w 0

w

*Т V

T

■ W N -

w =( w A W , ... Wr "n V n,0 n,i n, L-

i )T, i )T.

(7)

Выражение для сглаженной оценки мощности в этом случае примет вид:

px (t) = YiSpx (t), six (t) = px (t)- P (t),

l N-1 H

Px (t) = T7 Z XHn Xn, Px (t)=LM(Pmln, Pmax, (t)),

N n=0

Yi = ^i ((sPx (t) < 0) ? 2-l : l), L e (3...9), px(0) = x0HXo,

где Xn = (Xn,0 Xni, ... , Xn,L-i)H - вектор парциальных сигналов с выхода адаптивного КИХ -фильтра; LIM (хь х2, х) - функция ограничения входного аргумента x(t) значениями xi, х2; Pmln, Pmax - минимальное и максимальное значение

max

мощности АШП, используемое в системе обработки; Ад - интервал дискретизации входного сигнала (помехи); Px(t) - сглаженная оценка мощности помехи; Т00 - постоянная времени сходимости процесса самонастройки АКМИ; Тф - постоянная времени сглаживающего фильтра для мгновенной оценки мощности АШП.

На рис. 2 представлена упрощенная структурная схема АКМИ на базе адаптивных фильтров с КИХ в компенсационных каналах обработки. С целью упрощения на структурной схеме не показаны каналы фильтровой обработ-

Импульсная характеристика ФОС

8x10

6x10"

4x10

2x10

-2x10

Йк\

/ Г 1 \\\

4 'л// \

V Л ' л

10

20

30

Рис. 3. Импульсные характеристики ФОС с гауссовской амплитудно-частотной характеристикой в каналах обработки АКМИ

ки входных сигналов и помех. При этом полагается, что аддитивная смесь сигналов и помех, действующих на входе основного и компенсационных каналов, получена путем предварительной фильтрации в фильтре одиночного сигнала (ФОС) с гауссовской импульсной характеристикой (ИХ) в каждом канале обработки.

Моделировалась тактическая ситуация, в которой для прикрытия цели использовались три постановщика АШП-источника мешающих радиоизлучений (ИРИ) заданной интенсивности, генерирующих шумовую ХИП со случайным изменением периода повторения, длительности импульсов, а также линейной перестройкой частоты в пределах полосы пропускания ФОС. Многолучевое распространение учитывалось для первой АТТТП, число точек отражения от подстилающей поверхности принималось равным Np = 4. Для оставшихся 2 ПАП многолучевой канал распространения сигнала не моделировался.

Амплитудно-фазовые искажения (неидентичности) трактов обработки учитывались путем формирования предварительно искаженных комплексных импульсных характеристик ФОС, независимо в каждом из каналов обработки.

На рис. 3 приведены импульсные характеристики фильтра одиночного сигнала (ИХ ФОС) длиной Lf = 32 отсчета, используемые для предварительной частотной селекции в каналах обработки. Неидентичность фильтровых каналов моделировалась введением случайных отклонений (разброса) амплитуды дискретных

отсчетов ИХ ФОС относительно номинальных значений.

Угловое положение целей и ПАП в пределах сектора обзора представлено на рис. 4. По осям обозначено угловое отклонение относительно линии визирования РЛС - цель в радианах.

На рис. 5 представлены результаты моделирования работы устройства АКМИ при наличии трех источников АТТТП Слева приведено изображение индикатора угол - дальность при выключенном АКМИ, справа - при включенном АКМИ. По оси абсцисс отложен номер углового дискрета (0.127) в пределах сектора сканирования, по оси ординат - номер дискрета дальности (0.1630). Величина углового дискрета составляет 0.08 градуса, величина дискрета дальности - 25 метров, размер изображения в координатах дальность - угол составляет 42000 м х 10 град.

На рис. 6 приведены синфазная и квадратурная составляющие комплексной амплитуды сигналов в основном канале обработки на входе и выходе АКМИ.

Оценка интегрального коэффициент подавления составляет (7.11) дБ, что что свидетельствует о невысокой эффективности компенсации АШП.

Модель функционирования устройства адаптивной компенсации АШП на основе алгоритма последовательной регрессии

Алгоритм последовательной регрессии основан на использовании в рекуррентном алго-

Положение АШП относительно центра

Положение целей относительно центра

0.2

-0.2

0.2

«с 0

■ -----

-02 0 02 -02 0 02 Рис. 4. Угловое положение ПАП и целей в пределах сектора обзора наземной РЛС

-0 -1 Рис. 5. Экраны индикаторов дальность - угол при выключенном и включенном АКМИ для рекурсивного МНК -

алгоритма

ритме поиска оптимального ВВК (минимума многомерно оценки обра мех (КМП):

Классический метод получения оценки об-многомерной функции) методом Ньютона ратной КМП состоит в нахождении выбороч-оценки обратной ковариационной матрицы по- ной или рекурсивной оценки ковариационной

матрицы Л хх ^) с последующим ее обращением:

XV (() = Цо Л )У „ №, w)) = -ЦоК > )*« МО,

(8)

где Л-1 (0 - оценка обратной ковариационной матрицы помехи.

К

С) = -1 Кхх (0 - Кхх (0 ),

1 J-1 и Кхх «) = 7 X Xу (Г)ХН (Г), - ]=0

(9)

Рис. 6. Синфазная и квадратурная компоненты комплексной амплитуды сигнала в основном канале обработки в зависимости от дальности для фиксированного периода повторения 1пп = 65

где х ^) - ]-е выборочное значение вектора х(0; J - число выборочных значений вектора х(0, используемых для расчета оценки; Тф - постоянная времени сглаживания оценки Ы хх ) в рекурсивной фильтре первого порядка.

Оценка обратной КМП Ы-X ((), полученной согласно (9), требует повышенной производительности сигнальных процессоров, а также мер по снижению чувствительности оценки к степени обусловленности (отношению максимального и минимального собственных значений) матрицы, что влечет за собой увеличение разрядности и вычислительных затрат сигнального процессора. При плохо обусловленной КМП (отношения максимального и минимального собственных чисел составляет 100 и более) погрешность оценки возрастает, что приводит к снижению отношения сигнал - помеха на выходе АКМИ и уменьшению коэффициента подавления.

Одним из путей преодоления указанных трудностей является построение рекуррентного алгоритма оценки обратной ковариацион-

ной матрицы ЫхХх{1), исходя из уравнения:

Ы ^) =

(1f

V ТФ

J-1

^ Xу {I)хн (Г) - Ыхх «)

]=0

(10)

В дискретной форме уравнение примет вид:

Ы ж({п ) =

( (

Ы XX ^п-0 + а0

J-l Н ^

ы хх (^п-1)-х х у (()х н а) (11) ]=0

-1

А^

ап =-

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

ТФ

где А?д - интервал дискретизации входного сигнала (помехи).

Вместо переменной 1п в дальнейшем будем использовать подстановку ^ п .

Используя для оценки обратной КМП тождество Шермана-Моррисона-Вудбери (лемму об обращении матрицы) [10.13]:

= А-1 - А

(А + UCV) 1 =

-и (

С-1 + VA-1и) 1VA-1

(12)

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

Выполняя подстановку А = Ы хх (?), и = х(0, и = хн(0 и С = 1, получим:

Ы XX (п ) =

1

1 -а

Ы -хЧ-1) +

1»-х^М^ )хН (Гп )Ы-Ж-О

(13)

0

(1 -а0) + а0хн (гп)Ы-1 (^ №)

-1

В начальный момент времени можно принять

1 N-1

К-1(0) = ЦххI = "Г"1, =х хНх , (14)

]=0

где I - единичная матрица размером (пхп); схх -мощность помехи на входе; х] - вектор сигнальных отсчетов на входе соответствующих отводов адаптивного КИХ-фильтра; х у (?п) = h

н„

г ,■ -

] ]

выходной сигнал ФОС в ]-м канале обработки в момент времени п ,?((п), ]п) - сигнал на входе антенных систем основного и дополнительных каналов.

Рис. 7. Упрощенная структурная схема устройства АКМИ с использованием алгоритма последовательной регрессии для оценивания обратной КМП

Рис. 8. Экраны индикаторов дальность - угол при выключенном/включенном АКМИ для рекуррентного метода

оценки обратной КМП

На рис. 7 представлена упрощенная структурная схема устройства АКМИ с использованием алгоритма ПР для оценивания обратной КМП.

Проблема плохой обусловленности матрицы Л хх (п) в большинстве практически значимых случаев подавления АШП решается эври-

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

Рис. 9. Синфазная и квадратурная компоненты комплексной амплитуды сигнала в основном канале обработки в зависимости от дальности для фиксированного периода повторения 1пп = 65

Оценка интегрального коэффициент подавления составляет (12.19) дБ, что на (5.6) дБ превосходит эффективность работы АКМИ, реализованного на основе рекурсивного МНК -алгоритма.

Анализ результатов математического моделирования показывает:

- эффективность компенсации АТТТП в условиях многолучевого распространения сигналов, а также при неидентичных каналах обработки существенно зависит от используемого алгоритма оценки ВВК;

- для получения максимальной эффективности подавления АТТТП необходимо использовать методы, основанные на использовании обратной КМП, полученной непосредственным обращением исходной оценки матрицы помех либо путем последовательного приближения (регрессии);

- использовании для расчета ВВК рекуррентной оценки КМП по алгоритму последо-

вательной регрессии позволяет повысить эффективность работы АКМИ на (5.6) дБ, по сравнению с МНК - алгоритмом.

Заключение

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

Литература

1. Монзинго Р. А., Миллер Т. У. Адаптивные антенные решетки: Введение в теорию. Пер с англ. - М.: Радио и связь, 1986. -448 с., ил.

2. Адаптивные фильтры. Пер с англ. / Под ред. Коуэна К. Ф., Гранта П. М. - М.: Мир, 1988. - 392 с., ил.

3. Пистолькорс А. А., Литвинов О. С. Введение в теорию адаптивных антенн. - М.: Наука, 1991. - 200 с., ил.

4. Ширмам Я. Д., Манжос В. Н. Теория и техника обработки радиолокационной информации на фоне помех. - М.: Радио и связь, 1981. - 416 с., ил.

5. Адаптивная компенсация помех в каналах связи. / Под ред. Лосева В. И. - М.: Радио и связь, 1988. - 208 с.

6. Журавлев А. К., Лукошкин А. П., Поддубный С. С. Обработка сигналов в адаптивных антенных решетках. - Л.: Изд. ЛГУ, 1983. - 240 с., ил.

7. Сверхбольшие интегральные схемы и современная обработка сигналов: Пер. с англ. / Под ред. С Гуна, Х. Уайтхауса, Т. Кай-лата. - М.: Радио и связь, 1989. - 472 с., ил.

8. Ратынский М. В. Адаптация и сверхразрешение в антенных решетках. - М.: Радио и связь, 2003. - 200 с.

9. Дж. Голуб, Ч. Ван Лоун. Матричные вычисления. Пер с англ. - М.: Мир, 1999. - 548 с., ил.

10. Haykin S. Adaptive Filter Theory. Prentice - Hall, Inc., NJ, USA. 1996. - 992 p.

11. Woodbury matrix identity [Электронный ресурс]. - Режим доступа: https://en.wikipedia.org/wiki/Woodbury_ma-trix_identity. - Дата доступа: 05.12.2016.

12. Sherman-Morrison_formula [Электронный ресурс]. - Режим доступа: https://en.wikipedia.org/wiki/ Sherman-Morrison_formula. - Дата доступа: 05.12.2016.

13. Quasi-Newton method [Электронный ресурс]. - Режим доступа:: https://en.wikipedia.org/wiki/Quasi-Newton_ method. - Дата доступа: 05.12.2016.

References

1. Monzingo Р. А, Miller T. U. Adaptive antenna arrays: Introduction in the theory. Translate with English - М: Radio and communication, 1986. - 448 p.

2. The Adaptive filters. Translate with English / Under the editorship of Kouena K. F., the Grant of Item М. - М: the World, 1988. - 392 p.

3. Pistolkors A. A., Litvinov O. S. Introduction in the theory of adaptive antennas. - М: Science, 1991. - 200 p.

4. Shirman J. D., Manzhos V. N. The theory and technics of handling of the radar-tracking information against noises. -М: Radio and communication, 1981. - 416 p.

5. The Adaptive compensating of noises in communication paths. / Under the editorship of Losev V. I. - М: Radio and communication, 1988. - 208 p.

6. Zhuravlev A. K., Lukoshkin A. P., Poddubny S. S. Signal processing in the adaptive antenna arrays. - L: LSU, 1983. -240 p.

7. Very-large scale integrated circuits and the modern signal processing. Translate with English / Under the editorship of With Guna, H. Uajthausa, T. Kajlata. - М: Radio and communication, 1989. - 472 p.

8. Ratynsky M. V. Adaptation and super-resolution in antenna arrays. - М: Radio and communication, 2003. - 200 p.

9. J. Golub, C. Van Loun. Matrix calculations. Translate with English. - М: World, 1999. - 548 p.

10. Haykin S. Adaptive Filter Theory. Prentice Hall, Inc., NJ, USA. 1996. - 992 p.

11. Woodbury matrix identity. [The Electronic resource]. - The Access mode: https://en.wikipedia.org/wiki/Woodbury_ matrix_identity. - Access Date: 10.12.2016.

12. Sherman-Morrison formula [the Electronic resource]. - The Access mode: https://en.wikipedia.org/wiki/Sher-man-Morrison_formula. - Access Date: 03. 12.2016.

13. Quasi-Newton method [the Electronic resource]. - The Access mode: https://en.wikipedia.org/wiki/Quasi-Newton_ method. - Access Date: 23.11.2016.

Поступила После доработки Принята к печати

05.02.2017 22.02.2017 06.03.2017

Markevich V. E.

INCREASE OF EFFICIENCY OF SUPPRESSION OF THE ACTIVE NOISE JAMMING IN RADAR INTERCEPTION OF AIR TARGETS IN THE REVIEW

MODE

The modified algorithm of the adaptive compensating of the active noise jamming, based on usage of the recurrence estimation of inverse covariance matrix ofjamming in the course of an self-tuning is considered. The matrix estimation is formed by a method of a serial regression (SRM), a consequence application of a lemma about inversing of a matrix to differential - to the difference equation for a rounded estimation of direct covariance matrix of jamming (CMJ). The single estimation of straight line CMJ is formed as result dyad products of the instant vectors of noises in handling channels. In partial channels the adaptive transversal filters with the self-tuning on value of a signal of an aggregate error are implemented.

In article skeleton diagrams of devices of the compensating implemented by various methods are resulted. Mathematical modeling of operation of the algorithms, implemented by a method of the least squares (criterion of a minimum RMS error) and a method of a serial regression is performed.

Results of mathematical modeling of operation of the device of compensating of active noises jamming in the conditions of the active radio-electronic counteraction at multipath propagation of signals both non-identical amplitude and phase characteristics of channels of handling are resulted. Optimal values of technical parameters of devices of compensating of the active noises jamming providing obtaining of demanded values of coefficient of suppression are given. The comparative analysis of indexes of efficiency of classical and offered algorithms of compensating is made. Advantage of algorithm of compensating of the active noises jamming, using the recurrence estimation of a inverse covariance matrix ofjamming is shown.

The algorithm of functioning on the basis of SRM allows to lower essentially requirements to specialized VLSI (FPGA) at implementation is hardware - program complexes for preprocessing of the radar, radio communication and radio - navigation information in the conditions of influence of the active noises jamming in the presence of effects of multipath propagation of signals, in the presence of delay of signals on an antenna system aperture, and also at unremovable technological dispersion amplitude and phase characteristics of the microwave of channels of handling.

Keywords: the active noise jamming, the adaptive compensating of noises jamming, covariance matrix of jamming, a vector of weight coefficients, a method of a serial regression.

Маркевич Виталий Эдмундович - кандидат технических наук, ведущий научный сотрудник, ОАО «АЛЕВКУРП», e-mail: mark.vit@tut.by.

Markevich Vitaliy Edmundovich - Candidate of Engineering Sciences, leader scientist, OJSC «ALEVKURP». E-mail: mark.vit@tut.by

i Надоели баннеры? Вы всегда можете отключить рекламу.