快速傅立叶变换(FFT)创新.pptVIP

第6章 快速傅立叶变换(FFT) 6.1 引言 一.DFT的计算工作量 两者的差别仅在指数的符号和因子1/N. 通常x(n)和 都是复数,所以计算一个X(k)的值需要N次复数乘法运算,和 次复数加法运算.那么,计算全部N点的X(k)就要N2次复数乘法运算,N(N-1)次复数加法运算.一般来说,乘法运算要比相加运算复杂,为讨论简单起见,我们以复数乘法运算次数近似作为运算工作量的衡量标准.当N很大时,运算量将是惊人的,如N=1024,则要完成1048576 次(一百多万次)运算.这样,难以做到实时处理. 二.改进的途径 1. 的对称性和周期性 FFT算法正是基于这样的基本思想发展起来的。1965年,库利(cooley)和图基(Tukey)首先提出FFT算法.对于N点DFT,仅需(N/2)log2N 次复数乘法运算.例如N=1024=210 时，需要(1024/2)log2 210 =512*10=5120次。5120/1048576=4.88% ,速度提高20倍. 应用第三代数字信号处理器TMS320C30-40，具有50ns的单周期执行时间，每秒20000次浮点运算，完成乘法、加法及运算控制、数据存取共需约1s左右，而采用FFT算法，计算时间可缩减为2.366ms。 FFT有多种形式，但基本上可分为两类：时间抽取法和频率抽取法。 按“基数”分：基-2FFT算法；基-4FFT算法；混合基 FFT算法；分裂基FFT算法 其它方法：线性调频Z变换(CZT法) 1. 设输入序列长度为N=2L(L为正整数，将该序列按时间顺序的奇偶分解为越来越短的子序列，称为基2按时间抽取的FFT算法。也称为Coolkey-Tukey算法。 2. 其中基数2----N=2L，L为整数.若不满足这个条件，可以人为地加上若干零值（加零补长）使其达到 N=2L (1).N/2点DFT 1.先将 按n的奇偶分为两组作DFT,这样有: n为偶数时: n为奇数时: 由于: 所以,上式可表示为: 其中, 2.两点结论: (1) X (k),X (k)均为N/2点的DFT。 (2) X(k)=X (k)+W X (k)只能确定出 X(k)的k= 个； 即前一半的结果。 一个2点的DFT蝶形流图 另一个2点的DFT蝶形流图 4.蝶形运算所需系数 例子：求 N=23=8点DIF （1）先按N=8--N/2=4，做4点的DIF： 将N=8点分解成2个4点的DIF的信号流图 (2)N/2(4点)--N/4(2点)FFT(a)先将4点分解成2点的DIF： （b)一个2点的DIF蝶形流图 (c)另一个2点的DIF蝶形流图 (3)将N/4（2点）DFT再分解成2个1点的DFT(a)求2个一点的DFT (b)2个1点的DFT蝶形流图 (4)一个完整N=8的按频率抽取FFT的运算流图 利用FFT计算IFFT时在命名上应注意： (1)把FFT的时间抽取法，用于IDFT运算时，由于输入变量由时间序列x(n)改成频率序列X(k),原来按x(n)的奇、偶次序分组的时间抽取法FFT，现在就变成了按X(k)的奇偶次序抽取了。 (2)同样，频率抽取的FFT运算用于IDFT运算时，也应改变为时间抽取的IFFT。 二、改变FFT流图系数的方法 1、在IFFT的运算中，常常把1/N分解为(1/2)m，并且在M级运算中每一级运算都分别乘以1/2因子，就可得到IFFT的两种基本蝶形运算结构。 四、直接利用FFT流图方法的注意点 （1）FFT与IFFT连接应用时，注意输入输出序列的 排列顺序，是自然顺序还是倒序。 （2）FFT和IFFT共用同一个程序时，也应注意利用 FFT算法输入输出的排列顺序，即应注意自然 顺序还是例位序 （3）当把频率抽取FFT流图用于IDFT时，应改称时 间抽取IFFT流图。 （4）当把时间抽取FFT流图用于IDFT时，应改称频 率抽取IFFT流图。 6.6 任意基数的FFT算法 以P=3，Q=4, N=12为例 (1) 先将 x(n)通过 x(n1Q+n0)改写成 x(n1,n0)。因为 Q=4, n1=0,1,2, n0=0,1,2,3,故输入是按自然顺序的，即 x(0,0)=x(0) x(0,1)=x(1) x(0,2)=x(2) x(0,3)=x(3) x(1,0)=x(4