三、周期性矩形脉冲频谱.doc

山东理工大学备课纸 FFT快速傅式变换 3.2 直接计算DFT的问题及改进的途径 3.2.1直接计算DFT的运算量问题 两者的差别仅在指数的符号和因子1/N. 一个X(k)的值的工作量,如X(1) 通常x(n)和都是复数,所以计算一个X(k)的值需要N次复数乘法运算,和次复数加法运算.那么,所有的X(k)就要N2次复数乘法运算,N(N-1)次复 数加法运算. N点DFT运算量：复数乘法： N2 复数加法： N(N-1) 3.2.2 改进的途径 的对称性、周期性、可约性 对称性: 周期性： 可约性： 3.3 按时间抽取(DIT)的基－2FFT算法（库利-图基算法） 3.3.1 算法原理 (一)N/2点DFT 1.先将x(n) 按n的奇偶分为两组作DFT,设N=2M , 这样有: n为偶数时: n为奇数时: (n为偶数) (n为奇数) 2.两点结论: (1) X1 (k),X2 (k)均为N/2点的DFT。 (2) X(k)=X (k)+W X (k)只能确定出(k)的k=个；即前一半的结果。 3.X(k)的后一半的确定 由于 （周期性） 这就是说,X1(k),X2(k)的后一半,分别等于其前一半的值。 又由于 ,所以 可见,X(k)的后一半，也完全由X1(k), X2 (k)的前一半所确定。 *N点的DFT可由两个N/2点的DFT来计算。 4.蝶形运算 前一半 后一半 实现上式运算的流图称作蝶形运算（N/2个蝶形) X X(k) X(k+N/2) -1 WNnk X1(k) X2(k) 5.计算工作量分析 按奇、偶分组后的计算量： (1)N/2点的DFT运算量:复乘次数: 复加次数: (2)两个N/2点的DFT运算量:复乘次数: 复加次数: (3)N/2个蝶形运算的运算量:复乘次数: 复加次数: 总共运算量: 复乘次数: 复加次数: 可见分解后,计算工作点差不多减少一半 例如 N=8 时的DFT,可以分解为两个N/2=4点的DFT.具体方法如下: (1)n为偶数时,即 分别记作: (2) n为奇数时,分别记作: (3)对X (k)和X (k)进行蝶形运算,前半部为X(0) X(3),后半部分为X(4) X(7) 整个过程如下图所示: N/2 N/2点 DFT X1(0) X1 (1) X1 (2) X1 (3) X(0) X(1) X(2) X(3) x(0) x(2) x(4) x(6) N/2点 DFT X2 (0) X2 (1) X2 (2) X2 (3) X(4) X(5) X(6) X(7) x(1) x(3) x(5) x(7) WN1 WN2 WN3 -1 -1 -1 -1 WN0 (二) N/4点DFT 由于N=2 M 进行N/4点的DFT,得到 (偶中偶) (偶中奇) 前N/4的X1(k)： 后N/4的X1(k)为： 同样对n为奇数时，N/2点分为两个N/4点的序列进行DFT，则有: (奇中偶) (奇中奇) 例如,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) 。 这样,又