2019年B题露天矿生产的车辆安排南昌大学全国一等奖.doc

1、C题 SARS的传播摘要本文首先采用抽样检测法对SARS早期的模型的合理性及实用性进行了评 价，然后我们通过对传染病的共性及SARS的特性的分析。得出三个基本假设并 且把人群理想化为三类(S类，I类，R类)，建立起基本的SIR模型，再对SIR 模型中三类人群间的相互转化关系的分析，结合马氏链得岀三种人群间变化率的 矩阵T,由于SARS的特性,可知,SIR模型中的两个参数a(t),b(t)是以时间为 变量的函数。我们根据北京疫情的数据，通过多项式的数据拟合法分别得 a(t),b(t)的表达式，我们把a(t),b(t)及T结合，从而建立出模型。由于医疗 条件的逐步改善，必会硏制出其疫苗。于是我们在

2、不改变人群分类的情况下，增 加了一个系数c, (c表示疫苗日成功接种率，曲于在疫情期间，疫苗未能及时改 良，故c为常数。)进一步完善了我们的模型。本文利用数学软件(Mathematica, Matlab)很好的实现了模型运算，并结合 实际数据得岀了各类人群与时间的关系图。从图中可以很好的反映出各类人群的 变化规律，它们的变化规律与实际变化相吻合，从而证明了我们的模型基本符合 要求。1一问题的提出严重急性呼吸道综合症，简称SARS,是21世纪第一个在世界范围内传播的 传染病。它对全球的经济和生活造成巨大的破坏，尽管H前疫情已得到控制，但 对这种新冠状病毒及其流行规律的研究还刚刚开始，因此，有必要

3、根据SARS流 行的特点，建立数学模型预测其传染，从而釆取措施预防和控制其发展。而建立 该模型我们要综合各方面的因素才能使模型合理化。二问题的分析通过分析北京，香港和广东三地的受感染人数的变化规律，我们就可以对不 同地区预测流行病的变化趋势提出以下模型假设。模型的假设：1 将人群分为三类易感染者人数（疑似病例人用s表示；病人数（已受感染者，即确疹者）：用I表示；移出者人数（包括“被治愈者”和“死亡者”）：用R表示2 该地区人口不流动，疫情阶段无病原的输入和输出，设最初易感染者人 数为N,此时I, R均为0。3 被隔离人群完全断绝与外界接触，不再具有传染性。三模型的分析与建立问题一：对于附件1所

4、提供的早期模型即X （t）二N。（1+K）在疾病初期应该有它的 可参性。因为模型N （t）二门（1+K）是符合指数增长规律的。而对于这种传播 疾病，由于社会来不急防备以及群众的不重视，使得SARS疫情基本上呈自然规 律增长，这与香港的实际疫情拟合图基本一致。通过对香港和北京高峰前期的病 例与天数的对应关系计算出的数据与实际基本吻合。见下图表：香港(K=0 16204)病例数42080320实际天数10203040计算天数9.2319. 9529. 1838.41实际天数-计算天数0. 770. 050. 821.59北京（K二0. 13913）病例数13441 10470实际天数2030405

5、0计算天数19. 6929. 0537. 9347. 23实际天数-计 算天数0.310. 952.072. 77山上图可知在高峰期前即疾病早期基本符合指数增长规律，同时，我们也可以看到天数越接近高峰期，实际天数与计算天数就相差越大，这说明这个模型特 别适用于早期初，但从上表中的数据差来看，它们相差并不大，对整个早期阶段 也有较好的预测性。因此，该模型具有一定的合理性和实用性。问题二:模型的建立建立SIR模型川易感染者，感染者，移出者之和是个恒量即N二S+I+R。由于病人康复后具有 终生免疫力，人与人之间有相同的接触率。最终山如下两种假设决定状态之间的 转变率：(1)感染者的增长率是和感染者I

6、与易受感染者S的乘积成正比的；(2) 感染者I到移出者R的变化率是与感染者I成正比。基于以上假设得出模型的微 分方程：cIS “=aSIdt = aSIbI dtdR=bldt其中a, b都是以时间为变量的参数，a(t)为日感染率，b(t)为日移出率，于是 我们可以得出三类人的转换状态图：1SaXI.根据Markov Chain2理论，我们得出一个矩阵T:1 一 ，00T 二at1一$00b.1其中J就是当天易感染者S变为感染者I的日感染率,b就是当天感染者I转变为移出者R的移出率.(af, b分别为a(t), b(t)某一天的值)设初始值X产N, 0, 0,于是我们可以由X与T的转置矩阵相乘

