%中點放炮水平界面?波的時距曲線
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;