数字信号处理.pptxVIP

数字信号处理 7.有限离散傅氏变换（II） 2022-7-28 数字信号处理 第七章（2） 2 7.2快速傅氏变换（FFT) 有限离散傅氏变换 根据上式直接计算，工作量太大，特别当 比较大时， 以至于在实际中无法使用。因此，才出现了快速傅氏变 换（Fast Fourier Transform——FFT）。 A.FFT的原理与算法之一——时域分解FFT算法 2022-7-28 数字信号处理 第七章（2） 3 7.2快速傅氏变换（FFT) 设 ，有限离散信号 ，则其有限 离散频谱为 令 ，则 按上式直接计算 ，需要 次复数乘法和加法，计算 个 则需要 次复数乘法和加法，当 很大时，该计算量很大。 2022-7-28 数字信号处理 第七章（2） 4 7.2快速傅氏变换（FFT) 将信号 按下标偶数项和奇数项分成两 部分，即令 它们的有限离散频谱分别为 注意 和 是以 为周期的函数，这是由于 2022-7-28 数字信号处理 第七章（2） 5 7.2快速傅氏变换（FFT) 从而有 进一步有 若令 （当 ），则有 ，于是 2022-7-28 数字信号处理 第七章（2） 6 7.2快速傅氏变换（FFT) 上式告诉我们，计算有 项的有限离散频谱问题可以转 化为计算有 项的有限离散频谱问题。同样， 项的 计算又可以转化为 项的计算，一直下去， 最后化到 ，这就是快速傅氏变换的思想。 下面介绍计算的过程，取 ，即 则（见下页） 2022-7-28 数字信号处理 第七章（2） 7 7.2快速傅氏变换（FFT) 对于一般 ，对于 ，则 2022-7-28 数字信号处理 第七章（2） 8 7.2快速傅氏变换（FFT) 二进制逆序排列 表 原序信号和逆序信号 原序信号 二进 序号 二进逆序号 逆序信号 2022-7-28 数字信号处理 第七章（2） 9 7.2快速傅氏变换（FFT) 计算量： 大的循环共进行 次，根据公式 计算两个频谱值，需要2次加法和1次乘法。计算 个频 谱需 次加法和 次乘法。进行 次循环，共需 次 加法和 次乘法。注意到 ，因此乘法可以 减少，即每一组计算中，可减少1次乘法。若 次循环中，将 个数分成 组，共减少了 次乘法，所以 次循环，共可减少乘法数为 2022-7-28 数字信号处理 第七章（2） 10 7.2快速傅氏变换（FFT) 总计算量： 加法, 乘法 （ ） 时域分解FFT算法流程图（ ）见下页。 2022-7-28 数字信号处理 第七章（2） 11 7.2快速傅氏变换（FFT) 2022-7-28 数字信号处理 第七章（2） 12 7.2快速傅氏变换（FFT) B.FFT的原理与算法之二——频域分解FFT算法 原理：在频率域上，把一个有 项的离散频谱按奇、偶 序号分解为两个 项的离散频谱，并依次进行下去。 注意到 则有 2022-7-28 数字信号处理 第七章（2） 13 7.2快速傅氏变换（FFT) 再令 这样， 是 个信号 的频谱； 是 个信号 的频谱。 2022-7-28 数字信号处理 第七章（2） 14 7.2快速傅氏变换（FFT) 即将计算 项离散频谱的问题转化为计算 项离散频 谱的问题。转化公式是通过将离散频谱 按奇、偶序