地震资料数字处理程序(3)

main() { void conv(); void autocorr(); double temp,bo[73],ao[25],g[49],gr[49]; float z; double b[25],rbb[25],d[25],y[25],c[49]; int i; FILE *fp,*fp1,*fp2; fp=fopen("wavelet.dat","r"); for(i=0;i<25;i++) { fscanf(fp,"%f",&z); b[i]=z;

//printf("%f\n",b[i]); } autocorr(b,25,b,25,rbb); for(i=0;i<25;i++) { if(i==0)d[i]=1.0; else d[i]=0.0; } //求y(t)

tlvs(rbb,25,d,y); //printf("it is ok !\n"); fp1=fopen("yt.dat","w"); for(i=0;i<25;i++) { fprintf(fp1,"%f\n",y[i]); }

//求bo(0) conv(b,25,y,25,c,49);//c(t)=b(t)*g(t) temp=0.0; for(i=0;i<49;i++) { temp=temp+c[i]*c[i]; printf("%f\n",c[i]); }

bo[0]=1.0/sqrt(temp);