第4章 快速付里叶变换FFT.ppt

第四章快速付里叶变换(FFT)Fast FourierTransforming;§第1节 直接计算DFT算法存在的问题及改进途径;2.以DFT为例，计算DFT复数运算量;4.计算DFT需要的实数运算量;例子;二、改善DFT运算效率的基本途径;利用DFT运算的系数 的固有对称性和周期性，改善DFT的运算效率;例子;1 合并法：步骤1分解成虚实部;合并法：步骤3合并有些项;2、将长序列DFT利用对称性和周期性分解为短序列DFT--方法;结论; 一、算法原理 设输入序列长度为N=2M(M为正整数，将该序列按时间顺序的奇偶分解为越来越短的子序列，称为基2按时间抽取的FFT算法。也称为Coolkey-Tukey算法。 其中基数2----N=2M，M为整数.若不满足这个条件，可以人为地加上若干零值（加零补长）使其达到 N=2M;例子;二、算法步骤1.分组，变量置换;2.代入DFT中;3.求出子序列的DFT;4.结论1;5.求出后半部的表示式;6.结论2;三、蝶形结;例子：求 N=23=8点FFT变换 （1）先按N=8--N/2=4，做4点的DFT：;（b)求 一个蝶形结需要的运算量;（c）分解为两个N/2=4点的DFT的运算量;(d)用2个4点来求N=8点的FFT所需的运算量;(e)将N=8点分解成2个4点的DFT的信号流图 ;(2)N/2(4点)--N/4(2点)FFT(a)先将4点分解成2点的DFT：;(b)求2点的DFT;（c)一个2点的DFT蝶形流图;（d)另一个2点的DFT蝶形流图;(3)将N/4（2点）DFT再分解成2个1点的DFT(a)求2个一点的DFT;(b)2个1点的DFT蝶形流图 ;按时间抽取将一个N点DFT分解 为两个N/2点DFT（N=8） ;图 一个N=8点DFT分解为四个N/4点DFT ;;四、FFT算法中一些概念 （1）“级”概念; （2）“组”概念 每一级都有N/2个蝶形单元，例如：N=8,则每级都有4个蝶形单元。每一级的N/2个蝶形单元可以分成若干组，每一组具有相同的结构，相同的 因子分布，第m级的组数为：;(3) 因子的分布;（4）按时间抽取法;五、按时间抽取的FFT算法与直接计算DFT运算量的比较;例子;六、按时间抽取FFT算法的特点;1.原位运算（in-place);例子;2.码位倒读规则;例子;图倒位序的形成 ;整序重排子程序;七、按时间抽取的FFT算法的若干变体1.思路;(2)输入是自然顺序而输出是乱序;(3)输入和输出都是自然顺序;一、算法原理;二、算法步骤1.分组;2.代入DFT中;3. 变量置换--1;3. 变量置换--2;3. 变量置换--3;3. 变量置换--4;4.结论1;4.结论2;三、蝶形流图表示;例子：求 N=23=8点DIF （1）先按N=8--N/2=4，做4点的DIF：;将N=8点分解成2个4点的DIF的信号流图 ;(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(a)求2个一点的DFT;(b)2个1点的DFT蝶形流图 ;(4)一个完整N=8的按频率抽取FFT的运算流图 ;（6）DIF与DIT比较1;（6）DIF与DIT比较2;(4)一个完整N=8的按频率抽取FFT的运算流图 ;例子;第五节IFFT运算方法;一、利用FFT计算IFFT的思路1;利用FFT计算IFFT时在命名上应注意： (1)把FFT的时间抽取法，用于IDFT运算时，由于输入变量由时间序列x(n)改成频率序列X(k),原来按x(n)的奇、偶次序分组的时间抽取法FFT，现在就变成了按X(k)的奇偶次序抽取了。 （2）同样，频率抽取的FFT运算用于IDFT运算时，也应改变为时间抽取的IFFT。;二、改变FFT流图系数的方法1.思路;2.IFFT的基本蝶形运算;三.直接利用FFT流图的方法1.思路;2.直接利用FFT流图方法的推导;3.直接利用FFT流图方法的注意点;作业;第五节FFT的应用;一、快速卷积;1.利用圆周卷积代替线卷积;2、用FFT计算y(n)的步骤;2、用FFT计算y(n)与直接计算y(n)的运算量比较;3、分段卷积的方法;（1）重叠相加法1;（1）重叠相加法2;（1）重叠相加法例子;（2）重叠保留法1;（2）重叠保留法2;（2）重叠保留法例子;二、频谱分析;1.频谱分析仪的框图;2.部件说明;直接用DFT计算的运算量与用FFT计算的运算量比较，减少运算量的途径