传染病数学模型的建立及稳定性分析

《传染病数学模型的建立及稳定性分析》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、传染病数学模型的建立及稳定性分析数理医药学杂志990202数理医药学杂志JournalofMathematicalMedicine1999年　第12卷　第2期　Vol.12　No.2　1999传染病数学模型的建立及稳定性分析张　运　权摘要　利用疾病传播的一般规律及人口守恒统计法测建立起两室与三室的传染病模型［1］，再运用文献［2、3］的数学方法，重点对两室的传染病模型进行定性与稳定性分析，从而得出相应情况下的生态意义。　　关键词　稳定性　轨线　阈值定理1　预备知识1.1　两室的模型　把城市人口分为健康人与传染病人两个室(集合)，其人数分别记作：S(t),I(t)。1.2　三室的模型　把城市人口

2、分为健康人S(t)、传染病人I(t)及病愈免疫(包括死亡)的人R(t)。　　我们知道疾病传播一般服从下列法则：　　法则1　在所考虑的时期内，人口总数保持在固定水平N(即S(t)+I(t)+R(t)=N)。　　法则2　易受传染者S(t)人数的变化率正比于传染病患者I(t)与S(t)人数的乘积。　　法则3　由I(t)向R(t)转变的速率与I(t)成正比。2　两室的模型　　由上述疾病传播法则，不难得出传染病的数学模型(1)且初始状态为S(0)=S0>0,I(0)=I0>0其中常数λ、α称为传染率、移除率，其值均大于零。　　令σ=λ/α,1/σ=α/λ称为相对移除率，同时为了讨论问题的方便，不妨假设N

3、=1，即总体。　　定理1　(阈值定理)设S(t),I(t)是初值问题(1)的解，如果σS0<1,那么，当t→+∞时，I(t)单调减少趋于零。如果σS0>1，当t→+∞时，I(t)先增加达到最大值1-1/σ-1/σln万方数据file:///E

4、/qk/slyyxzz/slyy99/slyy9902/990202.htm（第1／5页）2010-3-237:21:16数理医药学杂志990202(σS0),此时S=1/σ,而后单调减少趋于零，S(t)是一个单调减少函数，并且其极限，是方程1-S+(ln(S/S0))/σ=0在(0，1/σ)内的根(见图1)。图12.1　在不考虑自然出生和死亡的前提下，

5、一种传染病发生时，如果易感染人的总数小于等于该病的相对移除率，此传染病不可能发生流行，将很快被消灭。如果易感染人的总数大于该病的相对移除率，此传染病可能发生流行，得病的人数将猛增，当易感染人数下降到S=1/σ时，得病人数I(t)达到最大值1-1/σ-1/σln(σS0)，而后得病人数逐渐减少，最终趋向于零，即传染病被消灭。在整个过程中易感染人数单调减少，最终并不是所有的易感染人都会得病。因此，我们说疾病不是因为缺少受传染者而停止传播，而是因为没有了传染者才停止传播。2.2　由于人类对传染病的认识提高以及现代医学水平的发展，对于许多传染病可以做到提前预防，使人群对许多种传染病具有免疫能力，例如打

6、预防针、进行免疫接种等，我们在模型(1)的基础上增加考虑直接进入消除类的因素，使其注射预防针的速率λ与I(t)成正比。　　经过调整，得到如下的数学模型：万方数据file:///E

7、/qk/slyyxzz/slyy99/slyy9902/990202.htm（第2／5页）2010-3-237:21:16数理医药学杂志990202(2)且S0>0、I0>0,其中λ,α，δ均大于零，令σ=λ/α,1/σ=α/λ为相对移除率。　　以下对模型(2)进行分析。　　我们在(S，I)相平面上考察轨线。　　首先由　　(dS)/(dt)=-λIS-δS=-(λIS+δS)<0得S(t)是单调减少的，所以，对于所有

8、t>0有S(t)0成立，此时I(t)单调减少。当S0>(1)/(σ)时由S(t)单调减少可得存在唯一的t1,使S(t1)=(1)/(σ),因此有：00，此时I(t)增加，t>t1后(dI)/(dt)<0,此时I(t)单调减少，所以I(t1)是最大值。即在相平面上的轨线I(r)在r=(1)/(σ)时，I((1)/(σ))为最大值。　　显然方程组(2)的轨线方程为：由此我们得：上式说明采取预防措施后，可以减少得病人数，并且I(t)的最大值小

9、于不采取预防措施时的最大值。　　令　D={(S,I)

10、0

11、/qk/slyyxzz/slyy99/slyy9902/990202.htm（第3／5页）2010-3-237:21:16数理医药学