SARS传播的数学模型.docx

SARS传播的数学模型.docx

《SARS传播的数学模型.docx》由会员分享，可在线阅读，更多相关《SARS传播的数学模型.docx（25页珍藏版）》请在冰豆网上搜索。

SARS传播的数学模型

SARS传播的数学模型

SARS传播的数学模型

摘要

通过对题目附件1的SARS模型进行分析和评价，加深了对SARS的认识和了解。

根据传染病的传播特点，建立了关于SARS病人率和疑似病人率两个常微分方程模型。

以所给数据为基本依据，用Matlab软件进行数值计算，与图形模拟方法求得模型中的有关参数。

当λ1=1.5和λ2=1时，理论图形与实际图形有良好的吻合，分别得到了SARS病人率和疑似病人率比较符合实际数据的变化图，能正确地预测它们的发展趋势。

他们对于模型中的参数有非常强的灵感性，λ1的值作微小的改变对于整个疫情的发展有很大的影响，所以政府采取对SARS疫情的有关措施是完全正确的。

本文重点分析了关于SARS病人率的模型一，根据求得的参数，利用相轨线理论对结果加以分析并对整个疫情作出预测，并推论出SARS病人率关于t的表达式i（t）,然后提出了对传染病的控制方案，同时列举了具体方法，并论证了方法的合理性和可行性，用其它地区的数据对模型进行检验，说明模型的参数有区域性。

关键词：

SARS微分方程曲线拟合数学模型相轨线

本文首先分析评价了附件1中SARS传播的数学模型，指出该模型可以对疫情走势进行预测,但同时也存在一定缺点，第一，混淆了累计病例数与累计确诊人数的概念；第二，对参数的确定缺乏根据；第三，预测时借助了其他地区的参数，偏差较大.本文针对其缺点建立了一个比较完善的传播模型.该传播模型按政府开始控制的时刻分为控制前与控制后两个模型，两个模型均以潜伏期5天为周期，以一个周期为整体建立差分方程模型.再结合5月15日以前北京疫情的公开数据，配合不同的政府监控力度，对整个北京的SARS疫情状况进行了预测.预计政府的监控力度一直保持在5月10日－5月15日的水平上时，6月10日－6月15日北京将会无新增病例，最后累积病例数为2993.对卫生部门采取的措施进行了评价：

若提前或延后5天采取严格的隔离措施最后累计病例数分别为1300多与5200左右.进一步通过对人群的不同分类，建立了两个微分方程组，可分别预测出实际发病人数、不可控/可控带菌者人数与当天疑似病例数、累计确诊人数、不可控/可控带菌者人数及治愈、死亡人数，结合两者的信息就可以得到足够的信息量.但模型中的部分参数无法确定给模型求解带来困难.可以通过搜集更多的数据和资料加以解决.本文同时就外国来京旅游人数受SARS的影响，建立了模型，估算出4、5、6、7四个月中北京地区入境旅游人数比往年同期减少了94.8万人，旅游经济损在4.74亿美元至9.48亿美元之间.并预测出在2003年10月上旬，旅游人数将恢复到正常水平.最后给报纸写了一篇短文，说明了建立传染病数学模型的必要性与重要性.一、问题的提出公元2003年春天，一种叫SARS的病毒从天而降，降到人类赖以生存的星球，降到中国人的头上.SARS究竟是什么，它为什么会代给人类这么

[4]寿纪麟，《数学建模—方法与范例》，西安：

西安交通大学出版社，1993年[5]赵静但琪等，《数学建模与数学实验》，北京：

高等教育出版社，2000年[6]姜启源，《数学模型》，北京：

高等教育出版社，1993年附件附件一给当地报刊的短文数学枪模型弹瞄准病毒一二三2002年末，一种叫SARS的病毒看中了我们耐以生存的蓝色星球，决定在这里繁衍生息。

它们站在这颗星球的主人－人类的头上大声叫嚣：

