- •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. Необходимым
	и достаточным условием устойчивости
	фильтра является следующее требование
	к его импульсной характеристике:
	
.
