fft#include <math.h> #define SWAP(a,b) tempr=(a) (a)=(b) (b)=tempr //当isign=1时,将输入数据data[1..2*nn]用它的离散Fourier变换系数代替 //当isign=-1时,将输入数据替换为其的离散IFFT系数乘以nn //data是一个保存nn个复数的数组 //Fourier系数的保存方法为F0的实部和虚部保存在data[1],data[2] //正的一次谐波系数保存在dat

快速傅里叶变换,一种快速算法，C语言编写-This is fast Fourier transform.The standard deviation of the FFT method is about 1/6 of that of the log-fit method and about a half of the standard deviation of the Prony method.