Из временной области в комплексную плоскость
Z-преобразование отображает последовательность x_n в функцию X(z) комплексной переменной:
X(z) = Σ_{n=−∞}^{∞} x_n · z^{−n}
Переменная z параметризует комплексную плоскость. Различные области этой плоскости соответствуют различным качественным поведениям фильтра.
Геометрические области
| Область | z | Поведение | ||
|---|---|---|---|---|
| Внутри единичного круга | < 1 | Устойчивые полюсы: затухающий отклик | ||
| Единичный круг | = 1 | Частотная ось: z = e^{i2πf} | ||
| Вне единичного круга | > 1 | Неустойчивые полюсы: растущий отклик |
Единичный круг играет в устойчивости дискретного времени ту же роль, что мнимая ось в устойчивости непрерывного времени (Лапласа).
Связь с преобразованием Лапласа
Для систем непрерывного времени преобразование Лапласа использует переменную s. Мнимая ось s = iω — это место, где живет частотная характеристика. Устойчивость: полюсы должны иметь Re(s) < 0 (левая полуплоскость).
Билинейное преобразование отображает s → z: z = (1 + Ts/2) / (1 − Ts/2). Это отображает левую полуплоскость во внутреннюю часть единичного круга — геометрический перевод 'устойчивой левой полуплоскости' в 'устойчивый внутренний круг'.
Билинейное преобразование как конформное отображение
Билинейное преобразование z = (1 + Ts/2) / (1 − Ts/2) является преобразованием Мёбиуса — конформным (сохраняющим углы) отображением комплексной плоскости.
Его ключевые геометрические свойства:
- Отображает s = iω (мнимую ось) на |z| = 1 (единичный круг)
- Отображает Re(s) < 0 (левую полуплоскость) на |z| < 1 (внутреннюю часть единичного круга)
- Отображает Re(s) > 0 (правую полуплоскость) на |z| > 1 (внешнюю часть единичного круга)
- Искажение частоты: отображение ω → f нелинейно — ω_analog = (2/T)·tan(πf_digital)
Это искажение сжимает высокие частоты к точке Найквиста. Проектировщики учитывают это путем предварительного искажения аналоговой спецификации перед применением билинейного преобразования.
Полюсы Баттерворта: расположение на окружности
Фильтры Баттерворта достигают максимально плоской полосы пропускания путем размещения аналоговых полюсов на окружности радиуса ω_c в s-плоскости.
Для фильтра Баттерворта N-го порядка полюсы расположены в:
s_k = ω_c · e^{iπ(2k+N−1)/(2N)} для k = 0, 1, …, N−1
Это размещает их равномерно на левой половине окружности радиуса ω_c. (Полюсы на правой половине были бы неустойчивыми; сохраняются только полюсы в левой полуплоскости.)
Почему расположение на окружности → максимально плоская полоса пропускания?
Полином Баттерворта |B_N(jω)|² = 1 + (ω/ω_c)^{2N} имеет все свои полюсы при |s| = ω_c. Ограничение равного радиуса означает, что все полюсы одинаково способствуют амплитудной характеристике при ω = ω_c. Теорема о максимальной плоскости: среди всех полиномов N-го порядка с полюсами на этой окружности полином Баттерворта имеет наибольшее количество производных, равных нулю при ω = 0.
Полюсы Чебышева: расположение на эллипсе
Полюсы Чебышева расположены на эллипсе в s-плоскости (не на окружности). Эллипс имеет большую и малую полуоси, определяемые параметром пульсаций ε. Полоса пропускания с равными пульсациями возникает из свойства равноколебательности полиномов Чебышева.
Полюсы эллиптических фильтров: расположение на эллиптической функции
Полюсы эллиптических (Кауэра) фильтров также расположены на эллипсе — но с ПОЛЮСАМИ И НУЛЯМИ, способствующими частотной характеристике. Нули расположены на мнимой оси (полюсы конечного ослабления в полосе задержания). Отображение эллиптической функции оптимально распределяет нули для достижения равных пульсаций в обеих полосах одновременно.
Вычисление расположения полюсов Баттерворта
Для фильтра Баттерворта 4-го порядка с ω_c = 1 (нормализованное), полюсы расположены в:
s_k = e^{iπ(2k+3)/8} для k = 0, 1, 2, 3
k=0: s₀ = e^{i3π/8} (в левой полуплоскости)
k=1: s₁ = e^{i5π/8} (в левой полуплоскости)
k=2: s₂ = e^{i7π/8} (в левой полуплоскости)
k=3: s₃ = e^{i9π/8} (в левой полуплоскости)
Эти четыре полюса расположены с равным угловым интервалом на единичной окружности, все с отрицательными действительными частями (левая полуплоскость).
Расстояние от полюсов до единичного круга
Теоретическая устойчивость требует |p| < 1. На практике возникают два дополнительных условия.
Запас устойчивости
Запас устойчивости IIR-фильтра — это минимальное расстояние от любого полюса до единичного круга: min_k (1 − |p_k|).
Полюс при |p| = 0.99 технически устойчив, но оставляет только 1% запаса. Арифметика с конечной точностью (округление в представлении коэффициентов & накопление ошибок округления) может эффективно передвигать полюсы. Если квантование коэффициентов сдвинет полюс с 0.99 на 1.001, фильтр становится неустойчивым.
Геометрическое следствие
Полюсы, очень близкие к единичному кругу, производят очень острые пики частотной характеристики — узкополосные резонаторы. Но узкополосные резонаторы требуют высокую точность: небольшие ошибки коэффициентов значительно сдвигают частоту пика.
Геометрический компромисс: острота пика ∝ 1 / (1 − |p|). При |p| → 1, острота → ∞, но запас устойчивости → 0 & чувствительность к ошибкам коэффициентов → ∞.
Секции второго порядка
Высокопорядковый IIR-фильтр, реализованный как один полином, численно чувствителен — округление одного коэффициента может передвинуть много полюсов. Стандартное решение: реализовать как каскад секций второго порядка (биквадратные звенья), каждая с только одной парой сопряженных полюсов & одной парой сопряженных нулей. Ошибки в одном разделе не могут возмутить полюсы в других.