第四章 快速傅里叶变换.docVIP

快速傅里叶变换 4．1引言 4．2直接计算DFT的问题及改进的途径 一．DFT的计算量 两者的差别仅在指数的符号和因子1/N。 一个X(k)的值的计算量,如X(1) 通常x(n)和 都是复数,所以计算一个X(k)的值需要N次复数乘法运算,和N-1次复数加法运算.那么,所有的X(k)就要N2次复数乘法运算,N(N-1)次复数加法运算.当N很大时,运算量将是惊人的,如N=1024,则要完成1048576 次(一百多万次)运算.这样,难以做到实时处理. 二.改进的途径 1.的对称性和周期性 对称性： 周期性： 得： 利用上述特性,可以将有些项合并,并将DFT分解为短序列,从而降低运算次数,提高运算速度.1965年,库利(cooley)和图基(Tukey)首先提出FFT算法.对于N点DFT,仅需(N/2)log2N 次复数乘法运算.例如N=1024=210 时，需要(1024/2)log2 210 =512*10=5120次。 5120/1048576=4.88% ,速度提高20倍 4．3按时间抽选(DIT)的基-2 FFT算法（库利-图基算法） 算法原理(基2FFT) (一)N/2点DFT 1．将按n的奇偶分为两组作DFT,设N=2L ,不足时,可补些零。这样有: n为偶数时: n为奇数时: 因此： 由于： 所以，上式可表示为： 其中： 2．两点结论: (1) X (k),X (k)均为N/2点的DFT。 (2) X(k)=X (k)+W X (k)只能确定出X(k)的k=个；即前一半的结果。 3．X(k)的后一半的确定 由于（周期性）, 所以: 同理： 这就是说,X1(k),X2(k)的后一半,分别等于其前一半的值。 又由于, 所以： 可见,X(k)的后一半，也完全由X1(k), X2 (k)的前一半所确定。 *N点的DFT可由两个N/2点的DFT来计算。 4.蝶形运算 由X1(k)、X 2(k)表示X(k)的运算是一种特殊的运算-碟形运算（N/2个蝶形） 前一半 后一半 5.计算工作量分析 按奇、偶分组后的计算量： (1)N/2点的DFT运算量：复乘次数： 复加次数： (2)两个N/2点的DFT运算量：复乘次数： 复加次数： (3)N/2个蝶形运算的运算量：复乘次数： 复加次数： *但是,N点DFT的复乘为N2 ;复加N(N-1);与分解后相比可知,计算工作点差不多减少 一半。 例如 N=8 时的DFT,可以分解为两个N/2=4点的DFT.具体方法如下: n为偶数时,即 分别记作: n为奇数时, 分别记作: (3)对X (k)和X (k)进行蝶形运算,前半部为X(0) X(3),后半部分为X(4) X(7) 整个过程如下图所示: (二) N/4点DFT 由于N=2 L ,所以 N/2仍为偶数,可以进一步把每个N/2点的序列再按其奇偶部分分解为两个N/4的子序列。 例如，n为偶数时的 N/2点分解为: 进行N/4点的DFT,得到 (偶中偶) (偶中奇) 从而可得到前N/4的X1(k) 后N/4的X1(k)为 例如,N=8时的DFT可分解为四个N/4的DFT,具体步骤如下: (1) 将原序列x(n)的“偶中偶”部分: 构成N/4点DFT，从而得到X3(0), X3(1)。 (2) 将原序列x(n)的“偶中奇”部分: 构成N/4点DFT，从而得到X4(0), X4(1)。 (3) 将原序列x(n)的“奇中偶”部分: 构成N/4点DFT，从而得到X5(0), X5(1)。 (4) 将原序列x(n)的“奇中奇”部分: 构成N/4点DFT，从而得到X6(0), X6(1)。 （5）由 X3(0), X3(1), X4(0), X4(1)进行碟形运算, 得到X1(0), X1(1), X1(2), X1(3)。 （6）由 X5(0), X5(1), X6(0), X6(1)进行碟形运算, 得到X2(0), X2(1), X2(2), X2(3)。 （7）由X1(0), X1(1), X1(2), X1(3)，X2(0), X2(1),X2(2),X2(3)进行碟形运算, 得到 X(0), X(1), X(2), X(3) X(4), X(5), X(6), X(7) 。 这样,又一次分解,得到四个N/4点DFT,两级蝶形运算,其运算量有大约减少一半，即为N点DFT的1/4。