7、，相乘一次得 到笫一天的易感染者，确诊病人及排除者的人数，再由该人数与T的转置矩阵相 乘一次得到笫二天各类人群的数H,依此类推，我们可以得到第t天的各类人群 的数目，于是我们可以得出任何一天各类人群的数目的初步模型即：X/X 严 T由于T是由a, b确定的，所以要建立模型还必须求出a(t),b(t).于是我们参照附件2中的数据，利用Matlab的多项式数据拟合，可以得到 它们相应的多项式。下面进一步讲拟合过程：笫一步：对附件2中的疑似病例进行20次数据拟合，得出每天的疑似病例 数的多项式：f s (t)二238. 2523-80. 1594t+480. 6525t 2 -325. 9497t3

8、+110. 1561t4 -22. 3320t5 +2. 9772t6-0. 275t7 +0. 0181ts-0. 0009t9 (1)笫二步：对附件2中已确疹的病例数据进行运算，新确诊的病例数m(t)就 是把当天已确诊病例数的累积n(t)减去前一天确诊病例数的累积n(t-l)即： m(t)=n(t)-n(t-l)对m(t)进行20次数据拟合，得出关于新增病例与时间的多项式：f / (t) =407. 5716-569. 7285t+476. 2872t2-227. 1678t3 +67. 6374t4-13. 3813t5+1. 8488t6-0.10854t7+0. 0139 t8-0.

