Устройство формирования весовых коэффициентов в условиях нестационарной помеховой обстановки

 

Использование: обеспечение сопровождения движущихся объектов в условиях быстрых изменений помеховой обстановки. Сущность: содержит генератор опорного сигнала (1), блок (2) вычитания, блок (3) анализа помеховой обстановки, блок (4) формирования матричных коэффициентов передачи калмановского фильтра, блок (5) формирования скалярного произведения, блок (6) суммирования, блок задержки (7), блок (8) формирования оценки экстраполяции , что позволяет повысить устойчивость адаптивного формирования весовых коэффициентов , 1 з.п. ф-лы, 4 ил,

СОЮЗ СОВЕТСКИХ

СОЦИАЛИСТИЧЕСКИХ

РЕСПУБЛИК (51)5 6 01 S 7/36

ГОСУДАРСТВЕННОЕ ПАТЕНТНОЕ

ВЕДОМСТВО СССР (ГОСПАТЕНТ СССР) ОПИСАНИЕ ИЗОБРЕТЕНИ

К АВТОРСКОМУ СВИДЕТЕЛЬСТВУ (21) 4902005/09 (22) 11.01.91 (46).15.04.93. Бюл, N 14 (72) А.А,Адаменко, B.À.Äoðoùóê, В.И.Шандрык и Ю,M,Ãâoýäü (5@ Монзинго P.À., Миллер Т.У. Адаптивные антенные решетки, Введение в теорию. M.:

Радио и связь, 1986, стр. 169 — 73. (54) УСТРОЙСТВО ФОРМИРОВАНИЯ ВЕСОВЫХ КОЭФФИЦИЕНТОВ B УСЛОВИЯХ

НЕСТАЦИОНАРНОЙ ПОМЕХОВОЙ ОБСТАНОВКИ

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

На фиг, 1 приведена структурная схема устройства формирования весовых коэффициентов в условиях нестационарной помеховой обстановки; на фиг. 2 — структурная схема блока анализа помеховой обстановки, Описываемое устройство формирования весовых коэффициентов содержит последовательно соединенные генератор

„„ Ы„„1809401 А1

2 (57) Использование: обеспечение сопровождения движущихся объектов в условиях быстрых изменений помеховой обстановки, Сущность . содержит генератор опорного сигнала (1), блок (2) вычитания, блок (3) анализа помеховой обстановки, блок(4) формирования матричных коэффициентов передачи калмановского фильтра, блок (5) формирования скалярного произведения, блок (6) суммирования, блок задержки (7), блок (8) формирования оценки экстраполяции, что позволяет повысить устойчивость адаптивного формирования весовых коэффицйе тов. 1 э.п. ф-лы, 4 ил. опорного сигнала 1, первый блок вычитания

2, блок формирования матричных козффи- Ф циентов передачи калмановского фильтра 4, блок суммирования 6, блок задержки 7, блок формирования оценки экстраполяции 8, а блок формирования скалярного проиэведе- "р ния 5, выход которого соединен со вторым входом первого блока вычитания 2, а выход блока формирования оценок экстраполяции

8 подключен также ко второму входу блока + суммирования 6, выход первого блока вычи- С) тания 2 подключен к первому входу 11 блока анализа помеховой обстановки 3, а его выход 13 ко второму входу блока формирования матричных. коэффициентов передачи калмановского фильтра 4, второй выход которого соединен со вторым входом 12 блока анализа помеховой обстановки 3, третий вход блока формирования матричных коэффициентов передачи калмановского фильтра 4 и второй вход блока формирования скалярного произведения 5 соединены между собой и являются входом устройства 9, 1809401 выходом устройства является выход 10 блока задержки 7.

Блок анализа помеховой обстановки 3 содержит последовательно соединенные квадратор 14, блок деления 15, M-тактовый блок задержки 16, второй блок вычитания

18, пороговое устройство 19, многораэрядHbIA ключ 20, матричный сумматор 21, выход

13 которого подключен ко второму входу блока формирования матричных коэффициентов передачи калмановского фильтра 4, а также накапливающий сумматор 17, вход которого подключен к выходу блока деления

15, а выход ко второму входу второго блока вычитания 18, постоянное запоминающее устройство 22, первый выход которого подключен ко второму входу порогового устройства 19, второй выход ко второму входу многоразрядного ключа 20, третий выход ко второму входу матричного сумматора 21, вход 11 квадратора подключен к выходу. первого блока вычитания 2, второй вход 12 блока деления подключен ко второму выходу блока формирования матричных коэффициентов передачи калмановского фильтра 4, Рассмотрим функционирование устройства формирования весовых коэффициентов в условиях нестационарных помех.

Уравнение состояния, отражающее динамические характеристики оптимального вектора весовых коэффициентов айтенной решетки

К...(+ 1) = ©(k+ 1/ ) don (k)+ 4(k) +

О(в) д (к + 1, m), 1Иопт(0) = Яо, где Ф (k+1/k) — переходная матрица состояния, ((К) — вектор белого шума с нулевым средним и ковариационной матрицей, Ej ф(k) g (j)) = Q(k) д (К!), д (К)) — символ

