微分方程及几个简单实例.ppt

微分方程模 型 浙江大学数学建模实践基地 §3.1 微分方程的几个简单实例 * * 在许多实际问题中，当直接导出变量之间的函数关系较为困难，但导出包含未知函数的导数或微分的关系式较为容易时，可用建立微分方程模型的方法来研究该问题， 本节将通过一些最简单的实例来说明微分方程建模的一般方法。在连续变量问题的研究中，微分方程是十分常用的数学工具之一。 例1 （理想单摆运动）建立理想单摆运动满足的微分方程，并得出理想单摆运动的周期公式。 从图3-1中不难看出，小球所受的合力为mgsinθ，根据牛顿第二定律可得： 从而得出两阶微分方程： （3.1） 这是理想单摆应满足的运动方程 （3.1）是一个两阶非线性方程，不易求解。当θ很小时，sinθ≈θ，此时，可考察（3.1）的近似线性方程： （3.2） 由此即可得出 （3.2）的解为: θ(t)= θ0cosωt 其中 当 时,θ(t)=0 故有 M Q P mg 图3-1 （3.1）的近似方程 例2 我方巡逻艇发现敌方潜水艇。与此同时敌方潜水艇也发现了我方巡逻艇，并迅速下潜逃逸。设两艇间距离为60哩，潜水艇最大航速为30节而巡逻艇最大航速为60节，问巡逻艇应如何追赶潜水艇。 这一问题属于对策问题，较为复杂。讨论以下简单情形： 敌潜艇发现自己目标已暴露后，立即下潜，并沿着直 线方向全速逃逸，逃逸方向我方不知。 设巡逻艇在A处发现位于B处的潜水艇，取极坐标，以B为极点，BA为极轴，设巡逻艇追赶路径在此极坐标下的方程为r=r(θ)，见图3-2。 B A A1 dr ds dθ θ 图3-2 由题意， ，故ds=2dr 图3-2可看出， 故有: 即: （3.3） 解为： （3.4） 先使自己到极点的距离等于潜艇到极点的距离然后按（3.4）对数螺线航行，即可追上潜艇。 追赶方法如下： 例3 一个半径为Rcm的半球形容器内开始时盛满了水，但由于其底部一个面积为Scm2的小孔在t=0时刻被打开，水被不断放出。问：容器中的水被放完总共需要多少时间？ 解: 以容器的底部O点为 原点，取坐标系如图3.3所示。令h(t)为t时刻容器中水的高度，现建立h(t)满足的微分方程。 设水从小孔流出的速度为v(t)，由力学定律，在不计水的内部磨擦力和表面张力的假定下，有： 因体积守衡，又可得： 易见： 故有： 即： 这是可分离变量的一阶微分方程，得 R x y S O 图3-3 h r 例4 一根长度为l的金属杆被水平地夹在两端垂直的支架上，一端的温度恒为T1，另一端温度恒为T2，（T1、T2为常数，T1 T2）。金属杆横截面积为A，截面的边界长度为B，它完全暴露在空气中，空气温度为T3，（T3 T2，T3为常数），导热系数为α，试求金属杆上的温度分布T(x)，（设金属杆的导热率为λ） 一般情况下，在同一截面上的各点处温度也不尽相同，如果这样来考虑问题，本题要建的数学模型当为一偏微分方程。 但由题意可以看出，因金属杆较细且金属杆导热系数又较大，为简便起见，不考虑这方面的差异，而建模求单变量函数T(x)。 热传导现象机理：当温差在一定范围内时，单位时间里由温度高的一侧向温度低的一侧通过单位面积的热量与两侧的温差成正比，比例系数与介质有关。 T1 T2 o x A B T3 l dt时间内通过距离O点x处截面的热量为： dt时间内通过距离O点x+dx处截面的热量为： 由泰勒公式： 金属杆的微元[x,x+dx]在dt内由获得热量为： 同时，微元向空气散发出的热量为： 系统处于热平衡状态，故有： 所以金属杆各处温度T(x)满足的微分方程: 这是一个两阶常系数线性方程，很容易求解 *