传染病模型-SI-SIR-SIS.doc

数学模型实验实验报告10学院： 专 业： 姓 名: 学号：_ _ 实验时间：_ _ 实验地点： 一、实验项目:传染病模型求解二、实验目得与要求a、求解微分方程得解析解b、求解微分方程得数值解三、实验内容问题得描述各种传染病给人类带来得巨大得灾难，长期以来,建立传染病得数学模型来描述传染病得得传播过程，分析受感染人数得变化规律，探索制止传染病蔓延得手段等,一直就是各国有关专家与官员关注得课题.不同类型传染病有各自不同得特点,在此以一般得传播机理建立几种3模型。分别对3种建立成功得模型进行模型分析,便可以了解到该传染病在人类间传播得大概情况.模型一（SI模型）：(）模型假设1、在疾病传播期内所考察地区得总人数N不变，人群分为健康人与病人，时刻t这两类人在总人数中所占比例为s()与（t)。2、每个病人每天有效接触得平均人数就是常数a，a成为日接触率，当病人与健康者有效接触时,可使其患病。(2)建立模型根据假设,每个病人每天可使a(t）个健康人变成病人，t时刻病人数为N()，所以每天共有aNs(）i（t)个健康者被感染,即病人得增加率为: Ndi/t=aNi又因为s（t）+i（）1再记时刻t=0时病人得比例为i0则建立好得模型为: (0)i(3）模型求解 （代码、计算结果或输出结果)sys a i 0 :日接触率，i：病人比例，s：健康人比例，i0:病人比例在=0时得值idle(i=i(1i)，i(0）0，t）； ysb(i，a,i0，0、3，、2);ezpl（y,00)fgueisr2duble（i）;i0：0、01：1；y=、3i、*(-i）；lot(i，) SI模型得it曲线 S模型得/dti曲线(4）结果分析由上图可知，在i=0：1内，d/dt总就是增大得,且在i、5时，取到最大值,即在tnf时，所有人都将患病.上述模型显然不符合实际，为修正上述结果,我们重新考虑模型假设,建立SIS 模型模型二（SI模型)（1） 模型假设假设条件1、2与S模型相同；3、每天被治愈得病人数占病人总数得比例为常数u，成为日治愈率,病人治愈后成为仍可被感染得健康者。显然1就是平均传染期.（2)模型建立 病人得增加率:d/dtaNui 且 i（)+s（t）1；则有: d/t=ai(1i)-ui 在此定义k=a/b，可知k就是整个传染传染期内每个病人有效接触得平均人数,成为接触数。则建立好得模型为:(）=i0;（2） 模型求解 (代码、计算结果或输出结果） ys a u i0 a：日接触率，i：病人比例，：日治愈率,i0:病人比例在t=0时得值 ove（Diai*（1-i)u*i,i（0）=i0，t） 求用u表示得it解析式sms k k:接触数 =a/; i=dsolv(i=aii+i（-1/k),i（0)=0,t) 求用k表示得it解析式 给k、a、i0指定特殊值，作出相关图像=sub（i,，a,i0,2,0、3,0、); %得情况,以k=2为例ezplot（y，0,10)use 作i图，分析随时间得增加, 得变化 gext（/k） legen(1 本例中k=)fiure i=stdoube（i）；i=0：0、01:;y=0、*i、*-12； plot(i,) 作didi得图像 get（1-1/k,在此图中为、) lgn(=2)sus(i，k,a，i0，0、8,、，0、02）; ki=str2duble(i)；i:0、0:1; y=、3i、*i(-1/0、)； pot(,) 作di/dt-i得图像leged（k=0、8）gtext(k1） I模型得it曲线（k1） S模型得d/dti曲线（k1) SIS模型得t曲线（k1时,i(t）得增减性取决于i0得大小,但其极限值i（)-1/k随k得增加而增加；当k=1时,病人比例i（t)越来越小,最终趋于0，这就是由于传染期内经有效解除从而使健康者变为得病人数不超过原来病人数得缘故.模型三。SR模型（1）模型假设1. 总人数N不变，人群分为健康者、病人与病愈免疫得移出者三类，称SIR模型。时刻三类人在总人数N中占得比例分别记作与。2. 病人得日接触率为,日治愈率为（与SI模型相同)，传染期接触数为。(2)模型建立由假设显然有 (1）对于病愈免疫得移出者而言应有 （2）再记初始时刻得健康者与病人得比例分别就是0（s0）与i(i00)(不妨设移出者得初始值r=0），则I模型得方程可以写作 (3)（）模型求解我们无法求出解析解,先做数值计算：设，用TLA软件编程：unin y=il（t,x）a=;b=0、3;=ax（1)x（2)*x(1）， ax（）x(2）;t=0：50；x0=、02，0、8；t,=oe45（11，t，x0);t，xpl(,x（:,1）,t，x（：，),rid,pauseplo(x(：，)，x（：,)）表1 得数值计算结果t235678i（t)0、020、03900、0320、12850、23、2750、3320、34440、24s(t）0、9800、9250、01、869、600、80、399、28390、027t045i（）0、26、410、070、0223、61、070、0050、00010s(t)0、1493、140、05430、034、0080、400、03990、09、398 得图形 is图形(相轨线)（4）结果分析得图形见左图, 得图形见右图，称为相轨线,随着得增加,沿轨线自右向左运动。由上图结合表1可知,由初值增长至约时达到最大值，然后减少,则单调减少.进行相轨线分析，可得：平面称为相平面,相轨线在相平面上得定义域为 在方程(3)中消去，并注意到得定义，可得， （4) 容易求出它得解为 (5) 在定义域D内，上式表示得曲线即为相轨线1. 不论初始条件如何,病人终将消失，即 （)其证明如下,首先,由（)，而故存在；由(2），而，故存在，再由(1),对于充分大得有，这将导致，与存在相矛盾。2. 最终未被感染得健康者得比例就是，在()式中令得到,就是方程 （）在内得根。在图形上，就是相轨线与轴在内交点得横坐标.3. 若，则先增加,当时,达到最大值 （8）然后减小且趋近于0，则单调减小至。4. 若,则单调减少至0,单调减少至。如果仅当病人比例有一段增长得时期才认为传染病在蔓延,那么就是一个阈值，当（即)时传染病就会蔓延。而减小传染期接触数,即提高阈值，使得(即),传染病就不会蔓延(健康者比例得初始值就是一定得，通常可认为接近）。并且，即使,从（7）,（8)式可以瞧出，减少时，增加(通过作图分析）,降低,也控制了蔓延得程度，我们注意到,在中，人们得卫生水平越高,日接触率越小;医疗水平越高,日治愈率越大，于就是越小，所以提高卫生水平与医疗水平有助于控制传染病得蔓延。从另一方面瞧，就是传染期内一个病人传染得健康者得平均数,称为交换数,其含义就是一个病人被个健康者交换,所以当即时，必有,既然交换数不超过1，病人比例绝不会增加,传染病不会蔓延。建模所得：1. 符号变量如何使用2. 如何求微分方程得解析解与数值解3. 对符号变量方程作图时，先将其中得符号变量赋值，再将其变成数值变量，这也就是一种有效得解决方法.