matlab编程合成地震记录(共2页).doc

《matlab编程合成地震记录(共2页).doc》由会员分享，可在线阅读，更多相关《matlab编程合成地震记录(共2页).doc（2页珍藏版）》请在淘文阁 - 分享文档赚钱的网站上搜索。

1、精选优质文档-倾情为你奉上clear;clc;%s1=input(请输入文件名: ,s);fid=fopen(yy-10.txt,r);c1=fscanf(fid,%f);N=length(c1);for i=2:3:N k=(i-2)/3+1; deltt(k)=c1(i); vv(k)=/deltt(k); rr(k)=0.31*vv(k)(1/4); Z(k)=vv(k)*rr(k);endn1=N/3;dp=360.2:0.2:2303for k=1:n1-1 R(k)=(Z(k+1)-Z(k)/(Z(k+1)+Z(k); endfigure(9);plot(dp,R);%=%对反射系

2、数序列进行低通滤波%=r1=fft(R);r1(1001:8716)=0.0;figure(10);plot(abs(r1);r2=ifft(r1);R1=real(r2);figure(11);plot(dp,R1);for i=1:n1-1 if(abs(R1(i)0.01) R1(i)=0.0; endendfigure(12);plot(dp,R1);f=30;wl=50;t=-wl:wl;deltt=0.002;b=(1-2*(pi*f*t*deltt).2).*exp(-(pi*f*t*deltt).2);figure(1);plot(b);x=conv(R1,b);xl=length(x);rl=wl;for i=rl:xl-rl-1 x1(i-rl+1)=x(i);endfigure(2);plot(dp,x1);fp=fopen(I:showDatawlgqgg30.Dat,w);for i=1:8 for j=1:n1-1 fwrite(fp,x1(j),float32); endendfclose(fp);for i=1:3:N k=(i-1)/3+1; depth(k)=c1(i);end专心-专注-专业