传染病传播的数学模型.docx

传染病传播的数学模型.docx

《传染病传播的数学模型.docx》由会员分享，可在线阅读，更多相关《传染病传播的数学模型.docx（12页珍藏版）》请在冰豆网上搜索。

传染病传播的数学模型

【下载本文档，可以自由复制内容或自由编辑修改内容，更多精彩文章，期待你的好评和关注，我将一如既往为您服务】

传染病传播的数学模型

很多医学工作者试图从医学的不同角度来解释传染病传播时的一种现象，这种现象就是在某一民族或地区，某种传染病传播时，每次所涉及的人数大体上是一常数。

结果都不能令人满意，后来由于数学工作者的参与，用建立数学模型来对这一现象进行模拟和论证，得到了较满意的解答。

一种疾病的传播过程是一种非常复杂的过程，它受很多社会因素的制约和影响，如传染病人的多少，易受传染者的多少，传染率的大小，排除率的大小，人口的出生和死亡，还有人员的迁入和迁出，潜伏期的长短，预防疾病的宣传以及人的个体差异等。

如何建立一个与实际比较吻合的数学模型，开始显然不能将所有因素都考虑进去。

为此，必须从诸多因素中，抓住主要因素，去掉次要因素。

先把问题简化，建立相应的数学模型。

将所得结果与实际比较，找出问题，修改原有假设，再建立一个与实际比较吻合的模型。

从而使模型逐步完善。

下面是一个由简单到复杂的建模过程，很有代表性，读者应从中体会这一建模过程的方法和思路。

1.最简单的模型

假设：

（1）每个病人在单位时间内传染的人数是常数k；

（2）一

个人得病后经久不愈，并在传染期内不会死亡。

以i（t）表示t时刻的病人数，ko表示每个病人单位时间内传染的人数，i（0）=io表示最初时有io个传染病人，则在.t时间内增加的病人数为

itt-it二koitt

两边除以上t,并令't-o得微分方程

k0t

其解为it-ioe

这表明传染病的转播是按指数函数增加的。

这结果与传染病传播初期比较吻合，传染病传播初期，传播很快，被传染人数按指数函数增长。

但由（2.1）的解可知，当t-K时，i⑴这显然不符合实际情况。

最多所有的人都传染上就是了。

那么问题在那里呢？

问题是就出在于两条假设对时间较长时不合理。

特别是假设

（1），每个病人

单位时间内传染的人数是常数与实际情况不符。

因为随着时间的推移，病人越来越多，而未被传染的人数却越来越少，因而不同时期的传播情况是不同的。

为了与实际情况较吻合，我们在原有的基础上修改假设建立新的模型。

2.模型的修改

将人群分成两类：

一类为传染病人，另一类为未被传染的人，分别用i（t）和s（t）表示t时刻这两类人的人数。

i（O）=io。

假设：

（1）每个病人单位时间内传染的人数与这时未被传染的人数成正比。

即k^kst；

（2）—人得病后，经久不愈，并在传染期内不会死亡。

由以上假设可得微分方程

dit

kstitdt

stit二n

i0-i

I

这是变量分离方程，用分离变量法可求得其解为

n

it二

1--1eknt……

io

（2.2）

（2.3）

其图形如下图2-1所示

模型（2.2）可以用来预报传染较快的疾病前期传染病高峰到来的时询。

医学上称dt_t为传染病曲线，它表示传染病人的增加率与时

dt

间的关系，如图2-2所示。

由（2.3）式可得

d2i（t）d2i（t）

再求二阶导数dt2，并令出2二0，可解得极大点为

In—-1It“0丿

（2.5）

t1

kn

从（2.5）式可以看出，当传染病强度k或人口总数n增加时，右

都将变小，即传染病高峰来得快。

这与实际情况吻合。

同时，如果知道了传染率k（k由统计数据得到），即可预报传染病高峰t1到来的时间,这对于预防传染病是有益处的。

模型（2.2）的缺点是：

当t-K时，由（2.3）式可知i（t）-n,即最后人人都要得病。

这显然与实袜情况不符。

造成这个结果的原因是假设

（2）中假设一人得病后经久不愈，也不会死亡。

为了得到与实际情况更吻合的模型，必须修改假设

（2）。

实际

上不是每个人得病后都会传染别人，因为其中一部份会被隔离，还有由于医治和人的身抵抗力会痊愈，有的人会死亡从而也就不再会传染给别人了。

因此必须对模型作进一步的修改，建立新的模型。

3.模型的进一步完善

从上面的分析我们看到模型（2.2）的假设

（2）是不合理的。

即不可能一人得病后会经久不愈，必有一部份人因医治或自身的免疫力，或是被隔离，或是死去而成为不会再继续传染给别人的第三类人。

因此我们把人群分成三类：

第一类由能够把疾病传染给别人的那些传染者组成的。

用l（t）

表示t时刻第一类人数。

第二类是由并非传染者但能够得病而成为传染者的那些人组成

的，用S（t）表示t时刻第二类人数。

第三类包括患病后死去的人，病愈后具有长期免疫力的人，以及在得病后被隔离起来的人。

