- 1、本文档共7页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 8、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
精品文档
数学模型实验—实验报告10
学院: 专 业: 姓 名:
学号:___ ____ 实验时间:__ ____ 实验地点:
一、实验项目:传染病模型求解
二、实验目的和要求
a.求解微分方程的解析解
b.求解微分方程的数值解
三、实验内容
问题的描述
各种传染病给人类带来的巨大的灾难,长期以来,建立传染病的数学模型来描述传染病的的传播过
程,分析受感染人数的变化规律,探索制止传染病蔓延的手段等,一直是各国有关专家和官员关注的课
题。
不同类型传染病有各自不同的特点,在此以一般的传播机理建立几种3 模型。分别对3 种建立成功
的模型进行模型分析,便可以了解到该传染病在人类间传播的大概情况。
模型一(SI模型):
(1)模型假设
1.在疾病传播期内所考察地区的总人数N 不变,人群分为健康人和病人,时刻t 这两类人在总人数中所
占比例为s (t )和i (t )。
2.每个病人每天有效接触的平均人数是常数 a,a 成为日接触率,当病人与健康者有效接触时,可使其
患病。
(2)建立模型
根据假设,每个病人每天可使as (t )个健康人变成病人,t 时刻病人数为Ni (t ),所以每天共有aNs (t )
i (t )个健康者被感染,即病人的增加率为:Ndi/dt=aNsi
又因为s (t )+i (t )=1
再记时刻t=0 时病人的比例为i0
则建立好的模型为:
di
ai(1i)
dt
i(0)=i0
(3)模型求解 (代码、计算结果或输出结果)
syms a i t i0 % a :日接触率,i :病人比例, s:健康人比例,i0 :病人比例在t=0 时的值
i=dsolve('Di=a*i*(1-i)','i(0)=i0','t');
y=subs(i,{a,i0},{0.3,0.02});
ezplot(y,[0,100])
figure
精品文档
精品文档
i=str2double(i);
i=0:0.01:1;
y=0.3*i.*(1-i);
plot(i,y)
SI 模型的i~t 曲线 SI 模型的di/dt~i 曲线
(4)结果分析
由上图可知,在i=0:1 内,di/dt 总是增大的,且在i=0.5 时,取到最大值,即在t->in 时,所有人
都将患病。
上述模型显然不符合实际,为修正上述结果,我们重新考虑模型假设,建立SIS 模型
模型二(SIS模型)
(1) 模型假设
假设条件1.2 与SI 模型相同;
3.每天被治愈的病人数占病人总数的比例为常数u ,成为日治愈率,病人治愈后成为仍可被感染的
健康者。显然1/u 是平均传染期。
(2)模型建立
病人的增加率:Ndi/dt=aNsi-uNi 且 i (t )+s(t)=1 ;
则有: di/dt=ai(1-i)-ui
在此定义k=a/b,可知k 是整个传染传染期内每个病人有效接触的平均人数,成为接触数。
则建立好的模型为:
di
ai[i(11/ k)]
dt
i(0)=i0;
(2) 模型求解 (代码、计算结果或输出结果)
>> syms a i u t i0 % a:日接触率,i :病人比例,u :日治愈率,i0 :病人比例在t=0 时的值
>> dsolve('Di=a*i*(1-i)-u*i','i(0)=i0','t') % 求用u 表示的i—t 解析式
>> syms k % k :接触数
>> k=a/u;
>> i=dsolve('Di=-a*i*i+a*i*(1-1/k)','i(0)=i0','t') % 求用k 表示的i—t 解析式
% 给k 、a、i0 指定特殊值,作出相关图像
精品文档
精品文档
>> y=subs(i,{k
您可能关注的文档
- 北师大版五年级数学下册第三单元分数乘法单元测试卷(含答案).pdf
- 北师大版六年级上册数学百分数的应用单元测试(含答案).pdf
- 北师大版八年级数学上册:第三章《位置与坐标》教案.pdf
- 北师大版六年级下册数学《圆柱与圆锥》单元测试卷(含答案).pdf
- 北师大版四年级语文下册《跳水》拔高练习.pdf
- 北师大版小学三年级语文测试上册第八单元练习题.pdf
- 北师大版小学数学一年级下册第一单元加与减(一)单元检测卷 (含答案).pdf
- 北师大版数学二年级下册第七、八单元测试卷(含答案).pdf
- 北师大版数学四年级下册期末测试卷(含答案).pdf
- 华东师大版2019-2020学年七年级数学下册第七章 一次方程组单元测试卷(含答案).pdf
文档评论(0)