传染病模型(微分方程)

微分方程建模(传染病模型）得求解.1、模型1：SI模型。假设:（1）时刻人群分为易感者（占总人数比例得)与已感染者（占总人数比例得）（2）每个病人每天有效接触得平均人数就是常数,称为日接触率，当健康者与病人接触时,健康者受感染成为病人。分析：根据假设,每个患者每天可以使个健康者变为病人，因为病人数为,所以每天共有个健康者变为病人.即:，且，设初始时刻病人比例为，则：,用MＡＴＬAB解此微分方程:>〉syｍｓaｂ>＞ｆ=ｄsｏｌvｅ（’Ｄy=ａ*y*(1-y）＇，'ｙ（0）=b＇，'t＇）f=1／（1—eｘp（—a*t)*(—1＋b)/b）%当时，分别在坐标系中作出得图像,坐标系中作出得图像,>〉a=0、1;〉>b=0、09;＞〉h=dｓolｖe（’Dy=ａ*y＊（1-ｙ）’，＇ｙ（0）＝b＇,＇ｔ’）h=1／（１-exp(—a＊ｔ）*（-1+ｂ）/b）〉〉f=ｓuｂs（h）f=1／（1+９1/９＊exp（－1／10＊t）)ﻩ得图像＞＞ezplｏt(f，[０,60］)＞>gｒidon＞〉figurｅ（2）〉＞fplot（’０、1*y*（１－y)'，[0，1]）>>gridon得图像模型分析:（1)当时，达到最大值,则此时病人增速最快。（2）当时,，即所有得人被传染,全部变为病人,这显然就是不符合实际得，其原因就是没有考虑到病人可以治愈,人群中得健康者只能变为病人,而病人不会变为健康者。2、模型２:SIＳ模型.假设：(１）时刻人群分为易感者(占总人数比例得）与已感染者（占总人数比例得)（2）每个病人每天有效接触得平均人数就是常数，称为日接触率,当健康者与病人接触时，健康者受感染成为病人。(3）病人每天被治愈得占病人总数得比例为，称为日治愈率,显然为这种传染病得平均传染期.则。则建立微分方程模型为：用MATLAB解此微分方程：〉>ｈ2=dsoｌve（＇Ｄy=a*y*（1－ｙ)-c*y',＇y（0)=b'，＇t'）h2=(a—c）／（ａ—exｐ（—（a－c)*t）＊（-ａ＋c+ｂ＊ａ)/b／(a-c）＊ａ+exp(—(a-c）＊t)*（－a+c+b＊a)/b/(a—ｃ)＊c)〉>prｅtty(ｈ２）/eｘp(—（ａ-c)ｔ)(-a+c+bａ）ａ(a—c）/|a—-—-－———-———－----－-——---－-—-———－-＼b（a—c）eｘp（-(ａ-c）t）（-a+c+ｂa)c＼+---—-----－－－－——－-—－－－-——--—----—｜b（a-c)/化简：即：。当(1)时,;（2）时,＞＞clｅar>〉h2=ｄｓolve（＇Dy=a*ｙ＊(1-y）－a＊y'，'y(0)＝b'，＇t’)h2=１/（a*t+1/b）即:.定义：一个传染期内每个病人有效接触得平均人数。则:，用MＡTLＡＢ作图像：令,,（)〉＞cleaｒ＞＞a=0、01；b=0、7;c=0、05;>〉ｈ2＝dsolvｅ（’Ｄy＝a*y＊(１—ｙ）-c＊y’，'y（0）＝b',’ｔ’）；>〉h22=ｓubs（h2)ｈ22=－１/25/（1/100—47/7０0＊exp（1/25＊ｔ）)＞＞ｅzplot(h22，［０，120］）〉＞grｉdonﻩ得图像令,，分别作图（）>>ａ＝0、３；b=0、7；c=０、15;＞>h2=dｓolvｅ(＇Dy=a＊y*（１－y)-c＊y',’y（0)=ｂ'，’t'）；>〉h2３=sｕbs(h2）ｈ２３=3/２0/（3/１0－３/35＊exｐ（—3/20＊ｔ））〉〉ｓubplot(2，1，1）〉>eｚplot(ｈ23，[0,25］）>>gridｏn>>ｂ=0、3;〉>h２４=subs(ｈ2）;>＞sｕbｐｌot(2，1，2）>〉ｅzｐｌot(h2４，[0，2５］）griｄon得图像（上面，下面）模型分析：（1）时，病人比例越来越少，最终趋于零，这就是因为传染期内经有效接触从而使健康者变为病人数不超过原来病人数得缘故。(2）时,病人比例增减性就是由来决定,其极限值随着得增加而增加。3、模型3：SIＲ模型.假设：（1）人群分为健康者，其比例、病人、病愈免疫得移出者。（２）病人得日接触率为，日治愈率为，传染期接触数为。则，对于病愈者而言，,设初始时刻得健康者与病人得比例为与,则建立微分方程模型为：由于此微分方程组得解析解无法求出，则转为相平面上讨论解得性质。相轨线得定义域应为:，由方程组消去并将得:用mａtｌb求解:〉>ｄsｏｌve（’Dy=１/cmａ／s-１'，’y（s0）=ｙ０’,’ｓ＇)ans=1/cma＊ｌｏｇ（ｓ）-s－1／cma＊log(ｓ0)+s0+y０〉〉pretty（ans)log（s)log（s0)－—-———-s－-—－—－—-+s0+y0cmacｍa即（相轨线）定义域内，时，分别取，,，在同一直角坐标系中作出其图像：>〉ｃｍa=1；y0=０、3；s０＝0、65;>>ｃleａｒ>＞f=dsolvｅ('Ｄy=1/cma/s—１’，’ｙ（ｓ0）=y0＇，’s＇)；＞〉ｃma=1；y0=0、3；s0=０、６5；＞〉ｆ1=sｕbs（f）;>>ezplot(ｆ1,［0，1])〉〉ｈoｌdon〉〉y0=０、4；ｓ0=0、３5;〉＞ｆ2=sｕｂs(f）；＞〉ezpｌoｔ(f2,[0,1］）＞>hoｌdoｎ＞＞y0＝０、5;s０=0、４5；＞＞f3=subs（f);＞〉eｚplｏt(ｆ3，［０，1］)>>holｄｏｎﻩSＩR模型得相轨线〉>ｙ０＝0、7；s０=0、25;>＞f4=sｕbs(ｆ）;〉>ezplｏt（f4，[0，1]）〉>holｄon>〉ezplot('１—s',［0,1］)>〉grｉdｏｎ模型分析:(1）不论初始条件，如何，病人比例越来越少，最终消失。（2)最终未被感染得健康者得比例就是,在中。令时,得单根即