微分方程建模课堂讨论之一-利用导数的几何及物理意义建模..docVIP

微分方程建模课堂讨论之一 利用导数的几何以及物理意义建模 1．将一个手表系在一条链子上，当沿直线拉此链的一端，在平面上，此手表所描述的轨迹是什么？ 【注】：微积分的发明人莱布尼兹讨论过此问题。 【模型】 HYPERLINK /retype/zoom/2fb98c88d0d233d4b14e6928?pn=1x=0y=0raww=69rawh=28o=jpg_6_0_______type=picaimh=28md5sum=76d92ac9b8ae8340cdb160sign=2111baad35zoom=png=0-741jpg=0-1888 \t _blank ：y(x)=- 【结果】 HYPERLINK /retype/zoom/2fb98c88d0d233d4b14e6928?pn=1x=0y=28raww=65rawh=24o=jpg_6_0_______type=picaimh=24md5sum=76d92ac9b8ae8340cdb160sign=2111baad35zoom=png=0-741jpg=0-1888 \t _blank HYPERLINK /retype/zoom/2fb98c88d0d233d4b14e6928?pn=1x=0y=52raww=67rawh=29o=jpg_6_0_______type=picaimh=29md5sum=76d92ac9b8ae8340cdb160sign=2111baad35zoom=png=0-741jpg=0-1888 \t _blank ：y(x)=x 【求解】：方法？换元法。 aa*arctan，此为解析解。 【预备知识】：（1）积分方法；（2）双曲函数。 2．【小孩和玩具问题】：一个小孩在平面上沿一曲线行走，设次曲线由两个世纪的函数X(t)和Y(t)确定。假设此小孩借助长度为a的硬棒，拉或推某玩具，设x(t)和y(t)是玩具的位置，试确定玩具的轨迹。 【建模分析】： （1）．(X(t)，Y(t))与(x(t)，y(t))之间的距离总是硬棒的长度，于是 (X-x)+(Y-y)=a222 （1） （2）．玩具总是在硬棒的方向上运动，因此，两个位置的差向量是玩具的速度方向的倍数，vT=(x,y)T： ?X-x??x?=λ ? ?,λ0 （2） ?Y-y??y? （3）玩具的速度依赖于小孩的速度向量Vc的方向。小孩的速度Vc在硬棒上的投影是玩具的速度VT的模。 【特例】：假设小孩在半径为a（硬棒的长）的圆上行走。在此特殊 情况下，玩具停留在此圆的圆心，根本不运动。 将式（2）代入式（1）可得， HYPERLINK /retype/zoom/2fb98c88d0d233d4b14e6928?pn=2x=0y=0raww=81rawh=31o=jpg_6_0_______type=picaimh=31md5sum=76d92ac9b8ae8340cdb160sign=2111baad35zoom=png=742-3202jpg=1889-3543 \t _blank λ= （3） 于是 HYPERLINK /retype/zoom/2fb98c88d0d233d4b14e6928?pn=2x=0y=31raww=81rawh=30o=jpg_6_0_______type=picaimh=30md5sum=76d92ac9b8ae8340cdb160sign=2111baad35zoom=png=742-3202jpg=1889-3543 \t _blank x??X-x??= ?y??Y-y? （4） （4）如何利用|VT |=|VC|cosα ？ T ● 标准化差向量(X● 确定V● C T =(X,Y) -x,Y-y) ，可得单位向量w； TC 在w生成的子空间上的投影：V w ； VT=(x,y)=(VCw)w TT （5）数值解法 （6）数值解Matlab程序： 程序1：小孩的运动： function [X,Xs,Y,Ys]=child(t) %小孩的运动轨迹 X=5*cos(t);Y=5*sin(t); Xs=-5*sin(t);Ys=5*cos(t); 程序3：主程序 %main1.m y0=[10,0]; [t,y]=ode45(toolandchild,[0 100],y0); clf;hold on; axis([-6 10 -6 10]); axis(square); plot(y(:,1),y(:,2)); t=0:0.05:6.3; [X,Xs,Y,Ys]=child(t); plot(X,Y,:); hold off; 程序