信号与系统复习资料快速傅里叶变换FFT.pptVIP

第四章快速傅里叶变换（FFT） 主要内容 DIT-FFT算法 DIF-FFT算法* IFFT算法 §4.1 引言 FFT: Fast Fourier Transform 1965年，Cooley-Turky 发表文章《机器计算傅里叶级数的一种算法》，提出FFT算法，解决DFT运算量太大，在实际使用中受限制的问题。 FFT的应用: 频谱分析、滤波器实现、实时信号处理等。 DSP芯片实现。TI公司的TMS 320C30，10MHz时钟，基2-FFT 1024点FFT时间15ms。 典型应用：信号频谱计算、系统分析等 §4.2 直接计算DFT的问题及改进途径 3、降低DFT运算量的考虑 §4.3 按时间抽取（DIT）的FFT算法 1、算法原理 设序列点数 N = 2L，L 为整数。 若不满足，则补零 将N点DFT定义式分解为两个长度为N/2的DFT 将上式表达的运算用 “蝶形”信流图表示。 用“蝶形结”表示上面运算的分解： N点DFT分解为四个N/4点的DFT 类似进一步分解 类似的分解一直继续下去，直到分解为最后的两点蝶形运算为止(2点DFT). 如上述N=8=23，N/4=2点中： FFT运算量与运算特点 1． N=2L时，共有L=log2N级运算；每一级有N/2个蝶形结。 2．每一级有N个输入数据和N个输出数据，且每级输出 只与输入有关，前级输出为后级输入 (即原位计算迭代，节省存储单元） 。 3．计算量： 每级N/2次复乘法，（每蝶形只乘一次）。共有L*N/2=(N/2)log2N 次复乘法。 按时间抽取FFT蝶形运算特点 1、关于FFT运算的顺序处理 由于时间序列按时间序位的奇偶抽取，故输入序列是混序的，为此需要先进行混序处理。 混序规律： x(n)按n的二进制码倒置排列 位倒序实现： DSP实现采用位倒序寻址. 例 计算 ， 。计算 点FFT。用时间抽取输入倒序算法，问倒序前寄存器的数 和倒序后 的数据值？ DIT FFT中最主要的蝶形运算实现 （1）参与蝶形运算的两个结点（信号）间“距离”（码地址）与其所处的第几级蝶形有关；第m级的“结距离”为 （2）每级迭形结构为 DIT算法的其他形式流图 输入倒位序输出自然序 输入自然序输出倒位序 输入输出均自然序 相同几何形状 输入倒位序输出自然序 输入自然序输出倒位序 §4.4 按频率抽取（DIF）的FFT算法 §4.5 IDFT的FFT算法（FFT应用一） 一、定义比较分析 离散时间信号与系统 Z变换与离散时间傅立叶变换(DTFT) 离散傅立叶变换(DFT) 快速傅立叶变换 数字滤波器的基本结构 无限单位冲激响（IIR)应数字滤波器设计 有限单位冲激响（FIR)应数字滤波器设计 频谱分析 数字滤波器 离散时间信号与系统 Z变换与离散时间傅立叶变换(DTFT) 离散傅立叶变换(DFT) 快速傅立叶变换 数字滤波器的基本结构 无限单位冲激响（IIR)应数字滤波器设计 有限单位冲激响（FIR)应数字滤波器设计 FFT * -*- 反回 系统分析 1、 DFT与IDFT 2N+2 (N – 1)= 2(2N – 1) 4N 一个X (k) 2N (2N – 1) 4N 2 N个X (k) (N点DFT) 2 一次复加 2 4 一次复乘 实数加法 实数乘法 2、DFT与IDFT运算特点 N (N – 1) N 2 N个X(k) (N点DFT) N – 1 N 一个X(k) 复数加法 复数乘法 同理：IDFT运算量与DFT相同。 1= 1= FFT算法分类: 时间抽算法 DIT: Decimation-In-Time 频率抽算法 DIF: Decimation-In-Frequency 将序列x(n)按n的奇偶分成两组： N为2的整数幂的FFT算法称基-2FFT算法。 记： ………（1） 可约性 再利用周期性求X(k)的后半部分 注：a. 上支路为加法，下支路为减法； b. 乘法运算的支路标箭头和系数。 0+8/2 3+8/2 K=0 K=3 分解后的运算量： N (N/2 –1) N2/2 两个N/2点DFT N N/2 N/2个蝶形 总计 2 1 一个蝶形 N/2 (N/2 –1) (N/2)2 一个N/2点DFT 复数加法 复数乘法 运算量减少了近一半 由于 ，