Последовательности и их представления.
Дискретные сигналы определяются для дискретных значений независимой переменной - времени.
Время квантуется обычно равномерно:
T=nT;
Представление произвольных последовательностей.
Рассмотрим последовательность
а(0), а(1), а(2), … , а(n);
где а(n) – величина n-го элемента.
Такую последовательность можно представить:
Дискретная система.
Алгоритм преобразования
y(n)=Ф[x(n)]
Т.О., чтобы вычислить выходную последовательность системы на заданную входную, необходимо вычислить цифровую свертку входной последовательности и импульсной характеристики системы.
Физическая реализуемость.
ЛПП – система физически реализуема, если отклик при n=n0 зависит только от отсчетов входной последовательности с номерами
Это означает, что импульсная характеристика h(n) системы равна нулю, при всех
Примеры систем физически нереализуемых.
- Идеальный ФНЧ
- Идеальный дифференциатор.
Устойчивость ЛПП.
Если при любой ограниченной входной последовательности выходная последовательность так же ограничена.
Т.о. последовательность y (0) не ограничена, так, что неравенство (*) является необходимым условием устойчивости.
Доказательство достаточности.
Пусть (*) выполняется. На вход поступает ограниченная последовательность
Тогда:
Последовательность y(n) ограничена, поэтому система устойчива, что и требовалось доказать.
Примеры импульсных характеристик.
При <1 соответствующая система является устойчивой.
Структурная схема цифровых фильтров.
Уравнение цифровой свертки:
Z-преобразование:
где Y(z) – z-преобразование выходной последовательности,
X(z) - z-преобразование входной последовательности,
H(z) - z-преобразование импульсной характеристики фильтра.
Тогда:
Выразим H(z) в виде дробно-рационального полинома от переменной z-1, т.е.:
причем b0 д.б. равно 1.
Приведем к общему знаменателю:
Или:
Взяв обратное z-преобразование от левой и правой части (обратное преобразование – не вводилось понятие) можно перейти к последовательностям.
т.к. b0=1, тогда:
- уравнение цифрового фильтра.
Простая форма 1 его реализации:
Запишем:
Обозначив первый и второй множители за H1(z) и H2(z) соответственно, можно получить два фильтра:
*Прямая форма 2 (неканоническая)
Очевидно, что можно использовать один набор элементов задержки, т.к. задерживается один и тот же сигнал:
*Прямая форма 2 или каноническая форма
Можно записать:
Получим третью структуру построения цифрового фильтра. Множители Hi(z) соответствуют либо блокам второго порядка:
либо блокам первого порядка, т.е.:
Последовательная (каскадная) форма.
Разложим на простые дроби:
где Hi(z) – соответствуют или блокам второго порядка:
или блокам первого порядка:
Нерекурсивные цифровые фильтры.
Порядок расчета фильтров. Свойства КИХ-фильтров. Характеристики КИХ-фильтров с ЛФХ. Частотные характеристики КИХ-фильтров с ЛФХ.
Квадрат амплитудной характеристики
При расчете БИХ фильтров с использованием аппроксимации только АЧХ (без учета фазовой характеристики_ удобнее всего оперировать с квадратом амплитудной характеристики:
Фазовая характеристика
Метод инвариантного преобразования импульсной характеристики.
Этот метод также называется методом стандартного преобразования.
В качестве импульсной характеристики рассчитываемого цифрового фильтра используется дискретизированная импульсная характеристика соответствующего аналогового фильтра.
Передаточная характеристика аналогового фильтра:
(*) Разложим на простые дроби,
где
Предположение, что M<N (* - разложение на простые дроби) – является в данном случае обязательным, в противном случае наложения в частотной области станут недопустимыми.
Импульсная характеристика h(t) аналогового фильтра:
Дискретизируя ее получим.
,
где Т – период дискритизации.
Найдем Z – преобразование.
Изменим порядок суммирования.
(**)
Сравним (*) (**), переход от H(s) к H(z) осуществляется с помощью отображения, при котором используется замена:
Частотная характеристика цифрового фильтра, рассчитанная методом инвариантного преобразования импульсной характеристики, образуется путем положения частотной характеристики дискретизуемого аналогового фильтра.
Т.о.
,
где
- угловая частота дискретизации цифрового фильтра
Каждая горизонтальная полоса, шириной из S – плоскости отображается на Z – плоскость.
Все смежные полосы из S – плоскости будут при отображении накладываться друг на друга в Z – плоскости.
Отсюда следует, что для того что бы частотные характеристики исходного аналогового фильтра и рассчитываемого методом инвариантного преобразования импульсной характеристики цифрового фильтра соответствовали друг другу, необходимо что бы полоса пропускания аналогового фильтра находилась в пределах диапазона.
Для выполнения этого условия необходимо до начала преобразования (т.е. до АЦП) вводить дополнительный фильтр нижних частот, гарантирующий соответствующее ограничение полосы пропускания аналогового фильтра.
Прямые методы рассчета цифровых фильтров.
Существуют так же прямые методы расчета цифровых фильтров в частотной или временной областях, без аналоговых прототипов.