各类地震波时距曲线的matlab编码

%中点放炮水平界面 波的时距曲线
clear all
dx=50;x=-500:dx:500;%炮检距
h=100;V1=1500;
theta=25*pi/180;
V2=V1/sin（theta）;
t1=sqrt（x.*x+4*h*h）/V1;%反射波时距曲线
t2=abs（x）./V1;%直达波时距曲线
%折射波时距曲线
xm=2*h*tan（theta）;%求盲区
k=1;
for i=1:length（x）
    if x（i）<-xm
        t3_1（k）=2*h*cos（theta）/V1+abs（x（i））/V2;
        x_1（k）=x（i）;
        k=k+1;
    end
end
k=1;