C语言实现FFT并与DFT比较（徐士良老师编写的c语言算法程序）

一 . FFT方法说明
计算n个采样点

P=p ₀, p ₁, ..........., p _n-1

的傅里叶变换，可以归结为计算多项式 :

f(x)= p ₀+p ₁ x+ p ₂ x ² +...... +p _n-1 x ^n-1

在各n次单位根1，w，w ²，…,w ^n-1上的值，即

f ₀=p ₀+p ₁+..........+p _n-1 f ₁=p ₀+p ₁w+...........+p _n-1w ^n-1 f ₂=p ₀+p ₁w ²++p ₂(w ²) ² +..........p _n-1(w ²) ^n-1 .......... . f _n-1=p ₀+p ₁w ^n-1++p ₂(w ^n-1) ² +..........p _n-1(w ^n-1) ^n-1

其中, w= e − j 2 π N {\rm e}^{-j\frac{2π}{N