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;
}
————————————————
版权声明:本文为CSDN博主「xiaotu2011」的原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/tuzanbo/article/details/75309647

 

Linux下FFTW库的安装

FFTW(the Fastest Fourier Transform in the West)库是由MIT(Massachusetts Institute of Technology)的Matteo Frigo和Steven G. Johnson开发的,用于一维和多维实数或复数的离散傅里叶变换。

1. 下载 fftw-3_2_2_tar.gz              (www.fftw.org, or www.rpmfind.net )


2. tar zxvf fftw-3_2_2_tar.gz            展开压缩文件


3. 在Linux中安装FFTW:
总共需要两次安装  
 第一次安装就按下面的 a 配置就行
第二次安装时为了 生成libfftwf.so的库文件   以f结尾的代表可以的数据类型是 float  加上 --enable-float   第一次安装的时候不能加这个选项  否则会出错

a.    
./configure --enable-type-prefix --prefix=/usr/local/fftw --with-gcc --disable-fortran --enable-i386-hacks  --enable-shared=yes
其中,
--enable-type-prefix    参数是为了同时使用single precision(单精度)和double precision(双精度),如果不使用它,最后只有以rfftw开头的文件被安装(real fftw);

--prefix=                      参数是设定安装目录;

--with-gcc                    使用gcc编译器;

--disable-fortran         参数为了不包含Fortran调用的机制;

--enable-i386-hacks   为Pentium和x86以后的CPU优化gcc的编译速度。
 --enable-shared=yes  为了生成动态库.so文件  否则生成的都是.a 文件  没法编译 (我就卡在这里)

第一次安装
1.
make                          编译
2.
make install               安装,这一次安装完后,在安装目录中存在以dfftw和drfftw开头文件,但没有sfftw开头的文件
3.
make clean              还需要安装一次,先清除

第二次安装
1.
./configure --enable-float --enable-type-prefix --prefix=/usr/local/fftw --with-gcc --disable-fortran --enable-i386-hacks --enable-shared=yes
其中,--enable-float      为了生成单精度计算的头文件和库文件,即以sfftw开头的文件。
2.
make                      重新编译
3.
make install            再一次安装,安装完后,目录中便会同时存在sfftw和dfftw开头的文件(用于复数函数/complex function的FFT变换)和srfftw与drfftw开头的文件(用于实数函数的FFT变换)
————————————————
版权声明:本文为CSDN博主「IT张先生」的原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接及本声明。
原文链接:https://blog.csdn.net/zwjzwj108108/article/details/24254067

posted @ 2020-03-16 10:31  苍月代表我  阅读(7753)  评论(0编辑  收藏  举报