БИХ-фильтр

БИХ-фильтр (фильтр с бесконечной импульсной характеристикой, англ. Infinite Impulse Response, IIR) — это тип цифрового фильтра, особенностью которого является наличие обратной связи: выходной сигнал фильтра зависит не только от входного сигнала, но и от предыдущих значений выходного сигнала. Благодаря обратной связи импульсная характеристика такого фильтра теоретически бесконечна (не затухает до нуля за конечное время), что и дало название классу. БИХ-фильтры широко применяются в цифровой обработке сигналов (ЦОС) для частотной селекции, подавления помех, коррекции сигналов и синтеза звука.

История

Развитие теории БИХ-фильтров связано с переходом от аналоговой схемотехники к цифровым методам обработки сигналов во второй половине XX века. Основы цифровой фильтрации были заложены в работах Р. Хэмминга, А. Оппенгейма, Р. Шафера и других исследователей. Первоначально для синтеза БИХ-фильтров использовались методы, основанные на преобразовании аналоговых фильтров-прототипов (например, фильтров Баттерворта, Чебышёва, Кауэра) в цифровую форму. Это стало возможным благодаря разработке инвариантного преобразования импульсной характеристики и билинейного преобразования. В 1960–1970-х годах с развитием микропроцессорной техники и появлением специализированных цифровых сигнальных процессоров (DSP) БИХ-фильтры стали активно внедряться в системы связи, аудиотехнику, радиолокацию и измерительные приборы.

Отличия от КИХ-фильтров

БИХ-фильтры противопоставляются фильтрам с конечной импульсной характеристикой (КИХ, англ. FIR). Основные различия:

Принцип работы и математическое описание

Работа БИХ-фильтра описывается разностным уравнением:

\[ y[n] = \sum_{k=0}^{M} b_k x[n-k] - \sum_{l=1}^{N} a_l y[n-l] \]

где:

• \(x[n]\) — входной сигнал,
• \(y[n]\) — выходной сигнал,
• \(b_k\) — коэффициенты прямой связи (нерекурсивная часть),
• \(a_l\) — коэффициенты обратной связи (рекурсивная часть),
• \(M\) и \(N\) — порядки нерекурсивной и рекурсивной частей соответственно.

Передаточная функция в z-области:

\[ H(z) = \frac{\sum_{k=0}^{M} b_k z^{-k}}{1 + \sum_{l=1}^{N} a_l z^{-l}} \]

Знаменатель передаточной функции определяет наличие полюсов, которые при неправильном выборе коэффициентов могут вывести фильтр за пределы устойчивости (попасть на единичную окружность или за её границы). Для устойчивости все полюсы должны находиться строго внутри единичного круга на z-плоскости.

Классификация и типы

БИХ-фильтры классифицируются по нескольким признакам.

По типу частотной характеристики

• Фильтры нижних частот (ФНЧ) — пропускают низкие частоты, подавляют высокие.
• Фильтры верхних частот (ФВЧ) — пропускают высокие, подавляют низкие.
• Полосовые фильтры (ПФ) — пропускают заданный диапазон частот.
• Режекторные (заграждающие) фильтры — подавляют узкий диапазон частот.
• Гребёнчатые фильтры — имеют периодическую АЧХ; используются в аудиоэффектах.

По структуре реализации

• Прямая форма I — простая, но чувствительна к ошибкам квантования.
• Прямая форма II (каноническая) — использует минимальное количество линий задержки.
• Транспонированная форма — модификация прямой формы, улучшающая численные свойства.
• Каскадная (последовательная) форма — разложение передаточной функции на биквадратные звенья (секции второго порядка). Наиболее распространена на практике, так как обеспечивает лучшую устойчивость и меньшую чувствительность к шумам квантования.
• Параллельная форма — разложение на сумму простых дробей.

По типу аппроксимации

• Фильтры Баттерворта — максимально гладкая АЧХ в полосе пропускания, без пульсаций.
• Фильтры Чебышёва I рода — равноволновая (пульсирующая) АЧХ в полосе пропускания, более крутой спад.
• Фильтры Чебышёва II рода — равноволновая АЧХ в полосе задерживания.
• Эллиптические (фильтры Кауэра) — равноволновая АЧХ и в полосе пропускания, и в полосе задерживания; обеспечивают максимальную крутизну спада при заданном порядке.
• Фильтры Бесселя — максимально линейная фазовая характеристика, но меньшая крутизна спада.

Методы синтеза

Билинейное преобразование

