2003年全国数学建模优秀论文北京SARS的传播研究

1、小组成员姓名专业电话北京SARS的传播研究摘要SARS从2003年陆续传入，期间先后感染6000多人其中北京感染2847，我国给我过经济社会带来严重额的影响，为减少疾病的危害，提高人们对疾病的ARS的认识，疫情分析及对北京疫情走势的预测研究也变得尤为重要。为改善现状并提高人们对疾病的是SARS的认识，我们对北京市的SARS传播问题建立数学模型。关键词： SARS 人群分类 微分模型 整体拟合 1、问题重述1.1问题的背景严重急性呼吸综合征（Severe Acute Respiratory Syndromes），又称传染性非典型肺炎，简称SARS，是一种因感染SARS冠状病毒引起的新的呼吸系统传

2、染性疾病。主要通过近距离空气飞沫传播，以发热，头痛，肌肉酸痛，乏力，干咳少痰等为主要临床表现，严重者可出现呼吸窘迫。本病具有较强的传染性，在家庭和医院有显著的聚集现象。首发病例，也是全球首例。于2002年11月出现在广东佛山，并迅速形成流行态势 1.2问题的叙述现阶段北京SARS的传播正处于高峰期。由于人们对该种疾病的传播机理还不太清楚，因此引起人们的恐慌，它关系社会的稳定和经济的发展。因此对该问题的研究非常有必要，我们把人口分成四类，即:健康人S(t)SARS病人I(t)病人免疫（包括死亡）的人R(t)及疑似病人P(t)四类人，利用现有数据着重从四类人口中：把该传染病进行统计学分析，归纳出主

3、要特征通过假设，参数以及它们的相互联系，进行数据判定，数据假设，数据处理，数据分析，建立模型，数据总结等得出较为科学的SARS问题的分析，相关信息（见附件1、2、3）附件1SARS疫情分析及对北京走势的预测 附件2北京市疫情的数据附件3北京市接待海外游客人数附件4相关编程1.3问题的提出问题一：对附件1所提供的一个早期的模型，评价其合理性和实用性。问题二：建立自己的模型，说明为什么优于附件1中的模型，对于卫生部门所采取的措施做出评论，如：提前或延后5天采取严格的隔离措施，对疫情传播所造成的影响做出估计。问题三：影响SARS传播因素以及对SARS疫情期间北京疫区人群进行分类。2、符号说明N: 北

4、京人口基数S：健康易感人群在人口基数中所占的比例，是与时间t有关的函数 ：I:未被隔离的SARS病毒感染者在人口基数中所占的比例;e:已被隔离人群在人口基数中的比例（确诊、疑似）： ：隔离人群中SARS病毒感染者在人口基数中的比例，排除了e(t)例中没有非典的那部分人：r:人口基数中所占的比例:死亡人数在人口基数中所占的比例: :每个未被隔离的SARS病毒感染者每天接触传染的人数： :每个已被隔离SARS病毒感染者每天接触传染的人数 : :隔离率，随着时间推移会发生变化 : :表示每天被隔离的人数 :非典患者的治愈率:非典患者的死亡率:3、问题分析我们在科学、客观数据的基础之上，对未来几个月北

5、京SARS病人发展形势做出科学、合理、简单的假设。不计影响较小的因素。在假设的条件之下，建立数学模型。然后要参考附件1，预测从2003年至2004年的北京SARS病人数量的变化情况。可以考虑的到拟合的方法。针对问题一：1)合理性：附件1所给出的模型为：N(t)=NO(1+K)它是基于现实中的自然状态，描述出了SARS传染病最核心最本质的变化趋势。K的取值采用半模拟循环计算方法，发展趋势由K值的变化体现。该模型的优点在于简单，易行，方便对数据采用拟合处理和利用取对数求方差估计与实际数据的误差，说明了该模型所具有的合理性。2)实用性：任何具有传染性的疾病大致都是会经历“发展(快速蔓延)期一相对稳定

6、期一逐渐消亡期”这样的一个过程，附件一模型准确地体现出了这点，因此它具有普遍实用性。3)模型的缺陷：此模型把实际问题过于简单化了，有不合理的地方：(1)模型中的K的取值只能根据已经有的数据拟合，因此模型的精确度严重地依赖与所给数据的准确度。实际中，统计所给的数据本身就有一定误差，拟合一个本身就包含偏差的数据势必造成与现实规律更大的背离。我们根据图直观的看出，模型只能给出接近的前期发展趋势，后期拟合与实际曲线有相当误差。(2)模型本身不具有预测性，它的K值是由数据拟合决定的。如果背离题目本意，我们让K按照某种规律变化，预测发展趋势，其产生的误差是很大的。(图略)(3)随着时间的推移，社会中存在各

