第四章 离散傅里叶变换.pptVIP

那么，X1(k)又可表示为第28页，共71页，星期日，2025年，2月5日X2(k)也可以进行相同的分解：注意：通常我们会把写成。第29页，共71页，星期日，2025年，2月5日N点DFT的第二次时域抽取分解图(N=8)2点DFT2点DFT2点DFT2点DFTx(0)x(4)x(2)x(6)x(1)x(5)x(3)x(7)X3(0)X3(1)X4(0)X4(1)X5(0)X5(1)X6(0)X6(1)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)4点DFT4点DFTx(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)第30页，共71页，星期日，2025年，2月5日88X3(0)X3(1)x(0)=x3(0)x(4)=x3(1)第31页，共71页，星期日，2025年，2月5日N点DIT―FFT运算流图(N=8)x(0)x(4)x(2)x(6)x(1)x(5)x(3)x(7)X(0)X(1)X(2)X(3)X(4)X(5)X(6)X(7)第32页，共71页，星期日，2025年，2月5日3、DIT―FFT算法与直接计算DFT运算量的比较1)、N=2M的DFT运算可分成M级，每一级有N/2个蝶形，每个蝶形有一次复乘两次复加。2)、所以M级共有次复乘和次复加。3)、若直接计算DFT，需N2次复乘和N(N-1)次复加。显然，当N较大时，有：例如，N=210=1024时第33页，共71页，星期日，2025年，2月5日FFT算法与直接计算DFT所需乘法次数的比较曲线第34页，共71页，星期日，2025年，2月5日4、DIT―FFT的运算规律及编程思想FFT的每级（列）计算都是由N个复数数据（输入）两两构成一个蝶型（共N/2个蝶形）运算而得到另外N个复数数据（输出）。当数据输入到存储器以后，每一组运算的结果，仍然存放在这同一组存储器中直到最后输出。例：将x(0)放在单元A(0)中，将x(4)放在单元A(1)中，W80放在一个暂存器中。将x(0)+W80x(4)→送回A(0)单元将x(0)-W80x(4)→送回A(1)单元X3(0)X3(1)x(0)x(4)1）原位运算(亦称同址计算)第35页，共71页，星期日，2025年，2月5日x(0)x(4)x(2)x(6)x(1)x(5)x(3)x(7)X(0)X(1)X(2)X(3)X(4)X(5)X(6)X(7)回顾：N点DIT―FFT运算流图(N=8)第36页，共71页，星期日，2025年，2月5日如上所述，N点DIT―FFT运算流图中，每级都有N/2个蝶形。每个蝶形都要乘以因子WNP，称其为旋转因子，p称为旋转因子的指数。2)旋转因子的变化规律观察FFT运算流图发现，第L级共有2L-1个不同的旋转因子。N=23=8时的各级旋转因子表示如下：L=1时，WNp=WN/4J，N/4=21=2L，J=0L=2时，WNp=WN/2J，N/2=22=2L，J=0，1L=3时，WNp=WNJ，N=23=2L，J=0，1，2，3第37页，共71页，星期日，2025年，2月5日对N=2M的一般情况，第L级的旋转因子为：第38页，共71页，星期日，2025年，2月5日设序列x(n)经时域抽选(倒序)后，存入数组X中。如果蝶形运算的两个输入数据相距B个点(B=2L-1)，应用原位计算，则蝶形运算可表示成如下形式：下标L表示第L级运算，XL(J)则表示第L级运算后数组元素X(J)的值。第39页，共71页，星期日，2025年，2月5日3)编程思想及流程图开始送入x(n)和N=2M调整输入x(n)的顺序for(L=1;L=M;L++)B=2L-1for(J=0;J=B-1;J++)p=J·2M-Lfor(k=J;