数字信号处理——第四章、快速傅里叶变换（FFT）

4.2 基-2FFT算法

4.2.1 直接计算DFT的特点及减少运算量的基本途径

有限长序列x(n)的N点DFT为：
$$X(k)=\sum _{-\infty }^{\infty }x(n)W_N^{kn}$$
计算一个k值需要N次乘法，N个k就需要N*N次乘法，计算次数巨大。所以我们可以将N点DFT分解为几个较短的DFT。另外，旋转因子$W_N^m$具有明显的周期性：
$$\begin{gathered} 周期性：W_N^{m+lN}=W_N^m \\ 对称性：W_N^{-m}=W_N^{N-M}或者[W_N^{N-m}]=W_N^m \end{gathered}$$
FFT算法就是不断地把长序列的DFT分解成几个短序列的DFT，并利用$W_N^{kn}$的周期性和对称性来煎炒DFT的运算次数

4.2.2 时域抽取发基-2FFT基本原理

基-2FFT算法分为两类：时域抽取法（DIIT-FFT）和频域抽取法（DIF-FFT）

序列x(n)可以按n的奇偶分解成N/2的子序列
$$\begin{gathered} x_1(r)=x(2r)r=0,1,\dots \dots ,N/2-1 \\ {x}_1(r)=x(2r+1)r=0,1,\dots \dots ,N/2-1 \end{gathered}$$
则x(n)的DFT为：
X ( k ) = ∑ r = 0 N / 2 − 1 x ( 2 r ) W N 2 k r +