LR_TsOS_2_BIKh
.pdfФедеральное агентство связи ордена Трудового Красного Знамени
федеральное государственное бюджетное образовательное учреждение высшего образования
«Московский технический университет связи и информатики»
Кафедра радиотехнических систем
Практикум по дисциплине
ЦИФРОВАЯ ОБРАБОТКА СИГНАЛОВ
Лабораторная работа № 3
СИНТЕЗ ЦИФРОВЫХ БИХ-ФИЛЬТРОВ МЕТОДОМ БИЛИНЕЙНОГО Z -ПРЕОБРАЗОВАНИЯ
Москва, 2016 г.
УДК 621.391:519.27 |
План подготовки УМД 2015/2016 уч. года |
Практикум по дисциплине
«ЦИФРОВАЯ ОБРАБОТКА СИГНАЛОВ»
СИНТЕЗ ЦИФРОВЫХ БИХ-ФИЛЬТРОВ МЕТОДОМ БИЛИНЕЙНОГО Z - ПРЕОБРАЗОВАНИЯ
В лабораторной работе № 3 изучается метод билинейного z - преобразования на примере синтеза цифровых ФНЧ и ФВЧ по аналоговым фильтрам прототипам Чебышёва и Баттерворта. Провести моделирование синтезированного цифрового фильтра в среде «Спектр-2».
Для студентов радиотехнических и телекоммуникационных специальностей.
Список лит. 11 назв., ил. 9, табл. 5.
Составители: Лобов Е.М., Смердова Е.О. Рецензент
Издание утверждено советом факультета Радио и Телевидения. Протокол № … от
_._.2016 г.
2
ЛАБОРАТОРНАЯ РАБОТА № 3
СИНТЕЗ ЦИФРОВЫХ БИХ-ФИЛЬТРОВ МЕТОДОМ БИЛИНЕЙНОГО Z -ПРЕОБРАЗОВАНИЯ
Цель работы:
Изучение метода билинейного z - преобразования на примере синтеза цифровых ФНЧ и ФВЧ по аналоговым фильтрам прототипам Чебышёва и Баттерворта. Моделирование синтезированного цифрового фильтра проводится в среде «Спектр-2».
Номер варианта выбирается студентом из следующей таблицы:
Условные |
|
Список требований |
Задаваемые значения |
||||||||||||
обозначения |
|
||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
F |
|
|
|
Частота |
F |
10000 200N |
|
|
|
|
||||
|
|
|
|
|
бр |
|
|
|
|||||||
|
|
д |
|
|
|
дискретизации |
д |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
F 1200 20N |
бр |
– Баттерворт |
||||||
|
|
1 |
|
|
|
Граничная частота ПП |
1 |
|
|
|
|
|
|
|
|
|
F |
|
|
|
|
|
|
|
|
– Чебышев |
|||||
|
|
|
|
|
|
|
F 1500 30N |
бр |
|||||||
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
N |
бр |
|
|
|
|
|
|
|
|
|
Fз |
3120 50Nбр |
10 |
|
– Баттерворт |
||||
|
|
|
|
|
|
|
5 |
||||||||
|
F |
|
|
|
Граничная частота ПЗ |
|
|
|
|
|
|
|
|||
|
|
з |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Fз |
2800 40Nбр |
|
– Чебышев |
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Максимально |
|
|
|
|
|
|
|
|
|
|
a |
max |
|
|
|
допустимое затухание |
a |
|
0.4455 |
|
|
|
|
|
|
|
|
|
|
|
max |
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
в ПП (дБ) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Минимально |
|
|
|
|
|
|
|
|
|
|
amin |
|
|
|
допустимое затухание |
amin 40 |
|
|
|
|
|
|
|||
|
|
|
|
|
|
в ПЗ (дБ) |
|
|
|
|
|
|
|
|
|
Nb |
1, 2,...40 |
, |
Nb |
- результат целочисленного деления номера бригады на 5 (с |
|||||||||||
|
|||||||||||||||
|
|
|
5 |
|
|
|
|
|
|
|
|
|
|
отбрасыванием цифр после запятой у результата деления).
Тип функции фильтрации определяется по чётности/нечётности последней цифры в номере студенческого билета. Если последняя цифра номера нечётная, то требуется провести синтез фильтра Чебышёва, если чётная, то рассчитывается фильтр Баттерворта.
По предпоследней цифре в номере студенческого билета определяется тип частотной избирательности фильтра. Если предпоследняя цифра номера
3
нечётная, то рассчитывается фильтр нижних частот, чётная – фильтр верхних частот.
1 Домашний расчёт.
1.1 Домашний расчёт состоит из следующих пунктов:
1.Формулировка требований к ЦФ.
2.Синтез передаточной функции АФП, включающий в себя: 2.1.оценку порядка ПФ АФП;
2.2.расчет полюсов и нулей, а также коэффициента усиления передаточной функции АФП;
2.3.запись аналитического выражения передаточной функции АФП.
3.Переход от передаточной функции АФП к передаточной функции ЦФ с помощью билинейного z -преобразования.
4.Запись аналитического выражения передаточной функции ЦФ в виде произведения передаточных функций биквадратных звеньев.
1.2 Пример расчёта домашнего задания для ФНЧ Чебышёва.
Спроектируем ФНЧ Чебышева I рода. Сформулируем требования к фильтру, - для ФНЧ это:
1.значения граничных частот полосы задержания и полосы пропускания ( Fз и F1 соответственно);
2.значение коэффициента максимального ослабления в полосе пропускания aр max , дБ;
3.значение коэффициента минимального ослабления в полосе
задержания aр min |
, дБ; |
4. частота дискретизации сигналов, обрабатываемых фильтром,
Fд
.
Примем значения граничных частот и коэффициентов
a |
р max |
|
,
a |
р min |
|
низкочастотного аналогового фильтра прототипа равными:
F 1900Гц, F 4940Гц, |
||||
|
1 |
з |
|
|
a |
p max |
0.4455дБ, a |
p min |
40дБ, |
|
|
|
Fд 17000 Гц.
Синтез АФП осуществляется в области нормированных частот. При этом граничная угловая частота полосы пропускания принимается равной 1 рад/с. Проведём соответствующее нормирование оси частот. Тогда
4
нормированные граничные частоты полосы пропускания и полосы задержания будут равны соответственно
|
2 F |
1, |
|
|
2 F |
2.6. |
|
1 |
|
з |
|
||||
|
|
|
|
|
|
||
1 |
2 F |
|
з |
|
2 F |
|
|
|
|
|
|
|
|
||
|
1 |
|
|
|
1 |
|
|
Следующим шагом синтеза является определение порядка фильтра и |
|||||||
коэффициента неравномерности ослабления в |
полосе пропускания |
. |
Напомним, что данные параметры необходимы для нахождения корней знаменателя модуля квадрата амплитудно-частотной характеристики (АЧХ) которая имеет вид:
Н (j ) |
2 |
|
|
|
|
|
1 |
1 |
2 |
ch |
2 |
nArch |
|
|
,
(1)
где 2 F модуля квадрата, ch(x) - гиперболический косинус, Arch(x) - гиперболический арккосинус.
С учетом (1) квадрат модуля передаточной функции АФП можно записать в форме:
Н ( p) |
|
2 |
H ( p)H ( p) |
2 |
|
2 |
1 |
|
|
||||||
|
|
|
nArch( jp) |
||||
|
|
|
1 |
ch |
|||
|
|
|
|
|
,
(2)
где |
H ( p) |
- операторная передаточная функция АФП, |
|
T ( ) ch n Arch |
- полином Чебышёва, являющийся |
p j .
функцией
фильтрации данного фильтра.
Коэффициент неравномерности рабочего ослабления найдём через коэффициент максимального подавления и значение граничной частоты полосы пропускания:
A |
|
|
|
|
|
2 |
2 |
( )) a |
|
, |
||
10lg(1 |
T |
|
p max |
|||||||||
p |
|
|
|
|
|
|
|
|
|
1 |
|
|
T ( ) T(1) 1 |
|
|
|
|||||||||
|
1 |
|
|
|
|
|
|
|
|
|
|
|
10lg(1 |
2 |
) a |
|
|
|
, |
|
|
||||
|
|
p max |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|||
|
10 |
0.1a |
p max |
1. |
|
|
|
|
|
|||
Подставив численное значение |
|
aр max , получаем 0.32867 . |
Порядок фильтра можно определить через значения граничной частоты полосы задержания и коэффициента минимального ослабления в полосе задержания.
5
|
|
|
10 |
0.1a |
p min |
1 |
|
|
|
|
|
|
|
||||
nArch |
3 |
Arch |
|
|
2 |
|
, |
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
10 |
0.1ap min |
1 |
|
|
|||
|
Arch |
|
|
|
|
|
|
||
|
|
|
2 |
|
|
||||
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||
n |
|
|
|
|
|
|
|
|
. |
Arch |
|
|
|
|
|
||||
|
3 |
|
|
|
|||||
|
|
|
|
|
|
|
|
Порядок фильтра может быть только целым числом,
A |
2 |
2 |
( |
)) a |
|
10lg(1 T |
|
p min |
|||
p |
|
|
з |
|
тогда справедливо:
,
|
2ch2 |
|
nArch |
3 |
|
|
|
10lg 1 |
|
|
|
|
a |
p min |
|
,
|
|
|
|
|
10 |
0.1a |
p min |
1 |
||||||
ch nArch |
|
|
|
|
|
|
|
|
|
|||||
3 |
|
|
|
|
|
2 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
10 |
0.1ap min |
|
1 |
|
|
|||||||
|
Arch |
|
|
|
|
|
|
|
|
|||||
|
|
|
2 |
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
Arch |
|
|
|
|
|
|
|
||||
|
|
|
3 |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
,
(3)
где x – операция округления в сторону ближайшего целого значения,
превышающего x .
Проведя арифметические операции с учетом заданных требований к фильтру, получим n 4 . Теперь мы имеем все данные для нахождения корней знаменателя квадрата модуля АЧХ.
Корни находятся по формуле:
pk
где
|
1 |
|
|
|
|
|
|
1 |
|
|
|
|
|
|
||||
|
Arsh |
|
|
2k 1 |
|
|
Arsh |
|
|
2k 1 |
|
|
||||||
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
|
||||||||||
sh |
|
sin |
|
|
|
jch |
|
cos |
|
|
|
, k 1,2,...,2n. |
||||||
|
|
|
|
|
|
|
|
|||||||||||
|
n |
|
|
|
2n |
|
|
n |
|
|
|
2n |
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
sh(x) |
- гиперболический синус. |
|
|
|
|
|
|
|
|
|
Из рассчитанных корней выбираем отвечающие соображениям устойчивости (система устойчива, если полюсы её передаточной функции лежат в левой полуплоскости её p -плоскости, т.е. имеют отрицательную вещественную часть).
С учетом изложенного, отобранные корни равны:
6
p -0.43804-j0.42352, |
|
1 |
|
p |
-0.18144-j1.0225, |
2 |
|
p |
-0.18144+j1.0225, |
3 |
|
p |
-0.43804+j0.42352. |
4 |
|
(4)
В результате расчетов получены 2 пары комплексно-сопряжённых корней, которые мы относим к H ( p) . Последним шагом в формировании низкочастотного АФП Чебышёва является запись аналитического выражения его операторной передаточной функции, которая имеет вид:
Н ( p) |
1 |
|
|
n 1 |
|
||
2 |
|
||
|
где |
|
1 |
|
– коэффициент усиления. |
|
2 |
n 1 |
|
|||
|
|
||||
|
|
|
В нашем случае H ( p) принимает вид:
1
n ( p
k 1
,
p |
k |
) |
|
|
Н ( p) |
0.38032 |
|
|
|
|
|
|
( p p )( p p )( p p )( p p ) |
|||
1 |
2 |
3 |
4 |
где |
1 |
0.38032 . |
3 |
||
|
2 |
|
,
(5)
Дальнейшие действия синтеза аналоговых фильтров рассматриваться не будут, так как в них нет необходимости при получении передаточной функции цифрового БИХ-фильтра.
Для перехода от передаточной функции АФП к передаточной функции ЦФ используется билинейное z - преобразование, которое заключается в подстановке в передаточную функцию выражения
где
Fд
|
1 z |
1 |
|
p 2F |
|
||
|
|
1 |
|
д |
1 |
z |
|
|
|
– частота дискретизации сигналов,
,
обрабатываемых фильтром,
(6)
z |
– |
комплексный параметр Z-преобразования. Перед использованием билинейного
передаточную функцию Н ( p) (5) в удобном
Z-преобразования, преобразуем для этого виде:
Н ( p) |
|
|
0.38032 |
|
|
|
|
|
|||
|
|
|
|||||||||
( p2 p( p p ) p p )( p2 |
p( p p ) p p ) |
||||||||||
|
|
1 |
4 |
1 |
4 |
2 |
3 |
2 |
3 |
(7) |
|
|
|
|
|
0.38032 |
|
|
|
|
|
||
|
|
|
|
|
|
|
|
. |
|
||
|
|
||||||||||
( p2 0.87608 p 0.37124)( p2 0.36288 p 1.0784) |
|
В (7) проведено перемножение множителей, содержащих комплексносопряжённые корни полинома знаменателя ПФ (5) (см. значение корней (4)).
7
Теперь можно проводить билинейное |
z |
- |
предварительно нормировав частоту дискретизации (примем
|
|
F |
|
|
|
Fs |
|
д |
1.424 |
, |
|
2 F |
|||||
|
|
|
|
||
|
|
1 |
|
|
тогда замена примет вид
преобразование,
Fд |
17000Гц ) |
p 2F 1 z
s 1 z
В результате подстановки в |
Н |
1 |
1 |
z 1 |
|
|
2.848 |
|
|
1 |
1 |
z 1 |
|
( p) выражения |
, |
(8) |
(8) получаем:
Н (z) |
|
|
|
|
|
|
|
|
|
|
|
|
|
0.38032 |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
1 |
|
2 |
|
|
|
|
|
|
|
1 |
|
|
|
|
||||||
|
|
|
|
1 z |
|
|
|
|
|
1 z |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
2.848 |
|
|
|
|
0.87608 |
2.848 |
|
|
|
0.37124 |
|
||||||||||||||
|
|
|
|
1 |
|
|
|
|
|
|
1 |
|
|
||||||||||||||
|
|
|
|
1 z |
|
|
|
|
|
1 z |
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
. |
|
|
|
|
|
|
|
|
|
1 |
|
|
2 |
|
|
|
|
|
|
|
|
1 |
|
|
|
|
||||
|
|
|
|
|
1 z |
|
|
|
|
|
|
1 z |
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
2.848 |
|
|
0.36288 |
2.848 |
|
1.0784 |
|
|
||||||||||||||||
|
|
|
|
|
1 |
|
|
|
|
|
1 |
|
|
|
|||||||||||||
|
|
|
|
|
1 z |
|
|
|
|
|
1 z |
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Умножим числитель и знаменатель каждой дроби на
|
|
1 z 1 |
2 |
|
и получим
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.38032 1 z 1 |
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||
Н (z) |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||
|
|
|
|
2.848 1 z 1 |
|
|
|
0.87608 2.848 1 |
z 1 |
|
|
|
1 z 1 |
|
|
0.37124 1 z 1 |
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 z 1 |
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
. |
|
||||||
|
|
|
|
|
|
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|||||||||||||
|
|
|
|
|
2.848 1 z 1 |
|
|
0.36288 2.848 |
1 z 1 |
|
|
1 z 1 |
|
1.0784 |
1 z 1 |
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||||||||||
Раскроем скобки |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2z 1 z 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
Н (z) |
|
|
|
|
|
|
|
|
|
|
|
0.38032 1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
z 2 |
|
|
|
|
z 2 |
|
|
|
|
|
|
|
|
2z 1 |
z 2 |
|
|
|
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||||||||||
|
|
|
|
|
|
|
8.1111 1 2z 1 |
|
|
2.4951 1 |
|
|
0.37124 1 |
|
|
|
|
|
|
||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2z 1 |
z 2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.0335 |
|
|
|
|
|
|
|
|
|
z 2 |
|
|
|
|
|
||||||||||||||
|
|
|
|
|
|
|
|
8.1111 1 2z 1 z 2 |
|
1 z 2 |
|
|
|
1.0784 1 2z 1 |
|
|
|
|
|
|
|
||||||||||||||||||||||||
и сложим слагаемые с одинаковыми степенями при z |
|
|
|
|
|
|
|
|
|
|
|
|
|
||||||||||||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.38032 1 |
|
2z 1 z 2 |
|
|
|
|
|
|
|
|
|
|
|
|
1 2z 1 |
z 2 |
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
15.4797z 1 5.9872z 2 |
|
|
|
|
|
|
|
|
|
|
|
|
8.1560z 2 |
|
|
||||||||||||||||||||||
Н (z) |
10.9774 |
|
|
|
10.2230 14.0654z 1 |
|
. |
Поделим числитель и знаменатель первой дроби на 10.9774 , а второй дроби - на 10.2230 , для того чтобы свободные члены в знаменателе стали равны единице, и получим
8
|
|
|
2z 1 |
z 2 |
|
|
|
|
2z 1 z 2 |
|
|
|
|||
Н (z) |
0.03464 1 |
|
0.09781 1 |
|
|
|
. |
||||||||
|
1.4101z 1 |
0.5454z 2 |
|
1.3759z 1 |
|
0.7978z 2 |
|
||||||||
|
|
||||||||||||||
|
1.0 |
1.0 |
|
|
Тем самым |
мы |
|
Н (z) 0.0033888 * H |
(z)H |
(z) |
1 |
2 |
|
привели |
наше |
выражение |
к |
виду |
и получили |
два |
звена второго |
порядка, где |
0.0033888 0.03464 0.09781 – общий коэффициент усиления. Выпишем отдельно выражения для каждого звена:
|
|
|
|
1 2z |
1 |
z |
2 |
|
|
b |
b z |
1 |
b z |
2 |
||||||
H |
(z) |
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
0 |
1 |
|
|
|
2 |
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
1 |
1 1.4101z |
1 |
0.5454z |
2 |
|
1 a z |
1 |
a z |
2 |
||||||||||
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|
|
|
|
|
1 2z 1 |
z 2 |
|
|
b b z 1 |
b z 2 |
|||||||||||
H |
|
(z) |
|
|
|
|
|
|
|
|
|
0 |
1 |
|
|
|
2 |
|
|
|
2 |
|
1.3759z 1 |
0.7978z 2 |
1 a z 1 |
|
|
z 2 |
|||||||||||||
|
1 |
|
a |
|||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
2 |
|
|
|
(9)
(10)
Выпишем отдельно коэффициенты каждого звена в таблицу, без учета общего коэффициента усиления фильтра.
Таблица 1. Коэффициенты звеньев ФНЧ Чебышева
Звено № |
b |
b |
b |
a |
a |
a |
|
0 |
1 |
2 |
0 |
1 |
2 |
1 |
1 |
2 |
1 |
1 |
-1.4101 |
0.54542 |
2 |
1 |
2 |
1 |
1 |
-1.3759 |
0.79781 |
Определим искажение частот, получившееся в
билинейного z |
-преобразования. Частоты АЧХ АФП |
соотношениями: |
|
ходе проведения и ЦФ связаны
где
ˆ
/
Fд
ˆ 2arctg |
|
|
2F |
||
|
||
|
д |
F |
F |
|
д |
||
|
||
|
|
- нормированная частота
|
2arctg |
F |
||
F |
||||
|
|
|
||
|
|
|
д |
|
tg |
ˆ |
, |
|
|
2 |
|
|||
|
|
|
АЧХ цифрового фильтра,
2
(11)
(12) f –
угловая частота АЧХ цифрового фильтра.
Подставляя в выражение для нахождения частоты АЧХ цифрового фильтра, значения граничных частот полосы пропускания и полосы задержания аналогового фильтра, получим граничные частоты ЦФ:
f |
|
|
|
Fд |
|
arctg |
F1 |
1827.2311 Гц, |
|
|
|
|
|
||||||
1 |
|
|
|
F |
|
|
|||
|
|
|
|
|
|
|
д |
|
(13) |
|
|
|
|
Fд |
arctg Fз |
||||
f |
|
|
4003.8081 Гц. |
||||||
з |
|
||||||||
|
|
|
|
Fд |
|
|
|||
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
9 |
|
На этом домашний расчёт БИХ-фильтра Чебышёва окончен.
1.3 Пример расчёта домашнего задания для ФНЧ Баттерворта.
Требования к фильтру и, соответственно, значения нормированных частот такие же, как и для ФНЧ Чебышёва (см. пункт 1.1 домашнего расчёта)
Квадрат модуля АЧХ фильтра Баттерворта имеет вид:
Н (j ) |
2 |
|
1 |
, |
|
|
|||
|
1 2 2n |
где ( ) 2n - функция фильтрации фильтра Баттерворта. Квадрат модуля передаточной функции АФП можно записать виде:
где
H (
p)
Н ( p) |
2 |
H ( p)H ( p) |
|
1 |
|
2 |
( jp |
||
|
|
|
||
|
|
|
1 |
- операторная передаточная функция АФП,
) |
2 |
|
|
p |
n |
, |
|
|
|
j
.
Коэффициент неравномерности рабочего ослабления в полосе пропускания и порядок фильтра для фильтра Баттерворта определяется по формулам:
|
|
|
|
|
10 |
0.1a |
p max |
1 |
, |
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
10 |
0.1ap min |
1 |
|
|
|
|
|
10 |
0.1a |
p min |
1 |
|
||||
|
lg |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
2 |
|
|
lg |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2 |
|
|
|
|
n |
|
|
|
|
|
|
|
|
|
|
|
|
, |
||||||
|
lg |
|
|
|
|
|
|
|
2lg |
|
|
||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
3 |
|
|
|
|
|
|
3 |
|
||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
где x – операция округления в сторону ближайшего целого значения,
превышающего x .
Проведя несложный арифметический расчёт, учитывающий выдвинутые к фильтру требования, получим значения 0.32867 , n 6 .
Все данные для нахождения корней знаменателя квадрата модуля АЧХ получены.
Корни найдём по формуле для чётных значений порядка фильтра Баттерворта:
p |
|
|
|
1 |
|
cos |
|
2k 1 |
|
|
j sin |
|
k |
n |
|
|
|
2n |
|
|
|||||
|
|
|
|
|
|
|||||||
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
Из рассчитанных корней, выбираем те, устойчивости:
2k 1 |
|
|
|
2n |
|
|
которые
,
k 1, 2,..., 2 n .
отвечают соображениям
10