C++ fftw3 实现fft

您好，关于如何在C++中使用fftw3库实现FFT，我可以提供以下示例代码和步骤：

1. 安装fftw3库。您可以从官网（https://www.fftw.org/download.html）下载适合您平台的版本，并按照官方文档进行安装。

2. 在C++代码中包含fftw3头文件，并链接fftw3库。

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

3. 定义输入和输出数组，并创建fftw3计算计划。

int n = 1024; // FFT长度
double* in = (double*)fftw_malloc(sizeof(double) * n); // 输入数组
fftw_complex* out = (fftw_complex*)fftw_malloc(sizeof(fftw_complex) * n); // 输出数组
fftw_plan plan = fftw_plan_dft_r2c_1d(n, in, out, FFTW_ESTIMATE); // 创建计算计划

在上面的代码中，我们定义了一个长度为n的输入数组in和输出数组out，并创建了一个fftw3计算计划plan。fftw_plan_dft_r2c_1d表示我们要计算一个实数序列的FFT，FFTW_ESTIMATE表示我们使用快速计算方法。

4. 填充输入数组，并执行计算计划。

for (int i = 0; i < n; i++) {
    in[i] = sin(2 * M_PI * i / n);
}

fftw_execute(plan);

在上面的代码中，我们填充了输入数组，然后使用fftw_execute函数执行计算计划。

5. 获取计算结果。

for (int i = 0; i < n; i++) {
    printf("%f + %fi\n", out[i][0], out[i][1]);
}

在上面的代码中，我们打印了输出数组的全部元素，每个元素是一个复数，由实部和虚部组成。

以上是一个基本的fftw3库的使用示例，您可以根据您的需求修改输入数组和计算计划的选项，以及处理输出数组的方式。希望能对您有所帮助！