高教社杯全国大学生数学建模竞赛题目[001]

1、个人资料整理，仅供个人学习使用1 / 72003 高教社杯全国大学生数学建模竞赛题目 （请先阅读 对论文格式的统一要求”）A 题 SARS 的传播SARS（ Severe Acute Respiratory Syndrome，严重急性呼吸道综合症，俗称：非典型肺炎） 是 21 世纪第一个 在世界范围内传播的传染病。SARS 的爆发和蔓延给我国的经济发展和人民生活带来了很大影响，我们从中得到了许多重要的经验和教训，认识到定量地研究传染病的传播规律、为预测和控制传染病蔓延创造条 件的重要性。请你们对 SARS 的传播建立数学模型，具体要求如下：矚慫润厲钐瘗睞枥庑赖賃軔。（1）对附件 1 所提供的一

2、个早期的模型，评价其合理性和实用性。（2） 建立你们自己的模型，说明为什么优于附件1 中的模型；特别要说明怎样才能建立一个真正能够预测以及能为预防和控制提供可靠、足够的信息的模型，这样做的困难在哪里？对于卫生部门所采取的措施做出评论，如：提前或延后5 天采取严格的隔离措施，对疫情传播所造成的影响做出估计。附件2 提供的数据供参考。聞創沟燴鐺險爱氇谴净祸測。（3） 收集 SARS 对经济某个方面影响的数据，建立相应的数学模型并进行预测。附件3 提供的数据供参考。（4）给当地报刊写一篇通俗短文，说明建立传染病数学模型的重要性。附件 1 :SARS 疫情分析及对北京疫情走势的预测2003 年 5 月

3、 8 日在病例数比较多的地区，用数理模型作分析有一定意义。前几天，XXX 老师用解析公式分析了北京SARS疫情前期的走势。在此基础上，我们加入了每个病人可以传染他人的期限（由于被严格隔离、治愈、死亡 等），并考虑在不同阶段社会条件下传染概率的变化，然后先分析香港和广东的情况以获得比较合理的参数，最后初步预测北京的疫情走势。希望这种分析能对认识疫情，安排后续的工作生活有帮助。残骛楼諍锩瀨濟溆塹籟婭骤。1 模型与参数假定初始时刻的病例数为 N0，平均每病人每天可传染 K 个人（K 一般为小数），平均每个病人可以直接感 染他人的时间为 L 天。则在 L 天之内，病例数目的增长随时间 t（单位天）的关

4、系是：酽锕极額閉镇桧猪訣锥顧荭。N（ t）= N0 （1+K）t如果不考虑对传染期的限制，则病例数将按照指数规律增长。考虑传染期限L 的作用后，变化将显著偏离指数律，增长速度会放慢。我们采用半模拟循环计算的办法，把到达L 天的病例从可以引发直接传染的基数中去掉。彈贸摄尔霁毙攬砖卤庑诒尔。参数 K 和 L 具有比较明显的实际意义。L 可理解为平均每个病人在被发现前后可以造成直接传染的期限，在此期限后他失去传染作用，可能的原因是被严格隔离、病愈不再传染或死去等等。从原理上讲，这个参 数主要与医疗机构隔离病人的时机和隔离的严格程度有关，只有医疗机构能有效缩短这个参数。但我们分 析广东、香港、北京现有

5、的数据后发现，不论对于疫情的爆发阶段，还是疫情的控制阶段，这个参数都不能用得太小，否则无法描写好各阶段的数据。该参数放在15-25 之间比较好，为了简单我们把它固定在20（天）上这个值有一定统计上的意义，至于有没有医学上的解释，需要其他专家分析。謀荞抟箧飆鐸怼类蒋薔點鉍。参数 K 显然代表某种社会环境下一个病人传染他人的平均概率，与全社会的警觉程度、政府和公众采取的各种措施有关。在疾病初发期，社会来不及防备，此时K 值比较大。为了简单起见，我们从开始至到高峰期间均采用同样的 K 值（从拟合这一阶段的数据定出），即假定这阶段社会的防范程度都比较低，感染率 比较高。到达高峰期后，我们在10 天的范

6、围内逐步调整 K 值到比较小，然后保持不变，拟合其后在控制阶段的全部数据，即认为社会在经过短期的剧烈调整之后，进入一个对疫情控制较好的常态。显然，如果 疫情出现失控或反复的状态，则K 值需要做更多的调整。厦礴恳蹒骈時盡继價骚卺癩。2 计算结果2.1 对香港疫情的计算和分析。香港的数据相对比较完整准确。但在初期，由于诊断标准等不确切，在3月 17 日之前，没有找到严格公布的数据。我们以报道的2 月 15 日作为发现第一例病人的起点，2 月 27日从报道推断为 7 例。3 月 17 日后则都是正式公布的数据。累积病例数在图 1 中用三角形表示。 我们然后用上述方法计算。4 月 1 日前后（从起点起

