信号分析与处理上机试验报告 .docVIP

信号分析与处理上机试验报告 电气化0501 李科 200502020209 实验一 快速傅利叶变换的软件实现 1. 设计程序源代码如下所示： #include stdio.h #include math.h #include conio.h /*清屏命令clrscr()需加载的预处理函数*/ #define swap(a,b) temp=(a);(a)=(b);(b)=temp /*宏定义的交换函数*/ /*快速傅利叶变换程序，数组A、B分别是带变换序列的实部和虚部ap是变换的类别，ap=1是做FFT; ap=-1是做IFFT */ void fft(float A[],float B[],unsigned M,int ap) {unsigned long N,I,J,K,L,LE,LE1,P,Q,R; int i; float Wr,Wi,W1r,W1i,WTr,WTi,theta,Tr,Ti,temp; N=pow(2,M); /*N=2M 是序列的总长度*/ J=0; for(I=0;IN-1;I++) /*码位倒置*/ {if(JI) {swap(A[I],A[J]); swap(B[I],B[J]);} K=N1; /*即k=k/2 */ while(K=2J=K) {J-=K; K=1;} J+=K;} for(L=1;L=M;L++) /*外层循环由级数L控制，执行M次*/ {LE=1L; /* LE=2L是群间隔*/ LE1=LE/2; /* LE1=2L-1是每个群的系数W数目*/ Wr=1.0; Wi=0.0; theta=(-1)*ap*3.1415926536/LE1; /*ap用于在IFFT时改变极性*/ W1r=cos(theta); W1i=sin(theta); for(R=0;RLE1;R++) /*中层循环由群系数控制，执行次*/ {for(P=R;PN-1;P+=LE) /*R是群系数的编号，P、Q是基本蝶形运 算两个输入数据在数组中的编号，循环每 次完成同一个系数的蝶形运算*/ {Q=P+LE1; Tr=Wr*A[Q]-Wi*B[Q]; Ti=Wr*B[Q]+Wi*A[Q]; /*Tr、Ti是的实部和虚部*/ A[P]+=Tr; /*即*/ B[P]+=Ti;} A[Q]=A[P]-Tr; /*即*/ B[Q]=B[P]-Ti; WTr=Wr; /*Wr、Wi是的实部和虚部*/ WTi=Wi; /*下面用和差化积公式求的实虚部Wr、Wi*/ Wr=WTr*W1r-WTi*W1i; /**/ Wi=WTr*W1i+WTi*W1i; /**/ } } if(ap==-1){ /*ap=-1是做IFFT，应使变换后的序列除以N*/ for(i=0;iN;i++){ A[i]/=N; B[i]/=N; } } return; } main(){ int i,t,ap; float A[100],B[100]; unsigned M; clrscr(); /*清屏*/ printf(\nplease input bian huan lei bie(FFT:1 or IFFT:-1):); scanf(%d,ap); /*输入序列变换的类别*/ printf(\nplease input jishu M:); /*M是级数*/ scanf(%d,M); t=pow(2,M); for(i=0;it;i++){ printf(\nplease input x(%d):(Re Im),i); scanf(%f %f,A[i],B[i]); /* 输入离散序列的实、虚部*/ } fft(A,B,M,ap); /*进行FFT(或IFFT)*/ printf(\n*************Result**************); for(i=0;it;i++) printf(\nX(%d): %.5f+j%.5f,i,A[i],B[i]); /*输出变换的结果*/ } 2. 实例验证： 2.1 用课本序列计算FFT得： please input bian huan lei bie(FFT:1 or IFFT:-1