第四章快速付里叶变换2.pptVIP

第四章 快速付里叶变换 FFT: Fast Fourier Transform 复习 什么是FFT？ 直接计算DFT的工作量有多大？ 的性质和特殊点？ 周期性、对称性、可约性。 时域抽取法基2FFT基本原理？ 时间抽取法基2FFT具体的运算步骤？ 基2FFT具体的运算步骤 时间抽取基2FFT流图绘制 N个输入，N个输出； 输入为倒序码，输出为顺序码； N＝2M，表示运算共M级数，取值为0～M－1； 蝶形运算两节点的距离: 2m，其中m表示第m级 每一级都有N/2个蝶形单元，可以分成若干组； 第m级的组数为： 每一组具有相同的结构，相同的 因子分布； 4.2.5 频域抽取法FFT(DIF-FFT) 2.代入DFT中 4.结论 三、蝶形流图表示 蝶形流图的另一种形式 例子：求 N=23=8点DIF （1）先按N=8--N/2=4，做4点的DIF： 将N=8点分解成2个4点的DIF的信号流图 N=8点分解成2个4点的另一种形式 (2)N/2(4点)--N/4(2点)FFT(a)先将4点分解成2点的DIF： （b）一个2点的DIF蝶形流图 (c)另一个2点的DIF蝶形流图 (3)将N/4（2点）DFT再分解成2个1点的DFT (b)2个1点的DFT蝶形流图 (4)一个完整N=8的按频率抽取FFT的运算流图 (5)DIF与DIT比较 1.相同点 (1)进行原位运算； (2)运算量相同,均为（N/2) Log2N次复乘,N Log2N次复加。 2.不同点 (1)DIT输入为倒位序,输出为自然顺序； DIF正好与此相反。 (2)DIF与DIT根本区别：在于蝶形结不同。 DIT的复数相乘出现在减法之前。 DIF的复数相乘出现在减法之后。 4.2.6 IDFT的高效算法 一、利用FFT计算IFFT的思路1 二、利用FFT计算IFFT的思路2 作业 课后习题1。 试画出N=16点的基-2按时间抽取的FFT流图（DIT）。 试画出N=16点的基-2按频率抽取的FFT流图（DIF）。 将下列两式进行比较 利用FFT计算IFFT时在命名上应注意： (1)把FFT的时间抽取法，用于IDFT运算时，由于输入变量由时间序列x(n)改成频率序列X(k),原来按x(n)的奇、偶次序分组的时间抽取法FFT，现在就变成了按X(k)的奇偶次序抽取了。 (2)同样，频率抽取的FFT运算用于IDFT运算时，也应改变为时间抽取的IFFT。 三、改变FFT流图系数的方法1.思路 在IFFT的运算中，常常把1/N分解为(1/2)m，并且在M级运算中每一级运算都分别乘以1/2因子，就可得到IFFT的两种基本蝶形运算结构。 2.IFFT的基本蝶形运算 A B A B (a)频率抽取IFFT的蝶形运算 (b)时间抽取IFFT的蝶形运算 四.不改(FFT)的程序直接实现IFFT * * X1(k) X2(k) 实现上式运算的流图称作蝶形运算 X1(k):偶中偶 X2(k):偶中奇 进行N/4点的DFT,得到 X3(k):偶中偶 X4(k):偶中奇 X5(k):奇中偶 X6(k):奇中奇 2点DFT 2点DFT x(0) x(4) x(2) x(6) X3(0) X3(1) X4(0) X4(1) X1(0) X1(1) X1(2) X1(3) 2点DFT 2点DFT x(1) x(5) x(3) x(7) X5(0) X5(1) X6(0) X6(1) X2(0) X2(1) X2(2) X2(3) (0) (4) (2) (6) (1) (5) (3) (7) W N 0 W N 0 W N 0 W 0 N -1 -1 -1 -1 X (0) X (1) X (0) X (1) X (0) X (1) X (0) X (1) 3 3 4 4 5 5 6 6 W N 0 W N 2 W N 0 W N 2 -1 -1 -1 -1 X (0) X (1) X (2) X (3) X (0) X (1) X (2) X (3) 1 1 1 2 1 2 2 2 W W W W N 0 N 1 N 2 N 3 -1 -1 -1 -1 X(0) X(1) X(2) X(3) X(4) X(5) X(6) X(7) 因此,8点DFT的FFT的运算流图 0 0 0 0 0 0 0 0 1 0 0 1 1 0 0 4 2 0 1 0 0 1 0 2