【數字信號處理】【傅里葉分析】【FFT】快速傅里葉變換的完整公式推導

給出離散傅里葉變換DFT的公式：
$\begin{aligned} DFT[x(n)]=X(m)&=\sum_{n=0}^{N-1}x(n)exp(-j\frac{2\pi mn}{N}),m\in[0,N-1]\\ &=x(0)exp(-j\frac{2\pi m0}{N})+x(1)exp(-j\frac{2\pi m1}{N})+x(2)exp(-j\frac{2\pi m2}{N})+x(3)exp(-j\frac{2\pi m3}{N})+... \end{aligned}$

FFT將DFT公式分成奇和偶兩部分，即，將含$x(0),x(2),...$的項取出，作爲偶數項，將含$x(1),x(3),...$的項取出，作爲奇數項，有
$\begin{aligned} X(m)&=\sum_{n=0}^{\frac{N}{2}-1}x(2n)exp(-j\frac{2\pi m2n}{N})+\sum_{n=0}^{\frac{N}{2}-1}x(2n+1)exp(-j\frac{2\pi m(2n+1)}{N}),m\in[0,N-1]\\ &=\sum_{n=0}^{\frac{N}{2}-1}x(2n)exp(-j\frac{2\pi mn}{\frac{N}{2}})+\exp(-j\frac{2\pi m}{N})\sum_{n=0}^{\frac{N}{2}-1}x(2n+1)\exp(-j\frac{2\pi mn}{\frac{N}{2}}),m\in[0,\frac{N}{2}-1]\\ \end{aligned}$
可以看到經過變換後，$m$的取值範圍變成了一半。現在考慮另一半的m，有
$\begin{aligned} X(m+\frac{N}2)&=\sum_{n=0}^{\frac{N}{2}-1}x(2n)\exp(-j\frac{2\pi n(m+\frac{N}2)}{\frac{N}{2}})+\exp(-j\frac{2\pi (m+\frac{N}2)}{N})\sum_{n=0}^{\frac{N}{2}-1}x(2n+1)\exp(-j\frac{2\pi n(m+\frac{N}2)}{\frac{N}{2}}),m\in[0,\frac{N}{2}-1]\\ &=\sum_{n=0}^{\frac{N}{2}-1}x(2n)\exp(-j\frac{2\pi mn}{\frac{N}{2}})\exp(-j2\pi n)+\exp(-j\frac{2\pi m}{N})\exp(-j\pi)\sum_{n=0}^{\frac{N}{2}-1}x(2n+1)\exp(-j\frac{2\pi nm}{\frac{N}{2}})\exp(-j2\pi n),m\in[0,\frac{N}{2}-1]\\ &=\sum_{n=0}^{\frac{N}{2}-1}x(2n)exp(-j\frac{2\pi mn}{\frac{N}{2}})-\exp(-j\frac{2\pi m}{N})\sum_{n=0}^{\frac{N}{2}-1}x(2n+1)\exp(-j\frac{2\pi mn}{\frac{N}{2}}),m\in[0,\frac{N}{2}-1] \end{aligned}$
令
$W^\theta_N=\exp(-j\frac{2\pi \theta}{N})$則
$X(m)=\sum_{n=0}^{\frac{N}{2}-1}x(2n)W_{\frac{N}2}^{mn}+W_N^m\sum_{n=0}^{\frac{N}{2}-1}x(2n+1)W_{\frac{N}2}^{mn},m\in[0,\frac{N}{2}-1]\\ X(m+\frac{N}2)=\sum_{n=0}^{\frac{N}{2}-1}x(2n)W_{\frac{N}2}^{mn}-W_N^m\sum_{n=0}^{\frac{N}{2}-1}x(2n+1)W_{\frac{N}2}^{mn},m\in[0,\frac{N}{2}-1]\\$
注意到
$\sum_{n=0}^{\frac{N}{2}-1}x(2n)W_{\frac{N}2}^{mn} ①\\ \sum_{n=0}^{\frac{N}{2}-1}x(2n+1)W_{\frac{N}2}^{mn} ②$
①和②分別是$N$點序列$x(n)$中對偶數項與奇數項的$\frac N 2$點的DFT。

上述兩個式子$X(m)$和$X(m+\frac N2)$爲FFT的核心公式，即一個$N$點DFT可以分成2個$\frac N2$點的DFT，並且前$\frac N2$點DFT與後$\frac N2$點DFT之間，只有正負號的不同，所以可以使用前$\frac N2$點DFT公式中的元素，計算後$\frac N2$點DFT中的元素，使計算量降低。

