4.1-4.2快速傅里叶变换.pptVIP

第四章 快速傅里叶变换（FFT） §4.1 基2FFT算法 运算量比较 DIT-FFT的运算规律及编程思想 算法的推导 DIF-FFT的蝶式运算流图 DIF-FFT的二次分解运算流图 DIF-FFT的运算流图 (N=8) DIF-FFT与DIT-FFT比较 IFFT算法的推导 利用 FFT 进行频谱分析 用FFT计算线性卷积 一、利用FFT进行频谱分析的基本方法 几个常用基本概念 典 型 例 题 用FFT进行频谱分析存在的两个问题 二、用FFT计算线性卷积 线性卷积 利用循环卷积计算线性卷积的条件 利用FFT进行线性卷积的步骤 线性调频Z变换(Chirp- Z 变换)算法 问题的引入 算法的基本原理 Chirp- Z 变换的方框图 说明 比较两式可知：只要将FFT中的旋转因子 改为 ， 再乘以 1/N 即可得到 IDFT 的快速算法 IFFT。 IFFT基本思想 ，∴还可将常数 1/N 分配到每级运算中，也就是每级蝶形运算均乘以 ? 。这样就实现了 FFT 与IFFT 运算的统一。 思考： 是否可以直接调用FFT子程序计算IFFT 直接调用FFT子程序计算IFFT的方法： §4.2 DFT 与 FFT 的应用 设x(n)为长为 N 的有限长序列，则： 利用 FFT 进行频谱分析的实现过程框图为： 1、数字频率分辨率： 2、模拟频率分辨率： 3、用于FFT的采样点数： 4、频率刻度值： 5、模拟信号长度(对信号的观察长度)： 6、分辨率： 例：用FFT来分析信号的频谱，若已知信号的最高频率为 ，要求频率分辨率为 ，试确定： 1、采样间隔 T ； 2、采用基-2FFT的最小样点数 N ，以及与此相对应的最小记录长度； 3、按您确定的参数所获得的实际分辨率。 解： 1、据采样定理，采样间隔 2、基-2FFT的最小样点数N 当采用基-2FFT算法时，要求 与此相对应的最小记录长度为： 3、按确定的参数所获得的实际分辨率 若想提高一倍的分辨率，怎么办？ 利用DFT分析模拟信号的频谱 用DFT分析模拟信号的频谱：计算信号的傅里叶变换 1、频谱泄漏 若序列 的长度为无限长，为了利用 FFT 进行频谱分析，首先必须将其截断为有限长序列 卷积定理 显然，两种频谱是有差别的，该现象就是频谱泄漏 2、栅栏效应 利用 FFT 进行频谱分析时，只知道离散频率点 的整数倍处的频谱。在两个谱线之间的情况就不知道，这如同通过一个栅栏观察景象一样，故称作栅栏效应。 解决办法：在序列后面补零点加大FFT点数 ，可使谱线间隔变小来提高分辨力，以减少栅栏效应。 注意：若需要加窗，则应先加窗再补零。 设 是 和 的线性卷积： 总运算量为： 思路：利用 FFT ,通过循环卷积来计算线性卷积 设 是x(n)和h(n)长为L的循环卷积： 其中 LMax[N,M], 上式表明：利用循环卷积计算线性卷积的条件为： 利用循环卷积计算线性卷积如下图 ①、将已知序列 （长为N）和 （长为M）补零延长，使它们的长度 。若采用基 -2 FFT算法，还应使 大于或等于 的 2 的最小整数次幂。 ②、做 和 的长为 点的 FFT 得到 和 ，并求它们的积 。 ③、求 的 IFFT并取前 点获得线性卷积的结果为 尽管采用 FFT 可计算出有限长序列的DFT ，但它要求序列长度 为2的整数幂或合数。实际中①有时只对信号的某一频段感兴趣，即只需要计算单位圆上某一段的频谱值，例如对窄带信号进行频谱分析时，总是要求在窄带范围内的抽样点足够密集，而窄带范围外则不需考虑。此时若依然采样以上方法，则需增加频域抽样点数，增加了窄带范围外的不需要的计算量；②在语声信号处理时，信号极点处的频谱十分关键，而极点位置往往离单位圆较远，此时不能采用FFT；③若 是大素数时也无法采用FFT。 上述几种情形可采用线性调频Z变换(Chirp- Z 变换)算法加以解决。 线性调频Z变换(Chirp- Z 变换)算法 设 是有限长序列， 沿Z平面上的一段螺线