用R（t）表示t时刻第三类人数。

假设疾病传染服从下列法则：

（1）在所考虑的时期内人口总数保持在固定水平N，即不考虑出

生及其他原因引起的死亡，以及人口的迁入迁出的情况。

⑵易受传染者人数S（t）的变化率正比于第一类的人数l（t）与第

二类人粉S（t）的乘积。

（3）由第一类向第三类转变的速度与第一类的人数成正比。

在这三条假设情况下可得如下微分方程：

其中r、入为比例常数，r为传染率，入为排除率

由方程（2.6）的三个方程相加得

d-n

—S（t）+I（t）+R（t）]=0dt」」

则StItRt常数二人口总数

故Rt二N-St-It

因此只要求出S（t）、l（t）即可求出R（t）。

方程组（2.6）的第一个和第二个方程与R（t）无关。

因此，由

dlrSI-I’

==—4r

dS-rSIrS

由初始条件：

当t二to时，Ito=Io,Sto=So并记

代入上式可确定常数C=10•S0-1nSo

最后得

IS=IoSo-STn

（2.9）

F面我们讨论积分曲线（2.9）的性质，由（2.8）知

S<

・>0

所以当Sp时，I（S）是S的减函数。

又有I（0）=―乂，IS0=1。

>0,由连续函数的中间值定理及单调性知，存在唯一点S：

：

0

：

=0，而当S0。

dScdlcgc

由（2.7）知1=0时，T二°，二=0，所以S,0为方程组

dtdt

（2.7）的平衡点。

当t-t0时，方程（2.9）的的图形如图2-3。

当t由to变到s时,点（S（t）,l（t））沿曲线（2.9）移动，并沿S减少的方向移动，因为S（t）随时间的增加而单调减少。

因此，如果So小于p，则|（t）单调减少到零，S（t）单调减少到S：

。

所以，如果为数不多的一群传染者Io分

|少一3

如果So」，则随着S（t）减少到P时，I⑴增加，且当S=p时，I（t）达到最大值。

当S（t）vP时I（t）才开始减少。

由上分析可以得出如不结论：

只有当居民中的易受传染者的人数超过阈值时传染

r

病才会蔓延。

用一般常识来检验上面的结论也是符合的。

当人口拥挤，密度高，缺少应有的科学文化知识，缺乏必要的医疗条件，隔离不良而排除率低时，传染病会很快蔓延；反之，人口密度低，社会条件好，有良好的医疗条件和较好的管理而排除率高时，则传染病在有限范围内

出现会很快被消灭。

丄

传染病学中的阈值定理设S。

二'•r,且假设7同1相比

是小量。

并设最初传染者人数Io很小，则最终患病人数为2r。

即是易受传染者的人数最初比阈值高多少，那么最终就会比阈值低多少。

这就是有名的传染病阈值定理。

生物数学家Kermack和Mekendrick在1927年首先证明了这个定理（证明从略）

根据阈值定理就可以由起初易受传染者的人数来估计最终患病的人数。

这定理解释了研究人员长期以来难以解释的为什么对于某一民族或地区，某种传染病传播时，每次所涉及的人数大体上是一常数的现象。

在传染病发生的过程中，不可能准确地调查每一天或每一星期的得病人数。

因为只有那些来医院就医者才能被人知道他们得了病，并把他们隔离起来防止传染。

因此，统计的记录是每一天或星期新排除者的人数，而不是新得病的人数。

所以，为了把数学模型所预示的结果同疾病的实际情况进行比较，必须解出（2.6）中的第三个方程。

dRdt

-\-N-R-S

dS

dSdRrS\rS

/S二-

dR

dtdt\

因为

dS

=

dR

S

p

R

所以

S

rre

fR

dR_mDc"p

从而有HjN—R—Soe（2.10）

方程（2.10）虽是可分离变量的方程，但是不能用显式求解，如果传染病不严重，则R/p是小量，取泰勒级数前三项有

其中

i

[0J^sz-s。

"2

=tanh'a

亠dR

方程

（2.11）在—前平面上定义了一条对称钟形曲线，称为疾病

传染曲线。

疾病传染曲线很好地说明了实际发生的传染病的情况：

每天报告的新病案的数目逐渐上升到峰值，然后又减少下来。

Kermak和Mekendrick把（2.11）得到的值，同取自1905年

下半年至1906年上半年在印度孟买发生的瘟疫资料进行比较，他们假设

dR

890sech20.2i-3.4

dt

其中t按星期计，在图2-4中的实际数字（图中用“表示）同理论曲线非常一致。

这就表明模型（2.6）是在固定居民中传染病传播的准确而可靠的数学模型。

对于传染病传播的数学模型还有人用随机模型，这

不是本章的内容，读者可参看有关的其他资料。

本节所介绍的传染病传播的数学模型的建模方法，是实际数学建模步骤和方法的典型例子。

在实际建模过程中往往都是从简单的开始得出数学模型，再和实

际比较逐步修改假设和模型，最终达到完善的地步。

这是值得大家仿效和学习的。