matlab傳染病模型

1、實徐目的：理解傳搽病的四類棧型，學會利用Matlab 件求解微分方程(組丿。賣絵題目：利用Matlab求解傳染病的SIS微分方程模型，并繪制教材P139頁圖3-圖6。SISM假設：(1) , F肘刻人群分為易威者(占總?cè)藬?shù)比例的$(/)和己威染者(占總?cè)藬?shù)比例的 /(/)(2) 、每個病人每夭有數(shù)接紐的平均人數(shù)是常數(shù)兄，幾稱為目接融率，當健康者與病人接紘肘，健康者受威染成為病人。(3) ,病人每天彼治念的&病人總數(shù)的比例為“，稱為可焙念率，顯然.丄為“這種傳捺病的平均傳染期。則建立微分方程棧型為：務=丸(1_但,?：() = %令療=令則模型可寫作分別作圖：當 sigma 1 肘Stepl:先

2、主義函救function y=prl(ijambda,sigma)y=-lambda.*i/f(i-(1-L/sigma)$tep2:作圖lambda=0.3;$igma=2: i=0:001:1：y=prl (ijambda.sigma) plot(Ly)0.020-0.02-0.04-0.06-0.08-0.1-0.12-0.14-0.1600.10.20.30.40.50.60.70.80.91生 sigma 1 肘Step*先定義函救function di=crb(t,ijambda,sigma)di=-lambda*!* (i- (1 -1/sigma)%$tep2:求解并作圖clc

3、clearlambda = 0.01:sigma=2;t JI = ode45(crb,0 J000,0.9.( Jambda,sigma);tJ2 = ode45(crb40J000.0.2,Jambda.$igma);plot(tJLtJ2jJ/2)legend (sigmar)0.9 -0.8 -0.7 -0.6 -0.5 _二一p0.4 -0.3 -0 2 11ccc01002003004005006007008009001000當 sigma = 1 肘Step*先定義函救function di=crb(t,ijambda,sigma)di=-lambda*!* (i- (1 -1/

4、sigma)%$tep2:求解并作圖clcclearlambda = 0.01:sigma = 1t JI = ode45(crb,0 J000,0.9.( Jambda,sigma);plot(tjl)lege nd(sigma =)0.9 - 0.8 -0.7 -0.6 0.50.40.30.20.1Q Lc c l c c01002003004005006007008009001000當 sigma 1 肘Step*先定義函救function di=crb(t,ijambda,sigma)di=-lambda*!* (i- (1 -1/sigma)%$tep2:求解并作圖clcclear

5、lambda = 0.01:$igma = 0.5;t JI = ode45(crb,0J000,0.9.Jambda,sigma); plot(tjl)legend (sigmar)a10.8 -0.6.0 2 Lttcc01002003004005006007008009001000創(chuàng)用matlab求鮮上面SIS模憂.提示(畫圖5程序)：Stepl:先定義函救function y=prl0Jambda,sigma)y=-lambda.*i/(i-(1-l./sigma)$tep2:作圖lambda = 0.3：sigma=2:i=0:0.01:l：y=prl (ijambda,sigma)plot(i.y)提示(畫圖6程序)：Stepl:丸岌義函救function di=crb(tjjambda,sigma)di=-lambda*i* (i- (1 -1/sigma)%step2:求解并作圖clcclearlambda = 0.01:sigma=0.2;tj = ode45(crb,0J00.0.94)Jambda,sigma); plot(tj)legend(sigmar)實驗指導書：利用matbb求下面定解問題、并作圖dy= 2y(3-y)