第5章快速傅立叶变换(FFT)讲解.pptVIP

第五章快速傅立叶变换（FFT） 主要内容 5.1 引言 5.2 基2FFT的时间抽取算法及蝶形图 5.3 基2FFT的频率抽取算法 5.4 进一步减少计算量的措施 5.5 基4FFT的频率抽取算法 5.6 分裂基FFT算法 5.7 多种形式的FFT算法 5.8 结束语 5.1 引言 (1)离散傅立叶变换（DFT） (2)意义 ①DFT是离散谱。②从理论上解决了离散谱的计算问题，具有重要的理论意义。 (3)计算量和存储量分析 ①计算量用复数加法和复数乘法次数表示。与N2成正比。 复数加法次数：N(N-1)次。复数乘法次数：N2次。 ②存储量用占用随机存取存储器和只读存储器的数量表示。 RAM单元：N+N=2N个。 ROM单元：N2个。 结论：计算量和存储量，随着数据量的增加都迅速增大。应用困难。 (4)解决方法－快速算法 J.W.Cooley (IBM), J.W.Tukey(MIT), “机器计算傅立叶级数的一种算法”，Math. Computation, Vol.19,1965.称为 Cooley and Turkey 算法。 以后，出现Sande Tukey算法，P.Dohamel and H.Hollmann的分裂基算法。 5.2 基2FFT的时间抽取算法及蝶形图（1） (1) Cooley and Turkey 算法的基本思想 W因子的特性： 基本思想：分解。 分解的唯一条件：数据长度满足N=2r。否则，补零解决。 (2)分解公式的推导(按奇偶序号分解，输入每2点抽1个值) 5.2 基2FFT的时间抽取算法及蝶形图（2） (3) FFT定理或分解公式 (4) 蝶形图和一次分解图：用蝶形图计算DFT。运算量:一次乘，两次加。 一次分解的计算量 5.2 基2FFT的时间抽取算法及蝶形图（3） (5) 8点数据时间抽取(Decimation-In-Time)算法的完整分解过程和蝶形流程图 蝶形图（蝶形运算流图）结构：N/2行，M=log2 N 列。 5.2 基2FFT的时间抽取算法及蝶形图（4） (6)时间抽取算法的特点 ①命名：Cooley-Turkey算法，时间抽取算法（简称DIT-FFT: Decimation In Time FFT） ② DIT-FFT的运算量（结构：N/2行，M=log2 N 列） ③ 旋转因子(twiddle factor)变化规律 5.2 基2FFT的时间抽取算法及蝶形图（5） ④原址计算（同址计算：In place computation ）：指输入数据、中间、最终结果可存放在同一存储单元。 ⑤序列的整序（倒序） 整序：字位倒置规则（bit reversed） 硬件电路和汇编程序用2进制倒序， 而高级语言则要换算成10进制。 5.2 基2FFT的时间抽取算法及蝶形图（6） (7) 16点数据时间抽取算法的完整信号流图。 5.3 基2FFT的频率抽取算法（1） (1) 基2FFT的频率抽取(Decimation-In-Frequency)算法的推导 5.3 基2FFT的频率抽取算法（2） (2)分解公式 (3)蝶形图、蝶形运算：两个短序列x1(n)、x2(n)可成对运算。 蝶形运算符：一次乘法，两次加法。先加减，后乘。而DIT-FFT：先乘，后加减。 5.3 基2FFT的频率抽取算法（3） (4)一次分解图(n=0,1,2,3) 一次分解的计算量与时间抽取算法相同。 (5)基本思想 长序列分解成短序列。从过程看，将输入数据对半分开为两个序列。从结果看，将DFT的输出X(k)按照奇、偶序号分解为两个序列，每2点抽1个值。 命名：基2FFT的频率抽取算法（Sand-Turkey算法，简称DIF-FFT: Decimation In Frequency FFT） 5.3 基2FFT的频率抽取算法（4） (6) 8点数据频率抽取算法的完整分解过程和蝶形流程图 ①完整分解过程 5.3 基2FFT的频率抽取算法（5） ②蝶形流程图（结构：N/2行，M=log2 N 列） (7)运算特点：运算量、旋转因子规律、同址运算、存储量、整序规则 整序规则：字位倒置规则。 (8)与时间抽取算法比较 分解思想，计算量存储量，整序规则相同。分解方