fft c语言

FFT是一种快速傅里叶变换算法，可以高效地计算数字信号的频谱。在C语言中，可以使用库函数进行FFT计算，例如使用FFTW库。

以下是一个使用FFTW库进行FFT计算的示例程序：

#include <stdio.h>
#include <stdlib.h>
#include <fftw3.h>

int main() {
    int n = 8; // FFT长度
    double *in, *out;
    fftw_plan plan;
    int i;

    // 分配内存空间
    in = (double*) fftw_malloc(sizeof(double) * n);
    out = (double*) fftw_malloc(sizeof(double) * n);

    // 初始化输入信号
    for (i = 0; i < n; i++) {
        in[i] = i;
    }

    // 创建FFT计算计划
    plan = fftw_plan_r2r_1d(n, in, out, FFTW_R2HC, FFTW_ESTIMATE);

    // 执行FFT计算
    fftw_execute(plan);

    // 输出结果
    for (i = 0; i < n; i++) {
        printf("%f\n", out[i]);
    }

    // 释放内存空间
    fftw_destroy_plan(plan);
    fftw_free(in);
    fftw_free(out);

    return 0;
}

这个示例程序演示了如何使用FFTW库进行FFT计算。其中，fftw_plan_r2r_1d函数创建了一个计算计划，FFTW_R2HC表示实数序列到复数序列的变换，FFTW_ESTIMATE表示使用默认的计算策略。fftw_execute函数执行计算计划，计算结果保存在out数组中。最后，需要释放内存空间和计算计划。