快速算法的研究.docVIP

目录 摘要………………………………………………………………..1 第一部分 FFT的回顾……………………………………………2 第一节 DFT与FFT的过渡…………………………………….2 第二节 按时间抽取以2为基的算法(FFT-DIT2)（库利-图基算法）………………………………………………………………...2 第三节 按频率抽取以2为基的算法(FFT-DIT2)（桑德-图基算法）...………………………………………………………………6 第二部分 混合基与分裂基的研究………………………………..7 第一节 混合基的相关知识………………………………………7 第二节 分裂基的相关知识………………………………………9 第三部分 FFT算法改进…….…………..……………………….12 第四部分 FFT实例运行结果及比较……………………………14 第五部分 总结与感想………………………………………….……18 参考文献 …..……………………………………………………...19 摘要 快速傅立叶变换FFT是离散傅立叶变换DFT的一种快速算法,我们知道有限长序列的主要特点是其频域也可以离散化成有限长序列，即能进行DFT变换。DFT广泛应用于信号的频谱分析、滤波器的设计以及系统的分析、设计和实现中。 而在相当长时间中由于N值大时，计算量大。无法达到实时处理的要求。1965年,库利(cooley)和图基(Tukey)首先提出FFT算法.对于N点DFT,仅需(N/2)log2N 次复数乘法运算。大大提高了计算效率。后来又有桑德快速算法以及一些算法改进。FFT是DFT的一种快速计算机算法，无任何近似，无精度损失。 本论文主要介绍基本FFT算法与实现，以及混合基、分裂基等一些先进算法，最后附有方案改进和算法程序设计。 第一部分 FFT的算法回顾 DFT与FFT的过渡 N点有限长序列x(n)的DFT为 二者仅有一符号之差和一个常熟因子1/N。由DFT可知完成整个DFT运算共需N2次复数乘法和N(N-1)次复数加法。由复数知识可知，此运算4N2次实数乘法和 2N(N-1)次实数加法运算。因此DFT 正比于N2，当N很大时，其运算量很可观。为此我们对此加以改进。 在DFT中我们知道WNnk具有一些特定性质： 可得出 由此可对DFT中的一些项合并，将其分解为短序列DFT，FFT就是由此而得出的。 第二节 按时间抽取以2为基的算法(FFT-DIT2)（库利-图基算法） 设序列点数N=2L，L∈Z。如条件不满足可人为加上若干个0，N=2L序列n的奇偶分两组 于是得 经过WNnk的性质转换使X（k）前后分离如下： 对X(k)的前一半（k=0,1, …,N/2-1）,由DFT定义: 对X(k)的后一半（k=N/2, N/2+1, … ,N-1） 由此我们只需求0——N/2－1区间的，，就可求出X（k）。其算法流图如下图一。 X1（0） x(0) X（0） X1（1） WN0 x(2) X（1） X1（2） WN1 x(4) X（2） X1（3） WN2 x(6) X（3）