fftw的使用

fftw 号称是最快的fft函数库

fftw_complex in[1024], out[1024];
    fftw_plan p;


    p=fftw_plan_dft_1d(1024,in,out,FFTW_FORWARD,FFTW_MEASURE);
    for(int i=0;i <1024;i ++) {
        in[i][0]=data[i];
        in[i][1]=0.0;
    }
   

    fftw_execute(p);

	for(int i = 1;i<=N/2;i++)
    {
        real = out[i][0];
		imag = out[i][1];
		result[i-1]=sqrt(XTRACT_SQ(real) + XTRACT_SQ(imag))/N;
		result[N/2+i-1]=i*q;
		max=result[i-1]>max ? result[i-1]: max;
    }
	fftw_destroy_plan(p)