DSP第五章快速付里叶变换FFT.ppt

第四章快速付里叶变换(FFT)Fast FourierTransforming 本章主要内容 1.直接计算DFT算法存在的问题及改进途径。 2.多种DFT算法(时间抽取算法DIT算法，频率抽取算法DIF算法，线性调频Z变换即CZT法) 3.FFT的应用 二、FFT产生故事 当时加文(Garwin)在自已的研究中极需要一个计算付里叶变换的快速方法。他注意到图基(J.W.Turkey)正在写有关付里叶变换的文章，因此详细询问了图基关于计算付里叶变换的技术知识。图基概括地对加文介绍了一种方法，它实质上就是后来的著名的库利(Cooley J.W)图基算法。在加文的迫切要求下，库利很快设计出一个计算机程序。1965年库利--图基在计算数学、Mathematic of Computation 杂志上发表了著名的“机器计算付里级数的一种算法”文章，提出一种快速计算DFT的方法和计算机程序--揭开了FFT发展史上的第一页，促使FFT算法产生原因还有1967年至1968年间FFT的数字硬件制成，电子数字计算机的条件， 使DFT的运算大简化了。 §第1节 直接计算DFT算法存在的问题及改进途径 一、直接计算DFT计算量 问题提出： 设有限长序列x(n),非零值长度为N,计算对x(n)进行一次DFT运算，共需多大的运算工作量? 2.以DFT为例，计算DFT复数运算量 计算一个X(k)(一个频率成分)值，例k=1则 运算量为N次复数乘法+(N-1)次复数加法。所以，要完成整个DFT运算，其计算量为： N*N次复数相乘和N*(N-1)次复数加法 4.计算DFT需要的实数运算量 每运算一个X(k)的值，需要进行4N次实数相乘和 2N+2(N-1)=2(2N-1)次实数相加. 整个DFT运算量为：4N2次实数相乘和2N(2N-1)次实数相加. 由此看出：直接计算DFT时，乘法次数与加法次数都是和N2成比例的。当N很大时，所需工作量非常可观。 例子 例1:当N=1024点时，直接计算DFT需要： N2=220=1048576次，即一百多万次的复乘运算 这对实时性很强的信号处理(如雷达信号处理)来讲，它对计算速度有十分苛刻的要求--迫切需要改进DFT的计算方法，以减少总的运算次数。 例2：石油勘探，24道记录，每道波形记录长度5秒，若每秒抽样500点/秒， 每道总抽样点数=500*5=2500点 24道总抽样点数=24*2500=6万点 DFT运算时间=N2=(60000)2=36*108次 二、改善DFT运算效率的基本途径 利用DFT运算的系数 的固有对称性和周期性，改善DFT的运算效率。 1.合并法：合并DFT运算中的某些项。 2. 分解法：将长序列DFT利用对称性和周期性，分解为短序列DFT。 利用DFT运算的系数 的固有对称性和周期性，改善DFT的运算效率 例子 例: 1 合并法：步骤1分解成虚实部 合并DFT运算中的有些项 对虚实部而言 所以代入DFT中： 合并法：步骤3合并有些项 2、将长序列DFT利用对称性和周期性分解为短序列DFT--方法 结论 快速付里时变换(FFT)就是在此特性基础上发展起来的，并产生了多种FFT算法，其基本上可分成两大类： 按抽取方法分： 时间抽取法(DIT);频率抽取法(DIF) 按“基数”分：基-2FFT算法；基-4FFT算法；混合基FFT算法；分裂基FFT算法 其它方法：线性调频Z变换(CZT法) 第二节基--2按时间抽取的FFT算(Coolkey-Tukey) 一、算法原理 设输入序列长度为N=2M(M为正整数，将该序列按时间顺序的奇偶分解为越来越短的子序列，称为基2按时间抽取的FFT算法。也称为Coolkey-Tukey算法。 其中基数2----N=2M，M为整数.若不满足这个条件，可以人为地加上若干零值（加零补长）使其达到 N=2M 例子 设一序列x(n)的长度为L=9,应加零补长为 N=24=16 应补7个零值。 二、算法步骤1.分组，变量置换 2.代入DFT中 3.求出子序列的DFT 4.结论1 一个N点的DFT被分解为两个N/2点DFT。X1(k),X2(k)这两个N/2点的DFT按照： 再应用W系数的周期性，求出用X1(k),X2(k)表达的后半部的X(k+N/2)的值。 5.求出后半部的表示式 6.结论2 三、蝶形结 例子：求 N=23=8点FFT变换 （1）先按N=8--N/2=4