学习笔记：快速傅里叶变换（FFT）

引言

  滴！解锁研究生关卡！博主入学已经一个学期了，万万没想到迎来的第一个课题就是与FFT有关的。要知道，本科的时候数字信号处理课程中的FFT、蝶形图就是本人的噩梦，虽然它将对称美发挥到了极致，但是！他喵的我就是学不懂嘛呜呜呜呜。怕什么还真就来什么，没辙，和他交个朋友吧。

一、DFT（Discrete Fouier Transform）

  FFT本质上是DFT的一种快速算法，所以想弄清楚其原理，就必须先明白DFT的内涵。DFT，即离散傅里叶变换，用来将时域序列转换为频域序列，一个N点的DFT变换具体公式如(1-1)所示：
$$X(k)=\sum _{n=0}^{N-1}x(n)W_N^{nk}n,k\in \left[0,N-1\right]\text{\hspace{1em}(1-1)}$$
  在FFT中，常取$N=2^M$。下面以8点FFT为例，讲解FFT时域抽取和频域抽取的两种算法的基本原理。

二、时域抽取（DIT）

2.1原理

  时域抽取，一般来讲是将序列x[n]按照奇偶划分为两部分，具体分解步骤如下：

step1. 将DFT变换公式按时域序列n的奇偶拆写成两个部分：
$$\begin{array}{rl}X(k)& =\sum _{n=0}^{N-1}x(n)W_N^{nk}\\ & =\sum _{n=0}^{N/2-1}x(2r)W_N^{2rk}+\sum _{n=0}^{N/2-1}x(2r+1)W_N^{2rk}{W}_N^k\\ & =\sum _{n=0}^{N/2-1}{x}_1(r)W_{N/2}^{rk}+\sum _{n=0}^{N/2-1}{x}_2(r)W_{N/2}^{rk}{W}_N^k\\ & =X_1(k)+X_2(k)W_N^{k}k=0,1,...,N-1\end{array}$$
其中：
$$\begin{array}{r}{X}_1(k)=\sum _{n=0}^{N/2-1}{x}_1(r)W_{N/2}^{rk}=DFT[x_1(r)]k=0,1,...,N-1\\ X_2(k)=\sum _{n=0}^{N/2-1}{x}_2(r)W_{N/2}^{rk}=DFT[x_2(r)]k=0,1,...,N-1\end{array}$$
$X_1(k)$和$X_2(k)$是以N/2为周期的DFT序列，$W_N^{k+N/2}=-W_N^k$，所以频域公式还可以写成：
$$\begin{array}{rl}X(k)& =X_1(k)+X_2(k)W_N^{k}k=0,1,...,N/2-1\\ X(k+N/2)& =X_1(k)-X_2(k)W_N^{k}k=0,1,...,N/2-1\end{array}$$
蝶形图：

step 2.将$x_1(r)$按r的奇偶进行分解：
$$\begin{array}{rl}{X}_1(k)& =\sum _{n=0}^{N/2-1}{x}_1(r)W_{N/2}^{nk}\\ & =\sum _{n=0}^{N/4-1}{x}_1(2l)W_{N/2}^{2lk}+\sum _{n=0}^{N/4-1}{x}_1(2l+1)W_{N/2}^{2lk}{W}_{N/2}^k\\ & =\sum _{n=0}^{N/4-1}{x}_3(r)W_{N/4}^{lk}+\sum _{n=0}^{N/4-1}{x}_4(r)W_{N/4}^{lk}{W}_{N/2}^k\\ & =X_3(k)+X_4(k)W_{N/2}^{k}k=0,1,...,N/2-1\end{array}$$
$X_3(k)$和$X_4(k)$都是以N/4为周期的序列，即$W_{N/2}^{k+N/4}=-W_{N/2}^k$，所以：
$$\begin{array}{rl}{X}_1(k)& =X_3(k)+X_4(k)W_{N/2}^{k}k=0,1,...,N/4-1\\ X_1(k+N/4)& =X_3(k)-X_4(k)W_{N/2}^{k}k=0,1,...,N/4-1\end{array}$$
同理：
$$\begin{array}{rl}{X}_2(k)& =X_5(k)+X_6(k)W_{N/2}^{k}k=0,1,...,N/4-1\\ X_2(k+N/4)& =X_5(k)-X_6(k)W_{N/2}^{k}k=0,1,...,N/4-1\end{array}$$

