|
Вычисление спектра выборки, когда выборка больше, чем требуемое разрешение |
|
|
|
Oct 20 2008, 14:36
|
Частый гость
 
Группа: Свой
Сообщений: 161
Регистрация: 19-01-08
Из: Питер
Пользователь №: 34 234

|
Допустим нужно нам получить спектр сигнала и допустим нужно сделать для этого именно N-точечное БПФ. Однако у нас есть выборка сигнала длинной kN, где k - целое число > 1. Тут имеется три варианта: 1) Взять любые подряд идущие N-точек из выборки и посчитать на них БПФ, но так мы теряем полезный сигнал, который можно было бы использовать. 2) Сделать k раз БПФ и результат усреднить. Вроде бы хорошо, но долго. 3) Каким-то образом усреднить выборку до N-точек и сосчитать один раз БПФ.
Собственно по 3-му пункту мне достался какой-то код, который не до конца понятно что делает. В чем я успел разобраться так это в том, что он умножает всю выборку на окно Кайзера длинной kN, затем ещё на какое-то окно, а затем хитрым образом складывает отсчеты сигнала. А именно так, если s - это исходная выборка умноженная на окна, а a - выборка, которая подается на БПФ, то: a_0 = s_0 + s_N + s_{2*N} + ... + s_{(k-1)*N} a_1 = s_1 + s_{N+1} + s_{2*N+} + ... + s_{(k-1)*N+1} ... a_{N-1} = s_{N-1} + s_{2*N-1} + ... + s_{k*N - 1}
Собственно, кто-нить пояснит что здесь за математика такая? Как это работает и что здесь делается.
|
|
|
|
|
 |
Ответов
|
Aug 22 2009, 06:00
|

Местный
  
Группа: Участник
Сообщений: 240
Регистрация: 20-09-08
Пользователь №: 40 347

|
Цитата(DMax @ Oct 20 2008, 18:36)  Допустим нужно нам получить спектр сигнала и допустим нужно сделать для этого именно N-точечное БПФ. Однако у нас есть выборка сигнала длинной kN, где k - целое число > 1. Тут имеется три варианта: 1) Взять любые подряд идущие N-точек из выборки и посчитать на них БПФ, но так мы теряем полезный сигнал, который можно было бы использовать. 2) Сделать k раз БПФ и результат усреднить. Вроде бы хорошо, но долго. 3) Каким-то образом усреднить выборку до N-точек и сосчитать один раз БПФ.
Собственно, кто-нить пояснит что здесь за математика такая? Как это работает и что здесь делается. Это называется полифазное БПФ (polyphase FFT). Читай что не ясно спрашивай.
|
|
|
|
|
Aug 23 2009, 14:24
|
Группа: Участник
Сообщений: 4
Регистрация: 9-08-09
Пользователь №: 51 803

|
Цитата(bahurin @ Aug 22 2009, 13:00)  Это называется полифазное БПФ (polyphase FFT). Читай что не ясно спрашивай. Не совсем понятно, что понимается под разрешением. Расстояние между пинами в обоих случаях одинаковое. Насколько я понимаю, окна нужно аккуратно подбирать в обоих вариантах, обычное БПФ тоже потеряло гармонику в последнем примере. Ну и в эксперименты шуму наверно стоило добавить, тогда бы преимущества полифазного БПФ были бы заметны.
|
|
|
|
Сообщений в этой теме
DMax Вычисление спектра выборки Oct 20 2008, 14:36 petrov Читайте про банки фильтров. а_0...a_{N-1} - это вы... Oct 21 2008, 20:03 ivan219 А можно по подробнее??? Nov 23 2008, 21:02 ivan219 Up Aug 21 2009, 19:06 _ea_ Над s выполняется БПФ размера kn, находится только... Aug 22 2009, 01:53 petrov Цитата(bahurin @ Aug 22 2009, 10:00) Алго... Aug 22 2009, 08:53  bahurin Цитата(petrov @ Aug 22 2009, 12:53) Это н... Aug 22 2009, 09:45   petrov Цитата(bahurin @ Aug 22 2009, 13:45) Все ... Aug 22 2009, 10:20    bahurin ЦитатаPerfect reconstruction банки фильтров для лю... Aug 23 2009, 09:59     petrov Цитата(bahurin @ Aug 23 2009, 13:59) Поли... Aug 23 2009, 10:15   Xenia Цитата(bahurin @ Aug 22 2009, 12:45) ТАК ... Aug 22 2009, 12:34  bahurin Цитата(_ea_ @ Aug 23 2009, 18:24) Не совс... Aug 23 2009, 16:30 Xenia Есть предложение получить суммы (или среднее) для ... Aug 22 2009, 09:19 ivan219 А причём тут частота Найквиста? Нужно к примеру то... Aug 23 2009, 09:59 bahurin Цитата(ivan219 @ Aug 23 2009, 13:59) А пр... Aug 23 2009, 10:14
1 чел. читают эту тему (гостей: 1, скрытых пользователей: 0)
Пользователей: 0
|
|
|