Отправлено fontp 14 июля 2005 г. 12:23 В ответ на: обновлять результаты FFT надо с частотой оцифровки тогда... я так понимаю. получится беда :) отправлено <font color=gray>2b|!2b?..</font> 13 июля 2005 г. 22:06

если все вычисления в целых числах

1. По всем каналам сигнал умножается на табличный SIN и COS на каждом отсчёте
C(i) = X(i)*COS(W*i)
S(i) = X(i)*SIN(W*i)
Эти штуки отправляются в линии задержки

2. По всем каналам DFT вычисляется рекурентно, например

D1(n) = D1(n-1) + C(n) - C(n-M)
D2(n) = D2(n-1) + S(n) - S(n-M)

М - длина окна. Для DFT кратна периоду всех частот

Хотя рекурсия находится на границе устойчивости, с целыми числами такие штуки с интеграторами обычно проходят поскольку вычисления не содержат ошибок

Получается 96 умножений на отсчёт и 192 сложения для 48 частот. Плёвое дело ;-)

• понятно это ещё не совсем DFT - c точностью до накрученого фазового — fontp   (14.07.2005 12:26 80.64.97.142, 69 байт)