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;
}