7、 45 天左右）是疫情高峰时期，在此之前我们取K=0.16204。此后的 10 天，根据数据的变化将 K 逐步调到 0.0273,然后保持 0.0273 算出后面控制期的结果。 短期内 K 调 整的幅度很大，反映社会的变化比较大。图中实心方黑点是计算的累积病例数。从计算累积病例数，很容 易算出每天新增病例数（当然只反映走向，实际状况有很大涨落）。可以看出，香港疫情从起始到高峰大2 / 7个人资料整理，仅供个人学习使用约 45 天，从高峰回落到 1/10 以下（每天几个病例）大约 40 天（5 月上中旬），到基本没有病例还要再经 过近一个月（至 U 6 月上中旬）。茕桢广鳓鯡选块网羈泪镀齐。22

8、对广东疫情的计算和分析。广东的起点是02 年 11 月 16 日 到今年 2 月下旬达到高峰，经过了约100天。在今年 2 月 10 日以前的数据查不到，分析比较困难。总体上看，广东持续的时间比香港长得多，但 累积的总病例数却少一些，这反映出广东的爆发和高峰都不强烈。但广东的回落也比较慢。从2 月下旬高峰期到现在经过了约 70 天，还维持着每天 10 来个新增病例，而同样过程香港只用了约40 天。这种缓慢上升和下降的过程也反映到K 值上。比较好的拟合结果是，在高峰期之前（t 101 天），K=0.0892 ;在随后的 10 天逐步调整到 0.031。用这组参数算出的后期日增病例数比实际公布的偏

9、小，说明实际上降低得更 慢。这种情况与疫情的社会控制状况有没有什么关系，需要更仔细的分析。鹅娅尽損鹌惨歷茏鴛賴縈诘。2.3 对北京疫情的分析与预测。北京的病例起点定在 3 月 1 日,经过大约 59 天在 4 月 29 日左右达到高峰。我们通过拟合起点和 4 月 20 日以后的数据定出高峰期以前的K=0.13913。这个值比香港的 0.16204 来得低，说明北京初期的爆发程度不如香港，但遗憾的是上升时间持续了近60 天，而香港是 45 天，这就造成了累积病例数大大超过香港。从图2 中还看出 4 月 20 日以前公布的数据大大低于计算值。而我们从对香港、广东情况的计算中，知道疫情前期我们的计算

10、还是比较可行的。从而可以大致判断出北京前期实际的病例 数。图中的公布数据截止到5 月 7 日（从起点起 67 天），其后的计算采用的是香港情况下获得的参数。按这种估算，北京最终累积病例数将达到3100 多。籟丛妈羥为贍债蛏练淨槠挞。3 / 7图 1 对香港疫情的拟合个人资料整理，仅供个人学习使用4 / 7图 2 对北京疫情的分析图 3 是计算的日增病例数。后期下降得较快的实心方黑点是采用香港参数获得的。这就是说，如果北京的疫情控制与香港相当或更好的话，就可以在高峰期后的40 天（从起点起 100 天）左右，即 6 月上中旬下降到日增几例。然后再经过约一个月，即7 月上中旬达到日增 0 病例。但

11、如果北京的新病例下降速度与广东类似的话，则要再多花至少一个月，才能达到上述的效果，且累积总病例数会到3800 左右。至于什么原因造成香港下降速度快而广东下降速度慢，需要有关方面作具体分析。預頌圣鉉儐歲龈讶骅籴買闥。II|1|1 rI110204060ao 100120口 3 厅 ill 厘型图 3 北京日增病例走势分析3 结论每个病人可以造成直接感染他人的期限平均在20 天左右，这个值在不同地区和不同疫情阶段似乎变化不大。病人的平均每天感染率与社会状况有关，在疫情爆发期较大，在疫情控制期要小很多。香港的初期爆 发情况比广东和北京都剧烈，但控制效果明显比较好。北京后期如果控制在香港后期的感染率水

12、平上，则 有望在 6 月上中旬下降到日增几例。然后再经过约一个月，即7 月上中旬达到日增 0 病例。而累积总病例数将达到 3100 多。但如果北京的新病例下降速度与广东类似的话，则要再多花至少一个月，才能达到上 述的效果，且累积总病例数会到3800 左右。渗釤呛俨匀谔鱉调硯錦鋇絨。附件 2 :北京市疫情的数据（据：）铙誅卧泻噦圣骋贶頂廡缝勵。4 月 20 日33940218334 月 21 日48261025434 月 22 日58866628464 月 23 日69378235554 月 24 日77486339644 月 25 日87795442734 月 26 日988109348764

13、 月 27 日1114125556784 月 28 日119912755978Q个人资料整理，仅供个人学习使用5 / 74 月 29 日1347135866834 月 30 日1440140875905 月 01 日15531415821005 月 02 日16361468911095 月 03 日17411493961155 月 04 日180315371001185 月 05 日189715101031215 月 06 日196015231071345 月 07 日204915141101415 月 08 日213614861121525 月 09 日217714251141685 月 1