Кронекера, включает также случайный вектор 0(m), определяющий характер измерения обстановки, Уравнение наблюдения определим установив соответствие между отсчетами опорного сигнала d(k) и значением выходного сигнала антенной решетки x (k)Won (k)

d(k) = xò (k)Wonò(k) + v(k), где 1 I — знак транспонирования, v(k)— белый гауссовский шум с нулевым средним и дисперсией E(v(k)v(j)} = п2 д (k,j) характеризует погрешности аппроксимации опорного сигнала, x(k) — вектор отсчетов xI(k) сигнала с выходов антенных элементов решетки.

Оптимальная оценка вектора весовых коэффициентов решетки по критерию минимума среднеквадратической ошибки в соответствии с результатами теории калмановской фильтрации при H(m) = 0 определяется соотношением;

Д

WonT(k/k)= Wont(k/k-1) + K(k)z(k), (3) 10 в котором и h

Wom(k/k-1) Ф(К/k-1)Wom(k-1/k-1)-(4) 15 оценка экстраполяции, A

z(k) = d(k) - х (k)W(k/k-1)" (5) 20 обновляющий процесс, K(k) = P(k/k-1)x(k)/cJz (k)(6) матричный.коэффициент усиления, Входящая в выражение (6) ковариационная матрица оценки экстраполяции

P(k/k-1) задается соотношением

P(k/k-1) = Ф(М/k-1)P(k-1/k-1) Ф (1с/k-1)+Q(k) зо р

° ковариационная матрица оценки весового вектора(Jz (k) = x (k)P(k/k-1)x(k) + oz (k) (9)

Процедура адаптации; описываемая выражениями (3)-(9) при стационарной помеховой обстановке О (m) = Отобеспечивает минимальные флюктуации оценки весового

45 вектора, однако обладает. плохой устойчивостью в случае несоответствии реальной и принятой модельной зависимости Won от времени. Последнее обусловлено в первую очередь скачкообразными изменениями по50 меховой обстановки. Обнаружение "скачка" оптимального вектора весовых коэффициентов s устройстве осуществляется на основе анализа статистических характеристик обновляющего процесса z(k). В стационар55 ном режиме обновляющий процесс z(k) в согласованном с моделью фильтре представляет собой белый гауссовский шум с нулевым средним и дисперсией, определяемой (9):

P(k/k) = P(k/k-1) - P(k/k-1) х (к)Р(1/k-1) х

35 хо (), (8) P(o/о) = Po дисперсия обновляющего процесса, 1809401 (12) Q(k) = Q(k) + y, нулевым средним и дисперсией, определяемой (9).

Если вектор Wppr(k) претерпевает "скачок, то статистические характеристики процесса z(k) изменяются. Наиболее простой метод обнаружения "скачка" состоит в формировании статистики вида

l(k) =, z (k) о,2(k), (10)

j =к — @+i имеющий распределение типау- квадрат с M степенями свободы. При этом решающее правило приобретает вид

< (< ) — о — стационарная поме <1

g(f<) = ховая обстановка

t(k) >1, — "скачок", . (.1 1) где Io — значение порога, Качественные характеристики такого правила обнаружения рассчитываются известным способом, используя таблицы g квадрат распределения.

Для предотвращения:потери устойчивостй и ускорения переднего процесса после обнаружения "скачка" необходимо увеличить чувствительность устройства к вновь поступающей информации и расширить его полосу пропускания. . Для этого преднамеренно увеличиваются диагональные элементы ковариационной матрицы возмущающего процесса где у» 1,! — единичная матрица.

Исходя из (3)-(12) устройство функцио-. нирует следующим образом. На первый вход первого блока вычитания 2 поступает отсчет опорного сигнала д(к), на вычитающий вход в виде параллельного двоичного

" кода поступает, оценка выходного сигнала антенной решетки, сформированного в результате умножения оценки вектора весоЛ вых коэффициентов%(!<-1/k-1) с выхода. блока задержки 7 на матрицу Ф (k/k-1) p.. блоке формирования оценки экстраполяции

8 и на вектор отсчетов выходных сигналов антенных элементов x(k) в блоке формирования скалярного произведения 5, Отсчет обновляющего процесса z(k),(5) поступает на . входы блока анализа помеховой обстановки

3 и блока формирования-матричных коэффи-. циентов передачи калмэновского фильтра 4.

В .блоке формирования матричных коэффициентов передачи калмановского фильтра 4 осуществляется формирование матричного коэффициента усиления K(k) в соответствии . с (6) — (9), расчетного значения дисперсии обновляющего процесса о . (k)(9), и произведег ния компонентов вектора K(k) на значения отсчета z(k). В блоке суммирования 6 суммируются значения оценки экстраполяции

5 Q(k/k 1)Wpqy(k 1/k 1) с выхода блока форми рования оценки экстраполяции 8 in произве- дения Цс)л(К) и формируется оптимальный вектор весовых коэффициентов W«T(k/k),,который через время Т поступает на выход

