intN=;// FFT窗口大小 fftwf_complex*in=(fftwf_complex*)fftwf_malloc(sizeof(fftwf_complex)*N); fftwf_complex*out=(fftwf_complex*)fftwf_malloc(sizeof(fftwf_complex)*N); fftwf_plan plan=fftwf_plan_dft_r2c_1d(N,(float*)in,out,FFTW_ESTIMATE); // 进行频域分析 for(inti=0;i<N;...
fftw_plan plan; intn=numSamples; double*in=newdouble[n]; fftw_complex*out=reinterpret_cast<fftw_complex*>(fftw_malloc(sizeof(fftw_complex)*(n/2+1))); plan=fftw_plan_dft_r2c_1d(n,in,out,FFTW_ESTIMATE); // 填充输入数组(将音频数据转换为double类型) for(inti=0;i<n;++i){ in[i]...
I try the FFTW c2r to go "backwards" (full source code at the end) plan_forward = fftw_plan_dft_r2c_1d(N, in, out, FFTW_ESTIMATE); plan_backward = fftw_plan_dft_c2r_1d(N, out, in, FFTW_ESTIMATE | FFTW_PRESERVE_INPUT); fftw_execute(plan_forward); printf("\tDONE! (1= %f...