Министерство цифрового развития, связи и массовых коммуникаций Российской Федерации
Ордена Трудового Красного Знамени федеральное государственное бюджетное образовательное учреждение высшего образования МОСКОВСКИЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
СВЯЗИ И ИНФОРМАТИКИ (МТУСИ)
Факультет "Радио и телевидение" Кафедра "Теория электрических цепей"
САМОСТОЯТЕЛЬНАЯ РАБОТА
по теме «Цифровые фильтры»
Дисциплина "Основы компьютерного анализа электрических цепей"
Москва, 2023
Оглавление |
|
|
1. ОПИСАНИЕ ЗАДАНИЯ ДЛЯ САМОСТОЯТЕЛЬНОЙ РАБОТЫ ПО |
|
|
ТЕМЕ «ЦИФРОВЫЕ ФИЛЬТРЫ» ................................................................... |
4 |
|
1.1. |
ЦЕЛЬ САМОСТОЯТЕЛЬНОЙ РАБОТЫ............................................... |
4 |
1.2. |
ИСХОДНЫЕ ДАННЫЕ (ПРИМЕР) ....................................................... |
4 |
1.3. |
ЗАДАНИЕ .................................................................................................. |
4 |
1.3.1. ПРОВЕРИТЬ УСТОЙЧИВОСТЬ ФИЛЬТРА ....................................... |
4 |
|
1.3.2. НАЙТИ ПОЛЮСА И НУЛИ ПЕРЕДАТОЧНЫХ ФУНКЦИЙ |
|
|
|
ФИЛЬТРА................................................................................................. |
4 |
1.3.3. ИЗОБРАЗИТЬ ФУНКЦИОНАЛЬНУЮ СХЕМУ ФИЛЬТРА ............. |
4 |
|
1.3.4. РАССЧИТАТЬ ИМПУЛЬСНУЮ РЕАКЦИЮ ФИЛЬТРА ................ |
4 |
|
1.3.5. РАССЧИТАТЬ АЧХ ФИЛЬТРА ............................................................ |
4 |
|
2. ПРИМЕР ВЫПОЛНЕНИЯ САМОСТОЯТЕЛЬНОЙ РАБОТЫ .................... |
5 |
|
2.1. |
АНАЛИЗ КИХ-ФИЛЬТРА ....................................................................... |
5 |
2.1.1. УСТОЙЧИВОСТЬ. .................................................................................. |
5 |
|
2.1.2. НАХОЖДЕНИЕ КОРНЕЙ...................................................................... |
5 |
|
2.1.3. ФУНКЦИОНАЛЬНАЯ СХЕМА ............................................................ |
7 |
|
2.1.4. ИМПУЛЬСНАЯ РЕАКЦИЯ ФИЛЬТРА ............................................... |
8 |
|
2.1.5. АЧХ ФИЛЬТРА ....................................................................................... |
9 |
|
2.2. |
АНАЛИЗ БИХ-ФИЛЬТРА ..................................................................... |
10 |
2.2.1. УСТОЙЧИВОСТЬ. ................................................................................ |
10 |
|
2.2.2. НАХОЖДЕНИЕ КОРНЕЙ.................................................................... |
11 |
|
2.2.3. ФУНКЦИОНАЛЬНАЯ СХЕМА .......................................................... |
11 |
|
2.2.4. ИМПУЛЬСНАЯ РЕАКЦИЯ ФИЛЬТРА ............................................. |
13 |
|
2.2.5. АЧХ ФИЛЬТРА ..................................................................................... |
15 |
|
|
2 |
|
2.3. АНАЛИЗ ПОСЛЕДОВАТЕЛЬНОГО СОЕДИНЕНИЯ КИХ- |
|
ФИЛЬТРА И БИХ-ФИЛЬТРА ............................................................... |
16 |
2.3.1. УСТОЙЧИВОСТЬ. ................................................................................ |
16 |
2.3.2. НАХОЖДЕНИЕ КОРНЕЙ.................................................................... |
17 |
2.3.3. ФУНКЦИОНАЛЬНАЯ СХЕМА .......................................................... |
17 |
2.3.4. ИМПУЛЬСНАЯ РЕАКЦИЯ ФИЛЬТРА ............................................. |
18 |
2.3.5. АЧХ ФИЛЬТРА ..................................................................................... |
20 |
3. ВЫВОДЫ .......................................................................................................... |
21 |
ПРИЛОЖЕНИЕ 1. ИНДИВИДУАЛЬНЫЕ ЗАДАНИЯ ............................... |
22 |
3
1. Описание задания для самостоятельной работы по теме
«Цифровые фильтры»
1.1.Цель самостоятельной работы
Цель самостоятельной работы на примере заданных передаточных функций КИХ-фильтра и БИХ-фильтра освоить методики проверки устойчивости фильтров, вычисления импульсных и амплитудно-частотных характеристик.
1.2.Исходные данные (пример)1
Заданы передаточные функции КИХ фильтра |
|
|
|
|
|
||
HКИХ z 1 0.9z 1 |
0.1z 2 1 0.1z 1 0.72z 2 |
|
(1) |
||||
и БИХ-фильтра |
|
|
|
|
|
|
|
HБИХ z |
|
1 |
|
|
|
. |
(2) |
|
|
|
|
|
|||
1 1.25z 1 |
0.25z 2 1 z 1 |
0.24z 2 |
|
||||
|
|
|
|
|
|
||
Период дискретизации T 1 мс
1.3.Задание
1.3.1.Проверить устойчивость фильтра
1.3.2.Найти полюса и нули передаточных функций фильтра
1.3.3.Изобразить функциональную схему фильтра
1.3.4.Рассчитать импульсную реакцию фильтра
1.3.5.Рассчитать АЧХ фильтра
Пункты 1.3.1-1.3.5 выполняются отдельно для КИХ-фильтра, БИХ-
фильтра и последовательно-соединѐнных КИХ и БИХ фильтров.
1 Все исходные данные - индивидуальные. Определяются номером в журнале и заданы в разделе Индивидуальные задания стр. 22.
4
2. Пример выполнения самостоятельной работы
2.1. |
Анализ КИХ-фильтра |
|
|
Передаточная функция КИХ-фильтра, |
согласно (1), |
задана в |
|
следующем виде |
|
|
|
|
HКИХ z 1 0.9z 1 0.1z 2 1 0.1z 1 0.72z 2 |
(3) |
|
2.1.1.Устойчивость.
Поскольку в передаточной функции КИХ фильтра нет знаменателя,
следовательно, у неѐ нет полюсов. Поэтому КИХ-фильтр является всегда устойчивым
2.1.2.Нахождение корней.
Нулями КИХ-фильтра являются все корни полиномов его передаточной функции. Вычисление корней может быть осуществлено с использованием математических программ, либо вручную.
Программно2. Для нахождения корней полиномов можно использовать функции пакетов MatLab – roots(a) или MathCad – polyroots(a), где a – вектор коэффициентов полинома.
Передаточная функция КИХ-фильтра задана в виде произведения двух полиномов, поэтому корнями передаточной функции будут все корни этих полиномов.
Для первого полинома 1 0.9z 1 0.1z 2 имеем
roots([1 -0.9 -0.1]) =
ans =
1.0000
-0.1000
2 При выполнение самостоятельной работы можно использовать любой из рассмотренных методов: программный или вручную.
5
Для второго полинома 1 0.1z 1 0.72z 2 имеем
roots([1 0.1 -0.72]) |
|
|
|
ans = |
|
|
|
-0.9000 |
|
|
|
0.8000 |
|
|
|
Таким образом, нулями передаточной функции будут значения |
|
||
z0 [1, 0.1, |
0.9, |
0.8] . |
(4) |
Вручную. Передаточная функция |
КИХ-фильтра задана в |
виде |
|
произведения полиномов второго порядка (квадратный трѐхчлен),
поэтому для нахождения корней можно вычислить вручную корни каждого квадратного трѐхчлена.
Для первого полинома 1 0.9z 1 0.1z 2 имеем уравнение
z2 0.9z 0.1 0
Решение
D 0.9 2 4 1 0.1 1.21 0
|
|
0.9 |
|
|
1 |
|
|
|
D |
||||
z1,2 |
|
|
|
|
|
|
|
|
|
||||
|
|
2 |
|
0.1 |
||
Для второго полинома 1 0.1z 1 0.72z 2 имеем уравнение
z 0.1z1 0.72 0 . |
(5) |
Решение
D 0.1 2 4 1 0.72 2.89 0
|
|
0.1 |
|
|
0.8 |
|
|
D |
|||
z1,2 |
|
|
|
|
|
|
|
||||
|
|
2 |
|
0.9 |
|
Таким образом, нулями передаточной функции будут значения
zКИХ [1, 0.1, |
0.8, 0.9] . |
(6) |
6
2.1.3.Функциональная схема
Для построения функциональной схемы нужно передаточную функцию
представить в виде одного полинома, т.е. перемножить два полинома
HКИХ z 1 0.9z 1 0.1z 2 1 0.1z 1 0.72z 21 0.9z 1 0.1z 2 0.1z 1 1 0.9z 1 0.1z 2
0.72z 2 1 0.9z 1 0.1z 2
1 0.9z 1 0.1z 2 0.1z 1 0.09z 2 0.01z 30.72z 2 0.648z 3 0.072z 4
1 0.8z 1 0.92z 2 0.647z 3 0.072z 4
Таким образом, результирующая передаточная функция КИХ-фильтра будет иметь вид
HКИХ z 1 0.8z 1 0.92z 2 0.647z 3 0.072z 4 . |
(7) |
Имеем фильтр четвѐртого порядка, который реализуется в виде 5-
звенного фильтра. Схема фильтра приведена на рисунке 1.
xn
Z-1 |
Z-1 |
Z-1 |
Z-1 |
|
|
|
|
0.072 |
|
|
|
|
0.647 |
yn |
|
|
|
|
|
|
|
0.92 |
|
|
|
|
|
|
|
|
0.8 |
|
|
|
1 |
|
|
|
|
Рисунок 1 – Функциональная схема КИХ-фильтра |
|
|||
Следует отметить, что вычисление результирующей передаточной функции можно выполнить программно с использованием, например, пакета
Matlab двумя способами.
1. Перемножением полиномов с помощью функции conv(a,b), где a и b
векторы коэффициентов перемножаемых полиномов:
conv([1 -0.9 -0.1],[1 +0.1 -0.72])
ans =
7
1.0000 -0.8000 -0.9100 0.6380 0.0720.
2. Реконструкция полинома по его корням с использованием функции poly(roots), где roots – вектор всех корней всех перемножаемых полиномов:
poly([1 -0.1 0.8 -0.9])
ans =
1.0000 |
-0.8000 |
-0.9100 |
0.6380 |
0.0720 |
2.1.4.Импульсная реакция фильтра
Импульсная реакция фильтра – это выходной сигнал фильтра при воздействии на входе единичного импульса. Выходной сигнал КИХ-фильтра с передаточной функцией вида (7) определяется следующим выражением:
yn a0 xn a1xn 1 a2 xn 2 |
a3xn 3 |
a4 xn 4 , |
(8) |
||||
где a0 1, a1 0.8, a2 |
0.92, |
a3 0.647, |
a4 0.72 |
|
|
||
Для вычисления передаточной функции нужно сформировать входной |
|||||||
сигнал |
|
|
|
|
|
|
|
|
|
1 |
для |
|
n 1 |
|
|
xn |
|
0 |
для всех |
n 1 |
, |
(9) |
|
|
|
|
|
||||
и вычислить yn для |
n 1,2,....N 1 , где N – число звеньев в фильтре. В |
||||||
данном случае N 5 . |
|
|
|
|
|
|
|
Для вычисления импульсной характеристики можно использовать |
|||||||
функции пакетов MatLab (или MathCad). |
Например, функцию filter(). Для |
||||||
КИХ-фильтра это выглядит следующим образом: |
|
|
|
||||
x=zeros(1,6);
x(1)=1;
y=filter(a,1,x) y =
1.0000 |
-0.8000 |
-0.9100 |
0.6380 |
0.0720 |
0 |
figure(1) stem(y,'LineWidth',2); grid
Результат расчѐта выводится на графике и приведѐн на рисунке 2.
8
Рисунок 2 – Импульсная реакция КИХ фильтра
2.1.5. |
АЧХ фильтра3 |
|
|
|
|
|
|
|
||
Для вычисления АЧХ фильтра нужно вычислить модуль передаточной |
||||||||||
функцию, |
как функцию от частоты f |
в диапазоне |
от 0 |
до |
1 |
, т.е. |
|
в |
||
|
||||||||||
|
|
|
|
|
|
T |
|
|
||
выражении для передаточной функции заменить z e j 2 fT (или |
z 1 e j 2 fT ) |
|||||||||
где T - заданный период дискретизации. |
Для функции (7) имеем |
|
|
|||||||
HКИХ z e j 2 fT |
|
1 0.8e j 2 fT 0.92e j 4 fT 0.647e j 6 fT 0.072e j8 fT |
|
.(10) |
||||||
|
|
|||||||||
АЧХ КИХ фильтра в диапазоне от 0 до 1000 Гц |
приведена на рисунке |
|||||||||
3. |
|
|
|
|
|
|
|
|
|
|
3 Для построения АЧХ фильтра число точек по частоте должно быть не менее 100.
9
Рисунок 3 – АЧХ КИХ-фильтра |
2.2.Анализ БИХ-фильтра
Передаточная функция КИХ-фильтр задана в следующем виде
HБИХ z |
|
1 |
|
|
|
(11) |
|
|
|
|
|
||
1 1.25z 1 |
0.25z 2 |
1 z 1 |
0.24z 2 |
|
||
|
|
|
|
|
|
2.2.1.Устойчивость.
Поскольку в заданной передаточной функции БИХ-фильтра числитель не зависит от z (в числителе полином нулевого порядка), то все корни знаменателя будут полюсами передаточной функции. Для определения устойчивости нужно будет вычислить корни знаменателя. Это сделано в следующем разделе.
10
