перестраиваемый полосно-заграждающий фильтр. Опции

[_sv_]

Добрый день всем.

Сразу: железо реализовано на дохлом 32х разрядном арм.
Поэтому ищу алгоритмы с целочисленной арифметикой.
Спектр сигнала мне уже известен (в программе уже реализован анализатор спектра).

Моя проблема стара как мир. Возбуждение в звуковом тракте. Микрофон-Динамик.
Что хочется попробовать: перестраиваемый полосно-заграждающий фильтр.
Частоту возбуждения - знаю. Полосу подавления - тоже. Степень подавления - тоже могу высчитать.

Проблемы:
1) Что-то не нахожу реализаций полосно-заграждающих фильтров.
Неужели надо лепить ФНЧ+ФВЧ?

2) Что-то не нахожу реализаций с целочисленной арифметикой.
Неужели таковых нет?

3) Вообще мало реализаций перестраиваемых полосно-заградительных фильтров.
Где искать?

Добавил:
Данные 16 бит, 16кгц дискретизация.

Сообщение отредактировал _sv_ - Sep 21 2014, 07:28

[Skryppy]

Можете заранее расчитать коэффиценты заграждающих фильтров в FDATools пакета MATLAB,там,кстати можете экспортировать фильтр в C header. "С Coder" матлаба также позволит вам реализовать различные фильтры, причем вы сами настраиваете арифметику-с плавающей или фиксированной точкой. А про перестраиваемые фильтры обсуждение было в соседней теме "Цифровой узкополосный фильтр с динамически изменяемой центральной частотой".

Сообщение отредактировал Skryppy - Sep 23 2014, 14:12

[_sv_]

Можете заранее расчитать коэффиценты заграждающих фильтров

Да. Могу. Заранее рассчитать - не проблема. Но мне надо чуточку иное.

Существует ли способ перевода коэфицинетов "плавающей точки" в фиксированную?
Т.е. посчитать коэфицинеты в варианте "плавающей точки" (а1, а2, в0, в1 и в2) получается.
Вроде бы получается. По крайней мере - фильтр делает то, что я ожидаю.
Но перевод в фиксированную точку пока не удался.

[des00]

Но перевод в фиксированную точку пока не удался.

round(pipa*2^N) не работает, а потом проверить на устойчивость ?