传染病数学模型的建立及稳定性分析_精品文档.pdf
- 文档编号:3177572
- 上传时间:2022-11-19
- 格式:PDF
- 页数:2
- 大小:85.48KB
传染病数学模型的建立及稳定性分析_精品文档.pdf
《传染病数学模型的建立及稳定性分析_精品文档.pdf》由会员分享,可在线阅读,更多相关《传染病数学模型的建立及稳定性分析_精品文档.pdf(2页珍藏版)》请在冰豆网上搜索。
传染病数学模型的建立及稳定性分析张运权(湖北工学院武汉430068)摘要利用疾病传播的一般规律及人口守恒统计法测建立起两室与三室的传染病模型1,再运用文献2、3的数学方法,重点对两室的传染病模型进行定性与稳定性分析,从而得出相应情况下的生态意义。
关键词稳定性轨线阈值定理1预备知识1.1两室的模型把城市人口分为健康人与传染病人两个室(集合),其人数分别记作:
S(t),I(t)。
1.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两室的模型由上述疾病传播法则,不难得出传染病的数学模型dSdt=-KISdIdt=KIS-AI
(1)且初始状态为S(0)=S00,I(0)=I00其中常数K、A称为传染率、移除率,其值均大于零。
令R=KA,1R=AK称为相对移除率,同时为了讨论问题的方便,不妨假设N=1,即总体。
定理1(阈值定理)设S(t),I(t)是初值问题
(1)的解,如果RS01,当t+时,I(t)先增加达到最大值1-1R-1Rln(RS0),此时S=1R,而后单调减少趋于零,S(t)是一个单调减少函数,并且其极限limt+S(t)S(+),是方程1-S+ln(S/S0)R=0在(0,1R)内的根(见图1)。
2.1在不考虑自然出生和死亡的前提下,一种传染病发生时,如果易感染人的总数小于等于该病的相对移除率,此传染病不可能发生流行,将很快被消灭。
如果易感染人的总数大于该病的相对移除率,此传染病可能发生流行,得病的人数将猛增,当易感染人数下降到S=1R时,得病人数I(t)达到最大值图11-1R-1Rln(RS0),而后得病人数逐渐减少,最终趋向于零,即传染病被消灭。
在整个过程中易感染人数单调减少,最终并不是所有的易感染人都会得病。
因此,我们说疾病不是因为缺少受传染者而停止传播,而是因为没有了传染者才停止传播。
2.2由于人类对传染病的认识提高以及现代医学水平的发展,对于许多传染病可以做到提前预防,使人群对许多种传染病具有免疫能力,例如打预防针、进行免疫接种等,我们在模型
(1)的基础上增加考虑直接进入消除类的因素,使其注射预防针的速率K与I(t)成正比。
经过调整,得到如下的数学模型:
dSdt=-KIS-DSdIdt=KIS-AI
(2)且S00、I00,其中K,A,D均大于零,令R=KA,1R=AK为相对移除率。
以下对模型
(2)进行分析。
我们在(S,I)相平面上考察轨线。
首先由99数理医药学杂志1999年第12卷第2期dSdt=-KIS-DS=-(KIS+DS)0有S(t)S0,再由方程dIdt=KIS-AI=KI(S-1R)可知,当S01R时,有dIdt0成立,此时I(t)单调减少。
当S01R时由S(t)单调减少可得存在唯一的t1,使S(t1)=1R,因此有:
0t0,此时I(t)增加,tt1后dIdt0,此时I(t)单调减少,所以I(t1)是最大值。
即在相平面上的轨线I(r)在r=1R时,I(1R)为最大值。
显然方程组
(2)的轨线方程为:
I(t)+lnI(t)I(0)=1-S+1Rln(SS0)由此我们得:
I(1R)=1-1R-1Rln(RS0)-lnI(1R)I01-1R-1Rln(S0R)上式说明采取预防措施后,可以减少得病人数,并且I(t)的最大值小于不采取预防措施时的最大值。
令D=(S,I)0S1,0I(t)1,S+I=1是一个由S轴到I转以及直线S+I=1所围成的三角形区域。
对于方程组
(2),其轨线为:
S=0,I=I0e-AtI=0,S=S0e-DtO(0,0)点在
(2)在D上唯一的平衡点,并且点(0,0)是局部渐近稳定的,这是因为特征根-D及-A均小于零。
对于在直线S+I=1上的所有解均有d(S+I)dt=-AI-DS0所以,在D内出发的轨线不会越出区域D。
令D-=(S,I)0S1,0I(t)1,S+I=1在D-上取Dulac函数B(S,I)=1I(t),由Dulac定理3知在D-上不存在极限环,所以由D-上出发的轨线当t+时,必趋于平衡点(0,0)。
综上所述,我们得出如下定理。
定理2对于初始问题
(2),区域D-是平衡点O(0,0)的渐近稳定区域(见图2)。
定理3(阈值定理)设S(t),I(t)是初值问题
(2)的解,如果RS01,当t+时I(t)增加到达最大值I(1R),而后单调减少趋于零。
S(t)同时单调减少趋于零。
2.3在传染病流行之前,对易感染的人群进行有效的预防可以使易感染的人数下降,从而达到防止传染病流行的目的。
2.4在传染病发生之后,立即对易感染的人群进行有效的预防,同样可以使易感染的人数下降,从而减少得病人数。
此种图2情况下,如果在发病初期易感染的人数S01R,那么疾病会很快被消灭。
如果在发病初期易感染人数S01R,那么得病人数先增加,当其达到最大值I(1R)后,得病人数逐渐减少而后疾病被消灭。
此种情况下的最大值I(1R)小于不作预防时的最大值1-1R-1Rln(RS0)。
2.5经过一段时间以后,整个人群将趋于对该疾病具有免疫力。
3三室的模型由前述疾病传播的一般法则及人口守恒定律,可得到三室的数学模型。
dSdt=-KISdIdt=KIS-AIdRdt=AI(3)且S00,I00,R00方程组(3)是三维的,采取与
(1)的同样讨论方法,只须把第一个和第二个方程联立即可,同样得到与定理1相类似的阈值定理。
参考文献1J.E.Nash等.TheemploymentofunithydrographstodetermineTheflowsofIrisharterialdrainagechannels,ProcInstn1975.2张锦炎.常微分方程几何理论与分支问题.北京大学出版社,1987.3张芷芬等.微分方程定性理论.科学出版社,1985.收稿日期:
1999-01-08100JournalofMathematicalMedicineVol.12NO.21999
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 传染病 数学模型 建立 稳定性 分析 精品 文档
![提示](https://static.bdocx.com/images/bang_tan.gif)