FFTW编译以及使用

FFTW官网www.fftw.org上可以下载fftw-3.3.4.tar.gz

解压tar -xzvf fftw-3.3.4.tar.gz

1 ./configure

     

2 make

3 make install

以编译float版本为例：

    

    ./configure --prefix=/home/tuzb/fftw --enable-shared --enable-float --disable-fortran

 make

 make install

1.FFTW有三个版本的数据类型float , double , long double。

2.都使用同样的头文件fftw3.h

3.所有以小写fftw_开头的函数替换为fftwf_或者fftwl_

4.函数参数中double替换成float或者long double

使用实例

#include <stdlib.h>

#include <stdio.h>

#include <math.h>

#include "fftw3.h"

int main(int argc,char** argv){

    printf("--------------start---------------\n");

    int len = 8;

    float *in = NULL;

    // 如果要使用float版本,需先引用float版本的lib库,然后在fftw后面加上f后缀即可.

    fftwf_complex *out = NULL; // fftwf_complex --> 即为float版本

    fftwf_plan p;

    in = (float *) fftwf_malloc(sizeof(float) * len);

    out = (fftwf_complex *) fftwf_malloc(sizeof(fftw_complex) * len);

    float dx = 1.0 / len;

    // 输入纯实数

    for (int i = 0; i < len; i++) {

        in[i] = sinf(2 * M_PI * dx * i) + sinf(4 * M_PI * dx * i);

        printf("%.2f ", in[i]);

    }

    printf("\n\n");

    // 傅里叶变换

    p = fftwf_plan_dft_r2c_1d(len, in, out, FFTW_ESTIMATE);

    fftwf_execute(p);

    // 输出幅度谱

    for (int i = 0; i < len; i++) {

        float len = sqrt(out[i][0] * out[i][0] + out[i][1] * out[i][1]);

        printf("%.2f ", len);

    }

    printf("\n");

    // 释放资源

    fftwf_destroy_plan(p);

    fftwf_free(in);

    fftwf_free(out);

    printf("-------------- end ---------------\n");

    return 0;

}