数学建模-传染病模型.ppt
《数学建模-传染病模型.ppt》由会员分享，可在线阅读，更多相关《数学建模-传染病模型.ppt（8页珍藏版）》请在冰豆网上搜索。
传染病模型传染病模型传染病是人类的大敌，通过疾病传播过程中若干重要因传染病是人类的大敌，通过疾病传播过程中若干重要因素之间的联系建立微分方程加以讨论，研究传染病流行的规素之间的联系建立微分方程加以讨论，研究传染病流行的规律并找出控制疾病流行的方法显然是一件十分有意义的工作。
律并找出控制疾病流行的方法显然是一件十分有意义的工作。
在本节中，我们将主要用多房室系统的观点来看待传染病的在本节中，我们将主要用多房室系统的观点来看待传染病的流行，并建立起相应的多房室模型。
流行，并建立起相应的多房室模型。
医生们发现，在一个民族或地区，当某种传染病流传时，医生们发现，在一个民族或地区，当某种传染病流传时，波及到的总人数大体上保持为一个常数。
即既非所有人都会波及到的总人数大体上保持为一个常数。
即既非所有人都会得病也非毫无规律，两次流行（同种疾病）的波及人数不会得病也非毫无规律，两次流行（同种疾病）的波及人数不会相差太大。
如何解释这一现象呢？
试用建模方法来加以证明。
相差太大。
如何解释这一现象呢？
试用建模方法来加以证明。
问题的提出：
问题的提出：
设某地区共有设某地区共有n+1人，最初时刻共有人，最初时刻共有i人得病，人得病，t时刻已时刻已感染（感染（infective）的病人数为的病人数为i（t），假定每一已感染者在单位假定每一已感染者在单位时间内将疾病传播给时间内将疾病传播给k个人（个人（k称为该疾病的传染强度），且称为该疾病的传染强度），且设此疾病既不导致死亡也不会康复设此疾病既不导致死亡也不会康复模型模型1此模型即此模型即MalthusMalthus模型，它大体上反映了传染病流行初期模型，它大体上反映了传染病流行初期的病人增长情况，在医学上有一定的参考价值，但随着时间的的病人增长情况，在医学上有一定的参考价值，但随着时间的推移，将越来越偏离实际情况。
推移，将越来越偏离实际情况。
已感染者与尚未感染者之间存在着明显的区别，有必要将已感染者与尚未感染者之间存在着明显的区别，有必要将人群划分成已感染者与尚未感染的易感染，对每一类中的个体人群划分成已感染者与尚未感染的易感染，对每一类中的个体则不加任何区分，来建立两房室系统。
则不加任何区分，来建立两房室系统。
则可导出：
则可导出：
故可得：
故可得：
（3.15）模型模型2记记t时刻的病人数与易感染人数时刻的病人数与易感染人数（susceptible）分别为分别为i（t）与与s（t），初始时刻的病人数为初始时刻的病人数为i。
根据病人不死也不会康复根据病人不死也不会康复的假设及（竞争项）统计筹算律，的假设及（竞争项）统计筹算律，其中：
其中：
解得：
解得：
（3.17）可得：
可得：
（3.16）统计结果显示，统计结果显示，（3.173.17）预报结果比预报结果比（3.153.15）更接近实际情况。
医学上称曲线更接近实际情况。
医学上称曲线为传染病为传染病曲线，并称曲线，并称最大值时刻最大值时刻t1为此传染病的流行为此传染病的流行高峰。
高峰。
令：
令：
得：
得：
此值与传染病的实际高峰期非常接近，可用作医学上的预报公式。
模型模型22仍有不足之处，它仍有不足之处，它无法解释医生们发现的现无法解释医生们发现的现象，且当时间趋与无穷时，象，且当时间趋与无穷时，模型预测最终所有人都得模型预测最终所有人都得病，与实际情况不符。
病，与实际情况不符。
为了使模型更精为了使模型更精确，有必要再将确，有必要再将人群细分，建立人群细分，建立多房室系统多房室系统infectiverecoveredsusceptiblekl（3.18）l称为传染病恢复系数求解过程如下：
求解过程如下：
对（对（3）式求导，由（）式求导，由
（1）、（）、
（2）得：
）得：
解解得：
得：
记：
记：
则：
则：
将人群划分为三类（见右图）：
易感染者、已感染将人群划分为三类（见右图）：
易感染者、已感染者和已恢复者（者和已恢复者（recovered）。
）。
分别记分别记t时刻的三类人数为时刻的三类人数为s（t）、i（t）和和r（t），则可建立下面的三房室模型：
则可建立下面的三房室模型：
模型模型3infectiverecoveredsusceptiblekl由由
（1）
（1）式可得：
式可得：
从而解得：
从而解得：
积分得：
积分得：
（3.19）不难验证，当t+时，r（t）趋向于一个常数，从而可以解释医生们发现的现象。
为揭示产生上述现象的原因（为揭示产生上述现象的原因（3.183.18）中）中的第（的第（11）式改写成：
）式改写成：
其中其中通常是一个与疾病种类有关的通常是一个与疾病种类有关的较大的常数。
较大的常数。
下面对下面对进行讨论，请参见右图进行讨论，请参见右图如果如果，则有则有，此疾病在该地区根本流行不起来。
，此疾病在该地区根本流行不起来。
如果如果，则开始时，则开始时，i（t）单增。
但在单增。
但在i（t）增加的同时，增加的同时，伴随地有伴随地有s（t）单减。
当单减。
当s（t）减少到小于等于减少到小于等于时，时，i（t）开始减小，开始减小，直至此疾病在该地区消失。
直至此疾病在该地区消失。
鉴于在本模型中的作用，鉴于在本模型中的作用，被被医生们称为此疾病在该地区医生们称为此疾病在该地区的阀值。
的阀值。
的引入解释了为什的引入解释了为什么此疾病没有波及到该地区么此疾病没有波及到该地区的所有人。
的所有人。
图3-14综上所述，模型综上所述，模型33指出了传染病的以下特征：
指出了传染病的以下特征：
（11）当当人人群群中中有有人人得得了了某某种种传传染染病病时时，此此疾疾病病并并不不一一定定流流传，仅当易受感染的人数与超过阀值时，疾病才会流传起来。
传，仅当易受感染的人数与超过阀值时，疾病才会流传起来。
（22）疾疾病病并并非非因因缺缺少少易易感感染染者者而而停停止止传传播播，相相反反，是是因因为为缺少传播者才停止传播的，否则将导致所有人得病。
缺少传播者才停止传播的，否则将导致所有人得病。
（33）种群不可能因为某种传染病而绝灭。
）种群不可能因为某种传染病而绝灭。
模型检验：
模型检验：
医疗机构一般依据医疗机构一般依据r（t）来统计疾病的波及人数来统计疾病的波及人数，从广义上，从广义上理解，理解，r（t）为为t时刻已就医而被隔离的人数，是康复还是死亡对时刻已就医而被隔离的人数，是康复还是死亡对模型并无影响。
模型并无影响。
及：
及：
注意到：
注意到：
可得可得：
（3.20）通常情况下，传染病波及的人数占总人数的百分比不会通常情况下，传染病波及的人数占总人数的百分比不会太大，故太大，故一般是小量。
利用泰勒公式展开取前三项，有：
一般是小量。
利用泰勒公式展开取前三项，有：
代入（代入（3.203.20）得近似方程：
）得近似方程：
积分得：
积分得：
其中：
其中：
这里双曲正切函数这里双曲正切函数：
而：
而：
对对r（t）求导求导：
（3.21）曲线曲线在医学上被称为疾病传染曲线。
在医学上被称为疾病传染曲线。
图图3-14（a）给出了（给出了（3.21）式曲线的图形，可用医疗单位式曲线的图形，可用医疗单位每天实际登录数进行比较拟合得最优曲线。
每天实际登录数进行比较拟合得最优曲线。
图3-14（a）图图3-14（b）记录了记录了1905年下半年至年下半年至1906年上半年印度孟买瘟年上半年印度孟买瘟疫大流行期间每周死亡人数，不难看出两者有较好的一致性。
疫大流行期间每周死亡人数，不难看出两者有较好的一致性。
升级会员 