传染病模型SI、SIS、SIR.doc

第5章 微分方程模型 （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([0 13 0 1.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); （2）SIS模型 利用MATLAB求解传染病模型中的SIS模型的解析解: 程序中a即λ, b即μ, y即i y=dsolve(Dy=a*(y-y^2)-b*y,y(0)=y0) y = (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模型改写为 画图：SIS模型的di/dt~i曲线（σ1） 程序中x即i，y即di/dt, λ=1, μ=0.3σ=3.33331. 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 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,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(i); SIR模型 利用MATLAB求解传染病模型中的SIR模型的数值解: 程序中a=λ=1, b=μ=0.3i(0)=0.02, s(0)=0.98 Ill.m文件中： 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 = 0 0.0200 0.9800 1.0000 0.0390 0.9525 2.0000 0.0732 0.9019 3.0000 0.1285 0.8169 4.0000 0.2033 0.6927 5.0000 0.2795 0.5438 6.0000 0.3312 0.3995 7.0000 0.3444 0.2839 8.0000 0.3247 0.2027 9.0000 0.2863 0.1493 10.0000 0.2418 0.1145 11.0000 0.1986 0.0917 12.0000 0.1599 0.0767 13.0000 0.1272 0.0665 14.0000