9、 0008 t9 (2)由第二步对f, (t)求导得出f；(t), f；(t)表示为每天已确诊病例的变化情况, 因为f5(o是当天疑似病例数，所以得出笫三步：我们把死亡病例和治愈者都归为同一组r类，于是把他们相加得到 每天的R类人群K (t)，再由笫二步的计算方法我们得到每天死亡和治愈的人数 1 (t)即1(t)=K (t) -K (t-1)然后对l(t)进行二次数据拟合，得岀关于每天死亡和治愈人数与时间的多项式: fK (t)=-0. 0323t2+2. 95235t-13. 1284(3)第四步：按照第二步的方法对第二步中的m(t)再进行一次二次拟合，得到 f7 (t)的一个二次多项式：f

10、, (t) =0. 056112-5. 7216t+141. 7628同理，由(3), (4)式可得 b(t)二以(t)/f/ (t)B|J：(前面(1), (2), (3) , (4)式的Mat lab拟合图及程序见附录I)通过上面的分析我们可以得到一个完整的模型:(D/(O四模型的改进由于我们逐步对SARS的研究和认识，那么在不久的将来预防它的疫苗也将 出现显然医疗卫生部门就会对群众进行疫苗接种,于是我们的模型能够通过接 种疫苗而进行改进，但是没有必要增加其他人群种类，我们可以简单而直接的把 接种疫苗的人群从S中分离到R中，所以我们假定一个系数c,它表示每天有c的 人成功接种疫苗.(例,c

11、=8%,表示每天有8%的人成功接种疫苗)，我们前面的T矩 阵可以改为：1 一 终 一 c 00T= at 1一0cht1X/X/T上面的模型就可以改进为：(II)(/)=b=(其中C为常数)五 模型的求解与结果分析对于上面所建立的模型可以通过mathemdtica编程实现（程序见附录II）。 现已山北京市疫情的数据，通过模型求得的各类人群的变化规律图如下：（说明: 山于现阶段SARS的疫苗还没正式用于预防，所以只用模型（I）求解）I感染人数与时间的变化图2治愈和死亡人数与时间的变化图3通过分析上图我们可知：对于图1中疑似病例是随时间的增加而逐渐下降， 这与附件2中的疑似病例的变化规律基本一致；

12、图2中感染人数先随时间的增加 达到高峰之后平缓下降，这于附件1中北京日增病例走势相吻合；图3中治愈和 死亡人数随时间的变化而递增，这显然与实际相符。下面我们给出卫生部们提前或者延后五天采取严格的隔离措施后疫情传播 的变化图如下：山下图容易看出，图4中提前五天采取隔离措施，疑似病例下降速度比延后 五天的速度快；图5中提前五天感染人数远小于延后五天，而且提前五天时疫情 时间持续更短；图6中提前五天变化速度比延后五天来得更快。疑似病例与时间的关系图4感染人数与时间的关系图5治愈和死亡人数与时间关系图6综上所述，提前采取严格隔离措施比延后效果更好，也使得对社会的影响更低。四 模型的优点及发展方向1优点

13、：（1）我们的模型可以对SARS的所有阶段进行预测，而附件1的模 型只适用于早期；（2）我们的模型操作性强，结果精确，计算完全可由相应的数学软 件mathematics,mntlab的编程来实现。它可以根据疫情a,b的值会按设定 的函数自行修正，这也符合疾病随地域的变化而改变；（3）在我们的改进模型中加了一个系数c,它使得模型更精确，也 符合传染病的发展规律；（4）此模型可以直观的反映疾病发展规律。我们只要通过数学软件 画出图形，在图形上就可以显示出SARS在每天各类人群的分布情况。2模型的发展方向：要建立一个真正能够预测以及能为预防和控制提供可靠，足够的信息 的模型。要考虑以下儿点因素：（1

14、）必须准确掌握疾病的传播规律，获得更多的直观信息，例如传染率, 死亡率，治愈率，传染期限，潜伏期等：（2）考虑各种社会因素，包括医疗机构的医疗程度，人的免疫程度，群 众对SARS的认识，以及政府机构对SARS的宣传等以上条件可用参数P 来反映；（3）还要考虑城市间人口的流动，此因素用B表示。解决以上因素的困难：a）III于SARS突然性发生以及它的传播速度快的特性，这就使得在初期 很难预防；b）要对因素P, B考虑是一个非常复杂的过程。而要达到控制疾病的程 度，困难就在于社会经济的落后，技术的不先进，人们素质的不高，显然 要想改变这一现象不是一朝一夕的过程，也是一个变数很大的因素。问题三建立传

15、染病数学模型的重要性2002年末开始，全国部分地区遭受了非典型性肺炎的袭击，全国经济出现 了前所未有的滑坡，使人们的生命受到威胁。2003年初我国在毫无准备的情况 下受到“非典”的突击，当时全民上下在一片迷茫之中，不知道病原体的来历， 医务人员也无从下手，致使非典病原体不断扩散，使得感染人数不断增加。此时 我们才意识到要对这个新的传染病去做深入研究。研究发现SARS是一种独特性传染病，它具有传播速度快，传播途径广等特 点，正因为如此，在人们还没有认清SARS时，人们不知道如何采取措施对它进 行控制，至使在SARS疫情前期，人们损失惨重。经过一段时间的传播，人们对 这种病原体有了一定的研究，基本

16、掌握了它的传播途径。但是这时没有找到任何 一种疫苗，没有能力完全消除它，只能釆取适当可行的措施减少它的发病率。我 国在弄清情况后就采取了严格的防预措施。例如，在有发现疫情的地区，封锁各 出现疫情点，封锁学校等；在没有发现疫情的地区也采取了封闭措施，例如，切 断外来人员，部分大学也实施封校。全国上下掀起了一场抗击非典的战争，上至 政府下至平明口姓无一不重视SARS的传播。但对于SARS的传播规律还没正确的 了解，这些做法只是治标而不能治本。关键是我们还不能够通过某种途径来找岀 它的发展规律。对上所述我们知道，SARS之初之所以那么肆虐是因为人们对它还没有很好 的认识，更没有有效的预防和控制措施。

17、那么如何去研究解决这个难题就变得迫 在眉睫，而这个关键就在于能否建立起合理的传染病模型。假若一个好的传染病模型出来后，就可以利用它对传染病进行及时的预测， 并对该疾病提供有效的信息，也可以给相关卫生部门研究接种疫苗提供可靠的数 据，从而可以达到控制传染病高病发率的发生。除此之外，它还使得政府能够根 据预测的信息做好足够的人力，物力准备，做到从容以对。从而就给人民的社会 生活和经济生活提供了强力的保障，也给社会稳定增加了祛码。所以无论是从自 身利益出发,还是从国家的稳定发展考虑，建立传染病的数学模型是至关重要的, 也是摆在我们面前的一个迫切的任务。所以，建立合理的数学模型是预防和控制传染病至关重要的一步，也是不可 缺少的环节。参考文献1 姜启源，数学模型第三版，北京，高等教育出版社，2003年7月2 http:/pespmcl. vub. ac. be/ASC/MARKOV_CHAIN. html 2003 年 9 月 23 日论文点评：本文首先根据题LI给出的条件和要求，用抽样检验法对早期模型的合理性和 实用性进行了评价，其合理性在于符合指数增长规律，实用性在于对早期阶段的 疫悄