内容发布更新时间 : 2024/12/27 12:05:47星期一 下面是文章的全部内容请认真阅读。

第五章 微分方程建模案例

微分方程作为数学科学的中心学科，已经有三百多年的发展历史，其解法和理论已日臻完善，可以为分析和求得方程的解（或数值解）提供足够的方法，使得微分方程模型具有极大的普遍性、有效性和非常丰富的数学内涵。微分方程建模包括常微分方程建模、偏微分方程建模、差分方程建模及其各种类型的方程组建模。微分方程建模对于许多实际问题的解决是一种极有效的数学手段，对于现实世界的变化，人们关注的往往是其变化速度、加速度以及所处位置随时间的发展规律，其规律一般可以用微分方程或方程组表示，微分方程建模适用的领域比较广，涉及到生活中的诸多行业，其中的连续模型适用于常微分方程和偏微分方程及其方程组建模，离散模型适用于差分方程及其方程组建模。本章主要介绍几个简单的用微分方程建立的模型，让读者一窥方程的应用。下面简要介绍利用方程知识建立数学模型的几种方法：

1．利用题目本身给出的或隐含的等量关系建立微分方程模型

这就需要我们仔细分析题目，明确题意，找出其中的等量关系，建立数学模型。

例如在光学里面，旋转抛物面能将放在焦点处的光源经镜面反射后成为平行光线，为了证明具有这一性质的曲线只有抛物线，我们就是利用了题目中隐含的条件——入射角等于反射角来建立微分方程模型的。

2．从一些已知的基本定律或基本公式出发建立微分方程模型

我们要熟悉一些常用的基本定律、基本公式。例如从几何观点看，曲线

y?y(x)上某点的切线斜率即函数y?y(x)在该点的导数；力学中的牛顿第二运

动定律：F?ma，其中加速度a就是位移对时间的二阶导数，也是速度对时间的一阶导数等等。从这些知识出发我们可以建立相应的微分方程模型。 例如在动力学中，如何保证高空跳伞者的安全问题。对于高空下落的物体，我们可以利用牛顿第二运动定律建立其微分方程模型，设物体质量为m，空气阻

209

力系数为k，在速度不太大的情况下，空气阻力近似与速度的平方成正比；设时刻t时物体的下落速度为v，初始条件：v(0)?0. 由牛顿第二运动定律建立其微分方程模型：

mdv?mg?kv2 dt求解模型可得：

mg(exp[2tv?kg]?1)m kgk(exp[2t]?1)m由上式可知，当t???时，物体具有极限速度：

v1?limv?t??mg， k其中，阻力系数k???s，?为与物体形状有关的常数，?为介质密度，s为物体在地面上的投影面积。根据极限速度求解式子，在m,?,?一定时，要求落地速度v1不是很大时，我们可以确定出s来，从而设计出保证跳伞者安全的降落伞的直径大小来。

3．利用导数的定义建立微分方程模型 导数是微积分中的一个重要概念，其定义为

f?(x)?lim?x?0f(x??x)?f(x)?y?lim， ?x?0?x?x商式

?y表示单位自变量的改变量对应的函数改变量，就是函数的瞬时平均变化?x率，因而其极限值就是函数的变化率。函数在某点的导数，就是函数在该点的变化率。由于一切事物都在不停地发展变化，变化就必然有变化率，也就是变化率是普遍存在的，因而导数也是普遍存在的。这就很容易将导数与实际联系起来，建立描述研究对象变化规律的微分方程模型。

例如在考古学中，为了测定某种文物的绝对年龄，我们可以考察其中的放射性物质（如镭、铀等），已经证明其裂变速度（单位时间裂变的质量，即其变化率）与其存余量成正比。我们假设时刻t时该放射性物质的存余量R是t的函数，

210

由裂变规律，我们可以建立微分方程模型:

dR??kR dt期中k是一正的比例常数，与放射性物质本身有关。求解该模型，我们解得：

R?Ce?kt，其中c是由初始条件确定的常数。从这个关系式出发，我们就可以测

定某文物的绝对年龄。（参考碳定年代法）

另外，在经济学领域中，导数概念有着广泛的应用，将各种函数的导函数（即函数变化率）称为该函数的边际函数，从而得到经济学中的边际分析理论。

4．利用微元法建立微分方程模型

一般的，如果某一实际问题中所求的变量p符合下列条件：p是与一个变量

t的变化区间[a,b]有关的量；p对于区间[a,b]具有可加性；部分量?pi的近似值可表示为f(?i)?ti。那么就可以考虑利用微元法来建立微分方程模型，其步骤是：首先根据问题的具体情况，选取一个变量例如t为自变量，并确定其变化区间

[a,b]；在区间[a,b]中随便选取一个任意小的区间并记作[t,t?dt]，求出相应于

这个区间的部分量?p的近似值。如果?p能近似的标示为[a,b]上的一个连续函数在t处的值f(t)与dt的乘积，我们就把f(t)dt称为量p的微元且记作dp.这样，我们就可以建立起该问题的微分方程模型：

dp?f(t)dt.

对于比较简单的模型，两边积分就可以求解该模型。

例如在几何上求曲线的弧长、平面图形的面积、旋转曲面的面积、旋转体体积、空间立体体积；代数方面求近似值以及流体混合问题；物理上求变力做功、压力、平均值、静力矩与重心；这些问题都可以先建立他们的微分方程模型，然后求解其模型。

5．熟悉一些经典的微分方程模型，对一些类似的问题，经过稍加改进或直接套用这些模型。

多年来，在各种领域里，人们已经建立起了一些经典的微分方程模型，熟悉这些模型对我们是大有裨益的。

211