数模实验第四版传染病模型求解

《数模实验第四版传染病模型求解》由会员分享，可在线阅读，更多相关《数模实验第四版传染病模型求解（8页珍藏版）》请在装配图网上搜索。

1、数学模型实验实验报告3学院： 河北大学工商学院 专 业： 电气七班 姓 名： 李青青 学号： 2012484098 实验时间： 2014/4/15 实验地点： B3-301 一、实验项目：传染病模型求解二、实验目的和要求a. 求解微分方程的解析解b. 求解微分方程的数值解三、实验内容1.用Matlab求解SI和SIS模型的解析解。 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/(

2、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);当a=2，同理得图a代表到达最大点的程度当y0=1时得由图知y代表变化曲线 SIS模型程序中程序中a=，b即，y=i(1)得出=0.51,=0.6，i(0)=0.02 (1与 1的图形方向相反。 2. 用Matlab求解SIR模型的数值解。 SIR模型程序中a=1,b=0.3，i(0)=0.02,s(0)=0.98.m文件中命令行中：绘制健康者,病人和移出者的曲线：取定参数值,绘制健康者,病人和移出者的曲线以及健康者-病人的相轨线.参考教材内容P140.