FFT海浪模拟——理论部分（二）

由于使用标准傅里叶变换计算海面计算量过大，所以使用离散形式的傅里叶变换（DFT）代替

1.离散傅里叶变换

1）定义式

$$X[k]=\sum _{n=0}^{N-1}x[n]e^{-j\frac{2\pi k}{N}n},0\le k\le N-1$$

为简化表达，令
$$W_N=e^{-j\frac{2\pi }{N}}$$
则定义式可化为
$$X[k]=\sum _{n=0}^{N-1}x[n]W_N^{kn}$$

2）性质

i.对于带有e的指数项W，有如下性质：
$$W_N^{k(2n)}=e^{-j\frac{2\pi k}{N}2n}=e^{-j\frac{2\pi k}{\frac{N}{2}}n}=W_{\frac{N}{2}}^{kn}(代换)$$

$$W_N^{kn}=W_N^{k(n+N)}=W_N^{(k+N)n}(周期性)$$

3）计算方式

采用了分治算法的思想，将一个长度为N的信号分为两个长度为N/2的信号进行计算，一个为偶数信号，另一个为奇数信号，表示如下：

首先将x[n]分解为奇数项和偶数项

偶数项：x[2n]

奇数项：x[2n+1]

此处的n取值范围为
n ∈ { 0 , 1 , . . . , N 2 − 1 } n \in \{0,1,...,\frac{N}{2}-1 \} n∈{0,1,...,2N​−1}

则可将DFT定义式分解为
$$\begin{gathered} X[k]=\sum _{n=0}^{\frac{N}{2}-1}x[2n]W_N^{k(2n)}+\sum _{n=0}^{\frac{N}{2}-1}x[2n+1]W_N^{k(2n+1)} \\ =\sum _{n=0}^{\frac{N}{2}-1}x[2n]W_N^{k(2n)}+W_N^k\sum _{n=0}^{\frac{N}{2}-1}x[2n+1]W_N^{k(2n)} \\ \sum _{n=0}^{\frac{N}{2}-1}x[2n]W_{\frac{N}{2}}^{kn}+W_N^k\sum _{n=0}^{\frac{N}{2}-1}x[2n+1]W_{\frac{N}{2}}^{kn} \end{gathered}$$
再令偶数项和奇数项分别为
$$\begin{gathered} Even:E[k]=\sum _{n=0}^{\frac{N}{2}-1}x[2n]W_{\frac{N}{2}}^{kn} \\ Odd:O[k]=\sum _{n=0}^{\frac{N}{2}-1}x[2n+1]W_{\frac{N}{2}}^{kn} \end{gathered}$$
根据W式的周期性，有
$$\begin{gathered} E[k]=E[k+\frac{N}{2}] \\ O[k]=O[k+\frac{N}{2}] \end{gathered}$$
这样我们就可以将定义式表示成
$$X[k]=E[k]+W_N^{k}O[k],0\le k\le N-1$$

2.一种计算示例（蝶形图表示）

取N=8，则计算过程如下（从最终结果反推初始输入）：

第1次分解（4上4下）：

第二次（2上2下）：

第三次（1上1下）：

以上公式推导过程来自

1.fft海面模拟（二）
2.如何快速理解傅里叶变换的计算和概念？