clear;close all;clc;
x=linspace(-1500,1500,151);
h=linspace(990,990,151);
d=sqrt(x.^2+h.^2);
vp1=1500;vp2=1800;
t=d./vp1;
plot(x,t,'-r')
hold on;
h11=linspace(1010,1010,151);
d11=sqrt(x.^2+h11.^2);
t11=d./vp1;
plot(x,t11)
hold on;
% x1=solve('25*x1^2*((0-2*x1)^2+990^2)=36*(0-2*x1)^2*(x1^2+2000^2)','x1');
% x1=vpa(x1,5);
% x1=real(x1);
% x1=double(x1);
%
% for j=1:4
% if x1(j,1)<750&&x1(j,1)>0
% x1=x1(j,1);
% break;
% end
% end
% fid=fopen('lc123.txt','at+');
% fprintf(fid,'%f\n',x1);
% fclose(fid);
a=load('lc123.txt');
a=a';
b=zeros(1,75);
for i=1:75
b(1,i)=a(1,76-i);
end
x1=[a,b];
h1=linspace(2000,2000,151);
d1=sqrt(x1.^2+h1.^2);
k=linspace(1500,0,76);
x0=k-2*a;
b1=zeros(1,75);
for i=1:75
b1(1,i)=x0(1,76-i);
end
x0=[x0,b1];
d2=sqrt(x0.^2+h.^2);
t1=d1./vp1+d2./vp2;
plot(x,t1,'-b*')
axis([-1500 1500 0 2.1]);
hold off;
title('海底检波器接收到的直达波、反射P波时距曲线');
xlabel('距离/m');
ylabel('时间/s');
grid on;
legend('直达波时距曲线','直达波炮点鬼波时距曲线','反射P波时距曲线');
clear;close all;clc;
x=linspace(-1500,1500,151);
h=linspace(990,990,151);
d=sqrt(x.^2+h.^2);
vp1=1500;vp2=1800;
t=d./vp1;
plot(x,t,'-r')
hold on;
h11=linspace(1010,1010,151);
d11=sqrt(x.^2+h11.^2);
t11=d./vp1;
plot(x,t11)
hold on;
% x1=solve('25*x1^2*((0-2*x1)^2+990^2)=36*(0-2*x1)^2*(x1^2+2000^2)','x1');
% x1=vpa(x1,5);
% x1=real(x1);
% x1=double(x1);
%
% for j=1:4
% if x1(j,1)<750&&x1(j,1)>0
% x1=x1(j,1);
% break;
% end
% end
% fid=fopen('lc123.txt','at+');
% fprintf(fid,'%f\n',x1);
% fclose(fid);
a=load('lc123.txt');
a=a';
b=zeros(1,75);
for i=1:75
b(1,i)=a(1,76-i);
end
x1=[a,b];
h1=linspace(2000,2000,151);
d1=sqrt(x1.^2+h1.^2);
k=linspace(1500,0,76);
x0=k-2*a;
b1=zeros(1,75);
for i=1:75
b1(1,i)=x0(1,76-i);
end
x0=[x0,b1];
d2=sqrt(x0.^2+h.^2);
t1=d1./vp1+d2./vp2;
plot(x,t1,'-b*')
axis([-1500 1500 0 2.1]);
hold off;
title('海底检波器接收到的直达波、反射P波时距曲线');
xlabel('距离/m');
ylabel('时间/s');
grid on;
legend('直达波时距曲线','直达波炮点鬼波时距曲线','反射P波时距曲线');
clear;close all;clc;
x=linspace(-1500,1500,151);
h=linspace(990,990,151);
d=sqrt(x.^2+h.^2);
vp1=1500;vp2=1800;
t=d./vp1;
plot(x,t,'-r')
hold on;
h11=linspace(1010,1010,151);
d11=sqrt(x.^2+h11.^2);
t11=d./vp1;
plot(x,t11)
hold on;
% x1=solve('25*x1^2*((0-2*x1)^2+990^2)=36*(0-2*x1)^2*(x1^2+2000^2)','x1');
% x1=vpa(x1,5);
% x1=real(x1);
% x1=double(x1);
%
% for j=1:4
% if x1(j,1)<750&&x1(j,1)>0
% x1=x1(j,1);
% break;
% end
% end
% fid=fopen('lc123.txt','at+');
% fprintf(fid,'%f\n',x1);
% fclose(fid);
a=load('lc123.txt');
a=a';
b=zeros(1,75);
for i=1:75
b(1,i)=a(1,76-i);
end
x1=[a,b];
h1=linspace(2000,2000,151);
d1=sqrt(x1.^2+h1.^2);
k=linspace(1500,0,76);
x0=k-2*a;
b1=zeros(1,75);
for i=1:75
b1(1,i)=x0(1,76-i);
end
x0=[x0,b1];
d2=sqrt(x0.^2+h.^2);
t1=d1./vp1+d2./vp2;
plot(x,t1,'-b*')
axis([-1500 1500 0 2.1]);
hold off;
title('海底检波器接收到的直达波、反射P波时距曲线');
xlabel('距离/m');
ylabel('时间/s');
grid on;
legend('直达波时距曲线','直达波炮点鬼波时距曲线','反射P波时距曲线');
clear;close all;clc;
x=linspace(-1500,1500,151);
h=linspace(990,990,151);
d=sqrt(x.^2+h.^2);
vp1=1500;vp2=1800;
t=d./vp1;
plot(x,t,'-r')
hold on;
h11=linspace(1010,1010,151);
d11=sqrt(x.^2+h11.^2);
t11=d./vp1;
plot(x,t11)
hold on;
% x1=solve('25*x1^2*((0-2*x1)^2+990^2)=36*(0-2*x1)^2*(x1^2+2000^2)','x1');
% x1=vpa(x1,5);
% x1=real(x1);
% x1=double(x1);
%
% for j=1:4
% if x1(j,1)<750&&x1(j,1)>0
% x1=x1(j,1);
% break;
% end
% end
% fid=fopen('lc123.txt','at+');
% fprintf(fid,'%f\n',x1);
% fclose(fid);
a=load('lc123.txt');
a=a';
b=zeros(1,75);
for i=1:75
b(1,i)=a(1,76-i);
end
x1=[a,b];
h1=linspace(2000,2000,151);
d1=sqrt(x1.^2+h1.^2);
k=linspace(1500,0,76);
x0=k-2*a;
b1=zeros(1,75);
for i=1:75
b1(1,i)=x0(1,76-i);
end
x0=[x0,b1];
d2=sqrt(x0.^2+h.^2);
t1=d1./vp1+d2./vp2;
plot(x,t1,'-b*')
axis([-1500 1500 0 2.1]);
hold off;
title('海底检波器接收到的直达波、反射P波时距曲线');
xlabel('距离/m');
ylabel('时间/s');
grid on;
legend('直达波时距曲线','直达波炮点鬼波时距曲线','反射P波时距曲线');clear;close all;clc;
x=linspace(-1500,1500,151);
h=linspace(990,990,151);
d=sqrt(x.^2+h.^2);
vp1=1500;vp2=1800;
t=d./vp1;
plot(x,t,'-r')
hold on;
h11=linspace(1010,1010,151);
d11=sqrt(x.^2+h11.^2);
t11=d./vp1;
plot(x,t11)
hold on;
% x1=solve('25*x1^2*((0-2*x1)^2+990^2)=36*(0-2*x1)^2*(x1^2+2000^2)','x1');
% x1=vpa(x1,5);
% x1=real(x1);
% x1=double(x1);
%
% for j=1:4
% if x1(j,1)<750&&x1(j,1)>0
% x1=x1(j,1);
% break;
% end
% end
% fid=fopen('lc123.txt','at+');
% fprintf(fid,'%f\n',x1);
% fclose(fid);
a=load('lc123.txt');
a=a';
b=zeros(1,75);
for i=1:75
b(1,i)=a(1,76-i);
end
x1=[a,b];
h1=linspace(2000,2000,151);
d1=sqrt(x1.^2+h1.^2);
k=linspace(1500,0,76);
x0=k-2*a;
b1=zeros(1,75);
for i=1:75
b1(1,i)=x0(1,76-i);
end
x0=[x0,b1];
d2=sqrt(x0.^2+h.^2);
t1=d1./vp1+d2./vp2;
plot(x,t1,'-b*')
axis([-1500 1500 0 2.1]);
hold off;
title('海底检波器接收到的直达波、反射P波时距曲线');
xlabel('距离/m');
ylabel('时间/s');
grid on;
legend('直达波时距曲线','直达波炮点鬼波时距曲线','反射P波时距曲线');
clear;close all;clc;
x=linspace(-1500,1500,151);
h=linspace(990,990,151);
d=sqrt(x.^2+h.^2);
vp1=1500;vp2=1800;
t=d./vp1;
plot(x,t,'-r')
hold on;
h11=linspace(1010,1010,151);
d11=sqrt(x.^2+h11.^2);
t11=d./vp1;
plot(x,t11)
hold on;
% x1=solve('25*x1^2*((0-2*x1)^2+990^2)=36*(0-2*x1)^2*(x1^2+2000^2)','x1');
% x1=vpa(x1,5);
% x1=real(x1);
% x1=double(x1);
%
% for j=1:4
% if x1(j,1)<750&&x1(j,1)>0
% x1=x1(j,1);
% break;
% end
% end
% fid=fopen('lc123.txt','at+');
% fprintf(fid,'%f\n',x1);
% fclose(fid);
a=load('lc123.txt');
a=a';
b=zeros(1,75);
for i=1:75
b(1,i)=a(1,76-i);
end
x1=[a,b];
h1=linspace(2000,2000,151);
d1=sqrt(x1.^2+h1.^2);
k=linspace(1500,0,76);
x0=k-2*a;
b1=zeros(1,75);
for i=1:75
b1(1,i)=x0(1,76-i);
end
x0=[x0,b1];
d2=sqrt(x0.^2+h.^2);
t1=d1./vp1+d2./vp2;
plot(x,t1,'-b*')
axis([-1500 1500 0 2.1]);
hold off;
title('海底检波器接收到的直达波、反射P波时距曲线');
xlabel('距离/m');
ylabel('时间/s');
grid on;
legend('直达波时距曲线','直达波炮点鬼波时距曲线','反射P波时距曲线');
clear;close all;clc;
x=linspace(-1500,1500,151);
h=linspace(990,990,151);
d=sqrt(x.^2+h.^2);
vp1=1500;vp2=1800;
t=d./vp1;
plot(x,t,'-r')
hold on;
h11=linspace(1010,1010,151);
d11=sqrt(x.^2+h11.^2);
t11=d./vp1;
plot(x,t11)
hold on;
% x1=solve('25*x1^2*((0-2*x1)^2+990^2)=36*(0-2*x1)^2*(x1^2+2000^2)','x1');
% x1=vpa(x1,5);
% x1=real(x1);
% x1=double(x1);
%
% for j=1:4
% if x1(j,1)<750&&x1(j,1)>0
% x1=x1(j,1);
% break;
% end
% end
% fid=fopen('lc123.txt','at+');
% fprintf(fid,'%f\n',x1);
% fclose(fid);
a=load('lc123.txt');
a=a';
b=zeros(1,75);
for i=1:75
b(1,i)=a(1,76-i);
end
x1=[a,b];
h1=linspace(2000,2000,151);
d1=sqrt(x1.^2+h1.^2);
k=linspace(1500,0,76);
x0=k-2*a;
b1=zeros(1,75);
for i=1:75
b1(1,i)=x0(1,76-i);
end
x0=[x0,b1];
d2=sqrt(x0.^2+h.^
评论0