数字信号处理第2版张小虹4数字信号处理2课件教学.pptVIP

以此类推，L≥3时，第L级的两个无关紧要的旋转因子 WN与WN 减少复数乘法的次数为2N / 2L = N/2L-1 。 N/4 0 从 L=3 到 L=M 共减少复数乘法的次数为 N 2L?1 L=3 ? M (1/2)L L=3 ? M =2N [(1/2)3 +(1/2)4 +???+ (1/2)M] L=3 ? M =2N =2N (1/2)3 [1+(1/2)4 +???+ (1/2)M?3] 1?(1/2)M?2 1?1/2 =2N (1/2)3 =2N (1/2)3?2?[1?(1/2)M?2] =(N/2) [1?(1/2)M?2] =(N/2) ?2 考虑以上这些减少的复数乘法后，这时的复数乘法运算 次数为 实现一次复数乘法，一般需要四次实数乘法，两次实数 2、减少运算的因子 相乘时，有 N N CM(2)= (M?2)?( ?2) 2 2 N 2 = (M?3)+2 数加法。但当任意复数 a + jb 与WN = (1?j) ? 2 /2 N/8 实数加法和两次实数乘法。 由式(4.4-6) 可见，任意复数与 WN 相乘只需要两次 N/8 (1?j) ? 2 (a + jb) 2 ? 2 (a+jb?ja+b) 2 = ? 2 [ (a+b)?j(a?b) ]=A+jB 2 = (b?a) B= ? 2 2 (4.4-6) ? 2 (a+b) 2 A = 式中 相同为 (N/2)?2次。 这样，所有含有WN 的蝶形，减少两次实数乘法。 N/8 从 L=3 到 L = M级，每级都包含WN ，第 L 级中， N/8 WN 对应2M-L = N/2L个蝶形运算，所以从第三级到 N/8 最后一级， WN 减少的实数乘法次数与 (4.4-4)式 N/8 3、改进后的实数乘法计算量 考虑所有相关的旋转因子，从实数运算考虑，计算 N = 2M点FFT的实数乘法次数为 13 2 =N (2M? ) +10 N RM(2)= 4 2 N 2 (M?3)+2 ? ?2 在基2FFT程序中，含有全部旋转因子的，该算法为 一类蝶形单元运算；去掉了 WN = ?1旋转因子的，为 r 二类蝶形单元运算；去掉了 WN = ?1 , ?j旋转因子的， r 为三类蝶形单元运算；去掉了 WN = ?1 , ?j旋转因 r 子，又对 WN 作了处理的，为四类蝶形单元运算。 N/8 类型越多，编程越复杂。当N较大时，乘法运算的减少 乘法次数是一类蝶形单元运算的75%。 量是相当可观的。例N=4096时，三类蝶形单元运算的 后三种运算也称为多类蝶形单元运算。显然，蝶形单元 的生成方法影响FFT运算速度。 实际应用中，旋转因子WN =cos(2?m/N) ?jcos(2?m/N) m WN ， 0≤m≤(N/2)?1，存放在数组中。在程序执行 m 产生旋转因子方法有两种，一种方法是每级运算时 直接计算产生，优点是实时运算占用内存少，缺点 是运算速度受限。另一种方法是事先计算出所有的 过程中查表得到所需系数，优点是运算速度快，缺 点是要占较多内存。 * 在线教务辅导网： 更多课程配套课件资源请访问在线教务辅导网 4.3 基2频率抽取 FFT算法 1、将N点DFT分解为两个N/2点DFT 序列基础上。 4.3 .1 基2频率抽取 FFT算法 础上；频选法FFT建立在把X(k)分解成越来越短的子 做法：将x(n)分成前后两段 时选FFT建立在把将x(n)分解成越来越短的子序列基 N = 2M 对第二项令 n?=n?N/2 n=n?+N/2 n? :0~(N/2)?1 n : N/2~N?1 X(k)= kn n=0 ? N?1 x(n)WN = kn n=0 ? N/2?1 x(n)WN + kn n=N/2 ? N?1 x(n)WN kn n=0 ? N/2?1 x(n)WN + n?=0 ? N/2?1 x(n?+N/2)WN X(k)= (n?+N/2)k n=0 ? N/2?1 [x(n) +(?1)k x(n+N/2)]WN = nk WN =(?1)k Nk/2 N/2点DFT(X(k)的偶次项) r=0,1,2,???, N/2?1 X(2r)=X1(r) X(2r+1)=X2(r) 令 X(k)= 则 n=0 ? N/2?1 [x(n) +x(n+N/2)]WN 2rn X(2r)= n=0 ? N/2?1 x1(n)WN/2 rn = =X1(r) x1(n)= x(n) +x(n+N/2) N/2点DFT(X(k)的奇次项) n=0 ? N/