14、0 日222713971161755 月 11 日226514111201865 月 12 日230413781292085 月 13 日234713381342445 月 14 日237013081392525 月 15 日238813171402575 月 16 日240512651412735 月 17 日242012501453075 月 18 日243412501473325 月 19 日243712491503495 月 20 日244412251543955 月 21 日244412211564475 月 22 日245612051585285 月 23 日24651179160

15、5825 月 24 日249011341636675 月 25 日249911051677045 月 26 日250410691687475 月 27 日251210051728285 月 28 日25149411758665 月 29 日25178031769285 月 30 日252076017710065 月 31 日252174718110876 月 01 日252273918111246 月 02 日25227341811157个人资料整理，仅供个人学习使用6 / 76 月 03 日252272418111896 月 04 日252271818112636 月 05 日2522716

16、18113216 月06日252271318314036 月 07 日252366818314466 月 08 日252255018415436 月 09 日252245118416536 月10日252235118617476 月11日252325718618216 月12日252315518718766 月 13 日25227118719446 月 14 日2522418919946 月 15 日2522318920156 月 16 日2521319020536 月 17 日2521519021206 月 18 日2521419121546 月 19 日2521319121716 月20日

17、2521319121896 月21日2521219122316 月22日2521219122576 月 23 日252121912277附件 3 :北京市接待海外旅游人数（单位：万人）年1 月月2 月 12月3 月4 月5 月6 月7 月8 月9 月10 月111997 9.411.316.819.820.318.820.924.924.724.319.41998 18.61999 9.611.715.819.919.517.817.823.321.424.520.12000 15.92001 10.112.917.721.021.020.421.925.829.329.823.62002 1

18、6.52003 11.426.019.625.927.624.323.027.827.328.532.818.511.526.420.426.128.928.025.230.828.728.122.220.713.729.723.128.929.027.426.032.231.432.629.222.915.417.123.511.61.782.618.816.2B题露天矿生产的车辆安排钢铁工业是国家工业的基础之一，铁矿是钢铁工业的主要原料基地。许多现代化铁矿是露天开采的，它的 生产主要是由电动铲车（以下简称电铲）装车、电动轮自卸卡车（以下简称卡车）运输来完成。提高这些 大型设备的利用率是增加

19、露天矿经济效益的首要任务。擁締凤袜备訊顎轮烂蔷報赢。露天矿里有若干个爆破生成的石料堆，每堆称为一个铲位，每个铲位已预先根据铁含量将石料分成矿石和7 / 7个人资料整理，仅供个人学习使用岩石。一般来说，平均铁含量不低于25%的为矿石，否则为岩石。每个铲位的矿石、岩石数量，以及矿石的平均铁含量（称为品位）都是已知的。每个铲位至多能安置一台电铲，电铲的平均装车时间为5 分钟。贓熱俣阃歲匱阊邺镓騷鯛汉。卸货地点（以下简称卸点）有卸矿石的矿石漏、2 个铁路倒装场（以下简称倒装场）和卸岩石的岩石漏、岩场等，每个卸点都有各自的产量要求。从保护国家资源的角度及矿山的经济效益考虑，应该尽量把矿石 按矿石卸点需要

20、的铁含量（假设要求都为29.5% 1%，称为品位限制）搭配起来送到卸点，搭配的量在一个班次（8 小时）内满足品位限制即可。从长远看，卸点可以移动，但一个班次内不变。卡车的平均卸车 时间为 3 分钟。坛搏乡囂忏蒌鍥铃氈淚跻馱。所用卡车载重量为 154 吨，平均时速 28 国。卡车的耗油量很大，每个班次每台车消耗近1 吨柴油。发动机点火时需要消耗相当多的电瓶能量，故一个班次中只在开始工作时点火一次。卡车在等待时所耗费的 能量也是相当可观的，原则上在安排时不应发生卡车等待的情况。电铲和卸点都不能同时为两辆及两辆以 上卡车服务。卡车每次都是满载运输。蜡變黲癟報伥铉锚鈰赘籜葦。每个铲位到每个卸点的道路都

21、是专用的宽60 的双向车道，不会出现堵车现象，每段道路的里程都是已知的。一个班次的生产计划应该包含以下内容：出动几台电铲，分别在哪些铲位上；出动几辆卡车，分别在哪些 路线上各运输多少次（因为随机因素影响，装卸时间与运输时间都不精确，所以排时计划无效，只求出各 条路线上的卡车数及安排即可）。一个合格的计划要在卡车不等待条件下满足产量和质量（品位）要求， 而一个好的计划还应该考虑下面两条原则之一：買鯛鴯譖昙膚遙闫撷凄届嬌。1.总运量（吨公里）最小，同时出动最少的卡车，从而运输成本最小；2利用现有车辆运输，获得最大的产量（岩石产量优先；在产量相同的情况下，取总运量最小的解）。请你就两条原则分别建立数学模型，并给出一个班次生产计划的快速算法。针对下面的实例，给出具体的 生产计划、相应的总运量及岩石和矿石产量。綾镝鯛駕櫬鹕踪韦辚糴飙铳。某露天矿有铲位 10 个，卸点 5 个，现有铲车 7 台