Теория систем автоматического регулирования

       

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


Рабочие файлы: [Аналоговые прототипы дискретных фильтров]

Примеры ЧХ дискретного фильтра: [z2_tf_fr.vsm][z2_tf_fr.mcd]

Передаточную функцию (ПФ) для дискретного фильтра любого типа можно получить из ПФ аналогового прототипа посредствам Билинейного преобразования (БЛП). Одновременно БЛП может быть использовано для масштабирующих операций в частотном домене. В частности для смещения полосы пропускания и регулировки ее ширины.

Дискретная ПФ фильтра второго порядка имеет вид:

(1)

b0 + b1*z^-1 + b2*z^-2 H(z) = ------------------------ a0 + a1*z^-1 + a2*z^-2

Для сокращения нагрузки на ЦВМ шесть коэффициентов z-ПФ могут быть заменены эквивалентными, один из которых будет равен единице (одна из шести операций умножения будет сокращена). Если проектируемая z-ПФ лишь одна из нескольких включенных последовательно, то можно привести к единице еще один коэффициент, за счет коррекции общего коэффициента передачи цепочки фильтров или системы. Тогда z-ПФ будет иметь вид:

(2)

(b0/a0) + (b1/a0)*z^-1 + (b2/a0)*z^-2 H(z) = --------------------------------------- 1 + (a1/a0)*z^-1 + (a2/a0)*z^-2

или

(3)

1 + (b1/b0)*z^-1 + (b2/b0)*z^-2 H(z) = (b0/a0) * --------------------------------- 1 + (a1/a0)*z^-1 + (a2/a0)*z^-2

Дальнейшее сокращение операций умножения возможно, если z-ПФ имеет пары равных коэффициентов.

Рекуррентный алгоритм, реализующий z-ПФ (2), должен рассчитывать следующее разностное уравнение (РУ):

(4)

y[n] = (b0/a0)*x[n] + (b1/a0)*x[n-1] + (b2/a0)*x[n-2] - (a1/a0)*y[n-1] - (a2/a0)*y[n-2]



Содержание раздела