Научная статья на тему 'Адаптивная подстройка весовых коэффициентов в антенной решетке'

Адаптивная подстройка весовых коэффициентов в антенной решетке Текст научной статьи по специальности «Компьютерные и информационные науки»

CC BY
528
123
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
АДАПТАЦИЯ / АНТЕННАЯ РЕШЕТКА / ЦИФРОВОЙ ФИЛЬТР / ВЕСОВЫЕ КОЭФФИЦИЕНТЫ / МАКРООПЕРАЦИЯ / АЛГОРИТМ / ADAPTATION / ARRAY / DIGITAL FILTER / WEIGHT FACTORS / MACRO / ALGORITHM

Аннотация научной статьи по компьютерным и информационным наукам, автор научной работы — Новиков Александр Валериевич

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

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

Похожие темы научных работ по компьютерным и информационным наукам , автор научной работы — Новиков Александр Валериевич

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

ADAPTIVE ADJUSTMENT OF WEIGHT FACTORS IN ARRAY

A conflict-free data structure for organizing a concurrent process of calculation of an adaptive array correlation matrix has been presented. A concurrent program of weight factor adjustment iterative process has been offered.

Текст научной работы на тему «Адаптивная подстройка весовых коэффициентов в антенной решетке»

УДК 621.396.6

АДАПТИВНАЯ ПОДСТРОЙКА ВЕСОВЫХ КОЭФФИЦИЕНТОВ

В АНТЕННОЙ РЕШЕТКЕ

А.В. Новиков

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

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

В условиях разнообразной и случайным образом изменяющейся внешней обстановки, обусловленной мешающими излучениями (активными помехами), для обеспечения требуемой эффективности систем связи, многофункциональных РЛС все более широкое применение получают адаптивные антенные решетки (ААР) [1]. Основным их элементом является цифровой адаптивный фильтр, предназначенный для формирования коэффициентов Wi,...,Wm, с которыми взвешиваются сигналы АР, предварительно преобразованные в цифровую форму в аналогово-цифровых приемных модулях.

Градиентный метод адаптации весовых коэффициентов в антенных решетках является широко известным в силу своей простоты и универсальности [1-3]. Рассмотрим ^-элементную антенную решетку, выходной сигнал которой определяется выражением:

y(t ) = WTx(t). (1)

Здесь WT - вектор комплексных весовых коэффициентов антенной решетки, x(t) - вектор-столбец входных сигналов x(t) = [x1(t), X2(t),... XN(t)]T. Обозначив через d(t) ожидаемый (опорный) сигнал, сигнал ошибки e(t) можно представить как:

e (t ) = d (t)-WTx (t). (2)

Для минимизации ошибки используется рекуррентный метод подстройки вектора весовых коэффициентов W (k), описываемый следующим выражением:

W (k +1)= W (k)- 2As (RxxW (k)- rTd ) (3)

где Yxd = [[x1(t)*d(t), X2(t)*d(t),..., XN(t)*d(t)]T ] - вектор коэффициентов корреляции входных сигналов решетки Xi(t) с опорным сигналом d(t), Rxx -корреляционная матрица входных сигналов антенной решетки:

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

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

Не

1

2 3

Х|

XI.L

X1.L+I xl,L+2 *1,L+3

X1,2L-I

XljiL+l xUL+2 Xl,nL+3

4,(n+l)L-l

*2 x2,2

x2,3

X2,L

x2,L+2 *2,L+3

X2,2L-I

x2,nL+l x2.nL+2 X2,nL+3

x2,(n+l)L-l

xN.I XN,2 *N,3

XN,L

XN.L+1 xN,L+2 XN.L+3

XN,2L-1

xN.nL< I xN.nL+2 XN,nL+3

xN,(n+])L-l

L+l L+2 L+3

2L

At;

Xl,

Xii

X1,L+I

X],L+2

Xi,2L-l

