- •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.4. Преобразование полосы для цифровых фильтров
По аналогии с фильтрами непрерывного времени существует несколько простых преобразований цифрового фильтра нижних частот (с частотой среза wc) в другой фильтр нижних частот (с другой частотой среза wu), а также в цифровой фильтр верхних частот, полосовой или режекторный. Далее приведены формулы для этих преобразований.
Фильтр нижних частот фильтр нижних частот:
,
где ;wu – заданная частота среза фильтра нижних частот.
Фильтр нижних частот фильтр верхних частот:
,
где wu – заданная частота среза фильтра верхних частот.
Фильтр нижних частот полосовой фильтр:
Здесь ;;w0 – центральная частота полосового фильтра.
Фильтр нижних частот режекторный фильтр:
Здесь ; ; w0 – центральная частота режекторного фильтра.
4. Методы реализации цифровых фильтров
Цифровые фильтры с заданной передаточной функцией можно построить различными способами. В любом реальном фильтре шумы и погрешности, появляющиеся при квантовании, существенно зависят от структуры фильтра. Прежде всего все фильтры можно разделить на два больших класса:
рекурсивные;
нерекурсивные.
Для рекурсивных фильтров соотношение между входной последовательностью {x(n)} и откликом фильтра {y(n)} может быть записано следующим образом: т. е. текущий отсчет отклика y(n) определяется не только текущим и предшествующим значениями входной последовательности, но и предшествующими отсчетами отклика.
В нерекурсивных фильтрах связь между входной последовательностью и откликом имеет вид , т. е. текущий отсчет отклика зависит от текущего и предшествующих значений входной последовательности.
Реализация может осуществляться на основе следующих форм построения схем фильтра:
прямой;
канонической прямой;
каскадной;
параллельной.
4.1. Прямая форма
Рассмотрим передаточную функцию N-го порядка вида
, (4.1)
причем b0 = 1. Приведя это равенство к общему знаменателю, получим:
или
.
Если рассматривать члены вида z–kY(z) как обратные z-преобразования последовательностей y(n – k), то, взяв обратные z-преобразования обеих частей последнего равенства, можно получить искомое разностное уравнение
.
Поскольку b0 = 1, разностное уравнение можно решить относительно y(n):
.
Простая структура реализации данного разностного уравнения, называемая прямой формой, показана на рис. 4.1.
Рис. 4.1
В ней для образования цепей используются раздельные элементы задержки. Характерными чертами этой структуры являются ее простота и непосредственная связь с z-преобразованием. Однако эта структура очень чувствительна к квантованию коэффициентов.
4.2. Прямая каноническая форма
Запишем формулу (4.1) в следующем виде:
.
Цифровой фильтр, соответствующий этой формуле, состоит из двух последовательно соединенных фильтров с коэффициентами передачи, соответственно, H1(z) и H2(z). Первый фильтр имеет только полюсы, а второй – только нули. Если записать
то получится пара разностных уравнений (в предположении, что b0 = 1)
которые можно реализовать, как показано на рис. 4.2.
Рис. 4.2
Рис. 4.3
Поскольку в цепях, соответствующих H1(z) и H2(z), сигнал w(n) задерживается одинаково, то для построения фильтра достаточно использовать один набор элементов задержки (рис. 4.3). Такую структуру называют канонической формой, так как в ней используется минимальное количество сумматоров, умножителей и элементов задержки.