- •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
Министерство образования РФ
______________________
Санкт-Петербургский государственный электротехнический
университет "ЛЭТИ"
______________________
А. В. Бадейкин В. В. Геппенер И. А. Корнеев
СИНТЕЗ ЦИФРОВЫХ ФИЛЬТРОВ
С ИСПОЛЬЗОВАНИЕМ ПАКЕТА ПРОГРАММ MATLAB
Учебное пособие
Санкт-Петербург
Издательство СПбГЭТУ "ЛЭТИ"
2001
УДК 621.372.542
ББК 3.973.233–018я7
Б15
Бадейкин А. В., Геппенер В. В., Корнеев И. А. Cинтез цифровых фильтров c использованием пакета программ MATLAB: Учеб. пособие. СПб.: Изд-во СПбГЭТУ "ЛЭТИ", 2001. 72 с.
Содержит основные теоретические сведения о методах синтеза цифровых фильтров. Рассмотрены вопросы использования пакета MATLAB для решения задач синтеза цифровых фильтров на ПЭВМ. Приведены основные сведения по базовым функциям, используемым в задачах синтеза цифровых фильтров. Приведены примеры расчета цифровых фильтров различного типа.
Предназначено для студентов специальности 220400 "Программное обеспечение вычислительной техники и автоматизированных систем".
Рецензенты: кафедра "Радиоэлектронные системы управления" (БГТУ "Военмех" им. Д. Ф. Устинова); канд. техн. наук Н. А. Федоров (ОАО "Радиоавионика").
Утверждено
редакционно-издательским советом университета
в качестве учебного пособия
ISBN 5-7629-0404-0 © СПбГЭТУ "ЛЭТИ", 2001
1. Общая характеристика цифровых фильтров
Различают два общих класса сигналов: аналоговые и дискретные. Аналоговым сигналом называется сигнал, определенный для каждого момента времени, дискретным сигналом – сигнал, определенный только в дискретные моменты времени. Как дискретный, так и аналоговый сигналы могут быть однозначно представлены некоторыми функциями частоты, которые называются их частотными спектрами.
Фильтрацией называется процесс изменения частотного спектра сигнала в некотором желаемом направлении. Этот процесс может привести к усилению или ослаблению частотных составляющих в некотором диапазоне частот, к подавлению или выделению какой-нибудь конкретной составляющей и т. п.
Цифровым фильтром называется цифровая система, которую можно использовать для фильтрации дискретных сигналов. Он может быть реализован программным методом или с помощью специальной аппаратуры, и в каждом из этих случаев цифровой фильтр можно применить для фильтрации сигналов как в реальном времени, так и предварительно записанных.
Цифровой фильтр можно представить структурной схемой, изображенной на рис. 1.1. На этой схеме x(n) и y(n) – соответственно, входное воздействие и реакция фильтра на это воздействие. Функционально они связаны соотношением
,
где вид оператора зависит от свойств конкретной системы.
Рис. 1.1
Реакцию цифрового фильтра на произвольное воздействие можно представить с помощью импульсной характеристики фильтра. Допустим, что x(n) – входная, а y(n) – выходная последовательности фильтра и пусть h(n) – отклик на единичный импульс, называемый импульсной характеристикой. Тогда
.
Таким образом, x(n) и y(n) связаны соотношением типа свертки. Частотная характеристика фильтра определяется следующим выражением:
. (1.1)
Поскольку частотная характеристика является периодической функцией частоты , равенство (1.1) можно рассматривать как разложение в ряд Фурье, причем коэффициенты являются одновременно отсчетами импульсной характеристики. Согласно теории рядов Фурье, коэффициентыh(n) могут быть выражены через :
.
Из этого соотношения видно, что h(n) по существу является суперпозицией синусоид с амплитудами, которые можно представить следующим образом:
.
Выражениеназывают амплитудной характеристикой фильтра, а– фазовой характеристикой фильтра.
1.1. Свойства цифровых фильтров
Дадим несколько определений, посвященных цифровым фильтрам.
Цифровой фильтр называется стационарным, если его параметры не изменяются во времени, т. е. предварительно невозбужденный фильтр, в котором x(n) = y(n) = 0 при всех n < 0, называют стационарным тогда и только тогда, когда для всех возможных воздействий.
Цифровой фильтр называют линейным тогда и только тогда, когда для всехa и b – произвольных постоянных и для всех допустимых воздействий x1(n) и x2(n).
Цифровой фильтр называют физически реализуемым, если величина отклика при n = n0 зависит только от значений входной последовательности с номерами n £ n0. Это означает, что импульсная характеристика h(n) равна нулю при n < 0.
Цифровой фильтр называется устойчивым тогда и только тогда, когда реакция на ограниченное воздействие ограничена, т. е. если из при всехn следует при всехn. Необходимым и достаточным условием устойчивости фильтра является следующее требование к его импульсной характеристике:
.