第三章离散傅里叶变换和快速付里叶变换FFT2课题.pptVIP

3.4 快速傅立叶变换 FFT产生 1965年，库利--图基在《计算数学》（Mathematic of Computation）杂志上发表了著名的“机器计算傅里级数的一种算法”文章，提出一种快速计算DFT的方法和计算机程序--揭开了FFT发展史上的第一页。 本节主要内容 ? 直接计算DFT算法存在的问题及改进途径。 ? FFT算法(基－2时间抽取算法） ? FFT的应用 一、直接计算DFT算法存在的问题及改进途径 直接计算DFT计算量 问题提出：设有限长序列x(n)，非零值长度为N，计算对x(n)进行一次DFT运算，共需多大的运算工作量? (1) 以DFT为例，计算DFT复数运算量 计算一个X(k)(一个频率成分)值，运算量为 例k=1则 要进行N次复数乘法+(N-1)次复数加法 所以，要完成整个DFT运算，其计算量为： N*N次复数相乘+N*(N-1)次复数加法 (2) 一次复数乘法换算成实数运算量 复数运算要比加法运算复杂，需要的运算时间长。 一个复数乘法包括4个实数乘法和2个实数相法。 (a+jb)(c+jd)=(ac-bd)+j(bc+ad) (3) 计算DFT需要的实数运算量 每运算一个X(k)的值，需要进行 4N次实数相乘和 2N+2(N-1)=2(2N-1)次实数相加. 整个DFT运算量为：4N2次实数相乘和2N(2N-1)次实数相加. 由此看出：直接计算DFT时，乘法次数与加法次数都是和N2成比例的。当N很大时，所需工作量非常可观。 当N=1024点时，直接计算DFT需要： N2=220=1048576次，即一百多万次的复乘运算。 这对实时性很强的信号处理(如雷达信号处理)来讲，它对计算速度有十分苛刻的要求 迫切需要改进DFT的计算方法，以减少总的运算次数。 (4) 比较DFT与IDFT之间的运算量 2. 改善DFT运算效率的基本途径 利用DFT运算的系数 的固有对称性和周期性，改善DFT的运算效率。 ? 合并法：合并DFT运算中的某些项。 ? 分解法：将长序列DFT利用对称性和周 期性，分解为短序列DFT。 利用DFT运算的系数 的固有对称性和周期性，改善DFT的运算效率 例: 将长序的DFT利用对称性和周期性分解为短序列DFT DFT的运算量与N2成正比 如果一个大点数N的DFT能分解为若干小点数DFT的组合，则显然可以达到减少运算工作量的效果。 方 法 结 论 快速付里时变换(FFT)就是在此特性基础上发展起来的，并产生了多种FFT算法，其基本上可分成两大类： 按抽取方法分： 时间抽取法(DIT)；频率抽取法(DIF) 按“基数”分：基-2FFT算法；基-4FFT算法；混合基FFT算法；分裂基FFT算法 其它方法：线性调频Z变换(CZT法) 二、基--2按时间抽取的FFT算法Decimation-in-Time(DIT) 1 算法原理 设输入序列长度为N=2M(M为正整数，将该序列按时间顺序的奇偶分解为越来越短的子序列，称为基2按时间抽取的FFT算法。也称为Coolkey-Tukey算法。 其中基数2----N=2M，M为整数。若不满足这个条件，可以人为地加上若干零值（加零补长）使其达到 N=2M。 例子 设一序列x(n)的长度为L=9，应加零补长为N=24=16 应补7个零值。 2 算法推导（1）分组 (2) 代入DFT中 (3) 求出子序列的DFT (4) 结论1 一个N点的DFT被分解为两个N/2点DFT。X1(k),X2(k)这两个N/2点的DFT按照： （5）求出后半部的表示式 (6) 结论2 (7）结论3 3 蝶形结（基本单元） 将N=8点分解成2个4点的DFT的信号流图 一个2点的DFT蝶形流图 另一个2点的DFT蝶形流图 将N/4（2点）DFT再分解成2个1点的DFT 2个1点的DFT蝶形流图 一个完整N=8的按时间抽取FFT的运算流图 4 FFT算法中一些概念 将N 点DFT先分成两个N/2点DFT，再是四个N/4点DFT…直至N/2个两点DFT.每分一次称为“一”级运算。 因为N=2M所以N点DFT可分成M级 如上图所示依次m=1,m=2….M共M级 （2）组 (3) 因子的分布 N=8点和N=1024点时直接计算DFT与用基2-按时间抽取法FFT的运算量。 5 按时间抽取FFT算法的特点 根据DIT基2-FFT算法原理，能得出任何