1.实验7-1传染病模型2
传染病模型
个比例就不会改变。而当阈值提高时， 减小，于
是这个比例就会降低。
ห้องสมุดไป่ตู้
群体免疫和预防
根据对模型的分析，当 s0 1/时，传染病不会
蔓延，因而制止传染病蔓延的途径有两条
1．提高卫生和医疗水平（使阈值变大）；
2．通过预防接种使群体得到免疫（降低 s）0
r0
现在有一个公司想承接这项业务，他们允诺：提供一 种新型的可持久的檐槽，它包括一个横截面为半圆形(半径 为7．5厘米)的水槽和一个竖直的排水管(直径为10厘米)， 并且不管天气情况如何，这种檐槽都能排掉房顶的雨水．
但是房管部门还在犹豫，考虑公司的承诺能否实现，于 是想请你用数学的方法给一个详细的分析，论证它这个方案 的可行性
疾病的传染高峰期
此时 d 2I 意义： dt2
0
ln( n 1)
计算高峰期得： t0
I0 kn
1、当传染系数k或n增大时，t0随之减少，表示传 染高峰随着传染系数与总人数的增加而更快
的来临，这与实际情况比较符合。
2、令λ=kn，表示每个病人每天有效接触的平均
人数，称日接触率。t0与 λ成反比。 λ表示该 地区的卫生水平， λ越小卫生水平越高。故
1
1
（＊）
只要通过群体免疫使初始时刻的移出者比例
（即免疫者比例）r0 满足（＊）式，就可以制
止传染病的蔓延．
课后任务
请各位同学进行一些调查，根据模型算一 算在广州，非典型肺炎爆发的高潮大概是在何 时，与实际情况相吻合吗？根据模型请给出你 的建议。
思考题1
设某城市共有n+1人，其中一人出于某种目 的编造了一个谣言。该城市具有初中以上文化 程度的人占总人数的一半，这些人只有1/4相信 这一谣言，而其他人约有1/3会相信。又设凡相 信此谣言的人每人在单位时间内传播的平均人 数正比于当时尚未听说此谣言的人数，而不相 信此谣言的人不传播谣言。试建立一个反映谣言 传播情况的微分方程模型。
传染病最简单模型
传染病最简单模型：已感染人数 (病人) x(t)，每个病人每天有效接触(足以使人致病)人数为λ 有()()()x t t x t x t t λ+∆-=∆ 又设()00x x =，得微分方程dxx dtλ= 解得0()t x t x e λ=SI 模型：区分已感染者(病人)和未感染者(健康人)。
总人数N 不变，λ为日接触率，病人和健康人的比例分别为i(t),s(t)。
则有di si dt λ=，又有s(t)+i(t)=1。
所以有0(1),(0)dii i i i dtλ=-=。
求解出01()11(1)ti t e i λ-=+- ，传染速度最快时刻为101ln(1)mt i λ-=-SIS 模型：传染病无免疫性。
总人数N 不变，病人的日接触率为λ，病人和健康人的比例分别为i(t),s(t)，接触数σ(感染期内每个病人的有效接触人数)。
病人日治愈率为μ，所以有diN Nsi Ni dtλμ=- ， 0(0)i i =。
由s(t)+i(t)=1，/σλμ=，就推出1[(1)]di i i dt λσ=---。
SIR 模型：传染病有免疫性。
总人数N 不变，病人、健康人和移出者的比例分别为i(t)，s(t)，r(t) ，病人的日接触率为λ，病人日治愈率为μ，接触数/σλμ=。
且有s(t)+i(t)+r(t)=1。
则有r(0)=r0很小，故000i s +≈。
推出00d ,(0)d d ,(0)d i si i i i ts si s s t λμλ⎧=-=⎪⎪⎨⎪=-=⎪⎩ 经济增长模型；1 ）道格拉斯(Douglas)生产函数 Q(t)，K(t)，L(t)，0f 分别表示某地区在t 时刻的产值、资金、劳动力和技术。
静态模型令z=Q/L ，y=K/L ，则z 是每个劳动力产值，y 是每个劳动力投资。
由于z 随y 增加而增长，但增速递减。
)(/0y g f L Q z ==，10,)(<<=ααy y g ，α)/(0L K L f Q =αα-=10),(L K f L K Q 此为Douglas 生产函数。
数学建模例题题
数学建模试题一、传染病模型医学科学的发展已经能够有效地预防和控制许多传染病，但是仍然有一些传染病暴发或流行，危害人们的健康和生命。
社会、经济、文化、风俗习惯等因素都会影响传染病的传播，而最直接的因素是：传染者的数量及其在人群中的分布、被传染者的数量、传播形式、传播能力、免疫能力等。
一般把传染病流行范围内的人群分成三类：S类，易感者(Susceptible)，指未得病者，但缺乏免疫能力，与感染者接触后容易受到感染；I类，感病者(Infective)，指染上传染病的人，它可以传播给S类成员；R类，移出者(Removal)，指被隔离或因病愈而具有免疫力的人。
要求：请建立传染病模型，并分析被传染的人数与哪些因素有关？如何预报传染病高潮的到来?为什么同一地区一种传染病每次流行时，被传染的人数大致不变？二、线性规划模型—销售计划问题某商店拟制定某种商品7—12月的进货、售货计划，已知商店仓库最大容量为1500件，6月底已存货300件，年底的库存以不少于300件为宜，以后每月初进货一次，假设各月份该商品买进、售出单价如下表。
要求：若每件每月的库存费用为0.5元，问各月进货、售货各为多少件，才能使净收益最多？建立数学模型，并用软件求解。
【注】线性规划在MATLAB的库函数为：linprog。
语法为：x = linprog(f,A,b)x = linprog(f,A,b,Aeq,beq,lb,ub,x0,options)[x,fval,exitflag,output,lambda] = linprog(...)例如：线性规划目标函数的系数：f = [-5; -4; -6]约束方程的系数及右端项：A = [1 -1 13 2 43 2 0];b = [20; 42; 30];lb = zeros(3,1);调用线性规划程序linprog求解，得：[x,fval,exitflag,output,lambda] = linprog(f,A,b,[],[],lb);x= 0.000015.00003.0000三、一阶常微分方程模型—人口模型与预测 下表列出了中国1982-1998年的人口统计数据，取1982年为起始年（0=t ），1016540=N 万人，200000=m N 万人。
传染病模型——精选推荐
3.12传染病模型摘要：本文是一个对传染病的研究问题。
通过把一般把传染病流行范围内的人群分成三类：S类，易感者(Susceptible)，指未得病者，但缺乏免疫能力，与感染者接触后容易受到感染；I类，感病者(Infective)，指染上传染病的人，它可以传播给S类成员；R类，移出者(Removal)，指被隔离或因病愈而具有免疫力的人。
建立数学模型用极限和微积分等数学方法对传染病传播规律进行研究。
关键词：传染病极限和微积分正文1 传染病〔Infectious Diseases〕是由各种病原体引起的能在人与人、动物与动物或人与动物之间相互传播的一类疾病。
病原体中大部分是微生物，小部分为寄生虫，寄生虫引起者又称寄生虫病。
有些传染病，防疫部门必须及时掌握其发病情况，及时采取对策，因此发现后应按规定时间及时向当地防疫部门报告，称为法定传染病。
中国目前的法定传染病有甲、乙、丙3类，共37种医学科学的发展已经能够有效地预防和控制许多传染病，天花在世界范围内被消灭，鼠疫、霍乱等传染病得到控制。
但是仍然有一些传染病暴发或流行，危害人们的健康和生命。
在发展中国家，传染病的流行仍十分严重；即使在发达国家，一些常见的传染病也未绝迹，而新的传染病还会出现，如爱滋病（AIDS）等。
有些传染病传染很快，导致很高的致残率，危害极大，因而对传染病在人群中传染过程的定量研究具有重要的现实意义。
传染病流行过程的研究与其他学科有所不同，不能通过在人群中实验的方式获得科学数据。
事实上，在人群中作传染病实验是极不人道的。
所以有关传染病的数据、资料只能从已有的传染病流行的报告中获取。
这些数据往往不够全面，难以根据这些数据来准确地确定某些参数，只能大概估计其范围。
基于上述原因，利用数学建模与计算机仿真便成为研究传染病流行过程的有效途径之一。
2问题提出上世纪初，瘟疫还经常在世界的某些地区流行，被传染的人数与哪些因素有关？如何预报传染病高潮的到来？为什么同一地区一种传染病每次流行时，被传染的人数大致不变？3 模型分析社会、经济、文化、风俗习惯等因素都会影响传染病的传播，而最直接的因素是：传染者的数量及其在人群中的分布、被传染者的数量、传播形式、传播能力、免疫能力等，在建立模型时不可能考虑所有因素，只能抓住关键的因素，采用合理的假设，进行简化。
传染病模型
模型2（SI模型） 假设条件为
1.在疾病传播期内所考察地区的总人 数N不变，即不考虑生死，也不考虑迁移。 人群分为易感染者（Susceptible）和已感 染者（Infective）两类（取两个词的第一 个字母，称之为SI模型），以下简称健康 者和病人。时刻t这两类人在总人数中所 占比例分别记作s(t)和i(t)。 2.每个病人每天有效接触的平均人数 是常数，称为日接触率。当病人与健康 者接触时，使健康者受感染变为病人。
再设t 0时有x0有个病人，即得微分方程
dx x , dt x ( 0) x
0
(1)
方程（1）的解为
x(t ) x e
0
t
(2)
结果表明，随着t的增加，病人人数x(t)无 限增长，这显然是不符合实际的。 建模失败的原因在于：在病人有效接 触的人群中，有健康人也有病人，而其中 只有健康人才可以被传染为病人，所以在 改进的模型中必须区别这两种人。
不同类型传染病的传播过程有其各自 不同的特点，弄清这些特点需要相当多的 病理知识，这里不可能从医学的角度一一 分析各种传染病的传播，而只是按照一般 的传播模型机理建立几种模型。 模型1 在这个最简单的模型中，设时 刻t的病人人数x(t)是连续、可微函数， 并且每天每个病人有效 接触(足使人致病) 的人数为常数 考察 t到 t t 病人人数的 增加，就有 x(t t ) x(t ) x(t )t
传染病模型
随着卫生设施的改善、医疗水平的提 高以及人类文明的不断发展，诸如霍乱、 天花等曾经肆虐全球的传染性疾病已经得 到有效的控制。但是一些新的、不断变异 着的传染病毒却悄悄向人类袭来。20世纪 80年代十分险恶的爱滋病毒开始肆虐全球， 至今带来极大的危害。长期以来，建立制 止传染病蔓延的手段等，一直是各国有关 专家和官员关注的课题。
传染病模型在不同场景下的应用(2)
传染病模型在不同场景下的应用在过去的几十年里，传染病模型在许多不同的场景中都得到了广泛的应用。
从疫情爆发的防控，到公共卫生政策的制定，再到疫苗的研发，传染病模型都发挥了重要的作用。
让我们来看一下传染病模型在疫情爆发时的应用。
当一种新的传染病爆发时，政府和卫生部门需要迅速采取措施来控制疫情的蔓延。
传染病模型可以帮助他们预测疫情的发展趋势，确定疫情的热点地区，以及评估不同防控措施的效果。
通过这些信息，政府可以制定出更加科学合理的防控策略，从而有效地控制疫情的蔓延。
除了在疫情爆发时的应用，传染病模型还在公共卫生政策的制定中扮演着重要的角色。
公共卫生政策的制定需要考虑到许多不同的因素，如人群的流动性，医疗资源的分布，以及人们的行为习惯等。
传染病模型可以帮助政策制定者更好地理解这些因素对疫情的影响，从而制定出更加有效的公共卫生政策。
另外，传染病模型还在疫苗的研发中发挥着重要的作用。
疫苗的研发需要进行大量的实验和临床试验，而这些实验和临床试验需要依据一定的理论模型来进行。
传染病模型可以帮助研究人员更好地理解疫苗的效果和安全性，从而为疫苗的研发提供有力的理论支持。
传染病模型在许多不同的场景中都得到了广泛的应用。
从疫情爆发的防控，到公共卫生政策的制定，再到疫苗的研发，传染病模型都发挥了重要的作用。
我相信，在未来的发展中，传染病模型将继续发挥其重要的作用，为人类的健康事业做出更大的贡献。
在处理疫情爆发时，传染病模型的应用是至关重要的。
以2003年非典型肺炎（SARS）为例，当时的疫情迅速蔓延至全球多个国家和地区。
通过建立传染病模型，研究人员和公共卫生官员能够预测疫情的发展趋势，评估防控措施的效果，从而指导政府和卫生部门采取有效的应对措施。
例如，模型可以帮助确定隔离措施和疫苗接种策略的重点地区，确保资源的有效分配。
在疫情期间，我和我的团队也参与了传染病模型的构建和应用，通过不断优化模型参数，我们能够更准确地预测疫情走势，为决策者提供关键信息。
数学建模传染病模型例题
数学建模传染病模型例题一、传染病模型简介传染病模型是数学建模的一个重要分支，主要用于描述传染病在人群中的传播规律。
通过构建合适的数学模型，可以研究传染病的传播动力学、预测疫情发展趋势以及评估防控措施的效果。
本文将重点介绍几种常见的传染病模型及其应用。
二、传染病模型的类型及应用1.SIR模型SIR模型是一种基于微分方程的传染病模型，其中S、I、R分别代表易感者（Susceptible）、感染者（Infected）和康复者（Recovered）。
该模型通过描述易感者感染、感染者康复以及康复者不再易感的动态过程，揭示了传染病在人群中的传播规律。
SIR模型在分析疫情爆发、研究防控措施等方面具有广泛应用。
2.SEIR模型SEIR模型是在SIR模型基础上发展的一种传染病模型，其中E代表潜伏者（Exposed）。
与SIR模型相比，SEIR模型增加了潜伏期这一概念，使得模型更加符合实际情况。
该模型可以用于研究传染病的传播速度、预测疫情发展趋势以及评估疫苗的效果。
3.SI模型SI模型是一种简化的传染病模型，仅包含易感者和感染者两个群体。
该模型适用于分析短期传染病，如流感等。
通过研究易感者与感染者的动态关系，可以预测疫情爆发的时间和规模。
三、传染病模型的参数估计与预测传染病模型的参数估计是数学建模的关键环节，通常采用最大似然估计、贝叶斯估计等方法。
此外，基于传染病模型的预测技术在疫情防控中也具有重要意义。
通过构建时间序列模型，如ARIMA、SVM等，可以预测未来一段时间内疫情的发展趋势。
四、数学建模在传染病防控中的实际应用数学建模在传染病防控中具有广泛应用，如疫情监测、防控措施评估、疫苗研究等。
通过对传染病模型的深入研究，可以为政府部门提供科学依据，协助制定针对性的防控策略。
五、案例分析本文将结合具体案例，如我国2003年非典疫情、2020年新冠肺炎疫情等，详细阐述传染病模型在实际应用中的重要作用。
通过分析案例，可以加深对传染病模型的理解，为今后疫情防控提供借鉴。
传染病模型精选推荐(一)2024
传染病模型精选推荐（一）引言：传染病模型是研究传染病传播方式和防控策略的重要工具。
本文将介绍5个精选的传染病模型，并探讨它们的特点和应用领域。
大点一：SIR模型1. SIR模型是传染病模型中最基本的一种，包括易感者(Susceptible)、感染者(Infected)和康复人群(Recovered)。
2. SIR模型适用于研究人群中的疾病传播情况，可以预测传染病的爆发和蔓延趋势。
3. SIR模型假设人群中没有出生死亡和迁移，并且感染后具有免疫力。
4. SIR模型可以通过改变参数来研究不同防控措施的效果，如隔离、疫苗接种等。
大点二：SEIR模型1. SEIR模型在SIR模型的基础上增加了潜伏期(Exposed)的状态，即潜伏期内已经感染但还未展现症状的人群。
2. SEIR模型适用于研究传染病的潜伏期和潜伏期内的传播方式。
3. SEIR模型可以更准确地描述疾病的传播过程，并提供更精确的防控策略。
4. SEIR模型可以通过添加接触率和潜伏期的参数来模拟不同传染性和潜伏期的疾病。
大点三：SEIRD模型1. SEIRD模型在SEIR模型的基础上增加了死亡者(Death)的状态，用于研究传染病的死亡率和致死风险。
2. SEIRD模型适用于研究死亡率高的传染病，如高致病性禽流感等。
3. SEIRD模型可以通过改变死亡率和康复率的参数来预测传染病的死亡数量和康复情况。
4. SEIRD模型有助于评估不同防控策略对死亡率的影响，如加强医疗资源、提高疫苗接种率等。
大点四：Agent-based模型1. Agent-based模型是一种基于个体行为和交互的传染病模型。
2. Agent-based模型可以模拟个体之间的接触和传播过程，更加现实和细致。
3. Agent-based模型适用于研究人口密集区域的传染病传播，如城市、机场等。
4. Agent-based模型能够考虑到不同个体的行为差异和健康状态，有助于制定个体化的防控策略。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
河北大学《数学模型》实验实验报告一、实验目的二、实验要求1.实验7-1 传染病模型2（ SI模型）——画di/dt~ i曲线图（参考教材 p137-138）传染病模型 2（ SI 模型）：;di/dt=ki(1-i),i(0)=i其中， 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 轴也被限制。
本题可用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,’颜色线型数据点图标’, 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)同一水平线上。
’字符串 1’为水平控制属性值，取三个值之一：‘left’，点(x,y)位于文本标识的左边。
‘center’，点(x,y)位于文本标识的中心点。
‘right’，点(x,y)位于文本标识的右边。
格式 2text(x,y, 文本标识内容,’VerticalAlignment’,’字符串 2’)x,y 给定标注文本在图中添加的位置。
’VerticalAlignment’为垂直控制属性，控制文本标识起点位于点(x,y)同一垂直线上。
’字符串 1’为垂直控制属性值，取四个值之一：‘middle’，’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模型di/dt~i曲线');xlabel('i');ylabel('di/dt');2.实验7-2 传染病模型2（ SI模型）——画i~t曲线图（参考教材 p137-138）传染病模型 2（ SI 模型）：di/dt=ki(1-i),i(0)=i;其中，i(t)是第t 天病人在总人数中所占的比例。
k 是每个病人每天有效接触的平均人数（日接触率）。
i0是初始时刻（t=0）病人的比例求出微分方程的解析解 i(t)，画出如下所示的 i~t 曲线（ i(0)=0.15, k=0.2,t=0~30）。
试编写一个 m 文件来实现。
（在图形窗口菜单选择 Edit/Copy Figure，复制图形）[提示]1）求解微分方程常微分方程符号解用函数 dsolve,调用格式如下：dsolve(‘equ1’,’equ2’,…,’变量名’)以代表微分方程及初始条件的符号方程为输入参数，多个方程或初始条件可在一个输入变量内联立输入，且以逗号分隔。
默认的独立变量为 t，也可把 t 变为其他的符号变量。
字符 D 代表对独立变量的微分，通常指 d/dt。
本题可用x=dsolve(‘Dx=k*x*(1-x)’,’x(0)=x0’)2) 画出 i~t 曲线（ i(0)=0.15, λ=0.2, t=0~30）用 for 循环，函数 length, eval, plot, axis, title, xlabel, ylabel3.实验7-3 传染病模型3（ SIS模型）——画di/dt~ i曲线图（参考教材 p138-139）已知传染病模型 3（ SIS 模型）：di/dt=- i[i-(1-1/ )],i(0)=i其中，i(t)是第t 天病人在总人数中所占的比例。
λ是每个病人每天有效接触的平均人数（日接触率）。
i0是初始时刻（t=0）病人的比例。
σ是整个传染期内每个病人有效接触的平均人数（接触数）。
取λ=0.1，σ =1.5，画出如下所示的di/dt~ i曲线图。
试编写一个 m 文件来实现。
（在图形窗口菜单选择 Edit/Copy Figure，复制图形）[提示]用fplot函数画出di/dt~ i曲线图；在上图上用plot函数画一条过原点的水平用title, xlabel, ylabel标注。
4.实验7-4 传染病模型3（ SIS模型）——画i~t曲线图（参考教材 p138-139）已知传染病模型 3（ SIS 模型）：di/dt=- i[i-(1-1/ )],i(0)=i其中，i(t)是第t 天病人在总人数中所占的比例。
λ是每个病人每天有效接触的平均人数（日接触率）。
i0是初始时刻（t=0）病人的比例。
σ是整个传染期内每个病人有效接触的平均人数（接触数）。
实验要求 :求出微分方程的解析解 i(t)。
取λ=0.2, σ =3, t=0~40，画出如下所示的图形。
试编写一个 m 文件来实现。
其中蓝色实线为 i(0)=0.2 时的 i~t 曲线（第 1 条）；黑色虚点线为过点（ 0, 1-1/σ）的水平线（第 2 条）；红色虚线为 i(0)=0.9 时的 i~t 曲线（第 3 条）。
[提示]图例标注可用legend('i(0)=0.2','1-1/¦σ','i(0)=0.9');5.实验7-5 传染病模型4（ SIR模型）（参考教材 p140-141）SIR 模型的方程 :di/dt= si- i i(0)=i0ds/dt=- si s(0)=s实验要求：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.3312 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 45三、实验内容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('fun1',[0 1.1 0 0.03]);x=fminbnd('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/dt~i曲线');xlabel('i');ylabel('di/dt');hold off2.实验7-2 传染病模型2（ SI模型）——画i~t曲线图在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);for i=1:1001t=tt(i);xx(i)=eval(x);endplot(tt,xx)axis([0,31,0,1.1]);title('图1 SI模型i~t曲线');xlabel('t(天)');ylabel('i（病人所占比例）');在命令行输入以下代码：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/dt~i曲线');xlabel('i');ylabel('di/dt');>> hold on>> plot([0,0.4],[0,0])4.实验7-4 传染病模型3（ SIS模型）——画i~t曲线图在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);for i=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('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);endplot(tt,xx,'-r');axis([0,40,0,1]);title('图1 SI模型i~t曲线(λ =0.2, σ =3)');xlabel('t(天)');ylabel('i（病人所占比例）');legend('i(0)=0.2','1-1/σ ','i(0)=0.9');在命令行输入以下代码：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)),grid,四、实验结果及其分析1.实验7-1 传染病模型2（ SI模型）——画di/dt~ i曲线图分析：,这时病人增加得在最快，可以认为是医院的门诊当i=1/2时di/dt达到最大值(di/dt)m量最大的一天，预示着传染病高潮的到来，是医疗卫生部门关注的时刻。