微分方程建模学习
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
微分方程建模
一般说来,微分方程建模的方法大致可以分为以下的几个步骤:
1•根据实际问题的要求确定要研究的量,包括自变量、未知函数、必要的参数等以及它 们各自的变化区间; 2•列方程。
可以在合理假设的前提下,禾U 用导数表示斜率、速度、变化率的实际意义, 根据一些基本定理(几何的、物理的、化学的或生物学的等等)或规律,找出未知函数的导 数(或微分)与相关各量之间的等量关系式,建立微分方程并确定定解条件(注:如果没有 现成的定理可供利用,也可以用微元分析法与模拟近似法列出微分方程)
;
3.解微分方程;
4•对模型的适用性作出评价,即用已知的数据检验微分方程的解是否与实际相符。
若结 果与实际存在一定的差距,则还要对方程进行修正和调整,直到得出较满意的结果为止。
下面,我们就通过一些实例说明微分方程建模的具体步骤。
一. 增长模型
在自然界和社会的经济活动中,许多量的变化都遵循着一个基本的规律: 任一单位时间
的增量都与该量自身当时的大小成正比。
运用这一基本规律,就可以建立起各种各样的增长
模型。
1. 马尔萨斯人口模型
严格地讲,讨论人口问题所建立的模型应属于离散型模型。
但在人口基数很大的情况下,
突然增加或减少的只是单一的个体或少数几个个体,
相对于全体数量而言, 这种改变量是极
其微小的,因此,我们可以近似地假设人口随时间连续变化甚至是可微的。
这样,我们就可
以采用微分方程的工具来研究这一问题。
最早研究人口问题的是英国的经济系家马尔萨斯 (Malthus ) ( 1766—1834)。
他根据百余
年的人口资料,经过潜心研究,在
1798年发表的《人口论》中首先提出了人口增长模型。
他的基本假设是:任一单位时刻人口的增长量与当时的人口总数成正比, 且比例系数为常数。
于是,设t 时刻的人口总数为 y (t ),则单位时间人口的增长量即为
y(t t) y(t)
t
令t 0 ,可得微分方程
根据基本假设,有
y(t t) y(t)
r y(t)
(r 为比例系数)
dy dt
这就是著名的马尔萨斯人口方程。
若假设t t o 时的人口总数为y o ,则不难求得该方程的特
即任一时刻的人口总数都遵循指数规律向上增长。
人们曾用这个公式对 六十余年世界的人口资料进行了检验,发现计算结果与人口的实际情况竟然是惊人的吻合!
2. 放射性元素衰变模型
放射性元素的质量随时间的推移而逐渐减少(负增长) ,这种现象称为衰变。
由物理学
定律知,放射性元素任一时刻的衰变速度与该时刻放射性元素的质量成正比。
根据这一原理,
我们也可以通过微分方程研究放射性元素衰变的规律。
设放射性元素t 时刻的质量m m(t),则其衰变速度就是 也,于是可得
dt
dm
m
(4.3)
dt
其中 0是比例常数,可由该元素的半衰期(质量蜕变到一半所需要的时间)确定; 前
置负号表明放射性元素的质量 m 是随时刻t 递减的。
如果在初始时刻(t 0)放射性元素的质量 m m 0,则可求得该方程的特解为
m(t) m ° e t
(4.4)
这说明放射性元素的质量也是随时刻按指数规律递减的。
为了能将求得的放射性元素衰变规律应用于实际, 还必须确定上式中的比例常数 。
这
时,我们可以假设放射性元素的半衰期为
T ,从而有
m °
T
m ° e 2
ln 2
解之,得
——,于是反映放射性元素衰变规律的( 4.4)式又可以表示为
T
m(t) m ° e T
(4.5)
并由此可解得
丄
T , m °
t In (4.6)
In 2 m(t)
它所反映的是放射性元素由初始时刻的质量 m 0衰减到m(t)所需要的时间。
放射性元素的衰变规律常被考古、 地质方面的专家用于测定文物和地质的年代, 其中最
(4.1)
解为
y o
e r(t
t o )
(4.2)
1700 — 1961 达二百
常用的是14C (碳-12的同位素)测定法。
这种方法的原理是:大气层在宇宙射线不断的轰
击下所产生的中子与氮气作用生成了具有放射性的 14
C ,并进一步氧化为二氧化碳,二氧化
碳被植物所吸收,而动物又以植物为食物, 于是放射性碳就被带到了各种动植物的体。
对于
具有放射性的
14
C 来说,不论是存在于空气中还是生物体,它都在不断地蜕变。
由于活着的
生物通过新代不断地摄取
14
C ,因而使得生物体的14C 与空气中的14
C 有相同的百分含量;
旦生物死亡之后,随着新代的停止,尸体的14
C 就会不断地蜕变而逐渐减少, 因此根据
14
C 蜕变减
少量的变化情况并利用(4.6 )式,就可以判定生物死亡的时间。
F 面,我们就来看一个运用 14
C 测定法确定年代的具体实例:
1972年8月,出土了马王堆一号墓(注:出土时因墓中女尸历经千年而未腐曾经轰动 世界)。
经测定,出土的木炭标本中 14
C 的平均原子蜕变速度为
29.78次/分,而新砍伐烧成
的木炭中
14
C 的平均原子蜕变速度为
38.37次/分;如果
14
C 的半衰期取为5568年(注:14
C
的半衰期在各种资料中说法不一,分别有 5568年、5580年和5730年不等),那么,怎样才
能根据以上数据确定这座墓葬的大致年代呢?
(4.6)式作进一步的修改:
从而有
上面两式相除,得
在确定衰变时间的公式(4.6)中,由于
m 0和m (t )表示的分别是该墓下葬时和出土时
木炭标本中
14
C 的含量,而测量到的是标本中 14
C 的平均原子蜕变速度,所有我们还要对
对(4.4)式求导,得
m (t)
t
m °e
m(t)
m(0)
m(0)
m o
代入(4.6),得
曰 是
m m(t)
(0) m 0 m(t)
t 丄a
In 2 m(t)
(4.7)
衰变时间由(4.6)式根据14C 含量的变化情况确定就转化为由(4.7)式根据14C 衰变速 度的变化情况来确定,这就给实际操作带来了很大的方便。
在本例中,T 5568年,m (t ) 29.78次/分,m (0)虽然表示的是下葬时所烧制的木 炭中
14
C 的衰变速度,但考虑到宇宙射线的强度在数千年的变化不会很大,
因而可以假设现
代生物体中
14
C 的衰变速度与马王堆墓葬时代生物体中
14
C 的衰变速度相同,即可以用新砍
伐烧成的木炭中
14
C 的平均原子蜕变速度 38.37次/分替代m (0)。
代入(4.7)可求得
t 工 ln^ 逢"酒 2036 (年) In 2 m (t) In 2 29.78
若以T 5580年或T 5730年计算,则可分别算得 t 2040年或t 2095年,即马王堆 一号墓大约是2000多年前我国汉代的墓葬。
(注:后经进一步考证,确定墓主人为汉代国
丞相利仓的夫人,名辛追。
)
3. 固定资产折旧模型
企业在进行成本核算的时候, 经常要计算固定资产的折旧。
一般说来,固定资产在任 单位时刻的折旧额与当时固定资产的价值都是成正比的。
函数关系。
假定某固定资产五年前购买时的价格是 试估算再过10年后的价值。
固定资产单位时刻的折旧额可表示为 一
t)
P(t)
,由题意可得
t
不难求得该方程的通解为
为方便计算,记五年前的时刻为
t 0,于是有初始条件
P(0) 10000
k t
P(t) 10000 e
试研究固定资产价值 P 与时间t 的
10000元,而现在的价值为 6000元,
首先我们可以假设t 时刻该固定资产的价值为
P P(t),则在[t, t
t ]这段时间该
P(t t) P(t)
t
kP(t)
(k 0)
dP kP
dt
P(t) Ce
kt
代入通解,可求得C
10000,故原方程的特解为
为确定比例常数
k ,可将另一个条件 P(5)
6000代入上式,
5 k
6000
10000e
从而有
P(t) 10000e 5 3
10000(f) 这就是价值P 与时间
t 之间的函数关系。
于是,再过
10 年(即 t 15)
该固定资产的价
值即为
P (15) 10000 e|)
3
2160 (元)
二.阻滞增长模型
与以上所讨论的增长模型不同,实际中存在着大量的另一类增长模型:
1.弗尔哈斯特人口模型
在人口模型的研究中,马尔萨斯得出了“任一时刻的人口总数都遵循指数规律向上增长” 的结论,并得到了一段时期人口数据的验证,然而,随着人口基数的增大,公式( 4.2)所
暴露的不足之处也越来越明显了。
根据公式(4.2)我们不难计算出,世界人口大约35年就要翻一番。
事实上,设某时刻
的世界人口数为y o,人口增长率为2%,且经过T年就要翻一番,则有
0.02T 0.02T
2y o y°e 即 e 2
解之,即得T 50ln 2 34.6 (年)。
于是,我们以1965年的世界人口33.4亿为基数进行计算,可以得到如下的一系列人口数据:
2515年200万亿
2625年1800万亿
2660年3600万亿
若按人均地球表面积(包括水面、船上)计算,2625年仅为0.09平方米/人,也就是说必须人挨着人站着才能挤得下;而35年后的2660年,人口又翻了一番,那就要人的肩上再
站人了。
而且随着时间的推移,我们有
lim y0e r(tt0)
这显然不符合人口发展的实际。
这说明,在人口基数不是很大的时候,马尔萨斯人口方程还
能比较精确地反映人口增长的实际情况,但当人口数量变得很大时,其精确程度就大大降低
了。
究其根源,是随着人口的迅速膨胀,资源短缺、环境恶化等问题越来越突出,这些都将限制人口的增长。
如果考虑到这些因素,就必须对马尔萨斯人口方程进行修改。
1845年,荷兰的数学、生物学家弗尔哈斯特(Verhulst)提出了一个修改方案,即将方
程修改为
dy
r y by2(0 b r)
dt
其中r、b称为“生命系数”。
由于b r,因此当y不太大时,by2这一项相对于r y可
以忽略不计;而当y很大时,by2这一项所起的作用就不容忽视了,它降低了人口的增长
速度。
慢增长期”。
在利用(4.9)式对人口的发展情况进行预测之前,还必须确定恰当的 以下方法来计算: 由方程(4.8)可得
dy dt
. r by y
于是,我们就有了下面的人口模型:
普 r y by
2
y
tt o
y
(4.8)
这是一个可分离变量的一阶微分方程。
解之,可得
r y 0
y
by。
(r by 。
)e r(t t0)
(4.9)
这就是人口 y 随时间t 的变化规律。
F 面,我们就对(4.9 )作进一步的讨论,并根据它对人
口的发展情况作一些预测。
首先,由于
lim y lim
t
t
r y °
r (t t 0)
by 0 (r by 0) e
即不论人口的基数如何,随着时间的推移,人口总量最终将趋于一个确定的极限值
其次,由理
dt
r y by 2
可得
y r y 2by y (r
2by)
r
,易知这正是函数(4.9)
2b
的图象 (称为 “人口增长曲线”或“
曲线”)拐点的纵坐标,它恰好位于人口总量极限值
半的位置(如图所示)
由于y —时y
2b
0,故d^是递增的,此时
dt
称为人口的“加速增长期”
r
;而当y 时y 0,
2b
故dy 是递减的,此时称为人口的“缓 dt
b 值,它可以按
dy
其中,dy
表示人口的理论增长率,而
旦 则表示人口的实际增长率。
如果我们以 1965年
dt y
率0.02代入上式,有
2.单种群动物模型
所谓单种群动物模型指的是某单一种群动物在一个相对封闭环境中生长衰亡规律的模 型。
这个模型可以运用微元分析法建立和求解。
假设某动物种群t 时刻的数量为 N N(t),出生率与死亡率分别为 m 、n ,且任一 时刻的出生数与死亡数都与当时的动物数量成正比,则在时间间隔 [t,t dt],该动物的出
生量与死亡量就分别为 mNdt 与nNdt ,由于
增量=出生量一死亡量
故有
dN mNdt nNdt (m n) Ndt
即
dN
(m n)N
(4.10)
dt
这也是一个可分离变量的微分方程,显然,它与马尔萨斯人口方程(
4.1)在本质上是完全
相同的。
如果初始时刻该动物种群的数量为 N 0,即N|t 0 N 0,则易求得该方程的特解为
N N 0e
(m n)t
这说明该动物种群的数量也是遵循指数规律增长的,其中
m n 称为“自然增长率”。
从以上的结果中可以看出,如果m n ,则该动物种群的数量随着时间的推移将会无限 制地增加,这显然是不可能的。
之所以产生这样的谬误,是因为我们把
m 、n 当成了常量,
即视出生率与死亡率一成不变,而事实却并非如此。
因为随着动物数量的增加,食品短缺、 环境恶化、瘟疫流行等问题都将导致出生率的降低和死亡率的升高,为此,可令
的人口数3.34 109
为初值,
并把某些生态学家估计的
r 的自然值0.029及人口的实际增长
即可求得 b 2.695 10 12
0.02 0.029 b (3.34
109
)。
于是,世界人口的极限值
0.029
12
b 2.695 10
107.6 (亿)
若以1965年的人口数
3.34 9
10为初值,则2000年的世界人口将达到
y
t 2000
9
0.029 3.34 10
0.009 0.02e
0029(2000 1965)
59.6
(
亿)
m a bN , n p qN
其中a 、b 、p 、q 均为正常数。
这说明出生率与死亡率不再是常数,而是
,前者随N 均匀减少,后者随 N 均匀增加。
于是方程(4.10)就被修 正为
分离变量并积分,得
Ce Akt
无限增加,而是逐渐接近于一个固定的上限
A ,因此它还是比较符合实际情况的。
事实上,
自然界各种动物(包括人类)的生长规律,以及传染病的发病人数、商品销售量的变化趋势 等问题,都可以用这个模型来描述。
形如(4.11)的微分方程称为逻辑斯蒂 (Logistic )方程,该方程的解(4.12)所表示的曲
线称为逻辑斯蒂曲线,即所谓的“
S 型曲线”。
N 的线性函数
(这是最简单的情形) dN “
[(a dt
[(a
bN) (P qN)] N P) (b q)N] N (b q)(a
若记k b q 、A
a p
,则上式可简记为
q
这就是描述动物种群数量
dN dt
kN(A
N)
(4.11)
N 与时间t
之间函数关系的微分方程。
将初始条件N
t0
N o 代入上式,可求得
—,故方程(4.11)
A N o
的特解即为
解出N 并化简, 可得
No A N
Akt
-e
A kt
AN 0 e A N o N o e Akt
A
A A kt
1 ( 1)e
A
1 Be t
(4.12)
A
其中B —
N 。
Ak 。
由于何
lim t 1 Be t
A ,这说明该种群中动物的数量
N 不再随时间的推移
三. 传染病模型
现代医学的发展已经能够有效地控制和预防许多传染病, 个世纪就已经在世界围被消灭,
鼠疫、霍乱等恶性传染病也得到了控制,
但是一些常见的传
染病并没有完全绝迹,一遇到合适的气候就会不断地爆发和流行, 至今仍然威胁着人类。
而
快,导致很高的致残率和死亡率, 危害极大,因而对传染病在人群中传播过程的定量研究并 应用于传染病的
防治就具有了重要的现实意义。
传染病流行过程的研究与其它学科的研究有所不同, 科学的数据只能通过以往传染病流
行时的报告来获取而不能在人群中实验, 而报告的数据往往难于满足研究的需要, 因而通过 建
立数学模型与计算机仿真就成为研究传染病流行过程的重要途径之一。
1.相对封闭环境中的传染病模型
假设N 个人共同生活在一个相对封闭的环境中,
如果其中的一个人感染了某种传染病,
而这种传染病又有一定的潜伏期,因而未发病时人们(包括患者自己)
过程是怎样的呢?
时间被感染的人数也比较少,
因而传播速度较慢;而在传染病流行的后期,
由于大多数人已
经被感染(I (t)较大),未被感染的人数 N I (t)已经不多了,所以此时单位时间被感染的 人数也不多,因而传播速度也很慢。
排除上述两种极端的情况,当有很多的感染者和很多的 未感染者时,传染病的传播速度是很快的。
因此,传染病的发病率一方面受感染人数的影响,
另一方面也受未感染人数的制约。
基于以上的分析,我们就可以建立微分方程如下:
其中k 是比例常数,可根据发病情况的统计数据来确定。
(注)
不难看出,这个方程与(4.11 )完全相同,因而它的通解为
古 Ce
Nkt
—丄 e
Nkt
如以前曾经肆虐一时的天花上
且随着时间的推移,还会出现一些新的传染病,如
AIDS 、SARS 等。
有些传染病传染得很
是不知道的。
一旦患 者发病被隔离时,
这种传染病实际上已经在传播了。
在这种情况下,
传染病在人群中的传播
如果我们以
I (t)表示发现首例病人后
t 时刻被感染的人数,则
N I (t)就表示此时刻
未被感染的人数。
在传染病流行的初期,由于 l(t)较小,能接触到感染者的人数少,单位
dl dt
kI(N I)
(4.13)
注意到I (0)
1,代入上式,可求得 C
故该方程的特解为
N (N 1)e
Nkt
定的距离,但在传染病流行的前期这个模型还是可用的, 因而传染病学者曾用它来预报过传
染病高潮到来的时刻。
解(4.14),即可求出比例常数 k 的值。
)
2.传染病的SIR 模型
在上一个模型中,我们实际上只是对人群进行了简单的划分: 感染者与健康人。
如果我
们把病愈免疫和死亡者也考虑在,情况就会有很大的不同。
首先,我们把人群分成三类:
S 类:称为易感类,该类成员没有染上传染病,但缺乏免疫力,可以被感染; I 类:称为传染类,该类成员已经染上了传染病,而且可以传播给 S 类成员;
R 类:称为恢复类或排除类,该类成员具有免疫能力或者已经死亡。
记以上三类人员在 t 时刻的人数分别为 S(t)、I (t)和R(t),总人数为 N ,显然在任 一时刻,都有
S(t) I(t) R(t) N
为了建立数学模型,我们作出以下的两点假设: (1)
单位时
间,一个病人传染的人数与当时健康者的人数成正比,比例系数为 k (称
为传染系数);
(2) 单位时间,病愈免疫(包括死亡)的人数与当时感染者的人数成正比,比例系数
为(称为恢复系数)。
由假设(1),可知S 类成员中在t 时刻的单位时间被感染的人数为 kS(t) I(t);由假 设(2),可知t 时刻的单位时间由I 类成员转化为 R 类的人数为 I (t)。
由于t 时刻单位时 间各类人员的人数就是该类人员总人数的变化率,故有
dS dt
如果初始时刻的感染人数为 I 0,即I (0) I 0,则S(0) N I 0 (此时R(0) 0 ),于 是我们就可以得到如
下的传染病模型:
(4.14)
这就是该传染病的传播规律。
显然, I (t)单调增,且当t 时,有I A ,这表明
如果任其发展而不采取积极有效的措施,
最终所有的人都将被感染。
尽管这与实际情况有
(注:假设发现首例病人
24小时后被感染的人数为
3人,则将t
24、 I
3代入特
dI dt
kSI
kSI
(4.15)
虽然微分方程组(4.15)的解析解难以求得,但是我们可以通过对 I 、S 之间函数关系
的研究讨论(4.15)解的性态。
为此,将(4.15)中的两个方程相除,得
dI kSI I
1 dS kSI kS
若将 记作 ,则称
为“特征指数”,且对于同一地区、同一种传染病,
为常数。
k
k
于是方程(4.17)可表示为
dI
1
( 4.18)
dS S
其初始条件为
I (S o ) N S o
方程(4.18)是一个可分离变量的微分方程,易求得其满足初始条件的解为
S
I (S ) In NS
件 19)
S o
(4.19)表示的只是患者数量|与健康人数量 S 之间的函数关系,为了建立它们与时刻
t 之间的函数关系以便作进一步的分析,我们可以利用(
I (t) (— 1) S(t)
时,I (t ) o ,这时l (t )单调减少。
这也就是说,仅当传染
病开始流行时健康人数 S o 超过 的情况下,传染病才会蔓延。
因此 是一个“阀值”(俗
称“门槛”)。
注意到S o N I o ,而通常I o 是很小的,所以可近似地认为 S o N ,于是
在总人数N 不变的情况下,提高门槛
的数值,对制止传染病的蔓延是十分有利的。
由于
匚,因此一方面要及时采取隔离措施,想方设法降低传染系数 大恢复系数 ,不
断提高该地区的医疗水平、免疫水平和卫生保健水平,以避免传染病的大 规模爆发和流行。
其初始条件为
1(0) I o ,
di
dt dS
dt
kSI I kSI S(0) N I o
(4.16)
(4.17)
4.19)式求I (S )对t 的导数:
dl dt
dl dS dS dt
1)
dS dt
由于S (t )是单调减函数,即
S (t ) o ,故由上式可以看出,当
时,I (t) o ,
从而I (t )单调增加;当S
k ,另一方面要努力增
(4.23)
四. 弱肉强食模型
在生态环境中,所有的动物都要通过食物来生存, 不是吃其它的动物,就是吃植物。
从 这种意义上讲,没有一种动物是完全孤立生活的,
而且彼此之间一定是相互影响的。
为了搞
清楚它们之间相互依存的关系,就必须把它们合在一起进行研究,这就产生了多种群模型。
弱肉强食模型就是其中较简单的一种。
假设有一个包含两个群体的系统,
其中的一个以另一个为食,
比如狐狸与兔子、 鲸鱼与
磷虾、天敌与害虫、正常细胞与肿瘤细胞……等等,我们把这种系统笼统地称为“捕食者- 被食者”系统。
设x(t)、y(t)分别表示t 时刻被食者与捕食者的数量。
如果它们各自单独生活,则被 食者的增长速率正比与当时的数量,即
dx
x
dt
而捕食者由于没有捕食对象,其数量减少的速率也正比与当时的数量,即
(4.22)中的两个方程相除,得
(0)
(4.20)
dy dt
(0)
(4.21)
现在两者生活在一起,被食者有一部分遭到了捕食者的捕食,增长速率要变缓,即 (4.20)式
中的 将减少,而且减少的量正比于捕食者的数量,所以
(4.20)式应修改为
dx dt
y) x
(0)
类似地,捕食者有了捕食对象,减少的速率也要变缓,即 (4.21)式中的
也减少,而且减少
的量正比于被食者的数量,所以
(4.21)式应修改为
牛
(x)y
将修改后的两个方程联立,可得方程组
赵(y) x dt
f
( x)y
称为Volterra -Lotka 方程,其初始条件是
x(0) X 。
,
y(0) y °
(0)
(4.22)
虽然这个方程组是非线性方程组, 不易直接求出它的解析解, 但是我们可以采用传染病 模型中使用过的方法,通过对
y 之间函数关系的研究讨论(
4.22 )解的性态。
为此,将
生态的平衡。
F 面我们再进一步讨论这个模型的另一个问题:假如由于某种外界的因素(比如天灾)
使得捕食者与被食者同时按比例消亡,那么一旦这种因素的作用停止,是被食者恢复得快, 还是捕食者恢复得快呢?
中,乘积项xy 对被食者来说是消亡项,对捕食者来说是生长项。
当 如说x 、y 都减
半时,则 x 、 y 也随之减半,但 xy 和 xy 却只减为原来的四分之一。
这说明数量减半以后,被食者增长速度的减少还不到一半, 而捕食者生长的速度却开始变缓, 因而被食者受到的损失比捕食者要小。
所以灾难过后,被食者的恢复要比捕食者恢复得快。
这个推断来自 Volterra 方程,因而也称之为 Volterra 原理。
这个原理在杀虫剂的施用方面得到了证实。
1968年,澳洲的吹棉蚧(一种昆虫)偶然 传入美国后,严重地威胁着美国的柑桔业,为此, 人们又引入了它的天敌澳洲瓢虫,才使吹 棉蚧的数量得到了控制。
后来,园艺家们使用了
DDT ,以期能进一步减少吹棉蚧。
然而,
分离变量并积分,可得通解
利用 y In e y
及 x In e
dy ( x )y dx
In y
y)x
In x In C
x
,可将上式化为
这就是方程(4.23)的隐式通解;当 x x 0、
y o 时,C
y o
e yo
X o
e 乂°
就是方程(4.23 )的一个隐式特解,它的图象是 要初始值X 。
和y o 不为零,这条闭曲线就永远不 通过零点。
当被食者较多时,由于食物丰富使捕 食者不断增加,于是被食者又开始减少;被食者 减少使捕食者也随之减少,从而被食者又会增多。
如此,两者的数量循环起伏,周而复始,维持着
y o
y o
xy 平面上的一条闭曲线(示意图如下)。
只
这个问题可以从方程(
4.22) 中找到答案。
在方程组(4.22)
dx dt dy dt
y) x x)y
dx
dt dy dt
xy xy
x (t )、y (t )都减小,譬
(4.23)
由于DDT 在杀灭吹棉蚧的同时也杀死了它的天敌,一段时间过后,人们发现吹棉蚧的数量 反而增多了。
上述原理在临床医学上也是有意义的。
癌症患者体的正常细胞与癌细胞就构成了一个 “捕食者-被食者” 系统,某些具有较高毒性的药物既能抑制癌细胞的生长,
同时也能抑制
正常细胞的生长。
癌症患者在进行化疗的时候,药物的毒性使作为捕食者的正常细胞也受到 了较大的伤害,一定停止用药,作为被食者的癌细胞可能会变本加厉地增长起来。
这始终是
癌症化疗中需要解决的一个难题。
五•溶液浓度模型
溶液浓度问题是工农业生产和治理环境污染中经常要碰到的问题,
描述为如下的实验室模型:一个容器有一个入口,
一个出口,里面盛满了某种溶液。
如果从
于是,在时间间隔
[t,t dt],有
从而含盐量的微元即为
t
In 2
0.6931 0.005 0.005
即经过约2小时18分37秒可使容器的含盐量达到 100公斤。
入口以不变的速率向容器注入一定浓度的相同溶液(或清水) ,搅拌均匀后以同样的速率从
出口排出,假设搅拌是在瞬间完成的,那么容器溶液浓度的变化规律是怎样的呢? 为了讨论的方便,我们来看一个具体的例子。
已知容器盛有1000公斤的清水,如果以每分钟5公斤的速率注入浓度为 0.2的盐水且
不停地搅拌,并以同样的速率排出搅拌后的盐水, 那么经过多少时间能使容器的含盐量达到 100公斤? 这个问题我们可以采用微元分析的方法建立微分方程求解。
为此,可以设 t 时刻容器 的含盐量为y(t),则此时溶液的浓度为 y(t)
1000
0.001 y(t),
进盐量:
0.2 5 dt dt
出盐量: 0.001y 5 dt
0.005ydt
此类问题通常都可以
dy dt 0.005ydt
业 0.005y
1
dt
这是一个一阶线性非齐次微分方程,易求得该方程满足初始条件为
y(0) 0的特解为
y 200(1 e
0.005t
这就是容器的含盐量 y 随时刻t 变化的规律。
将 y 100代入上式,可求得
138.62 (分)。