最新SI.SIR.SIS 模型资料.pdfVIP

  1. 1、本文档共7页,可阅读全部内容。
  2. 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
  3. 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载
  4. 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
  5. 5、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
  6. 6、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们
  7. 7、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
  8. 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(1i) 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(11/ 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

文档评论(0)

qiuzhimin2019 + 关注
实名认证
文档贡献者

该用户很懒,什么也没介绍

1亿VIP精品文档

相关文档