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

* 在线教务辅导网： 教材其余课件及动画素材请查阅在线教务辅导网 QQ:349134187 或者直接输入下面地址： 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/2?1 [x(n)?x(n+N/2)]WN (2r+1)n X(2r+1)= n=0 ? N/2?1 [x(n)?x(n+N/2)]WN WN 2rn = n n=0 ? N/2?1 x2(n)WN/2 rn = =X2(r) x2(n) =[x(n)?x(n+N/2)]WN n 上式的运算关系可以用蝶形表示 运算量：一次复乘；两次复加。 n=0,1,2, …(N/2)?1 (4.3-3) x2(n) =[x(n)?x(n+N/2)]WN n 即有 x1(n)= x(n) +x(n+N/2) x(n) +x(n+N/2) [x(n)?x(n+N/2)]WN n x(n) x(n+N/2) WN n 合成N点X(k) 均为N/2点DFT 同样是利用两个N/2点DFT，合成一个N点DFT。 (4.3-4) n=0 ? N/2?1 x2(n)WN/2 rn X(2r+1) = =X2(r) 即 n=0 ? N/2?1 x1(n)WN/2 rn X(2r) = =X1(r) x1(n) ? X1(k) N/2 x2(n) ? X2(k) N/2 x(n) = 例 N=8 x(0) x(1) x(2) x(3) x(4) x(7) x(5) x(6) X(0) X(4) X(2) X(6) X(1) X(7) X(5) X(3) N/2点 DFT N/2点 DFT WN 0 WN 1 WN 2 WN 3 复乘计算量： 2、继续将N/2点DFT分解为两个N/4点DFT。 对第二项令 n?=n?N/4 n=n?+N/4 n? :0~(N/4) ?1 n : N/4~(N/2) ?1 2 N 2 @ mF = 2? N 2 N 2 2 + 2 N 2 = N 2 + n=0 ? N/2?1 x1(n)WN/2 rn X1(r)= n=0 ? N/4?1 x1(n)WN/2 rn n=N/4 ? N/2?1 x1(n)WN/2 rn + = rn n=0 ? N/4?1 x1(n)WN/2 + n?=0 ? N/4?1 x1(n?+N/4)WN/2 (n?+N/4)r X1(r)= l = 0,1,2, …(N/4)?1 rn n=0 ? N/4?1 x1(n)WN/2 + n?=0 ? N/4?1 x1(n?+N/4)WN/2 (n?+N/4)r X1(r)= WN/2 =(?1)r Nr/4 n=0 ? N/4?1 [x1(n) +(?1)r x1(n+N/4)]WN/2 = nr X1(2l)=X3(l) X1(2l+1)=X4 (l) 令 X1(r)= n=0 ? N/4?1 [x1(n) +x1(n+N/4)]WN/2 2ln X1(2l)=X3(l)= n=0 ? N/4?1 x3(n)WN/4 ln = N/4点DFT(X1(k)的偶次项) N/4点DFT(X1(k)的奇次项) n=0 ? N/4?1 [x1(n) ?x1(n+N/4)]WN/2 (2l+1)n X1(2l+1)=X4(l)= n=0 ? N