c语言fft

FFT是快速傅里叶变换（Fast Fourier Transform）的缩写，是一种高效的计算离散傅里叶变换（DFT）的算法。在C语言中，可以使用库函数来实现FFT，如fftw库、Intel MKL库等。

以下是一个使用fftw库实现FFT的例子：

#include <fftw3.h>

void fft(double *input, double *output, int n) {
    fftw_complex *in, *out;
    fftw_plan plan;

    in = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * n);
    out = (fftw_complex*) fftw_malloc(sizeof(fftw_complex) * n);
    plan = fftw_plan_dft_r2c_1d(n, input, out, FFTW_ESTIMATE);

    for (int i = 0; i < n; i++) {
        in[i][0] = input[i];
        in[i][1] = 0;
    }

    fftw_execute(plan);

    for (int i = 0; i < n / 2 + 1; i++) {
        output[i] = sqrt(out[i][0] * out[i][0] + out[i][1] * out[i][1]);
    }

    fftw_destroy_plan(plan);
    fftw_free(in);
    fftw_free(out);
}

该函数接受一个输入数组`input`，一个输出数组`output`和一个数组长度`n`，并计算出输入数组的FFT结果，存储在输出数组中。这里使用了fftw3库的`fftw_plan_dft_r2c_1d`函数来创建一个FFT计划，并使用`fftw_execute`函数执行计划。

值得注意的是，由于FFT计算结果包含实部和虚部，因此在输出结果中，需要将它们的模值计算出来并存储在输出数组中。