16、快速傅里叶变换与频域离散卷积

快速傅里叶变换与频域离散卷积

1. 快速傅里叶变换（FFT）

1.1 FFT 概述

快速傅里叶变换（FFT）是一种用于快速计算离散傅里叶变换（DFT）的数值方法。除了舍入噪声（即数字计算机在表示数字和执行操作时因精度有限而产生的噪声）外，其结果与 DFT 相同。这种技术的有效性不容小觑，特别是在处理代表信号的样本数量相对较大时，能在实际应用中产生显著差异。由于 FFT 能大幅节省时间，因此在实际中很少使用方程（2.56）中的 DFT 矩阵计算方法。

1.2 FFT 的优势

FFT 方法的主要优势在于它消除了 DFT 矩阵计算中的许多冗余。实际上，FFT 算法涉及的基本运算（即加法和乘法）数量级为 O(NLogN)，而标准 DFT 的数量级为 O(N²)。

1.3 FFT 算法的基数要求

一些 FFT 算法，如经典的 Cooley 和 Tukey 方法，要求 N 的值为 2 的整数次幂，如 32、64、128 等，这种方法被称为基数为 2 的方法。也可以使用其他基数，如 4、5、6 等。在一般情况下，N 必须是基数的整数次幂。在实践中，可以通过使用最小的 k 值，使得向量大小小于 2^k，并将 DFT 输入向量 h 中未使用的位置用零填充来满足这一要求。然而，这种做法可能会导致不连续性，并由于吉布斯效应在恢复信号中引入振荡。为缓解这一问题，可以将向量未使用部分的前半部分填充为原始函数的最后一个值，后半部分填充为原始函数的第一个值（因为 DFT 意味着该向量表示的函数是周期性的）。

1.4 FFT 与数字信号处理器

FFT 具有简单性和较少的数值运算（主要涉及复数乘法）