课件第4章 快速傅立叶变换.pptVIP

自然顺序n 二进制码表示 码位倒读 码位倒置顺序n’ 以N=8为例： 0 1 2 3 4 5 6 7 000 001 010 011 100 101 110 111 000 100 010 110 001 101 011 111 0 4 2 6 1 5 3 7 看出：码位倒读后的顺序刚好是数据送入计算机内的顺序。 四、按频域抽取的基2-FFT算法 设序列x(n)长度为N=2M，首先将x(n)前后对半分开，得到两个子序列，其DFT可表示为如下形式 -1 DIF―FFT一次分解运算流图(N=8) -1 -1 -1 -1 W N 0 W N 1 W N 2 W N 3 两点DFT及其蝶形图 -1 考虑到 上图可变形为: -1 DIF―FFT运算流图(N=8) 0 N W 1 N W 2 N W 3 N W -1 -1 -1 -1 x(0) x(3) x(1) x(2) x(4) x(5) x(6) x(7) 0 N W 2 N W 2 N W 0 N W X(0) X(6) X(4) X(2) X(1) X(5) X(3) X(7) 0 N W 0 N W 0 N W 0 N W -1 -1 -1 -1 -1 -1 -1 -1 时间抽取算法与频率抽取算法的比较 1) 频率抽选法和时间抽选法总的计算量是相同的 复乘： 复加： 2) 频率抽取法和时间抽取法一样，都适用于原位运 算， 即蝶形的输入和输出占用同一个存储单元。 3) 均存在码位倒序问题。 4) 频率抽选法和时间抽选法一样，基本运算也是蝶形 运算。但两者的蝶形形式略有不同。 * * 第四章 快速傅立叶变换 Fast Fourier Transform FFT产生故事 虽然频谱分析和DFT运算很重要，但在很长一段时间里，并没有得到真正的运用。而频谱分析仍大多采用模拟信号滤波的方法解决。 Cooley J W, Tukey J W. An algorithm for the machine computation of complex Fourier series. Mathematics of Computation, 1965, pp297~301 提出一种快速计算DFT的方法，揭开了FFT发展史上的第一页。 1984年，法国的杜梅尔和霍尔曼将基2分解和基4分解糅合在一起，提出了分裂基FFT算法,使运算效率进上步提高。 一、直接计算DFT的问题 1、问题的提出 设有限长序列x(n)，非零值长度为N，若对x(n)进行一次DFT运算，共需多大的运算工作量？ 2. DFT的运算量 1）x(n)为复数， 也为复数。 2）DFT与IDFT的计算量相当。 注意： 计算机运算时（编程实现）： N次复乘，N-1次复加 N个点 以DFT为例： N (N – 1) N 2 N个X(k) (N点DFT) N – 1 N 一个X(k) 复数加法 复数乘法 2N+2 (N – 1)=2 (2N – 1) 4N 一个X (k) 2N (2N – 1) 4N 2 N个X (k) (N点DFT) 2 一次复加 2 4 一次复乘 实数加法 实数乘法 运算量 (a+jb)(c+jd)=(ac-bd)+j(bc+ad) 例：设做一次复乘用时1μs。石油勘探，有24个通道的记录 ，每通道波形记录长度为5秒，若每秒抽样500点/秒， 1）每道总抽样点数：500*5=2500点 2）24道总抽样点数：24*2500=6万点 3）DFT复乘运算时间：N2=(60000)2=36*108次 所需时间为 由于计算量大，且要求相当大的内存，难以实现实时处理，限制了DFT的应用。长期以来，人们一直在寻求一种能提高DFT运算速度的方法。 1、利用DFT运算的系数的固有对称性和周期 性，合并DFT中的某些项。 1）对称性 2）周期性 3）可约性 二、改善DFT运算效率的基本途径 2、将长序列DFT利用对称性和周期性分解为短 序列DFT来减少计算量。 因为DFT的运算量与N2成正比的，如果一个大点数N的DFT能分解为若干小点数DFT的组合，则显然可以达到减少运算工作量的效果。 N点 DFT N/2点 DFT N/2点 DFT N/4点 DFT N/4点 DFT N/4点 DFT N/4点 DFT ……. 复乘： FFT算法分类: 时域抽取法 DIT: Decimation-