基 于基 于2时间抽取 时间抽取 时间抽取 FFT FFT流图设计.pdf

班级：电子131 姓名：ZWQ 学号：1319200035 基于基2 时间抽取FFT 流图设计的用于计算1024 点序列的FFT 源程序 #includestdio.h #includemath.h #define pi 3.1415926 #define N 10 /*设置二进制位数*/ #define K 10/*设置总级数,级数应该是小于二进制位数*/ #define L 1024 /*设置点数*/ /****************************************/ struct Fushu{ float R;float I; };/***定义结构体*****/ /*计算2 的i 次方*/ int bcheng( int i ) { int j; if(i==0) j=1; else for(j=1;i;i--)j=j*2; return j; } /****产生WN 数组上标********/ void WM(int e[],int k)/****k 是级数********/ { int i,j,B,b; for(i=0;iL;) { B= bcheng(k); b= bcheng(k-1); if((i%B)==b) {for(j=0;jb;j++,i++)e[i]=j*bcheng(N-k);} else {e[i]=0;i++;} } } /********************************/ /***计算e 的次方***/ void WnmRI(int n,int e[],struct Fushu Wn[])/*****(总点数,产生上标的数组,Wn 数组)*********/ { int i; for(i=0;iL;i++) {Wn[i].R=cos(2.0*pi*(e[i])/n); Wn[i].I=sin(-2.0*pi*(e[i])/n); } } /********************************/ /*******复数乘积***************/ void Fushucheng(struct Fushu a[],struct Fushu b[]) { int i;float R,I; for(i=0;iL;i++) { R=a[i].R*b[i].R-a[i].I*b[i].I; I=a[i].R*b[i].I+a[i].I*b[i].R; a[i].R=R; a[i].I=I; } } /********************************/ void nextarray(struct Fushu x[],int k)/***(操作数组,级数)******/ { int i,j,B,b; float R,I; for(i=0;iL;) { B= bcheng(k); b= bcheng(k-1); if((i%B)==0) { for(j=0;jb;j++,i++) { R=x[i].R+x[i+bcheng(k-1)].R; I=x[i].I+x[i+bcheng(k-1)].I; x[i+bcheng(k-1)].R=x[i].R-x[i+bcheng(k-1)].R; x[i+bcheng(k-1)].I=x[i].I-x[i+bcheng(k-1)].I; x[i].R=R; x[i].I=I; } } else { i++;} } } /********************************/ /****************序列号转换***********************/ int change( int k ) {