LIBS = -L/usr/local/fftw/lib

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              （www.fftw.org, or www.rpmfind.net ）
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;
}