Министерство цифрового развития, связи и массовых коммуникаций Российской Федерации
Ордена Трудового Красного Знамени федеральное государственное бюджетное образовательное учреждение высшего образования
Московский технический университет связи и информатики
___________________________________________________________________
Кафедра «Теории электрических цепей»
Лабораторная работа №15 по дисциплине ОКАЭЦ «Исследование БИХ-фильтров»
Выполнил |
|
|
Студент группы БИК2205 |
_________________________ |
|
Проверил |
|
|
Доцент кафедры ТЭЦ |
_________________________ |
Бакулин М.Г. |
Москва 2023
1.Цель работы
Спомощью программы Micro-Cap получить основные временные и частотные характеристики фильтров с бесконечной импульсной характеристикой (БИХ-фильтров).
2.Предварительный расчёт
Исходные данные:
FT 0,10-4 1 – используемый в опыте диапазон;
T – интервал дискретизации по времени;
j
-1 – мнимая единица;
3.Расчётные формулы
Для БИХ-фильтра первого порядка.
Рисунок 1 – Пример схемы типового звена БИХ-фильтра первого порядка
yi a0 xi+a1 xi-1+b1 yi-1 – алгоритм работы цифрового фильтра первого порядка;
a0 0, a1 2, b1_1 1, b1_2 -1 – коэффициенты;
Передаточная функция БИХ-фильтра описывается выражением:
m
H(z) |
∑ak Z-k |
|
k=0 |
. |
|
|
n |
|
|
1-∑bk z-k |
|
|
k=1 |
|
Для БИХ-фильтра первого порядка (m=n=1) передаточная функция будет иметь вид:
H1_1(z) 1a-0+b1a_11 zz--11 , H1_2(z) 1a-0+b1a_12 zz--11 .
Для проверки устойчивости БИХ-фильтра необходимо решить уравнение
n
1-∑bk z-k 0. Фильтр будет устойчив, если полюсы его передаточной
k=1
функции лежат внутри единичной окружности комплексной плоскости. Для БИХ-фильтра первого порядка получается: если |b1|<1, то фильтр будет
устойчив.
Рассматриваемый БИХ-фильтр первого порядка не устойчив для обоих b1, так как единица не меньше единицы.
Импульсная характеристика h – это реакция на воздействие в виде функции Дирарка (дельта-функции). Так как z-преобразование дельта-функции равно единице только для нулевого отсчёта, то все остальные отсчёты будут равны нулю. Функцию Дирарка можно описать следующим соотношением:
i 0,1 20, k 2,3 20;
σ ‖if i<1 |
||. |
||
i |
‖ |
‖ |
|| |
|
‖ ‖1 |
|| |
|
|
‖else if i≥1|| |
||
|
‖ |
‖ |
|| |
|
‖ |
‖0 |
|| |
Импульсная характеристика рассматриваемого фильтра будет иметь вид:
|
h1_1 |
0 |
a0 |
σ +a1 0+b1_1 0=0 |
, h1_2 |
0 |
a0 |
σ +a1 0+b1_2 0=0 ; |
||||||||
|
|
|
|
|
0 |
|
|
|
|
|
0 |
|
|
|
||
h1_1 |
|
a0 σ +a1 σ +b1_1 h1_1 |
=2 |
, h1_2 |
1 |
a0 |
σ +a1 |
σ +b1_2 h1_1 |
=2 ; |
|||||||
|
1 |
|
|
1 |
|
0 |
0 |
|
|
|
1 |
0 |
0 |
|
||
h1_1 |
a0 |
σ +a1 |
σ +b1_1 h1_1 |
k-1 |
, h1_2 |
k |
a0 |
σ +a1 |
σ +b1_2 h1_2 |
. |
||||||
|
|
k |
|
|
k |
|
k-1 |
|
|
|
k |
k-1 |
|
k-1 |
||
|
|
2 |
|
|
1.8 |
|
|
1.6 |
|
|
1.4 |
|
|
1.2 |
h1_1 |
|
1 |
i |
|
|
|
0.8 |
|
|
|
0.6
0.4
0.2
0
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
i
Рисунок 2 – Импульсная характеристика БИХ-фильтра первого порядка с b1 1
|
|
2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
h |
1_2i |
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
15 |
16 |
17 |
18 |
19 |
20 |
||
|
-0.4 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-0.8 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-1.2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-1.6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-2 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
Рисунок 3 – Импульсная характеристика БИХ-фильтра первого порядка с |
||||||||||||||||||||||
|
|
|
|
|
|
|
|
|
b1 -1 |
|
|
|
|
|
|
|
|
|
|
|
||
Частотный коэффициент передачи (комплексная передаточная функция)
БИХ-фильтра первого порядка выглядит следующим образом: |
|
|||||
|
|
|
substitute,z ej 2 π FT |
|
2 |
|
H1_1 |
(FT) H1_1 |
(z)――――――→ |
|
; |
||
|
|
|||||
|
|
|
e2i |
FT π-1 |
||
|
|
|
substitute,z ej 2 π FT |
|
2 |
|
H1_2(FT) H1_2 |
(z)――――――→ |
|
; |
|||
|
FT π |
|||||
|
|
|
e2i |
+1 |
||
АЧХ рассматриваемого БИХ-фильтра первого порядка:
H1_1_АЧХ(FT) abs H1_1(FT) ;
H1_2_АЧХ(FT) abs H1_2(FT) .
0.5
31
28
25
22
19
16
H1_1_АЧХ(FT) 13
10
7
4
1
1
0 |
0.1 |
0.2 |
0.3 |
0.4FT0.5 |
0.6 |
0.7 |
0.8 |
0.9 |
1 |
||
|
|
|
|
|
|
|
|
|
|
|
|
Рисунок 4 – АЧХ БИХ-фильтра первого порядка с b1 1
|
|
|
0.5 |
|
|
|
31 |
|
|
|
|
28 |
|
|
|
|
25 |
|
|
|
|
22 |
|
|
|
|
19 |
|
|
H1_2_АЧХ |
(FT) |
16 |
|
|
13 |
|
|
||
|
|
|
|
|
|
|
10 |
|
|
|
|
7 |
|
|
|
|
4 |
|
|
|
|
1 |
|
|
|
|
1 |
|
|
|
|
0 |
0.1 0.2 0.3 0.4FT0.5 0.6 0.7 0.8 0.9 |
1 |
Рисунок 5 – АЧХ БИХ-фильтра первого порядка с b1 -1
Для БИХ-фильтра второго порядка
Рисунок 6 – Пример схемы типового звена БИХ-фильтра второго порядка
yi a0 xi+a1 xi-1+a2 xi-2+b1 yi-1+b2 yi-2 – алгоритм работы цифрового фильтра второго порядка;
a0 1, a1 1, a2 -2, b1 0.5, b2 -1 – коэффициенты;
Передаточная функция БИХ-фильтра описывается выражением:
m
∑ak Z-k
H(z) k=0 |
n |
. |
|
|
1-∑bk z-k
k=1
Для БИХ-фильтра первого порядка (m=n=2) передаточная функция будет иметь вид:
H2(z) a0+a1 z--11+a2 z--22
1-b1 z -b2 z
Для проверки устойчивости БИХ-фильтра необходимо решить уравнение
n
1-∑bk z-k 0. Фильтр будет устойчив, если полюсы его передаточной
k=1
функции лежат внутри единичной окружности комплексной плоскости.
1-b1 z |
-1 |
-b2 z |
-2 |
solve,z |
0.25 |
+0.96824583655185422129i |
|
, |
|
|
0 ――→ |
0.25 |
-0.96824583655185422129i |
|
|||
|
|
|
|
|
|
z1 abs(0.25+0.96824583655185422129i)=1 , z2 abs(0.25-0.96824583655185422129i)=1 .
Следовательно, рассматриваемый БИХ-фильтр второго порядка не устойчив для, так как единица не меньше единицы (корни вышеописанного уравнения лежат на единичной окружности, а не внутри).
Импульсная характеристика h – это реакция на воздействие в виде функции Дирарка (дельта-функции). Так как z-преобразование дельта-функции равно единице только для нулевого отсчёта, то все остальные отсчёты будут равны нулю. Функция Дирарка описана в пункте отчёта для БИХ-фильтра первого порядка.
Импульсная характеристика рассматриваемого фильтра будет иметь вид:
k 3,4 20;
h20 a0 σ0+a1 0+a2 0+b1 0+b2 0=1 , h21 a0 σ1+a1 σ0+a2 0+b1 h20+b2 0=1.5 , h22 a0 σ2+a1 σ1+a2 σ0+b1 h21+b2 h20=-2.25 , h2k a0 σk+a1 σk-1+a2 σk-2+b1 h2k-1+b2 h2k-2.
|
|
4.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3.75 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
2.25 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
h2 |
|
0.75 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
0 |
1 |
2 |
3 |
4 |
5 |
6 |
7 |
8 |
9 |
10 |
11 |
12 |
13 |
14 |
|
15 |
16 |
17 |
18 |
19 |
20 |
|
|
|
0 |
|
|
|||||||||||||||||||||
|
|
-0.75 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-1.5 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-2.25 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
-3.75 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
i |
|
|
|
|
|
|
|
|
|
|
|
|
Рисунок 7 – Импульсная характеристика БИХ-фильтра второго порядка |
|
|||||||||||||||||||||||
Частотный коэффициент передачи (комплексная передаточная функция) |
|
|||||||||||||||||||||||
БИХ-фильтра второго порядка выглядит следующим образом: |
|
|
|
|
|
|||||||||||||||||||
H2(FT) H2 |
substitute,z ej 2 π FT |
|
|
|
|
|
1.5 |
e |
6.2831853071795864769i FT |
-... |
||||||||||||||
(z)――――――→ |
|
|
|
|
|
|
|
|
|
|
|
|||||||||||||
|
|
|
|
|
|
|
|
|
|
e6.2831853071795864769i FT 2 -0.5 e6.283185307 |
||||||||||||||
АЧХ рассматриваемого БИХ-фильтра второго порядка: |
|
|
|
|
|
|
|
|
|
|||||||||||||||
|
|
|
|
|
H2_АЧХ(FT) abs H2(FT) . |
|
|
|
|
|
|
|
|
|
||||||||||
|
|
|
|
|
|
|
|
|
0.211 |
|
|
|
0.5 |
|
|
|
|
0.791 |
|
|
|
|
||
|
|
|
|
|
27 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
24 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
21 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
18 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
H2_АЧХ(FT) |
|
|
15 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
12 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
9 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
6 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0.8 |
3 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
0 |
|
0.1 |
|
0.2 |
0.3 |
|
0.4 |
0.5 |
|
0.6 |
|
0.7 |
0.8 |
|
0.9 |
|
1 |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
FT |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Рисунок 8 – АЧХ БИХ-фильтра второго порядка |
|
|
|
|
|
|||||||||||||||||
4. Исследование БИХ-фильтра на ЭВМ
Рисунок 9 – Схема исследуемого БИХ-фильтра первого порядка
Рисунок 10 – АЧХ БИХ-фильтра первого порядка с b1 1
Рисунок 11 – АЧХ БИХ-фильтра первого порядка с b1 -1
Рисунок 12 – Схема исследуемого БИХ-фильтра второго порядка
