科学技术当中的微分方程模型.ppt

问题分析与求解： 因为圆桶下沉满足牛顿第二定律：运动物体受到的力等于该物体的质量与其运动加速度的乘积，即F=ma。在本问题中有 ，F=W-B-D=W-B-kv，y(0)=0，v(0)=0， 于是可以得到如下微分方程： (1) 重力 浮力 阻力 ①　求解问题（１），Matlab命令为 用Matlab软件求解。 求解出函数y(t)后，求出圆桶下落到水深为91米海底的时间，即求满足y(t)=91的时间t1,然后求出在时间t1的速度即可以得到问题的解答。这里的速度关系由直接求函数y(t)对t的导数得到。 syms t; syms y; y=dsolve(m*D2y+k*Dy-w+b=0,y(0)=0,Dy(0)=0) 运行结果为 y =-m/k^2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k^2 即 运行结果如图 b=1025.94*0.208; k=0.12; w=239.456; m=w/9.8; % 这里重力加速度取为9.8米/秒2 t=0?:0.1?:20?; y=-m/k^2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k^2-91?; plot(t,y) ②　求解方程y(t)=91，先画出图形，命令为 可见方程的根在t=13附近， 为求解方程，建立M文件： function y=my0(t) b=1025.94*0.208; k=0.12; w=239.456; m=w/9.8; y=-m/k^2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k^2-91?; 结果为t =13.204223061992241 t=fsolve(‘my0, 13) 用如下的命令求解： ③　求y(x)的导数，Matlab代码为 b=1025.94*0.208; k=0.12; w=239.456; m=w/9.8; syms t; y1=diff(-m/k^2*exp(-k/m*t)*(-w+b)-(-w+b)/k*t+m*(-w+b)/k^2-91)? 运行结果为 y1 = -3190720108293556007631/14692224126156800000*exp(- 21/4276*t)+81439/375 即 　计算结果表明圆桶下沉到水深约为91米海底时的速度约为13.6361米/秒，这个速度大于12.2米/秒的速度，因此，计算结果说明美国原子能委员会以往处理浓缩的放射性废料方法是不安全的。 故t =13.204223061992241时的速度为v =13.636101129315506 问题： 案例四　 弱肉强食微分方程模型 　　生活在同一环境中的各类生物之间，进行着残酷的生存竞争。设想一海岛，居住着狐狸与野兔，狐吃兔，兔吃草，青草如此之丰富，兔子们无无食之忧，于是大量繁殖；兔子一多，狐易得食，狐量亦增，而由于狐狸数量增加吃掉大量兔子，狐群又进入饥饿状态而使其总数下降，这时兔子相对安全，于是兔子总数回升。就这样，狐兔数目交替地增减，无休止的循环，遂形成生态的动态平衡。如何用建立数学模型描述并预测下一阶段情况？ 变量说明： 模型各项意义 其中a, b, c, d 均为正常数。 （１） x(t)―――t时刻兔子数目 y(t)―――t时刻狐狸数目 数学模型： ax表示兔子的繁殖速度与现存兔子数成正比； -bxy表示狐兔相遇，兔子被吃掉的速度； -cy表示狐狸因同类争食造成的死亡速度与狐狸总数成正 比； dxy表示狐兔相遇，对狐狸有好处而使狐狸繁殖增加的速 度。 这是意大利著名物理学著名生物数学家Volterra建立的微分方程模型，这是达尔文主义的数学表达，称为Volterra-Lotka模型。 （２） 设该问题的初始条件为 令 模型求解与分析: 得系统的两个稳定点（奇点）为 和 ①　取x(t)=0，求解问题（１）与（２），得其解为 　取y(t)=0，求解问题（１）与（２），得其解为 说明x轴和y轴都是方程组（１）的轨线，它们的实际意义是：如果没有兔子，狐狸因为没有食物将趋于死亡；而如果没有狐狸，兔子将无限制增长。 ② 当x , y 0 时，方程组（１）等价于 这是一个可分离变量的方程，积分得通解 两边取指数得 （３） 令 ③　证明：当初值 时，方程（３）确定了xoy面上位于第一象限的一族封闭曲线。 ，由于 所以当 时，方程无 x 0 , y