FFTW3库的应用。_fftw 复数运算

号称最快的FFT算法。是个C语言程序集。

下面是具体应用方法和注意事项：

以复数傅里叶变换为例了。我用到的就是复数的计算。

1.定义输入输出变量

fftw_complex * fttInput; //输入需要变换的复数数据
fftw_complex * fttOutput;//输出变换后的复数结果

2.定义变换方案

fftw_plan  p;//创建变换方案。例如我做的是对图像的傅里叶变换 就选择p = fftw_plan_dft_2d。对一维数据做变换就选择fftw_plan_dft_1d.

3.实施变换

fftw_execute(p1);//实施变换

4.注意事项：一个fftw_plan只能对固定的in/out进行变换， 但可以在变换后改变in的内容（大小不变）以用同一个方案执行新的变换。这个尤其重要。在我的试验中发现第一次对某个图像做傅里叶变换时，性能并没有提升反而慢于opencv，然而对同样大小图像内容不同做多次变换（不重新创建plan）时，性能提升明显。所以提速关键应该在此吧。还有FFTW_FORWARD是正变换，FFTW_BACKWARD是逆变换。

上完整代码如下：

fftw_complex * fttinput;
fftw_complex * fttoutput;
fftw_plan  p;

fttinput = (fftw_complex *)fftw_malloc(sizeof(fftw_complex) * w* h);

fttoutput = (fftw_complex *)fftw_malloc(sizeof(fftw_complex) * w* h);

p = fftw_plan_dft_2d(h, w,fttinput, fttoutput,FFTW_FORWARD, FFTW_ESTIMATE);

for（int i = 0； i < 100; i ++ ）

{

     //对某组图像做傅里叶变换 不断改变fttinput的内容 用fftw_execute执行傅里叶变换。不断刷新结果

//结果不断刷新在fttoutput 中。

     fftw_execute(p);

}

//这样就完成了对100张图像的复数傅里叶变换了。