matlab编程合成地震记录

clear;clc;%s1=input('请输入文件名:','s');fid=fopen('yy-10.txt','r');c1=fscanf(fid,'%f');N=length(c1);fori=2:3:Nk=(i-2)/3+1;deltt(k)=c1(i);vv(k)=1000000/deltt(k);rr(k)=0.31*vv(k)^(1/4);Z(k)=vv(k)*rr(k);endn1=N/3;dp=360.2:0.2:2303fork=1:n1-1R(k)=(Z(k+1)-Z(k))/(Z(k+1)+Z(k));endfigure(9);plot(dp,R);%=============================================================%对反射系数序列进行低通滤波%==============================================================r1=fft(R);r1(1001:8716)=0.0;figure(10);plot(abs(r1));r2=ifft(r1);R1=real(r2);figure(11);plot(dp,R1);fori=1:n1-1if(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;fori=rl:xl-rl-1x1(i-rl+1)=x(i);endfigure(2);plot(dp,x1);fp=fopen('I:\showData\wlgqgg30.Dat','w');fori=1:8forj=1:n1-1fwrite(fp,x1(j),'float32');endendfclose(fp);fori=1:3:Nk=(i-1)/3+1;depth(k)=c1(i);end