快速傅立叶变换FFT程序设计.docVIP

快速傅立叶变换（FFT）程序设计 DFT(Discrete Fourier Transformation)是数字信号分析与处理，如图形、语音及图像等领域的重要变换工具，直接计算DFT的计算量与变换区间长度N的平方成正比。当N较大时，因计算量太大，直接用DFT算法进行谱分析和信号的实时处理是不切实际的。快速傅立叶变换(Fast Fourier Transformation，简称FFT)使DFT运算效率提高1—2个数量级，其原因是当N较大时，对DFT进行了基4和基2分解运算。FFT算法利用DFT系数的对称性、周期性和可约性等性质将长序列的DFT分解为若干个短序列的DFT计算，然后再按一定规则将其合并，从而得到整个的DFT。按时间抽取的基2 FFT算法是FFT最常用的一种表现形式。 设（=2）点序列的离散傅立叶变换为，则有： (1) (2) 式中：=0,1，……,；=0,1,……,。 将(=0,1,……, )按的奇偶分成以下两组： 将上式带入(1)式并化简，就得到(3)和(4)式。这样N点序列就变为两个点序列和离散傅立叶变换的组合： (3) (4) 式中：=0,1,……, ， 其中 (3)式和(4)式可以用下面的碟形节信号流图图1来表示。 图1 2点碟形节 对于和可以继续利