.
…… 66.5 66.51 66.52 66.53 66.54 66.55 66.56 66.57 66.58 66.59 66.6 66.61 66.62 66.63 66.64 66.65 66.66 66.67 …… 0.16685 0.15685 0.14685 0.13685 0.12685 0.11685 0.10685 0.096853 0.086853 0.076853 0.066853 0.056853 0.046853 0.036853 0.026853 0.016853 0.0068532 6.7389e-011 …… -0.0011133 -0.00098384 -0.00086239 -0.000748 -0.00064349 -0.00054603 -0.00045657 -0.00037511 -0.00030164 -0.00023618 -0.00017871 -0.00012924 -8.7773e-005 -5.4301e-005 -2.8827e-005 -1.1351e-005 -1.8739e-006 0 【小结】
小船起初速度较快，相当于在顺水被冲下去。后来快到达终点时较慢，相当于逆水而行。
6.4解析解的图像
先建立函数
function x=Guohe2(y,k)
x=0.5*(-0.01)^(-k)*y.^(1-k)-0.5*(-0.01)^k*y.^(1+k); end
y=0:-0.01:-100; x=Guohe2(y,0.5); plot(x,y)，grid;
title('图4.过河轨迹，解析解') xlabel('x/m') ylabel('y/m')
11页
.
可以看出，由数值解法得到的数据绘成的图与解析解的结果几乎一样。
6.5改变流速v1=0
ts=0:0.01:80; x0=[0,-100];
option=odeset('reltol',1e-6,'abstol',1e-9); [t,x]=ode23s(@Guohe,ts,x0,option,0,2); plot(t,x),grid
title('图1.分位移-时间') xlabel('t/s') ylabel('x,y/m') gtext('x(t)') gtext('y(t)') pause
plot(x(:,1),x(:,2)),grid, title('图2.过河轨迹') xlabel('x/m') ylabel('y/m')
12页
.
6.6 改变流速v1=0.5m/s
程序基本相同，在此略去，只给出图像。
13页
.
6.7 改变流速v1=1.5m/s
14页
.
6.8 改变流速v1=2m/s
15页