C语言写一个DSP的FFT变换算法

以下是一个简单的C语言实现FFT(快速傅里叶变换)算法的示例代码：

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

#define PI 3.14159265358979323846

// 计算离散傅里叶变换(DFT)
void dft(double *in, double *out, int n)
{
    int i, j;
    double re, im;

    for (i