Узкополосное преобразование Фурье

[sgk]

Задача вроде несложная, но с заковыкой.
Есть звук один канал оцифрованный на 44Кгц.
Нужно посмотреть спектр с хорошим частотным разрешением в полосе предположим от 500 до 530Гц....

Если посмотреть в буквальном смысле, то Спектралаб с 1 М точек. Будет 24 точки на Герц.
Пример на картинке, можно взять набор данных для дальнейшей обработки, сравнить с тем что у Вас получится после написания собственной программы.
Сергей sgk

[expert_elk]

Если посмотреть в буквальном смысле, то Спектралаб с 1 М точек. Будет 24 точки на Герц.
Пример на картинке, можно взять набор данных для дальнейшей обработки, сравнить с тем что у Вас получится после написания собственной программы.
Сергей sgk

Похоже тема уже ушла в небытие. Однако интересно, решилась ли проблема и каким образом?

[Леонид3]

expert_elk, Виталий, эта "проблема" решается просто: коэффициенты разложения по Фурье для любой частоты считается как сумма произведений дискретных значений выборки на соответствующие значения SIN (COS) -- синусы можно (и нужно) вычислить заранее. Если длина выборки кратна четырём, то число умножений можно сократить в четыре раза, а при длине равной 2^n - могут встретиться частоты, на которых число умножений меньше и в 8, 16 ... раз -- обычная алгебра.

[UR5ZQV]

Вот такая фигня получилась:
-Тупо устанавливаю частоту дискретизации 1500 Гц.
-выборка 1024 (2048).
-оконная функция (лучше всего получилось с Блэкманом-Харрисом).
-БПФ.
Разрешение ок.7Гц (4 Гц).
Но считает зараза ок. 0.5 сек. (1 сек).

P.S. Relayer, Если можно алгоритм децимации (или ссылочку).