FFT не для всех частот? Опции

[ataradov]

Есть задача посторить спектр от сигнала (комлексного) оцифрованного с частотой 10 МГц (длинна выборки - 4096-32768 отсчетов). Но от этого спектра нужно только +/- 10 кГц. Но эти 10 кГц нужно иметь с достаточно большим разрешением (нужно хотя-бы 512 дискретов в частотной области).

Эта задача легко решается с помощью неполного DFT, но уж больно долго считается.

Я если честно не очень понимаю алгоритмы БПФ. Скажите пожалуйста можно-ли как-то аналогично БПФ ускорить вычисление вышеизложенной задачи или нет? Может предложите что-то свое?

И если это поможет - в дальнейшем нужно будет вычислять это преобразование непрерывно бегущим окном, длинна корого указана выше.

Спасибо.

[vladr]

Есть задача посторить спектр от сигнала (комлексного) оцифрованного с частотой 10 МГц (длинна выборки - 4096-32768 отсчетов). Но от этого спектра нужно только +/- 10 кГц. Но эти 10 кГц нужно иметь с достаточно большим разрешением (нужно хотя-бы 512 дискретов в частотной области).

Эта задача легко решается с помощью неполного DFT, но уж больно долго считается.

Я если честно не очень понимаю алгоритмы БПФ. Скажите пожалуйста можно-ли как-то аналогично БПФ ускорить вычисление вышеизложенной задачи или нет? Может предложите что-то свое?

И если это поможет - в дальнейшем нужно будет вычислять это преобразование непрерывно бегущим окном, длинна корого указана выше.

Спасибо.

Надеюсь, что Вам поможет http://www.exstrom.com/journal/sigproc/specmag.pdf.
Статья с С кодом: "A C program for spectrum magnification: When an FFT is not enough". Обещают эффективность, гораздо лучшую, чем вставление нулей.