fft c语言程序

fft c语言程序
以下是一个基于C语言的FFT（快速傅里叶变换）程序示例：
```c
#include <stdio.h>
#include <math.h>
#define PI 3.14159265
void FFT(double complex *x, int n ) {
if ( n <= 1 ) return;
double complex even ( n /2 );
double complex odd ( n /2 );
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++) {
double complex t = cexp(-I *2 *PI * i / n ) * odd (i); x[i] = even (i) + t ;
x[i + n /2 ] = even (i) - t ;
}
}
int main() {
double complex x[] = {1 , 2 , 3 , 4 , 5 , 6 , 7 , 8 };
int n = sizeof(x) / sizeof(double complex );
FFT(x, n );
for ( int i =0 ; i < n ; i++) {
printf( "%f + %fi\n" , creal(x[i]), cimag(x[i]));
}
return 0 ;
}
```
这个程序使用了递归的方法来实现快速傅里叶变换。

它接受一个复数数组`x`作为输入，并计算其傅里叶变换。

在主函数中，定义了一个复数数组`x`，并调用`FFT`函数对其进行变换。

最后，输出变换后的结果。

请注意，这个程序只是一个简单的示例，实际的FFT程序可能会更加复杂，并且可能需要根据具体的应用场景进行调整。