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

①　求解问题（１），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，求解问题（１）与（２），得其解为 * 西安理工大学应用数学系 西安理工大学应用数学系 * 实际问题当中的微分方程模型 　　对于现实世界的变化，人们关注的往往是变量之间的变化率，或变化速度、加速度以及所处的位置随时间的发展规律，之中的规律一般可以写成一个（偏）微分方程或方程组。所以实际问题中，有大批的问题可以用微分方程来建立数学模型，涉及的领域包括物理学、化学、天文学、生物学、力学、政治、经济、军事、人口、资源等等（有传统领域，有新的非传统领域）。本专题首先归纳出微分方程模型所涉及的领域以及常见模型类型，然后列举几个领域的典型例子。 一、不同领域中的微分方程模型 1 社会及市场经济中的微分方程模型 （1）综合国力的微分方程模型 （2）诱发投资与加速发展的微分方程模型 （3）经济调整的微分方程模型 （4）广告的微分方程模型 （5）价格的微分方程模型 2 战争中的微分方程模型 （1）军备竞赛的