以下是实现FFT算法的C语言代码:
#include <math.h>
#include <complex.h>
#define N 128 // 数组长度
void FFT(double complex *x, int n)
{
int i, j, k, m;
double complex w, t;
for (m = n; m >= 2; m >>= 1)
{
for (i = 0; i < m / 2; i++)
{
w = cexp(-2 * I * M_PI * i / m);