文档介绍：该【传染病模型 】是由【泰山小桥流水】上传分享，文档一共【8】页，该文档可以免费在线阅读，需要了解更多关于【传染病模型 】的内容，可以使用淘豆网的站内搜索功能，选择自己适合的文档，以下文字是截取该文章内的部分文字，如需要获得完整电子版，请下载此文档到您的设备，方便您编辑和打印。传生病模型(微分方程)传生病模型(微分方程)1/8传生病模型(微分方程).微分方程建模(传生病模型)的求解。1、模型1:SI模型。假定:(1)t时刻人群分为易感者(占总人数比率的s(t))和已感染者(占总人数比率的y(t))(2)每个病人每日有效接触的平均人数是常数,称为日接触率,当健康者与病人接触时,健康者受感染成为病人。剖析:依照假定,每个患者每日能够使s(t)个健康者变为病人,由于病人数为Ny(t),所以每日共有Ns(t)y(t)个健康者变为病人。即:NdyNsy,且s(t)y(t)1,设初始时刻病人比率为b,则:dtdyy(1y),用MATLAB解此微分方程:dty(0)bsymsabf=dsolve('Dy=a*y*(1-y)','y(0)=b','t')f=1/(1-exp(-a*t)*(-1+b)/b)%y(t)11b1et1(11),,分别在坐标系oty中作出y(t)的图像,坐标系oyy中作出yy(1y)的图像,1/(1+91/9exp(-1/10t))>>a=;>>b=;>>h=dsolve('Dy=a*y*(1-y)','y(0)=b','t')=/(1-exp(-a*t)*(-1+b)/b)>>f=subs(h)=0t01020405060301/(1+91/9*exp(-1/10*t))y(t)的图像>>ezplot(f,[0,60])>>>>figure(2)>>fplot('*y*(1-y)',[0,1])(微分方程)传生病模型(微分方程)2/8传生病模型(微分方程)>>gridon模型剖析:(1)当(2)当yy(1y)的图像1时,dy达到最大值,则此时病人增速最快。2dt时,y(t)1,即全部的人被传染,全部变为病人,这显然是不传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)符合实质的,其原因是没有考虑到病人能够治愈,人群中的健康者只能变为病人,而病人不会变为健康者。2、模型2:SIS模型。假定:(1)t时刻人群分为易感者(占总人数比率的s(t))和已感染者(占总人数比率的y(t))(2)每个病人每日有效接触的平均人数是常数,称为日接触率,当健康者与病人接触时,健康者受感染成为病人。(3)病人每日被治愈的占病人总数的比率为,称为日治愈率,显然1为这种传生病的平均传染期。则NdyNsyNy。则成立微分方程模型为:传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)dydtdty(1y)y传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)y(0)b用MATLAB解此微分方程:h2=dsolve('Dy=a*y*(1-y)-c*y','y(0)=b','t')h2=(a-c)/(a-exp(-(a-c)*t)*(-a+c+b*a)/b/(a-c)*a+exp(-(a-c)*t)*(-a+c+b*a)/b/(a-c)*c)pretty(h2)exp(-(a-c)t)(-a+c+ba)a(a-c)/|a---------------------------------\b(a-c)exp(-(a-c)t)(-a+c+ba)c+--------------------------------|b(a-c)/传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)优选传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程).化简:ace(ac)t.(acba)ae(ac)t.(acba)cab(ac)b(ac)b(ac)2ab(ac)e(ac)t.(acba)e(ac)t.(acba)cb(ac)2ab(ac)(ca)e(ac)t.(acba)ab(ac)(ca)e(ac)t.(acba)1b(ac)21(1a)e(ac)tacbac(11即:y(t))e()t。b(11当(1)时,y(t))e()t;b(2)时,clearh2=dsolve('Dy=a*y*(1-y)-a*y','y(0)=b','t')h2=1/(a*t+1/b)传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)1即:y(t)tb1。传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)定义:一个传染期内每个病人有效接触的平均人数。传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)1则:y()1,(1),用MATLAB作图像:0(1),,()=;b=;c=;h2=dsolve('Dy=a*y*(1-y)-c*y','y(0)=b','t');h22=subs(h2)优选-1/25/(1/100-47/700exp(1/25t))(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)t传生病模型(微分方程)传生病模型(微分方程)12/8传生病模型(微分方程).h22=-1/25/(1/100-47/700*exp(1/25*t))ezplot(h22,[0,120]),,(a=;b=;c=;h2=dsolve('Dy=a*y*(1-y)-c*y','y(0)=b','t');h23=subs(h2)h23=3/20/(3/10-3/35*exp(-3/20*t))>>subplot(2,1,1)>>ezplot(h23,[0,25])>>>>b=;>>h24=subs(h2);>>subplot(2,1,2)>>ezplot(h24,[0,25]))3/20/(3/10-3/35exp(-3/20t))10152025t3/20/(3/10+1/5exp(-3/20t))10152025t传生病模型(微分方程)传生病模型(微分方程)13/8传生病模型(微分方程)1的图像(,)模型剖析:(1)1时,病人比率越来越少,最后趋于零,这是由于传染期内经有效接触进而使健康者变为病人数不高出原来病人数的缘故。(2)11时,病人比率y(t)增减性是由b来决定,其极限值y()1随着的增加而增加。传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)优选传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程).3、模型3:SIR模型。假定:(1)人群分为健康者,其比率s(t)、病人i(t)、病愈免疫的移出者r(t)。(2)病人的日接触率为,日治愈率为,传染期接触数为。则s(t)i(t)r(t)1,dr关于病愈者而言,y,dt设初始时刻的健康者和病人的比率为s0和y0,则成立微分方程模型为:dysyy传生病模型(微分方程)传生病模型(微分方程)16/8传生病模型(微分方程)dtdsdtsy传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)y(0)y0,s(0)s0由于此微分方程组的剖析解无法求出,则转为相平面sy上讨论解的性质。相轨线的定义域(s,y)D应为:D(s,y)s0,y0,sy1,由方程组消去dt并将得:传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)dydsyss0用matlb11sy0求解:传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)dsolve('Dy=1/cma/s-1','y(s0)=y0','s')ans=1/cma*log(s)-s-1/cma*log(s0)+s0+y0pretty(ans)log(s)log(s0)-------s--------+s0+y0cmacma即y(s)(s0y0)s1lns(相轨线)s0定义域内,1时,(y0,s0)分别取(,),(,),(,),(,)传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)优选传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程).在同素来角坐标系中作出其图像:cma=1;y0=;s0=;clearf=dsolve('Dy=1/cma/s-1','y(s0)=y0','s');cma=1;y0=;s0=;>>f1=subs(f);1-s>>ezplot(f1,[0,1])>>>>y0=;s0=;>>f2=subs(f);>>ezplot(f2,[0,1])>>>>y0=;s0=;>>f3=subs(f);0>>ezplot(f3,[0,1])>>holdonSIR模型的相轨线y0=;s0=;f4=subs(f);ezplot(f4,[0,1])holdonezplot('1-s',[0,1])gridon模型剖析:(1)无论初始条件s0,y0怎样,病人比率越来越少,最后消失。(2)最后未被感染的健康者的比率是s,在y(s)(s0y0)1ssln中。s0令y(s)0时,(s0y0)s1lns0的单根即为s:最后未被感染的健康者的比s0例。在图像上:相轨线与s轴在(0,1)内交点的横坐标。(3)当s01时传生病不会延伸,(如最左边的曲线,随着t的增加,病人数y(t)在减小)。所以提高医疗卫生水平(使日接触率减小或许使日治愈率增大),进而使1变大,也可降低s0(设y00,则s0r01),则s01r01,r011,即便免疫者比率增大。这实质上是比较困难的。如5,r080%。如有侵权请联系见告删除,感谢你们的配合!传生病模型(微分方程)传生病模型(微分方程)22/8传生病模型(微分方程)优选传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程).传生病模型(微分方程)传生病模型(微分方程)24/8传生病模型(微分方程)优选传生病模型(微分方程)传生病模型(微分方程)8/8传生病模型(微分方程)