Информационные технологии
УДК 621.391
СИНТЕЗ ЦИФРОВЫХ КИХ-ФИЛЬТРОВ ДЛЯ РЕШЕНИЯ ЗАДАЧ ВОССТАНОВЛЕНИЯ СИГНАЛОВ С ИСПОЛЬЗОВАНИЕМ КРИТЕРИЯ МОМЕНТОВ
В.И. Батищев, И.И. Волков, А.Г. Золин
Самарский государственный технический университет 443100, г. Самара, ул. Молодогвардейская, 244
Е-майл: vib@list.ru
Рассмотрен метод построения обратных цифровых фильтров для решения обратных задач восстановления сигналов на основе аппроксимационного подхода. В качестве критерия адекватности используется критерий моментов.
Ключевые слова: критерий моментов, КИХ-фильтр, весовая функция, обратный фильтр.
В практике обработки и интерпретации экспериментальных данных часто возникает необходимость рассмотрения обратной задачи, заключающейся в восстановлении неизвестного входного воздействия по результатам регистрации откликов на выходе средств измерения. В большинстве случаев это задача компенсации искажающего действия аппаратной функции, обеспечивающая улучшение разрешающей способности различного рода измерительных приборов и систем.
На сегодняшний день методология синтеза оптимальных алгоритмов восстановления сигналов разработана достаточно полно. Однако существующие методы либо требуют для своей реализации не всегда доступной априорной информации, либо сталкиваются с вычислительными проблемами, связанными с некорректностью обратных задач и необходимостью использования регуляризующих процедур [1, 2].
В статье предлагаются алгоритмы синтеза цифровых фильтров для решения задач восстановления сигналов, основанные на использовании метода моментов, свойства и достоинства которого известны [3, 4].
Наблюдается временной ряд у(тп), который можно представить в виде преобразованного КИХ-фильтром с весовой функцией й0(г) временного ряда л(ш) . Ставится задача синтеза обратного КИХ-фильтра с такой функцией Н(Г), чтобы при подаче на его вход временного ряда у(ш) его выходной сигнал был оценкой сигнала л(ш) .
Виталий Иванович Батищев (д.т.н., проф.), заведующий кафедрой «Информационные технологии».
Игорь Иванович Волков (к.т.н., доц.), доцент кафедры «Информационные технологии». Алексей Георгиевич Золин (к.т.н., доц.), доцент кафедры «Информационные технологии».
N 0-1
В этом случае
У(т) = 2 h0(i)x(m - (1)
i=0
N-1
X(m) = 2h(i) y(m - i). (2)
i=0
Подставив в (2) y(m) из (1), получим
N0+N-2
X(m) = 2 H (i) x(m - i). (3)
i=0
Здесь H(i) - весовая функция КИХ-фильтра, преобразующего исходный временной ряд x(m) в его оценку.
Из соотношения (3) следует, что для идеального решения задачи (X(m) = x(m)) функция H (i) должна представлять собой символ Кронекера
H (i ) = S(i). (4)
Так как вид функции H(i) целиком и полностью определяется видом весовой функции обратного фильтра h(i), то эта функция должна выбираться такой, чтобы H (i) была как можно ближе к символу Кронекера.
Это приближение предлагается осуществлять по критерию моментов.
Пусть выбрана любая система линейно независимых функций (i),k = 0, p,
таких, чтобы функция ^ (i + j) была функцией с разделяющимися переменными i и j. Такими функциями, в частности, являются полиномы, экспоненты, синусы, косинусы и т. п.
Будем называть К-м моментом символа Кронекера величины
N0+N-2
у к = 2 *(0^ (i) = ^к (0), (5)
i=0
а К-м моментом весовой функции H
N 0 + N - 2
Мк = 2 h(o^k (i). (6)
i=0
Функцию ht будем выбирать такой, чтобы выполнялось условие
Ук = Мк, к = 0, p. (7)
В этом случае
lim H(i) = S(i). (8)
Условие (8) с учетом (6) и (5) сводится к такому:
N 0 + N - 2
2 H (i)4k (0 = ^к (0), к = 0, р, (9)
i=0
получаем
N0 + N - 2 N _^о-1
^ (г)Н(г) = XI ^)Н(г)^к (г + V).
г=0 г=0 V=0
Подстановка полученного результата в (9) дает систему уравнений, из которых должна определиться весовая функция обратного фильтра:
^ _^-1 _
I Хк0^)Н(г)'¥к(г + v) = (0), к = 0,р. (10)
V=0 !=0
Теперь решим систему уравнений (10) для различных видов базисных функций. Полиномиальные базисные функции.
ВД = гк • (11)
В этом случае соотношение (10) принимает вид
N0 -1N-1 I Iho(v)h(i)(i + V)к = 8к,
V=0 !=0
или, после преобразований,
к
1сЫк _ Ч)ад = 8к • (12)
д=0
Здесь
No -1
«0^) = I *4(0; (13)
=0
N-1
«V = Х'ЖО; (14)
г=0
V - е - моменты исходного и обратного фильтров.
Так как ^(О известна, известны и моменты «0^ )• Из системы уравнений (12) могут быть определены моменты ц обратного фильтра:
1
«0 =
0>(0)
1 к-1
а =---'
к /л\ к 0
0>(0) ?=0
1сХ(к - Ч)«9, к = 1, р.
(15)
у=0
Теперь задача сводится к синтезу обратного фильтра по критерию моментов
V-1 ^
ХД(г)гv = а„, V = 0,1. (16)
=0
Исследования показали, что для получения наилучшего результата наиболее эффективным способом весовую функцию обратного фильтра следует выбирать как взвешенную сумму исходных базисных функций 100
Н(г ) = 1СкС,к+к. (17)
к=0
В этом случае система уравнений (12) примет вид
р N-1
1Ск1С+/' = «V, V = 0, р. (18)
к=0 г=0
Из этой системы уравнений определяется величина Ск :
С = 1 (-1)к+т (р + к +1)!(р + т +1)!(N - т - 1)!^т+1 к т=0 к!т!т!(р - к)!(р - т)!(т + к +1)^ + к)!,
где
^ = а
X 0 ^0
X, =2 (-і)'"-” Р.,ам
V-
=1
»0-1
00
і=0
N-1
і=0 °с0
а0(к )
Опять пришли к синтезу фильтра по критерию моментов. Выбираем
!,(і) = 2с,А<”'" . (23)
V-
0
Подставим к(і) из (23) в (21), будем иметь
Я(к )
Решение данной системы находится по алгоритму:
(19)
Экспоненциальные базисные функции.
^ (г) = А(к+1)!';
Х = е«\
Сразу обратим внимание на то, что для этих параметрических функций
«0(к)Мк =1, к = 0,р, (20)
где
(21)
цк = 1*(0А<* *'>'.
г =0
Из (20) и (21) получим следующую систему уравнений:
N-1 1 „__
ХА(0А(*+1)г =-----------------------------------------, к = 0, р. (22)
Р 1 +г+ 2) N 1
2Cv 1 Я+У+2 = ’к = °’Р- (24)
^0 1 - Я А0(к)
1 Я к+"^+2)N V-1
к ,v
1 - Я
=
= 0, р, = 0, к
к+v+2
-2^,, ^к,д
д=0
к(ш,т) =1
"
к(т, д) = - 2 к("’к)ск,,
к = д + 1 " =0, р
д =0," - 1
с = 2^"^ 2А(",к), V=р0.
v £ ^"," к=0 а0(к) ^
(25)
Г армонические базисные функции с периодом N.
Данные базисные функции имеют вид
^к (і) =
2кіп , ;--
sm-------------, к = 1, р;
cos-
N
2кіж
N
(26)
к = 0, р.
В этом случае система уравнений (10) примет вид
^( ) ■ (2к(г + V>^) 0 1-
I ) = °, к = 1,р;
у=0 г=0 N
I Ih(i>ho (V) со8(2к(г+-^) =1, к = 0, р.
у=0 г=0
N
После тригонометрических преобразований приведем эту систему к виду
а0(к )Д + А)(к )ак = 0, к = 1, р;
А>(0)А>=1 _
А>(к )Рк-а0(к )ак =1, к = 1, Р.
(27)
Здесь
No 1 2кіж
а0(к) = 2 Л0(г>т
і=0
No -1
А>(к) = 2 Л0(і)с^
і=0
N
2кіп
N
N-1
а
=2*(і)
БІП-
і=0
N-1
& = 2А(і)біп
2 кіл
~н~;
2кіл
=0 N
Величины «0(к) и Д)(к) известны, т. к. определяются по известной весовой функции ^0 (0. Величины сск и Рк находим из системы уравнений (27)
1
$0(0)
А)(к )
а02(к) + $02(к)
(29)
а0(к )
,к = 1, р.
ак = - 2 2
к а0(к) + 00 (к)
Из (29) и (28) приходим к следующей системе уравнений:
'N-1 1
2*(і)
і=0
N-1
2^(0 БІП
і=0
N-1
2^(і) СОБ
00(0)
2 кіл
N
2 кіл
а0(к )
а0 (к) + 002(к)
(30)
00(к) ,к = 1, р.
і=0 N а0 (к) + 002 (к)
Выберем весовую функцию исследуемого фильтра в виде
2кіл ^ . 2кіл
ко=с0+2скС0Б +2Ак біп —— ■
к=1 N к=1 N
(31)
Подставив Н(і) в систему уравнений (31) и приняв во внимание ортогональность базисных функций, после ряда преобразований будем иметь:
с =
1
с =
N00(0)
2 00(к)
N а0 (к) + 00(к)
(32)
Л, = -
2 а°(к) ,к = 1, р
N «0 (к) + 00 (к)
И, наконец, подстановка значений Ск и вк в (31) дает такую весовую функцию обратного фильтра:
к(і) = —
N
1
00(0)
2кіл . 2кіл ^
Р 00(к)С0^^^ -а0(к)біп 77 -22— N N
к=1
ао(к) + 0(2 (к)
0
к
Рассмотрим восстановление числовой последовательности из случайных значений, которая была обработана прямоугольным фильтром с весовой функцией Н0=1/Ы0. Значения исходного и двух восстановленных рядов показаны на рисунке.
Исходный и восстановленные ряды:
1 - исходный ряд;
2 - восстановленный ряд, N0 = 3, ширина обратного фильтра 6, количество моментов 4;
3 - восстановленный ряд, N0 = 5, ширина обратного фильтра 6, количество моментов 4
Относительная средняя квадратическая погрешность (ОСП) восстановления вычислялась по формуле Л=|| х - Х||/|| х ||. Полученные результаты при использовании полиномиальных базисных функций приведены в таблице.
Результаты восстановления последовательности
№ п/п N0 Ширина обратного фильтра Кол-во моментов ОСП № п/п N0 ОСП
1 3 3 2 0,05 11 5 0,06
2 4 2 0,07 12 0,07
3 5 2 0,10 13 0,08
4 4 3 0,055 14 0,07
5 5 3 0,053 15 0,06
6 6 3 0,08 16 0,06
7 5 4 0,06 17 0,07
8 6 4 0,04 18 0,07
9 7 4 0,06 19 0,07
10 9 4 0,11 10 0,07
1. Батищев В.И., Золин А.Г., Косарев Д.Н., Романеев А.Е. Аппроксимационный подход к решению обратных задач анализа и интерпретации экспериментальных данных // Вестник СамГТУ. Сер. Технические науки. - 2006. - Вып. 40. - С. 57-65.
2. Леонов А. С. Решение некорректно поставленных обратных задач. - М.: Либроком, 2010.
3. Батищев В.И., Мелентьев В.С. Аппроксимационные методы и системы промышленных измерений, контроля, испытаний, диагностики. - М.: Машиностроение, 2007. - 393 с.
4. Батищев В.И., Волков И.И., Золин А.Г. Построение и оптимизация ортогональных базисных систем для аппроксимации спектрально-корреляционного анализа и идентификации линейных динамических объектов // Вестник СамГТУ. Сер. Технические науки. - 2007. - Вып. 40. - С. 47-52.
Статья поступила в редакцию 5 сентября 2012 г.
SYNTHESIS OF FIR DIGITAL FILTERS FOR SOLVING SIGNAL RECONSTRUCTION USING CRITERIA OF MOMENTS
V.I. Batishchev, 1.1. Volkov, A. G. Zolin
Samara State Technical University
244, Molodogvardeyskaya st., Samara, 443100
The method of constructing the inverse digital filters for solving inverse problems of reconstruction of signals based on an approximation approach. The criterion used criterion of moments.
Keywords: criterion points FIR filter, the weight function, the inverse filter.
Vitaly I. Batishchev (Dr. Sci. (Techn.)), Professor.
Igor I. Volkov (Ph.D. (Techn.)), Associate Professor. Aleksey G. Zolin (Ph.D. (Techn.)), Associate Professor.