7、种控制的综合作用，用一个单纯笼统的K 的变化已很难刻画出复杂因素的影响，因为各种因素对SARS的影响不尽相同，有的可能抑制传播，有的则可能促进流行，致使模型的一致性在后期变差，误差越来越大。因此，至少应设为某种函数形式，引入一些参量因子进行考虑。(4)此模型单单刻画出了传染病的一般性，那么SARS和其它的传染病也就没什么本质上的区别了，缺乏对其SARS的，特征进行具体深入分析。针对问题二：对于附件1的模型建立优于它的模型。根据定义与假设列出相应的所需的方程组，由直接拟合推导各个参数存在较大的困难，因此采用整体拟合。再通过相应式子计算预测每日新增的隔离的SARS病毒感染者，整理相应的数据；最后预

8、测北京最终的累计感染非典人数。据此，在后标题“模型的评论与改进”中阐述对卫生部门锁采取的措施的评论。针对问题三通过对早期模型和实际情况的分析，我们认为影响SARS传播因素众多，大致可分为时域因素和地域因素。列举如下：(1)时域因素a媒体宣传：初期疫情较轻，媒体宣传强度很弱，导致民众的自我保护意识不足，容易感染；后期疫情较重，媒体宣传强度很大，民众的自我保护意识大大加强。b政府干预：初期疫情较轻，政府介入不足，后期疫情较重，政府加强干预（如：强行隔离，公共场所消毒等行为）。c认识程度：当一种新的传染病出现时，初期由于人们的认识程度不足，无法采取有效的预防和治疗措施，但随着研究的深入，认识程度会越

9、来越高。(2)地域因素a经济水平和医学水平：经济水平和医学水平高的地区的疫情控制情况会明显比水平低的地方好。b人口密度和人口流动：人口密度和人口流动大的城市若爆发传染病，疫情程度会比人口密度和人口流动小的城市大。c气候：SARS适合在春秋两季传播，且各城市的气候会疫情程度。综上我们认为一个较好的传染病传播模型因该具有如下功能： a能较好的描述疫情的大致走势。 b能较精确的给出关键时间（初期爆发时刻；中期稳定时刻；高峰期；0病例增长的时刻），以便政府和卫生部门针对不同作出及时而正确的措施。 c能给出描述疫情的指标，以便政府和卫生部门决定其各项工作的力度。4、模型的相关假设1、所获得的数据由权威部

10、门提供的全国疫情统计数据真实可信；2、将SARS所有传播途径都视为与病源的接触 3、在疾病传播期内所考察地区的总人数视为常数N，即认为本地区流入的人数与流出的人数均相等，时间以天为单位：4、假设每个病人单位时间有效接触人数r(所谓“有效接触”是指病人与健康这接触时，足以使健康者受到感染而成为病人)为常数：5、根据国家卫生部门资料可知处于潜伏期的SARS病人不具有传染性：6、根据国家卫生部资料，SARS康复者二度感染的概率为0，他们已经退出传染体系，因此将他们归为“退出者”。7、将人群分为五类： 1. 健康者：(易感人群) 2已被隔离的SARS感染患者 疑似病人：被隔离但没有确诊或排除的人员。确

11、诊SARS感染患者3. 未被隔离的SARS感染患者。 4. 治愈人群（不会传染SARS的人群）死亡人群治愈人群5、模型的建立与求解5.1 SARS传播模型的建立1 根据之前的定义和假设，我们知道每个未被隔离的SARS病毒感染者每天可以使S个健康易感染者，假设未被感染者的SARS感染人数为。因而知道每天被未隔离感染者感染的健康易感染者共有S。同理，每天被已隔离的感染者感染的健康者共有。此外还知道每天新增被隔离的感染SARS患者总数为：。综合可得出，每天新增的、未被隔离病人感染的数量总数为：+ （1）在每天新增被隔离的SARS感染者中减去每天治愈的人数和死亡的人数，可以得到有效的每天感染者的数量：

