Адаптивный цифровой фильтр

 

Изобретение относится к области вычислительной техники и может быть использовано для фильтрации цифровых сигнаTffjrm 5л. лов в различных системах обработки информации . Цель изобретения - увеличение точности за счет подавления коррелированных помех улучшением устойчивости процедуры вычисления весового вектора путем отыскания квазирешения системы линейных уравнений на компактном множестве значений весового вектора. Поставленная цель достигается за счет того, что в состав устройства входят ключи 1, 2, коррелятор 3, блок 4 определения весовых коэффициентов , элементы 5 задержки, умножители 6, сумматор 7, синхронизатор 8. Блок определения весовых коэффициентов содержит узел 9вычисления собственных значений, узел 10вычисления вектора, делители 11, узлы 12 селекции, узел 13 перехода. 9 ил. & (Л 00 оо со

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

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

РЕСПУБЛИН

А1

ОПИСАНИЕ ИЗОБРЕТЕНИ

К А ВТОРСКОМУ СВИДЕТЕЛЬСТВУ

ГОСУДАРСТВЕННЫЙ НОМИТЕТ СССР

ПО ДЕЛАМ ИЗОБРЕТЕНИЙ И ОТНРЫТИЙ (21) 4138757/24-24 (22) 27.10.86 (46) 07.04.88. Бюл. № 13 (71) Московский авиационный институт им. Серго Орджоникидзе (72) В. Я. Плекин и М. М. Леднев (53) 681.32 (088.8) (56) Авторское свидетельство СССР № 800995, кл. G 06 F 15/353, 1981.

Гольденберг Л. М., Матюшина Б. Д., Поляк М. Н. Цифровая обработка сигналов. — М.: Радио и связь, 1985, с. 165, рис. 3.4 (54) АДАПТИВНЫЙ ЦИФРОВОЙ

ФИЛЬТР (57) Изобретение относится к области вычислительной техники и может быть использовано для фильтрации цифровых сигнаTwm &.

„„SU„„1387173 (5И 4 H 03 Н 17 06 G 06 F 15/353 лов в различных системах обработки информации. Цель изобретения — увеличение точности за счет подавления коррелированных помех улучшением устойчивости процедуры вычисления весового вектора путем отыскания квазирешения системы линейных уравнений на компактном множестве значений весового вектора. Поставленная цель достигается за счет того, что в состав устройства входят ключи 1, 2, коррелятор 3, блок 4 определения весовых коэффициентов, элементы 5 задержки, умножители 6, сумматор 7, синхронизатор 8. Блок определения весовых коэффициентов содержит узел

9 вычисления собственных значений, узел

10 вычисления вектора, делители 11, узлы 12 селекции, узел 13 перехода. 9 ил.

1387173

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

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

На фиг. 1 изображена структурная схема адаптивного цифрового фильтра; на фиг. 2 — структурная схема узла вычисления вектора; на фиг. 3 — структурная схема весового узла; на фиг. 4 — структурная схема узла вычисления собственных значений; на фиг. 5 — структурная схема узла селекции; на фиг. 6 — структурная схема узла перехода; на фиг. 7 — схема коррелятора; на фиг. 8 — схема накопителя; на фиг. 9 — управляющие сигналы на выходах синхронизатора.

Предлагаемый фильтр содержит ключи 1 и 2, коррелятор 3, блок 4 определения весовых коэффициентов, элементы 5 задержки, умножители 6, сумматор 7, синхронизатор 8.

Блок 4 определения весовых коэффициентов содержит узел 9 вычисления собственных значений, узел 10 вычисления вектора, делители 11, узлы 12 селекции, узел 13 перехода.

Узел 9 вычисления вектора (фиг. 2) содержит умножитель 14, весовой узел 15, делители 16 и 17, сумматор 18, вычитатель

19, сумматор 20. Весовой 15 узел (фиг. 3) содержит умножитель 21 и регистр 22 числа.

Узел 9 вычисления собственных значений (фиг. 4) содержит квадратор 23, весовой узел 24, сумматор 25, элемент 26 вычислений -/Х (квадратного корня), квадратор 27, вычитатель 28, весовой узел 29, сумматоры 30 и 31, весовой узел 32, сумматор 33, регистр 34 числа, вычитатель 35.

Узел 12 (фиг. 5) содержит элемент 36 вычисления модуля, компаратор 37, регистр

38 (хранения верхней границы), элементы

И 39 и 40, элемент НЕ 41, элемент ИЛИ 42, регистр 43 (хранения нижней границы), элементы И 44 и 45, элемент НЕ 46, компаратор 47 и элемент ИЛИ 48.

Узел 13 перехода (фиг. 6) содержит группу элементов НЕ 49, весовой узел 50, делители 51 и 52, умножители 53 и 54, сумматоры 55 — 57.

Коррелятор 3 (фиг. 7) содержит регистр

58 сдвига, первый 59, второй 60 и третий 61 умножители, квадратор 62, накопители 63 — 66 и три делителя 67 — 69.

Каждый из накопителей (фиг. 8) содержит регистр 70 сдвига, сумматор 71, счетчик 72, компаратор 73, регистр 74 числа и элемент И 75.

Цифровой фильтр работает следующим образом.

Адаптивный режекторный фильтр, синтезированный в соответствии с адаптивным байесовским подходом по критерию максимума коэффициента подавления помехи, использует для вычисления весового вектора Н = (hl, h,..., hl ) следующий алгоритм:, и

RH= — В; Н= — R В. и где R= (рк е1 — оценка коррелнцнонной матрицы помехи;

B= р, — оценка корреляционного вектора; оценка коэффициентов корреляции, 15 1, j, 1 = 1,N (N — порядок фильтра).

Так как процесс измерения параметров помехи сопровождается ошибками (т.е. A =

= К + оК, В = В + 6В, где R,  — истинные значения корреляционных характеристик, 6Р, 6 — матрица и вектор оши2О бок), то вектор Н, полученный из решения системы уравнений (1), отличается от своего оптимального значения Но = — R В:

Н= Н +6Н.

При этом относительное по норме возмущение весового вектора Л Н = II 6H 11/ ll H II можно определить при помощи следующего соотношения:

ЛН = cond R (ЛК+ЛВ) /(1 — ЛКсоп(1К), (3) где ЛК= 11оК!(/IIRII, AB= II6BII/II BII — относительные возмущения матрицы и вектора;

cond R = Aмаке./р.мнн — число обусловленности матрицы; амин, рма«с — минимальное и MBKcHMBJlbное собственные значения матрицы К;

li ° II — означает норму матрицы или век35 тора.

Из анализа (3) следует, что в случае сильно коррелированных помех, характеризующихся плохо обусловленной корреляционной матрицей (cond R — oo), вычисленное значение вектора Н может значительно отли4О чаться от своего оптимального значения Н11, что приводит к существенному снижению эффективности подавления помех.

Получение устойчивых решений системы

4 уравнений (1) возможно лишь при условии, когда значения вектора Н образуют компактное множество. В этом случае применение метода квазирешений плохо обусловленных систем линейных уравнений приводит к следующему алгоритму. Пусть все значения весового векторра Н принадлежат компактному множеству F: HgF. Тогда в качестве решения системы уравнений (1) берется такой вектор Н*gF, что д л дк

IIRH* + Bll = ш1 II КН + В Н. (4

55 Рассмотрим расширенный весовой вектор

Й = (h11, Н), который определяется из решения расширенной системы уравнений

RH= 0, () 73 р/Xh;= О. (б) 10

30 — 1+а; 50

13871 и. где К = (рд; j — расширенная корреляционная матрица, k, m = О, N.

Из формулы (5) при 1зо = 1 следует система уравнений (1). Допустим, что р,„ „„— вЂ: 1 (сильно коррелированная помеха). Тогда решением системы (5) может служить такой вектор Й = (hp,...,hN), составляющие которого удовлетворяют условию

Соотношению (6) отвечает вектор с биномиальными коэффициентами:

Ф

;hp = 1, h;= (— 1) CN, i = 1, N, где С вЂ” число сочетаний из N no i. 15

Отсюда следует, что каждая координата вектора Н не превосходит некоторой константы

h; и Рг. Ис- 25 ходя из этого, для каждой координаты вектора определяются нижняя и верхняя гра1.

6 ницы Ьг и Ьг с учетом полученных границ и соотношений (4) и (7) определяется алгоритм для вычисления весового вектора:

h;= h;,еслиЬгаh;аЬ,;

h; = Ьг если h; О Ьг (8)

h;= Ьг,если Ь;(Ьг.

Алгоритм (8) обладает свойством устойчивости, так как решение плохо обусловленной системы линейных уравнений ищется на компактном множестве F.

С использованием алгоритма (8) синтезирован нерекурсивный режекторный фильтр третьего порядка.

Для избежания трудоемкой операции 40 обращения матрицы используют ортогональное преобразование (характеризующееся матрицей С) такое, что

С RC=Л, где Л= (Л"...Л»- ) — диагональная матри- 45 ца из собственных значений.

Собственные числа корреляционной матрицы R определяются соотношениями

Ai= 1+

Лг= 1 — рг, г /г ()0) з = 1 -+- рг+(рг+8р ) 1+ в, с учетом которых матрица С принимает вид

Преобразование (9) переводит систему уравнений (1) в базис из собственных векторов матрицы:

RH = — В = ) ЛН = — D, (12) где D =С В вЂ” корреляционный вектор помехи в базисе из собственных векторов, составляющие которого без учета знака в соотношении (12) определяются следующим образом:

di = pi + 2pip2/а + рз, pl + p3 (13) с1з = p> + 2р рг/в + рз.

Для перехода в старый базис используется преобразование Н = СН, которое в координатной форме принимает вид = hi юг+1 3

h = 2р h>/à + 2р hp/B; (14 )

ha = Ь + "г + "з.

С выхода синхронизатора 8 на управляющий вход ключа 1 начинают поступать импульсы То оценки (фиг. 9а), в результате чего цифровые отсчеты входного сигнала поступают в коррелятор 3 (фиг. 7), который производит оценку коэффициентов корреляции помехи в соответствии с алгоритмом м м р; = ZU (п) U (и — i)/XU (и), (15) где U (n) — отсчеты входного сигнала;

М вЂ” объем выборки оценивания.

С выходов умножителей 59 — 61 и квадратора 62 отсчеты сигнала поступают в накопители 63 — 66 (фиг. 8). Счетчик 72 производит подсчет числа импульсов То, при этом его выходной код поступает в компаратор 73, на другой вход которого с выхода регистра 74 числа подается код числа М. В момент равенства кодов на выходе компаратора 73 формируется единичный сигнал (уровень) и через элемент И 75 накопленные значения сигналов поступают в делители 67 — 69, на выходах которых в соответствии с (15) формируются оценки коэффициентов корреляции.

Полученные таким образом значения

A р;, i = 1, 3, поступают в узел 9 вычисления. собственных значений матрицы R, который в соответствии с алгоритмом (10) определяет величины а, в, Ai, Лг и Лз. Узел 10 вычисления вектора D по оценкам р;, i = 1, 3, и значениям а и в, которые поступают на его входы с выходов узла 9, определяет в соответствии с (13) компоненты вектора D .

После чего на выходах делителей 11 формируются составляющие весового вектора в базисе из собственных векторов

1 — 1 1

С = 2p i/а О 2р./1г (11)

1 1 1

h; — d;/A;, i — 1,3, поступающие в узлы 12, реализующие операцию (8) .

1387173

Через элемент 36, вычисления модуля абсолютные значения Ц поступают в компаратор 37, в котором сравниваются со своими верхними границами hq, При выполнении

I ( условия ф; ) Ь2 на выходе компаратора формируется единичный уровень, благодаря чему через элемент И 39 значение hB2 поступает в компаратор 47. Если выполняется усI t ловие Ц (Ь2, то на выходе компаратора 37 формируется нулевой уровень, поступающий через элемент НЕ 41 на вход элемента И 40. В результате этого значение

ih,.i через элемент ИЛИ 42 подается в компаратор 47. Верхняя (фиг. 5) часть схемы работает аналогичным образом.

С выходов логических узлов значения

/ Ц поступают в узел 13 перехода, который в соответствии с алгоритмом (14) осуществляет переход в старый базис. Причем группа элементов НЕ 49 необходима для учета знака минус в соотношениях (12), а сумматор 57 работает как вычитающее устройство.

С выходов блока 4 определения весовых коэффициентов отсчеты весового вектора

h;, i = 1, 3, поступают на входы умножителей 6. После чего с выхода синхронизатора 8 на управляющий вход ключа 2 и тактовые входы элементов 5 задержки начинают поступать импульсы Тоб обработки (фиг. 9б) . В результате этого отсчеты входного сигнала проходят в настроенный нерекурсивный фильтр, при этом на выходе сумматора 7 формируется выходной сигнал в виде взвешенной суммы задержанных отсчетов.

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

Адаптивный цифровой фильтр, содержащий блок определения весовых коэффициентов, три элемента задержки, три умножителя и сумматор, выход которого является информационным выходом фильтра, выход первого элемента задержки подключен к первому входу первого умножителя и информационному входу второго элемента задержки, выход которого подключен к первому входу второго умножителя и информационному входу третьего элемента задержки, выход которого подключен к первом у входу третьего умножителя, i-й (i = 1; 3) выход блока определения весовых коэффициентов подключен к второму входу i-го умножителя, выход которого подключен к i-му входу сумматора, отличающийся тем, что, с целью повышения точности, в него введены два ключа, коррелятор и синхронизатор, первый выход которого подключен к тактовым входам первого, второго и третьего элементов задержки и управляющему входу первого ключа, выход которого подключен к информационному входу первого элемента задержки и четвертому входу сумматора, второй выход синхронизатора подключен к тактовому входу коррелятора и управляющему входу второго ключа, выход которого под5

55 ключен к информационному входу коррелятора, 1-й выход которого подключен к 1-му входу блока определения весовых коэффициентов, информационные входы первого и второго ключей соединены и являются информационным входом фильтра, причем блок определения весовых коэффициентов содержит три делителя, три узла селекции, узел перехода, узел вычисления собственных значений и узел вычисления вектора, i-й выход которого подключен к первому входу i-го делителя, выход которого подключен к входу

i-го узла селекции, выход которого подключен к i-му входу узла перехода, второй вход i-ro делителя подключен к i-му выходу узла вычисления собственных значений, первый и второй входы которых соединены соответственно с первым и вторым входами узла вычисления вектора и являются соответственно первым и вторым входами блока определения весовых коэффициентов, третьим входом которого является третий вход узла вычисления вектора, четвертый и пятый выходы узла вычисления собственных значений подключены соответственно к четвертому и пятому входам узла вычисления вектора, и узла перехода, шестой вход которого подключен к первому входу блока определения весовых коэффициентов, i-м выходом которого является i-й выход узла перехода, причем узел вычисления вектора содержит два делителя, два сумматора, вычитатель, два умножителя и регистр числа, выход которого подключен к первому входу первого умножителя, выход которого подключен к первым входам первого и второго делителей, выходы которых подключены к первым входам соответственно первого и второго сумматоров, вторые входы которых соединены с первым входом вычитателя, первым входом второго умножителя и являются первым входом узла вычисления вектора, вторым входом которого является второй вход второго умножителя, выход которого подключен к второму входу первого умножителя, третьи входы первого и второго сумматоров соединены с вторым входом вычитателя и являются третьим входом узла вычисления вектора, четвертым и пятым входами которого являются вторые входы соответственно первого и второго делителей, выход первого сумматора, выход вычитателя и выход второго сумматора являются соответственно первым, вторым и третьим выходами узла вычисления вектора, при этом узел вычисления собственных значений содержит четыре сумматора, два вычитателя, элемент вычисления - /Г, четыре регистра числа, три умножителя и два квадратора, вы ходы которых подкл ючен ы соответственно к первому входу первого умножителя и первому входу первого сумматора, выходы которых подключены соответственно к второму входу первого сумматора и входу элемента вычисления - l выход которого под1387173 ключен к первым входам первого вычислителя и второго сумматора, выходы которых подключены к первым входам соответственно второго и третьего умножителей, выходы которых подключены к первым входам соответственно третьего и четвертого сумматоров, вторые входы которых соединены с первым входом второго вычитателя и подключены к выходу первого регистра числа, вторые входы первого, второго и третьего умножителей подключены к выходам соответственно второго, третьего и четвертого регистров числа, вход первого квадратора является первым входом узла вычисления собственных значений, вторым входом которого являются соединенные между собой вход второго квадратора, вторые входы первого и второго вычитателей и второй вход второго сумматора, выходы второго и третьего умножителей являются соответственно четвертым и пятым выходами узла вычисления собственных значений, первым, вторым и третьим выходами которого являются выходы соответственно третьего сумматора, второго вычитателя и четвертого сумматора, причем каждый узел селекции содержит два регистра, четыре элемента И, два элемента

НЕ, два компаратора, два элемента ИЛИ и элемент вычисления модуля, выход которого подключен к первому входу первого элемента И и первому входу первого компаратора, выход которого подключен к первому входу второго элемента И и входу первого элемента НЕ, выход которого подключен к второму входу первого элемента И, выход которого подключен к первому входу первого элемента ИЛИ, выход которого подключен к первым входам третьего элемента И и второго компаратора, выход которого подключен к второму входу третьего элемента И и входу второго элемента НЕ, выход которого подключен к первому входу четвертого элемента И, выход которого подключен к первому входу второго элемента ИЛИ, второй вход которого подключен к выходу третьего элемен5 та И, выход первого регистра подключен к вторым входам первого компаратора и второго элемента И, выход которого подключен к второму входу первого элемента ИЛИ, выход второго регистра подключен к вторым входам второго компаратора и четвертого элемента И, вход элемента вычисления модуля является входом узла селекции, выходом которого является выход второго элемента ИЛИ, при этом узел перехода содержит три сумматора, два делителя, три умножителя, регистр числа и группу элементов НЕ, выход первого элемента НЕ группы подключен к первым входам первого сумматора и первого умножителя, выход которого подключен к первому входу второго сумматора, второй вход

20 которого подключен к выходу второго умножителя, первый вход которого соединен с вторым входом первого сумматора и подключен к выходу второго элемента НЕ группы, выход второго элемента которой подключен к третьему входу первого суммато25 ра, выход третьего умножителя подключен к первым входам первого и второго делителей, выходы которых подключены к вторым входам соответственно первого и второго умножителей, i-й (i=1, 3) вход третьего сумматора соединен с входом i-го элемента НЕ группы и является i м входом узла перехода, i-м выходом которого является выход i-го сумматора, вторые входы первого и второго делителей являются соответственно четвертым и пятым входами узла перехода, шестым входом которого является первый вход третьего умножителя, второй вход которого подключен к выходу регистра числа.

1387173

Щг7

VE

Составитель A. Баранов

Редактор М. Бланар Техред И. Верес Корректор A. Тяско

Заказ 1230/56 Тираж 928 Подписное

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

l! 3035. Москва. Ж вЂ” 35, Раушская наб., д. 4/5

Производственно-полиграфическое предприятие, г. Ужгород, ул. Проектная, 4