实验4常微分方程数值解 下载本文

内容发布更新时间 : 2024/12/23 10:36:41星期一 下面是文章的全部内容请认真阅读。

.

…… 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页