《刘学江传染病接种问题模型》由会员分享，可在线阅读，更多相关《刘学江传染病接种问题模型（23页珍藏版）》请在金锄头文库上搜索。
1、2009-2010第二学期数学模型 期末考试承 诺 书我完全明白，在期末考试不能以任何方式（包括电话、电子邮件、网上咨询等）与队外的任何人研究、讨论与题有关的问题。我知道，抄袭别人的成果是违反考试规则的, 如果引用别人的成果或其他公开的资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参考文献中明确列出。我郑重承诺，严格遵守考试规则，以保证考试的公正、公平性。如有违反考试规则的行为，我将受到严肃处理。专业名称 ： 1. 数学与应用数学 考试题目 ： 2. 流感疫苗接种问题 姓名 (打印并签名) ：3. 刘学江 班级： 4. 2008级1班 学号： 5. 20080601
2、22 成绩： 6. 日期： 2010 年 06月 05 日 目录零 封面. . . .1一流感疫苗接种的数学模型. .2 1. 摘要 . .2 2. 问题的重述. .2 3. 问题的分析. .2 4. 详细设计. .4 4.1 初等数学模型. . .44.2 房室模型模型的建立. . . . . .9 5. 模型检验. .17二 总结. .22三 参考文献. . . . . . . . . . . . . . .22 一流感疫苗接种的数学模型1摘要本文通过初等数学模型，病毒指数增长模型，Logistic模型，SIS模型，SIR模型，随机房室模型分析菌种1、2的感染人数，计算相应的疫苗接种人数和
3、患病人数，通过模型间的对比找出用尽可能少的疫苗预防尽可能多的人数。对病毒指数模型只作定性分析，重点讨论:初等数学模型，随机房室模型，SIS模型，SIR模型。关键词：流感病毒菌种指数生长模型 随机房室模型 SIS模型 SIR模型 Logistic模型2问题的重述流感病毒有两种菌种，现已研制成两种疫苗。疫苗1对菌种1有85的预防效果，对菌种2有70的预防效果；疫苗2对菌种1有60的预防效果，对菌种2有90的预防效果。两种疫苗不能在一个人身上同时使用。(1)为使尽可能多的居民具有免疫力，需要进一步了解那些信息？(2)为使尽可能多的居民具有免疫力，应采取何种接种疫苗的策略？(3)在采取你所推荐的策略的
4、情况下，估计有多少居民具有免疫力(平均的估计和最坏情况的估计）。3问题的分析问题1：针对问题1主要从以下3方面进行：（1） 菌种角度：1、菌种的传播方式2、菌种的传播能力3、菌种的首次感染率4、菌种的生长方式5、菌种的致死率6、菌种的交叉感染率7、菌种的再次感染能力（2） 疫苗角度：1、疫苗的预防能力菌种12疫苗10.850.720.60.9 2.、疫苗的数量3、疫苗的接种限制4、疫苗的价格5、疫苗的有效期（3） 地区角度：1、该地区的总人口数2、初始时刻感染人数3、地区的预防效率4、每人感染后的接触人数5、每种菌种的治愈可能问题2:“民具有免疫力，应采取何种接种疫苗的策略”采用模型对比法找出
5、最佳接种策略，免疫人数P=P1接种疫苗人数-P2接种疫苗后患病人数，为使P尽可能的大可以使P1尽可能的大P2尽可能的小，但在实际接种过程中考虑人力和财力因素不可能使接种疫苗的人数无限制的多，反而是在高效的前提下尽可能得节约疫苗对人力财力和居民健康也有好处，这里采取的接种方式是对感染者周围的一定人群接种，对其余人群均不接种疫苗，所以对疫苗的接种数和有效性是问题2的关键。问题3： “在采取你所推荐的策略的情况下，估计有多少居民具有免疫力(平均的估计和最坏情况的估计” 。通过相关参数的不同设置拟合出平均和最坏情况4详细设计符号说明：P表示感染或患病人数 N表示该地区总人口 n表示菌种感染者接触的人数
6、 t表示时间 K表示某些比例系数 c表示菌种致死率 m表示疫苗不能接种的人数比例4.1初等数学模型 模型1 设t时刻的病人数X（t）是连续可微函数，设每人接触的人数是n, 从 t到t+t时刻病人人数的增加，X(t+t)- X（t）=n *X（t）*t 设初始时刻有X。,得微分方程dx/dt=n*x,x(0)= X。解得X（t）=n*exp(n*t),病人数呈指数增长不符合实际情况，根据 SIS模型 SIR模型 Logistic模型 理论感染人数不可能无限制的增长，如下图所示。针对疫苗接种情形通常在病情初发期进行的，所以可以将病人数看做指数型增长对菌种1、2建立不同的指数接触人数增长模型：P1t
7、、P2t表示t时刻人数，在从 t到t+t时刻可能感染菌种1的病人人数的增加为P1(t+t)- P1t=n1* P1t*t , P1t(0)=P10得 解得P1t=P01*exp(n1*t) 除去原始感染菌种1人数剩下的就是被接触人数，也就是需要接种疫苗的目标人数将P1t修正为P1t=P01*exp(n1*t)-P01在从 t到t+t时刻可能感染菌种2的病人人数的增加为P2(t+t)- P2t=n2* P2t*t , P2t(0)=P20得 解得P2t=P02*exp(n2*t) 除去原始感染菌种2人数剩下的就是被接触人数，也就是需要接种疫苗的目标人数将P2t修正为P2t=P02*exp(n2*
8、t)-P02由于菌种1、2均有致死率c1、c2，还要在感染的基础上除去死亡人数则菌种1在t时刻的感染人数：P1t=P01*exp(n1*t) （1-c1）。.菌种2在t时刻的感染人数为：P2t=P02*exp(n2*t)(1-c2)。.若采用疫苗1接种具有免疫力的人数P1=P1t*0.85+P2t*0.7。 若采用疫苗2接种具有免疫力的人数P2=P1t*0.6+p2t*0.9。. 比较P1 、P2的大小就可以得到采用哪种疫苗作为接种策略。在实际接种中常常采用不同疫苗针对不同菌种的接种方法，因此模型的接种策略需要修正。模型2对模型进行修正：菌种12疫苗10.850.720.60.9疫苗1、2针对
9、菌种的预防能力由上表可知对菌种1采用疫苗1接种最高效啊，对菌种2采用疫苗2接种最高效。因此对模型一的进行修正得到新的免疫人数：P=P1t*0.85+P2t*0.9。在实际疫苗接种工作每种疫苗都有不能接种的人群，在这里假设疫苗不能接种的人数比例为m1,疫苗2不能接种的人数比例为m2因此对还需啊哟修正：P= P1t(1-m10*0.85+P2t（1-m2）*0.9。带入各项数据得：P= P01*exp(n1*t)（1-c1）(1-m1)*0.85+ P02*exp(n2*t)(1-c2)（1-m2）*0.9。采用C语言编程：#include(math.h)#include(stdio.h)Void mian()double n1,n2,c1,c2,m1,m2;/*定义疫苗感染者的接触人数n，菌种的致死率c，疫苗不能接种的人群比例m*/ Int t； /*定义时间*/Double P1t,P2t;/*t时刻防疫菌种1、2的人数*/Double p;/*最后的防疫总人数*/Int P01，P02;/*定义初始时刻的菌种感染人数*/Double P1t,P2t;；/*定义菌种1。2在
