浙江大学远程教育2015年绩效管理在线作业

1传染病模型摘要当今社会，人们开始意识到通过定量地研究传染病的传播规律，建立传染病的传播模型，可以为预测和控制传染病提供可靠、足够的信息。本文利用微分方程稳定性理论对传统传染病动力学建模方式进行综述，且针对甲流，SARS等新生传染病模型进行建模和分析。不同类型的传染病的传播过程有其各自不同的特点，我们不是从医学的角度一一分析各种传染病的传播，而是从一般的传播机理分析建立各种模型，如简单模型，SI模型，SIS模型，SIR模型等。本文中，我们应用传染病动力学模型来描述疾病发展变化的过程和传播规律，运用联立微分方程组体现疫情发展过程中各类人的内在因果联系，并在此基础上建立方程求解算法。然后，通过借助MATLAB程序拟合出与实际较为符合的曲线并进行了疫情预测，评估各种控制措施的效果，从而不断完善文中的模型。本文由简到难、全面地评价了该模型的合理性与实用性，而后对模型和数据也做了较为扼要的分析，进一步改进了模型的不妥之处。同时，在对问题进行较为全面评价的基础上又引入更为全面合理的假设，运用双线性函数模型对卫生部的措施进行了评价并给出建议，做好模型的完善与优化工作。关键词传染病模型,简单模型，SI，SIS,SIR，微分方程，MATLAB。2一、问题重述有一种传染病（如SARS、甲型H1N1）正在流行，现在希望建立适当的数学模型，利用已经掌握的一些数据资料对该传染病进行有效地研究，以期对其传播蔓延进行必要的控制，减少人民生命财产的损失。考虑如下的几个问题，建立适当的数学模型，并进行一定的比较分析和评价展望。1、不考虑环境的限制，设单位时间内感染人数的增长率是常数，建立模型求T时刻的感染人数。2、假设单位时间内感染人数的增长率是感染人数的线性函数，最大感染时的增长率为零。建立模型求T时刻的感染人数。3、假设总人口可分为传染病患者和易感染者，易感染者因与患病者接触而得病，而患病者会因治愈而减少且对该传染病具有很强的免疫功能，建立模型分析T时刻患病者与易感染者的关系，并对传染情况（如流行趋势，是否最终消灭）进行预测。二、问题分析1、这是一个涉及传染病传播情况的实际问题，其中涉及传染病感染人数随时间的变化情况及一些初始资料，可通过建立相应的微分方程模型加以解决。2、问题表述中已给出了各子问题的一些相应的假设。3、在实际中，感染人数是离散变量，不具有连续可微性，不利于建立微分方程模型。但由于短时间内改变的是少数人口，这种变化与整体人口相比是微小的。因此，为了利用数学工具建立微分方程模型，我们还需要一个基本假设感染人数是时间的连续可微函数。3三、模型假设模型二和模型三的假设条件假设一在疾病传播期内所考察地区的总人数N不变，即不考虑生死，也不考虑迁移。人群分为易感染者（SUSCEPTIBLE）和已感染者（INFECTIVE）两类（取两个词的第一个字母，称之为SI模型），以下简称健康者和病人。时刻T这两类人在总人数中所占比例分别记作ST和IT。假设二每个病人每天有效接触的平均人数是常数，称为日接触率。当病人与健康者接触时，使健康者受感染变为病人。假设三模型三在假设一和假设二的基础上进行考虑，然后设病人每天治愈的比例为，称为日治愈率。病人治愈后成为仍可被感染的健康者，显然1/是这种传染病的平均传染期。模型四的假设条件假设四总人数N不变。人群分为健康者、病人和病愈免疫的移出者（REMOVED）三类，称SIR模型。三类人在总数N中占的比例分别记作ST,IT和RT。假设五病人的日接触率为，日治愈率为（与SI模型相同），传染期接触为/。4四、符号说明T某一具体时刻XT病人人数每天每个病人有效接触的人数N总人数ST健康者总人数IT病人总人数I初始时0刻病人的比例T病人的最大M值日治愈率1/平均传染率接触率RT移出者S初始时刻0健康者的比例56五、模型的建立与求解模型1在这个最简单的模型中，设时刻T的病人人数XT是连续、可微函数，并且每天每个病人有效接触（足以使人致病的接触）的人数为常数，考察T到病人人数的增加，就有TXTTX程有个病人，即得微分方时有再设0XT10,DXT方程（1）的解为20TETX结果表明，随着T的增加，病人人数XT无限增长，这显然是不符合实际的。建模失败的原因在于在病人有效接触的人群中，有健康人也有病人，而其中只有健康人才可以被传染为病人，所以在改进的模型中必须区别这两种人。模型2（SI模型）的增加率，即就是病人数个健康者被感染，于是有，所以每天共为病人数为个健康者变为病人，因天可使根据假设，每个病人每NISITINSTNIT3DSITIN又因为41TIS7，则病人的比例为再记初始时刻00IT5,1D0ITI方程（5）是LOGISTIC模型。它的解为610TEI所示。和图的图形如图和2DITTI，这个时刻为到达最大值时可知，第一，当式及图由MTITIIDD2/116,57LN01ITM这时病人增加的最快，可以认为是医院的门诊量最大的一天，预示着传染病高潮的到来，是医疗卫生部门关注的时刻。况。，这显然不符合实际情将被传染，全变为病人即所有人终时到来。第二，当可以推迟传染病高潮的健设施、提高卫生水平以改善保越小卫生水平越高。所，表示该地区的卫生水平成反比，因为日接触率与,1ITTM其原因是模型中没有考虑到病人可以治愈，人群中的健康者只能变成病人，病人不会再变成健康者。8模型3（SIS模型）有些传染病如伤风、痢疾等愈合后免疫力很低，可以假定无免疫性，于是病人被治愈后变成健康者，健康者还可以被感染再变成病人，所以这个模型成为SIS模型。考虑到这一模型的假设条件，于是有（8）NITTIITITNS可得微分方程0（9）II1IDT定义（10）其中是整个传染期内每个病人有效接触的平均人数，称为接触数。得到（11）1/DIIT模型4（SIR模型）大多数传染者如天花流感肝炎麻疹等治愈后均有很强的免疫力，所以病愈的人既非健康者（易感染者），也非病人（已感染者），因此他们将被移除传染系统，我们称之为移除者，记为R类。SIR模型是指易感染者被传染后变为感染住，感病者可以被治愈，并会产生免疫力，变为移除者。人员流动图为SIR。1模型构成在假设1中显然有9STITRT1（12）对于病愈免疫的移出者的数量应为（13）RTDNI不妨设初始时刻的易感染者、染病者、恢复者的比例分别为（0），（S0I0），0，则SIR基础模型用微分方程组表示如下I0R（14）DITSDRTIIST，IT的求解极度困难，在此我们先做数值计算来预估计ST，IT的一般变化规律。2数值计算在方程（3）中设1，03，I（0）002，S（0）098，用MATLAB软件编程FUNCTIONYILLT,XA1B03YAX1X2BX1AX1X2TS050X0020,098T,XODE45ILL,TS,X0PLOTT,X,1,T,X,2PAUSE10PLOTX,2,X,1输出的简明计算结果列入表1。IT,ST的图形以下两个图形，IS图形称为相轨线,初值I0002,S0098相当于图2中的P0点,随着T的增,S,I沿轨线自右向左运动由表1、图1、图2可以看出,IT由初值增长至约T7时达到最大值,然后减少,T,I0,ST则单调减少,T,S00398并分析IT,ST的一般变化规律表1IT,ST的数值计算结果T012345678IT002000039000732012850203302795033120344403247ST098000952509019081690692705438039950283902027T91015202530354045IT02863024180078700223000610001700005000010ST0149301145005430043400408004010039900399003981113相轨线分析我们在数值计算和图形观察的基础上，利用相轨线讨论解I（T）,S（T）的性质。D（S，I）|S0，I0，SI1（15）在方程（14）中消去并注意到的定义，可得TD（16）1IS0|SI所以1ISDD00I1SSD利用积分特性容易求出方程5的解为（17）001LNSISI12在定义域D内,17式表示的曲线即为相轨线,如图3所示其中箭头表示了随着时间T的增加ST和IT的变化趋向。图3下面根据（14），（17）式和图3分析ST,IT和RT的变化情况T时它们的极限值分别记作,和SIR1不论初始条件S0,I0如何,病人将消失,即0,TI2最终未被感染的健康者的比例是,在7式中令I0得到,是方001LNSSI在0,1/内的根在图形上是相轨线与S轴在0,1/内交点的横坐标3若1/,则开始有，IT先增加,令0，可0S1ISDO1ISD得当S1/时,IT达到最大值001LNMISIS（13然后S1/即1/S0时传染病就会蔓延而减小传染0S期接触数,即提高阈值1/使得1/即1/,传染病就不会蔓延0S0S健康者比例的初始值是一定的,通常可认为接近1。0S0S并且,即使1/,减小时,增加通过作图分析,降低,也控制MI了蔓延的程度我们注意到在中,人们的卫生水平越高,日接触率越小医疗水平越高,日治愈率越大,于是越小,所以提高卫生水平和医疗水平有助于控制传染病的蔓延从另一方面看,是传染期内一个病人传染的健康者的平均数,1/S称为交换数,其含义是一病人被个健康者交换所以当即时01/S01S必有既然交换数不超过1,病人比例IT绝不会增加,传染病不会蔓延。5群体免疫和预防根据对SIR模型的分析,当时传染病不会蔓延所以为制止蔓延,01/S除了提高卫生和医疗水平,使阈值1/变大以外,另一个途径是降低,这可以0S通过比如预防接种使群体免疫的办法做到。14忽略病人比例的初始值有,于是传染病不会蔓延的条件0I01SR01/S可以表为0R这就是说，只要通过群体免疫使初始时刻的移出者比例（即免疫比例）就可以制止传染病的蔓延。这种办法生效的前提条件是免疫者要均匀分布在全体人口中，实际上这是很难做到的。据估计当时印度等国天花传染病的接触数5，至少要有80的人接受免疫才行。据世界卫生组织报告，即使花费大量资金提高，也因很难0R做到免疫者的均匀分布，使得天花直到1977年才在全世界根除。而有些传染病的更高，根除就更加困难。6模型验证上世纪初在印度孟买发生的一次瘟疫中几乎所有病人都死亡了。死亡相当于移出传染系统，有关部门记录了每天移出者的人数，即有了的实际数据，RTDKERMACK等人用这组数据对SIR模型作了验证。首先，由方程（12），（14）可以得到SRTDDISIT，两边积分得1SRT上式两边同时乘以可00SRD0LN|SR0RSE所以80RTSTE15再901RRTDIRSSE当时，取（13）式右端TAYLOR展开式的前3项得1/RRE102001RTSDS在初始值0下解高阶常微分方程得0R110212TTSH其中，从而容易由（10）式得出200SI01T20RTDTSCH然后取定参数S0,等，画出（11）式的图形，如图4中的曲线，实际数据在图中用圆点表示，可以看出，理论曲线与实际数据吻合得相当不错。16六、模型评价与推广根据传染病的模型建立研究进而推广产生了传染病动力学模型。传染病动力学1是对进行理论性定量研究的一种重要方法,是根据种群生长的特性,疾病的发生及在种群内的传播,发展规律,以及与之有关的社会等因素,建立能反映传染病动力学特性的数学模型,通过对模型动力学性态的定性,定量分析和数值模拟,来分析疾病的发展过程,揭示流行规律,预测变化趋势,分析疾病流行的原因和关键。对于2003年发生的SARS疫情,国内外学者建立了大量的动力学模型研究其传播规律和趋势,研究各种隔离预防措施的强度对控制流行的作用,为决策部门提供参考有关SARS传播动力学研究多数采用