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