传染病传播的数学模型---上课

1、微分方程模型学习目的1,加深对微分方程概念的理解，掌握针对一些问题通过建立微分方程的方法及微分方程的求解过程；2,了解微分方程模型解决问题思维方法及技巧；3,领会建立微分方程模型的逐步改进法的核心及优点，并掌握该方法；4,理解微分方程的解的稳定性的意义，会用稳定性判定模型的解是否有效；5.体会微分方程建摸的艺术性。在自然学科（如物理、化学、生物、天文）以及在工程、经济、军事、社会等学科中大量的问题可以用微分方程来描述。正如列宁所说：“自然界的统一性显示在关于各种现象领域的微分方程式的惊人的类似中."（列宁选集第二卷，人民出版社1972年版第295页）;要建立微分方程模型，读者必须掌握

2、元素法（有关元素法，在高等数学中已有介绍）。所谓元素法，从某种角度上讲，就是分析的方法，它是以自然规律的普遍性为根据并且以局部规律的独立的假定为基础。在解决各种实际问题时，微分方程用得极其广泛。读者通过下面的几个不同领域中的模型介绍便有所体会，要想掌握好它，在这方面应作大量的练习。§ 17.1 传染病传播的数学模型学习目标1,通过学习建立传染病传播的数学模型的思维方法，能归纳出该类建模的关键性步骤及思维方法；并能指出求解传染病传播的数学模型的方法技巧；2,能用已知的传染病传播的数学模型，预报某种传染病的传播；3.学会从简单到复杂的处理问题的方法。由于人体的疾病难以控制和变化莫测，因此

3、医学中的数学模型较为复杂。生物医学中的数学模型分为两大类：传染病传播的数学模型和疾病数学模型。以下仅讨论传染病的传播问题。人们将传染病的统计数据进行处理和分析，发现在某一民族或地区，某种传染病传播时，每次所涉及的人数大体上是一常数。这一现象如何解释呢？关于这个问题，医学工作者试图从医学的不同角度进行解释都得不到令人满意的解释。最后由于数学工作者的参与，在理论上对上述结论进行了严格的证明。同时又由于传染病数学模型的建立，分析所得结果与实际过程比较吻合，这个现象才得到了比较满意的解释。传染病传播所涉及的因素很多，如传染病人的多少，易受传染者的多少，传染率的大小，排除率的大小，人口的出生和死亡等。如

4、果还要考虑人员的迁入与迁出，潜伏期的长短以及预防疾病的传播等因素的影响，那么传染病的传播就变得非常复杂。如果一开始就把所有的因素考虑在内，那么将陷入多如乱麻的头绪中不能自拔，倒不如舍去众多的次要因素，抓住主要因素，把问题简化，建立相应的数学模型。将所得结果与实际比较，找出问题，修改原有假设，再建立一个与实际比较吻合的模型。下面由简单到复杂将建模的思考过程作一示范，读者可以从中得到很好的启发。模型一、考虑最简单的情形：欢迎共阅假设(1),每个病人在单位时间内传染的人数是常数K。；假设(2),一人得病后，经久不愈，并在传染期内不会死亡。记i(t)表示t时刻病人数，K。表示每个病人单位时间内传染的人

5、数，i(0)=i。，即最初有认个传染病人。则在At时间内增加的病人数为'di(t)于是得微分方程=Koi(1),其解为i(t)=ioekoti(0)=i。结果表明：传染病的传播是按指数函数增加的。这个结果与传染病传播初期比较吻合，传染病传播初期，传播快，被传染人数按指数函数增长。但由方程(1)的解可以推出，当tT8时，i(t)T8,这显然是不符合实际情况的。问题在于两条假设均不合理。特别是假设(1),每个病人在单位时间内传染的人数是常数与实际不符。因为在传播初期，传染病人少，未被传染者多；而在传染病传播中期和后期，传染病人逐渐增多，未被传染者逐渐减少，因而在不同时期的传染情况是不同的。

6、为了与实际情况吻合，我们在原有基础上修改假设建立新的模型。模型二、用i(t),s(t)表示t时刻传染病人数和未被传染人数，i(0)=io。假设(1),每个病人单位时间内传染的人数与这时未被传染的人数成正比，即Ko=Ks(t)；假设(2),一人得病后，经久不愈，并在传染期内不会死亡；假设(3),总人数为n,即s(t)+i(t)=n.由以上假设得微分方程T=Ks(t)i(t)dt，s(t)+i(t)=n(2)i(o)=i0.'.l1L|用分离变量法求得其解为i(t)=n(3)n-Knt1+I-1eUo)其图形如图17.1所示.模型二可以用来预报传染较快的疾病前期传染病高峰到来的时间。图17

