2019年微分方程详解.doc

第二章 微分方程 本章学习目的： 本章的主要目的在于：学习微分方程模型的建立、求解方法、分析结果及解决实际问题的全过程。 1．知道求解微分方程的解析法、数值解法以及图形表示解的方法； 2．理解求微分方程数值解的欧拉方法，了解龙格——库塔方法的思想； 3．熟练掌握使用MATLAB软件的函数求微分方程的解析解、数值解和图形解； 4．通过范例学习怎样建立微分方程模型和分析问题的思想。 §2.1 引例 在《大学物理》中，我们曾学习过简谐振动（如：弹簧振子、单摆）2xd2?x?0?，那就是一个典型的二阶常微分方程的模型。 2dt这里我们讨论“倒葫芦形状容 器壁上的刻度问题”。 1 x 对于圆柱形状容器壁上的容积刻度，可以利用圆柱体体积公式：2?与相对于容器底部的任为常数，体积V，其中容器的直径D4V?HD/ H成正比，因此在容器壁上可以方便地标出容积刻度。意高度而对于几何形状不规则的容器，比如“倒葫芦形状”的容器壁上如何 标出容积刻度呢？ 如图所示，建立坐标系，由微元法分析可知：12?dx(dV?x)DD表示高度，直径是高度的函数，记为，其中x 41dV2?)?xD( x）。可得微分方程：（4dx0)?V(0的部分测试数无解析表达式，只给出D(x)D(x)如果该方程中的函数 据，如何求解此微分方程呢？h=0.2; d=[0.04,0.11,0.26,0.56,1.04,1.17]; 2 x(1)=0;v(1)=0; for k=1:5 x(k+1)=x(k)+h; v(k+1)=v(k)+(h/2)*(pi/4)*(d(k)^2+d(k+1)^2); end x=x(1:6),v=v(1:6), plot(x,v) x = Columns 1 through 5 0 0.2000 0.4000 0.6000 0.8000 Column 6 1.0000 v = Columns 1 through 5 0 0.0011 0.0073 0.0373 0.1469 Column 6 0.3393 3 0.350.30.20.0.10.0.00.0.0.0.0.0.0.0.0. 微分方程模型的建立2.2 §等关键词提”、“减少”、“增加”“在工程实际问题中，“改变”、变化”“边际的、、“衰变”、示我们注意什么量在变化，关键词“速率”“增长” 等常涉及到导数。dyv? 就是一个简单的一阶微分方程。我们熟悉的速度公式： dt 微分方程是指含有导数或微分的等式。 一般形式：)(n?0,y)?yF(x,,y,? )?(n(n)1?).x或：yf?(,,y,y,y?4 常用的建立微分方程的方法有：运用已知物理定律；利用平衡与增长式；运用微元法；应用分析法。 2.2.1 运用已知物理定律 建立微分方程模型时应用已知物理定律，可事半功倍。 0C18的房间内，该物体最初的例2.1 一个较热的物体置于室温为000C 需要多C。想知道它的温度降到30C，3分钟以后降到5060温度是少时间？10分钟以后它的温度是多少？ 牛顿冷却（加热）定律：将温度为T的物体放入处于常温 m 的介质中时，T的变化速率正比于T与周围介质的温度差。 分析：假设房间足够大，放入温度较低或较高的物体时，室内温度基本不受影响，即室温分布均衡，保持为m，采用牛顿冷却定律是一个相当好的近似。 建立模型：设物体在冷却过程中的温度为T(t)，t≥0， dT与（T?m）成正比，建立微分方根据牛顿加热（冷却）定律： dtdT????k(T?m) 程 （2.1） ?dt?T(0)?60?，m=18。其中参数k >0 2.2.2 利用平衡与增长式 许多研究对象在数量上常常表现出某种不变的特性，如封闭区域内的能量、货币量等。利用变量间的平衡与增长特性，可分析和建立有关变量间的相互关系。 5 此类建模方法的关键是分析并正确描述基本模型的右端，使平衡式成立。 例2.2 战斗模型：两方军队交战，希望为这场战斗建立一个数学模型，应用这个模型达到如下目的： 1. 预测哪一方将获胜？ 2. 估计获胜的一方最后剩下多少士兵？ 3. 计算失败的一方开始时必须投入多少士兵才能赢得这场战斗？ 解：模型建立： 设 x(t)： t 时刻X方存活的士兵数 y(t)： t 时刻Y方存活的士兵数 假设： 1）双方所有士兵不是战死就是活着参加战斗，x(t)与y(t)都是连续变量； 2）Y方军队