Эффективная децимация в два раза Опции

[Alex65111]

Задача понизить частоту дискретизации в два раза. Соответственно надо сигнал подфильтровать, а потом децимировать. Сейчас я рассчитываю фильтр с большой крутизной, малыми (не более 0.1дБ) пульсациями в полосе пропускания, и подавлением не менее 96дБ методом Ремеза. В принципе работает. Но мне не дает покоя мысль, что где-то я встречал какой-то "хитрый" подход к решению проблемы децимации в два раза, в котором получается достаточно вычислительно эффективный фильтр. То ли малого порядка, то ли с кучей нулей в ИХ. Подскажите, на самом ли деле есть более лучшее решение?

[Alex65111]

Данные вещественные. Входной антиалисный аналоговый фильтр обеспечивает нужный уровень подавления примерно на 60кГц, хотя полезная полоса - до 30кГц. Соответственно частота дискретизации выбрана 120кГц. Далее необходимо цифровым фильтром подогнать полосу к 30кГц и продецимировав довести частоту дискретизации до 60кГц (опять же для вещественного сигнала).

Подход к проблеме в лоб в принципе задачу решает - Ремезом рассчитываю фильтр с переходной полосой от 29 до 30кГц, с пульсациями в полосе пропускания не более 0.1дБ и подавлением 96дБ. Фильтр получается где-то около 500 порядка.

Задача - попытаться найти решение не в лоб. По подсказке попробовал halfband.
В fdatool задал следующие параметры:
тип - halfband lowpass
filter order - 400
Fs=120000кГц
Fpass=29000кГц

Вроде у него каждый второй тап получился равен нулю. Да и пульсации и крутизна подходит.

При попытке применить опцию расчета децимирующего фильтра с коэффициентом децимации 2 и фильтром найквиста получился фильтр у которого порядок просто изумительный всего 47, но у него нужное подавление возникает только на 36кГц, т.е. у меня кусок от 24 до 30кГц будет искажен (а разрешается быть искажено только 29-30кГц) .

Так что пока я склоняюсь к halfband lowpass.

Кстати, я что-то не понял - в чем существенная разница между halfband lowpass расчетом фильтра и децимирующего фильтра Найквиста?

Сообщение отредактировал Alex65111 - Nov 11 2007, 15:06

[zhorro]

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