分裂基 c语言算法,分裂基快速傅里叶变换

一、功能

计算复序列的分裂基快速傅里叶变换。

二、方法简介

序列\(x(n)(n=0,1,...,N-1)\)的离散傅里叶变换定义为

\[

X(k)=\sum_{n=0}^{N-1}x(n)W_{N}^{nk}, \qquad k=0,1,...,N-1

\]

其中\(W_{N}^{nk}=e^{-j\frac{2\pi nk}{N}}\)，将\(X(k)\)按序号\(k\)的奇偶分成两组。当\(k\)为偶数时，进行基2频率抽取分解， \(X(k)\)可表示为

\[

X(2k)=\sum_{n=0}^{N/2-1}[x(n)+x(n+\frac{N}{2})]W_{N}^{2nk} \ , \ k=0,1,...,\frac{N}{2}-1

\]

当\(k\)为奇数时进行基4 频率抽取分解，$ X(k)$可表示为

\[

\left\{\begin{matrix}X(4k+1)=\sum_{n=0}^{N/4-1}{[x(n)-x(n+\frac{N}{2})]-j[x(n+\frac{N}{4})-x(n+\frac{3N}{4})]}W_{N}^{n}W_{N}^{4nk}\\ X(4k+3)=\sum_{n=0}^{N/4-1}{[x(n)-x(n+\frac{N}{2})]+j[x(n+\frac{N}{4})-x(n+\frac{3N}{4})]}W_{N}^{n}W_{N}^{4nk}\end{matrix}\right.\\k = 0,1,...,\frac{N}{4}-1

\]

由此可见，一个\(N\)点的DFT 可以分解为一个\(N/2\)点的DFT 和两个\