【数字信号处理】一文讲清FFT(快速傅里叶变换)

快速傅里叶变换（Fast Fourier Transform，FFT）

FFT的背景

• 1、为什么要时域→频域频率？50Hz+频率120Hz正弦波信号叠加，再叠加一些随时间变换的干扰信号（噪声），信号时域上的波形如图2- 1（左）所示。
  显然，从图中很难看出真正的有效信息，更无法直接过滤干扰信号。若对叠加信号进行离散傅里叶变换处理，将时域波形→频域图，如图2- 1（右）为变换后的频谱图。此时便很容易分辨出真正有效的频率段是在50Hz和120Hz，而其它频率都是干扰信息。
  
• 2、那么问题来了，怎样将时域→频域呢？
  于是出现了大神傅里叶。他提出了——离散傅里叶变换（Discrete Fourier Transform，DFT）是一种基础的数字信号处理算法，它将一个离散的信号从时域转换到频域，使得我们可以分析复杂信号中的频率成分。
  
• ？？？傅里叶变换是啥玩意？用人话来说，“把信号比喻成一件衣服，任何一块布都可以由不同颜色的线织成，远看这张比布我们不知道是由哪些颜色的线织成的，但如果把布料拆散成一条条线，是不是就可以清晰地分析出所用线的颜色”
• 同理，傅里叶认为任何一个信号都可以用很多个不同幅度、频率的正弦波/余弦波叠加而成，如果把这些“正弦波/余弦波”分离出来，就可以分析这些波的频率。如果再把这些频率画在如图2- 1（右）的频谱分布图上，就直观地实现了时域→频域

傅里叶变换（法语：Transformation de Fourier，英语：Fourier transform，缩写：FT）是一种线性变换，通常定义为一种积分变换。
其基本思想是一个函数可以用（可数或不可数，可数的情况对应于傅里叶级数）无穷多个周期函数的线性组合来逼近，从而这些组合系数在保有原函数的几乎全部信息的同时，还直接地反映了该函数的“频域特征”。——维基百科

• 3、说的简单，那我算起来该怎么算呢？
  在现实世界中，一条连续的信号需要经过采样，变成一个个离散的信号后，才能用计算机计算。
  用到的算法就是“离散傅里叶变换（Discrete Fourier Transform，DFT）”。
  但是离散傅里叶变换的计算量极为庞大，能算，但是太慢！记住就行，因为会用到很多的乘法、加法运算，计算机实现乘法很麻烦！

（有兴趣可以看看为什么计算量大）一个N点复数序列x[n]的离散傅里叶变换定义为公式（3-1）所示，在计算N点序列x[n]的DFT时，每计算一次频率分量应进行N次复数乘法运算以及N-1次复数加法运算，故需要N2次复数乘法运算和(N-1)N次复数加法运算才能获取N个频率分量X[k][28]。由此看来，随着N值增加，信号处理中离散傅里叶变换的运算量会变得极为庞大，从而难以满足实时性的需求。

快速傅里叶变换（Fast Fourier Transform，FFT）

于是便出现了快速傅里叶变换（Fast Fourier Transform，FFT），说白了就是一类快速的计算离散傅里叶变换的方法。代表者是1965年Cooley(库利)和Tukey(图基)提出的FFT，可以将计算N点序列DFT的复乘次数降低到(N/2)log2N。

• FFT假设输入信号是周期性的，即它会将一个有限长度的序列视为一个无限长的周期信号。当我们对一个有限长度的信号进行FFT时，实际上是在假设这个信号会无限重复。这个过程称为周期延拓。对于一个正弦波信号，如果你截取了一个完整的周期并进行FFT，得到的频谱将是该正弦波的特征频率。
• 傅立叶变换（FFT）是一种算法，用于计算离散傅立叶变换（DFT）和其逆变换。
• DFT的核心思想是任何时域信号都可以表示为一系列不同频率的正弦波和余弦波的加权和

DFT的数学表达

DFT将一个长度为 (N) 的时域序列 (x(n)) 转换为一个等长的频域序列 (X(k))。每个频域样本 (X(k)) 是通过以下公式计算得到的：

$$X(k)=\sum _{n=0}^{N-1}$$