
- •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
2.1. Прямоугольное окно
N-точечное прямоугольное окно соответствует простому усечению ряда Фурье. Оно описывается весовой функцией:
Частотная характеристика прямоугольного окна описывается соотношением
.
Таким образом,
.
2.2. Обобщенное окно Хэмминга
Формула обобщенного окна Хэмминга имеет вид
причем,
лежит в пределах 0£
a £
1. Случай
=
0.5 соответствует окну Ханна, а случай
=
0.54 – окну Хэмминга.
Обобщенное окно Хэмминга для всех n может быть представлено в виде формулы
.
Следовательно, частотная
характеристика обобщенного окна Хэмминга
равна круговой свертке частотной
характеристики прямоугольного окна
с последовательностью импульсов и может
быть записана в виде
,
откуда
.
Ширина главного лепестка
частотной характеристики окна Хэмминга
в два раза больше, чем для прямоугольного
окна. Уровень боковых лепестков в случае
окна Хэмминга значительно ниже, чем у
характеристики прямоугольного окна.
Боковые лепестки функций
находятся в противофазе с боковыми
лепестками
,
поэтому общий уровень боковых лепестков
значительно уменьшается.
2.3. Окно Блэкмана
Функция окна Блэкмана подобна функциям окон Ханна и Хэмминга и определяется выражением
Дополнительный косинусоидальный компонент обеспечивает еще более значительное уменьшение пульсаций Гиббса. Ширина главного лепестка увеличивается и оказывается в три раза больше, чем у прямоугольного окна.
2.4. Окно Кайзера
Задача расчета хороших окон фактически сводится к математической задаче отыскания ограниченных во времени функций, преобразование Фурье которых наилучшим образом аппроксимируют функции, ограниченные по частоте, т. е. имеющие минимальную энергию за пределами заданного интервала частот. При решении этой задачи в замкнутой форме для непрерывных функций времени был выведен класс так называемых вытянутых сфероидальных волновых функций. Эти функции имеют достаточно сложный вид, поэтому Кайзер в качестве наилучшего окна предложил относительно простую аппроксимацию этих функций. Эта аппроксимация, названная окном Кайзера, имеет вид
,
где b – константа, определяющая компромисс между максимальным уровнем боковых лепестков и шириной главного лепестка (или долей энергии в главном лепестке) частотной характеристики окна; I0(x) – функция Бесселя нулевого порядка.
Частотная характеристика дискретного окна Кайзера в замкнутой форме не получена, но Кайзер показал, что для непрерывной функции окна частотная характеристика пропорциональна
,
где величина wb приблизительно равна ширине главного лепестка частотной характеристики.
Окно Кайзера является, по существу, оптимальным окном в том смысле, что оно представляет собой последовательность конечной длины, которая имеет минимум энергии спектра за пределами некоторой заданной частоты.
3. Бих-фильтры. Методы синтеза
Фильтром с бесконечной импульсной характеристикой (БИХ-фильтром) называют фильтр, длина импульсной характеристики которого не ограничена справа или слева.
Будем рассматривать БИХ-фильтры при условии, что они являются физически реализуемыми и устойчивыми. Для импульсных характеристик таких фильтров h(n) справедливы следующие ограничения:
Наиболее общая форма записи z-преобразования импульсной характеристики БИХ-фильтров имеет вид
.
Предположим, что M £ N. Системы, удовлетворяющие этому условию, называют системами N-го порядка.
Решение задачи расчета фильтров сводится к нахождению значений его коэффициентов ai и bi, обеспечивающих аппроксимацию заданных характеристик фильтра. Таким образом, задача расчета фильтра в значительной степени сводится к задаче аппроксимации и может быть решена чисто математическими методами.
Наиболее распространенным методом расчета цифровых БИХ-фильтров является метод дискретизации аналогового фильтра, удовлетворяющего заданным требованиям. При расчете цифровых фильтров верхних частот, полосовых и режекторных, используются два подхода, представленные на рис. 3.1
Рис. 3.1
В первом случае нормализованный аналоговый фильтр предварительно преобразуется в другой аналоговый фильтр, из которого путем дискретизации рассчитывается фильтр с заданными характеристиками. Во втором случае нормализованный фильтр нижних частот дискретизуется сразу же, а затем преобразованием его полосы частот формируется цифровой фильтр с заданными характеристиками.