Передаточные функции фильтров

Z-преобразование. Удобным методом решения разностных уравнений линейных систем является z-преобразование. Применяя z-преобразование к обеим частям равенства (4.1), c учетом сдвига функций (y(k-m) ó z^m Y(z)), получаем:

Y(z) a_mz^m = X(z) b_nzⁿ, (8.1)

где X(z),Y(z)- соответствующие z-образы входного и выходного сигнала. Отсюда, полагая a_o = 1, получаем в общей форме функцию связи выхода фильтра с его входом - уравнение передаточной функции системы в z-области:

H(z) = Y(z)/X(z) = b_nzⁿ (1+ a_mz^m). (8.2)

Для НЦФ:

H(z) = b_nzⁿ. (8.3)

При проектировании фильтров исходной, как правило, является частотная передаточная функция фильтра H(ω), по которой вычисляется ее Z-образ H(z) и обратным переходом в пространство сигналов определяется алгоритм обработки данных. В общей форме для выходных сигналов фильтра:

Y(z) = H(z)·X(z).

Y(z)·(1+ a_m·z^m) = X(z)· b_n·zⁿ

Y(z) = X(z)· b_n·zⁿ – Y(z)· a_m·z^m. (8.4)

После обратного Z-преобразования выражения (8.4):

y(k) = b_n·x(k-n) – a_m·y(k-m). (8.5)

При подаче на вход фильтра единичного импульса Кронекера _о, имеющего z-образ (z) = zⁿ = 1, сигнал на выходе фильтра будет представлять собой импульсную реакцию фильтра y(k) ≡ h(k), при этом:

H(z) = Y(z)/(z) = Y(z) = TZ[y(k)] = h(k)z^k, (8.6)

т.е. передаточная функция фильтра является z-образом ее импульсной реакции. При обратном z-преобразовании передаточной функции соответственно получаем импульсную характеристику фильтра:

h(k) ó H(z). (8.7)

Если функция H(z) представлена конечным степенным полиномом, что, как правило, характерно для НЦФ, то обратное z-преобразование осуществляется элементарно идентификацией коэффициентов по степеням z. Передаточная функция РЦФ также может быть представлена степенным полиномом прямым делением числителя на знаменатель правой части выражения (8.2), однако результат при этом может оказаться как конечным, так и бесконечным, т.е. система может иметь либо конечную, либо бесконечную импульсную характеристику. Фильтры с бесконечной импульсной характеристикой получили название БИХ-фильтров, с конечной импульсной характеристикой соответственно КИХ-фильтров. Нерекурсивные фильтры всегда являются КИХ-фильтрами, т.к. длительность импульсной реакции НЦФ определяется окном фильтра.

Примеры.

1. Передаточная функция РЦФ: H(z) = (1-z⁵)/(1-z).

Прямым делением числителя на знаменатель получаем: H(z) = 1+z+z²+z³+z⁴.

H(z) ó h(n) = {1,1,1,1,1}. Фильтр РЦФ является КИХ-фильтром.

2. Передаточная функция: H(z) = 1/(1-2z).

Методом обратного z-преобразования: h(n) = 2ⁿ.

Фильтр РЦФ является БИХ-фильтром.

Устойчивость фильтров. Фильтр называется устойчивым, если при любых начальных условиях реакция фильтра на любое ограниченное воздействие также ограничена. Критерием устойчивости фильтра является абсолютная сходимость отсчетов его импульсного отклика:

|h(n)| < ¥. (8.8)

Анализ устойчивости может быть проведен по передаточной функции. В устойчивой системе значение H(z) должно быть конечным во всех точках z-плоскости, где |z| £ 1, а, следовательно, передаточная функция не должна иметь особых точек (полюсов) на и внутри единичного круга на z-плоскости. Полюсы H(z) определяются корнями многочлена знаменателя передаточной функции (8.2).

Пример.

Передаточная функция фильтра рис. 6.1: H(z) = b₀/(1-a₁z). При а₁= 0.5 полюс знаменателя: z_р= 2. |z_р|>1. Фильтр устойчив.

Передаточная функция фильтра рис. 6.2: H(z) = b₀/(1+a₁z). При а₁= 1.1 полюс знаменателя: z_р= -0.909. |z_р| < 1. Фильтр неустойчив, что и подтверждает пример фильтрации.

Передаточная функция фильтра рис. 6.3: H(z) = 0.5(1+z)/(1-z). Полюс знаменателя: z_р= 1. В принципе, фильтр неустойчив, но эта неустойчивость проявляется только при k = ∞. Импульсный отклик фильтра h(n) = {0.5,1,1,1, ….}, сумма которого равна ∞ только при n = ∞, т.е. при интегрировании бесконечно больших массивов. При интегрировании конечных массивов результат всегда конечен.

Приведенный критерий устойчивости относится к несократимой дроби, т.к. в противном случае возможна компенсация полюса нулем передаточной функции, и следует проверить наличие однозначных нулей и полюсов.

Проверка на устойчивость требуется только для рекурсивных цифровых фильтров (систем с обратной связью), нерекурсивные системы всегда устойчивы.