Implement FFT algorithm with C language

"用C语言实现快速傅里叶变换算法 /*****************fft程序***************/#包括“typedef.h”#包括“math.h” 结构compx EE（结构compx b1，结构compx b2）{结构compx b3；b3.real=b1.real*b2.real-b1.imag*b2.imag；b3.imag=b1.real*b2.imag+b1.imag*b2.real；返回（b3）；} void FFT（struct compx*xin，int N）{int f，m，nv2，nm1，i，k，j=1，l；/*int f，m，nv2，nm1，i，k，j=N/2，l*/结构组件v、w、t；nv2=N/2；f=N；对于（m=1(f=f/2）=1.m++{；}nm1=N-1； /*变址运算*/对于（i=1；i<=nm1；i++{if（i{t=xin[j]；xin[j]=xin[i]；xin[i]=t；}k=nv2；而（k{j=j-k；k=k/2；}j=j+k；}{int le，lei，ip；浮点数；对于（l=1；l<=m；l++{le=pow（2，l）； // 这里用的是L而不是1 !!!!lei=le/2；pi=3.14159；v.real=1.0；v.imag=0.0；w.real=cos（pi/lei）；w.imag=-sin（pi/lei）；对于（j=1；j<=lei；j++{/*double p=pow（2，m-l）*j；双ps=2*pi/N*p；w.real=cos（ps）；w.imag=-sin（ps）*/对于（i=j；i<=N；i=i+le）{/*w.real=cos（ps）；w.imag=-sin（ps）*/ip=i+lei；t=EE（xin[ip]，v）；xin[ip].real=xin[i].real-t.real；xin[ip].imag=xin[i].imag-t.imag；xin[i].real=xin[i].real+t.real；xin[i].imag=xin[i].imag+t.imag；}v=EE（v，w）；}}}返回；} /*****************主程序********************/ #include#include#include#include“typedef.h” 浮动结果[257]；结构复合材料[257]；int Num=256；常量浮点pp=3.14159； main（）{inti=1；对于（；i<0x101；i++{s[i].real=sin（pp*i/32）；s[i].imag=0；}FFT（s，Num）；对于（i=1；i<0x101；i++{result[i]=sqrt（pow（s[i].real，2）+pow（s[i].imag，2））； }}"

Our other product:

Professional lightning protection system for Tower	Professional FM Radio Station Equipment Package	FMUSER FU-1000C 1000W FM Transmitter With USB Audio Input Interface
FMUSER 200W FM Radio Broadcasting Package	Hotel IPTV Solution	FMUSER FSN-1000T V5.0 1KW FM Radio Broadcast Transmitter