12、 （2）每日治愈人数 = （3）每日死亡人数 (4)并且，根据所占比例。应有： (5)上述式（1）（5）构成了求解SARS模型所需的常微分方程组，即SARS模型。5.2问题的求解。考虑到直接拟合推导各个比率参数, 存在很大的难度, 我们采取了整体拟合的策略, 以避免求解、时遇到相关数据缺乏所造成的困难。5.2.1预测北京最终的累计非典感染人数根据假设，每日新增的SARS病毒的感染者人数为 ，其中有来自新收治病例的20%，来自疑似病例中的 80%，根据已知数据有当天新隔离感染者为： =（当天确诊人数-前一天确诊人数）20%+（后一天确诊人数-当天确诊人数）80%以2003年4月21日作为第一天，

13、据此计算从4.21-5.17,每日新增隔离的病毒感染者人数如表1所示：日期123456789113.4105.285.898.6109.412393.2135.4104日期10111213141516171810989100.670.687.669.283.887.450.2日期19202122232425262748.240.438.842.2271917.215414.2根据这27组数据，对表1所得计算结果进行简单的描点（如图1折线所示）其分布与指数曲线 近似(如图1曲线所示)，所以对其进行整体回归拟合得到的表达式为： 通过拟合的曲线能够看到，在t=76的时候，新增隔离病毒感染者人数已近趋

14、于0。5.2.2 预测北京SARS疫情结束的大致时间当北京当天需要被隔离的人数降到0时，说明SARS病毒已经不再传播，疫情已经解除。每日处在隔离中的SARS感染的总人数为。 =当天已经确诊的病例-当天死亡人数-当天治愈出院人数t123456789113.4105.285.598.6109.412393.2135.4104103999320573743683796.499.269.585.697.411483.2132.492t10111213141516171810989100.670.687.669.283.887.450.2710963313711979543432937282.659.6

15、80.663.266.877.437.2t19202122232425262748.240.438.842.2271917.215.414.26711223685163422495511440.231.423.811.2-14611.2-1.4-23.8对其进行简单的描点，其图像大致和高斯函数图像较为吻合，如图2所示当=0时，也就是说在7月10或7月11左右，北京地区当天隔离的SARS感染人数为0，我们预测此时北京的非典疫情结束。6、模型评价及改进1、评价模型对北京地区中期的计算值与实际值基本吻合，说明该模型有一定的实用性。但对后期预测与后来的实际情况却有一定差距，而实际上，各地区的政策及人们

16、生活习惯各有所不同，因此用一个地区所获得的参数去预测另一地区，其结果只具有参考性，而不具备很强的可靠性。所以该模型的实用性有一定局限。经过计算，以4月20 日作为严格隔离开始时间，那么在外未被隔离的病人大概为1500人，如果严格隔离往后推5天，那么开始严格隔离时的在外患者人数约为2500人，那么这样的结果时SARS的疫情大大加深，当然，如果严格隔离提前5天起，那个疫情将会减轻。为了简化计算，我们的模型中没有考虑SARS病的潜伏期，SARS病毒的潜伏期一般为2-7天，并在潜伏期内，该病人不具有传染性，因此对北京最终SARS病毒感染者的总数没有太大影响，但它明显一定程度会影响高峰期和结束期，这是该

17、预测 模型需要改进的地方。模型采用微分方程本身就有一定的缺限，其计算结果的准确性、可靠性将受到限制，再加之数值解的不确定性，模型对长时间的预测有它的局限性。因时间限制模型没能更多考虑交叉分类进行。2、改进若能建立以随机偏微分方程组为基础的数学模型，将大大提高计算的准确性与可靠性，使得预测更加准确，但这样做将遇到模型求解，数据准确收集和数值求解的不精确性等诸多困难。参考文献1中华人民共和国卫生部网站，http：wwwm0hgovcn2百度百科：http：3王树禾著常微分方程模型与混沌M 4朱道元编著数学建模精品案例M东南大学出版社，1999年8月第1版5王兵团编著.数学建模基础。6吴建国编著.数

18、学建模案例精编。附录附件1：北京市疫情的数据( 据： )已确诊病例累计现有疑似病例死亡累计治愈出院累计当天退出数当天病人数当天病例退出率治愈率3394021833174310.0394430.076566482610254365200.0115380.0826925886662846166190.0258480.0743136937823555136840.0190060.0804097748633964127740.0155040.082687877954427398730.0103090.0836298810934876109900.0101010.07676811141255567831

