发布时间 : 星期三 文章2019年整理数学建模微分方程的应用举例55830更新完毕开始阅读6144e77177eeaeaad1f34693daef5ef7bb0d120a
第八节 数学建模——微分方程的应用举例
微分方程在物理学、力学、经济学和管理科学等实际问题中具有广泛的应用，本节我们将集中讨论微分方程的实际应用，读者可从中感受到应用数学建模的理论和方法解决实际问题的魅力.
内容分布
★衰变问题 ★逻辑斯谛方程 ★价格调整问题 ★人才分配问题模型 ★追迹问题
内容要点：
一、衰变问题
镭、铀等放射性元素因不断放射出各种射线而逐渐减少其质量, 这种现象称为放射性物质的衰变. 根据实验得知, 衰变速度与现存物质的质量成正比, 求放射性元素在时刻t的质量.
用x表示该放射性物质在时刻t的质量, 则速度与现存的质量成正比”可表示为
dx表示x在时刻t的衰变速度, 于是“衰变dtdx??kx. (8.1) dt这是一个以x为未知函数的一阶方程, 它就是放射性元素衰变的数学模型, 其中k?0是比例常数, 称为衰变常数, 因元素的不同而异. 方程右端的负号表示当时间t增加时, 质量x减少.
解方程(8.1)得通解x?Ce?kt.若已知当t?t0时, x?x0,代入通解x?Ce?kt中可得
C?x0ekt0, 则可得到方程(8.1)特解
x?x0e?k(t?t0),
它反映了某种放射性元素衰变的规律.
注: 物理学中, 我们称放射性物质从最初的质量到衰变为该质量自身的一半所花费的时间为半衰期, 不同物质的半衰期差别极大. 如铀的普通同位素(通常的镭(
226238U)的半衰期约为50亿年;
Ra)的半衰期是1600年.半衰期是上述放射性物质的特征, 然而半衰期却不依
226赖于该物质的初始量, 一克
Ra衰变成半克所需要的时间与一吨226Ra衰变成半吨所需要
的时间同样都是1600年, 正是这种事实才构成了确定考古发现日期时使用的著名的碳-14测验的基础.
二、 逻辑斯谛（Logistic）方程:
逻辑斯谛方程是一种在许多领域有着广泛应用的数学模型, 下面我们借助树的增长来建立该模型.
一棵小树刚栽下去的时候长得比较慢, 渐渐地, 小树长高了而且长得越来越快, 几年不见, 绿荫底下已经可乘凉了; 但长到某一高度后, 它的生长速度趋于稳定, 然后再慢慢降下来. 这一现象很具有普遍性. 现在我们来建立这种现象的数学模型.
如果假设树的生长速度与它目前的高度成正比, 则显然不符合两头尤其是后期的生长情形, 因为树不可能越长越快; 但如果假设树的生长速度正比于最大高度与目前高度的差, 则又明显不符合中间一段的生长过程. 折衷一下, 我们假定它的生长速度既与目前的高度,又与最大高度与目前高度之差成正比.
设树生长的最大高度为H(m), 在t(年)时的高度为h(t), 则有
dh(t)?kh(t)[H?h(t)] (8.2) dt其中k?0是比例常数. 这个方程为Logistic方程. 它是可分离变量的一阶常数微分方程.
下面来求解方程(8.2). 分离变量得
dh?kdt,
h(H?h)两边积分
dh?h(H?h)??kdt,
得 或
1[lnh?ln(H?h)]?kt?C1, Hh?ekH?t1CH?C2ekHt,
H?h故所求通解为
C2HekHtHh(t)??,
1?C2ekHt1?Ce?kHt其中的C?C????1?e?C1H?0?是正常数. C2?函数h(t)的图象称为Logistic曲线. 图8-8-1所示的是一条典型的Logistic曲线, 由于它的形状, 一般也称为S曲线. 可以看到, 它基本符合我们描述的树的生长情形. 另外还可以
算得
t???limh(t)?H.
这说明树的生长有一个限制, 因此也称为限制性增长模式（阻滞增长模型）.
注: Logistic的中文音译名是“逻辑斯谛”. “逻辑”在字典中的解释是“客观事物发展的规律性”, 因此许多现象本质上都符合这种S规律. 除了生物种群的繁殖外, 还有信息的传播、新技术的推广、传染病的扩散以及某些商品的销售等. 例如流感的传染、在任其自然发展(例如初期未引起人们注意)的阶段, 可以设想它的速度既正比于得病的人数又正比于未传染到的人数. 开始时患病的人不多因而传染速度较慢; 但随着健康人与患者接触, 受传染
的人越来越多, 传染的速度也越来越快; 最后, 传染速度自然而然地渐渐降低, 因为已经没有多少人可被传染了.
下面举两个例子说明逻辑斯谛的应用.
人口阻滞增长模型 1837年, 荷兰生物学家Verhulst提出一个人口模型
dy?y(k?by),dt其中k,b的称为生命系数.
y(t0)?y0 (8.3)
我们不详细讨论这个模型, 只提应用它预测世界人口数的两个有趣的结果.
有生态学家估计k的自然值是0.029. 利用本世纪60年代世界人口年平均增长率为2%以及1965年人口总数33.4亿这两个数据, 计算得b?2,从而估计得:
(1)世界人口总数将趋于极限107.6亿. (2)到2000年时世界人口总数为59.6亿.
后一个数字很接近2000年时的实际人口数, 世界人口在1999年刚进入60亿. 新产品的推广模型 设有某种新产品要推向市场, t时刻的销量为x(t),由于产品性能良好, 每个产品都是一个宣传品, 因此, t时刻产品销售的增长率虑到产品销售存在一定的市场容量N, 统计表明
dx与x(t)成正比, 同时, 考dtdx与尚未购买该产品的潜在顾客的数量dtN?x(t)也成正比, 于是有
dx?kx(N?x) dt其中k为比例系数. 分离变量积分, 可以解得
(8.4)
x(t)?由
N
1?Ce?kNt (8.5)
dxdtx的图像知，当x(t*)?N时, 则有
dxN?0,即销量x(t)单调增加. 当x(t*)?时, dt2NNd2xd2xd2x**?0;当x(t)?时, 2?0;当x(t)?时, 2?0;即当销量达到最大需求量222dtdtdtN的一半时, 产品最为畅销, 当销量不足N一半时, 销售速度不断增大, 当销量超过一半时,
销售速度逐渐减少.
国内外许多经济学家调查表明. 许多产品的销售曲线与公式(8.5)的曲线(逻辑斯谛曲线)十分接近. 根据对曲线性状的分析, 许多分析家认为, 在新产品推出的初期, 应采用小批量生产并加强广告宣传, 而在产品用户达到20%到80%期间, 产品应大批量生产; 在产品用户超过80%时, 应适时转产, 可以达到最大的经济效益.
三、价格调整模型
假设某种商品的价格变化主要服从市场供求关系. 一般情况下,商品供给量S是价格P的单调递增函数, 商品需求量Q是价格P的单调递减函数, 为简单起见, 分别设该商品的供
给函数与需求函数分别为
S(P)?a?bP,Q(P)????P (8.6)
其中a,b,?,?均为常数, 且b?0,??0.
当供给量与需求量相等时, 由(8.6)可得供求平衡时的价格
Pe?并称Pe为均衡价格.
??a ??b一般地说, 当某种商品供不应求, 即S?Q时, 该商品价格要涨, 当供大于求, 即
S?Q时, 该商品价格要落. 因此, 假设t时刻的价格P(t)的变化率与超额需求量Q?S成
正比, 于是有方程
dP?k[Q(P)?S(P)] dt其中k?0,用来反映价格的调整速度.
将(8.6)代入方程, 可得
dP??(Pe?P) (8.7) dt其中常数??(b??)k?0,方程(8.7)的通解为
??t P(t)?Pe?Ce假设初始价格P(0)?P0,代入上式, 得C?P0?Pe,于是上述价格调整模型的解为
??t P(t)?P?(P?P)ee0e由于??0知, t???时, P(t)?Pe.说明随着时间不断推延, 实际价格P(t)将逐渐趋近均衡价格Pe.
四、人才分配问题模型
每年大学毕业生中都要有一定比例的人员留在学校充实教师队伍, 其余人员将分配到国民经济其他部门从事经济和管理工作. 设t年教师人数为x1(t),科学技术和管理人员数目为x2(t),又设1个教员每年平均培养?个毕业生, 每年从教育、科技和经济管理岗位退休、死亡或调出人员的比率为?(0???1),?表示每年大学生毕业生中从事教师职业所占比率