Какая конструкция полосового фильтра даст самый короткий импульсный отклик?


14

Разработка простого полосового фильтра IIR II-го порядка по Бутеруорту с центральной частотой 500 Гц и шириной полосы в 1 октаву дает мне следующую частотную характеристику ...

введите описание изображения здесь

Теперь, если я возьму импульсный отклик и нормализую его и переведу в дБ, мы можем наблюдать затухание импульсного отклика.

введите описание изображения здесь

Затухание импульсного отклика приблизительно линейно со временем при построении графика на этой шкале, что позволяет нам определять статистику времени затухания (как в акустике помещения, где вы можете определить время реверберации). Чтобы импульсный отклик этого фильтра упал ниже 30 дБ, требуется около 11 мс.

Мы пытаемся минимизировать это время затухания, сохраняя следующую константу:

  • Полоса пропускания -3 дБ
  • Порядок фильтров

Я рад принять (в определенных пределах) полосу пропускания и пульсацию полосы пропускания и / или компромисс по крутизне полосы перехода для достижения этой цели. Кто-нибудь может предложить метод фильтрации с минимальной длительностью импульсного отклика, как определено выше?


2
Пожалуйста, включите частоту дискретизации, чтобы придать этим значениям некоторое значение.
Juancho

4
Поляки в фильтре будут давать экспоненциально затухающие члены в импульсной характеристике, которая, как вы показали, при построении графика в логарифмическом масштабе дает линейное затухание. Скорость распада связана с расстоянием полюсов до единичного круга; чем ближе они, тем медленнее затухание. Крутизна полосы перехода также связана с тем, насколько близко полюса находятся к единичному кругу. Я не знаю каких-либо методов проектирования, которые позволили бы вам расставить приоритеты именно этой характеристике.
Джейсон Р

@ Juancho Частота дискретизации была опущена, так как я думал, что это совершенно не имеет значения: использование 5 кГц или 500 кГц не меняет частоту затухания импульсной характеристики. Я нацеливаюсь на 44,1 кГц, если вам интересно. Спасибо за просмотр :)
learnvst

2
@JimClay да, я уверен, что вы можете, но я хочу, чтобы стоимость вычислений была очень низкой. Чтобы эффективно использовать Fir, мне нужно было бы использовать метод, основанный на fft, и это привело бы к задержке алгоритма, пока буфер FFT заполнен выборками. Да нет?
Learnvst

3
@JimClay, почему законы физики всегда мешают моим планам мирового господства! Ворчать ворчать
учись

Ответы:


8

В частотной характеристике фильтра Баттерворт в является результатом конкретных формул и его характеристикой является плоским частотной полосой пропускания отклика. Следовательно, если коэффициенты фильтра БИХ изменены каким-либо образом, фильтр может не поддерживать характеристики «Баттерворта».

В дополнение к ответам "Hilmar" и "Jason R", возможно, вы могли бы рассматривать это как проблему оптимизации фитнес-функции, которая фиксирует ваши спецификации.

Например, вы можете начать с некоторого дизайна (например, фильтра Баттерворта), а затем использовать технику оптимизации, чтобы переместить нули и полюса вокруг их местоположений (или изменить конструкцию путем добавления / удаления полюсов и нулей), пытаясь достичь ваших спецификаций (более четко) спад во временной области, поддержание пропускной способности и порядка фильтрации).

В этом направлении была проделана большая работа по созданию фильтров с использованием генетических алгоритмовздесь ) и имитации отжигаздесь ), которые могут оказаться полезными.


6

Боюсь, волшебной пули нет. Вы можете использовать эллиптический фильтр для независимого управления пульсациями в полосе пропускания и затуханием в полосе пропускания, однако вы обнаружите, что скорость затухания тесно связана с крутизной и общей полосой пропускания фильтра. Вы можете значительно уменьшить скорость фильтра, уменьшив порядок фильтров до 1, но с другой стороны, фильтр будет намного менее крутым.


Невозможно уменьшить фильтр до первого порядка, так как он используется в сети типа Linkwitz-Riley, но спасибо за то, что
нашли
Используя наш сайт, вы подтверждаете, что прочитали и поняли нашу Политику в отношении файлов cookie и Политику конфиденциальности.
Licensed under cc by-sa 3.0 with attribution required.