7、.1图17.2医学上称di/(dt-t)为传染病曲线，它表示传染病人增加率与时间的关系，如图17.2所示。由(3)式可得欢迎共阅Kn2didtnJi-1ei0Knt_Knt(4)1令.。dt2,得极大点为tiln(n1)10(5)Kn由此可见，当传染病强度K或总人数n增加时，匕都将变小，即传染病高峰来得快，这与实际情况吻合。同时，如果知道了传染强度K(K由统计数据得出)，即可预报传染病高峰匕到来的时间,这对于防治传染病是有益处的。模型二的缺点是：当tT吧时，由(3)式可知，i(t)Tn,即最后人人都要生病，这显然是不符合实际情况的。造成该问题的原因是假设(2)中假设了人得病后经久不愈。为了与实

8、际问题更加吻合，对上面的数学模型再进一步修改，这就要考虑人得了病后有的会死亡；另外不是每个人被传染后都会传染别人，因为其中一部分会被隔离。还要考虑人得了传染病由于医治和人的自身抵抗力会痊愈，并非象前面假设的那样，人得病后经久不愈。为此作出新的假设,建立新的模型。模型三在此模型中，虽然要考虑比前面两个模型复杂得多的因素，但仍要把问题简单化。设患过传染病而完全痊愈的任何人具有长期免疫力，不考虑反复受传染的情形，并设传染病的潜伏期很短，可以忽略不计，即一个人患了病之后立即成为传染者。在这种情况下，把居民分成三类：第一类是由能够把疾病传染给别人的那些传染者组成的，用I(t)表示t时刻第一类人数；第二类

9、是由并非传染者但能够得病而成为传染者的那些人组成的，用s(t)表示t时刻第二类人数；第三类包括患病死去的人，病愈后具有长期免疫力的人，以及在病愈并出现长期免疫力以前被隔离起来的人，用R(t)表示t时刻第三类人数。假设疾病传染服从下列法则:(1)在所考虑的时期内人口总数保持在固定水平N,即不考虑出生及其它原因引起的死亡，以及迁入迁出等情况；(2)易受传染者人数s(t)的变化率正比于第一类的人数I(t)与第二类人数s(t)的乘积；(3)由第一类向第三类转变的速率与第一类的人数成正比。由(1)、(2)、(3)条得微分方程组ds=rsIdtdI"=rsl-3(6)dt岖=3dt其中r、尸为两

10、个比例常数，r为传染率，尸为排除率。由(6)式的三个方程相加得则s(t)+I+R(t)=N(人口总数)故R(t)=N-s(t)-I(t)欢迎共阅由此可知，只要知道了s和I，即可求出R(t)o而(6)式的第一和第二个方程与R(t)无关。因此，由dsdtdl.dt-rslrsl一I(8)得M*一右,、，尸，I(s)=s+Ins+c。rV_当t=10时，1(I)=10,s(t0)=s0,记P=，有r.一sI(s)=10+s0s+Pin(9)s。下面讨论积分曲线(9)的性质。由(8)式知所以当s<P时，I(s)是s的增函数，s>P时，I(s)是s的减函数。由连续函数中间值定理及单调性知，存

11、在唯一点s»0cse<.%使得I(s*)=0。而当s3C<ss0时，I(s)>0。由知I=0时，ds/dt=0,dI/dt=00/qII1IJ"所以(s*0)为方程组(7)的平衡点.图173当t之t01寸，方程(9)的图形如图17.3如果s0>P,则随着s减小到P明呵眺到<1s店跚叫假党金。丁s(t)<P时，才开始减小。由以上分析可以得出女郎缚花方阴缈当属困用啊恻传麋啊的瞰超过阈值p=Z时二I因此，如果s0小于P,则I(t)单调减小到零，s(t)单调r传染病才会蔓延。用一般的常识来用涌嘤沙卮也抑必4gB力逃缺乏应有的科学文化知识，缺乏必要

