y=37448305898981035/25966066601558016*tan(7489661179796207/2814749767106560*t-7923735992701864684306921345231/5070602400912917605986812821504)+6.085*2;

%计算竖直距离 trapz(t,y) t=0:330; %水平速度

y =-40./(369./5.*t+1447837414805163./4503599627370496);

%计算水平距离 trapz(t,y)

附录6

%给出数据

p1=1025; %水的密度 p2=740; %物体密度 g=10;

c2=0.1; %上举力系数

s=0.09; %黑匣子的迎水面积 v=0.027; %黑匣子的体积 a=(p1-p2)*g/p1; b=c2*p2*s/2/p1/v; v2=177.7179;

c=(sqrt(a/b)+v2)/(sqrt(a/b)-v2);%竖直方向初速度

%当时间趋于正无穷时的速度

w=sqrt(2*(p1-p2)*v*g/c2/p2/s);

%计算竖直距离 t=0:0.001:1000;

y=w*(t-1/sqrt(a*b)*log((1+exp(-c))/(1+exp(-2*sqrt(a*b)-c)))); t(find(abs(y-1000)<0.01))

%ans =209.2530 209.2540 209.2550 209.2560 y=1000

%ans = 417.2790 417.2800 417.2810 417.2820 y=2000

%ans=625.3050 625.3060 625.3070 625.3080 y=3000

附录7

%给出数据

c1=2.2; %拖曳力系数 p1=1025; p2=740; s=0.09; v=0.027;

u=0.25; %水流速度

v1=124.4228; %水平方向初速度; aa=p2*c1*s/2/p1/v; bb=1/(u-v1);

%计算水平距离

t=209.2545 %此处的t由附录5中计算所得的四个时间取平均所得 x=u*t-1/aa*log(abs(aa*t+bb))+1/aa*log(abs(bb)); %x=48.1058 y=1000 %x=99.8516 y=2000 %x=151.7053 y=3000