根據該思路，可以將①和②繼續分解，分別求①中奇數項與偶數項的$\frac N 4$點DFT，和②中奇數項和偶數項的$\frac N 4$點DFT。這樣一直分解下去，直到分解成最簡單的2點DFT爲止。注意到，在FFT中，N點DFT會一分爲二，二分爲四，四分爲八，最後會分解成多個2點DFT，所以N必須是2的冪次。

接下來，本文用一個$N=8$點FFT進行說明。

首先寫出8點DFT原始公式
$DFT[x(n)]=X(m)=\sum_{n=0}^{7}x(n)exp(-j\frac{2\pi mn}{N}),m\in[0,7]$
根據上文的FFT核心公式，8點DFT可以分爲兩個4點DFT，即
$X(m)=\sum_{n=0}^{3}x(2n)W_4^{mn}+W_8^m\sum_{n=0}^3x(2n+1)W_4^{mn},m\in[0,3]\\ X(m+4)=\sum_{n=0}^3x(2n)W_4^{mn}-W_8^m\sum_{n=0}^3(2n+1)W_4^{mn},m\in[0,3]\\$
令
$X_1(m)=\sum_{n=0}^{3}x(2n)W_4^{mn}\\ X_2(m)=\sum_{n=0}^{3}x(2n+1)W_4^{mn}\\$
則
$X(m)=X_1(m)+W_8^mX_2(m),m\in[0,3]\\ X(m+4)=X_1(m)-W_8^mX_2(m),m\in[0,3]\\$
1、對$X_1(m)$進行奇偶分解，有
$\begin{aligned} X_1(m)&=\sum_{n=0}^{1}x(2(2n))W_4^{m2n}+\sum_{n=0}^1x(2(2n+1))W_4^{m(2n+1)},m\in[0,3]\\ &=\sum_{n=0}^{1}x(4n)\exp(-\frac{2\pi m2n}{4})+W_4^m\sum_{n=0}^1x(4n+2)\exp(-j\frac{2\pi m2n}4),m\in[0,3]\\ &=\sum_{n=0}^{1}x(4n)W_2^{mn}+W_4^m\sum_{n=0}^1x(4n+2)W_2^{mn},m\in[0,1]\\ &=x(0)W_2^0+x(4)W_2^m+W_4^m[x(2)W_2^0+x(6)W_2^m],m\in[0,1]\\ &=x(0)+x(4)W_2^m+W_4^m[x(2)+x(6)W_2^m],m\in[0,1]\\ \end{aligned}$
即
$X_1(0)=x(0)+x(4)+x(2)+x(6)\\ X_1(1)=x(0)-x(4)+W_4^1[x(2)-x(6)]$
考慮$m+2\in[2,3]$的情況
$\begin{aligned} X_1(m+2)&=\sum_{n=0}^{1}x(4n)W_2^{(m+2)n}+W_4^{m+2}\sum_{n=0}^1x(4n+2)W_2^{(m+2)n},m\in[0,1]\\ &=\sum_{n=0}^{1}x(4n)W_2^{mn}\exp(-j2\pi n)+W_4^{m}\exp(-j\pi)\sum_{n=0}^1x(4n+2)W_2^{mn}\exp(-j2\pi n),m\in[0,1]\\ &=\sum_{n=0}^{1}x(4n)W_2^{mn}-W_4^{m}\sum_{n=0}^1x(4n+2)W_2^{(m+2)n},m\in[0,1]\\ &=x(0)W_2^0+x(4)W_2^m-W_4^m[x(2)W_2^0+x(6)W_2^m],m\in[0,1]\\ &=x(0)+x(4)W_2^m-W_4^m[x(2)+x(6)W_2^m],m\in[0,1]\\ \end{aligned}$
即
$X_1(2)=x(0)+x(4)-[x(2)+x(6)]\\ X_1(3)=x(0)-x(4)-W_4^1[x(2)-x(6)]$
2、對$X_2(m)$進行奇偶分解，有
$\begin{aligned} X_2(m)&=\sum_{n=0}^{1}x(2(2n)+1)W_4^{m2n}+\sum_{n=0}^1x(2(2n+1)+1)W_4^{m(2n+1)},m\in[0,3]\\ &=\sum_{n=0}^{1}x(4n+1)\exp(-\frac{2\pi m2n}{4})+W_4^m\sum_{n=0}^1x(4n+3)\exp(-j\frac{2\pi m2n}4),m\in[0,3]\\ &=\sum_{n=0}^{1}x(4n+1)W_2^{mn}+W_4^m\sum_{n=0}^1x(4n+3)W_2^{mn},m\in[0,1]\\ &=x(1)W_2^0+x(5)W_2^m+W_4^m[x(3)W_2^0+x(7)W_2^m],m\in[0,1]\\ &=x(1)+x(5)W_2^m+W_4^m[x(3)+x(7)W_2^m],m\in[0,1]\\ \end{aligned}$
即
$X_2(0)=x(1)+x(5)+x(3)+x(7)\\ X_2(1)=x(1)-x(5)+W_4^1[x(3)-x(7)]$
考慮$m+2\in[2,3]$的情況
$\begin{aligned} X_2(m+2)&=\sum_{n=0}^{1}x(4n+1)W_2^{(m+2)n}+W_4^{m+2}\sum_{n=0}^1x(4n+3)W_2^{(m+2)n},m\in[0,1]\\ &=\sum_{n=0}^{1}x(4n+1)W_2^{mn}\exp(-j2\pi n)+W_4^{m}\exp(-j\pi)\sum_{n=0}^1x(4n+3)W_2^{mn}\exp(-j2\pi n),m\in[0,1]\\ &=\sum_{n=0}^{1}x(4n+1)W_2^{mn}-W_4^{m}\sum_{n=0}^1x(4n+3)W_2^{mn},m\in[0,1]\\ &=x(1)W_2^0+x(5)W_2^m-W_4^m[x(3)W_2^0+x(7)W_2^m],m\in[0,1]\\ &=x(1)+x(5)W_2^m-W_4^m[x(3)+x(7)W_2^m],m\in[0,1]\\ \end{aligned}$
即
$X_2(2)=x(1)+x(5)-[x(3)+x(7)]\\ X_2(3)=x(1)-x(5)-W_4^1[x(3)-x(7)]$
將$X_1(m)$和$X_2(m)$分解後的表達式代回$X(m)$和$X(m+4)$中，有
$\begin{aligned} X(m)&=X_1(m)+W_8^mX_2(m)\\ X(m+2)&=X_1(m+2)+W_8^{m+2}X_2(m+2)\\ X(m+4)&=X_1(m)-W_8^mX_2(m)\\ X((m+2)+4)&=X_1(m+2)-W_8^{m+2}X_2(m+2)\\ m&\in[0,1] \end{aligned}$
將$m\in[0,1]$代入，得到$X(0),X(1),...,X(7)$的表達式
$\begin{aligned} X(0)&=x(0)+x(4)+x(2)+x(6)+x(1)+x(5)+x(3)+x(7)\\ X(1)&=x(0)-x(4)+W_4^1[x(2)-x(6)]+W_8^1\{x(1)-x(5)+W_4^1[x(3)-x(7)]\}\\ X(2)&=x(0)+x(4)-[x(2)+x(6)]+W_8^2\{x(1)+x(5)-[x(3)+x(7)]\}\\ X(3)&=x(0)-x(4)-W_4^1[x(2)-x(6)]+W_8^3\{x(1)-x(5)-W_4^1[x(3)-x(7)]\}\\ X(4)&=x(0)+x(4)+x(2)+x(6)-[x(1)+x(5)+x(3)+x(7)]\\ X(5)&=x(0)-x(4)+W_4^1[x(2)-x(6)]-W_8^1\{x(1)-x(5)+W_4^1[x(3)-x(7)]\}\\ X(6)&=x(0)+x(4)-[x(2)+x(6)]-W_8^2\{x(1)+x(5)-[x(3)+x(7)]\}\\ X(7)&=x(0)-x(4)-W_4^1[x(2)-x(6)]-W_8^3\{x(1)-x(5)-W_4^1[x(3)-x(7)]\}\\ \end{aligned}$
即
$\begin{aligned} X(0)&=x(0)+x(4)+x(2)+x(6)+x(1)+x(5)+x(3)+x(7)\\ X(1)&=x(0)-x(4)-j[x(2)-x(6)]+\exp(-j\frac{\pi}4)\{x(1)-x(5)-j[x(3)-x(7)]\}\\ X(2)&=x(0)+x(4)-[x(2)+x(6)]-j\{x(1)+x(5)-[x(3)+x(7)]\}\\ X(3)&=x(0)-x(4)+j[x(2)-x(6)]+\exp(-j\frac{3\pi}4)\{x(1)-x(5)+j[x(3)-x(7)]\}\\ X(4)&=x(0)+x(4)+x(2)+x(6)-[x(1)+x(5)+x(3)+x(7)]\\ X(5)&=x(0)-x(4)-j[x(2)-x(6)]-\exp(-j\frac{\pi}4)\{x(1)-x(5)-j[x(3)-x(7)]\}\\ X(6)&=x(0)+x(4)-[x(2)+x(6)]+j\{x(1)+x(5)-[x(3)+x(7)]\}\\ X(7)&=x(0)-x(4)+j[x(2)-x(6)]-\exp(-j\frac{3\pi}4)\{x(1)-x(5)+j[x(3)-x(7)]\}\\ \end{aligned}$