VS2017配置FFTW库及FFTW库lib生成

是世界范围内最有名的计算fft的库，平时熟悉的Origin及Matlab计算fft采用的也是fftw。fftw库支持计算一维及多维的任意长度数据，反正用就行了，毕竟人家成熟商业软件计算fft用的都是这个库，更多详情请查看其官网fftw库官网。

FFTW官网提供了源码及头文件、dll、def文件，一般我们需要根据def文件生成lib文件，也可以用源码+Cmake生成dll和lib文件。网上得到lib文件的方式基本都是根据vs自带的lib.exe将.def文件生成.lib文件。如下图所示（利用vs2017生成64位的fftw lib文件），将ide文件中的4个文件msobj140.dll、mspdb140.dll、mspdbcore.dll、mspdbsrv.exe复制到x64文件夹下，并且还要将三个以def为后缀名的fftw文件拷到该x64文件夹下，最后利用cmd进到x64文件下执行lib.exe，然后分别执行：
lib /machine:x64 /def:libfftw3-3.def
lib /machine:x64 /def:libfftw3f-3.def
lib /machine:x64 /def:libfftw3l-3.def
以生成lib文件，最后还要在VS里配置一下。


用这个方法，出现如下错误，提示不能加载mspdcore.dll文件，最后发现可能是版本不一致的问题。于是，后来便在网上找相同版本的dll文件，没找到，放弃。为了得到lib文件，后来又尝试了其他几种方法，例如利用32位的def生成lib+采用64位的dll文件、卸载vs重新安装vs2017，但都无果。今晚，和师兄聊天，得知他那有利用Cmake生成的lib文件，于是便拷了一份，最终结束lib的生成之路，并能成功运行。

最后，给出福利链接：fftw lib文件（提取码：7385）
，里面有32位、64位的bin、lib、include文件，可以直接下载用，不用再自己生成lib文件。

VS配置：
（1）将文件夹fftw_X64和dll文件拷贝到项目文件夹下
（2）包含目录

./fftw_X64/include

（3）库目录

./fftw_X64/lib

（4）附加依赖项

libfftw3-3.lib
libfftw3f-3.lib

测试代码：

#include <iostream>
#include "fftw3.h"

using namespace std;

#define REAL 0
#define IMAG 1

int main()
{
	const double PI = 3.14159265358979323846;

	const int L = 10000;
	double Fs = 500;
	double T = 1.0 / Fs;

	fftw_complex x[L];
	fftw_complex y[L];

	double time = 0.0;

	for (int i = 0; i < L; i++) {
		time = i * T;		
		x[i][REAL] = 0.7*sin(2 * PI * 50 * time) + sin(2 * PI * 120 * time);
		x[i][IMAG] = 0;		
	}
	fftw_plan plan = fftw_plan_dft_1d(L, x, y, FFTW_FORWARD, FFTW_ESTIMATE);
	fftw_execute(plan);

	double *freq = new double[L];
	double *Amp = new double[L];

	cout<<"Freq"<< "  " << "Amp"<< endl;
	int i = 0;
	freq[0] = 0;
	Amp[0] = sqrt(y[i][REAL] * y[i][REAL] + y[i][IMAG] * y[i][IMAG]) / L;
	cout << freq[0] << "  " << Amp[0] << endl;
	for (int i = 1; i < L/2; i++) {
		freq[i] = i*Fs/L;
		Amp[i] = sqrt(y[i][REAL] * y[i][REAL] + y[i][IMAG] * y[i][IMAG])*2 / L;
		cout << freq[i] << "  " << Amp[i] << endl;
	}
	freq[i] = Fs/2;
	Amp[i] = sqrt(y[L/2][REAL] * y[L/2][REAL] + y[L/2][IMAG] * y[L/2][IMAG]) / L;
	cout << freq[i] << "  " << Amp[i] << endl;


	fftw_destroy_plan(plan);
	delete []freq;
	delete []Amp;

}

测试结果：