fft算法代码注释及流程框图.docVIP

基2的DIT蝶形算法源代码及注释如下： /************FFT***********/ //整个程序输入和输出利用同一个空间x[N]，节约空间 #include stdio.h #include math.h #include stdlib.h #define N 1000 //定义输入或者输出空间的最大长度 typedef struct { double real; double img; }complex; //定义复数型变量的结构体 void fft(); //快速傅里叶变换函数声明 void initW(); //计算W(0)~W(size_x-1)的值函数声明 void change(); //码元位置倒置函数函数声明 void add(complex,complex,complex *); /*复数加法*/ void mul(complex,complex,complex *); /*复数乘法*/ void sub(complex,complex,complex *); /*复数减法*/ void divi(complex,complex,complex *); /*复数除法*/ void output(); /*输出结果*/ complex x[N],*W; /*输出序列的值*/ int size_x=0; /*输入序列的长度，只限2的N次方*/ double PI; //pi的值 int main() { int i; system(cls); PI=atan(1)*4; printf(Please input the size of x:\n); /*输入序列的长度*/ scanf(%d,size_x); printf(Please input the data in x[N]:(such as:5 6)\n); /*输入序列对应的值*/ for(i=0;isize_x;i++) scanf(%lf %lf,x[i].real,x[i].img); initW(); //计算W(0)~W(size_x-1)的值 fft(); //利用fft快速算法进行DFT变化 output(); //顺序输出size_x个fft的结果 return 0; } /*进行基-2 FFT运算,蝶形算法。这个算法的思路就是， 先把计算过程分为log(size_x)/log(2)-1级（用i控制级数）； 然后把每一级蝶形单元分组（用j控制组的第一个元素起始下标）； 最后算出某一级某一组每一个蝶形单元（用k控制个数，共l个）。*/ void fft() { int i=0,j=0,k=0,l=0; complex up,down,product; change(); //实现对码位的倒置 for(i=0;ilog(size_x)/log(2);i++) //循环算出fft的结果 {