快速傅里叶变换FFTppt课件.ppt

45 第七讲 快速傅里叶变换 (FFT) 由于 故 因此 X(k)可表示为 2.N点DFT按k的奇偶分组可分为两个N/2的DFT 当k为偶数,即 k=2r时,(-1)k =1； 当k为奇数,即 k=2r+1 时 (-1)k =-1 。 这时 X(k) 可分为两部分： k为偶数时： 可见,上面两式均为N/2的DFT。 k为奇数时： -1 3.蝶形运算 -1 -1 -1 -1 W W W W N N N N 0 1 2 3 4.N=8时,按k的奇偶分解过程 先蝶形运算，后DFT: 45-* * 45-* 45-* QA 办公室:手 机Email: tangzhyguo@163.com 本讲在分析直接计算DFT的特点的基础上介绍DFT的快速算法-快速傅里叶变换(FFT)；同时简要介绍了FFT算法的发展历程；此外还要介绍FFT的两种最常用的算法－－基于时间抽取的FFT（DIT：库利－图基算法）和基于频率抽取的FFT（DIF：桑德－图基算法）。 一、直接计算DFT存在的问题 N点序列x(n)的DFT变换定义为： 计算X(k)的运算量：需要N2次复数乘法，N(N－1)次复数加法。在N较大时计算量很大。 例如：N＝1024时, 需要1,048,576次复数乘法, 即4,194,304次实数乘法 对于象雷达、通信、声纳等需要实时处理的信号，因为其运算量更大，所以无法满足信号处理的实时性要求。迫切需要有新的算法。 二、DFT运算的特点 实际上，DFT运算中包含有大量的重复运算。在WN 矩阵中，虽然其中有N2个元素，但由于WN的周期性，其中只有N个独立的值，即 ，且这N个值也有一些对称关系。总之，WN因子具有如下所述周期性及对称性： 1.对称性 2.周期性 由上述特性还可得出： 利用上述对称特性，可使DFT运算中有些项可以合并，这样，可使乘法次数减少大约一半；利用WN矩阵的对称性及周期性，可以将长序列的DFT分解为短序列的DFT，N越小，运算量能够减少。 例如，对于四点的DFT，直接计算需要16次复数乘法，根据上述特性可以有以下形式的算法： 第二列和第三列交换，则： 则有： 由此得出： 从上例可知，通过应用对称性和周期性，4点的DFT实际上只需要进行一次复数乘法。 三、FFT发展简介 FFT的实质：快速傅里叶变换(FFT)并不是一种新的变换，是为了改进和提高离散傅里叶变换(DFT)运算速度基于DFT运算特点而发展起来的DFT快速算法，其实质还是DFT。 FFT发展的原因：DFT是信号分析与处理中的一种重要变换，广泛应用于通信、图像处理、雷达及声纳等领域，由于其计算量与变换区间长度N的平方成正比，在N较大时，计算量很大，使得直接应用DFT进行实时处理信号是不现实的。 FFT的发展历程： 1965年，J. W. Cooley和J. W. Tukey巧妙应用DFT中W因子的周期性及对称性提出了最早的FFT，这是基于时间抽取的FFT。具有里程碑式的贡献(运算量缩短两个数量级) 1966年，G. Sand提出了基于频率抽取的FFT算法 1975年，Winogard提出WFTA法；1977年Kolha和Parks提出素因子算法（PFA） 1984年，P. Dohamel和H. Hollmann提出分裂基快速算法，进一步减少了计算量，提高了计算速度（目前最理想的算法） FFT的各种算法 纵观FFT的发展历程，FFT算法分成两大类： (1) 针对N等于2的整数次幂的算法，如基2算法、基4算法、实因子算法和分裂基算法 (2)针对N不等于2的整数次幂的算法，其以Winograd为代表的一类算法(素因子法PFA、Winograd算法WFTA) 简要介绍库利-图基算法和桑得-图基算法 设序列x(n)的长度为N，且满足N＝2M，M为自然数，按n的奇偶将x(n)分解为两个N/2的子序列： x1(r)=x(2r), r=0,1,2,…,N/2-1 x2(r)=x(2r+1) r=0,1,2,…,N/2-1 则x(n)的DFT为： 1.基本原理 四、按时间抽取(DIT)的FFT－－库利-图基算法 因为 ，所以： 上式中X1(k)和X2(