第四章习题与答案.docVIP

解: ⑴ 直接计算: 复乘所需时间: 复加所需时间: ⑵用FFT计算: 复乘所需时间: 复加所需时间: 由（1）式可得的路径，如下表所示: k 0 1 2 3 4 5 6 7 8 9 0.8 0.67 0.56 0.46 0.39 0.32 0.27 0.22 0.19 0.16 arg [] 8. 请用C语言编写程序: (1) 按频率抽取的FFT算法 (2) 分裂基FFT算法 基- 2 FFT ( 频率抽取DIF 法) 算法程序 /*Free_Copy*/ /* C语言编写的频率抽取FFT算法(最大计算64点）*/ /* 输入: 序列点数、序列值 * / /* 输出: 序列FFT变换后的数值及反变换(应与原序列相同 ) */ #include conio.h #include math.h #include stdio.h #define N 64 #define PI 3.1415926 #define w0 (0.125*PI) #define Cmul(a,b,c) a.x=b.x*c.x-b.y*c.y;a.y=b.x*c.y+b.y*c.x; #define Cequal(a,b) a.x=b.x;a.y=b.y; #define Cadd(a,b,c) a.x=b.x+c.x;a.y=b.y+c.y; #define Csub(a,b,c) a.x=b.x-c.x;a.y=b.y-c.y; #define Wn(w,r) w.x=cos(2.0*PI*r/n);w.y=-sin(2.0*PI*r/n); struct comp { float x; float y; }; void main() { int i,j,nu2,nm1,n,m,le,le1,k,ip,z; int flag,f,n1; struct comp a[N],t,t1,w,d; float a_ipx,m1; printf(\nThis program is about FFT by DIF way. ); printf(\nplease enter N : ); scanf(%d,n1); n=n1; m1=log(n1)/log(2); m=log(n1)/log(2); if (m!=m1) n=pow(2,m+1); for(i=0;in;i++) {a[i].x=a[i].y=0.0;} printf(\n); for(i=0;in1;i++) { printf(\nplease enter data(%d)_[Re]: ,i); scanf(%f,a[i].x); printf(\nplease enter data(%d)_[Im]: ,i); scanf(%f,a[i].y); } for(z=0;z=1;z++) { flag=-1; for (m=(log(n)/log(2));m=1;m--) { le=pow(2,m); flag++; le1=le/2; for( j=0;jle1;j++) { for (i=j;i=(n-1);i+=le) { ip=i+le1; Cequal(t,a[i]); Cequal(t1,a[ip]); f=(int) (i*pow(2,flag))%n; Wn(w,f); Cadd(a[i],t,t1); Csub(a[ip],t,t1); a_ipx=a[ip].x; if (z==1) { w.y*=-1; } a[ip].x=a[ip].x*w.x-a[ip].y*w.y; a[ip].y=