微分方程鸭子过河

1、鸭子过河设河边点O的正对岸为点A,河宽OA=h(图)，水流速度为a,有一鸭子从点A游向点O,设鸭子(在静水中)的游速为b (b>a),且鸭子游动的方向始终朝着点Oo设h=10m, a=1m/s, b=2m/s,用数值法求渡河所需时间、任意时刻鸭子的位置及游动曲线。建立任意时刻鸭子的位置和鸭子游动的数 学模型，并求其解析解。1 .模型的假设为了使问题确定和简化，实际上已经作了如下假设：假设河宽固定，设为h,且两岸为平行直线；鸭子游速为 b及水流速度a均为常数；鸭子游动的方向始终指向Oo2 .模型的建立和求解取O为坐标原点，河岸朝顺水方向为x轴，y轴指向对岸，如图所示。设时刻t鸭子位于点 P

2、 (x, y),设起点坐标(x, y) = (0, h),终点坐标(0, 0),设0为鸭子 速度方向与x轴正向间的夹角，b (bcos ,bsin )bOPa (a,0) , v a b于是鸭子游动的迹线满足：x(0)=0 , y(0)=h(1)模型的数值解实际上，从上述方程不能求得x(t) , y(t)的解析式，但在参数确定的情况下，可以通过数值解得到任意时刻x(2)=y ,编鸭子的位置。设 x=(x(1) , x(2) , x(1)=x 写如下的函数M文件：%鸟子过河、渡河function dx=duhe(t,x) %为duhe的函数M文件a=1;b=2;s=sqrt(x(1)A2+x(2

3、)A2);dx=a-b*x/s;-b*x/s;%在编写运行程序时，须设定时间t(可见，鸭子的渡河时间在 7s之间)ts=0:7;x0=0,10;t,x=ode45(duhe,ts,x0); %t,xplot(t,x),gridgtext('x(t)'),gtext('y(t)'),pause % plot(x(:,1),x(:,2),grid,建立名以向量形式表示方程组的起点及终点步长，可大致估计静水中的渡河时间，并作试探。%* y的初始值调用ode45计算姗出 t , x(t) , y(t)%按照数值输出作x(t) , y(t)的图形利用鼠标确定字符串位置%乍

4、y(t)的图形gtext('x'),gtext('y')得到的数彳t结果x(t) , y(t)为鸭子的位置列入表。x(t) , y(t)及y(x)的图形见图(a)和(b)表h=10 , a=1, b=2时的数值解tx(t)y(t)tx(t)y(t)图(a)和图(b)(2)模型的解析解为了得到更精确的运动轨迹，还必须对模型作进一步分析以得到其解析解。鸭子运动速度为:故有：由此得到微分方程：dx Vxdy Vy22a x ybyx一,x(h)=0y求解此齐次微分方程得到鸭子游动的轨迹方程为: a a1 1 -y b y b-,0<y<h (具体求解参见附

5、录(1)h h采用下列Matlab程序，我们可以画出鸭子运动的轨迹(图)。h=10;a=1;b=2;y=h:0;x=h/2*(y./h).A(1-a/b)-(y./h).A(1+a/b);plot(x,y,'bO-') legend('duck') xlabel('X');ylabel('Y');图鸭子运动的轨迹鸭子游动曲线轨迹的弧长可以用公式dsJix2dy求出,也可以用数值方法求解。3 .对解以及问题的进一步讨论关于解可以作进一步分析：如果 b<a,由上述轨迹方程当y-0,得到x-8。因此，这中情况下鸭子是不可能到达对岸

6、的，这与鸭子运动的力学分析结果是一致的。syms y;limit(10/2*(y/10)A(1-2)-(y/10)A(1+2),y,0,'left')syms y;limit(10/2*(y/10)A(1-2)-(y/10)A(1+2),y,0,'right') 结果分别为-Inf和Inf。很自然地，还可以探讨如下问题：如果鸭子上岸的地点不超过和对岸下游一定位置(比如与正对 岸距离为l ),鸭子的速度大小与方向不变，问鸭子以怎样的游动方向才能以最少的时间到达上岸地点？鸭子能够按要求到达对岸速度应满足什么条件？如果水流速度变化，进一步可研究2003年全国数学建模竞

7、赛D题：强渡长江。4 .建模过程总结这是一个微分方程应用题，整个解题过程已经包含了建立数学模型的基本内容，即根据问题背景和建模问题作出必要的简化假设一一鸭子速度和水流速度均为常数；用字母和符号表示有关变量(如鸭子速度、水流速度、时间及位置坐标等)；利用相应的物理（或其他）规律一一牛顿力学有关规律，列出微分方程；求解微分方程得到鸭子游动轨迹曲线解析解，此处我们还采用了数值解法得到了任意时刻鸭子的位置（坐标）；解的讨论及推广应用等。参考文献20011李志林，欧宜贵，数学建模及典型案例分析，北京：化学工业出版社，2同济大学应用数学系，高等数学（本科少学时类型） 上册（第二版），北京：高等教育出版社，

8、附录：（1）鸭子游动轨迹方程的求解将得到的微分方程dxaX22yx一化成齐次方程dxdybyydydxaH 1 x yydyb;x-的形式，得y（i-i）令ux一,贝1 x=yu , ydx dyduu y，代入上述方程，得 dydua 2u y vu1 u(1-2)dyb化简并分离变量得(1-3)两端积分，得ln(u Ju2 1) aln y G (其中。为常数)(1-4)b au 1 1 y b C2(其中 C2 eC1)(1-5)x将u一代入上式，得y(1-6)-J- 1 y2C2(其中 C2eC1)y . ya_a由 x(h)=0 # y=h, x=0 代入上式，得 1 h b C2 ,求得 C2 hba将C2hb代入式(1-5),得 a aJu2 1 hby b u(1-7)将上式平方并化简，得a aa