王兰美《数字信号处理》第14节.pptVIP

第四 章 离散傅里叶变换快速算法（FFT) 1、DFT运算的特点： FFT算法的基本思想： 2、按时间抽取的FFT（N点DFT运算的分解）:DIT—FFT 将DFT运算也相应分为两组： 蝶形信号流图 通过N点FFT运算可得到： Y(k)=X1(k)+jX2(k) ，N点 根据前面的讨论，得到 时间抽取法FFT的运算特点： 第一次分偶、奇，根据最低位n0的0、1状态来分，若n0=0，则为偶序列；n0=1则为奇序列，得到两组序列： 000 010 100 110 ? 001 011 101 111 第二次对这两个偶、奇序列再分一次偶、奇序列，这就要根据n1的０、１状态。若n1=0，则为偶序列；n1=1则为奇序列，得到四组序列： 000 100 ? 010 110 ? 001 101 ? 011 111 同理，再根据n2的０、１状态来分偶、奇序列，直到不能再分偶、奇时为止。对于N=8, n2已是最高位，最后一次分得结果为 000 ? 100 ? 010 ? 110 ? 001 ? 101 ? 011 ? 111 2.4 FFT应用中的几个问题 IDFT与DFT的差别： 1）把DFT中的每一个系数 改为 ， 2）再乘以常数 1/N ， 则以上所讨论的时间抽取或频率抽取的FFT运算均可直接用于IDFT运算，当然，蝶形中的系数 应改为 。 对 的一般情况，第L级的旋转因子为 从而可以确定第L级运算的旋转因子。 1)IDFT的运算方法 以上所讨论的FFT算法可用于IDFT运算——简称为IFFT 比较IDFT的定义式： IDFT： DFT：X(k)=DFT[x(n)] = ? b）第二种方法，完全不需要改动FFT程序，而是直接利用它作 IFFT。 考虑到 故 IFFT计算分三步： ① 将X（k）取共轭（虚部乘以-1） ② 对 直接作FFT ③ 对FFT的结果取共轭并乘以1/N，得x（n）。 * DFT在时域和频域都是离散的,在数字技术中占有很重要的地位。由于DFT运算复杂，并没有得到真正的运用，而频谱分析仍大多采用模拟信号滤波的方法解决，直到1965年Tukey, Cookey提出DFT运算的一种快速算法以后，情况才发生了根本变化，人们开始认识到DFT运算的一些内在规律，从而很快地发展和完善了一套高速有效的运算方法——快速付里变换（FFT）算法。FFT的出现，使DFT的运算大大简化，运算时间缩短一～二个数量级，使DFT的运算在实际中得到广泛应用。 首先分析有限长序列 x（n）进行一次DFT运算所需的运算量。 一般，x（n）和wnkN都是复数 一个X（k）值， N次复数相乘， N-1次复数相加 N个X（k）点， N2次复数相乘， N（N-1）次复数相加， 从上面的分析看到，在DFT计算中，不论是乘法和加法，运算量均与N2成正比。因此，N较大时，运算量十分可观。例，计算N=10点的DFT，需要100次复数相乘，而N=1024点时，需要1048576（一百多万）次复数乘法，如果要求实时处理，则要求有很高的计算速度才能完成上述计算量。 利用以下两个特性可减少运算量： 1）系数 是一个周期函数，它的周期性和对称性可用来改进运算，提高计算效率。 例 又如 因此 利用这些周期性和对称性，使DFT运算中有些项可合并； 2）利用 的周期性和对称性，把长度为N点的大点数的DFT运算依次分解为若干个小点数的DFT。因为DFT的计算量正比于N2，N小，计算量也就小。 FFT算法正是基于这样的基本思想发展起来的。它有多种形式，但基本上可分为两类：时间抽取法和频率抽取法。 基2FFT N=2M，M：正整数 首先将序列x（n）分解为两组，一组为偶数项，一组为奇数项，