Фильтр ФНЧ, Переходная характеристика Опции

[DDR4]

Здравствуйте. Подскажите мне пожалуйста - можно ли создать цифровой ФНЧ (желательно FIR), со следующими характеристиками: 35Гц -6дБ, 50Гц -40 дБ. Неравномерность в частоте пропускания - 1дБ. Параметры сигнала - полезный сигнал в диапазоне 0.5 - 25Гц, количество выборок в секунду - 1900. Большое количество коэффициентов или большой объём вычислений - не проблема (т.к считаем на ПК). Важна неизменность формы сигнала в полезной полосе. Проблема заключается в чём - в плохой переходной характеристике. Т.е улучшаем частотную - получаем реакцию на ступеньку в виде выброса. Пробовал разные окна - но в итоге всё равно про нормальной частотной хар-ке имеем выбросы (Всё это в MatLab). А хочется гладкую переходную характеристику (без выбросов с единичного скачка) и ослабление на частоте задержки в -40(-50 ещё лучше) дБ (в тех пределах, которые возможны с учётом того, что считать машине можно много и долго ). Подскажите, возможно ли это?

Сообщение отредактировал DDR4 - Oct 19 2008, 11:57

[RadioJunior]

Вопрос: хочу програмным образом изменять у ФНЧ (БИХ или КИХ и порядком пока не определился) частоту среза и крутизну завала АЧХ (вид ФЧХ при этом не интересует). На входе только ц.сигнал. Можно ли обойтись без fft при расчете ИХ (вектора коэффициентов)? Может кто подскажет.

[Самурай]

Вопрос: хочу програмным образом изменять у ФНЧ (БИХ или КИХ и порядком пока не определился) частоту среза и крутизну завала АЧХ (вид ФЧХ при этом не интересует). На входе только ц.сигнал. Можно ли обойтись без fft при расчете ИХ (вектора коэффициентов)? Может кто подскажет.

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

Если уж хотите совсем по простому, и нет жестких требований к оптимальности и характеристикам (неравномерность и подавление) фильтра, тогда классический вариант: рисуете на бумажке нужную АЧХ и ручками, на той же бумажке, выводите аналитическое уравнение для коэфф. фильтра (через преобразование Фурье от АЧХ). В этом самом уравнении характеристики фильтра (частота среза и крутизна) задаются парой параметров. Добавить умножение коэфф. фильтра на оконную ф-ю (по вкусу) и вот вам Ваш фильтр но, правда, с не самыми выдающимися характеристиками.

Если не ошибся, то вот, например коэфф. для фильтра с трапециидальной (трапециевидной? ) АЧХ (формулу окна найдете сами):

h(n) = (cos(Fp*n) - cos(Fs*n))/[(Fs-Fp)*pi*n*n],

n = -m...+m - порядок фильтра,
Fp - частота среза, от 0 до pi,
Fs - частота подавления, от 0 до pi

В принципе, этот алгоритм не сильно отличается от расчета через FFT, разница только в стационарности. Что проще решайте сами.