12、的医疗条纥隔离乖而排除率裕日%传染病会很快他；反之，人口密度低，社会条件好，有良好的公共3111曲汁较1的斡疝而翟耀肃很快被的较病在有限范围内出现十.炉金盆新国以得出如下结论：却俱一被出火。如果起初易受传染者的人数s0大于但接近于阈值P,即如果(80P)与P相比是小量，则最终患病的人数近似于2(s0-P),这就是着名的传染病学中的阈值定理。生物数学家Kermack和.,Mekendrick在1927年首先证明了这个定理。定理(传染病学中的阈值定理)：设s0=P+r,且假设r/P同1相比是小量。并设最初传染者人数I0很小，则最终患病的人数为2r。即易受传染者的人数最初比阈值高多少，那最终就会比阈

13、值低多少。证明略。根据阈值定理就可以由起初易受传染者的人数来估计最终患病的人数。这个定理解释了研究人员长期以来难以解释的为什么对于某一民族或地区，某种传染病传播时，每次所波及的人数大体上是一常数的现象。欢迎共阅在传染病发生过程中，不可能准确的调查每一天或每一星期得病的人数。因为只有那些来医院就医者才能被人知道他们得了病，并把他们隔离起来防止传染。因此，统计的记录是每一天或每一星期新排除者的人数，而不是新得病的人数。所以，为了把数学模型所预示的结果同疾病的实际情况进行比较，必须解出(6)式中的第三个方程：dsdsdR-rsI因为=dRdtdtI所以s(R);s0e®drssdsdR有空

14、=(N一R-s0e-R/：J(10)dt方程(10)虽是可分离变量的，但是不能用显式求解。如果传染病不严重，则R/P是小量，取泰勒级数e"=1.R1:2十一的前三项，取近似值得其解为R(t)=：2-1I),2s(Ns0FP因止匕dRdt:''sec'12soa¥t邛i(11)方程(11)在t-dR/dt平面上定义了一条对称钟形曲线，称为疾病传染曲线.疾病传染曲线很好的说明了实际发死亡人数生的传染病。每天报告的新病案的数目逐渐上升到峰值，然后又减少下来。dR一-1000dt八800600400200Xo12tp510152030星期数ar图17.4图1

15、7.5Kermack和Mekendrick把(11)得到的dRdt的值，同取自1905年下半年至1906年上半年在孟买发生的瘟疫资料进行比较，他们设其中t按星期计，在图17.5中，dR/dt的实际数字(图上用一表示)同理论曲线非常一致。这就表明了模型三是在固定的居民中传染病传播的准确而可靠的数学模型。对于同一事物，可用不同的数学工具来描述它。下面介绍一般随机传染病模型。模型四、一般随机传染病的数学模型:以上建立的常微分方程描述的传染病的传播是确定性的模型。但人生病是随机的，因而建立随机的传染病的数学模型才能更实际的反映传染病的传播。n表示易受传染者总数,设X(t)表示t时刻易受传染者人数，Y(

16、t)表示t时刻已受传染者人数,欢迎共阅又设t时刻有i(i>0)个易受传染者移入已受传染者中来。这种传染病传播的机制如下：(1)在群体中个体均匀的混和；在区间(t,t+担)内，一个新传染病例出现的概率为ZxyAt+o(At),其中八(九>0)是传染率;(3)在区间(t,t+At)内，排除一个个体的概率为RyAt+o(At),其中N(N0)是排除率；(4)在区间(t,t+&)内，有多次转移(即多个传染或排除)发生的概率为o(At);(5)在区间(t,t+&)内,无变化的概率为1(Xx+H)yAt+o(At),令Pxy(t)=6X(t)=x,Y(t)=y,x>0,y

17、之0,t之0。从(2)和(3)知有两种可能的转移(xtx-1,y->y+1)>(y-)y-1)0因此，表征随机传染病流行过程的差分方程为：'dPxy(t)=(x+1)(y-1)Px书,y/t)-y(x+P)Pxy(t)+P(y+1)Px,y书(t)出(12)dPni(t)门'J=_i(n+P)%(t)、dt其中0Wx+yWn+i,0Mx<n,0<y<n+i0在(12)式中P=N/K为相对排除率，并已对时间变量作了变换，使方程对传染率九是无量纲的。'方程(12)的初始条件为；/(6)。方程(12)所描述的随机传染病流行数学模型可以推得确定性模型1942年Wilson和Burker讨论了潜伏期的重要性，他们用微分差分方程来描述传染病传播的数学模型.其中A是易受感染者乘以恢复健康的比率，仃是易受感染者变成传染者的潜伏期，r是易受感染者与病人的接触率。Cooke除了考虑潜伏期，还引进了阻尼阈值的概念，这个概念说明个别成员从易受感染者成为传染者之前可能反复发病，Cooke模型