[an error occurred while processing this directive]
|
если все вычисления в целых числах
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 частот. Плёвое дело ;-)
E-mail: info@telesys.ru