БПФ для оценки низкочастотной составляющей Опции

[roman73]

Доброго времени суток!

Необходимо выделить низкочастотную составляющую временного ряда.
Алгоритм по ТЗ такой:
1.Делается БПФ.
2.Обнуляются все частоты, кроме некоторого числа низких.
3.Делается обратное БПФ.

Вопросы:
Если применять классическое БПФ, что делать с мнимой частью сигнала?
Допустим, мнимая часть входного сигнала устанавливается в 0.
А что делать с мнимой частью результирующего сигнала?
Или можно применить другое преобразование вместо Фурье, работающее только с действительным сигналом?
Косинус-преобразование? Хартли? Насколько я знаю, преобразование Хартли напрямую не даёт частоты, а что-то другое.
А как с этим у косинус-преобразования?
Применимо ли тут RealFFT?
Какой из перечисленных алгоритмов лучше для данного случая?
Где взять эффективные реализации?

Число отсчетов ~50000-100000

Сообщение отредактировал roman73 - Aug 22 2011, 14:55

[V_G]

1. Во входном сигнале мнимой компоненты нет.
2. Вычисленное ДПФ из такого сигнала автоматом имеет симметричную действительную и антисимметричную мнимую части.
3. Если при обнулении частей спектра вы сохраните свойства симметрии действительной и мнимой частей, ОДПФ автоматом даст в результате нулевую мнимую часть.
4. Откуда уверенность, что после применения предложенного в начале алгоритма на выходе будет строго периодический сигнал (не учитывая естественной периодичности, связанной с длиной выборки)?

ЗЫ. Что Вы имели в виду, когда писали, что фильтр - программный? Все цифровые фильтры вычислительные (и в этом смысле программные) по определению.

Сообщение отредактировал V_G - Aug 22 2011, 14:56

[roman73]

4. Откуда уверенность, что после применения предложенного в начале алгоритма на выходе будет строго периодический сигнал (не учитывая естественной периодичности, связанной с длиной выборки)?

Нет такой уверенности.
Так в ТЗ написано.
Переименовал тему.

ЗЫ. Что Вы имели в виду, когда писали, что фильтр - программный? Все цифровые фильтры программные по определению.

Ну теоретически цифровой фильтр тоже можно схемно реализовать
Имел в виду что интересует реализация на обычном компе с виндой.
Не интересует реализация на микроконтроллерах/плисах и т.д.