C FFT 快速傅里叶库- fftw-3.3.8
Posted 姚家湾
tags:
篇首语:本文由小常识网(cha138.com)小编为大家整理,主要介绍了C FFT 快速傅里叶库- fftw-3.3.8相关的知识,希望对你有一定的参考价值。
在编写C++ 数据采集软件时,会使用到快速傅里叶变换算法,例如振动分析。网络上介绍的比较多的是fftw3。
FFTW—Fastest Fourier Transform in the West,是由 MIT 的 Matteo Frigo 博士和 Steven G. Johnson 博士开发的一个完全免费的软件包。FFTW 最初的 release 版本于 1997 年发布,最新的 release 版本 fftw-3.3.8。正如它的名字声称的,号称是最快的FFT软件包。
fftw支持任意大小的、任意维数的数据的离散傅里叶变换(DFT),并且还支持离散余弦变换(DCT)、离散正弦变换(DST)和离散哈特莱变换(DHT)
安装
从官网(http://fftw.org/index.html)下载源码压缩包 fftw-3.3.8.tar.gz,然后解压到linuxwork 目录中。按照http://www.fftw.org/fftw3_doc/Installation-on-Unix.html的指导安装:
cd fftw-3.3.8
/configure CC=clang
make
sudo make install
./configure的格式为 ./configure CC="<the name of your C compiler>",我使用的是clang编译器。
测试程序
#include<fftw3.h>
int main(void)
int N;
N=1024;
fftw_complex *in, *out;
fftw_plan my_plan;
in = (fftw_complex*) fftw_malloc(sizeof(fftw_complex)*N);
out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex)*N);
my_plan = fftw_plan_dft_1d(N, in, out, FFTW_FORWARD, FFTW_ESTIMATE);
fftw_execute(my_plan);
fftw_destroy_plan(my_plan);
fftw_free(in);
fftw_free(out);
return 0;
编译
clang test.cpp -o test -lfftw3 –lm
-lm 是指使用math 库。
例2
下面的程序引用了 https://blog.csdn.net/cyh706510441/article/details/46676123 中的代码。
#include <stdio.h>
#include <stdlib.h>
#include <math.h>
#include "fftw3.h"
#define PI 3.1415926
int main()
int len = 128;
double *in = NULL;
// 如果要使用float版本,需先引用float版本的lib库,然后在fftw后面加上f后缀即可.
fftw_complex *out = NULL;// fftwf_complex --> 即为float版本
fftw_plan p;
in = (double *)fftw_malloc(sizeof(double) * len);
out = (fftw_complex *)fftw_malloc(sizeof(fftw_complex) * len);
double dx = 1.0 / len;
// 输入纯实数
for (int i = 0; i < len; i++)
in[i] = sin(2*PI * dx*i) + sin(4*PI * dx*i);
printf("%.2f ", in[i]);
printf("\\n\\n");
// 傅里叶变换
p = fftw_plan_dft_r2c_1d(len, in, out, FFTW_ESTIMATE);
fftw_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");
// 释放资源
fftw_destroy_plan(p);
fftw_free(in);
fftw_free(out);
return 0;
编译
clang test2.c -o test2 -lfftw3 –lm
结束语
好了,搞定!
以上是关于C FFT 快速傅里叶库- fftw-3.3.8的主要内容,如果未能解决你的问题,请参考以下文章