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

1、河北大学数学模型实验 实验报告班级专业15计科2班姓名张宇轩学号20151101006实验地点C1-229指导老师司建辉成绩实验项目1 .实验7-1传染病模型2 （ SI模型） 画di/dt i曲线图2 .实验7-2传染病模型2 （ SI模型）一一画it曲线图3 .实验7-3传染病模型3 （ SIS模型） 画di/dt i曲线图4 .实验7-4传染病模型3 （ SIS模型） 一一画it曲线图5 .实验7-5传染病模型4 （ SIR模型）、实验目的二、实验要求1.实验7-1传染病模型2 （ SI模型）画di/dti曲线图（参考教材 p137-138 ）传染病模型2 （ SI模型）：di/dt=k

2、i（1-i）,i（0）=i0；其中，i是第t大病人在总人数中所占的比例。人是每个病人每天有效接触的平均人数（日接触率）。i0是初始时刻（t=0 ）病人的比例。取k=0.1 ,画出di/dti曲线图，求i为何值时di/dt达到最大值，并在曲线图上标注。试编写一个m文件来实现。参考程序运行结果（在图形窗口菜单选择Edit/Copy Figure ,复制图形）：提示1)画曲线图用fplot函数，调用格式如下：fplot(fun,lims)fun必须为一个M文件的函数名或对变量 x的可执行字符串。若lims取xmin xmax,则x轴被限制在此区间上。若 lims 取xmin xmax ymin ym

3、ax ,贝U y 轴也被限制。本题可用fplot('0.1*x*(1-x)',0 1.1 0 0.03);2)求最大值用求解边界约束条件下的非线性最小化函数fminbnd ,调用格式如下:x=fminbnd( 'fun ' ,x1,x2)fun必须为一个M文件的函数名或对变量 x的可执行字符串。返回自变量x在区间x1<x<x2 上函数取最小值时的x值。本题可用x=fminbnd('-0.1*x*(1-x)',0,1)y=0.1*x*(1-x)4)指示最大值坐标用线性绘图函数plot ,调用格式如下：plot(x1,y1,'颜色

4、 线型 数据点图标，x2,y2,'颜色 线型 数据点图标，)说明参见数学实验p225本题可用hold on; % 在上面的同一张图上画线(同坐标系)plot(0,x,y,y, ':',x,x,0,y, ':');3)图形的标注使用文本标注函数text,调用格式如下：格式1text(x,y,文本标识内容:HorizontalAlignment ','字符串 1')x,y给定标注文本在图中添加的位置。'HorizontalAlignment '为水平控制属性，控制文本标识起点位于点(x,y)同一水平线上。'字符

5、串1'为水平控制属性值，取三个值之一：'left',点(x,y)位于文本标识的左边。'center '，点(x,y)位于文本标识的中心点。'right '，点(x,y)位于文本标识的右边。格式2text(x,y,文本标识内容:VerticalAlignment ','字符串 2')x,y给定标注文本在图中添加的位置。'VerticalAlignment '为垂直控制属性，控制文本标识起点位于点(x,y)同一垂直线上。'字符串1'为垂直控制属性值，取四个值之一：'middle

6、' ,' top ' ,' cap' ,' baseline ' ,' bottom '。(对应位置可在命令窗口应用确定)本题可用text(0,y,'(di/dt)m','VerticalAlignment','bottom');text(x,-0.001,num2str(x),'HorizontalAlignment','center');4)坐标轴标注调用函数 xlabel , ylabel 和 title本题可用title('SI

7、 模型 di/dti 曲线')；xlabel('i');ylabel('di/dt');2 .实验7-2传染病模型2 ( SI模型)一一画it曲线图(参考教材 p137-138 )传染病模型2 ( SI模型)：di/dt=ki(1-i),i(0)=i0;其中，i（t）是第t大病人在总人数中所占的比例。k是每个病人每天有效接触的平均人数（日接触率）。i0是初始时刻（t=0）病人的比例求出微分方程的解析解i（t）,画出如下所示的it曲线（i（0）=0.15, k=0.2,t=030 ）。试编写一个 m文件来实现。（在图形窗口菜单选择Edit/Copy Fig

8、ure ,复制图形）提示1）求解微分方程常微分方程符号解用函数 dsolve,调用格式如下：dsolve（ 'equ1 ' : equ2 '，变量名'）以代表微分方程及初始条件的符号方程为输入参数，多个方程或初始条件可在一个输入 变量内联立输入，且以逗号分隔。默认的独立变量为t,也可把t变为其他的符号变量。字符D代表对独立变量的微分，通常指 d/dt。本题可用x=dsolve（ ，Dx=k*x*（1-x） ，,，x（0）=x0，）2）画出 it 曲线（i（0）=0.15,入 =0.2, t=030 ）用 for 循环，函数 length, eval, plot,

