前言

本文介绍如何使用FFTW来进行快速傅里叶。
我是用的是vs2019。

1. 下载FFTW，并进行编译

下载地址 http://fftw.org/
下载完成后进行解压。

生成lib和dll文件

使用vs2019自带的工具，打开图中所示工具，建议使用管理员模式打开。

然后 cd 到刚刚解压的fftw文件夹下。

然后分别键入
32位

lib /machine:x86 /def:libfftw3-3.def
lib /machine:x86 /def:libfftw3l-3.def
lib /machine:x86 /def:libfftw3f-3.def

64位

lib /machine:x64 /def:libfftw3-3.def
lib /machine:x64 /def:libfftw3l-3.def
lib /machine:x64 /def:libfftw3f-3.def

就会生成三个对应的.lib和.exp文件文件
如下图所示：在这里插入图片描述
2. 将上图中三个标记处的文件移动到C:\Windows\SysWOW64中（重要）

2. 创建c++工程使用fftw

创建空白c++windows控制台工程。

略去过程

在工程中，创建main.cpp

下面为测试代码：

#include <iostream>

#include "fftw3.h"
#pragma comment(lib, "libfftw3-3.lib")

//实部与虚部
#define REAL 0
#define IMAG 1

using namespace std;

int main() {
	/*
	*fftw_complex 是FFTW自定义的复数类
	*引入<complex>则会使用STL的复数类
	*/
	fftw_complex x[5];
	fftw_complex y[5];

	for (int i = 0; i < 5; i++) {
		x[i][REAL] = i;
		x[i][IMAG] = 0;
	}

	//定义plan，包含序列长度、输入序列、输出序列、变换方向、变换模式
	fftw_plan plan = fftw_plan_dft_1d(5, x, y, FFTW_FORWARD, FFTW_ESTIMATE);

	//对于每个plan，应当"一次定义 多次使用"，同一plan的运算速度极快
	fftw_execute(plan);

	for (int i = 0; i < 5; i++) {
		cout << y[i][REAL] << "  " << y[i][IMAG] << endl;
	}

	//销毁plan
	fftw_destroy_plan(plan);

	cout << "\nPress Enter to exit..." << endl;
	cin.get();

	return 0;
}