10 10 блока задержки 7.

Оценка помеховой обстановки (обнаружение "скачка" вектора весовых коэффици,ентов) и формирование корреляционной матрицы возмущающего процесса Q осуще15 ствляется в блоке анализа помеховой обстановки 3.

Отсчет обновляющего процесса z(k) поступает на вход 11 квадратора 14 и затем на вход делимого блока деления 15, На вход

20 делителя 12 блока деления 15 поступает значение дисперсии обновляющего процесса о, сформированное в блоке формирова2 ния матричных коэффициентов передачи калмановского фильтра. С выхода блока де25 ления 15 значение частного z (k)! az (k) поступает на накапливающий сумматор 17 и затем на второй вход второго блока вычитания 18. Это же значение частного через

М-тактовый блок задержки 16 с. запаздыва30 нием на М тактов поступает на вычитающий вход второго блока вычитания 18, на выходе которого формируется значение достаточной статистики т(к) =,> z (k) о (k), посту2 -2

35 пающее на первый вход порогового устройства 19, На второй вход порогового устройства 19 поступает значение порога (о из постоянного запоминающего устройства

40 (ПЗУ) 22 . В том случае, если значение достаточной статистики l(k) больше значения порога 1о, на восходе порогового устройства 19 формируется разрешающий потенциал, по.ступающий на первый вход многоразрядно45 го ключа 20 и на.его выход передаются значения элементов матрицы jl из ПЗУ22 со второго выхода, поступающее на второй вход многоразрядного ключа 20. В матричном сумматоре 21 производится поэлемен50 тное суммирование матрицы jl и матрицы

0(k), поступающей на второй вход матричного сумматора,21 с третьего выхода ПЗУ

22. В случае, если значение достаточной статистики Як) меньше значения порога

55 1, запрещающий потенциал с порогового устройства 19 запирает многоразрядный ключ 20 и элементы матрицы, поступающие на первый вход матричного сумматора 21, приобретают нулевые значения, и значения элементов матрицы Q(k) на выходе матрич1809401 ного сумматора 21 равны значениям элементов матрицы 0(k).

Блок формирования матричных коэффициентов передачи калмановского фильтра 4 представляет собой специализированный вычислительный блок, осуществляющий формирование матричнОго коэффициента усиления K(k) в соответствии с (6) — (9), вычисление значения дисперсии обновляющего процесса о (k) в соответствии с (9), а также произведения компонентов вектора K(k) на значения отсчета z(k), На фиг. 3 представлена структурная схема блока формирования оценки экстраполяции 8. Он осуществляет перемножение оценки вектора весовых коэффициентов

W(k-1/k-1) с выхода блока задержки 7 на матрицу ©(k/k-1) и представляет собой N блоков формирования скалярного произведения 23,. на первые входы которых подаются значения элементов матрицы Ф(к/k-1), э на вторые входы отсчеты вектора весовых коэффициентов W(k-1/k-1).

На фиг. 4 представлена структурная схема блока формирования скалярного произведения 5. Он осуществляет перемножение.значений весовых коэффициент тов W(k/k-1), полученных на выходе лока оценки экстраполяции 8, на вектор отсчетов выходных сигналов антенных элементов x(k) в соответствии с (5), Блок представляет со- бой N перемножителей 24., на первые входы которых подаются отсчеты выходных сигналов антенных элементов x(k), на вторые входы значения весовых коэффициентов

W(k/k-1), а также сумматор 25, где суммируются сигналы с выходов перемножителей.

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

Формула изобретения

1. Устройство формирования весовых коэффициентов в условиях нестационарной помеховой обстановки, содержащее после45 ход — к второму входу матричного сумматора. при этом вход квадрэтора является

50 ховой обстановки.

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

БФМКПКФ и второй вход, блока формирования скалярного произведения соединены между. собой и являются входом устройства, выход линии задержки является выходом устройства.

2, Устройство формирования весовых коэффициентов по п.1, о т л и ч а ю щ е е с я тем, что блок анализа помеховой обстановки содержит последовательно соединенные квадратор, блок деления, М-тактовый блок задержки, блок вычитания, пороговый блок, многоразрядный ключ, матричный сумматор; выход которого является выходом 6лока анализа помеховой обстановки, а также накапливающий сумматор, вход которого подключен к выходу блока деления; а выход — к второму входу блока вычитания, постоянный запоминающий блок. первый выход которого подключен к второму входу порогового блока, второй выход — к второму входу многоразрядного ключа, третий Bblпервым входом блока анализа помеховой обстановки,,второй вход блока деления является вторым входом блока анализа поме1809401

1809401 ж срюе Мвмуф Вумай

Аж@/

Фиа 4

Составитель А. Адаменко

Техред М.Моргентал

Корректор И. Шмакова

Редактор B. Трубченко

Производственно-издательский комбинат "Патент", г. Ужгород, ул.Гагарина, 101

Заказ 12&4 Тираж Подписное

ВНИИПИ Государственного комитета по изобретениям и открытиям при ГКНТ СССР

113035, Москва, Ж-35, Раушская наб., 4/5