第4章 FFT变换.ppt

§4-1引言 一.DFT的计算工作量 两者的差别仅在指数的符号和因子1/N. 通常x(n)和 都是复数,所以计算一个 X(k)的值需要N次复数乘法运算,和 次 复数加法运算.那么,所有的X(k)就要N2次复 数乘法运算,N(N-1)次复数加法运算.当N很 大时,运算量将是惊人的,如N=1024,则要完 成1048576 次(一百多万次)运算.这样,难以做到实时处理. 二.改进的途径 1. 的对称性和周期性 利用上述特性,可以将有些项合并,并 将DFT分解为短序列,从而降低运算次数,提 高运算速度.1965年,库利(cooley)和图基 (Tukey)首先提出FFT算法.对于N点DFT,仅需 (N/2)log2N 次复数乘法运算.例如N=1024=210 时， 需要(1024/2)log2 210 =512*10=5120次。 5120/1048576=4.88% ,速度提高20倍 §4-2 按时间抽取(DIT)的FFT算法 —库利-图基算法 一.算法原理(基2FFT) (一)N/2点DFT 1.先将 按n的奇偶分为两组作DFT,设N=2L ,不足时,可补些零。这样有: n为偶数时: n为奇数时: 由于: 所以,上式可表示为: 其中, 2.N/2点DFT的结论: (1) X (k),X (k)均为N/2点的DFT。 (2) X(k)=X (k)+W X (k)只能确定出 X(k)的k= 个； 即前一半的结果。 $4.5 DFT和 FFT的应用 FFT频谱分析 FFT线性卷积计算 $4.5.1 用FFT 进行频谱分析 一、方法 有限长序列x(n)的在 之间的频谱： 幅度谱： 相位谱： 频率分辨率： 二、存在问题 1.混叠现象 2.频谱泄露 3.栅栏效应 $4.5.3 用FFT计算线性相关 1. 循环相关与线性相关的关系 线性相关： rxy取非零值长度为： 规律： 1、长度分别为M和N的序列的卷积序列的非零值长度为：M+L-1 2、长度分别为M和N的序列的相关序列的非零值长度为：M+L-1 3、用函数conv(x,y)计算卷积，产生卷积序列的长度：M+L-1 4、用函数xcorr(x,y)计算相关，产生相关序列的长度：2max(M,L)-1 $4.6 N为合数的DFT快速算法 $4.7 线性调频Z变换（Chirp-Z） 2.Chirp-Z实现步骤 第4章 FFT总结 一、DIT-FFT 将x(n)按n的奇偶分为两组作DFT 五、（一）用FFT 进行频谱分析 有限长序列x(n)的在 之间的频谱： 幅度谱： 相位谱： 频率分辨率： 本章结束 1.混叠现象 为避免混叠,由抽样定理可知，须满足 其中, 为抽样频率; 为信号的最高频率分量； 或者 其中,T为抽样间隔。 2.频谱泄漏 在实际应用中，通常将所观测的信号 限制在一定的时间间隔内,也 就是说，在时域对信号进行截断操作,或 称作加时间窗,亦即用时间窗函数乘以信号,由卷积定理可知，时域相乘,频域为卷积,这就造成拖尾现象,称之为频谱泄漏. 0 n 0 n n 3.栅栏效应 用DFT计算频谱时，只是知道为频率 的整数倍处的频谱。 在两个谱线之间的情况就不知道,这相当通过一个栅栏观察景象一样，故称作栅栏效应。 补零点加大周期，可使F变小来提高 辨力，以减少栅栏效应。 §4.5.2 线性卷积的FFT算法 一.线性卷积的长度 设一离散线性移不变系统的冲激响 应为 ,其输入信号为 .其输出 为 .并且 的长度为L点, 的 长度为M点,则： 以实例说明： 0 1 2 3 1 2 3 . 。 1 。 0 1 1 2 3 2 3 。 。 0 1 2 3 。 0 1 1 2 3 2 3 。 0 1 2 3 4 0 1 2 3 4 5 0 1 1 2 3 2 3 。 。 0 1 2 3 4 5 1 3 3 5 6 6 。 二.用FFT算 的步骤： FFT FFT IFFT x x(n) h(n) y(n) X(k