課程設(shè)計(jì)：信號(hào)分析與處理c語言編程(編輯修改稿)
2025-06-12 20:46 本頁面
【文章內(nèi)容簡介】 0}。 for(i=0。i200。i++) { s[i]=A1*exp(B*pow((i*data),2)) * sin(2*PI*f1*i*data) + A2*exp(B*pow((i*data),2)) * sin(2*PI*f2*i*data)。 n[i]=A3*(sin(2*PI*f3*i*data)+cos(2*PI*f4*i*data))。 x[i]=s[i]+n[i]。 } //導(dǎo)出 x[t] fp=fopen(D:\\sh\\py\\,w)。 for(i=0。iG。i++) fprintf(fp,%f\n,x[i])。 fclose(fp)。 //x[t]頻譜計(jì)算 fft(x,x1,Fx,G,1)。 //導(dǎo)出 X(k) fp=fopen(D:\\sh\\py\\,w)。 for(i=0。iG。i++) fprintf(fp,%f\n,Fx[i])。 fclose(fp)。 //處理 h[t] double dataf。 double m1,m2,m3,m4。 dataf=(data*G)。 m1=,m2=。 m3=Gm2,m4=Gm1。 //計(jì)算 ,導(dǎo)出單門 for(i=0。iG。i++) { if(i=int(m1)amp。amp。i=int(m2)) Fh1[i]=1。 else Fh1[i]=0。 } //計(jì)算 ,導(dǎo)出雙門 for(i=0。iG。i++) { if((i=int(m1)amp。amp。i=int(m2))||(i=int(m3)amp。amp。i=int(m4))) Fh2[i]=1。 else Fh2[i]=0。 } fp=fopen(D:\\sh\\py\\,w)。 for(i=0。i256。i++) fprintf(fp,%f\t%f\n,Fh1[i],Fh2[i])。 fclose(fp)。 for(i=0。iG。i++)//單門濾波 { Y1[i]=x[i]*Fh1[i]。 Y1i[i]=x1[i]*Fh1[i]。 Fy1[i]=sqrt(pow(Y1[i],2)+pow(Y1i[i],2))。 } for(i=0。iG。i++)//開雙門 { Y2[i]=x[i]*Fh2[i]。 Y2i[i]=x1[i]*Fh2[i]。 Fy2[i]=sqrt(pow(Y2[i],2)+pow(Y2i[i],2))。} fp=fopen(D:\\sh\\py\\,w)。//導(dǎo)出單雙門 for(i=0。iG。i++) fprintf(fp,%f\t%f\n,Fy1[i],Fy2[i])。 fclose(fp)。 //Y[k]作反變換 IFFT,并且導(dǎo)出實(shí)部，虛步 ifft(Y1,Y1i,G,1)。 ifft(Y2,Y2i,G,1)。 fp=fopen(D:\\sh\\py\\,w)。 for(i=0。iG。i++) fprintf(fp,%f\t%f\t%f\t%f\n,Y1[i],Y1i[i],Y2[i],Y2i[i])。 fclose(fp)。 printf(\nover\n\n)。} 相關(guān)程序代碼 include include include define PI void main() { int i,j。 double A1=,A2=。 double f1=30,f2=50。 double data=,p=100。 double x[100]={0},y[150]={0},xy[249]={0},yx[249]={0}。 double xx[199],yy[299],xyft[249]。 for(i=0。i100。i++) x[i]=A1*sin(2*PI*f1*i*data)。 for(j=0。j150。j++)//計(jì)算 y(t) y[j]=A2 * exp(p*j*data*j*data)*sin(2*PI*f2*j*data)。 correl(x,y,100,150,xy)。 correl(x,x,100,100,xx)。 correl(y,y,150,150,yy)。 correl(y,x,150,100,yx)。 FILE *fp,*fp1,*fp2,*fp3,*fp4,*fp5,*fp6。//導(dǎo)出數(shù)據(jù) int i1,j1,k,l,m,n。 fp1=fopen(D:\\sh\\data\\,w)。 fp2=fopen(D:\\sh\\data\\,w)。 fp3=fopen(D:\\sh\\data\\,w)。 fp4=fopen(D:\\sh\\data\\,w)。 fp5=fopen(D:\\sh\\data\\,w)。 fp6=fopen(D:\\sh\\data\\,w)。 for(i=0。i249。i++) xyft[i]=yx[i]。 fp=fopen(D:\\sh\\data\\,w)。 for(i=0。i249。i++) fprintf(fp,%f\n,xyft[i])。 fclose(fp)。 for(i1=0。i1100。i1++) fprintf(fp1,%f\n,x[i1])。 fclose(fp1)。 for(j1=0。j1150。j1++) fprintf(fp2,%f\n,y[j1])。 fclose(fp2)。 for(k=0。k249。k++) fprintf(fp3,%f\n,xy[k])。 fclose(fp3)。 for(l=0。l199。l++) fprintf(fp4,%f\n,xx[l])。 fclose(fp4)。 for(m=0。m249。m++) fprintf(fp5,%f\n,yx[m])。 fclose(fp5)。 for(n=0。n299。n++) fprintf(fp6,%f\n,yy
環(huán)評(píng)公示相關(guān)推薦
鄂ICP備17016276號(hào)-1