快速傅里叶变换FFT算法及其应用_[当文网提供]..doc

快速傅里叶变换FFT算法及其应用 摘 要 本文较为系统地阐述了快速傅里叶变换的算法原理及其在数字信号处理等工程技术中的应用。根据抽取方法的不同，一维基2 FFT算法分为两种：频域抽取的FFT算法和时频域抽取的FFT算法。第1节阐述了这两种FFT算法的原理。第2节给出了两种算法的编程思想和步骤。第3节阐述了一维非基2 FFT的两种算法：Cooley-tukey FFT算法和素因子算法(Prime Factor Algorithm)的思想原理，给出了在把一维非基2 DFT的多层分解式转化为二层分解的过程中，如何综合运用这两种算法以达到总运算次数最少的方案；并以20点DFT为例描述了非基2 FFT算法实现的一般步骤。第4节介绍了一维FFT算法在计算连续时间信号的傅里叶变换、离散信号的线性卷积、离散信号压缩和滤波等数字信号处理中的典型应用。第5节把一维FFT变换推广到二维FFT变换，并在一维FFT算法的基础上，给出了二维FFT算法的原理和实现过程。最后在附录中给出了一维DFT的基2 FFT 算法（包括频域抽取的FFT和IFFT算法、时域抽取的FFT和IFFT算法），一维任意非基2 FFT算法，二维DFT的基2 FFT 算法以及二维DFT的任意非基2 FFT 算法的详细的Visual C++程序。 本文通过各种流程图和表格，较为深入系统地阐述了FFT的算法原理；运用Matlab编程，通过大量生动的实例，图文并茂地列举出了FFT算法的各种应用，并在每个实例中都附上了完整的Matlab程序，可供读者参考。由于篇幅所限，本文未涉及FFT变换以及其应用的数学理论背景知识。 关键词：FFT算法的应用，一维基2 FFT算法，频域抽取，时域抽取，非基2 FFT算法，Cooley-Tukey算法，素因子算法，线形卷积，信号压缩和滤波，二维FFT算法 目 录 摘 要 1 目 录 2 1 一维DFT的快速算法—FFT 4 1.1 频域抽取的基2算法 4 1.1.1 正变换的计算 4 1.1.2 逆变换的计算 7 1.2 时域抽取的基2算法 8 2 一维基2 FFT算法编程 10 3 一维任意非基2 FFT算法 14 3.1 Cooley-Tukey FFT算法 14 3.2 素因子算法（PFA） 14 3.3 一维任意非基2 FFT算法 16 4 一维FFT算法的应用 20 4.1 利用FFT计算连续时间信号的傅里叶变换 20 4.2 利用FFT计算离散信号的线性卷积 23 4.3 利用FFT进行离散信号压缩 25 4.4 利用FFT对离散信号进行滤波 28 4.5 利用FFT提取离散信号中的最强正弦分量 31 5 二维DFT的快速变换算法及应用简介 36 5.1 二维FFT变换及其算法介绍 36 5.2 二维FFT变换算法的应用 37 参考文献 38 附 录 39 1．一维DFT的基2 FFT 算法Visual C++程序 39 （1） 频域抽取的FFT和IFFT算法 39 （2） 时域抽取的FFT和IFFT算法 44 2．一维任意非基2 FFT算法Visual C++程序 49 3．二维DFT的基2 FFT 算法Visual C++程序 54 4．二维DFT的任意非基2 FFT 算法Visual C++程序 62 1 一维DFT的快速算法—FFT 当序列的点数不超过时，它的点定义为 （1） 反变换定义为 （2） 二者形式相似，快速算法的原理一样，这里先就其正变换进行讨论。令，当依次取为时，可表示为如下的方程组： （3） 由上式可见，直接按照定义计算点序列的点时，每行含个复乘和个加，从而直接按定义计算点的总计算量为个复乘和个加。当较大时，很大，计算量过大不仅耗时长，还会因字长有限而产生较大的误差，甚至造成计算结果不收敛。所谓快速傅里叶变换就是能大大减少计算量而完成全部点计算的算法。下面介绍两种经典的的快速算法：频域抽取的算法和时域抽取的算法。 1.1 频域抽取的基2算法 1.1.1 正变换的计算 这里仅介绍基2算法，即是2的整次幂的情况。由定义 （4） 把分成两半，即和，代入（4）式得 （5） 由于 （5）式两项又可合并为 （6） 当为偶数时，注意到，，（6）式变为 （7） 当为奇数时， ，（6）式变为