06第6章 微分方程建模习题解答.docVIP

第 6章 微分方程建模习题解答 6.1 解 设导弹运行曲线的参数方程为 ( ( ( x x(t), y ( y(t), ( 即在时刻t ，导弹的位置在点 (x(t), y(t)) ，这时乙舰的位置在(1, 0t) v 。由于导弹始终对 准乙舰，而导弹运行方向是沿曲线的切线方向，所以有 dy v t ( y ( 0 ， dx 1( x 整理得 dy v0 ( y ( (1( ) ， t x dx 两边对 x 求导，得 dt dy d y dy 2 v0 ( ( (1( x) ( ， dx dx dx dx 2 即 dt d y 2 v0 ( (1( ) . （6.1） x dx dx 2 已知导弹是速度为 5v ，即 0 ( ( ( 2 2 dx ( ( ( ( dy ( ( ( ( dt ( ( ( dt 5v 0 ， dx 由于 ( 0 dt ，所以 2 dx ( ( ( dy ( 1( ( 5v ，即 0 dt ( ( dx dt dx ( 1 5v 0 1( ( ( ( dy dx 2 ( ( ( ， （6.2） 代入式（6.1），得到运动曲线满足的微分方程为 ( 2 ( dy ( ( ( 1 ( ( ( d2 y ( ( dx ( ( , 0 ( x ( ( dx 5(1( x) 1, 2 ( y(0) ( 0, y(0) ( 0. ( 可以求出上述微分方程的解析解，并求得当 x (1时， y ( 0.2083。导弹运行的轨 迹见图 6.1。 -35-(5 (-1)1/5 (-(-1)3/5)1/2)/8 -...- (5 (1/(-1)2/5) (x - 1)6/5)/(12 (-(-1)3/5)1/2) 0.2 0.15 0.1 0.05 0 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 x 图 6.1 导弹运行的轨迹图 符号求解的 Matlab 程序如下 clc, clear y=dsolve(D2y=sqrt(1+(Dy)^2)/5/(1-x),y(0)=0,Dy(0)=0,x) ezplot(y(1),[0,0.9999]) %符号求解时，得到两个分支，这里画出第一个分支 yy=subs(y(1),x,1) %求击中时乙舰行驶的距离 也可以利用 Matlab 求数值解，为了利用 Matlab 软件求数值解，需要做变量替换， 把上述二阶非线性常微分方程转化为一阶常微分方程组的初值问题，令 y ( y 1 ， dy y ( 2 ，则得到一阶微分方程组 dx (dy 1 ( y , ( 2 dx ( ( ( 2 dy 1 y ( 2 ( , 2 dx 5(1( x) ( ( y (0) ( 0, y (0) ( 0. ( ( 1 2 求解的 Matlab 程序如下 clc, clear dyy=@(x,yy)[yy(2); sqrt(1+yy(2)^2)/5/(1-x)]; %定义微分方程右端项的匿名函数 yy0=[0,0]; %初值条件 [x,yy]=ode45(dyy,[0,1-eps],yy0) %为避免奇异点 x=1,右端点取为 1-eps plot(x,yy(:,1)) %画出轨迹曲线 yys=yy(end,1) %求击中时乙舰行驶的距离 数值解的结果和符号解的结果一致。 6.2 解 设法定行车速度为 v ，交叉路口的宽度为 a ，典型的车身长度为b 。考虑到 0 -36-a ( b 车通过路口实际上指的是车的尾部必须通过路口，因此，通过路口的时间为 。 v 0 现在计算刹车距离。设 w 为汽车重量， ( 为摩擦系数，显然，地面对汽车的摩擦 力为 (w ，其方向与运动方向相反。汽车在停车过程中，行驶的距离 x 与时间t 的关系 可由下面的微分方程 2 x w ( d ( (( w ， （6.3） g dt 2 求得，其中 g 是重力加速度。 给出方程（6.3）的初始条件 dx x 0 ( 0 ， ( v ， （6.4） t( dt 0 t(0 于是，刹车距离就是直到速度v ( 0时汽车驶过的距离。 首先，求解二阶微分方程（6.3），对（6.3）式从 0 到t 积分，再利用初始条件 （6.4），得到 dx ( (( ( ， （6.5） gt v 0 dt dx 在 的条件下对（6.5）式从 0 到t 积分，得 ( v 0 dt t(0 1 ( x ( ( (gt2 v t . （6.6） 0 2 dx 在式（6.5）式中令 ( 0 ，得到刹车所用的时间 dt v t0 ( ， 0 (g 从而得到 v 2 x 0 ) ( 。 2(g (t 0 计算黄灯状态的时间 x(t0 ) ( a ( b y ( (T ， v 0 其中T