Наиболее распространённый метод. Аналоговый фильтр-прототип (например, Баттерворта) преобразуется в цифровой с помощью замены:

\[ s = \frac{2}{T} \cdot \frac{1 - z^{-1}}{1 + z^{-1}} \]

где \(T\) — период дискретизации. Метод сохраняет устойчивость и отображает всю ось частот в единичный круг, но вносит нелинейное искажение шкалы частот (frequency warping), которое компенсируется предварительным искажением (prewarping) частоты среза.

Инвариантное преобразование импульсной характеристики

Импульсная характеристика аналогового фильтра дискретизируется и используется как коэффициенты цифрового фильтра. Метод прост, но может приводить к наложению спектров (aliasing) и не гарантирует сохранения формы АЧХ.

Метод согласованного z-преобразования

Полюсы и нули аналогового фильтра отображаются на z-плоскость с помощью экспоненциальной замены \(z = e^{sT}\). Применяется реже из-за сложности и возможных проблем с устойчивостью.

Применение

БИХ-фильтры используются в областях, где важна низкая вычислительная сложность и не требуется строго линейная фаза:

• Аудиообработка — эквалайзеры, кроссоверы, фильтры низких частот для сабвуферов, эффекты (реверберация, хорус, фленжер).
• Связь — формирование спектра сигналов, подавление внеполосных помех, согласованные фильтры.
• Биомедицинская инженерия — фильтрация ЭКГ, ЭЭГ, пульсоксиметрии для удаления шумов и артефактов.
• Радиолокация и гидролокация — обработка сигналов, подавление помех, выделение целей.
• Управление — цифровые регуляторы (ПИД-регуляторы), фильтры в системах автоматического регулирования.
• Измерительная техника — сглаживание данных, подавление высокочастотных помех в датчиках.

Преимущества и недостатки

Преимущества

• Низкий порядок фильтра (обычно 2–10) для достижения высокой крутизны АЧХ.
• Меньшее количество операций умножения и сложения на отсчёт по сравнению с КИХ-фильтром аналогичной избирательности.
• Экономия памяти и вычислительных ресурсов, что критично для встраиваемых систем.

Недостатки

• Нелинейная фазовая характеристика, что может искажать форму сигнала (особенно для импульсных и переходных процессов).
• Потенциальная неустойчивость при неправильном выборе коэффициентов или при ошибках квантования.
• Чувствительность к шумам квантования коэффициентов и округлению промежуточных результатов, что может приводить к паразитным колебаниям (limit cycles).
• Сложность проектирования для фильтров с линейной фазой.

Реализация

БИХ-фильтры реализуются как программно (на микроконтроллерах, DSP, CPU), так и аппаратно (в FPGA, ASIC). В программной реализации часто используются библиотеки ЦОС, такие как CMSIS-DSP (ARM), Intel IPP, библиотеки для MATLAB, GNU Octave, SciPy (Python). Для уменьшения ошибок квантования применяются:

• Использование чисел с плавающей запятой (float/double) вместо фиксированной точки.
• Каскадная реализация в виде биквадратных звеньев.
• Увеличение разрядности промежуточных накоплений (аккумуляторов).
• Использование двойной точности на критических этапах.

Интересные факты

• БИХ-фильтры являются цифровым аналогом аналоговых фильтров с обратной связью (например, фильтров на операционных усилителях).
• Самый простой БИХ-фильтр — фильтр первого порядка, описываемый одним коэффициентом обратной связи и одним коэффициентом прямой связи. Он используется, например, в простых сглаживающих фильтрах (экспоненциальное скользящее среднее).
• В некоторых приложениях, где требуется минимальная задержка (например, в слуховых аппаратах), БИХ-фильтры предпочтительнее КИХ-фильтров из-за меньшей вычислительной задержки.
• Эллиптические фильтры (Кауэра) обеспечивают самую крутую переходную полосу среди всех БИХ-фильтров при заданном порядке, но имеют наибольшую неравномерность АЧХ.

Источники

• Оппенгейм А. В., Шафер Р. В. Цифровая обработка сигналов. — М.: Техносфера, 2006.
• Хэмминг Р. В. Цифровые фильтры. — М.: Советское радио, 1980.
• Лайонс Р. Цифровая обработка сигналов. — М.: Бином-Пресс, 2006.
• Proakis J. G., Manolakis D. G. Digital Signal Processing. — 4th ed. — Pearson, 2006.
• Smith S. W. The Scientist and Engineer's Guide to Digital Signal Processing. — California Technical Publishing, 1997.