傅里叶变换的c语言代码

以下是一个简单的使用C语言实现傅里叶变换的代码示例：

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

#define PI 3.14159265358979323846

// 定义复数结构体
typedef struct {
    double real;
    double imag;
} complex;

// 傅里叶变换函数
void fft(complex *x, int n)
{
    if (n <= 1) {
        return;
    }
    complex *even = (complex*)malloc(n/2 * sizeof(complex));
    complex *odd = (complex*)malloc(n/2 * sizeof(complex));
    for (int i = 0; i < n/2; i++) {
        even[i] = x[2*i];
        odd[i] = x[2*i+1];
    }
    fft(even, n/2);
    fft(odd, n/2);
    for (int i = 0; i < n/2; i++) {
        complex t = {cos(2*PI*i/n), -sin(2*PI*i/n)};
        t.real *= odd[i].real;
        t.imag *= odd[i].imag;
        x[i].real = even[i].real + t.real;
        x[i].imag = even[i].imag + t.imag;
        x[i+n/2].real = even[i].real - t.real;
        x[i+n/2].imag = even[i].imag - t.imag;
    }
    free(even);
    free(odd);
}

int main()
{
    int n = 8;
    complex x[] = {{1, 0}, {2, 0}, {3, 0}, {4, 0}, {5, 0}, {6, 0}, {7, 0}, {8, 0}};
    fft(x, n);
    for (int i = 0; i < n; i++) {
        printf("%f + %fi\n", x[i].real, x[i].imag);
    }
    return 0;
}

这段代码中，我们定义了一个复数结构体，用于表示复数。我们实现了一个递归的FFT算法，用于对输入向量进行傅里叶变换。在main函数中，我们定义了一个长度为8的复数向量，并对其进行傅里叶变换。最后，我们打印了变换后的结果。