Отправлено SM 30 ноября 2004 г. 09:33

Суть сего в том, что каждый отсчет ПФ можно рассмотреть как цифровой фильтр. Как свертку входного сигнала с N периодами комплексной экспоненты. X(n) = sum(k=0, N-1, x(k)*exp(j*2*pi*(n/N)(N-k)). Слегка помудрив с математикой, получается H(z)=(1-Z^-N)/(1-exp(j*2*pi*n/N)*z^-1). Выбросив числитель (вместо этого можно просто перед каждлй пачкой из N отсчетов занулять состояние фильтра), и преобразовав знаменатель так, чтобы выкинуть оттуда все мнимое, получится H(z)=(1-exp(-j*2*pi*n/N)*z^-1)/(1-2*cos(2*pi*n/N)*z^-1+z^-2). Итого надо вычислить N раз рекурсивную часть фильтра, и потом один раз нерекурсивную, и получить интересующий отсчет ДПФ. Ели интересует всего несколько точек ДПФ, этот алгоритм эффективнее остальных.

P.S. Если уж хотите засунуть исходник в сообщение, то не забывайте, что знаки больше и меньше это спец-символы в HTML. А то Вам даже и не ответишь как положено

• Это тоже ответ для AlexeyV. — -=ВН=-   (30.11.2004 10:41, 2549 байт)
  • В апликатионах Аналога для 21хх семейства есть хорошее объяснение. — Valera_by   (30.11.2004 12:05, 51 байт)
    • Самое лучшее объяснение - в любом учебнике по ЦОС. Хотя оно мало отличается от поста ВН. — SM   (30.11.2004 12:18, пустое)