传染病模型SI、SIS、SIR

1、WORD格式标准实用SI 模型利用 MATLAB求解传染病模型中的SI 模型的解析解:程序中 a 即 ,y 即 i y=dsolve(Dy=a*(y-y2),y(0)=y0)y = 1/(1-exp(-a*t)*(-1+y0)/y0)画图： SI 模型的 it曲线设 =1, i(0)=0.1 y=dsolve(Dy=y-y2,y(0)=0.1)y =1/(1+9*exp(-t) x=0:0.01:13; y=1./(1+9.*exp(-x); plot(x,y)title(SI模型的 it曲线 );xlabel(t);ylabel(i);axis(0 13 0 1.1);专业资料整理WORD格

2、式文案大全专业资料整理WORD格式标准实用画图： SI 模型的 di/dti曲线程序中 x 即 i ，y 即 di/dt, =1 x=0:0.01:1;y=x-x.*x; plot(x,y)title(SI模型的 di/dti曲线 );xlabel(i);ylabel(di/dt);专业资料整理WORD格式文案大全专业资料整理WORD格式标准实用SIS 模型利用 MATLAB求解传染病模型中的SIS 模型的解析解:程序中 a 即 ,b 即 ,y 即 i y=dsolve(Dy=a*(y-y2)-b*y,y(0)=y0)y =(a-b)/(a-exp(-(a-b)*t)*(-a+b+y0*a)/

3、y0/(a-b)*a+exp(-(a-b)*t)*(-a+b+y0*a)/y0/(a-b)*b)画图： SIS 模型的 di/dti曲线 1程序中 x 即 i ，y 即 di/dt, =1, =0.3 x=0:0.01:1; y=0.7.*x-x.2; plot(x,y)title(SIS模型的 di/dti曲线 );xlabel(i);ylabel(di/dt);专业资料整理WORD格式文案大全专业资料整理WORD格式标准实用画图： SIS 模型的 it曲线 1设 =1, =0.3,i(0)=0.02 y=dsolve(Dy=0.7*y-y2,y(0)=0.02)y = 7/(10+340*

4、exp(-7/10*t) x=0:1:16; y=7./(10+340.*exp(-7./10.*x); plot(x,y)title(SIS模型的 it曲线 );xlabel(t);ylabel(i);专业资料整理WORD格式文案大全专业资料整理WORD格式标准实用画图： SIS 模型的 di/dti曲线 1程序中 x 即 i ，y 即 di/dt, =0.5, =0.6 x=0:0.01:1; y=-0.5.*x.2-0.1.*x; plot(x,y)title(SIS模型的 di/dti曲线 );xlabel(i);ylabel(di/dt);专业资料整理WORD格式文案大全专业资料整理

5、WORD格式标准实用画图： SIS 模型的 it曲线 1设 =0.5, =0.6,i(0)=0.02 y=dsolve(Dy=-0.5*y2-0.1*y,y(0)=0.02)y = 1/(-5+55*exp(1/10*t) x=0:1:40; y=1./(-5+55.*exp(1./10.*x); plot(x,y)title(SIS模型的 it曲线 );xlabel(t);ylabel(i);专业资料整理WORD格式文案大全专业资料整理WORD格式标准实用SIR 模型利用 MATLAB求解传染病模型中的SIR 模型的数值解:程序中 a= =1, b= =0.3 ， i(0)=0.02,s(0

6、)=0.98M文件中：function y=ill(t,x)a=1;b=0.3;y=a*x(1)*x(2)-b*x(1),-a*x(1)*x(2);命令窗口中： t,x=ode45(ill,0:50,0.02,0.98);t,x ans =00.02000.98001.00000.03900.95252.00000.07320.90193.00000.12850.81694.00000.20330.69275.00000.27950.54386.00000.33120.39957.00000.34440.28398.00000.32470.20279.00000.28630.149310.00

7、000.24180.114511.00000.19860.091712.00000.15990.076713.00000.12720.066514.00000.10040.059315.00000.07870.054316.00000.06140.050717.00000.04780.048018.00000.03710.046019.00000.02870.044520.00000.02230.043421.00000.01720.042622.00000.01330.041923.00000.01030.041524.00000.00790.041125.00000.00610.04082

8、6.00000.00470.040627.00000.00360.040428.00000.00280.040329.00000.00220.040230.00000.00170.0401专业资料整理WORD格式文案大全专业资料整理WORD格式标准实用31.00000.00130.040032.00000.00100.040033.00000.00080.040034.00000.00060.039935.00000.00050.039936.00000.00040.039937.00000.00030.039938.00000.00020.039939.00000.00020.039940.

9、00000.00010.039941.00000.00010.039942.00000.00010.039943.00000.00010.039944.00000.00000.039845.00000.00000.039846.00000.00000.039847.00000.00000.039848.00000.00000.039849.00000.00000.039850.00000.00000.0398 plot(t,x(:,1),t,x(:,2),grid,pausei(t),s(t)图形如下：专业资料整理WORD格式文案大全专业资料整理WORD格式标准实用 plot(x(:,2),x(:,1),grid,pause is 图形相轨线如下：专业资料整理WORD格式文案大全专业资料整理WORD格式标准实用画图： SIR 模型的相轨线程序中 y 即 i, x即 s, =1, =0.3 s(0)=0.32; s(0)=0.58; s(0)=0.73; s(0)=0.85 x=0:0.01:1; y=1-x; y1=1-x+0.3.*(log(x)-log(0.32); y2=1-x+0.3.*(log(x)-log(0.58); y3=1-x+0.3.*(log(