地震数据处理快速傅里叶变换的c语言代码

以下是一个简单的 C 语言程序，用于执行快速傅里叶变换（FFT）以处理地震数据：

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

#define PI 3.14159265358979323846

// Function to perform the FFT
void FFT(double complex *x, int n)
{
    if (n == 1)
        return;

    double complex *even = malloc(n / 2 * sizeof(double complex));
    double complex *odd = malloc(n / 2 * sizeof(double complex));

    for (int i = 0; i < n / 2; i++)
        even[i] = x[2 * i];
    for (int i = 0; i < n / 2; i++)
        odd[i] = x[2 * i + 1];

    FFT(even, n / 2);
    FFT(odd, n / 2);

    for (int k = 0; k < n / 2; k++)
    {
        double complex t = cexp(-I * PI * k / (double)n) * odd[k];
        x[k] = even[k] + t;
        x[k + n / 2] = even[k] - t;
    }

    free(even);
    free(odd);
}

int main()
{
    // Sample input data
    double complex x[] = {1, 2, 3, 4, 5, 6, 7, 8};

    // Number of data points
    int n = sizeof(x) / sizeof(double complex);

    // Perform the FFT
    FFT(x, n);

    // Print the output
    for (int i = 0; i < n; i++)
        printf("x[%d] = %f + %fi\n", i, creal(x[i]), cimag(x[i]));

    return 0;
}

该程序使用递归实现快速傅里叶变换，并且是基于公式 `e^(-2*pi*i*k/N)` 来计算旋转因子。该代码可以很容易地进行修改，以适应不同的输入数据格式和变换参数。