19、0650.0028170.0732391199127559781212100.0099170.0644631347135866831612910.0123930.0642911440140875901713880.0122480.06484115531415821001814540.012380.06877616361468911091115410.0071380.0707331741149396115715920.0043970.07223618031537100118616790.0035740.07028189715101031211717360.0097930.069719601523

20、1071341018080.0055310.074115204915141101411318850.0068970.074801213614861121521819130.0094090.07945621771425114168919450.0046270.086375222713971161751519740.0075990.088652226514111201863119980.0155160.093093230413781292084120100.0203980.103483234713381342441319920.0065260.1224923701308139252619970.0

21、030050.126189238813171402571720080.0084660.127988240512651412733820060.0189430.136092242012501453072719820.0136230.154894243412501473322019580.0102150.169561243712491503495019450.0257070.179434244412251543955418950.0284960.208443244412211564478318530.0447920.24123245612051585285617790.0314780.296796

22、246511791605828817480.0503430.332952249011341636674116690.0245660.399641249911051677044416330.0269440.431108250410691687478515970.0532250.467752251210051728284115140.0270810.54689625149411758666314760.0426830.58672125178031769287914160.0557910.655367252076017710068513380.0635280.75186825217471811087

23、97512530.7781320.86751825217391902053672780.2410077.38489225217341902120352110.16587710.0473925217241912154171760.09659112.2386425217181912171181590.11320813.6540925217161912189421410.29787215.524822521713191223126990.26262622.535352521550191225720730.27397330.9178125214511912277-116354-21.53742.166

24、67252235118111243312170.0271160.92358325227118111573211840.0270270.9771962522418111897411520.0642361.0321182522318112635810780.0538031.171614252266818113218410200.0823531.295098252225718314031419360.1506411.498932252215518415431107950.1383651.940881252231841653966850.1401462.413139252251861747198589

25、0.3361632.966044252241871944523910.1329924.971867252231891994213390.0619475.882006252231892015-575319-1.802516.3166142523218314463788940.4228191.61745252321861821565160.1085273.52907252321871876-0.289883.118497附件3：北京市接待海外旅游人数（单位：万人）年1月 2月 3月 4月 5月 6月 7月 8月 9月 10月 11月 12月19971998199920002001200220039

26、.4 11.3 16.8 19.8 20.3 18.8 20.9 24.9 24.7 24.3 19.4 18.69.6 11.7 15.8 19.9 19.5 17.8 17.8 23.3 21.4 24.5 20.1 15.910.1 12.9 17.7 21.0 21.0 20.4 21.9 25.8 29.3 29.8 23.6 16.511.4 26.0 19.6 25.9 27.6 24.3 23.0 27.8 27.3 28.5 32.8 18.511.5 26.4 20.4 26.1 28.9 28.0 25.2 30.8 28.7 28.1 22.2 20.713.7 29.

27、7 23.1 28.9 29.0 27.4 26.0 32.2 31.4 32.6 29.2 22.915.4 17.1 23.5 11.6 1.78 2.61 8.8 16.2附4t123456789113.4105.285.598.6109.412393.2135.4104103999320573743683796.499.269.585.697.411483.2132.492t10111213141516171810989100.670.687.669.283.887.450.2710963313711979543432937282.659.680.663.266.877.437.2t1

28、9202122232425262748.240.438.842.2271917.215.414.26711223685163422495511440.231.423.811.2-14611.2-1.4-23.8113.4105.285.898.6109.4 123 93.2 135.4104日期10111213141516171810989100.670.687.669.283.887.450.2日期19202122232425262748.240.438.842.2271917.215414.2附5 x=1:1:27;y=113.4 , 105.2, 85.8 , 98.6 , 109.4

29、, 123 , 93.2 , 135.4 , 104 , 109 , 89 , 100.6 , 70.6 , 87.6 , 69.2 , 83.8 , 87.4 , 50.2 , 43.2 , 40.4 , 38.8 , 42.2 , 27 , 19 , 17.2 , 15.4 , 14.2;plot(x,y);clf;x=1:1:27;y=113.4 , 105.2, 85.8 , 98.6 , 109.4 , 123 , 93.2 , 135.4 , 104 , 109 , 89 , 100.6 , 70.6 , 87.6 , 69.2 , 83.8 , 87.4 , 50.2 , 43.2 , 40.4 , 38.8 , 42.2 , 27 ,