xI,nL+! X l,n L+2 X[,nLf3

Xl,(n+1)L-I

X2.l

x2,2 *2,3

X2.L

X2,L+! X2,L+2 Xi.L+3

X2.2L-I

x2,nL+l

X2,nL+2 X2,nL+3

X2,(n+I)L-I

Xn,! *N,2 *n,3

xn,l

XN.L+I

Xn,L+2

X-N.L+3

XN.2L-1

XN.JIL+I

XN,nL+2

XN,iiL+3

XM,(n+l)L-1

1 2.

К

K+l K.+ 2........... 2 К

NK+1 NK+2......(N+1)K

Рис. 1. Структура данных для вычисления корреляционной матрицы

Параллельная программа для вычисления корреляционной матрицы имеет следующий вид:

DCL X,Y(L VEC,T/L STREAM,N STREAM), c REG; for i=1 to N for j=1 to N; cadr init; c=0 ;

cadr corr;for k=1 to T/L do;c=c+pirsum (Y(*,k,i)*X(*,k,J); end; D(i,j)=c;end;end corr;end;

В соответствии с предложенной структурой данных из массивов X и Y выполняются циклы чтения данных и запись вычисленных коэффициентов корреляции в квадратный массив D. При структурно-процедурной реализации двойной цикл чтения будет реализован следующим образом. В цикле по i данные будут читаться из обоих массивов по адресам от 1 до A=T/L, где T - длительность радиолокационного сигнала, принимаемого антенной решеткой. Далее, в следующем цикле чтения в массиве X произойдет наращивание начального адреса на 1, а в массиве Y - возвращение к начальному адресу, и цикл чтения повторится снова. Процедура чтения повторяется до тех пор, пока не будет прочитан весь массив X. В цикле по J начальный адрес чтения массива Y последовательно наращивается на величину T/L. Запись вычисленных коэффициентов корреляции производится также двойным циклом. В цикле по i формируются i-t элементы строки, в цикле по j - формируется j-я строка массива. В данной задаче квадратная конфигурация массива D необходима для того, чтобы обеспечить максимальное распараллеливание вычислений, поскольку требования к быстродействию системы при решении задачи адаптивной настройки коэффициентов решетки являются чрезвычайно жесткими. Это связано с тем, что процесс адаптации занимает до нескольких сотен итераций, которые необходимо выполнять, чтобы итерационный процесс завершился.

Ниже приводится параллельная программа итерационного процесса подстройки весовых коэффициентов решетки:

DCL RXX(N VEC,N STREAM) ,W1,W2(N VEC,N stream) MEM, RXD(N STREAM) d(1, VEC) MEM; C,Dd(N VEC) COM; j=0

q: j=j+i;

cadr init_1;for k=1 to N by 1 do; c(k)=W1(k,k); end; end init_1;cadr adapt_odd;

for i=1 to N by1 do; for k=1 to N by1 do;

dd(k)=c(k)* RXX(k,i); W2(k,i) = (pirsum(dd(k)) - RXD(i))*delta -c(k)); end; end; end adapt_odd;

cadr init_2; for k=1 to N by 1 do; c(k)=W2(k,k) ; end; end init; cadr adapt_even ;

for i=1 to N by1 do; for k=1 to N by1 do; dd(k)=c(k)* RXX(k,i); W1(k,i) = (pirsum(dd(k)) - RXD(i))*delta -c(k)) ; end; end; end adapt_even;

cadr init_3;for k=1 to N by 1 do; c(k)=W2(k,k) ; end; end init_3; cadr convergence;

for k=1 to N; f= conv(c(k),W2(k,k); d(1)=pir_or(f); end; end convergence;

if d(1)=1 then go to q else end;

В кадре тй_1 в регистры макроопераций adapt_w1 из массива '1 записываются константы - веса антенной решетки. В кадре adapt_odd, реализующем выполнение итерации с четным номером, из массива ЯХХ читаются векторы элементов корреляционной матрицы, а из массива ЯХБ -элементы функции взаимной корреляции. Обработка информации выполняется на векторе макроопераций adapt_w1, а также макрооперации adapt_w2 и р^Биш. С выходов у2 макроопераций adapt_w1 результаты вычислений записываются в массив '2, при этом новые вычисленные веса будут расположены по диагонали данного массива. Аналогичным образом выполняется итерация с нечетным номером. На рис. 2 представлены граф-схемы макроопераций, используемых в процедуре вычисления весовых коэффициентов.

Рис 2. Граф-схемы макроопераций, используемых в процедуре вычисления весовых коэффициентов

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

Из результатов следует существенное преимущество параллельной реализации алгоритма рекуррентного вычисления вектора весовых коэффициентов при адаптации АР.

В заключение отметим, что при параллельной обработке необходимо увеличить на М 2 числа ячеек оперативной памяти для хранения результатов промежуточных вычислений.

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

1. Вендик О.Г., Парнес М.Д. Антенны с электрическим сканированием. Введение в теорию // под ред. Л. Д. Бахраха. 2001. 250 с.

2. Воскресенский Д.И. Антенны с обработкой сигнала. М.: Сайнс-Пресс, 2002. 80 с.

3. Курбатский С.А., Новиков А.В., Хомяков Д.А. Алгоритм вычисления вектора весовых коэффициентов в адаптивной антенной решетке //

Известия ТулГУ. Технические науки. Вып.11. Ч. 1. Тула: Изд-во ТулГУ, 2012. С. 261-268.

4. Акиншин Р.Н., Хомяков А.В., Новиков А.В. Оценки погрешностей при выборе математических моделей диаграмм направленности и рассеяния антенны // Журнал «Антенны». №12. 2013. С. 14-19.

Новиков Александр Валериевич, начальник сектора, rts@cdbae.ru, Россия, Тула, ОАО Центральное конструкторское бюро аппаратостроения

ADAPTIVE ADJUSTMENT OF WEIGHT FACTORS IN ARRAY

A.V. Novikov

A conflict-free data structure for organizing a concurrent process of calculation of an adaptive array correlation matrix has been presented. A concurrent program of weight factor adjustment iterative process has been offered.

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

Key words: adaptation, array, digital filter, weight factors, macro, algorithm.

Novikov Alexandr Valeriyevich, head of sector, rts@,cdbae.ru, Russia, Tula, JSC Central Design Bureau of Apparatus Engineering

УДК 621.391

ПСЕВДОМЕДИАННАЯ ФИЛЬТРАЦИЯ ДЛЯ ОБРАБОТКИ ОПТИЧЕСКИХ ИЗОБРАЖЕНИЙ

Н.С. Акиншин, А.В. Андреев, А.И. Кальной

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

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

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

у(ш, п) = F{x(m - и п - j)}, (1)

где i = 0, ..., p, j = 0, ..., p, k = 1, ..., Q, l = 1, ..., Q, m = 1, ..., n п = 1, ..., N.

Для подавления спекл-шумов на оптических изображениях в системах с синтезированием апертуры антенны в качестве операции Б использу-

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