1.实验7-1 传染病模型2

1、河北大学数学模型实验报告班里的专家十五课两班姓名张宇轩学号20151101006实验场所C1-229战斗机指导老师司建辉成绩实验项目1 .实验7-1描绘传染病模型2(SI模型) di/dtI的图表2 .实验7-2描绘传染病模型2(SI模型) It的图表3 .实验7-3描绘传染病模型3(SIS模型) di/dtI图表4 .实验7-4描绘传染病模型3(SIS模型) It的图表5 .实验7-5传染病模型4(SIR模型)1 .实验目的2 .实验要求1 .实验7-1描绘传染病模型2(SI模型) di/dtI的图表(参考教材p137-138 )传染病模型2(SI模型):di/dt=ki(1-i )、i(0

2、)=i0；其中，i(t )是第t天患者占总人数的比例。是一名患者每天有效接触的平均人数(日接触率)。i0是初始时刻(t=0)的患者的比例。设定k=0.1，绘制di/dt i的图表，求出I为什么取值时，di/dt达到最大值，绘制在图表上。 试试看创建并实现m文件。请参阅程序的执行结果(从图形窗口的菜单中选择“Edit/Copy Figure”，复制图形)提示。1 )画图表fplot函数以以下形式调用fplot(fun，lims )fun必须是m文件的函数名称或变量x的可执行字符串。如果lims取“xmin xmax”，则x轴限定于该区间。lims取“xmin xmax ymin ymax”时，y

3、轴也受到限制。这个问题可以利用fplot(0.1*x*(1-x )、0.10.03 )；2 )求最大值使用求解边界约束的非线性最小化函数fminbnd，调用形式如下x=fmin bnd (fun、x1、x2)fun必须是m文件的函数名称或变量x的可执行字符串。返回参数x为区间x1 hold on打印( 0，0.4 ， 0，0 )4 .实验7-4描绘传染病模型3(SIS模型) It的图表在matlab中创建m文件fun4.m代码如下所示：函数y=函数(x )x=d解算(dx=-0.2 * x * (x-(1-1/3) )、x(0)=0.2；TT=linspace (0，41，1001 )for

4、i=1:1001t=tt(i )xx(i)=eval(x )结束打印(TT，xx )hold on；打印( 0，40 ，1-1/3，1-1/3 ，-k )；x=d解算(dx=-0.2 * x * (x-(1-1/3) )、x(0)=0.9；TT=linspace (0，41，1001 )for i=1:1001t=tt(i )xx(i)=eval(x )结束打印(TT，xx，-r )axis ( 0，40，0，1 )；标题(图1 SI型号it曲线(=0.2、=3) )；xlabel(t (日) )；ylabel(i (患者的比例)长度(I (0)=0. 2，1-1 /，i(0)=0.9；在命令行中输入以下代码fun4；5 .实验7-5传染病模型4(SIR模型)在matlab中创建m文件fun5.m代码如下所示：函数y=函数(t，x )a=1；b=0.3；y=a*x(1)*x(2)-b*x(1)、-a*x(1)*x(2)；在命令行中输入以下代码ts=0:50；x