Цифровой фильтр в звуковом тракте это правило, по которому новый выходной отсчёт собирается из прошлых входных отсчётов и, в рекурсивном случае, из прошлых выходных. За этой простой формулировкой скрываются два принципиально разных семейства: фильтры с конечной импульсной характеристикой FIR и фильтры с бесконечной импульсной характеристикой IIR. Инженер, проектирующий эквалайзер, кроссовер для активной акустики или коррекцию помещения, обязан понимать, чем эти семейства различаются по памяти, устойчивости, фазовому поведению и вычислительной цене. Ниже разобраны математические основы обоих подходов, практические числа для типичной частоты дискретизации 48 кГц и критерии выбора структуры под конкретную задачу обработки сигнала.
Свёртка и линия задержки как физический смысл FIR фильтра
FIR фильтр описывается прямой свёрткой входного сигнала с набором коэффициентов: y[n] = Σ b_k x[n−k], где суммирование идёт по k от нуля до N−1. Каждый коэффициент b_k умножает входной отсчёт, задержанный ровно на k тактов дискретизации. Аппаратно это реализуется линией задержки, то есть сдвиговым буфером, который хранит N последних значений входа. Линия задержки и есть память фильтра: она буквально помнит прошлое сигнала длиной N отсчётов, и ровно эта память определяет, насколько низкие частоты и крутые переходные полосы фильтр способен различить. Импульсная характеристика совпадает с самими коэффициентами и заканчивается через N отсчётов, отсюда название - конечная импульсная характеристика.
Цена FIR следует из теории окон: чтобы получить узкую переходную полосу, нужна длинная импульсная характеристика, потому что ширина главного лепестка спектра обратно пропорциональна его длине. Пример: линейно-фазовый FIR фильтр нижних частот длиной 1024 отсчёта на 48 кГц. Групповая задержка такого фильтра равна половине длины, то есть 512 отсчётам, что при 48000 отсчётах в секунду составляет примерно 10,7 мс. Задержка одинакова для всех частот, форма волны не искажается, но в живом мониторинге десять миллисекунд уже ощутимы. Вычислительная нагрузка составляет 1024 умножения с накоплением на каждый выходной отсчёт, поэтому длинные FIR часто считают через быструю свёртку в частотной области: спектр блока умножается на спектр импульсной характеристики, сложность падает до O(log N) на отсчёт ценой блочной обработки с перекрытием и сохранением хвостов.
Хранение тоже ощутимо: тысяча коэффициентов в одинарной точности это четыре килобайта на фильтр, и в многоканальном устройстве таблицы превращаются в заметную статью памяти плюс буфер на 1024 входных отсчёта. Тем не менее предсказуемость окупает всё: поведение полностью описывается импульсной характеристикой, скрытых состояний нет, и смена таблицы коэффициентов не вызывает сюрпризов в динамике тракта.
Разностное уравнение с обратной связью и природа IIR фильтра
IIR фильтр добавляет к прямой ветви рекурсию: y[n] = Σ b_k x[n−k] − Σ a_k y[n−k], где коэффициенты a_k образуют обратную связь по прошлым выходам. Именно она делает импульсную характеристику бесконечной: единичный импульс запускает затухающую экспоненциальную реакцию, которая теоретически не заканчивается никогда. Фильтр помнит прошлое не через длинный буфер, а через внутреннее состояние из нескольких элементов, что и объясняет его экономность по вычислениям и памяти.
Z-преобразование переводит уравнение в передаточную функцию H(z) = B(z)/A(z) = (Σ b_k z^−k)/(1 + Σ a_k z^−k). Знаменатель порождает полюса, числитель нули, и карта полюсов на комплексной z-плоскости полностью определяет поведение. Условие устойчивости предельно ясное: все полюса должны лежать строго внутри единичного круга |z| = 1. Полюс на окружности даёт незатухающие колебания, а полюс снаружи превращает малейшую погрешность в растущий выброс, упирающийся в насыщение арифметики и превращающий тракт в генератор.
Практическая сила IIR в наследовании теории аналоговых фильтров. Аппроксимации Butterworth с максимально плоской характеристикой, Chebyshev первого и второго рода с пульсациями в полосе пропускания либо подавления и эллиптические фильтры с пульсациями в обеих полосах переносятся в цифру билинейным преобразованием. Фильтр четвёртого порядка на срезе 100 Гц при 48 кГц, для которого FIR потребовал бы сотен отводов, реализуется парой звеньев второго порядка, то есть примерно десятком коэффициентов. Разница в вычислительной цене достигает двух порядков, и поэтому почти весь параметрический эквалайзер построен на рекурсивных звеньях.
Обратная сторона - чувствительность к квантованию. Когда коэффициенты a_k округляются до сетки float или fixed point, полюса смещаются; для узких резонансных фильтров с полюсами у самой окружности сдвиг на доли процента выбрасывает полюс наружу, и фильтр становится неустойчивым уже после расчёта. Поэтому коэффициенты пересчитывают в двойной точности, а округление в рабочую сетку сопровождают проверкой радиуса полюсов и контрольным моделированием отклика.
Почему FIR всегда устойчив и даёт строго линейную фазу
У FIR передаточная функция H(z) = Σ b_k z^−k имеет только нули и полюса в нуле, не влияющие на устойчивость. Обратной связи нет, состояние конечно, любой ограниченный вход даёт ограниченный выход при любых конечных коэффициентах, включая жёстко квантованные. Инженер может округлять коэффициенты до шестнадцати бит, и худшее последствие - лёгкая деформация характеристики и рост уровня стоп-полосы, но никогда самовозбуждение. Это делает FIR безопасным выбором для автоматической коррекции по измерению помещения, где коэффициенты вычисляются алгоритмом и никто не проверяет их глазами.
Второе достоинство - точная линейная фаза. Если коэффициенты симметричны, b_k = b_{N−1−k}, либо антисимметричны, фазовая характеристика становится строго линейной функцией частоты, а групповая задержка константой (N−1)/2 отсчётов. Все спектральные компоненты задерживаются одинаково, огибающая транзиентов сохраняет форму, атака ударных не размазывается асимметрично. Симметрия к тому же экономит умножители: парные отводы складываются перед умножением, и число умножений сокращается почти вдвое.
У причинного IIR линейной фазы не бывает в принципе: обратная связь закручивает фазу тем сильнее, чем ближе частота к срезу и чем выше крутизна. Нулевую фазу даёт лишь офлайн двухпроходная фильтрация вперёд и назад. Тонкость есть и у FIR: линейно-фазовые фильтры делятся на четыре типа по чётности длины и знаку симметрии, и, например, антисимметричный фильтр чётной длины имеет принудительный нуль на частоте Найквиста и не годится для полной полосы пропускания, поэтому проектировщик выбирает тип заранее под требуемую форму отклика.
Цена крутого среза и сравнение вычислительных бюджетов
Задача: разделительный фильтр для двухполосной активной системы, срез 2 кГц, дискретизация 48 кГц, спад 48 дБ на октаву. IIR решение - каскад из четырёх biquad звеньев по топологии Linkwitz Riley, около двадцати умножений с накоплением на отсчёт. Линейно-фазовый FIR той же крутизны потребует двух-трёх сотен отводов, а при срезе 200 Гц счёт пойдёт на тысячи. Ширина переходной полосы FIR обратно пропорциональна длине: Δf примерно равна требуемому подавлению в децибелах, умноженному на частоту дискретизации и делённому на 22 и на N, поэтому низкие частоты среза - самый дорогой случай для свёртки и одновременно самый дешёвый для рекурсии.
Сравнивать надо не только умножения. FIR даёт постоянную задержку (N−1)/2, и для 1024 отсчётов это упомянутые 10,7 мс, тогда как у IIR задержка частотнозависима, мала в середине полосы и растёт у среза локальными горбами в десятки отсчётов. Для мониторного тракта музыканта средняя задержка IIR в единицы микросекунд решает всё; для сведения и кроссовера, где каналы суммируются, фазовая линейность FIR важнее. Гибридный приём - многоступенчатая децимация: дешёвые фильтры на высокой частоте снимают основную задачу, а крутой срез выполняется на пониженной частоте, где требуемая длина свёртки падает пропорционально коэффициенту децимации, и суммарный бюджет вычислений сокращается в разы.
Biquad звенья, порядок фильтра и численная точность
IIR высоких порядков никогда не реализуют одним уравнением: уже выше четвёртого порядка прямая форма разваливается от округления. Стандарт индустрии - каскад звеньев второго порядка, biquad, каждое с пятью коэффициентами: b0, b1, b2 для прямой ветви и a1, a2 для обратной при нормированном a0 равном единице. Пять коэффициентов и два элемента задержки на звено - и из таких кирпичиков собираются полки, пики, режекторы и кроссоверы любой сложности. Порядок полного фильтра равен удвоенному числу каскадов, полюса размещаются парами внутри единичного круга под контролем каждого звена отдельно.
Структура влияет на шум не меньше коэффициентов. Прямая форма I требует четырёх ячеек задержки и прилично ведёт себя в плавающей точке; транспонированная прямая форма II хранит два регистра состояния и является фактическим стандартом аудио библиотек; прямая форма II концентрирует ошибки округления в одной точке и проигрывает для сигналов с высоким динамическим диапазоном. В одинарной точности ошибка накапливается в регистрах состояния резонансных звеньев, а на срезах ниже нескольких сотен герц, где полюса жмутся к точке z = 1, теряются значащие биты разностей коэффициентов. Поэтому студийные плагины исполняют biquad в double, а расчёт коэффициентов ведут в двойной точности и квантуют в рабочий формат в последнюю очередь.
Практические правила укладываются в небольшой перечень:
- для одинарной точности и низких срезов применять транспонированную форму II с состоянием в double;
- каскадировать фильтры выше четвёртого порядка только через biquad секции с сортировкой по добротности;
- после квантования коэффициентов проверять, что полюса остались внутри единичного круга;
- в fixed point выбирать Q-формат с запасом знаков и насыщением в аккумуляторе;
- плавно интерполировать коэффициенты при изменении параметров, чтобы исключить щелчки от скачка состояния.
Групповая задержка, эквалайзеры и кроссоверы как поле применения
Групповая задержка есть минус производная фазы по круговой частоте и показывает, насколько разные компоненты проходят тракт за разное время. У линейно-фазового FIR она постоянна и равна (N−1)/2 отсчётов, что делает его эталоном для измерительных систем и суммируемых каналов. У IIR групповая задержка имеет горбы у срезов и резонансов: пиковый фильтр с добротностью 3 локально задерживает звук на миллисекунды, и на слух это проявляется как лёгкое смазывание атаки в обработанной полосе, хотя на низких частотах эффект маскируется естественным поведением самого помещения.
Параметрический эквалайзер почти всегда IIR: ручки частоты усиления и добротности требуют мгновенного пересчёта нескольких коэффициентов, минимальной задержки и копеечной стоимости вычислений, а фазовая кривизна даже благоприятна, потому что соответствует минимально-фазовому характеру самих акустических резонансов. Графический эквалайзер на тридцать одну полосу - это просто тридцать одно biquad звено в каскаде или параллельной сумме. Кроссовер иллюстрирует компромисс лучше всего: Linkwitz Riley четвёртого порядка даёт плоскую сумму полос за считаные такты на отсчёт, а линейно-фазовый FIR убирает фазовый поворот и ровнее сшивает полосы по направленности, но платит задержкой 5–10 мс, которую приходится выравнивать в соседних ветвях отдельной линией задержки. На срезах ниже 100 Гц длина FIR раздувается ещё сильнее, и процессоры делают гибрид: сабвуферная часть на IIR, основная полоса на FIR коррекции. Точность тоже разделяет лагеря: FIR спокойно живёт в одинарной точности, а biquad кроссовера на 40 Гц при 192 кГц в float уже граничит с потерей устойчивости и требует double.
Инженерная логика проста: где критичны фаза, надёжность и произвольная форма коррекции, берут FIR и платят длиной и задержкой; где критичны стоимость вычислений и мгновенная реакция на ручки управления, берут каскады biquad и платят фазовой кривизной и заботой о полюсах. Профессиональный тракт почти всегда содержит обе технологии, и мастерство разработчика состоит в том, чтобы каждой задаче назначить ту структуру, чьи слабости в данном месте тракта не заняты.
Отдельной строкой стоит культура проверки результатов. Перед выпуском кроссовера в продукт инженер свершает три вещи: смотрит АЧХ суммы полос с реальными драйверами, слушает вертикальную направленность на оси между излучателями и обязательно контролирует импульсный отклик на ступеньку - именно она выдаёт как фазовые горбы IIR-цепочек, так и колокольный «пре-рингинг» длинного FIR у среза. Если система проверена этими тремя способами и прошла, она переживёт и полевую работу; кажущиеся академическими различия структур на сцене оборачиваются вполне материальной картиной: устойчивый вокал в центре, уверенный бас без гудения и минимум усталости у слушателя через час работы.