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