生成一份C语言的FFT代码

生成C语言的FFT代码需要用到FFT算法的具体实现方法，以下是一份基于Cooley-Tukey算法的C语言FFT代码示例：

#include <stdio.h>
#include <math.h>

// Cooley-Tukey FFT algorithm
void fft(double* x, double* y, int n, int inv)
{
    if (n == 1)
        return;

    // Divide
    double* xe = new double[n / 2];
    double* ye = new double[n / 2];
    double* xo = new double[n / 2];
    double* yo = new double[n / 2];
    for (int i = 0; i < n / 2; i++) {
        xe[i] = x[2 * i];
        ye[i] = y[2 * i];
        xo[i] = x[2 * i + 1];
        yo[i] = y[2 * i + 1];
    }

    // Conquer
    fft(xe, ye, n / 2, inv);
    fft(xo, yo, n / 2, inv);

    // Combine
    for (int k = 0; k < n / 2; k++) {
        double t_re = xo[k] * cos(-2 * M_PI * k / n) - yo[k] * sin(-2 * M_PI * k / n);
        double t_im = xo[k] * sin(-2 * M_PI * k / n) + yo[k] * cos(-2 * M_PI * k / n);
        x[k] = xe[k] + t_re;
        y[k] = ye[k] + t_im;
        x[k + n / 2] = xe[k] - t_re;
        y[k + n / 2] = ye[k] - t_im;
    }

    // Inverse FFT
    if (inv == 1) {
        for (int i = 0; i < n; i++) {
            x[i] /= n;
            y[i] /= n;
        }
    }

    delete[] xe;
    delete[] ye;
    delete[] xo;
    delete[] yo;
}

int main()
{
    // Input data
    int n = 8;
    double x[] = { 1, 2, 3, 4, 4, 3, 2, 1 };
    double y[] = { 0 };

    // Compute FFT
    fft(x, y, n, 0);

    // Print results
    printf("FFT result:\n");
    for (int i = 0; i < n; i++)
        printf("%f + %fi\n", x[i], y[i]);

    return 0;
}

该示例代码使用了递归的Cooley-Tukey算法来实现FFT，可以计算长度为2的幂次方的输入序列的FFT。在main函数中，输入数据通过数组x和y给出，计算结果保存在这两个数组中。可以根据自己的需求修改输入数据和输出格式。