1、图形1中的横坐标数字表示时间的天数，如15即4月20日之后的第15天，40即4月20日之后的第40天。在45天之后的时间里，模型对k的值进行了调整，k=0.0273，我们再将k=0.0273代入模型 N（t）= （1+K），在45天之后的时间里，我们取了30天的时间，分别算出每天的病人累计数，如图3-3所示：图3-33.3对指数模型的验证和评价在图形3-3中的横坐标的数值表示图形1中所表示的天数之后的天数，如1即表示4月15日之后的45天之后的有第六天，也就是4月15日之后的第51天，即表示4月15日之后的第67天。首先在图形3-3结合图形3-1可以看出，图形3-1中的第45天与图形2中的第一
2、天（相隔一天）的人数统计是相差比较大的，存在这种情况的原因是在我们在计算第61天，数据值发生了改变，从0.16204到0.0273是一个很大的变化，而在实际的生活中的情况是k值每天都在进行数值在减小的改变，但改变的没有这么大，也正是因为k有了跳跃，N（t）的值才会发生这么大的变化，这是可以理解的。我们对图形2的整个曲线来与附件1中的图形1进行比较，可以发现，在整个阶段的数值曲线图形都是很接近的。我们在对全国在前期和后期k分别取k=0.16204和k=0.0273的值来代入所给的模型来计算并画出的图形,与实际的数据和图形进行了比较，是有着很好的吻合，同样我们也可以对k取值一个定值来对全国进行计算
3、和画图，同样也是合理的。因此我们就认为题目中给我们的那个模型N（t）= 是合理的。通过这个模型我们可以根据某一地区的疫情从爆发到高潮或某一阶段的时间的长短来拟合得到一个与该地区这种疫情的感染率，就可以用该模型来计算或预测该地区现在及以后的病人的累计数， 这也就是该模型的实用性所在。4建立新模型4.1模型假设模型假设与指数模型假设一致不在赘述。4.2模型分析与建立4.2.1模型分析初期由于疫情初期政府控制力度不够，大众的对SARS的防范意识不强，造成病情迅速蔓延。而当政府采取有力措施，人们的防患意识增强，疫情则趋于缓和，病患者人数迅速下降。所以SARS传播大体上可分为两个阶段：1）控制前期：即认
4、为病毒传播方式是自然传播。2）控制后期：政府强力介入之后的病毒传播模型。4.2.2 模型建立根据对指数模型的分析和4.2的分析疫情走势的微分方程如下；N（t） = K N（t） N（t L） . （1）4.3模型的求解如果假定有一个初始爆发时间，最初有N0 个病人突然出现，在L 天之内（t L）则N（t-L）=0 。在这个初发期间内，方程（1） 给出的发病人数呈指数增长N（t）=N （ 0tL） （2）当Lt2L的时候，N（t-L）这部分人就已经没有传播能力了，因此我们推算出了下列模型N（t）= N（1+K） （t-L）K（1+K） （Lt2L） （3）当2Lt3L的时候又有下列模型N（t-L
5、） （2L60时，在这之前失去传播能力的只占了少部分，因此规定当t60时也可用N（t）= NN（t-L）的模型。K的值其实是一个变量，它每天的值都在发生变化。疫情刚开始的时候，K的值大，原因可能有刚可能是政府部门还没有足够重视起来，人们也还没有重视，医疗部门也还没有比较好的设备，医生们对病情也还没有很了解，技术上可能也还有不足。但随着病情的日益加重，来自各个方面的重视程度都有很大的提高，这是K的值就比较小了。在此模型中，我们认为 感染率（K）在 数值上与病例的增长率是相等的，疫情患者他传播在传播给健康人的时候，健康人他可能是带病毒了，但健康热处于潜伏期状态，据“全国“非典” 科技攻关组公布七大
6、科研进展”与于2003-06-03日报道中指出潜伏期患者传染的可能很小。有关部门对非典暴发过程中两例传播链进行了细致的调查和分析，这两个案例中共追查到潜伏期密切接触者158人，无一人死亡。因此我们在模型中说的感染率只为疫情患者传染给他人，而且他人发病，若他人不发病则不为感染率。增长率在数值上即为感染率。我们对全国所提供的所有数据中的已确诊病例累计进行了分析计算，得出感染率K的变化数据并画出了曲线图。如图4-1所示：图4-1K（感染率）是一条跟t的值有关的曲线，我们通过回归法K的公式为：K = 7E-13t- 4E-10t+ 8E-08t- 1E-05t+ 0.0006t- 0.0191t +0
7、.2325 （5）图4-1中R=0.6988为曲线回归的均方差，可见存在的误差并不大。t为疫情流行的天数。4.4模型检验通过该公式可预测疫情开始时或以后的累计病人总数。例如 要预测某一天病人的累计总数，将时间t的天数代入方程（5）即可求得K（感染率）的大小，因为L的值定在20天，所以当0t20时，将K代入（2）；当20t40时，将K代入（3）；当40t60时，将K代入（4）。当t=10时，我们根据方程（5），可求得K=0.0923,我们再将K=0.0923代入（2）得到N8。当t=50时，我们根据方程（5），可求得K=0.0614,我们再将K=0.0614代入（2）得到N308。这与实际给出的
8、数据非常接近。可以说明我们的模型是一个比较能够预测以及能为预防和控制提供信息的模型。4.5模型的应用与推广此模型可以作为预测以及能为预防和控制提供可靠、足够的信息的模型。4.6与指数模型的比较1）我们对不同阶段的疫情的计算和预测建立了不同的模型，这样来分析比附件1所提供的早期模型更加的精确。2）对感染率K求出了方程，可以知道每一天的疫情感染率，可以更加有效的计算与预测有关数据。3）该模型实用性更强，能更加准确的反映实情。5 建立模型的关键和困难建立模型的关键在于对模型进行动态的分析，当传染病发展到一定阶段在政府的控传染率下降。此时还用之前的误差会很大。在建立模型过程中有以下几个方面的困难：1）
9、对不同地区SARS的卫生知识的宣传的多少的不同，K的值就不一样；2）对某一地区的不同地方的强化管理也不一样（如公交、商场、餐厅、娱乐场所等），K的值也就不一样；3）还有保护工具的使用、建筑物的通风条件、居住的卫生条件等等的不同，都会有有不同的K的取值。6对于卫生部门采取的措施的评价对于卫生部门提前或延后5天采取严格的隔离措施的影响，我们可以建立下面的模型进行辅助分析估计：1） 模型参数定义：S（t）t时刻易感人群总数I（t） t时刻出现的新增患者 患者从患病起经过时间，仍为患者的概率患者距发病时间，具有传染性的概率患者与易感人群接触率近断时间的医学研究表明，从正式发病到治愈一般需714天或更长
10、时间，假定平均治愈时间为12天。2） 基本条件假设：新患者出现的数量与现有患者的数量成正比，也与现有易感者的数量成正比，即发病率是患者人数和易感者人数的双线性函数。由基本假设条件可得：S（t+1）=S（t）-I（t+1）I（t+1）= S（t）经整理后得：S（t+1）=S（t）-S（t+1）=S（t）（1-S（t+1）/S（t）=（1- （5）虽然不能具体知道的数值，那么我可以根据较为理想的均匀平均递减概率参数，可得下表：12345678910111211/1210/129/128/127/126/125/124/123/122/121/12表6-1如果病人发病后5天才开始隔离，并且的值在疫情
11、初期又较大的话，那么由上表可知病人已经分别以11/12、10/12、9/12、8/12、7/12的大概率在社会上与易感人群接触和传染。由（5）式得:1 （6）也就是S（t+1） 7/12。当处于潜伏期时，传染性几乎为0，因而同理我们有理由相信： （7）即S（t+1）近似于S（t）。所以，如果在病人发病前提前5天隔离的话，新增病人数将变得很小。7 建立传染病数学模型的重要性随着卫生设施的改善、医疗水平的提高以及人类文明的不断发展，诸如霍乱、天花等曾经肆虐全球的传染病已经得到有效的控制，但是在世界的某些地区，特别是贫穷的发展中国家，还不时出现传染病流行的情况，与次同时，一些鲜为人知的险恶传染病则跨
12、国越界在既包括发达国家也包括发展中国家的更大范围内蔓延。一直以来，建立传染病的数学模型来描述传染病的传播过程，分析受感染人数的变化规律，预报传染病高潮的到来等等，有着重要的作用。以最近突发性的险恶传染病SARS为例。从2002年11月16日在中国广东佛山市首例发生家族聚集性发病至2003年5月，疾病呈迅速蔓延趋势。目前全世界30多个国家和 地区有病例报告。中国大陆、香港和台湾发病人数占全球的90%以上。世界卫生组织（WHO）总干事Brundtlard博士指出，SARS已威胁到全球人类的健康。由于目前对SARS尚无可靠的病理学诊断，所以只能根据医疗卫生部门提供的可靠数据统计资料，建立模型来描述S
13、ARS病毒的宏观传播过程，有助于从量的方面来分析受感染人数的变化趋势，掌握SARS的流行规律，从而及时对疫情进行控制，提供科学的数据，认清传染的基本要素，为防病提供必要的依据 。例如，5月8日，西安交通大学医学院紧急启动“建立非典流行趋势预测与控制策略数学模型”研究项目。于5月19日初步完成了第一批成果，这一数学模型利用实际数据拟合参数，并对全国和北京、山西等地的疫情进行了计算仿真。结果指出，将患者及时隔离对于抗击非典至关重要。分析报告说，就全国而论，若非典病人延迟隔离1天，就医人数将增加1000人左右，推迟两天约增加2100人左右；若外界输入1000人中包含一个病人和一个潜伏病人，将增加患病
14、人数100人左右；若4月21日以后，政府未采取隔离措施，则高峰期病人人数将达60万人。同时美国科学杂志网站5月23日发表的两份最新研究报告显示，如果对非典采取严密的公共卫生防治措施，这种新型疾病是能够得到控制的。而采取这种措施需要有一个预见性，这就需要人们通过模型的建立对SARS的发病周期、发病人数的变化趋势、疑似人数的变化趋势等来分析和预测。并为政府和医疗卫生部门进行决策和资料调配提供直接的服务，为相关的研究部门提供科学的数据。SARS作为新发传染病之一，虽有着其特殊性，但也符合一般传染病的传播规律。从SARS对人民身体健康造成严重危害可以看出及时对传染病建立模型并进行分析和预测对人类的生命
15、健康有着至关重要的作用。8 参考文献1 姜启源 数学模型 高等教育出版社 1993.82 云舟工作室 数学建模基础教程 人民邮电出版社 2001.73 刘双等 SARS临床病例及影像学分析 中国医药科技出版社2003.5附录A疫情的数据表日期已确诊病例累计现有疑似病例死亡累计治愈出院累计4月20日33940218334月21日48261025434月22日58866628464月23日69378235554月24日77486339644月25日87795442734月26日988109348764月27日1114125556784月28日1199127559 4月29日134713586683
16、4月30日1440140875905月1日15531415821005月2日16361468911095月3日17411493961155月4日180315371185月5日189715101031215月6日196015231071345月7日204915141101415月8日213614861121525月9日217714251141685月10日222713971161755月11日226514111201865月12日230413781292085月13日234713382445月14日237013081392525月15日238813171402575月16日24051265273
17、5月17日242012501453075月18日24341473325月19日243712491503495月20日244412251543955月21日12211564475月22日245612051585285月23日246511791605825月24日249011341636675月25日249911051677045月26日250410697475月27日251210051728285月28日25149418665月29日25178031769285月30日252076017710065月31日252118110876月16日19020536月17日21206月18日19121546
18、月19日21716月20日21896月21日22316月22日22576月23日22776月1日252273911246月2日73411576月3日72411896月4日71812636月5日71613216月6日71318314036月8日55018415436月9日45116536月10日35117476月13日7118719446月14日18919946月15日20156月7日252366814466月11日18216月12日1551876附录Bfunction E=BJ（theta） format long;x=;C= %出院人数 ;A= %累计个案 B=;死亡人数 y=A-B-C;s
19、0=theta（4）;n0=theta（5）;t,s=ode45（feidian,0:64,s0 n0 339 51,theta）;E=sum（y-s（:,3）.2+（B+C-s（:,4）.2）%+（C-s（:,5）.2） X,FVAL,EDITFLAG=fminsearch（bj,0.0002 0.2 0.05 1800 180）;options=odeset（reltol,10（-14）,abstol,10（-15）;theta=X IMIN=FVAL y=zeros（65,1）;B= %死亡人数 100,theta（4） theta（5） 339 51,options,theta）;plot（t,s（:,3） hold on T=0:64;plot（T,y,o） xlabel（时间）;yla
