《传染病模型》由会员分享，可在线阅读，更多相关《传染病模型（13页珍藏版）》请在金锄头文库上搜索。
1、病毒扩散与传播的控制模型摘要随着科技的发展，病毒扩散与传播越来越受到人们的关注。本文通过建立微分方程模型，描述了病毒扩散与传播的过程，最后通过分析， 得到了控制病毒扩散与传播的方法。对问题一，我们通过分析影响变量的因素，建立微分方程模型。对问题二，我们通过把增加的影响因素加入到问题一的微分方程模型中，改 善后得到了新的微分方程模型。最后把变量代入，求解微分方程模型，得到结果。 当t=13天时，确诊患者人数达到峰值6793000人；t=150天时，确诊人数减少到 116800人。对问题三、问题四、问题五，通过把改变后的条件代入到问题二中的微分方 程模型中，可以得到其对应结果。问题三的结果，当 t
2、=13 天时，确诊人数达到 峰值 6769000 人， t=150 天时，确诊人数减少到 108400 人。问题四的结果，当 t=13天时，确诊人数达到峰值6795000人，t=150天时，确诊人数减少到116200 人。问题五的结果，当t=12天时，确诊人数达到峰值6793000人，当t=150天 时，确诊人数减少到 113500人对问题六，结合前面所得到的结果，我们分析在其它因素都不变的情况下只 改变一种因素，分析得到该种因素的灵敏度，最后得出各个因素的灵敏度。可以 得到，尽快开始隔离、治愈时间t 、t ，加强隔离强度p、减少人均日接触率12r 都可以改善病情。对问题七，在问题六的基础上，
3、可以得出相应的减轻病情的方法和建议。关键词：微分方程模型 微分方程组求解（MATLAB）一、问题的重述已知某种不完全确知的具有传染性病毒的潜伏期为dld2天，病患者的治愈 时间为d3天。该病毒可通过直接接触、口腔飞沫进行传播、扩散，该人群的人 均每天接触人数为r。为了控制病毒的扩散与传播将该人群分为五类：确诊患者、 疑似患者、治愈者、死亡和正常人，可控制参数是隔离措施强度P （潜伏期内的 患者被隔离的百分数）。要求：1. 在合理的假设下试建立该病毒扩散与传播的控制模型；2. 利用你所建立的模型针对如下数据进行模拟条件 1： dl=l, d2=ll, d3=30, r=10,条件2：已经知道的初
4、始发病人数为890、疑似患者为2000条件3：隔离措施强度p=60%条件4：患者2天后入院治疗，疑似患者2天后被隔离,试给出患者人数随时 间变化的曲线图，并明确标识图中的一些特殊点的具体数据，分析结果的合理性。3. 若将2中的条件4改为条件：患者1.5天后入院治疗，疑似患者1.5天后被隔离,模拟结果有何变化？4. 若仅将2中的条件3改为条件：隔离措施强度p=40%，模拟结果有何 变化？5. 若仅将2中的条件1改为条件：d1=1, d2=11, d3=30, r=250,模拟 结果有何变化？6. 分析问题中的参数对计算结果的敏感性。7. 针对如上数据给政府部门写一个不超过400字的建议报告。二、
5、问题的分析2.1 问题一的分析问题一的解决，在于对疑似患者、确诊患者、治愈者、正常人、死亡者的理 解，在理解的基础上，我们分析影响它们的因素有哪些，最后通过建立微分方程 模型来解释这些影响关系。2.2 问题二的分析问题二的解决基于问题一，在问题已的基础上，我们对于增加的影响因素进项分析，并且改进问题一中的微分方程模型。我们会得到考虑更加全面的微分方 程模型，最后把已知的数据带入方程组中，最后用matlab进行求解，可以得 到结果。2.3 问题三的分析问题三的解决，即把已知的条件代入问题二的微分方程组中，最后用 matlab求解，可以得到结果。2.4 问题四的分析问题四的解决，即把已知的条件代入
6、问题二的微分方程组中，最后用 matlab求解，可以得到结果。2.5 问题五的分析问题五的解决，即把已知的条件代入问题二的微分方程组中，最后用 matlab求解，可以得到结果。2.6 问题六的分析通过对问题二、问题三、问题四、问题五结果的分析，通过相互之间的对比， 我们会发现不同影响因素对病情的影响，最终得出结论。2.7 问题七的分析 通过问题六的结果，我们会得到影响病情的因素以及影响关系，通过这些关 系，提出减轻病情的建议和方法。三、模型的假设1总人数 N 不变，人群分为确诊患者、疑似患者、治愈者、死亡和正常人； 2. 疑似患者是被病毒感染，但是没有发病的人；3确诊患者自动被隔离，不具备传染
7、病毒的能力，所以人群中只有未被隔离的 疑似患者能够传染疾病；4治愈者具有了免疫能力，不会再被该病毒感染； 5平均潜伏期为1d d125. 在问题二中，假设人群总数为 1000 0000四、符号说明S疑似患者在人群中的比例I确诊患者在人群中的比例R治愈者在人群中的比例H正常人在人群中的比例D死亡人数在人群中的比例N人群的总人数（假设为1000 0000）d1d2病毒的潜伏期d3治愈时间r人群的人均日接触人数P隔离强度d4平均潜伏期五、模型的建立与求解5.1 问题一的模型建立与求解易知显然有S+I+R+H+D=1 根据已知条件以及假设，分析可得： 疑似患者在人群中的比例决定于未被隔离的疑似患者和疑
8、似患者发病的比 例有关。确诊患者在人群中的比例决定于疑似患者发病和患病者得到治愈的比例有 关。治愈者在人群中的比例决定于确诊患者被治愈的比例有关。 正常人的比例决定于未被隔离的疑似患者的比例有关。死亡者的比例则是人群中每一时刻人群原来整数减去疑似患者、确诊患者 治愈者、正常人人数的比例。有以上分析以及题目中所数据和假设中的数据可将此问题满足的微分方程 模型建立，如下：dSSN二 Nr (1 p)SH - Ndtd4dISIN = N Ndtdd43dRIN = Ndtd3dHN=Nr (1 p) SH dt到此，我们建立该病毒扩散与传播的控制模型；5.2 问题二的模型建立与求解问题二模型的建立
9、是建立在问题一的基础上，由题中给出的额外条件患者t1 天后得到治疗和疑似患者t天后被隔离，所以我们基于问题一的模型。2考虑到此事疑似患者感染人群的天数增加 t 天，所以此时人群中每天疑似2rSH患者的增加人数比问题一多N 竺 ；相反人群中正常人每天比问题一减少t2t2考虑到此时患病者t天后开始治疗，所以这时治愈时间变为t + d天，此时1 1 3人群中每天治愈者的增加人数变为N- ；相反人群中每天确诊患者的变为d +131SIN Nd d +14 3 1所以得到问题二的微分方程模型，如下：dS rSHSN 二 N+ Nr (1 p)SH Ndt td24dI S IN = N N dt d d
10、 +1dRdt4 3 1二 Nd +131 NrSHdtNr (1 p) SH t把问题二中已知的条件 1、条件 2、条件 3、条件4可得：d =1，d = 11 ，d =30，r=10，p=60% ，t =t =2 ，12312I=0.000089，S =0.00002 ， R =0000由假设可得H =0.99970将这些数据代入上述方程，可得微分方程组，如下：dS = 9 x SHdtdI S Idt 6 32dRdt32dHdt=-9 x SHI =0.000089， S =0.00002 ， 00R =0， H =0.999700用 MATLAB 求解此微分方程模型，可以得到患者人数
11、随时间变化的曲线图如下 图所示：KV:患!s若人*I債恋亠V常亡韭歯泄正无三二X I问题二的结果由于死亡者的比例几乎解决于 0，甚至是一个很小的负数，所以与横坐标轴重合)在图中标出三个特殊点，(0,0.000089)(13,0.6793)(150,116800)。分析可 得，起初即 t=0 时，人群中确诊患者的比例为0.000089，即890 人；当 t=13 时， 确诊患者人数比例达到最多0.6793，即6793000 人。并且013天，确诊患者人 数递增；13150 天，确诊患者人数递减，最后在 150 天时，确诊患者比例为 0.01168 即 116800 人。5.3 问题三的模型建立与
12、求解由已知可得，改变问题二中的条件四，t =t二1.5可得微分方程组：12dSS=10.67 x SH 一dt6dI S Idt6 31.5dR = Idt 31.5dH= 10.67 x SH dtD = 1 S I R HI =0.000089 ， S =0.00002 ， 00R =0， H =0.999700用MATLAB求解此微分方程模型，可以得到患者人数随时间变化的曲线图如下图所示：问题三的结果在图中标出俩个特殊点，(0,0.000089)(13,0.6769)(150,108400)。分析可 得，起初即 t=0 时，人群中确诊患者的比例为0.000089，即890 人；当 t=1
13、3 时， 确诊患者人数比例达到最多0.6769，即6769000 人。并且013天，确诊患者人 数递增；13150天，确诊患者人数递减，最后在150天时，确诊患者比例为 0.01084即 108400 人。5.4 问题四的模型建立与求解由已知可得，改变问题二中的条件三改为p=40%，同理可得微分方程:dS = 11X SH - Sdt6dI S Idt 6 32dRdt 32dHdt=llx SHI =0.000089， S =0.00002 ， 00R =0， H =0.9997 00用MATLAB求解此微分方程模型，可以得到患者人数随时间变化的曲线图如下图所示:燎出昔人昔mii:-和議亡能
14、圄诸正死-=0SO1WISO问题四的结果在图中标出俩个特殊点，(0,0.000089)(13,0.6795)(150,116200)。分析可 得，起初即 t=0 时，人群中确诊患者的比例为0.000089，即890 人；当 t=13 时， 确诊患者人数比例达到最多 0.668，即 6795000 人。并且 013 天，确诊患者人数 递增；13150 天，确诊患者人数递减，最后在 150 天时，确诊患者比例为 0.01162， 即 116200 人。5.5 问题五的模型建立与求解由已知可得，改变问题二中的条件一改为d 二 1, d 二 11, d 二 30, r 二 250123同理可得微分方程：dS = 275 x SH - Sdt6dI S Idt 6 32dRdt32dH = -275 x SHdtD=1-S-I -R-HI =0.000089， S =0.00002 ， 00R =0， H =0.999700用 MATLAB 求解此微分方程模型，可以得到患者人数随时间变化的曲线图如下图所示：问题五的结果在图中标出俩个特殊点，(0,0.000089)和(12,0.6793)。分析可得，