9、 axis, title, xlabel, ylabel3 .实验7-3传染病模型3 （ SIS模型） 画di/dt i曲线图（参考教材 p138-139 ）已知传染病模型3 （ SIS模型）：di/dt=- ii-（1-1/）,i（0）=i 0其中，i（t）是第t大病人在总人数中所占的比例。人是每个病人每天有效接触的平均人数（日接触率）。i0是初始时刻（t=0）病人的比例。是整个传染期内每个病人有效接触的平均人数（接触数）。取入=0.1 , （t=1.5,画出如下所示的di/dt i曲线图。试编写一个 m文件来实现。（在图形窗口菜单选择 Edit/Copy Figure ,复制图形）提示用f

10、plot函数画出di/dti曲线图；在上图上用plot函数画一条过原点的水平用 title, xlabel, ylabel 标注。4 .实验7-4传染病模型3 （ SIS模型） 一一画it曲线图（参考教材 p138-139 ）已知传染病模型 3 （ SIS模型）：di/dt=- ii-（1-1/）,i（0）=i 0其中，i（t）是第t天病人在总人数中所占的比例。入是每个病人每天有效接触的平均人数（日接触率）。i0是初始时刻（t=0）病人的比例。b是整个传染期内每个病人有效接触的平均人数(接触数)。实验要求：求出微分方程的解析解i(t)。取入=0.2, (r=3, t=040，画出如下所示的图形

11、。试编写一个m文件来实现。图 1 型I曲线-=0 2, fT=3)05 1C 152025303540t(天】1 09 0.8 07 06 0.5 0.4 0.3 0.2 0.1D其中蓝色实线为i(0)=0.2 时的it曲线(第1条)；黑色虚点线为过点(0,1-1/ b )的水平线(第 2条)；红色虚线为i(0)=0.9 时的it曲线(第3条)。提示图例标注可用legend('i(0)=0.2','1-1/|','i(0)=0.9');5 .实验7-5传染病模型4 ( SIR模型)(参考教材 p140-141 )SIR 模型的方程：di/dt=s

12、i- i i(0)=i 0ds/dt=- si s(0)=s 0实验要求：1 .设 入=1 ,(1 =0.3 , i(0)=0.02 , s(0)=0.98 。输入p139 的程序，并修改程序中的t,x,使得输出的数据格式如下(提示：取 4位小数，使用四舍五入取整函数round ,矩阵剪裁和拼接)：ans =Columns 1 through 60 1 2 3 4 50.02 0.039 0.0732 0.1285 0.2033 0.27950.98 0.9525 0.9019 0.8169 0.6927 0.5438Columns 7 through 126 7 8 9 10 150.331

13、2 0.3444 0.3247 0.2863 0.2418 0.07870.3995 0.2839 0.2027 0.1493 0.1145 0.0543Columns 13 through 1820 25 30 35 40 450.0223 0.0061 0.0017 0.0005 0.0001 00.0434 0.0408 0.04010.0399 0.0399 0.03982 .运行结果与教材p140的内容比较。提示1）求解微分方程的数值解函数 ode45 ,格式如下：t,x=ode45（'fun',ts,x0）fun是由一个或多个待解方程写成的函数式m文件；ts=t0,

14、tf表示此微分方程的积分限是从 t0至U tf,也可以是一些离散的点，形式为 ts=to,ti,；X0为初值条件。2）等待用户反应命令pause :程序执行到该命令时暂停，直到用户按任意键后继续（处在命令窗口有效）。三、实验内容1 .实验7-1传染病模型2 （ SI模型） 画di/dt i曲线图在matlab 中建立 M 文件fun1.m代码如下：function y=fun（x）k=0.1;y=k*x*1-x;Fun2.m代码如下：function y=fun(x) k=0.1;y=-k*x*1-x;在命令行输入以下代码：fplot('funT,0 1.1 0 0.03);x=fmi

15、nbnd('fun2',0,1);y=0.1*x*(1-x);hold on;plot(0,x,y,y,'-',x,x,0,y");text(0,y,'(di/dt)m','VerticalAlignment','bottom');text(x,-0.001,num2str(x),'HorizontalAlignment','center');title('SI 模型 di/dti 曲线')；xlabel('i');ylabel('d

16、i/dt');hold off2 .实验7-2传染病模型2 ( SI模型)一一画it曲线图在matlab 中建立 M 文件fun22.m代码如下：k=0.2;x0=0.15;x=dsolve('Dx=k*x*(1-x)','x(0)=x0');tt=linspace(0,31,1001);fo门=1:1001t=tt(i)；xx(i)=eval(x);endplot(tt,xx)axis(0,31,0,1.1);title('图 1 SI 模型 it 曲线')；xlabel('t(天)')；ylabel('i (病

17、人所占比例)');在命令行输入以下代码：fun22;3 .实验7-3传染病模型3 ( SIS模型) 画di/dt i曲线图在matlab 中建立 M 文件fun3.m代码如下：function y=fun(x)a=0.1;b=1.5;y=-a*x*x-(1-1/b);在命令行输入以下代码：fplot('fun3',0 0.4 -0.0005 0.003);x=fminbnd('fun3',0,1);title('SIS 模型 di/dti 曲线')；xlabel('i');ylabel('di/dt');&

18、gt;> hold on>> plot(0,0.4,0,0)4 .实验7-4传染病模型3 ( SIS模型)一一画it曲线图在matlab 中建立 M 文件fun4.m代码如下：function y=fun(x)x=dsolve('Dx=-0.2*x*(x-(1-1/3)','x(0)=0.2');tt=linspace(0,41,1001);fo门=1:1001t=tt(i);xx(i)=eval(x);endplot(tt,xx);hold on;plot(0,40,1-1/3,1-1/3,'-k');x=dsolve(

19、9;Dx=-0.2*x*(x-(1-1/3)','x(0)=0.9');tt=linspace(0,41,1001);fo门=1:1001t=tt(i);xx(i)=eval(x);endplot(tt,xx,'-r');axis(0,40,0,1);title('图 1 SI 模型 it 曲线(入=0.2, (T =3)');xlabel('t(天)')；ylabel('i (病人所占比例)');legend('i(0)=0.2','1-1/6','i(0)=0.9

20、');在命令行输入以下代码：fun4;5 .实验7-5传染病模型4 ( SIR模型)在matlab 中建立 M 文件fun5.m代码如下：function y=fun(t,x)a=1;b=0.3;y=a*x(1)*x(2)-b*x(1),-a*x(1)*x(2)'在命令行输入以下代码：>> ts=0:50;>> x0=0.02,0.98;>> t,x=ode45('fun5',ts,x0);>> plot(t,x(:,1),t,x(:,2),grid,pause>> plot(x(:,2),x(:,1)

21、,grid,四、实验结果及其分析1 .实验7-1传染病模型2 （ SI模型）一一画di/dti曲线图分析:当i=1/2时di/dt达到最大值（di/dt） m，这时病人增加得在最快，可以认为是医院的门诊量最大的一天，预示着传染病高潮的到来，是医疗卫生部门关注的时刻。当t趋近于无穷时i趋近于1 ,即所有人终将被传染，全部变成病人，着显然不符合实际。原因是模型中没有考虑到病人可以治愈，人群中的健康者只能变成病人，病人不会再变成健康者。2 .实验7-2传染病模型2 （ SI模型）一一画it曲线图t（天）军三七左Y黑二分析：当i=1/2时di/dt达到最大值（di/dt） m，这时病人增加得在最快，可

22、以认为是医院的门诊量最大的一天，预示着传染病高潮的到来，是医疗卫生部门关注的时刻。当t趋近于无穷时i趋近于1 ,即所有人终将被传染，全部变成病人，着显然不符合实际。原因是模型中没有考虑到病人可以治愈，人群中的健康者只能变成病人，病人不会再变成健康者。3 .实验7-3传染病模型3 （ SIS模型） 画di/dt i曲线图分析：口是一个阈值，当u >1时,i的增减性取决于i0的大小，但其极限值i（无穷）=1-1/I。随。的增加而增加（试从I 口的含义给予解释）；当。<二1时，病人的比例i（t）越来越小。最终趋近于0,这是由于传染期内经有效接触从而使健康者变成病人数不超过原来的病人数的缘

23、4 .实验7-4传染病模型3 （ SIS模型） 一一画it曲线图10.908D.7S1 £1 模型i7曲线0 =0.2, a =3）06060.40.30.20.1510152025303540t（天）分析：口是一个阈值，当。1时,i的增减性取决于i0的大小，但其极限值i（无穷）=1-1/I。随。的增加而增加（试从1口的含义给予解释）；当。二1时，病人的比例i（t）越来越小。最 终趋近于0,这是由于传染期内经有效接触从而使健康者变成病人数不超过原来的病人数的缘 故。5.实验7-5传染病模型4 （ SIR模型）ans =00.02000.98001.00000.03900.95252.

24、00000.07320.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.050717.00000.04780.048018.00000.03710.046019.00000.0