Как достичь резких полос перехода фильтра с помощью методов проектирования фильтров Iir
Фильтры с бесконечным импульсным откликом (IIR) являются краеугольным камнем современной цифровой обработки сигналов, ценятся за их вычислительную эффективность и способность реализовывать крутые частотные отклики с небольшими коэффициентами. Достижение четкой полосы перехода фильтра - узкой области между полосой пропускания и полосой остановки - является критическим требованием в таких приложениях, как сети аудио-трансфера, выравнивание каналов связи, уравнивание биомедицинского сигнала и формирование импульсов радара. В то время как фильтры с конечным импульсным откликом (FIR) предлагают линейную фазу и присущую стабильность, фильтры IIR могут соответствовать или превышать характеристики отката конструкций FIR со значительно более низким порядком фильтра, что делает их незаменимыми для систем реального времени и ограниченных ресурсами. Эта статья обеспечивает всестороннее исследование методов проектирования фильтров IIR с бритвенно-острыми переходными полосами, охватывающих классические аналитические подходы, передовые структурные реализации и практические соображения стабильности.
Основы фильтров IIR и переходных полос
Определение полосы перехода
Переходная полоса фильтра — это диапазон частот, по которому отклик величины падает от края полосы пропускания до края полосы остановки. В идеальных фильтрах эта полоса имеет нулевую ширину, но для практических реализаций требуется конечная переходная область. Крутизна скатывания количественно определяется порядком фильтра и выбранной функцией приближения. Для фильтров IIR ширина полосы перехода обратно пропорциональна порядку фильтра, но удвоение порядка не обязательно вдвое сокращает полосу перехода — нелинейные эффекты возникают из нулевых взаимодействий полюса и рекурсивной природы фильтра.
Почему IIR Excels в Шарп Катоффс
Фильтры IIR достигают резких переходов с помощью обратной связи — выход зависит как от входных, так и от прошлых выходов. Эта обратная связь позволяет функции передачи фильтра размещать полюса вблизи единичного круга, создавая высокий коэффициент качества (Q) , который производит быстрое падение величины вокруг частоты среза. Один участок IIR второго порядка может реализовать резонансный пик, который потребует десятков кранов FIR. Следовательно, для данной спецификации полосы перехода фильтр IIR обычно требует порядка в пять-десять раз ниже, чем его аналог FIR, переводя на меньшее умножение на образец и более низкое использование памяти.
Ключевые методы проектирования для полос перехода Sharp
Фильтры IIR высшего порядка
Увеличение порядка фильтров является наиболее простым способом погружения в полосу перехода. Фильтр Баттерворта порядка N имеет откат 20NdB/decade. Однако прямое осознание фильтра IIR высокого порядка как единой передаточной функции может привести к сильной численной чувствительности и нестабильности. Две общие стратегии реализации смягчают эти проблемы:
Каскадные секции второго порядка (SOS)
Путем факторизации полинома высокого порядка в продукт секций второго порядка (биквады), каждый со своей парой сложносопряженных полюсов, фильтр становится гораздо более надежным против квантования коэффициента. Каждый биквад может быть спроектирован с его полюсами хорошо в пределах единичного круга, а общий ответ является продуктом отдельных секций. Это также упрощает настройку и отладку. Для резких переходов полюса каждой секции расположены ближе к единичному кругу, поэтому тщательное упорядочивание и спаривание полюсов с нулевым полюсом необходимы для предотвращения перелива и ограничения колебаний цикла. Стандартная практика заключается в том, чтобы упорядочить секции за счет увеличения частоты полюсов и спаривания ближайших полюсов и нулей, чтобы минимизировать пиковое усиление.
Параллельное внедрение формы
Альтернативой каскадной является параллельная форма, где фильтр представлен как сумма секций второго порядка (расширение частичной фракции). Такой подход снижает эффект ошибок коэффициента, поскольку каждый раздел вносит самостоятельный вклад в общий отклик. Параллельные реализации особенно привлекательны, когда фильтр имеет широкий динамический диапазон, поскольку они избегают кумулятивных ошибок усиления, присущих каскадным структурам. Однако параллельная форма может быть более вычислительно дорогостоящей для очень высоких порядков, поскольку число секций растет менее эффективно.
Эллиптические (Cauer) фильтры
Эллиптические фильтры являются чемпионами крутых переходных полос для данного порядка. Они достигают самой узкой ширины перехода, позволяя эквирипплу как в полосе пропускания, так и в полосе остановки. Стоп-пол имеет конечный пол затухания, к которому фильтр неоднократно приближается по мере увеличения частоты. Конструкция основана на якобианской эллиптической функции, и полученный фильтр имеет отклик величины, который соответствует заданной полосе пропускания (δ ]p ) и затуханию полосы пропускания (δ s ) с минимальным возможным порядком фильтра.
Продажа и продажа Stopband Ripple
Резкость переходной полосы эллиптического фильтра напрямую зависит от величины пульсации. Увеличение допустимой пульсации полосы пропускания позволяет фильтру подтолкнуть полосу остановки ближе к краю полосы пропускания, сужая область перехода. Аналогично, ослабление затухания полосы остановки позволяет еще более резкое начальное падение. Например, эллиптический фильтр с пульсацией полосы пропускания 0,1 дБ и затуханием полосы остановки 60 дБ может достичь коэффициента перехода (край полосы остановки / край полосы пропускания) 1,1, тогда как фильтр Баттерворта потребует порядка почти в три раза выше для того же соотношения.
Чебышев Тип I и Тип II Фильтры
Фильтры типа I Чебышева имеют равномерную полосу пропускания и монотонную полосу остановки, а фильтры типа II (обратный Чебышев) имеют монотонную полосу пропускания и равномерную полосу остановки. Оба предлагают более резкие переходы, чем Butterworth для того же порядка. Тип I помещает полюса в эллипс, который охватывает круг Баттерворта, позволяя более крутой откат вблизи отсечки. Тип II достигает аналогичной селективности, вводя нули на воображаемой оси для создания выемок, которые быстро снижают коэффициент усиления полосы пропускания. Для приложений, требующих нулевой пульсации полосы пропускания (например, овладение звуком), тип II предпочтительнее, хотя его нулевые полосы остановки могут производить фазовые искажения, которые должны быть оценены.
Butterworth и Bessel Filters — оригинальное название
Фильтры Баттерворта максимально плоские в полосе пропускания, но имеют самый медленный откат для данного порядка — 20N dB/decade. Фильтры Бесселя предназначены для постоянной задержки группы, жертвуя наклоном отката для линейной фазы. Ни один из них не является оптимальным для резких переходных полос, но оба являются надежными отправными точками, когда требуется стабильная фазовая реакция или минимальный перерасход. Для резких переходов дизайнеры обычно избегают Баттерворта, если порядок не может быть сделан очень высоким, и вычислительные ресурсы многочисленны. Фильтры Бесселя почти никогда не выбираются для крутых отсечек; вместо этого они используются в приложениях временной области, где сохранение импульса имеет решающее значение.
Использование билинейного преобразования с предварением
При проектировании фильтров IIR из аналоговых прототипов (наиболее распространенный подход) билинейное преобразование отображает переменную частоты непрерывного времени в дискретно-временной домен z. Билинарное преобразование сжимает бесконечную ось аналоговой частоты в конечную полосу от 0 до π, вызывая деформацию частоты. Критические частоты (полосы полосы пропускания и стоп-диапазона) должны быть предварительно деформированы до преобразования, чтобы гарантировать, что цифровой фильтр отвечает желаемым местам среза. Для резких полос перехода даже небольшие ошибки предварительного искажения могут сместить срез достаточно, чтобы нарушить спецификации, поэтому формула предварительного искажения Ωaa = (2/T) tan(ωd/2]
Расширенные возможности для Sharp Transition Bands
Размещение полюс-ноль и чувствительность
Резкие переходные полосы требуют полюсов, чрезвычайно близких к единичному кругу. Например, узкий фильтр нижних частот с отсечкой 0,1π и шириной полосы перехода 0,01π может иметь радиус полюса 0,998 или выше. Такая близость к единичному кругу делает фильтр очень чувствительным к квантованию коэффициента — ошибка округления даже 2—16 может вытолкнуть полюс за пределы единичного круга, вызывая нестабильность. Поэтому дизайнеры должны использовать двухточную арифметику во время фазы проектирования и рассмотреть масштабирование коэффициента и обратную связь ошибок (например, структуры дельта-сигма) в реализациях с фиксированной точкой. Управление квантованием коэффициента является критическим навыком для достижения надежных резких переходных полос.
Численные эффекты стабильности и квантирования
Помимо размещения полюсов, фактическая реализация фильтра IIR высокого Q страдает от предельных циклов (колебания из-за округления) и переливающихся колебаний. Структуры, оптимизированные для минимизации этих эффектов, включают форму Латтиса-Лестницы и фильтр State-Variable (Chamberlin). Перенесенная структура II прямой формы часто рекомендуется для ее более низкой чувствительности, но она все еще может накапливать шум. Для очень резких фильтров исследователи предложили использовать волновые цифровые фильтры (WDF) , которые моделируют физические пассивные сети и проявляют отличные числовые свойства даже тогда, когда полюса чрезвычайно близки к единичному кругу. Другой подход - структура сцепленной формы , которая изолирует косинусные и синусовые термины для снижения чувствительности.
Алгоритмы оптимизации для индивидуализированных переходных полос
Когда классические конструкции (Баттерворт, Чебышев, Эллиптик) не отвечают определенным нестандартным ограничениям, таким как полоса перехода, которая должна быть резкой только в определенном диапазоне частот или требование минимальной вариации задержки группы, могут быть использованы методы численной оптимизации. Итеративные перевзвешенные наименее квадратные (IRLS) и генетические алгоритмы могут регулировать нулевые местоположения полюсов, чтобы минимизировать максимальное отклонение между фактическим ответом величины и идеальной формой кирпичной стены. Эти методы часто начинаются с эллиптического дизайна и тонкой настройки коэффициентов. Однако оптимизация должна включать ограничения стабильности (поли внутри единичного круга) и часто добавляет штраф за радиус полюсов, превышающий 0,999, чтобы избежать чрезмерной хрупкости числа.
Практические шаги дизайна для фильтров Sharp IIR
Определение спецификации
Начните с указания полосы пропускания (δ]p в dB), затухания полосы остановки (δs в dB), частоты края полосы пропускания (ωp) и частоты края полосы остановки (ωs. Ширина полосы перехода Δω = |ωsp. Для данного типа фильтра эти спецификации непосредственно определяют минимальный порядок. Используя стандартные формулы или таблицы проектирования (например, ]ellipord в MATLAB Signal Processing Toolbox или в SciPy), вычислите порядок и естественные частоты, которые отвечают требованиям
Оценка фильтра заказов
Порядок N для эллиптического фильтра может быть аппроксимирован известной формулой с использованием коэффициента селективности k = ωp/ωs и отношения стоп-диапазона к пульсирующим энергиям.Как правило, каждое удвоение порядка фильтра уменьшает ширину полосы перехода вдвое при использовании эллиптических фильтров. Для типа I Чебышева порядок примерно . Мудро моделировать полученный порядок и проверять запас; если полоса перехода предельно слишком широка, увеличить порядок на единицу и повторно имитировать.
Моделирование и проверка
После проектирования коэффициентов фильтра, выполнить анализ частотной реакции с использованием FFT высокого разрешения (например, 8192 точки) для измерения фактической ширины перехода. Проверьте шаг отклика для пересечения и звонка - острые IIR фильтры часто проявляют значительные артефакты домена времени из-за сохраненной энергии в их петлях обратной связи. Используйте частотный размах, чтобы визуально подтвердить, что частотная реакция попадает в заданную полосу пропускания ряби и достигает необходимого затухания полосы остановки на краю полосы остановки. Для реализации в реальном времени запустите фильтр на тренажере DSP с фиксированной арифметикой для оценки воздействия квантования коэффициента и шума округления.
Сравнение с FIR для резких переходов
Никакое обсуждение резких переходов IIR не является полным без признания альтернативы: FIR-фильтры со многими кранами. FIR-фильтр линейной фазы порядка M имеет полосу пропускания перехода, пропорциональную 1/M (приблизительно 4π/M для конструкции окна Kaiser) (приблизительно 4π/]М. Для ширины перехода 0,01π FIR может достигать того же самого. Для IIR-фильтра требуется около 24 коэффициентов (6 биквадов × 4 коэффициента каждый) против 400 для FIR — экономия более 15× в памяти и 10× в многократно накапливаемых операциях. Однако FIR-фильтры могут быть реализованы с эффективными многофазными структурами для многократных приложений, и они гарантируют стабильность и линейную фазу. Для приложений, где нелинейная фаза приемлема, IIR-фильтры являются окончательным выбором для резки
Заключение
Проектирование IIR-фильтров с резкими переходными полосами является одновременно искусством и наукой, требующей глубокого понимания классических аналоговых приближений, дискретного отображения времени и численной стабильности. Эллиптические фильтры предлагают самый крутой скачок для данного порядка, в то время как чебышевские и Баттерворт служат конкретным компромиссам между плоскостью и селективностью. Фильтры высокого порядка требуют надежных структур реализации, таких как каскадные секции второго порядка и тщательное сопряжение полюсов-ноль. Передовые методы, такие как предварение, оптимизация коэффициентов и волновые цифровые фильтры, раздвигают границы достижимой селективности в практических системах. Следуя систематическим шагам проектирования, изложенным в этой статье - определение спецификации, оценка порядка, моделирование и проверка - инженеры могут уверенно создавать IIR-фильтры, которые отвечают самым требовательным требованиям переходной полосы, сохраняя стабильность и вычислительную эффективность.