首页 > 数学建模实验答案微分方程模型

数学建模实验答案微分方程模型

发布时间：2023-03-27 13:04:29 来源：文档文库

小中大

字号：

手机查看

实验07微分方程模型（2学时）
（第5章微分方程模型）
1.（验证）传染病模型2（SI模型）p136~138传染病模型2（SI模型）：
diki(1i,i(0i0dt其中，
i(t是第t天病人在总人数中所占的比例。
k是每个病人每天有效接触的平均人数（日接触率）。i0是初始时刻（t=0）病人的比例。
di1.1画~i曲线图p136~138dt取k=0.1，画出上标注。
参考程序：
%传染病模型2（SI模型）的di/dt~i曲线图%文件名：p137fig2.m%λ=0.1
clear;clc;fplot('0.1*x*(1-x',[01.100.03];x=fminbnd('-0.1*x*(1-x',0,1y=0.1*x*(1-xholdon
plot([0,x],[y,y],':',[x,x],[0,y],':';text(0,y,'(di/dtm','VerticalAlignment','bottom';text(x,-0.001,num2str(x,'HorizontalAlignment','center';didi~i的曲线图，求i为何值时达到最大值，并在曲线图dtdt
1
title('SI模型的di/dt~i曲线';xlabel('i';ylabel('di/dt';holdoff;提示：fplot,fminbnd,plot,text,title,xlabel1）画曲线图
用fplot函数，调用格式如下：
fplot(fun,limsfun必须为一个M文件的函数名或对变量x的可执行字符串。若lims取[xminxmax]，则x轴被限制在此区间上。若lims取[xminxmaxyminymax]，则y轴也被限制。本题可用
fplot('0.1*x*(1-x',[01.100.03];
2）求最大值
用求解边界约束条件下的非线性最小化函数fminbnd，调用格式如下：x=fminbnd('fun',x1,x2fun必须为一个M文件的函数名或对变量x的可执行字符串。返回自变量x在区间x1上函数取最小值时的x值。本题可用
x=fminbnd('-0.1*x*(1-x',0,1y=0.1*x*(1-x
3）指示最大值坐标
用线性绘图函数plot，调用格式如下：
plot(x1,y1,'颜色线型数据点图标',x2,y2,'颜色线型数据点图标',…本题可用

2

本文来源：https://www.2haoxitong.net/k/doc/5cc3d8a1aff8941ea76e58fafab069dc5122474c.html