数字信号处理—第二章.pptVIP

数字信号处理—第二章 第一节 直接计算DFT的问题及改进途径1、问题的提出 设有限长序列x(n)，非零值长度为N，若对x(n)进行一次DFT运算，共需多大的运算工作量？计算成本?计算速度? 2. DFT的运算量 回忆DFT和IDFT的变换式： 1）x(n)为复数， 也为复数。2）DFT与IDFT的计算量相当。注意： 计算机运算时（编程实现）： N次复乘，N-1次复加 N个点 以DFT为例： 复数乘法复数加法一个X(k)NN – 1N个X(k)(N点DFT)N 2N (N – 1)实数乘法实数加法一次复乘42一次复加2一个X (k)4N2N+2 (N – 1)=2 (2N – 1)N个X (k)(N点DFT)4N 22N (2N – 1)运算量(a+jb)(c+jd)=(ac-bd)+j(bc+ad) 例：计算一个 N点DFT ，共需N2次复乘。以做一次 复乘1μs计，若N =4096，所需时间为例：石油勘探，有24个通道的记录，每通道波形记 录长度为5秒，若每秒抽样500点/秒， 1）每道总抽样点数：500*5=2500点 2）24道总抽样点数：24*2500=6万点 3）DFT复乘运算时间：N2=(60000)2=36*108次 由于计算量大，且要求相当大的内存，难以实现实时处理，限制了DFT的应用。长期以来，人们一直在寻求一种能提高DFT运算速度的方法。 FFT便是 Cooley Tukey 在1965 年提出的的快速算法，它可以使运算速度提高几百倍，从而使数字信号处理学科成为一个新兴的应用学科。 第二节 改善DFT运算效率的基本途径 1、利用DFT运算的系数 的固有对称性和周期 性，改善DFT的运算效率。 1）对称性 2）周期性 3）可约性 2、将长序列DFT利用对称性和周期性分解为短 序列DFT的思路 因为DFT的运算量与N2成正比的，如果一个大点数N的DFT能分解为若干小点数DFT的组合，则显然可以达到减少运算工作量的效果。 N点DFTN/2点DFTN/2点DFTN/4点DFTN/4点DFTN/4点DFTN/4点DFT…….复乘： FFT算法的基本思想： 利用DFT系数的特性，合并DFT运算中的某些项 把长序列DFT→短序列DFT，从而减少运算量。FFT算法分类:时间抽选法 DIT: Decimation-In-Time频率抽选法 DIF: Decimation-In-Frequency 第三节 按时间抽选的基2-FFT算法1、算法原理 设输入序列长度为N=2M(M为正整数，将该序列按时间顺序的奇偶分解为越来越短的子序列，称为基2按时间抽取的FFT算法。也称为Coolkey-Tukey算法。 其中基2表示：N=2M，M为整数.若不满足这个条件，可以人为地加上若干零值（加零补长）使其达到 N=2M。 先将x(n)按n的奇偶分为两组，作变量置换: 当n=偶数时，令n=2r; 当n=奇数时，令n=2r+1; 分组，变量置换2、算法步骤得到： 带入DFT中 所以 由于? X1(k)、X2(k)只有N/2个点，以N/2为周期；而X (k)却有N个点，以N为周期。要用X1(k)、X2(k)表达全部的X (k) 值，还必须利用WN系数的周期特性。 后半部分前半部分又考虑到 的对称性：有： 后半部分前半部分蝶形运算流图符号说明： (1) 左边两路为输入 (2) 右边两路为输出 (3) 中间以一个小圆表示加、 减运算（右上路为相加 输出、右下路为相减输 出)1个蝶形运算需要1次复乘，2次复加 复数乘法复数加法一个N 点DFTN 2N (N–1)一个N / 2点DFT(N / 2)2N / 2 (N / 2 –1)两个N / 2点DFTN 2 / 2N (N / 2 –1)一个蝶形12N / 2个蝶形N / 2N总计N2/2 + N/2≈ N2/2N(N/2-1) + N≈ N2/2运算量减少了近一半 分解后的运算量： 先将N=8点的DFT分解成2个4点DFT：可知：时域上：x(0),x(2),x(4),x(6)为偶子序列 x(1),x(3),x(5),x(7)为奇子序列 频域上：X(0)~X(3)，由X(k)给出 X(4)~X(7)，由X(