- •1. Общая характеристика цифровых фильтров
- •1.1. Свойства цифровых фильтров
- •1.2. Представление цифрового фильтра в виде разностного уравнения
- •2. Ких-фильтры. Методы синтеза
- •2.1. Прямоугольное окно
- •2.2. Обобщенное окно Хэмминга
- •2.3. Окно Блэкмана
- •2.4. Окно Кайзера
- •3. Бих-фильтры. Методы синтеза
- •3.1. Аналоговые фильтры-прототипы
- •3.1.1. Фильтры Баттерворта
- •3.1.2. Фильтры Чебышева
- •3.1.3. Эллиптические фильтры
- •3.1.4. Фильтры Бесселя
- •3.2. Методы дискретизации аналогового фильтра
- •3.2.1. Метод инвариантного преобразования импульсной характеристики
- •3.2.2. Метод билинейного преобразования
- •3.3. Преобразования полосы частот для аналоговых фильтров
- •3.4. Преобразование полосы для цифровых фильтров
- •4. Методы реализации цифровых фильтров
- •4.1. Прямая форма
- •4.2. Прямая каноническая форма
- •4.3. Каскадная форма
- •4.4. Параллельная форма
- •5. Работа в среде matlab
- •5.1. Ввод матриц
- •5.2. Матричные операции
- •5.3. Операторы, выражения и переменные
- •5.4. Функции построения матриц
- •5.5. Операторы for, while, if
- •5.6. Скалярные функции
- •5.7. Векторные функции
- •5.8. Mатричные функции
- •5.9. Подматрицы и форма записи с двоеточием
- •5.10. M-файлы
- •5.11. Текстовые строки, сообщения об ошибках, ввод данных
- •5.12. Сравнение эффективности алгоритмов
- •5.13. Формат вывода
- •5.14. Протокол
- •5.15. Графика
- •6. Записная книжка matlab для Microsoft Word
- •6.1. Начало работы
- •6.2. Краткий пример использования Записной книжки
- •7. Проектирование цифровых фильтров в среде matlab
- •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. Преобразования czt–z-преобразование по спиральному контуру:
- •Dct– дискретное косинусное преобразование:
- •Impinvar – метод инвариантной импульсной характеристики для перевода аналогового фильтра в цифровой:
- •8. Примеры проектирования цифровых фильтров
- •8.1. Генерация входной последовательности сигнала
- •8.2. Реализация спектрального анализа с использованием бпф
- •8.3. Синтез цифрового ких-фильтра
- •8.4. Реализация цифровой фильтрации
- •Список литературы
- •Содержание
- •197376, С.-Петербург, ул. Проф. Попова, 5
3.2.2. Метод билинейного преобразования
Билинейное z-преобразование использует следующую замену:
(3.5)
На рис.3.3 показано, каким образом s-плоскость отображается в z-плоскость.
Рис. 3.3
Вся ось jW из s-плоскости отображается в единичную окружность на z-плоскости, левая полуплоскость s отображается в единичный круг, а правая полуплоскость s – в область, расположенную вне единичного круга на z-плоскости.
Из формулы (3.5) найдем выражение для z:
.
При s = jW получим:
.
Отсюда видно, что . При= 0 имеемz = 1 и при W = ¥ z = –1, в промежутке z монотонно меняется от 0 до . Подставив в формулу s = s + jW, получим:
.
При s < 0 (для левой полуплоскости) , т. е. точки расположены внутри единичной окружности.
При билинейном преобразовании передаточная функция цифрового фильтра H(z) рассчитывается с помощью алгебраической подстановки, т. е.
.
Из этого соотношения видно, что порядки знаменателей функций H(z) и H(S) совпадают, но порядки числителей могут различаться. Например, передаточная функция
имеет числитель нулевого порядка, а знаменатель – первого порядка. В то же время, получаемая методом билинейного преобразования функция равна
где и числитель и знаменатель 1-го порядка. Причиной этого является то, что H(S) имеет нуль на бесконечности (S = ¥), который при билинейном преобразовании отображается в точку z = –1.
Так как в единичную окружность на z-плоскости отображается вся ось jW из s-плоскости, то эффекты, связанные с наложениями в частотной характеристике цифрового фильтра, характерные для метода инвариантного преобразования импульсной характеристики, в данном случае наблюдаться не будут. Однако соотношение между частотами аналогового фильтра и цифрового фильтра оказывается нелинейным.
Рассмотрим характер этой нелинейности. Пусть в (3.5) и S = jW тогда
или
откуда
;
При небольших w отображение почти линейно, однако для основной части частотной шкалы оно существенно нелинейно и значительно ограничивает область применения билинейного преобразования. Существует, правда, довольно большой класс фильтров, для которых частотная деформация может быть скомпенсирована. К ним относятся фильтры нижних частот, высоких частот, полосовые и режекторные. Эффекты нелинейности соотношения между частотными шкалами аналогового и цифрового фильтров удается учесть лишь в том случае, когда частотная характеристика аналогового фильтра имеет вид ступенчатой функции. Кроме того, при билинейном преобразовании ни импульсная, ни фазовая характеристики аналогового фильтров не будут совпадать.
3.3. Преобразования полосы частот для аналоговых фильтров
Существует много различных методов преобразования фильтров нижних частот с частотой среза, равной 1 рад/с, в другой фильтр нижних частот (имеющий другую частоту среза), а также в фильтр верхних частот, полосовой или режекторный. Перечислим наиболее простые преобразования.
Фильтр нижних частот фильтр нижних частот:
.
Фильтр нижних частот фильтр верхних частот:
.
Фильтр нижних частот полосовой фильтр:
.
Фильтр нижних частот режекторный фильтр:
.
Здесь Wl – нижняя частота среза, Wu – верхняя частота среза.