数字信号处理-第四章快速傅立叶变换(new).ppt

由于 K取偶数 K取奇数 令 则 当N=2时， －1 当N=4时， －1 2点DFT,一个蝶形运算 －1 而 是一对蝶形运算输出 而 是一对蝶形运算输出 －1 －1 －1 －1 数字信号处理－第四章 快速傅立叶变换 4.2 直接计算DFT的问题及改进的途径 设 为N点有限长序列，其DFT为 一般： 为复数，计算一个 需要N次复数乘法及N-1次复数加法。由于 一共有N个点，因此完成整个DFT运算需要 次复数乘法及 复数加法。 （1） （2） 由于 所以 一次复数乘法需用四次实数乘法，一次复数加法则需二次实数加法。 每运算一个 需4N次实数乘法及4（N-1)+2=4N-2次实数加法。所以这个DFT运算总共需要 次实数乘法和 N（4N-2)次实数加法。 K一定：N-1次加法 K一定：N-1次加法 K一定：N-1次加法 K一定：N-1次加法 W因子性质： N，k均以N为周期 4.3 按时间抽选（DIT)的基－2FFT算法（库利－图基算法） 一：算法原理 基－2FFT：序列点数N为2的整数幂的FFT。 将 的序列 （n=0,1,…,N-1)先按n的奇偶分成以下两组： 蝶形运算 令： N点 DFT . . . N/2点 DFT . . . 蝶形 运算 N/2点 DFT . . . 蝶形运算 … N/2点 DFT . . . N/4点 DFT . . . 蝶形 运算 N/4点 DFT . . . N/2点 DFT . . . N/4点 DFT . . . 蝶形 运算 N/4点 DFT . . . 令 (4) (3) 如下 K周期为N/2 将式（3）、式（4）用蝶形信号流图符号表示。当支路未标出系数时，则该支路的传输系数为1 图4-1时间抽选法蝶形运算流图符号 －1 如果N=2 －1 N=2 如果N=4 －1 N=4 －1 －1 －1 N=8点的DFT可以分解成两个4点的DFT及蝶形运算 上半部分 下半部分 序号为偶次 序号为奇次 然后将 的DFT继续分解，分解成 的DFT及蝶形运算。 上半部分 下半部 分 序号为偶次 序号为奇次 DIT-FFT算法特点： （1）同址运算：每一级的每个蝶形的输入与输出数据，在运算前后可以存储在同一内存单元上称为同址运算。 （2）倒位序规律（基2）：DIT的基2FFT运算，其输出X(k)是自然顺序，但输入x(n)是码位倒置的顺序。在实际运算中，先按自然顺序输入x(n)存入存储单位，再通过变址运算，将自然顺序的存储转换为码位倒置顺序的存储。 （3）蝶距及W因子的变换规律：各类蝶形运算两个点相距的“距离”称蝶距。 蝶距规律：最后一级的蝶距为N/2,依次向左为N/4,N/8,… W因子规律：最后一级有一组，每组有N/2种，顺序为 向左一级有两组，每组N/4种，顺序为 其余类推。 m：蝶形运算的第m级 二、运算量 按时间抽选法FFT的流图可见，当 ,共有 L级碟形运算，每级都有 个碟形运算组成，每个碟形有一次复乘（4次 实乘）、二次复加（4次实加） ，这样FFT 运算总共需要复乘数： 运算总共需要实乘数： 复加数： 实加数： 1、如果一台通用计算机的速度为平均每次复乘 ，每次复加 比较1024点DFT[x(n)]与FFT[x(n)]复乘次数，复加次数 解: ⑴ 直接DFT计算:  复乘所需时间:  复加所需时间: ⑵用FFT计算:  复乘所需时间:  复加所需时间: 4.4 按频率抽选（DIF）的基－2 FFT算法（桑德－图基算法） 按频率抽选（DIF）的基－2 FFT算法：把输出序列 （也是N点序列）按的其顺序的奇偶分解为越来越短的序列。 IFFT IFFT 蝶形运算 N点 DFT . . . N/2点 DFT . . . 蝶形 运算 N/2点 DFT . . . . . . . . . . . . IFFT IFFT 蝶形运算 IF