Fir digital filters synthesis for solving signal reconstruction using two criteria

Abstract


The design method of the inverse nonlinear time-dependent digital filters for the solution of inverse problems of the recovery signals is given. Recovery method shown possible, in case of application to the original signal of the filter with known finite impulse response. Finding the impulse response of the inverse filter is based on the approximation approach. As a criterion for the adequacy of the solutions two criteria are considered: the criterion of complete coincidence of the original with the reconstructed signal and the criterion of the minimum weighted square error. Algorithm tests by restoring the signals in a variety of analytic functions are implemented .

Full Text

В настоящее время при решении задач обработки и интерпретации экспериментальных данных часто возникает необходимость рассмотрения обратной задачи, заключающейся в восстановлении неизвестного входного воздействия по результатам регистрации откликов на выходе средств измерения. В большинстве случаев это задача компенсации искажающего действия аппаратной функции, обеспечивающая улучшение разрешающей способности различного рода измерительных приборов и систем. В случае, когда для обработки доступна только часть искаженного сигнала, без начальных условий, задача становится недоопределенной и, соответственно, некорректной. На сегодняшний день методология синтеза оптимальных алгоритмов восстановления сигналов разработана достаточно полно. Однако существующие методы либо требуют для своей реализации не всегда доступной априорной информации, либо сталкиваются с вычислительными проблемами, связанными с некорректностью обратных задач и необходимостью использования регуляризующих процедур [1, 2]. В статье предлагаются алгоритмы синтеза нелинейных нестационарных (адаптивных) цифровых фильтров для решения задач восстановления сигналов, основанные на построении модели весовой функции обратного фильтра с использованием двух подходов: критерия полного соответствия и критерия минимума взвешенной среднеквадратической погрешности [3, 4]. Пусть наблюдается ряд , (1) представляющий собой ряд x(j), обработанный фильтром с конечной импульсной характеристикой (КИХ-фильтром) с известной весовой функцией h0(i). Требуется восстановить ряд x(j) по ряду y(j). Начальные условия неизвестны. Оценку ряда х в момент времени m будем искать по алгоритму (2) где р - величина участка ряда у, по которой определяется оценка х. Подставив в (1) вместо x(j) x’(j) из (2), получим . (3) Значения y’(j) при идеальной обратной фильтрации будут совпадать со значениями y(j), а при реальной будут отличаться от них. Весовая функция h(v) обратного фильтра должна выбираться такой, чтобы это отличие было наименьшим. Для простоты перепишем соотношение (3) в следующем виде: . Введем обозначение (4) Тогда (5) Отметим, что значения D(k) будут зависеть от значения текущего момента времени. Синтезируем алгоритмы вычисления весовой функции h(v) обратного фильтра по различным критериям. Критерий совпадения или . (6) Из (6) с учетом (5) получим систему уравнений для вычисления значений h(v) . (7) Значения h(v), найденные из этой системы уравнений, будут зависеть от значений текущего момента времени m и от y(j). Решение системы уравнений (7) может быть осуществлено по одному из следующих алгоритмов. Общая часть: Алгоритм 1: . (8) Алгоритм 2: (9) Критерий минимума взвешенной квадратичной погрешности Взвешенная квадратичная погрешность определяется соотношением (10) Для обеспечения минимума этой погрешности должно быть выполнено условие . Это условие с учетом (10) сводится к системе уравнений Для простоты введем обозначения: (11) Тогда получим такую систему уравнений: (12) Для ее решения получаем следующий алгоритм: (13) Апробация алгоритмов Исследуя алгоритмы (8) и (13) на примерах, для простых случаев можно получить аналитические выражения весовой функции обратного фильтра. Рассмотрим в качестве исходного сигнала х различные виды элементарных функций. В качестве прямого фильтра возьмем фильтр с весовой функцией h0(k)=1/3, k={0,1,2}. Показательная функция: . Используя алгоритм (8), для данного случая можно аналитически получить значение весовой функции обратного фильтра что дает нам абсолютно точное решение задачи. Полиномиальная функция: Для данного случая получим значение весовой функции обратного фильтра: Применение обратного фильтра с такой весовой функцией даст нам абсолютно точное решение задачи. Результаты восстановления x(m) Тригонометрическая функция: Используя алгоритм (8), получим значение весовой функции обратного фильтра: И в этом случае мы полиучаем абсолютно точное решение задачи. Рассмотрим более сложный пример: Весовая функция прямого фильтра h0={ 0.1,0.3,0.5,1,1,1,1,1,1,0.5,0.3,0.1}. Результат восстановления по алгоритму (8) представлен на рисунке. Относительная среднеквадратическая погрешность восстановления, вычисляемая по формуле , в случае использования алгоритмов (8) и (13) без учета переходного процесса в начале восстановления составила 0,25%.

About the authors

Vitaly I Batishchev

Samara State Technical University

244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
(Dr. Sci. (Techn.)), Professor

Igor I Volkov

Samara State Technical University

244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
(Ph.D. (Techn.)), Associate Professor

Aleksey G Zolin

Samara State Technical University

Email: zolin.a.g@gmail.com
244, Molodogvardeyskaya st., Samara, 443100, Russian Federation
(Ph.D. (Techn.)), Associate Professor

References

  1. Батищев В.И., Золин А.Г., Косарев Д.Н., Романеев А.Е. Аппроксимационный подход к решению обратных задач анализа и интерпретации экспериментальных данных // Вестник Самарского государственного технического университета. Сер. Технические науки. - 2006. - № 40. - С. 57-65.
  2. Леонов А.С. Решение некорректно поставленных обратных задач. - М.: Либроком, 2010.
  3. Батищев В.И., Мелентьев В.С. Аппроксимационные методы и системы промышленных измерений, контроля, испытаний, диагностики. - М.: Машиностроение, 2007. - 393 с.
  4. Батищев В.И., Волков И.И., Золин А.Г. Построение и оптимизация ортогональных базисных систем для аппроксимации спектрально-корреляционного анализа и идентификации линейных динамических объектов // Вестник Самарского государственного технического университета. Сер. Технические науки. - 2007. - № 40. - С. 47-52.

Statistics

Views

Abstract - 42

PDF (Russian) - 7

Cited-By


Article Metrics

Metrics Loading ...

PlumX

Dimensions

Refbacks

  • There are currently no refbacks.

Copyright (c) 2014 Samara State Technical University

Creative Commons License
This work is licensed under a Creative Commons Attribution 4.0 International License.

This website uses cookies

You consent to our cookies if you continue to use our website.

About Cookies