资源描述：
《传染病模型SI、SIS、SIR讲解学习.doc》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
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');ylabel('i');axis([01
2、301.1]);仅供学习与交流，如有侵权请联系网站删除谢谢11精品好文档，推荐学习交流画图：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');>>仅供学习与交流，如有侵权请联系网站删除谢谢11精品好文档，推荐学习交流SIS模型利用MATLAB求解传染病模型中的SIS模型的解析解:程序中a即λ,b即μ,y即i>>y=dsolve('Dy=a*(y-y^2)-b*y','y(0)=y0')y=(
3、a-b)/(a-exp(-(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');>>仅供学习与交流，如有侵权请联系网站删除谢谢11精品好文档，推荐学习交流画图：SIS模型的i~t曲线（δ>1）设λ=1,μ=0.3,
4、i(0)=0.02>>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）仅供学习与交流，如有侵权请联系网站删除谢谢11精品好文档，推荐学习交流程序中x即i，y即di/dt,λ=0.5,μ=0.6>>x=0:0.01:1;>>y=-0.5.*x.^2-0.1.
5、*x;>>plot(x,y)title('SIS模型的di/dt~i曲线');xlabel('i');ylabel('di/dt');>>画图：SIS模型的i~t曲线（δ≤1）设λ=0.5,μ=0.6,i(0)=0.02仅供学习与交流，如有侵权请联系网站删除谢谢11精品好文档，推荐学习交流>>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曲线');xlabe
6、l('t');ylabel('i');>>SIR模型仅供学习与交流，如有侵权请联系网站删除谢谢11精品好文档，推荐学习交流利用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.0
7、7320.90193.00000.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.05
