[工学]第四章数字信号处理FFT.ppt

§4-1引言 一.DFT的计算工作量 通常x(n)和 都是复数,所以计算一个 X(k)的值需要N次复数乘法运算,和 次 复数加法运算.那么,所有的X(k)就要N2次复 数乘法运算,N(N-1)次复数加法运算. 乘法比加法复杂，需要的运算时间多，尤其是复数相乘，每个复数相乘包括4个实数相乘和2个实数相加，例 又每个复数相加包括2个实数相加，所以，每计算一个 X（k）要进行4N次实数相乘和2N+2（N-1）=2（2N-1）次实数相加，因此，整个DFT运算需要4N2实数相乘和 2N（2N-1）次实数相加。 二.改进的途径 周期性也可以这样表达： 　　　　　　具有以下性质： 1.周期性 2.共轭对称性 3.可约性 4.正交性 对于N点DFT,仅需(N/2)log2N 次复数乘法运算.例如N=1024 时，只需要5120次。5120/1048576=4.88% ,速度提高20倍 §4-2 按时间抽取(DIT)的FFT算法 —库利-图基算法 一.算法原理(基2FFT) 只考虑N=2L 的情况，L：正整数 将 按n的奇偶分为两组， 　 　n为偶数时: n为奇数时: 由于: 两点结论: (1) X (k),X (k)均为N/2点的DFT。 (2) X(k)=X (k)+W X (k)只能确定出 X(k)的k= 个； 即前一半的结果。 将 4点DFT分解成2点DFT流程图： 对于N=4的情况，只进行了一次奇偶分解，全部运算过程分为两级蝶形运算。 对于 任意情况，每分一次，称为一级运算。 §4-5 线性卷积的FFT算法 利用循环卷积计算线性卷积如下图： (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的运算流图如下 ５.WNr 的确定 每一级WNr因子分布的规律是： (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算法中 因子分布的规律 ６.存储单元 输入序列x(n),n=0,1… N-1, 计N个单元; 　系数 , r=0,1…（N/2）-1, 　需N/2个存储单元； 　　共计(N+N/2)个存储单元。 §4-3 频率抽取(DIF)的基-2FFT算法 (桑德—图基算法) 一.算法原理 1.N点DFT的另一种表达形式 由于： 因此： 2.N点DFT按k的奇偶分组可分为两个N/2的DFT 当k为偶数,即 k=2r时,(-1)k =1； 当k为奇数,即 k=2r+1 时 (-1)k =-1 。 所以 X(k) 可分为两部分： X(k) 可分为两部分： k为偶数时： 注意：上面两式均为