linux 和 Windows下FFTW庫的安裝
main: test_fftw.o
g++ -o test_fftw test_fftw.cpp $(INCLUDESH) $(INCLUDESLIB) $(LIBS)
test_fftw.o:
g++ -c test_fftw.cpp $(INCLUDESH)
clean:
rm -f *.o
.PHONY:clean
也是可以成功的。
最後附上文章:
FFTW(the Fastest Fourier Transform in the West)庫是由MIT(Massachusetts Institute of Technology)的Matteo Frigo和Steven G. Johnson開發的,用於一維和多維實數或複數的離散傅立葉變換。
1. 下載 fftw-2_1_3_tar.gz (
2. tar zxvf fftw-2_1_3_tar.gz 展開壓縮檔案
3. 在Linux中安裝FFTW:
a. ./configure --enable-type-prefix --prefix=/usr/local/fftw --with-gcc --disable-fortran --enable-i386-hacks
其中,--enable-type-prefix 引數是為了同時使用single precision(單精度)和double precision(雙精度),如果不使用它,最後只有以rfftw開頭的檔案被安裝(real fftw);
--prefix= 引數是設定安裝目錄;
--with-gcc 使用gcc編譯器;
--disable-fortran 引數為了不包含Fortran呼叫的機制;
--enable-i386-hacks 為Pentium和x86以後的CPU優化gcc的編譯速度。
b.make 編譯
c.make install 安裝,這一次安裝完後,在安裝目錄中存在以dfftw和drfftw開標頭檔案,但沒有sfftw開頭的檔案
d.make clean 還需要安裝一次,先清除
e../configure --enable-float --enable-type-prefix --prefix=/usr/local/fftw --with-gcc --disable-fortran --enable-i386-hacks
其中,--enable-float 為了生成單精度計算的標頭檔案和庫檔案,即以sfftw開頭的檔案。
f.make 重新編譯
g.make install 再一次安裝,安裝完後,目錄中便會同時存在sfftw和dfftw開頭的檔案(用於複數函式/complex function的FFT變換)和srfftw與drfftw開頭的檔案(用於實數函式的FFT變換)
例如,如果需要用到雙精度的實數FFT變換/FFTs,那麼在編譯的連結命令中需要按如下順序加入
-ldrfftw -ldfftw引數
下面的是如何使用的一個例子
#include <complex> #include <fftw3.h> #include <math.h> #include <iostream> #define N 10 using namespace std; int main(int argc, char * argv[]){ fftw_complex in[N], out[N]; fftw_plan p; p=fftw_plan_dft_1d(N,in,out,FFTW_FORWARD,FFTW_MEASURE); for(int i=0;i <N;i ++) { in[i][0]=i; in[i][1]=0.0; } fftw_execute(p); for(int i=0;i <N;i ++){ cout<<out[i][0]<<" "<<out[i][1]<<endl; } complex<double> temp = 0.0; for(int k =0; k < N; k ++){ double pi = 4*atan(1.0); temp += exp(complex<double>(0.0,-2.0*pi*3*k/N))*complex<double>(in[k][0],in[k][1]); } cout<<"out[3] is "<<temp<<endl; fftw_complex out1[N]; fftw_plan p1; p1=fftw_plan_dft_1d(N,out1,in,FFTW_BACKWARD,FFTW_MEASURE); for(int i=0;i <N;i ++){ out1[i][0]=out[i][0]; out1[i][1]=out[i][1]; } fftw_execute(p1); for(int i=0;i <N;i ++){ cout<<in[i][0]<<" "<<in[i][1]<<endl; } fftw_destroy_plan(p); fftw_destroy_plan(p1); return 1; }結果是: