4 第四节 快速傅里叶变换.pptVIP

第四章　快速傅里叶变换CHARPTER 4 FAST FOURIER TRANSFORM 1.　引言 2.　基－2 DIT FFT算法 3. 基－2 DIF FFT算法 4.　IDFT快速计算方法 5.　进一步减少FFT计算量的方法 6. 分裂基FFT算法 7.　线性卷积与线性相关的FFT算法 第一节 引言 FFT的提出 DFT最大的缺陷：计算复杂度太大 Fast Fourier Transform (FFT) 是一种高效地计算DFT的方法。 军事上，工程上都需要这种高效性。 J. W. Cooley and J. W. Tukey, An algorithm for machine calculation of complex Fourier series,1965 FFT开始大规模应用于数字信号处理，数值分析等领域 For example: In 1969, the 2048 point analysis of a seismic trace took 13 ? hours. Using the FFT, the same task on the same machine took 2.4 seconds! 关于Gauss FFT的正式提出 直到1965年，Cooley和Tukey在《计算机科学 》发表著名的《机器计算傅立叶级数的一种算法》论文，FFT才开始大规模应用。 那个年代，有个肯尼迪总统科学咨询委员会。其中有项研究主题是，对苏联核测试进行检测，Tukey就是其中一员。 美国/苏联核测试提案的批准，主要取决于不实地访问核测试设施而做出检测的方法的发展。 其中一个想法是，分析离海岸的地震计情况，这种计算需要快速算法来计算DFT。其它应用是国家安全，如用声学探测远距离的核潜艇。 所以在军事上，迫切需要一种快速的傅立叶变换算法，这也促进了FFT的正式提出。 FFT的继续改进 桑德（G.Sand）、图基等快速算法相继出现，几经改进，很快形成现在的快速傅立叶变换（FFT）。这种算法使DFT的运算效率提高1到2个数量级，为数字信号处理技术应用于各种信号的实时处理创造了良好的条件，大大推进了数学信号处理技术。 1984年，法国的杜哈梅（P.Dohamel）和霍尔曼（H.Hollamann）提出的分裂基块快速算法，使运算效率进一步提高。 Advanced FFT is still developing… Principal Discoveries of Efficient Methods of Computing the DFT（部分） 如前所述，N点DFT的复乘次数等于N2。显然，把N点DFT分解为几个较短的DFT，可使乘法次数大大减少。 旋转因子WmN具有明显的周期性和对称性。 其周期性表现为 FFT算法基本上分为两大类：时域抽取法FFT(Decimation In Time FFT)和频域抽取法FFT(Decimation In Frequency FFT)。 1、DIT-FFT 设序列x(n)的长为N且满足 则x(n)的DFT为 X1(k)和X2(k)分别为x1(r)和x2(r)的N/2点DFT，即: 与第一次分解相同，将x1(r)按奇偶分解成两个N/4长的子序列x3(l)和x4(l)，即 用同样的方法可计算出 每一级运算都需要N/2次复数乘和N次复数加(每个蝶形需要两次复数加法)。所以，M级运算总共需要的复数乘次数为 1.原位计算 由上面分析可知，DIT-FFT运算很有规律。 点的FFT共进行M级运算，每级由N/2个蝶形运算组成。 2.旋转因子的变化规律 如上所述，N点DIT-FFT运算流图中，每级都有N/2个蝶形。每个蝶形都要乘以因子WpN，称其为旋转因子，p称为旋转因子的指数。 由分析可知，第L级有2 L-1个旋转因子。 N=23=8时的各级旋转因子表示如下： L=1时，WpN=WJ N/4=WJ2L， J=0 L=2时， WpN =WJ N/2=WJ2L， J=0，1 L=3时， WpN =WJN=WJ2L， J=0，1，2，3 对N=2M的一般情况，第L级的旋转因子为 3. 蝶形运算规律 设序列x(n)经时域抽选(倒序)后，存入数组X中。如果蝶形运算的两个输入数据相距B个点，应用原位计算，