数字信号处理_9new.pptVIP

* * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * (3) 将原序列x[k]的“奇中偶”部分: 构成N/4点DFT，从而得到X5[0], X5[1]。 (4) 将原序列x[k]的“奇中奇”部分: 构成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]。 N/4点 DFT x5[0]= x2[0] =x[1] x5[1]= x2[2] =x[5] X5[0] X5[1] x6[0]= x2[1] =x[3] x6[1]= x2[3] =x[7] N/4点 DFT X6[0] X6[1] -1 W40 X2[0] X2[2] -1 W41 X2[1] 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。 对于N=8时DFT，N/4点即为两点DFT,因此 4点DFT 4点DFT x[0] x[2] x[4] x[6] x[1] x[3] x[5] x[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] -1 -1 -1 -1 8点基2时间抽取FFT算法流图 基2时间抽取FFT算法流图 第一级 第二级 第三级 这种FFT算法,是在时间上对输入序列的次序是属于偶数还是属于奇数来进行分解的，所以称作按时间抽取的算法(DIT)。 基2时间抽取FFT算法流图 N=2 x[k]={x[0], x[1]} 4点基2时间抽取FFT算法流图 x[0] x[2] x[1] x[3] X1[0] X1[1] X2[0] X2[1] 2点DFT 2点DFT -1 -1 -1 -1 X [0] X [1] X [2] X [3] 4点基2时间抽取FFT算法流图 8点基2时间抽取FFT算法流图 4点DFT 4点DFT x[0] x[2] x[4] x[6] x[1] x[3] x[5] x[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] -1 -1 -1 -1 4点DFT 4点DFT x[0] x[2] x[4] x[6] x[1] x[3] x[5] x[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] -1 -1 -1 -1 8点基2时间抽取FFT算法流图 第一级 第二级 第三级 8点基2时间抽取FFT算法流图 二、运算量 当N = 2L时，共有L级蝶形，每级N/2个蝶形，每个蝶形有1次复数乘法2次复数加法。 当N=8 =23时， 由3级蝶形运算组成； 每级由4个蝶形运算组成；共有12个蝶形运算； 每个蝶形运算有1次复乘，2次复加； 共有12次复乘，24次复加； x[0] x[4] WN0 -1 X3(0) X3(1) x[2] x[6] WN0 -1 X4(0) X4(1) x[1] x[5] WN0 -1 X5(0) X5(1) x[3] x[7] WN0 -1 X6(0) X6(1) -1 WN0 X1(0) X1(2) -1 WN2 X1(1) X1(3) -1 WN0 X2(0) X2(2) -1 WN2 X2(1) X2(3) -1 WN0 X[0] X[4] -1 WN1 X[1] X[5] -1 WN2 X[2] X[6] -1 WN3 X[3] X[7] 复数乘法： 复数加法： 复数乘法次数