[an error occurred while processing this directive]
|
Можно экономию довольно существенную получить, если вспомнить (из алгоритма БПФ с прореж. по частоте по основ. 4), что
X[4k]=FFT(x[n]+x[n+N/4]+x[n+N/2]+x[n+3N/4]).
FFT размером N/4. Вычисленное любым способом. n,k - 0-N/4-1. x - сигнал, умноженный на окно, длительностью N. X[4k] - каждый 4-ый спектральный отсчет:-)
Эквивалентно тому, что Вы хотите. Но Фурье в 4 раза меньшего размера.
А вообще вся процедура эквив. фильтрации гребенкой полосовых, комплексных в общем случае, КИХ фильтров с последующим прореживанием во времени. В пределе - без временного прореживания. Имп. х-ка каждого фильтра в гребенке=ОКНО(m)*exp(j*2pi*m*k/N). Можно учесть и центрирование на середину окна.
Отличаются ИХ только частотой комплексной экспоненты,k. Она дает центральную частоту фильтра. Все подавительные св-ва определяются окном. Окно по сути ИХ НЧ фильтра, с такими же полосами, подавлениями, неравномерностями, как и полосовой, но с нулевой центральной частотой.
Соответственно не обязательно использовать класс. окна. Можно попробвать сгенерировать НЧ фильтр с требуемой полосой проп., неравномерностями, подавлениями, перех. полосой с помощью алг. Ремеза, например. И использовать ИХ этого фильтра в качестве окна.
Хотя, в силу принципа, лежащего в алгоритме Ремеза, полученное окно будет неуловимо похоже, скорее всего, на окно Дольфа-Чебышева:-)
Но все-таки.
Если Вам нужна только небольшая часть сп. отсчетов, то взгляд через фильтры поможет еще большую экономию получить, выплывет Герцель:-)
Ну и необязательно использовать КИХи. Можно и БИХи. Хотя ну их:-).
E-mail: info@telesys.ru