“我是病毒我怕谁？

”SARS确实厉害，在人类可以为自己在科技、医学等各方面的巨大成就而骄傲的时候，它给足了人类颜色。

由于我们品种齐全的“杀虫剂”中没有“灭SARS牌”的药剂，所以它们把整个世界搞得天昏地暗所做的工作，仅仅是轻轻推倒多米诺骨牌的第一张。

在这种危急时刻，人类的政治、科技领袖们带领大家扬着科学的旗帜，同舟共济，赢得了非典狙击战的巨大胜利。

而这其中，传染病数学模型的建立起到了至关重要的作用。

首先，正确的传染病模型客观反应了病毒的传播规律，为预测和控制传染病提供了必要前提；为我们的生产生活提供了科学的依据与正确的指导方向。

例如：

中科院非典模型研究小组建立的非典预测模型，就为我国的非典控制和消除起到了十分有效的指导作用！

其实，我们也正是通过这些传染病模型，积极阻断非典的传播途径，加大监控力度，以致最终取得这场战争的胜利。

其次，传染病模型并非仅仅可以描述病毒的传播规律。

一种传染病的流行，还会使国民经济，政治，文化受到冲击。

例如SARS对旅游业的负面影响巨大，有了SARS传播的数学模型，就可以预测SARS对旅游业的影响，估算损失。

例如：

通过SARS传播的数学模型，我们建立了非典对北京市接待海外旅游者人数影响的数学模型，并利用建立的模型预测出十月上旬左右北京市的海外旅游者的接客数量将会恢复到正常水平；在非典时期该旅游行业的总损失在474050000美元到948100000美元之间！

海外旅游接待人数将会减少948100人！

通过现有资料的调查，可以说明我们的预测结果是合理可行的。

当然亦可以分析SARS对其它方面的影响。

这样，各相关部门可以根据模型所提供的病情预测，更好地制定在传染病流行期间，其防止措施和生产计划，使疾病造成的损失降到最小！

这足以说明建立传染病模型的重要性！

知己知彼是我们在这场与SARS的斗争中取胜的法宝。

传染病数学模型的建立，不仅是我们在这场战役中使用的强大武器也是一支指挥棒。

一个正确的数学模型，能客观反应事物内在规律，有科学的预见作用。

而在进行了科学预测之后，我们才能对症下药，对传染病的传播加以正确控制，并达到事半功倍的效果。

马克思指出：

一种科学只有成功地应用数学时，才算达到了真正完美的地步。

有了传染病数学模型这个强有力的科学武器做指导，再面目狰狞的传染病我们也有科学的决策去应对它！

数学枪模型弹，瞄准病毒，一二三，开枪！

附件二“SAS的传播”原模型的半模拟循环计数法程序N0=1;%北京的初始非典病例k=0.13913;%开始阶段的每个病人平均传染概率t=1:

59;%高峰前的非典传播时间N1

（1）=N0*k;%该循环统计疫情高峰前的累积病例fori=1:

59ifi<=20N（i）=N0*（1+k）^i;%控制前累计病例ifi>=2N1（i）=N（i）-N（i-1）;%每天新增病例endendifi==20N（i）=（N（i-1）-N0）*k+N（i-1）;N1（i）=N（i）-N（i-1）;endifi>20N（i）=（N（i-1）-N（i-20））*k+N（i-1）;N1（i）=N（i）-N（i-1）;endendk1=0.0273;%10天过渡后的每个病人平均传染概率%在10天的过渡期，分四次调整k值

k2=k*5/12;ik=k-k2forj=1:

3i=i+1;N（i）=（N（i-1）-N（i-20））*k+N（i-1）;N1（i）=N（i）-N（i-1）;endk2=k/4;k=k-k2forj=4:

6i=i+1;N（i）=（N（i-1）-N（i-20））*k+N（i-1）;N1（i）=N（i）-N（i-1）;endk2=k/2.75;k=k-k2forj=7:

