C++实现快速傅里叶变换

基础公式

傅里叶变换（FT）：
$$F(\omega )=F[f(t)]={\int }_{-\infty }^{+\infty }f(t)e^{-j\omega t}dt$$
离散傅里叶变换（DFT）：
$$X(k)=\sum _0^{N-1}x(n)W_N^{kn}(k=0,1,2,3\cdots N-1)$$
其中$W_N^{kn}=e^{-j\frac{2\pi }{N}kn}$，称为旋转因子。

\quadFFT是基于DFT的一种算法，将长序列的DFT分解为短序列的DFT，利用旋转因子的周期性、对称性、可约性，减少了重复运算。

FFT原理

对于一个多项式
$$A(x)=\sum _0^{n-1}{a}_i{x}^i=a_0+a_{1}x+a_2{x}^2+\cdots +a_{n-1}{x}^{n-1}$$
按照下标的奇偶性把A(x)分成两部分
A(x)=(a0+a2x2+a4x4+⋯+an−2xn−2)+(a1x+a3x3+a5x5+an−1xn−1)  =(a0+a2x2+a4x4+⋯+an−2xn−2)+x(a1+a3x2+a5x4+an−1xn−2)A(x)=(a_0+a_2x^2+a_4x^4+\cdots+a_{n-2}x^{n-2})+(a_1x+a_3x^3+a_5x^5+a_{n-1}x^{n-1}) \\ \quad\quad\; =(a_0+a_2x^2+a_4x^4+\cdots+a_{n-2}x^{n-2})+x(a_1+a_3x^2+a_5x^4+a_{n-1}x^{n-2}) A(x)=(a0​+a2​x2+a4​x4+⋯+an−2​x