资源描述：
《传染病模型si、sis、sir》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。
1、标准实用SI模型利用MATLAB求解传染病模型中的SI模型的解析解:程序中a即λ,y即i>>y=dsolve('Dy=a*(y-y^2)','y(0)=y0')y=1/(1-exp(-a*t)*(-1+y0)/y0)画图：SI模型的i~t曲线设λ=1,i(0)=0.1>>y=dsolve('Dy=y-y^2','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模型的i~t曲线');xlabel('t');yl
2、abel('i');axis([01301.1]);文案大全标准实用画图：SI模型的di/dt~i曲线程序中x即i，y即di/dt,λ=1>>x=0:0.01:1;y=x-x.*x;>>plot(x,y)title('SI模型的di/dt~i曲线');xlabel('i');ylabel('di/dt');>>文案大全标准实用SIS模型利用MATLAB求解传染病模型中的SIS模型的解析解:程序中a即λ,b即μ,y即i>>y=dsolve('Dy=a*(y-y^2)-b*y','y(0)=y0')y=(a-b)/(a-ex
3、p(-(a-b)*t)*(-a+b+y0*a)/y0/(a-b)*a+exp(-(a-b)*t)*(-a+b+y0*a)/y0/(a-b)*b)画图：SIS模型的di/dt~i曲线（δ>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/dt~i曲线');xlabel('i');ylabel('di/dt');>>文案大全标准实用画图：SIS模型的i~t曲线（δ>1）设λ=1,μ=0.3,i(0)=0.02>>
4、y=dsolve('Dy=0.7*y-y^2','y(0)=0.02')y=7/(10+340*exp(-7/10*t))>>x=0:1:16;>>y=7./(10+340.*exp(-7./10.*x));>>plot(x,y)title('SIS模型的i~t曲线');xlabel('t');ylabel('i');>>画图：SIS模型的di/dt~i曲线（δ≤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,
5、y)title('SIS模型的di/dt~i曲线');xlabel('i');ylabel('di/dt');>>画图：SIS模型的i~t曲线（δ≤1）设λ=0.5,μ=0.6,i(0)=0.02文案大全标准实用>>y=dsolve('Dy=-0.5*y^2-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模型的i~t曲线');xlabel('t');ylabel
6、('i');>>SIR模型文案大全标准实用利用MATLAB求解传染病模型中的SIR模型的数值解:程序中a=λ=1,b=μ=0.3，i(0)=0.02,s(0)=0.98M文件中：functiony=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.000
7、00.12850.81694.00000.20330.69275.00000.27950.54386.00000.33120.39957.00000.34440.28398.00000.32470.20279.00000.28630.149310.00000.24180.114511.00000.19860.091712.00000.15990.076713.00000.12720.066514.00000.10040.059315.00000.07870.054316.00000.06140.050717.00000.
8、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.0