2.2 8点FFT举例

将DFT变换公式按时域序列n的奇偶拆写成两个部分：
x ( n ) = { x ( 0 ) , x ( 1 ) , x ( 2 ) , x ( 3 ) , x ( 4 ) , x ( 5 ) , x ( 6 ) , x ( 7 ) } x(n) =\{x(0),x(1),x(2),x(3),x(4),x(5),x(6),x(7)\} x(n)={x(0),x(1),x(2),x(3),x(4),x(5),x(6),x(7)}

x 1 ( n ) = { x 1 ( 0 ) , x 1 ( 1 ) , x 1 ( 2 ) , x 1 ( 3 ) } = { x ( 0 ) , x ( 2 ) , x ( 4 ) , x ( 6 ) } x_1(n) =\{x_1(0),x_1(1),x_1(2),x_1(3)\}=\{x(0),x(2),x(4),x(6)\} x1​(n)={x1​(0),x1​(1),x1​(2),x1​(3)}={x(0),x(2),x(4),x(6)}
x 2 ( n ) = { x 2 ( 0 ) , x 2 ( 1 ) , x 2 ( 2 ) , x 2 ( 3 ) } = { x ( 1 ) , x ( 3 ) , x ( 5 ) , x ( 7 ) } x_2(n) =\{x_2(0),x_2(1),x_2(2),x_2(3)\}=\{x(1),x(3),x(5),x(7)\} x2​(n)={x2​(0),x2​(1),x2​(2),x2​(3)}={x(1),x(3),x(5),x(7)}

x 1 ( n ) = { x 1 ( 0 ) , x 1 ( 1 ) , x 1 ( 2 ) , x 1 ( 3 ) } = { x ( 0 ) , x ( 2 ) , x ( 4 ) , x ( 6 ) } x_1(n) =\{x_1(0),x_1(1),x_1(2),x_1(3)\}=\{x(0),x(2),x(4),x(6)\} x1​(n)={x1​(0),x1​(1),x1​(2),x1​(3)}={x(0),x(2),x(4),x(6)}

• x 3 ( n ) = { x 3 ( 0 ) , x 3 ( 1 ) } = { x ( 0 ) , x ( 4 ) } x_3(n) =\{x_3(0),x_3(1)\}=\{x(0),x(4)\} x3​(n)={x3​(0),x3​(1)}={x(0),x(4)}
  • x 7 ( n ) = { x ( 0 ) } x_7(n) =\{x(0)\} x7​(n)={x(0)}
  • x 8 ( n ) = { x ( 4 ) } x_8(n) =\{x(4)\} x8​(n)={x(4)}

$X_7(0)=x_3(0)=x(0)$
$X_8(0)=x_3(1)=x(4)$
$X_3(0)=X_7(0)=x_4(0)=x(2)$
$X_3(1)=X_8(0)=x_4(1)=x(6)$

• x 4 ( n ) = { x 4 ( 0 ) , x 4 ( 1 ) } = { x ( 2 ) , x ( 6 ) } x_4(n) =\{x_4(0),x_4(1)\}=\{x(2),x(6)\} x4​(n)={x4​(0),x4​(1)}={x(2),x(6)}
  • x 9 ( n ) = { x ( 2 ) } x_9(n) =\{x(2)\} x9​(n)={x(2)}
  • x 10 ( n ) = { x ( 6 ) } x_{10}(n) =\{x(6)\} x10​(n)={x(6)}

$X_9(0)=x_4(0)=x(2)$
$X_{10}(0)=x_4(1)=x(6)$
$X_4(0)=X_9(0)=x_4(0)=x(2)$
$X_4(1)=X_{10}(0)=x_4(1)=x(6)$

x 2 ( n ) = { x 2 ( 0 ) , x 2 ( 1 ) , x 2 ( 2 ) , x 2 ( 3 ) } = { x ( 1 ) , x ( 3 ) , x ( 5 ) , x ( 7 ) } x_2(n) =\{x_2(0),x_2(1),x_2(2),x_2(3)\}=\{x(1),x(3),x(5),x(7)\} x2​(n)={x2​(0),x2​(1),x2​(2),x2​(3)}={x(1),x(3),x(5),x(7)}

