- •Расчет аналоговых и цифровых фильтров
- •Содержание
- •1. Общая теория
- •1.1. Свойства цифровых фильтров
- •1.2. Представление цифрового фильтра в виде разностного уравнения
- •1.3 Аналоговые фильтры-прототипы
- •1.3.1. Фильтры Баттерворта
- •1.3.2. Фильтры Чебышева
- •1.3.3. Эллиптические фильтры
- •1.3.4. Фильтры Бесселя
- •1.4. Преобразования полосы частот для аналоговых фильтров
- •1.5. Преобразование полосы для цифровых фильтров
- •1.6. Методы дискретизации аналогового фильтра
- •1.6.1. Метод инвариантного преобразования импульсной характеристики
- •1.6.2. Метод билинейного преобразования
- •1.7. Методы синтеза ких-фильтров.
- •1.7.1. Прямоугольное окно
- •1.7.2. Обобщенное окно Хэмминга
- •1.7.3. Окно Блэкмана
- •1.7.4. Окно Кайзера
- •1.8. Методы синтеза бих-фильтров.
- •1.9. Методы реализации цифровых фильтров
- •1.9.1. Прямая форма
- •1.9.2. Прямая каноническая форма
- •1.9.3. Каскадная форма
- •1.9.4. Параллельная форма
- •2. Задание на курсовую работу в таблице 2.1. Представлены варианты на курсовую работу. Вариант выбрается в соответсвии с указанием преподавателя.
- •3. Порядок выполнения задания №1. Расчет аналогового фильтра.
- •Порядок фильтра.
- •Аппроксимация ачх фильтров - общие замечания.
- •Аппроксимация по Баттерворту.
- •Аппроксимация по Чебышеву первого рода.
- •Аппроксимация по Чебышеву второго рода.
- •Аппроксимация по Кауэру. Эллиптический фильтр.
- •Сравнение порядков фильтров при различных способах аппроксимации ачх. Решение уравнения порядка фильтра.
- •Порядок расчет первого задания.
- •4. Порядок выполнения задания №2. Расчет цифрового фильтра методом частотной выборки с использованием окон.
- •Порядок расчет второго задания.
- •5. Порядок выполнения задания №3. Реализация цифрового фильтра в среде Simulink.
- •6. Проектирование цифровых фильтров в среде matlab
- •6.1. Анализ фильтров и выполнение
- •Filtic – cоздание начального состояния для функции filter:
- •Freqs – частотная характеристика аналогового фильтра:
- •Freqspace – формирование последовательности отсчетов частоты:
- •Freqz – частотная характеристика цифрового фильтра.
- •Grpdelay – групповая задержка цифрового фильтра:
- •Impz – импульсный отклик цифрового фильтра:
- •Unwrap – корректировка фазовых углов:
- •Zplane – отображение нулей и полюсов цифрового фильтра:
- •7.2. Проектирование цифровых бих-фильтров besself – проектирование аналогового фильтра Бесселя:
- •Butter – проектирование цифрового и аналогового фильтров Баттерворта:
- •Cheby1 – проектирование цифрового и аналогового фильтров Чебышева – первого типа:
- •Cheby2 – проектирование цифрового и аналогового фильтров Чебышева второго типа:
- •Ellip – проектирование эллиптического цифрового и аналогового фильтров:
- •Yulewalk – проектирование рекурсивного фильтра с использованием метода наименьших квадратов по заданной амплитудно-частотной характеристике:
- •7.3. Выбор порядка бих-фильтра buttord – выбор порядка фильтра Баттерворта:
- •Cheb1ord – выбор порядка для фильтра Чебышева первого порядка:
- •Cheb2ord – выбор порядка для фильтра Чебышева первого порядка:
- •Ellipord – выбор порядка эллиптического фильтра:
- •7.4. Проектирование ких-фильтров fir1 – фильтр fir проектируется с использованием метода окна:
- •Fir2 – проектирование фильтра fir с использованием оконного метода для произвольной формы фильтра:
- •Firls – проектирование ких-фильтра с использованием минимизации ошибок методом наименьших квадратов (мнк):
- •Intfilt – расчет интерполирующего ких-фильтра:
- •Remez – синтез оптимального fir-фильтра с равномерной (чебышевской) аппроксимацией на основе алгоритма Паркса – Мак-Клелана:
- •7.5. Преобразования
- •Dftmtx – матрица дискретного преобразования Фурье (дпф):
- •Impinvar – метод инвариантной импульсной характеристики для перевода аналогового фильтра в цифровой:
- •Список литературы
1.1. Свойства цифровых фильтров
Цифровой фильтр называется стационарным, если его параметры не изменяются во времени, т. е. предварительно невозбужденный фильтр, в котором = = 0 при всех
< 0, называют стационарным
тогда и только тогда, когда
для всех возможных воздействий.Цифровой фильтр называют линейным тогда и только тогда, когда
для всех
и
– произвольных постоянных и для всех
допустимых воздействий
и
.Цифровой фильтр называют физически реализуемым, если величина отклика при
зависит только от значений входной
последовательности с номерами
.
Это означает, что импульсная характеристика
равна нулю при
.Цифровой фильтр называется устойчивым тогда и только тогда, когда реакция на ограниченное воздействие ограничена, т. е. если из
при всех
следует
при всех
.
Необходимым и достаточным условием
устойчивости фильтра является следующее
требование к его импульсной характеристике:
1.2. Представление цифрового фильтра в виде разностного уравнения
Цифровой фильтр в общем виде представляется следующим образом как разностное уравнение:
(1.2)
где
и
–
вещественные или комплексные коэффициенты.
Цифровые фильтры принято делить на два класса: нерекурсивные (НФ) и рекурсивные (РФ). Если в (1.2) все коэффициенты = 0, что соответствует отсутствию обратной связи, то фильтр является нерекурсивным и описывается уравнением
. (1.3)
Если в (1.3) хотя бы один из
коэффициентов
,
то фильтр является рекурсивным и
представляет собой устройство с обратной
связью.
Для анализа систем, описываемых
разностными уравнениями, широко
применяется
-преобразование.
Прямое
-преобразование
последовательности
определяется формулой
. (1.5)
В разностных уравнениях
существенной операцией является
единичная задержка, описываемая
оператором
,
или
(т. е. для последовательности
-преобразование
будет иметь вид
Передаточной (системной)
функцией
цифрового фильтра называется отношение
-преобразований
выходного
и входного
сигналов фильтра. Для
рекурсивного и нерекурсивного фильтров
из (1.3) и (1.4), используя (1.5), получаем:
Комплексная частотная
характеристика цифрового фильтра,
представленного в виде разностного
уравнения (1.2), может быть получена
подстановкой в выражение для передаточной
функции значения
.
Для рекурсивного фильтра общего вида
частотная характеристика будет иметь
вид
Аналогично, для нерекурсивного фильтра имеем:
1.3 Аналоговые фильтры-прототипы
Приведем расчетные формулы
для нескольких стандартных типов
аналоговых фильтров. Допустим, нужно
рассчитать аналоговый фильтр нижних
частот с частотой среза
рад/с. В качестве аппроксимируемой
функции будет использоваться квадрат
амплитудной характеристики (исключением
является фильтр Бесселя).
Будем считать, что передаточная функция аналогового фильтра является рациональной функцией переменной s следующего вида:
1.3.1. Фильтры Баттерворта
Фильтры Баттерворта нижних частот характеризуются тем, что имеют максимально гладкую амплитудную характеристику в начале координат в s-плоскости. Квадрат амплитудной характеристики нормированного (т. е. имеющего частоту среза 1 рад/с) фильтра Баттерворта равен
где
– порядок фильтра.
Аналитически продолжая функцию на всю
-плоскость,
получим:
Все полюсы этой функции
находятся на единичной окружности на
одинаковом расстоянии друг от друга в
-плоскости.
Выразим передаточную функцию
через полюсы, располагающиеся в левой
полуплоскости s:
где
,
– константа нормирования.
Можно сформулировать несколько свойств фильтров Баттерворта нижних частот:
Фильтры Баттерворта имеют только полюсы (все нули передаточных функций этих фильтров расположены на бесконечности).
На частоте
рад/с коэффициент передачи фильтра
равен
(т. е. на частоте среза их амплитудная
характеристика снижается на 3 дБ).Порядок фильтра n полностью определяет весь фильтр.
