DSP_Chapter5_快速傅立叶变换.pptVIP

Chapter 5The Fast Fourier Transform 罗新龙 信息工程学院 北京邮电大学 Overview DFT的计算 快速傅立叶变换算法 短时傅立叶变换 1. DFT的计算 由于时域处理更加经济, 因此滤波器设计已经转向时域处理. 但是频域处理使得有些问题变得更加简单: 需要有效的方法来计算DFT. 离散傅立叶变换 DFT: 离散序列的离散变换. 矩阵形式: 具有许多特殊结构?可以设计高效算法. 计算复杂性 对于每一个点需要N个复数乘法和N-1个复数加法, k = 0, 1,2, …, N-1. 复数乘法: (a+bj)(c+dj) = (ac-bd)+j(ad+bc) 需要4个实数乘法 + 2个实数加法. 复数加法 = 2 个实数加法. 总共需要4N2 个实数乘法和 4N2 –2 N 个实数加法. Goertzel 算法 把X(k)化为: 以上式子具有卷积的形式, 令 Goertzel 算法 因此把X(k)表示为卷积的形式: X(k) = yk(N) 而 yk(n) = xe(n)* hk(n) 多项式可以表示为 Goertzel 算法 对于每一个X(k)分别滤波处理: 即使只有一些采样点也能计算. 不需要很大的缓存和系数表. 与完全傅立叶变换有同样的复杂度 (4N2 个实数乘法和 4N2