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

Еще можно посмотреть в сторону интерполированных КИХ-фильтров.
За счет включения M-элементов задержек удается получить узкую полосу при гораздо меньшей ( в 5..7 ) раз длине фильтра.
А прототип рассчитывать "оконным" методом, да хоть того же Кайзера, что довольно быстро.

Ну так я собственно и хочу. Токо формулы лень выводить для разных типов фильтров. Все надеюсь что есть готовые На дворе начало ХХI века, а эта тема должна была бы быть исчерпаной, уже к концу 60х

Поставьте банк КИХ фильтров на 128, 256,...1024,...8192,... полосы и не мучайтесь. Сразу исчезнут проблемы с обеспечением устойчивости и не нужны будут никакие формулы. Банк фильтров легко реализуется через БПФ. БПФ на 8192 точки для сигналов звукового диапазона потянет наверно любой современный PC в реальном времени без всяких проблем.

Эта тема неисчерпаема, как и Мир.

Вам сюда за реализацией:
http://musicdsp.org/archive.php?classid=3#225

Найдите статью
Regalia P.A., Mitra S.K. and Vaidyanathan P.P. The Digital All-Pass Filter: A Versatile Signal Processing Building Block // Proc. IEEE. 1988.V.76. №1. Р.1937.
Она же есть в переводе в ТИИЭР.
Там есть полосовые фильтры с формулами настройки по частоте и ширине,
есть фильтр - эквалайзер.
Фильтры замечательные тем, что не возбуждаются при любой комбинации коэффициентов.

Ха-ха! Спасибо конечно, но я там давно был и неоднократно. И там такого нет. В смысле есть обрезанный код с МИТ-овской лицензией, но он вычисляет все прямым способом, для риалтайма не годящимся. А вот это уже слегка похоже на то, что надо: http://www.abvolt.com/research/publications2.htm, токо книгу "Digital Filters: Basics and Design"
Dietrich Schlichtharle, (народную) не могу нарыть. Насчет статьи "Regalia P.A., Mitra S.K. and Vaidyanathan P.P. The Digital All-Pass Filter: A Versatile Signal Processing Building Block", нашел - весьма интересна, но пока не совсем понятны механизмы реализации.

Какого "такого"?

Параметрические эквалайзеры 2-го порядка в обозримом пространстве используются 2-х видов:

1 Регалии-митры (менее популярны)
2 Бристоу-джонсона (наиболее популярны)

Различаются определением полосы.

Вам-то что нужно?

По ссылке
http://musicdsp.org/archive.php?classid=3#225
хотя бы Delphi Class implementation of the RBJ filters
чем не устраивает ?

Или Вы чего-то не договариваете

А где там выбрать х-ку Баттерворта, Чебышева, Бесселя?

Нигде. Баттерворт, чебышев 1/2 ,бессель, эллиптический это все способы аппроксимации чх стандартных типов, т е

фнч фвч пф рф.

Эквалайзер не относится к стандартному типу фильтров. Вам что надо?

Если честно, то после прохождения сигнала через узкополосные БИХ-фильтры, Вы маму родную не узнаете
Я за КИХ-варианты.

А мне надо стандартные типы фильтров, только и всего. "Эквалайзер не относится к стандартному типу фильтров. " - если их не всегда делают такими, то это еще ничего не значит... Иногда все же делают.

Тада эквалайзеры (по-научному - корректоры ачх - фильтр с произвольной ачх, синтезируется только методом оптимизации) тут при чем?



Дауш... (ц)

Фазовый фильтр имеет модуль передаточной функции H(z), равный |H(z)| = 1. Его ФЧХ на частоте fR имеет фазовый сдвиг, равный 180°.
Если складываются сигналы от двух фазовых фильтров, выходной сигнал подавляется на частотах, для которых разность фаз равна 180°.
Результирующая передаточная функция:
HS = (H1(z) ± H2(z))/2,
соответствует различным фильтрам ФНЧ,ФВЧ, полосовому или режекторному фильтру в зависимости знака суммирования (ФНЧ или ФВЧ) и от порядка фазовых фильтров.
Так, если H1(z) -фазовое звено 2-го порядка, а H2(z)=1, и вычесть их выходы, то на всех частотах, кроме fR, результат будет нулевой,
а на частоте fR - он усилится в 2 раза.

Передаточная функция фазового фильтра 2-го порядка равна:

H1(z)= (a+b(1+a)Z^-1 +Z^-2)/ (aZ^-2+b(1+a)Z^-1 +1),

где b = cos(2πfR) – регулирует частоту среза fR, a = (1-t)/(1+t) - задает ширину полосы пропускания Δf, t = tg(πΔf).
Таким образом, изменяя b , центральная частота пропускания регулируется в пределах (0,05 - 0,45) fsampling с подавлением в полосе задержания до 50 дб.

Аппаратная реализация фильтра на VHDL описана в книжке А.Сергиенко "VHDL для...".

Не превращайте тему во флейм. В названии абсолютно четко указано, что мне надо, и другие разные способы (КИХ, формулы из аудиокукбукс и тд ) меня не устраивают.