Структура обратной связи
Фильтр КИХ (конечная импульсная характеристика) вычисляет каждый выходной отсчет как взвешенную сумму текущих и прошлых входов только. Без обратной связи. Импульсная характеристика имеет конечную длительность.
Фильтр БИХ (бесконечная импульсная характеристика) подает предыдущие выходы обратно в вычисление:
y_n = Σ b_k · x_{n−k} − Σ a_k · y_{n−k}
Рекурсивный член −Σ a_k · y_{n−k} создает обратную связь. Один импульс на входе будет циркулировать по циклу обратной связи бесконечно долго (геометрически затухая, если система устойчива).
Зачем использовать обратную связь?
Фильтр БИХ может достичь резкой частотной селективности с гораздо меньшим числом коэффициентов, чем фильтр КИХ. 2-полюсный БИХ может обеспечить такое же ослабление в полосе затухания, какое требует 50-коэффициентный КИХ.
Цена: потенциальная неустойчивость. Полюсы H(z) определяют устойчивость. Все полюсы должны лежать строго внутри единичного круга.
История Хэмминга о душе и обратной связи
Хэмминг использовал яркую личную историю для иллюстрации неустойчивости обратной связи.
Он всегда останавливался в одной и той же комнате отеля, потому что знакомство помогало ему ориентироваться, когда он был усталым. Сантехник установил трубы горячей воды большого диаметра в душе. Они создали значительную задержку между регулировкой крана и ощущением изменения температуры воды.
Каждое утро Хэмминг следовал одной и той же схеме: вода слишком холодная → добавить горячей → все еще холодная → добавить еще → вдруг кипящая → прыгнуть → убавить → повторить.
Задержка в цикле обратной связи означала, что его коррекции всегда приводили к перерегулированию. Он не мог адаптироваться к задержке, даже после многих повторений.
Урок инженерии: неустойчивость возникает либо из-за чрезмерного усиления в цикле обратной связи, либо из-за чрезмерной задержки в цикле обратной связи. Оба проявляются как одно и то же колебательное поведение. В терминах фильтров: полюсы на единичном круге или вне его производят именно такой колебательный или расходящийся отклик.
Характеристика неустойчивости
Хэмминг заметил, что одну и ту же неустойчивость в душе можно проанализировать двумя способами:
1. Его ответ был слишком сильным (чрезмерное усиление в действии коррекции).
2. Его обнаружение было слишком задержано (он поспешно входил в ванну, не дав системе установиться).
Оба описания дают один и тот же математический результат: полюс цикла обратной связи переместился вне единичного круга.
Четыре классических семейства
Теория аналоговых фильтров развивалась вокруг четырех классических семейств проектирования, каждое представляющее различный компромисс. Эти семейства преобразуются в дискретное время через билинейное преобразование или инвариантность импульса.
Баттерворт (максимально плоский)
Отклик в полосе пропускания: |H(jω)|² = 1 / (1 + (ω/ω_c)^{2N}). Монотонно убывающий. Без ряби в полосе пропускания или полосе затухания. Полюсы лежат на круге радиуса ω_c в s-плоскости (или преобразованном круге в z-плоскости). Самая плоская возможная полоса пропускания для заданного порядка N.
Чебышев тип I
Равная рябь в полосе пропускания, монотонный в полосе затухания. Для заданного порядка N и уровня ряби достигает более резкого срезания, чем Баттерворт. Полюсы лежат на эллипсе (в s-плоскости).
Чебышев тип II
Равная рябь в полосе затухания, монотонный в полосе пропускания. Зеркальное отражение типа I в частотной области.
Эллиптический (Кауэр)
Равная рябь в ОБЕИХ полосе пропускания и полосе затухания. Для заданного порядка N и уровней ряби достигает наиболее резкого возможного перехода от полосы пропускания к полосе затухания. Использует эллиптические функции для оптимального размещения полюсов и нулей. Хэмминг: название происходит от того, что эллиптические функции используются в выводе.
Фундаментальный компромисс
Все четыре семейства достигают одного и того же базового компромисса по-разному: более высокий порядок N дает более резкий переход. Допущение ряби (Чебышев, эллиптический) достигает более резкого перехода для того же N. Эллиптический достигает абсолютно самого резкого перехода для любых заданных N и спецификаций ряби.
Выбор между семействами фильтров
Выбор между семействами зависит от того, что может допустить приложение.
Сомнение в требовании эксперта
Хэмминг вспомнил, что определенные эксперты утверждали, что все фильтры БИХ (рекурсивные) обладают определенным свойством. Он спросил себя, действительно ли это было правдой — и нашел контрпример.
Его точка: эксперты часто несут утверждения, которые они поглотили в школе, не переэкзаменуя их в контексте текущих проблем. Если вы спросите себя, действительно ли то, что вам говорят, является правдой, удивительно, насколько много вы можете найти, что это либо, либо граничит с ложью, даже в хорошо развитой области.
Контрпример не был типом фильтра, который вы обычно проектировали, но он доказывал, что утверждение было поверхностным. Одного контрпримера достаточно для опровержения универсального утверждения.
Проектирование БИХ на практике
Хэмминг отметил, что он независимо разработал большую часть теории фильтров БИХ, решая другую проблему: вывод стабильных формул корректора для численного решения обыкновенных дифференциальных уравнений.
Форма формулы корректора: y_n = Σ a_k · y_{n−k} + Σ b_k · f(y_{n−k})
Обратная связь появляется и в членах y (линейная обратная связь), и в членах f(y) (нелинейная обратная связь через дифференциальное уравнение). Устойчивость для фильтров БИХ является частным случаем более общей проблемы устойчивости для численных интеграторов ОДУ.
Связь обратной связи через различные области
Одна и та же математическая структура — обратная связь, полюсы, граница устойчивости — появляется в цифровых фильтрах, численных решателях ОДУ, системах управления, биологических ритмах и экономических моделях.
В каждой области: цикл обратной связи вычисляет новое состояние из предыдущих состояний. Устойчивость требует, чтобы обратная связь не усиливала возмущения бесконечно.
Граница устойчивости единичного круга в Z-плоскости соответствует: мнимой оси в плоскости Лапласа s (непрерывное время), условию спектрального радиуса ρ(A) < 1 для линейных итераций, и условию показателя Ляпунова λ < 0 для нелинейных систем.