CC BY
67
УДК 621.372
DOI: 10.17586/0021-3454-2016-59-1-69-72
СИНТЕЗ ФИЛЬТРА С ЗАДАННОЙ ЧАСТОТНОЙ ПЕРЕДАТОЧНОЙ ФУНКЦИЕЙ
С. И. Зиатдинов
Санкт-Петербургский государственный университет аэрокосмического приборостроения,
190000, Санкт-Петербург, Россия E-mail: kaJ53@guap.ru
Предложена методика синтеза фильтра с независимыми и практически любыми амплитудно- и фазочастотными характеристиками. Получено выражение для импульсной характеристики фильтра. Приведен конкретный пример.
Ключевые слова: фильтр, частотная передаточная функция, импульсная характеристика.
При обработке информации наиболее часто используется процедура фильтрации сигнала с помощью различных фильтров — фильтров нижних и верхних частот, полосовых и ре-жекторных фильтров. В ряде случаев применяются адаптивные фильтры, параметры которых могут изменяться в зависимости от спектральных характеристик обрабатываемых сигналов.
Вместе с тем управление частотной передаточной функцией фильтра является достаточно сложной технической задачей. Далеко не всегда удается согласовать параметры фильтра с характеристиками сигналов на его входе. При этом существует достаточно жесткая связь между амплитудно-частотной характеристикой (АЧХ) фильтра и его фазочастотной характеристикой (ФЧХ).
Несомненно представляет практический интерес разработка фильтра с заранее заданными АЧХ и ФЧХ, которые в то же время были бы не связаны друг с другом. Решение этой задачи является основным содержанием настоящей статьи.
Пусть желаемые АЧХ и ФЧХ фильтра имеют вид, показанный соответственно на рис. 1, а, б жирными линиями. При этом АЧХ фильтра в области отрицательных частот является зеркальной копией АЧХ в области положительных частот относительно оси W(ю) , а ФЧХ также зеркальна относительно оси ф(ю), но в области отрицательных частот имеет отрицательный знак. а) б)
W(ra)
W(ro,) ф(ю)
Дю
Рис. 1
Для синтеза фильтра с заданными произвольными и независимыми АЧХ и ФЧХ проведем их дискретизацию по частоте с шагом Дю, в пределах которого не происходит существенных изменений Ж(ю) и ф(ю). В результате АЧХ и ФЧХ фильтра можно представить в
виде последовательности прямоугольников высотой Ж(юг-) и ф(юг-) с основанием Дю (см.
рис. 1, тонкие линии).
ю
ю
Очевидно, что чем меньше шаг дискретизации по частоте Дю, тем меньше дискретизи-рованные АЧХ и ФЧХ будут отличаться от исходных АЧХ и ФЧХ фильтра. В результате синтезированный фильтр можно представить в виде параллельно соединенных элементарных фильтров с АЧХ (ю) = Ж (юг) и ФЧХ фг (ю) = ф(юг), где г — номер элементарного фильтра. АЧХ и ФЧХ элементарного фильтра приведены на рис. 2, а, б соответственно.
а) б)
Ж(ю) ф(ю)
Щ(ю,)
ф, ••
0 Ю2 Ю1
Ю
-Ю2 -Ю1
-ф,
Ю2 Ю1
Рис. 2
В целом синтезированный фильтр будет иметь следующую частотную передаточную функцию:
N
W2 О) = Е Щ (ю)е
,=1
]ф , (ю)
(1)
где N — количество элементарных фильтров.
Из данного соотношения следует, что посредством независимого выбора необходимых значений Ж (ю) и фг- (ю) можно получить практически любую частотную передаточную
функцию (ую) .
Импульсная характеристика такого фильтра принимает вид [см. лит.]
1 ю ■ , 1 ю N / \ .
^ (0) = — | (= — I £ Ж (ю>ф(ю)в]Шс1ю . (2)
2п -ю 2п -ю г =1
Поменяв интеграл и сумму в выражении (2) местами, получим
ИеС) = £ 2- ю Ж (юУ[ш+фг (ю)]^ю = £ И (), (3)
г=12п -ю г=1
1 ю
где И (0) =— | (ю)в^ш+фг — импульсная характеристика элементарного фильтра.
2п -ю
Найдем импульсную характеристику элементарного фильтра, считая, что в пределах его полосы пропускания Дю величины Ж(ю) и ф,(ю) остаются постоянными:
1 -ю1 1 ю2 ■
И (0 =— I Щв-]фгв]Шс1ю +— I Цв^в^Ло =
2п
2п
1 _®1 1 ®2 = -L j we'^, d+— j Ще](ш+ф, d.
2п.
2л,
(4)
Используя в формуле (4) преобразование Лапласа e]ф = cos ф + j sin ф , после несложных тригонометрических преобразований получаем
W sin Г11 cos ÍЮ2±Ю1
1 ю2
h, (t) = — j Щ cos(rot + ф, )dro = —L sin I ——— t | cos I ——— t + ф, |. п ю nt V 2 J V 2
(5)
Преобразуем выражение (5) к виду
Ю
0
, , ч 2Wj . Äot
hi (t) = ~r sin — cos^T1 + ф,), (6)
nt 2
Ю2 + Ю1 A TT^ J.
где ю, = —±2—— средняя частота АЧХ элементарного фильтра с полосой пропускания
Дю = ю2 - о»!.
В частном случае можно записать:
- 1п , ,, , 2W, . Дюt .
ф, = + n /2, h- (t) = ± —- sin-sin o, t;
nt 2
- , „ , , ч >/2-W- • Дюt, , • л
ф, = + n /4, hi (t) =-- sin-(cos o,t ± sin o,t);
nt 2
_ . ^ , , ч W, . Дюt ( r-
Ф, = + n /6, h,-(t) = — sin-^v3 - coso,t ± sin ю,^;
nt 2
Л , , Л 2W, . Дюt
Ф, = 0, h, (t) = —- sin-cos o,t.
nt 2
В качестве примера рассмотрим синтез фильтра нижних частот первого порядка с частотной передаточной функцией вида
W (jo) = -j-,
1 + jOT
где T — постоянная времени фильтра.
Данной частотной передаточной функции соответствуют следующие АЧХ и ФЧХ:
W(ю) = , 1 =, ф(ю) = -arctg ют . (7)
V1 + ю2 т2
При этом импульсная характеристика фильтра определяется выражением
h(t) = - е-t/т . т
В соответствии с выражениями (3), (6) и (7) импульсная характеристика синтезируемого фильтра нижних частот принимает вид
N 2 Дю h£ (t) = Z-1 sin-(cos ю^ - arctg ю,т).
,=1 nt^/1 + ю2 т2 2 При расчетах полагалось т = 0,01 с, fmax = ю,max /2n = 10 кГц.
Качество синтезированного фильтра оценивалось по относительной среднеквадратиче-ской ошибке вычисления импульсной характеристики фильтра:
I i M 2
Oh =J M {[h(t) - h£ (t)]/h(t )}
ЦЩ) - пЕ )]/ п(г У2
к=1
где М — число учитываемых отсчетов импульсной характеристики, которое было принято равным 100; отсчеты брались с периодом 0,001 с.
Были получены следующие результаты: при ширине полосы пропускания элементарного фильтра Д/ = Дю /2п = 10 Гц нормированная среднеквадратическая ошибка Оп = 2,32 %; при Д/ = 5 Гц — оп = 0,91 % и при Д/ = 1 Гц — оп = 0,8 %.
Таким образом, с уменьшением полосы пропускания элементарного фильтра частотная передаточная функция синтезированного фильтра стремится к желаемой частотной передаточной функции.
Рассмотренная методика синтеза позволяет получить фильтр практически с любой частотной передаточной функцией.
ЛИТЕРАТУРА
Гоноровский И. С. Радиотехнические цепи и сигналы. М.: Радио и связь, 1986. 512 с.
Сведения об авторе
Сергей Ильич Зиатдинов — д-р техн. наук, профессор; СПбГУАП, кафедра информационно-сетевых технологий; E-mail: kaf53@guap.ru
Ссылка для цитирования: Зиатдинов С. И. Синтез фильтра с заданной частотной передаточной функцией // Изв. вузов. Приборостроение. 2016. Т. 59, № 1. С. 69—72.
SYNTHESIS OF FILTER WITH SPECIFIED FREQUENCY TRANSFER FUNCTION
S. I. Ziatdinov
Saint Petersburg State University of Aerospace Instrumentation, 190000, St. Petersburg, Russia
E-mail: kaf53@guap.ru
A method is proposed for synthesis of filter with independent and practical arbitrary amplitude and phase frequency characteristics. An expression describing impulse characteristic of the filter is derived. An example of filter synthesis with the use of the proposed method is presented.
Keywords: filter, frequency transfer function, impulse characteristic.
Data on author
Sergey I. Ziatdinov — Dr. Sci., Professor; St. Petersburg State University of Aerospace Instrumentation, Department of Information Networks Technologies; E-mail: kaf53@guap.ru
For citation: Ziatdinov S. I. Synthesis of filter with specified frequency transfer function // Izvestiya Vysshikh Uchebnykh Zavedeniy. Priborostroenie. 2016. Vol. 59, N 1. P. 69—72 (in Russian).
DOI: 10.17586/0021-3454-2016-59-1-69-72
Рекомендована кафедрой информационно-сетевых технологий
Поступила в редакцию 11.03.14 г.