10i=i+1;N（i）=（N（i-1）-N（i-20））*k+N（i-1）;N1（i）=N（i）-N（i-1）;endik=k1%该循环求疫情消失后的累积病例whileN（i）-N（i-1）>1i=i+1;N（i）=（N（i-1）-N（i-20））*k1+N（i-1）;endN附件三非典对北京海外旅游接客数的影响的模型和预测程序%A1—A7为北京市接待海外旅游人数1997—2003年的数据A1=[9.411.316.819.820.318.820.924.924.724.319.418.6]'A2=[9.611.715.819.919.517.817.823.321.424.5

20.115.9]'A3=[10.112.917.721.021.020.421.925.829.329.823.616.5]'A4=[11.426.019.625.927.624.323.027.827.328.532.818.5]'A5=[11.526.420.426.128.928.025.230.828.728.122.220.7]'A6=[13.729.723.128.929.027.426.032.231.432.629.222.9]'A7=[15.417.123.511.61.782.618.816.20000];t=1:

12x=[A1,A2,A3,A4];%采用时间序列法求没有非典的情况2003年的实际预测值fori=1:

12m（i）=mean（[A4（i）,A5（i）,A6（i）]）;endm1=（15.4-12.2）+m;%没有非典的情况2003年的实际预测值plot（[1:

12],m1）plot（[1:

84],[A1'A2'A3'A4'A5'A6'm1],'*',[1:

84],[A1'A2'A3'A4'A5'A6'm1]）;x1=[15.417.123.511.61.782.618.816.2];c=（m1（[1:

8]）-x1）./m1（[1:

8]）%这里，可以通过第3月的真实值检验该模型的预测精确性c（3）=[];[p1,s1]=polyfit（[4:

7],c（[4:

7]）,2）;%用二次曲线进行最小二乘拟合，预测在实际情况下的数据Y1=polyconf（p1,[4:

10],s1）;%预测值plot（[4:

10],Y1,'r'）;holdon;c=[c,[0.21462232228424-0.14273284947445-0.56883671831862]];plot（[1:

10],c,'*',[1:

4],c（[1:

4]））;holdon;plot（[1:

10],zeros（1,10）,'k'）

一、问题的提出

SARS俗称非典型肺炎，是21世纪第一个在世界范围内传播的传染病。

我国作为发展中大国深受其害：

SARS的爆发和蔓延给我国的经济发展和人民生活带来了很大影响。

在党和政府的统一领导下，全国人民与SARS顽强抗争，取得了可喜的阶段性胜利，并从中得到了许多重要的经验和教训，认识到在没有找出真正病因和有效治愈方法前，政府采取的强制性政策对抑制SARS自然发展最有效办法。

而本题的目的就是要建立一个适当的模型对SARS传播规律进行定量地分析、研究，为预测和控制SARS蔓延提供可靠、足够的信息，无论对现在还是将来都有其重要的现实意义。

二、模型的假设

1．地总人数N可视为常数，即流入人口等于流出人口。

2．据人口所处的健康状态，将人群分为：

健康者，SARS病人，退出者（被治愈者、免疫者和死亡者）。

3．在政府的强制措施下，人口基本不流动，故无病源的流入和流出，避免了交叉感染，降低了感染基数。

4．隔离的人断绝了与外界的联系，不具有传染性。

5．SARS康复者二度感染的概率为0。

6．国家完善了监控手段，加强了对SARS病毒监控的力度，故可假设所有感染SARS病毒的人群都进入了SARS病人类和疑似类。

7．由于对SARS病原体的研究不够深入，无有效药物可以使人体免疫，同时SARS病毒感染后，大量繁殖，破坏免疫系统，故不可免疫。

三、模型的建立

（一）参数的设定和符号说明

s（t）：

t时刻健康者在总体人群中的比例

i（t）：

t时刻SARS病人在总体人群中