КРАТКИЕ СООБЩЕНИЯ X
УДК 621.396:681.323 с1о1:10.15217/1ззп1684-8853.2016.5.98
СИНТЕЗ НЕРЕКУРСИВНЫХ ДИСКРЕТНЫХ ФИЛЬТРОВ ВО ВРЕМЕННОЙ ОБЛАСТИ
С. И. Зиатдинова, доктор техн. наук, профессор
аСанкт-Петербургский государственный университет аэрокосмического приборостроения, Санкт-Петербург, РФ
Введение: переход от непрерывных фильтров к нерекурсивным дискретным фильтрам с использованием отсчетов импульсных характеристик непрерывных фильтров-аналогов приводит к ошибкам в формировании переходных процессов. Вместе с тем динамические свойства фильтров в полном объеме определяются их переходными характеристиками. Синтез нерекурсивных дискретных фильтров на базе переходных характеристик позволяет обеспечить их динамические свойства, совпадающие с динамическими свойствами непрерывных фильтров. Цель: разработка методики синтеза нерекурсивных дискретных фильтров с использованием переходных характеристик. Методы: представление импульсной характеристики дискретного фильтра в виде первой конечной разности его переходной характеристики и нахождение с использованием интеграла наложения выходного сигнала. Результаты: получена новая методика синтеза нерекурсивного дискретного фильтра, позволяющая при заданном входном сигнале и известной переходной характеристике фильтра найти его выходной сигнал. Выдвинутые теоретические положения подтверждены конкретным примером. Практическая значимость: полученная методика позволяет при решении задач фильтрации сигналов в помехах и шумах использовать нерекурсивные дискретные фильтры, динамические свойства которых совпадают с динамическими свойствами непрерывных фильтров-аналогов.
Ключевые слова — частотная передаточная функция, импульсная и переходная характеристики, нерекурсивные дискретные фильтры, разностные уравнения, коэффициенты.
Введение
При обработке сигналов самое широкое распространение получили разнообразные фильтры, с помощью которых решаются задачи фильтрации, дифференцирования, интегрирования, экстраполяции и т. д. В каждом конкретном случае фильтр должен обладать определенными частотными свойствами.
В современных условиях все чаще приходится сталкиваться с цифровыми методами обработки на базе персональных компьютеров или специализированных вычислителей. При этом стоит задача преобразования непрерывных фильтров в дискретные. В настоящее время достаточно хорошо отработана методика синтеза дискретных фильтров по их непрерывным аналогам.
Можно выделить два основных метода синтеза: синтез дискретных фильтров в частотной области и синтез дискретных фильтров во временной области. При синтезе дискретных фильтров в частотной области [1] должны воспроизводиться амплитудно-частотные и фазо-частотные характеристики непрерывных фильтров с минимальными погрешностями. Преобразование частотной передаточной функции непрерывного фильтра в частотную передаточную функцию дискретного фильтра осуществляется на базе билинейного преобразования, которое используется для создания в основном фильтров верхних частот и режекторных фильтров.
В случае синтеза дискретных фильтров во временной области применяется метод инвариант-
ной импульсной характеристики, при котором отсчеты импульсной характеристики непрерывного фильтра используются для вычисления коэффициентов линейного разностного уравнения дискретного фильтра.
Для дискретных фильтров импульсная характеристика представляется последовательностью масштабированных отсчетов непрерывной импульсной характеристики Н(£) [2]:
ЛИ = тн(г),
где I = 0, 1, 2, ... — номер отсчета импульсной характеристики; Т — период следования отсчетов входного и выходного сигналов дискретного фильтра; ^ = 1Т — текущее дискретное время.
В результате алгоритм работы дискретного фильтра определяется дискретной сверткой отсчетов обрабатываемого сигнала х(£) и импульсной характеристики Н(£) фильтра [2, 3] к
у[к] = X х[Ь]Н[к - ¿], & = 0, 1, 2, ... .
1=0
Выходной сигнал дискретного фильтра в общем виде можно описать линейным разностным уравнением с постоянными коэффициентами:
у[к] = аох[к] + а-х[к -1] +... + апх[к - п] -- Ь-у[к -1] - Ь2у[к - 2] -... - Ьпу[к -п] =
п п
= Х а1х[к - ¿] -X Ьу[к - ¿], (1)
1=0 1=1
где п — порядок разностного уравнения; ар Ъ1 — постоянные коэффициенты.
КРАТКИЕ СООБЩЕНИЯ X
При этом синтез дискретного фильтра для заданного порядка п заключается в выборе постоянных коэффициентов аь, Ь, которые определяют вид частотной характеристики фильтра.
На практике широко используются как нерекурсивные, так и рекурсивные дискретные фильтры.
Нерекурсивные дискретные фильтры
Для нерекурсивных дискретных фильтров разностное уравнение (1) записывается следующим образом:
у[к] = а0 х[к] + а^[к -1] +... + апх[к - п] =
= ^ а^х[к - ¿].
1=0
(2)
При этом порядок уравнения п рассчитывается так, чтобы произведение пТ было больше длительности переходного процесса в фильтре пТ > а коэффициенты аь принимаются равными отсчетам импульсной характеристики Лр]:
а0 = А[0], а1 = А[1], ..., ап = А[п].
Тогда выражение (2) для выходного сигнала нерекурсивного дискретного фильтра принимает хорошо известный вид
у[к] = Н[0]х[к] + Н[1]х[к -1] +... + Н[п]х[к - п] =
(3)
= ^ Н[Ь]х[к - £].
1=0
Покажем, что в общем случае данное соотношение является неверным, т. е. не отражает точно истинные физические процессы, протекающие в дискретном фильтре.
В качестве примера рассмотрим фильтр нижних частот (ФНЧ) первого порядка с частотной передаточной функцией в непрерывном варианте
W (]ю) = 1
1 + ]'юх
где х — постоянная времени фильтра.
Данной частотной передаточной функции соответствуют импульсная и переходная характеристики вида [1]
2п 3 1 + ]Юх ср
= 1 - е ®ср*,
(4)
где ю = — частота среза фильтра. х
В соответствии с (3) переходная характеристика рассматриваемого нерекурсивного дискретного фильтра, как его реакция на единичное
входное воздействие, записывается следующим образом:
g[k] = Н[0]1[к] + Н[1]1[к -1] +... + Н[п]1х[к - п] =
= Ё ^[¿]1[к -1],
1=0
(5)
где 1[к] — единичный решетчатый входной сигнал.
Синтезированный нерекурсивный дискретный фильтр должен по своим параметрам соответствовать непрерывному фильтру, т. е. переходная характеристика дискретного фильтра (5) должна с точностью до постоянного множителя равняться переходной характеристике непрерывного фильтра (4) в дискретные моменты времени
г1 = ьт.
На рисунке показаны переходные характеристики рассматриваемого непрерывного ФНЧ £н[к] и нерекурсивного дискретного ФНЧ £д[к] для случая Тюср = 0,4. Здесь же представлено отклонение переходных характеристик непрерывного и дискретного фильтров Дg[k] = £дИ - £н[к].
Анализируя представленные результаты расчетов, можно отметить, что для ФНЧ первого порядка на протяжении длительности переходного процесса ^ » 8Т отклонение переходной характеристики дискретного фильтра от переходной характеристики непрерывного фильтра изменяется от 40 % практически до нуля относительно установившегося значения = 1. Расчеты, проведенные для ФНЧ второго порядка, показывают, что указанное отклонение переходных характеристик составляет 10 %.
Таким образом, существующая методика определения коэффициентов разностного уравнения нерекурсивных дискретных фильтров на базе импульсной характеристики является неточной и в результате не обеспечивает требуемого качества работы дискретных фильтров в пределах переходного процесса.
1
0,8 0,6
0,2 0
.-V /
/7
к ' / г / /
/
-/ /
I
I \
\
8нМ
Д*[Л]
02
6 8 10 12 14
■ Переходные характеристики ФНЧ и отклонение переходных характеристик
к
№ 5, 2016 X ИНФОРМАЦИОННО-УПРАВЛЯЮЩИЕ СИСТЕМЫ X 99
КРАТКИЕ СООБЩЕНИЯ
В полном объеме физические процессы, протекающие в фильтрах, описываются их переходными характеристиками. В связи с этим рассмотрим синтез нерекурсивных дискретных фильтров на базе их переходных характеристик.
В общем виде выходной сигнал нерекурсивного фильтра определяется интегралом наложения t
y(t) = J x(t -x)h(x)dx. (6)
0
При этом импульсная характеристика фильтра h(t) связана с его переходной характеристикой g(t) соотношением
h(t) = Ml = lim ¿(t) - g(t -At). (7)
dt At^o At
После подстановки (7) в (6) получим
y(t) = lim fx(t-x)^(T)-g(T-At)
At^Oo 0
At
dx.
(8)
Будем считать, что за время Дг не происходит заметных изменений переходной характеристики. Тогда, положив Дг = Т и заменяя интеграл в (8) суммой, получим
у( к) =Хх [ к -1] {я [I ] - я[ I -1]} х [ к -1 ] Ая [I ], (9)
1=0 1=0
где Д^р] = g[i] - g[i - 1]; g[i] — отсчеты переходной характеристики фильтра.
Для & = 0, 1, 2, ... на основании (9) можно записать следующую систему уравнений:
у[0] = х[0]Д^[0]; у[1] = х[1]Д^[0] + х[0]Д^[1]; у[2] = х[2]Д^[0] + х[1]Д£[1] + х[0]Д^[2];
у[Щ = х[&]Д£[0] + х[к - 1]Д^[1] + ... + х[0]Д^&]. Представим данную систему в виде у[0] = а0х[0]; У[1] = а0х[1] + а1х[0]; у[2] = а0х[2] + а1х[1] + а2х[0];
............................................................... (10)
у[&] = а0х[&] + а1х[& - 1] + ... + а&х[0],
где коэффициенты а0 = Д^[0], а1 = Д^[1], ..., =
= Дg[Щ.
Расчет переходных характеристик нерекурсивных дискретных фильтров как нижних, так и верхних частот первого порядка с коэффициентами, вычисленными с использованием соотношений (10), показал совпадение с переходными характеристиками соответствующих непрерывных фильтров на всем протяжении переходных характеристик. Для фильтров более высоких порядков при Тюср = 0,4 отклонение переходных характеристик непрерывных и дискретных фильтров не превышает десятых долей процентов и резко уменьшается с уменьшением произведения Тюср.
Заключение
Известная методика расчетов коэффициентов разностного уравнения нерекурсивных дискретных фильтров на базе отсчетов импульсной характеристики не обеспечивает правильной реализации переходных процессов.
Методика расчетов коэффициентов с использованием отсчетов переходной характеристики позволяет получить динамические свойства нерекурсивных дискретных фильтров, совпадающие с динамическими свойствами непрерывных фильтров-аналогов как в переходном, так и в установившемся режимах.
Предложенная методика является общей и может быть распространена на нерекурсивные и рекурсивные дискретные фильтры как нижних, так и верхних частот практически любых порядков.
Литература
1. Воробьев С. Н. Цифровая обработка сигналов. — М.: Академия, 2013. — 318 с.
2. Гоноровский И. С. Радиотехнические цепи и сигналы. — М.: Радио и связь, 1986. — 512 с.
3. Alan V. Oppenheim, Ronald W. Schafer. Discrete-Time Signal Processnig. — Prentice Hall, 1989. — 1120 p.
^ КРАТКИЕ СООБЩЕНИЯ X
UDC 621.396:681.323 doi:10.15217/issn1684-8853.2016.5.98
Synthesis of Non-Recursive Discrete Filters in Temporal Range
Ziatdinov S. I.a, Dr. Sc., Tech., Professor, kaf53@guap.ru
aSaint-Petersburg State University of Aerospace Instrumentation, 67, B. Morskaia St., Saint-Petersburg, Russian Federation
Introduction: The transition from continuous filters to non-recursive discrete filters using the readings of impulsive characteristics of continuous filters-analogs leads to mistakes in forming transient processes. At the same time, the dynamic properties of the filters are fully determined by their transient characteristics. Synthesis of non-recursive discrete filters based on transient characteristics provides that their dynamic properties coincide with the dynamic properties of continuous filters. Purpose: The goal is to create a method for the synthesis of non-recursive discrete filters using transient characteristics. Method: We represent the impulsive characteristic of a discrete filter as the first final difference of its transient characteristic, and find the output signal using an imposition integral. Results: A new technique has been obtained for the synthesis of non-recursive discrete filter, which allows you to calculate the output signal of a filter when its input signal is given and the transient characteristic known. The theoretical results are substantiated by an example. Practical relevance: The obtained technique allows you to solve the problems of filtering signals in noises by using non-recursive discrete filters whose dynamic properties coincide with the dynamic properties of continuous filters-analogs.
Keywords — Frequency Transmission Function, Impulse and Transient Characteristic, Non-Recursive Discrete Filters, Difference Equation, Coefficients.
References
1. Vorobiev C. N. Tsifrovaia obrabotka signalov [Digital Signal Processing]. Moscow, Akademiia Publ., 2013. 318 p. (In Russian).
2. Gonorovskii I. S. Radiotekhnicheskie tsepi i signaly [Radio Circuits and Signals]. Moscow, Radio i sviaz' Publ., 1986. 512 p. (In Russian).
3. Alan V. Oppenheim, Ronald W. Schafer. Discrete-Time Signal Processnig. Prentice Hall, 1989. 1120 p.
УВАЖАЕМЫЕ АВТОРЫ!
Научная электронная библиотека (НЭБ) продолжает работу по реализации проекта SCIENCE INDEX. После того как Вы зарегистрируетесь на сайте НЭБ (http://elibrary.ru/ defaultx.asp), будет создана Ваша личная страничка, содержание которой составят не только Ваши персональные данные, но и перечень всех Ваших печатных трудов, имеющихся в базе данных НЭБ, включая диссертации, патенты и тезисы к конференциям, а также сравнительные индексы цитирования: РИНЦ (Российский индекс научного цитирования), h (индекс Хирша) от Web of Science и h от Scopus. После создания базового варианта Вашей персональной страницы Вы получите код доступа, который позволит Вам редактировать информацию, помогая создавать максимально объективную картину Вашей научной активности и цитирования Ваших трудов.
№ 5, 2016 ^
ИНФОРМАЦИОННО-УПРАВЛЯЮЩИЕ СИСТЕМЫ
101