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]

//正的一次谐波系数保存在data[3],data[4]

//负的一次谐波系数保存在data[2*nn-1],data[2*nn]

//最高次谐波系数保存在data[nn+1],data[nn+2]

//nn必须为2的幂-fft