• x 5 ( n ) = { x 5 ( 0 ) , x 5 ( 1 ) } = { x ( 1 ) , x ( 5 ) } x_5(n) =\{x_5(0),x_5(1)\}=\{x(1),x(5)\} x5​(n)={x5​(0),x5​(1)}={x(1),x(5)}
  • x 11 ( n ) = { x ( 1 ) } x_{11}(n) =\{x(1)\} x11​(n)={x(1)}
  • x 12 ( n ) = { x ( 5 ) } x_{12}(n) =\{x(5)\} x12​(n)={x(5)}

$X_{11}(0)=x_5(0)=x(1)$
$X_{12}(0)=x_5(1)=x(5)$
$X_5(0)=X_{11}(0)=x_5(0)=x(1)$
$X_5(1)=X_{12}(0)=x_5(1)=x(5)$

• x 6 ( n ) = { x 6 ( 0 ) , x 6 ( 1 ) } = { x ( 3 ) , x ( 7 ) } x_6(n) =\{x_6(0),x_6(1)\}=\{x(3),x(7)\} x6​(n)={x6​(0),x6​(1)}={x(3),x(7)}
  • x 13 ( n ) = { x ( 3 ) } x_{13}(n) =\{x(3)\} x13​(n)={x(3)}
  • x 14 ( n ) = { x ( 7 ) } x_{14}(n) =\{x(7)\} x14​(n)={x(7)}

$X_{13}(0)=x_6(0)=x(3)$
$X_{14}(0)=x_6(1)=x(7)$
$X_6(0)=X_{13}(0)=x_6(0)=x(3)$
$X_6(1)=X_{14}(0)=x_6(1)=x(7)$
由于博主电脑上没有绘制蝶形图的软件，所以直接用了之前一篇博文的配图，原博的连接如下link，这篇文章讲述的十分详细，作者的博客还有专栏来介绍FFT和蝶形图，推荐大家阅读。

三、频域抽取(DIF)

1. 将DFT变换公式的时域序列从中间分成前、后两个部分：
   $$\begin{array}{rl}X(k)& =\sum _{n=0}^{N-1}x(n)W_N^{nk}\\ & =\sum _{n=0}^{N/2-1}x(n)W_N^{nk}+\sum _{n=N/2}^{N}x(n)W_N^{nk}\\ & =\sum _{n=0}^{N/2-1}x(n)W_N^{nk}+\sum _{n=0}^{N/2-1}x(n+N/2)W_N^{(n+N/2)k}\\ & =\sum _{n=0}^{N/2-1}\left[x(n)+x(n+N/2)W_N^{(N/2)k}\right]W_N^{nk}\\ & =\sum _{n=0}^{N/2-1}\left[x(n)+(-1{)}^{k}x(n+N/2)\right]W_N^{nk}\end{array}$$
   蝶形图:
   

2. 频域按k的奇偶划分成两个部分，可以写作下式：
   $$\begin{array}{rl}X(2m)& =\sum _{n=0}^{N/2-1}x\left[(n)+(-1{)}^{2m}x(n+N/2)\right]W_N^{2mn}\\ & =\sum _{n=0}^{N/2-1}\left[x(n)+x(n+N/2)\right]W_{N/2}^{mn}\\ & =\sum _{n=0}^{N/2-1}{x}_1(n)W_{N/2}^{mn}m=0,1,...,N/2\\ & =DFT[x_1(n)]=X_1(k)\\ X(2m+1)& =\sum _{n=0}^{N/2-1}\left[x(n)+(-1{)}^{(2m+1)}x(n+N/2)\right]W_N^{n(2m+1)}\\ & =\sum _{n=0}^{N/2-1}\left[x(n)-x(n+N/2)\right]W_N^n{W}_{N/2}^{mn}\\ & =\sum _{n=0}^{N/2-1}{x}_2(n)W_{N/2}^{mn}m=0,1,...,N/2\\ & =DFT[x_1(n)]=X_1(k)\end{array}$$
   其中：
   $$\begin{array}{rl}{x}_1(n)& =\left[x(n)+x(n+N/2)\right]\\ x_2(n)& =\left[x(n)-x(n+N/2)\right]W_N^n\end{array}$$
   3.重复上述步骤即可。