лекции / TsOS_Uchebnoe_posobie_2018
.pdfфункции цифрового фильтра. |
Операторная передаточная функция |
||||
аналогового фильтра имеет вид: |
H p |
1 |
|
. Нормированное значение |
|
|
|
||||
p 1 |
|||||
|
|
|
|||
частоты дискретизации принять |
Fs 0.25 . |
|
|||
4.Частота фильтра прототипа равна 500 кГц, частота дискретизации Fд 250 кГц. Определите нормированную частоту цифрового фильтра,
если известно, что переход от фильтра прототипа к цифровому фильтру был осуществлён с помощью билинейного z - преобразования.
5.Значение нормированной частоты некоторой точки АЧХ цифрового фильтра равна ˆ / 2 рад, частота дискретизации Fд 100 кГц. Определите соответствующую частоту фильтра прототипа, если известно, что переход от фильтра прототипа к цифровому фильтру был осуществлён с помощью билинейного z - преобразования.
1.2.7Примеры решения задач
1.Методом билинейного z - преобразования перейти от
операторной передаточной функции аналогового фильтра прототипа к передаточной функции цифрового фильтра. Операторная передаточная функция аналогового фильтра имеет
|
вид: H p |
1 |
|
. Частоту дискретизации принять Fд 1 |
Гц. |
|||||||||||||
|
|
|
|
|
|
|||||||||||||
|
p 10 |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
Решение |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
H p |
1 |
|
|
|
|
|
|
|||||
Подставим |
|
в |
|
|
|
выражение |
для |
билинейного |
|
z - |
||||||||
|
p 10 |
|||||||||||||||||
преобразования: |
p 2F |
1 z 1 |
: |
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
s 1 z 1 |
|
|
|
|
|
|
|
|
|||||
H z |
|
1 |
|
|
|
|
|
|
|
|
|
1 z 1 |
|
1 z 1 |
|
|
||
|
|
|
|
|
|
|
|
|
. |
|||||||||
|
1 z 1 |
|
|
|
2Fд 1 z 1 10 1 z 1 |
10 2Fд z 1 10 2Fд |
||||||||||||
|
2Fд |
|
|
|
10 |
|
|
|
|
|
|
|
|
|||||
|
1 z |
1 |
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Подставляя |
значение |
|
|
Fд 1 Гц |
в вышенаписанное |
выражение |
для |
|||||||||||
передаточной функции цифрового фильтра H z , получим:
H z |
1 z 1 |
. |
|
||
|
8z 1 12 |
|
2.Частота фильтра прототипа равна 1500 кГц, частота дискретизации Fд 750 кГц. Определите нормированную
51
частоту цифрового фильтра, если известно, что переход от фильтра прототипа к цифровому фильтру был осуществлён с помощью билинейного z - преобразования.
Решение
Для решения данной задачи подставим заданные значения и Fд в выражение (67) и произведём расчёт:
ˆ |
|
|
1500 |
|
2arctg 1 2 |
|
|
|
рад. |
|
2arctg |
|
|
|
|
|
|||
2arctg |
|
|
|
|
|
||||
|
2Fд |
2 750 |
|
|
4 |
|
2 |
|
|
52
1.3 Расчет цифровых БИХ-фильтров с внесением предыскажений
1.3.1Внесение предыскажений
Вобщем случае, билинейное z -преобразование, заключается в подстановке в передаточную функцию аналогового фильтра прототипа выражения
p |
2 1 |
z 1 |
(70) |
|||
|
|
|
|
|||
T 1 |
z 1 |
|||||
|
|
|||||
где T – интервал дискретизации сигналов, обрабатываемых фильтром, z – комплексный параметр Z-преобразования.
В ходе проведения билинейного z - преобразования происходит искажение (деформация) частотной оси. Частоты АФП и ЦФ связаны соотношениями:
|
2 |
|
|
|
T |
, |
(71) |
|
|
|
arctg |
|
|
||||
|
|
|||||||
|
T |
|
|
|
2 |
|
|
|
|
2 |
T |
, |
|
(72) |
|||
|
|
tg |
|
|
||||
|
|
|
||||||
|
|
|
T |
|
2 |
|
|
|
где - частота ЧХ аналогового фильтра, - частота ЧХ цифрового фильтра. Из-за этой деформации для получения передаточной функции цифрового фильтра с требуемой частотой среза (граничной частотой полосы пропускания) 1 необходимо использовать передаточную функцию аналогового фильтра, имеющего специально подобранную граничную частоту полосы пропускания 1 . Однако мы располагаем таблицами нормализованных передаточных функций аналоговых фильтров, т.е. имеющих 1 1 рад/с. Кайзер и Гоулден преодолевают эту трудность, используя преобразование нижних частот (т.е. заменяя нормализованную переменную p на p / 1 перед
выполнением преобразования (70)).
Относительно билинейного преобразования (70) можно сделать два следующих замечания:
1.действительный постоянный множитель 2 / T , как отметили Рейдер и Голд, можно опустить;
2.процесс двукратного преобразования (от нормализованного аналогового фильтра нижних частот к аналоговому фильтру нижних частот с частотой среза 1 и затем билинейное преобразование)
удваивает объём вычислений.
Рейдер и Голд учли, что константу 2 / T можно опустить и использовали преобразование в виде
53
p |
1 z 1 |
, |
(73) |
|||
|
|
|
||||
1 |
z 1 |
|||||
|
|
|
||||
для которого второе из сделанных Поскольку любой постоянный преобразования, можно использовать
выше замечаний также справедливо. множитель не изменяет формы преобразование вида
|
p k |
1 z 1 |
, |
|
(74) |
|
|
1 z 1 |
|
||||
|
|
|
|
|
||
где k - действительная константа, определяемая из условия |
|
|||||
1 |
|
|
T |
|
|
|
k tg |
1 |
|
, |
(75) |
||
|
|
|
2 |
|
|
|
где в свою очередь 1 - частота среза аналогового фильтра (граничная частота полосы пропускания), 1 - требуемая частота среза цифрового фильтра. Из (75) получим
|
T |
|
|
|
|
|
(76) |
|||
k ctg |
1 |
|
|
. |
|
|
|
|
||
|
|
2 |
|
|
|
|
|
|
|
|
Следовательно, преобразование (74) принимает форму |
|
|||||||||
|
T 1 z 1 |
|
||||||||
p 1ctg |
1 |
|
|
|
|
|
|
, |
(77) |
|
2 |
|
|
|
1 |
||||||
|
1 z |
|
|
|
|
|||||
а при 1 1 |
|
|
|
|
|
|
|
|
|
|
T |
1 z 1 |
|
||||||||
p ctg |
1 |
|
|
|
|
. |
(78) |
|||
2 |
|
z |
1 |
|||||||
|
1 |
|
|
|
|
|
||||
Таким образом, введя в преобразование константу общего вида k , мы избавились от необходимости выполнять его дважды.
Итак, преобразование (77) даёт возможность рассчитать цифровой фильтр нижних частот с заданной частотой среза по данным аналогового фильтра нижних частот.
При синтезе цифрового фильтра верхних частот, полосового или режекторного фильтров Кайзер предлагает до выполнения билинейного преобразования хорошо известными в аналоговой фильтрации методами преобразовать нормализованный аналоговый фильтр нижних частот. Эта методика потребует дополнительных вычислений, однако можно ввести новое преобразование, позволяющее рассчитывать цифровой фильтр (ФНЧ, ФВЧ, ПФ, РФ) непосредственно по передаточной функции аналогового фильтра нижних частот.
Рассмотрим синтез фильтра верхних частот, для которого можно непосредственно воспользоваться уже полученным преобразованием (73). При замене аргумента p передаточной функции данного аналогового фильтра
нижних частот на 1 / p получим передаточную функцию фильтра верхних
54
частот. С точки зрения преобразования (73) это означает, что если аргумент p
аналогового фильтра нижних частот заменить на 1 z 1 / 1 z 1 , то получим
выражение для передаточной функции цифрового фильтра верхних частот. Итак, в общем случае преобразование аналогового фильтра нижних
частот в цифровой фильтр верхних частот описывается формулой
p k |
1 z 1 |
, |
(79) |
|||
|
|
|
||||
1 |
z 1 |
|||||
|
|
|
||||
где k - действительная положительная константа, значение которой таково, что частоты среза аналогового фильтра нижних частот 1 и цифрового фильтра верхних частот 1 связаны соотношением
|
T |
|
|
|
|
|
|
|||
k 1tg |
|
1 |
|
. |
|
|
|
(80) |
||
|
|
|
2 |
|
|
|
|
|
|
|
При этом преобразование принимает вид |
|
|
|
|
|
|
|
|
||
|
T |
1 z 1 |
|
|||||||
p 1tg |
1 |
|
|
|
|
|
|
|
. |
(81) |
2 |
|
|
z |
1 |
||||||
|
|
1 |
|
|
|
|
||||
Если аналоговый фильтр нижних частот нормирован, т.е. если |
1 1 , то |
|||||||||
T |
|
1 z 1 |
|
|||||||
p tg |
1 |
|
|
|
|
|
. |
(82) |
||
|
1 z |
1 |
||||||||
|
2 |
|
|
|
|
|
|
|||
В каждом из двух рассмотренных случаев преобразования имеют форму, исключающую необходимость:
1. предварительной деформации (эта операция включена в действительную константу k );
2. преобразование передаточной функции аналогового фильтра нижних частот в передаточную функцию аналогового фильтра верхних частот.
Таким образом, преимущества этого преобразования с точки зрения минимизации вычислений очевидны.
1.3.2 Примеры расчета цифровых БИХ-фильтров с внесением предыскажений
1.3.2.1 Пример расчёта домашнего задания для ФВЧ Чебышёва и Баттерворта. Требования к фильтрам Чебышёва и Баттерворта остались такими же, как
впредыдущей лабораторной работе, - для ФНЧ это:
1.значения граничных частот полосы задержания и полосы пропускания ( Fз и F1 соответственно);
2.значение коэффициента максимального ослабления в полосе пропускания aр max , дБ;
3.значение коэффициента минимального ослабления в полосе задержания aр min , дБ;
55
4. частота дискретизации сигналов, обрабатываемых фильтром, Fд . Примем значения граничных частот и коэффициентов aр max , aр min
низкочастотного аналогового фильтра прототипа равными:
|
F1 1900Гц, Fз 4940Гц |
|||
|
ap max 0.4455дБ, ap min |
40дБ |
||
|
|
Fд 17000 Гц |
|
|
Значения |
граничных |
частот |
и |
коэффициентов aр max , aр min |
высокочастотного аналогового фильтра равны:
F1 4940Гц, Fз 1900Гц
ap max 0.4455дБ, ap min 40дБ
Fд 17000 Гц
Синтез фильтров верхних частот производится на основе рассчитанных в предыдущей лабораторной работе фильтров нижних частот Чебышёва и Баттерворта, путём проведения преобразования Константинидиса. Благодаря
этому преобразованию возможен переход от аналогового ФНЧ сразу к цифровому ФВЧ с внесением предыскажений.
Преобразование Константинидиса состоит в подстановке в операторную передаточную функцию выражения:
|
T |
1 |
z 1 |
|
|
|||
p tg |
1 |
|
|
|
|
|
, |
(83) |
2 |
|
|
z |
1 |
||||
|
1 |
|
|
|
||||
где 1 - граничная частота полосы пропускания ФВЧ в рад/с, T - интервал дискретизации в с. В рассматриваемом случае
|
T 1 z 1 |
1 z 1 |
|
|||||||
p tg |
1 |
|
|
|
|
1.2941 |
|
|
|
(84) |
|
|
z |
1 |
|
z |
1 |
||||
|
2 |
1 |
|
1 |
|
|
||||
Проведём данную подстановку сначала в операторную передаточную функцию ФНЧ Чебышёва ( H ( p) взята с домашнего расчёта лабораторной работы № 3 «Синтез цифровых БИХ – фильтров методом билинейного z -
преобразования»). |
|
|
|
|
0.38032 |
|
|
|
|
Н ( p) |
|
|
|
(85) |
( p2 0.87608 p 0.37124)( p2 |
0.36288 p 1.0784) |
|||
После подстановки в Н ( p) выражения |
(84) и проведения |
простых |
||
математических операций получаем передаточную (системную) функцию цифрового ФВЧ Чебышёва:
H (z) |
0.0371(1 z 1)4 |
(86) |
(1 0.8199z 1 0.2869z 2 )(1 0.3701z 1 0.7086z 2 ) |
||
|
56 |
|
Разобьём операторную функцию на множители не выше 2 порядка, т.е. получим Н (z) в форме биквадратных звеньев и выпишем отдельно выражения
для каждого звена: |
|
|
|
|
|
|
||||
|
|
|
|
0.0371(1 z 1)2 |
0.0371(1 2z 1 z 2 ) |
|
||||
H1(z) |
|
|
|
|
|
|
|
|
(87) |
|
|
|
0.8199z 1 |
|
|
(1 0.8199z 1 0.2869z 2 ) |
|||||
|
(1 |
0.2869z 2 ) |
|
|||||||
|
|
|
|
(1 z 1)2 |
(1 2z 1 z 2 ) |
|
||||
H2 |
(z) |
|
|
|
|
|
|
|
(88) |
|
(1 0.3701z 1 |
|
|
|
(1 0.3701z 1 0.7086z 2 ) |
||||||
|
|
|
0.7086z 2 ) |
|
||||||
Тем самым мы привели наше выражение к виду Н (z) H1 (z)H2 (z) и получили два звена второго порядка.
Выпишем отдельно коэффициенты каждого звена в таблицу, без учета общего коэффициента усиления фильтра.
Таблица 5. Коэффициенты звеньев ФВЧ Чебышева
Звено № |
|
b0 |
|
b1 |
|
b2 |
|
a0 |
a1 |
a2 |
||
1 |
1 |
|
-2 |
|
1 |
|
1 |
0.8199 |
0. 2869 |
|
||
2 |
1 |
|
-2 |
|
1 |
|
1 |
0.3701 |
0.7086 |
|
||
Аналогичный расчёт проводится и для фильтра Баттерворта, операторная |
||||||||||||
передаточная функция, которого равна: |
|
|
|
|
|
|
||||||
Н (p) |
|
|
|
|
3.0893 |
|
|
|
|
|||
|
|
|
||||||||||
(p2 2.3314p 1.4564)(p2 1.7067p 1.4564)(p2 +0.6247p 1.4564) |
||||||||||||
|
|
|
|
|
|
(89) |
|
|
|
|
|
|
Передаточная (системная) функция ФВЧ Баттерворта после |
||||||||||||
преобразования Константинидиса: |
|
|
|
|
|
|
|
|||||
H (z) |
|
|
|
|
|
0.0239(1 z 1 )6 |
|
|
|
|
||
(1 0.071z 1 |
0.0186z 2 )(1 0.0818z 1 |
0.1728z 2 )(1 0.1109z 1 0.5896z 2 ) |
||||||||||
Запишем Н (z) в форме биквадратных звеньев (в данном случае биквадратных звена, так как порядок фильтра шестой):
Н (z) H1 (z)H2 (z)H3 (z)
|
|
|
|
|
|
0.0239(1 z 1)2 |
||||
H1 |
(z) |
|
|
|
|
|
|
|||
(1 0.071z 1 |
|
|
||||||||
|
|
|
|
|
0.0186z 2 ) |
|||||
|
|
|
|
|
|
(1 z 1)2 |
||||
H2 |
(z) |
|
|
|
|
|
|
|||
|
0.0818z 1 |
|
|
|
||||||
|
|
(1 |
0.1728z 2 ) |
|||||||
|
|
|
|
|
|
(1 z 1)2 |
||||
H3 |
(z) |
|
|
|
|
|
|
|||
|
0.1109z 1 |
|
|
|
||||||
|
|
(1 |
0.5896z 2 ) |
|||||||
|
|
|
|
|
|
|
|
57 |
|
|
0.0239(1 2z 1 z 2 )
(1 0.071z 1 0.0186z 2 )
(1 2z 1 z 2 )
(1 0.0818z 1 0.1728z 2 ) (1 2z 1 z 2 )
(1 0.1109z 1 0.5896z 2 )
(90)
три
(91)
(92)
(93)
Таблица 6. Коэффициенты звеньев ФВЧ Баттерворта
Звено № |
b0 |
b1 |
b2 |
a0 |
a1 |
a2 |
1 |
1 |
-2 |
1 |
1 |
0.071 |
0.0186 |
2 |
1 |
-2 |
1 |
1 |
0.0818 |
0.1728 |
3 |
1 |
-2 |
1 |
1 |
0.1109 |
0.5896 |
На этом домашний расчёт БИХ-фильтров верхних частот Чебышёва и Баттерворта окончен.
1.3.2.2 Пример расчёта домашнего задания для ФВЧ Чебышёва и Баттерворта. Требования к фильтрам Чебышёва и Баттерворта остались такими же, как
впредыдущей лабораторной работе, - для ФНЧ это:
1.значения граничных частот полосы задержания и полосы пропускания ( Fз и F1 соответственно);
2.значение коэффициента максимального ослабления в полосе пропускания aр max , дБ;
3.значение коэффициента минимального ослабления в полосе задержания aр min , дБ;
4.частота дискретизации сигналов, обрабатываемых фильтром, Fд .
Примем значения граничных частот и коэффициентов aр max , aр min низкочастотного аналогового фильтра прототипа равными:
|
F1 1900Гц, Fз 4940Гц |
|||
|
ap max 0.4455дБ, ap min |
40дБ |
||
|
|
Fд 17000 Гц |
|
|
Значения |
граничных |
частот |
и |
коэффициентов aр max , aр min |
высокочастотного аналогового фильтра равны:
F1 4940Гц, Fз 1900Гц
ap max 0.4455дБ, ap min 40дБ
Fд 17000 Гц
Синтез фильтров верхних частот производится на основе рассчитанных в предыдущей лабораторной работе фильтров нижних частот Чебышёва и Баттерворта, путём проведения преобразования Константинидиса. Благодаря
этому преобразованию возможен переход от аналогового ФНЧ сразу к цифровому ФВЧ с внесением предыскажений.
Преобразование Константинидиса состоит в подстановке в операторную передаточную функцию выражения:
|
T |
1 |
z 1 |
|
|
||
p tg |
1 |
|
|
|
|
, |
(94) |
2 |
|
z |
1 |
||||
|
1 |
|
|
|
|||
|
58 |
|
|
|
|
|
|
где 1 - граничная частота полосы пропускания ФВЧ в рад/с, T - интервал дискретизации в с. В рассматриваемом случае
|
T |
1 z 1 |
1 z 1 |
|
||||||
p tg |
1 |
|
|
|
|
1.2941 |
|
|
|
(95) |
|
|
1 z |
1 |
|
z |
1 |
||||
|
2 |
|
|
|
1 |
|
|
|||
Проведём данную подстановку сначала в операторную передаточную |
||||||||||
функцию ФНЧ Чебышёва ( H ( p) |
|
взята с домашнего расчёта лабораторной |
||||||||
работы № 3 «Синтез цифровых БИХ – фильтров методом билинейного z - преобразования»).
Н ( p) |
0.38032 |
|
|
(96) |
|
|
|
||
( p2 0.87608 p 0.37124)( p2 |
0.36288 p 1.0784) |
|||
После подстановки в Н ( p) выражения |
(84) и проведения |
простых |
||
математических операций получаем передаточную (системную) функцию цифрового ФВЧ Чебышёва:
H (z) |
0.0371(1 z 1 )4 |
(97) |
(1 0.8199z 1 0.2869z 2 )(1 0.3701z 1 0.7086z 2 ) |
Разобьём операторную функцию на множители не выше 2 порядка, т.е. получим Н (z) в форме биквадратных звеньев и выпишем отдельно выражения
для каждого звена: |
|
|
|
|
|
|
||||
|
|
|
|
0.0371(1 z 1)2 |
0.0371(1 2z 1 z 2 ) |
|
||||
H1(z) |
|
|
|
|
|
|
|
|
(98) |
|
|
|
0.8199z 1 |
|
|
(1 0.8199z 1 0.2869z 2 ) |
|||||
|
(1 |
0.2869z 2 ) |
|
|||||||
|
|
|
|
(1 z 1)2 |
(1 2z 1 z 2 ) |
|
||||
H2 |
(z) |
|
|
|
|
|
|
|
(99) |
|
(1 0.3701z 1 |
|
|
|
(1 0.3701z 1 0.7086z 2 ) |
||||||
|
|
|
0.7086z 2 ) |
|
||||||
Тем самым мы привели наше выражение к виду Н (z) H1 (z)H2 (z) и получили два звена второго порядка.
Выпишем отдельно коэффициенты каждого звена в таблицу, без учета общего коэффициента усиления фильтра.
Таблица 7. Коэффициенты звеньев ФВЧ Чебышева
Звено № |
|
b0 |
b1 |
b2 |
|
a0 |
a1 |
a2 |
|
1 |
|
1 |
-2 |
1 |
|
1 |
0.8199 |
0. 2869 |
|
2 |
|
1 |
-2 |
1 |
|
1 |
0.3701 |
0.7086 |
|
Аналогичный расчёт проводится и для фильтра Баттерворта, операторная |
|||||||||
передаточная функция, которого равна: |
|
|
|
|
|||||
Н (p) |
|
|
3.0893 |
|
|
|
|
||
|
|
||||||||
(p2 2.3314p 1.4564)(p2 1.7067p 1.4564)(p2 +0.6247p 1.4564) |
|||||||||
|
|
|
|
(100) |
|
|
|
|
|
|
|
|
|
59 |
|
|
|
|
|
Передаточная (системная) функция ФВЧ Баттерворта после преобразования Константинидиса:
0.0239(1 z 1 )6
H (z) (1 0.071z 1 0.0186z 2 )(1 0.0818z 1 0.1728z 2 )(1 0.1109z 1 0.5896z 2 )
(101)
Запишем Н (z) в форме биквадратных звеньев (в данном случае три
биквадратных звена, так как порядок фильтра шестой):
Н (z) H1 (z)H2 (z)H3 (z)
|
|
|
|
|
|
0.0239(1 z 1)2 |
||||
H1 |
(z) |
|
|
|
|
|
|
|||
(1 0.071z 1 |
|
|
||||||||
|
|
|
|
|
0.0186z 2 ) |
|||||
|
|
|
|
|
|
(1 z 1)2 |
||||
H2 |
(z) |
|
|
|
|
|
|
|||
|
0.0818z 1 |
|
|
|
||||||
|
|
(1 |
0.1728z 2 ) |
|||||||
|
|
|
|
|
|
(1 z 1)2 |
||||
H3 |
(z) |
|
|
|
|
|
|
|||
|
0.1109z 1 |
|
|
|
||||||
|
|
(1 |
0.5896z 2 ) |
|||||||
0.0239(1 2z 1 z 2 )
(1 0.071z 1 0.0186z 2 )
(1 2z 1 z 2 )
(1 0.0818z 1 0.1728z 2 ) (1 2z 1 z 2 )
(1 0.1109z 1 0.5896z 2 )
(102)
(103)
(104)
Таблица 8. Коэффициенты звеньев ФВЧ Баттерворта
Звено № |
b0 |
b1 |
b2 |
a0 |
a1 |
a2 |
1 |
1 |
-2 |
1 |
1 |
0.071 |
0.0186 |
2 |
1 |
-2 |
1 |
1 |
0.0818 |
0.1728 |
3 |
1 |
-2 |
1 |
1 |
0.1109 |
0.5896 |
На этом домашний расчёт БИХ-фильтров верхних частот Чебышёва и Баттерворта окончен.
1.3.3Контрольные вопросы
1.Почему возникла необходимость вводить предыскажения при использовании метода билинейного z - преобразования?
2.Как кроме внесения предыскажений можно получить цифровой фильтр с заданной граничной частотой полосы пропускания 1 ? Чем это неудобно?
3.В чём заключается преобразование Константинидиса?
4.Поясните преобразование Константинидиса на примере фильтра нижних частот (ФНЧ).
5.Как выглядит преобразование Константинидиса для фильтра верхних частот (ФВЧ).
60
