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

[ataradov]

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

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

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

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

Спасибо.

[Alex255]

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

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

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

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

Спасибо.

БПФ в Вашем случае поможет. ДПФ это 512*4096 условных отсчета, в БПФ это будет всего (4096/2)*log2(4096) - разность почти на 2 порядка. При этом на выходе в первом случае будет 512 отсчетов, а во втором - 4096. Никакой магии или потери информации здесь нет. Суть БПФ в том, что ядро преобразования (которая экспонента) циклично и одно его значение умножается на разные отсчеты. Остается лишь грамотно организовать процесс комбинирования.

Почитать о БПФ можно например здесь:
http://www.analog.spb.ru/

[ataradov]

Никакой магии или потери информации здесь нет.

Но дело в том, что нужно получить эти 4096 отсчетов на первые 10 кГц. А в случае простого БПФ эти отсчеты получатся на весь сигнал и расстояние между ними будет 10e6/4096 ~ 2.4 кГц. то-есть в интересующий диапазон попадет всего 5-6 отсчетов.

Производить децимацию нежелательно, так как децимационные фильтры приведут к потере части информации, а ее в этом случае и так почти нет.

Делается это все для определения допплеровской частоты в очень плохом принятом сигнале.

[Alex255]

Но дело в том, что нужно получить эти 4096 отсчетов на первые 10 кГц. А в случае простого БПФ эти отсчеты получатся на весь сигнал и расстояние между ними будет 10e6/4096 ~ 2.4 кГц. то-есть в интересующий диапазон попадет всего 5-6 отсчетов.

Производить децимацию нежелательно, так как децимационные фильтры приведут к потере части информации, а ее в этом случае и так почти нет.

Делается это все для определения допплеровской частоты в очень плохом принятом сигнале.

Шаг (разрешение) определяется длиной (по времени) вашего окна, по любому БПФ даст Вам все что можно вытянуть из выборки. Нужно лучшее разрешение - делайте более длинную выборку (по времени!), в конечном итоге частоту отсчетов можно существенно снизить, зачем 10М если интересует 10к. Децимация актуальна для устранения шумов - нужна ли она Вам?