全国大学生数学建模竞赛全国一等奖论文资料设计

1、2012高教社杯全国大学生数学建模竞赛承诺书我们仔细阅读了中国大学生数学建模竞赛的竞赛规则 .我们完全明白，在竞赛开始后参赛队员不能以任何方式（包括电话、电子邮件、网 上咨询等）与队外的任何人（包括指导教师）研究、讨论与赛题有关的问题。我们知道，抄袭别人的成果是违反竞赛规则的，如果引用别人的成果或其他公开的 资料（包括网上查到的资料），必须按照规定的参考文献的表述方式在正文引用处和参 考文献中明确列出。我们郑重承诺，严格遵守竞赛规则，以保证竞赛的公正、公平性。如有违反竞赛规 则的行为，我们将受到严肃处理。我们参赛选择的题号是（从 A/B/C/D中选择一项填写）： 我们的参赛报名号为（如果赛区设

2、置报名号的话）：所属学校（请填写完整的全名）：参赛队员（打印并签名）：1.2.3.指导教师或指导教师组负责人（打印并签名广日期：年月日赛区评阅编号（由赛区组委会评阅前进行编号）：2012高教社杯全国大学生数学建模竞赛编号专用页赛区评阅编号（由赛区组委会评阅前进行编号）：赛区评阅记录（可供赛区评阅时使用）：评阅人评分备注全国统一编号（由赛区组委会送交全国前编号）：全国评阅编号（由全国组委会评阅前进行编号）甲型H1N1流感的预测、控制和影响模型摘要甲型H1N1流感是全国乃至全球人们最受关注的传染病，它的传播速度快，对人们 的身体健康危害极大。本文根据香港甲流疫情数据进行分析，对其传播的预测与控制进

3、 行研究并建出模型，并提出模型建立的关键和困难以及对卫生部门所采取的预防措施作 出评定估计。针对问题一，为了了解甲流的传播情况，先作出已确诊的病例散点图。根据散点图的情况，分别建立了马尔萨斯模型：xt1107.8e0.0175t ,阻滞增长模型:1,SIS模型:1- 1 e ti0di dti i (1 -) , SIR 模型:ds dtdidtSiNi案P s（小）g-sSi ,以及SIR模型的改进模型：工 qidtdti 0ioS 0S0di dt d dt(g ) qiP s()g从SIR模型的改进模型中，可以得出控制传染源、切断传播途径、保护易感人群、隔离 等措施进行预防和控制H1N1

4、甲流的传播。针对问题二，考虑H1N1对旅游经济的影响，对近几年香港接待海外游客的数据进行拟合，得出2009年后三个月的游客数目y1 25.5199,y2 26.7907 , y3 18.1984 ,进dx 1而建立灰色预测模型：dtx k 1-1-0.0124X250.7669,并对其模型(x(0)(1) b)ea b 20468e 0.0124 20468 a a进行了残差检验和关联度检验,从而较为准确的预测出2010的旅客人数为274.9568万【关键词】H1N1流感马尔萨斯模型Logistic 模型SIR模型灰色预测法、问题重述2009年3月底至4月中旬，由墨西哥、美国等地相继发生甲型H

5、1N1流感（A/H1N1influenza）疫情逐步迅速地蔓延到世界各地。甲型 H1N1流感（简称甲流）是一种新型 甲型流感病毒引起的急性呼吸道传染病。去年爆发期间全球数千万人染病，死亡人数超 过16000人。截至去年12月21日，我国内地确诊110590例，死亡442人。由于甲流 的传播速度快，对人们的身体健康危害大，因此得到世界卫生组织的重视和人们广泛的 关注。附件1是香港流感疫情的模拟数据；附件 2是香港接待海外旅游人数的模拟数据。 收集和阅读有关甲流的相关数据及文章，建立数学模型，解决如下问题：问题一：对甲流的传播数学模型进行分析，特别地说明怎样才能建立一个真正能够 预测以及能为预防和

6、控制提供可靠、足够的信息的模型，这样做的困难在哪里？同时， 对于卫生部门所采取的措施做出评论，如：提前或延后5天采取严格的隔离措施，对疫情传播所造成的影响做出估计（附件 1提供的数据可供参考）。问题二：收集甲流对经济某个方面影响的数据，建立相应的数学模型并进行预测（附 件2提供的数据可供参考）。、问题分析根据附件1香港疫情数据分析，我们初步观察到在对 65天甲流传播情况包含了对 已确诊病例、疑似病例、死亡人数累计量以及治愈出院人数累计量。依据这些数据，首 先我们对香港疫情中的已确诊病例情况做出定量分析，运用Mtlab7.1编程得出了甲流传播速度情况的散点图。针对传染病的传播过程，首先，我们用x

7、t表示时刻t的病人人数，用 表示每天每个病人有效接触的人数，考虑t至ht时刻病人人数的增加，建dx_,立微分万程 x , x 0 x0 ,通过马尔萨斯模型求解得：xt xgeto接着在病人的有效接触人群中只有病人方可被传染为病人，因此要区分健康人和病 人。那么我们再次对这些数据进行分析，用常数 表示日接触率；st表示健康者；it表示病人；用Ni t表示病人数。那么由此可知每天共有Nst i t个健康者被感染。建dt立模型Nd1 Nsi, st i t 1 ,通过阻滞增长模型求解得：it 1/ 1 - 1 e t接着我们考虑当治愈后的健康者还可被感染变成病人的情况，我们用 表示日治愈率，1表示平

8、均传染期，建立模型N di-Nsi Ni。dt对于问题二，首先我们利用 2003年至2008年后7至9月份各个月份的平均值与 2009年做差值，利用其差值进行拟合，利用 Mtlab7.1求得2003年至2008年与2009 年后三个月的差值为2.4468,-2.2407 , 0.6516,从而得到2009年后三个月香港海外旅 游人数。接着同样运用 Mtlab7.1编程对2003年到2009年香港海外旅游总人数进行了 处理并假设 X 0 k 229.2,217.3,250,292.7,297,326.1,196.69 ,再对其作一次性累加生成运算得到新的生成数列 X 1 k 229.2 ,446

9、.5,696.5,989.2 ,1286.2,1612.3 ,1809，紧接着对X 1 k作紧邻均值生成得出数据阵B和数据向量Yn,再对参数列a,bT进行最小二乘估计最后建立出了灰色模型 (GM(1,1)模型)。我们又经过对GM(1,1)模型的残差 检验和关联度检验，最终得出了预测结果。三、符号说明符号含 义单位备注日接触率人常量日治愈率人常量N疾病传播内所考察地区的总人数人常量整个传染病期内每个病人有效接触的平均人数人常量M易感人群总数人g隔离人数比例常量w未隔离人数比例常量P接触后没后及时隔离治疗的人数人Z新增病人数人病人仍患病的概率四、模型假设1、假设已确诊人数作为主要的预测模型的指标，

10、对于甲流感病情的预测没有影响。2、假设所有的统计数据真实，没有遗漏现象。3、假设与患者有效接触的易感染者（即未患过该病的健康者）均会被传染。4、假设所考查人群的总数恒定，没有其他病源的输入和输出，不考虑总人口的出生率 和自然死亡率。五、模型的建立与求解5.1 对问题一建立模型与求解5.1.1 已确诊病例散点图根据问题一，由附件1 （香港疫情数据）中的已确诊病例数据，用Mtlab7.1作出如 下散点图（程序参见附件3）:图1散点图从图1可看出，前25天（即5月20日至6月15日）,甲流的传播速度增长幅度较 大，而后四十天，甲流的传播速度持续增长，但增长速度趋于平缓。5.1.2 马尔萨斯模型（Ma

11、lthusian 模型）甲流传播预测 模型类 似于人口增长 的预测 模型，故首先 采用马 尔萨斯模型（Malthusian模型）进行建模。设时刻t的病人人数xt是连续、可微函数，并且每天每个病人有效接触（足以使人致病的接触）的人数为常数，考察t到t+ t病人人数的增加，则有再设t 0时有x0个病人，即得微分方程dx x dt x 0 X0解之可得：x t x0e t其中，xo，为常数。根据香港疫情数据中的已确诊的病例数据散点图（图1）,考虑利用马尔萨斯模型xtx0e t来预测甲流的传播情况。用 matlab7.1求得x0 1107.80.0175。即得马尔萨斯模型如下（程序参见附件 4）: 0

12、.0175tx t 1107.8e模型 ifile Edit Ifi av ItJE4rl Tods Qasktoi口齿。昌玲，鼻踮 Me 丹口图2马尔萨斯拟合及预测图形结果表明，随着t的增加，病人人数xt无限增长。即马尔萨斯拟合及预测图线与香港疫情中的已确诊病例数据图线拟合程度较差，且对未来预测情况跟实际显然是不太相符合的，因此暂不考虑用该模型进行数据预测。对模型I的结果分析：马尔萨斯模型是关于人口或种群增长的模型，它发现人口或 种群成指数增长。即在模型I中可引意为，患病人数随着时间得增长呈指数增长变化。 但现实生活中，由于病人在有效接触的人群中，包含健康人和病人，而其中只有健康人 才可以被

13、传染为病人，因此在改进的模型中必须避免将健康人和病人混为一体这种情 况，即要区别病人和健康人进行建模。5.1.3 阻滞增长模型（Logistic 模型）在疾病传播期内所考察地区的总人数 N不变，即不考虑生死，也不考虑迁移。人群 分为易感染者和已感染者两类，以下简称健康者和病人。时刻t这两类人在总人数中占得比例分别记作st和i t。假设病人每天的有效接触的平均人数是常数，成为日接触率。当病人与健康者有效接触时，使健康者受感染变为病人。根据假设，每个病人每天可使st个健康者变为病人，因为病人数为 Ni t ,所以每天共有Nst i t个健康者被感染，于是 Nsi就是病人数Ni的增加率，即有N di

14、 Nsi dt 又因为s t i t 1再记初始时刻t 0病人的比例为i0 ,则di dti 0i 1 iio用Mtlab7.1作出i tt和包i的图形如下 dt图3 Logistic 模型i t t曲线（程序参见附件5）:模型II结果分析：由图4可知，当i0.5时更达到最大值 dt生,这个时刻为出m解之得:11- 1 eio图4 Logistic 模型电i曲线dt,11dtm ln 1io此时病人数增加得最快，预示着传染病的高潮的到来。tm与 成反比，由于日接触率 反应了该地区的卫生水平，越小卫生水平越高。所以改善保健设施、提高卫生水平可以延缓传染病高潮的到来。而当t时i 1,即所有人终究将

15、被传染，全变为病人，这显然与实际情况不符相。其中的原因是模型中没有考虑到病人是可以治愈的，人 群中的健康者只能变成病人，而病人不会再变成健康者。下面模型中将讨论病人可以治 愈的情况。5.1.4 SIS 模型由于病人被治愈后变成健康者，健康者还可以被感染再变成病人，那么由此得到需增加的条件为：每天被治愈的病人数占病人总数的比例为常数，称为日治愈率。病人治愈后成为仍可被感染的健康者。显然-是这种传染病的平均传染期。N 电 Nsi Ni dts i(t) 1记初始时刻t 0病人的比例为八，则dt i(1 i)ii(0) i0设 一，则 可表示整个传染病期内每个病人有效接触的平均人数，称为接触数。利

16、用，可得如下模型：史 i i (1 -)模型田dt根据模型田，利用Mtlab7.1作出一 i的图形，如下(程序参见附件6)： dt图5 SIS模型的曳一曲图6 SIS模型的it曲dt1是一个阈值。由图5可知道，随着病人模型田结果分析：不难看出，接触数所占的人数越多，那么在时间t内病人的增长率就越大。当 1时it的增减性取决于io1的大小（见图6）,单具极限值i 1 ,随着 的增加而增加；当1时病人比例it越来越小，最终趋于0,这是由于传染期内经有接触从而使健康者变成的病人数不超过 原来病人数的缘故。5.1.5 SIR 模型由于病人在治愈后有一定的免疫力，所以病愈的人既非健康者（易感染者），也非

17、病人（已感染者），他们己经退出传染系统。人群分为健康者、病人、病愈与免疫的移出者三类，即SIR模型。这三类人在总人数 N中占得比例分别记作s（t）,i和r（t）。s（t） i（t） r（t） 1NdNidt记初始时刻的健康者和病人的比例分别是So So 0和io io 0 （设移出者的初始值di dt dsdt i o soro 。），则可得SIR模型:SiNi模型IVSiioSo由于模型IV无法直接求出S和i（t）的值，故先作数值运算。设 2,o.5,i（o） o.o2,S（o） o.96,用Mtlab7.1求解可得如下i（o）,S（t）图形和s i图形（程序参见附件7）:图 7 i(o),

18、s(t)图形图8 si图形（相轨图）模型IV结果分析：is平面称为相平面，相轨线在相平面上的定义域为s,is 0,i Qs消去dt可得:dids1,s S0i0利用积分特性可解得：isoio1 s 1.-ln,在定义域D内，该式表小的曲线即为相so轨线，如图9所示.其中箭头表示了随着时间t的增加s和i(t)的变化趋向.图9 SIR模型的相轨线根据图9,可分析s(t),i(t)和r(t)的变化情况如下:1 .不论初始条件s0,i。如何，病人消失将消失，即：i 02 .最终未被感染的健康者的比例是s .3 .若s0 1/ ，则i(t)先增加，后减小.4 .若s0 1/ ，则i(t)单调减小至0.5

19、.1.6 模型IV (SIR模型)的改进模型由于在H1N1流行的过程中，各个地方(包括香港)都采取了一定的措施，一般是 采取了隔离的制度，所以在模型模型IV ( SIR模型)的基础上进行改进。考虑到隔离人 数比例g和未隔离人人数比例w ,以及接触后没有及时隔离治疗的人数 p,从而建立如 下改进模型:dsdtdg gtdrdtqi模型Vdidt (g ) qid dtP s(由于该模型的分析过程过于复杂，所以该模型在这里将不多做讨论。但从该模型中, 可以看出为预防和控制提供可靠的信息，比如：控制传染源、切断传播途径、保护易感 人群、隔离等。5.1.7 建立模型的关键和困难建立模型的关键在于对模型

20、进行动态的分析，当传染病发展到一定阶段时，在医疗 水平提高、人员流动、出生率和死亡率以及以及采取防御传播措施等方面的影响促使传 染率下降。此时仍用之前模型的误差会很大。在建立模型过程中有以下几个方面的困难： 对不同地区H1NI的卫生知识的宣传程度，K值取值不同；对某一地区的不同地方的 强化管理也不一样，K值也就不一样；防护措施不同、卫生条件不一等，都会影响到 K的取值。另外，本文模型大多假设种群总数为常数，且考虑的影响因素较少，但在实 际问题中，由于疾病的复杂性往往涉及变动人口、年龄结构、隔离等多种因素的影响， 致使模型的建立错综复杂。5.1.8 对卫生部门采取的措施评价经上网查询得知医学研究

21、表明，从正式发病到治愈一般需一至两周，假定平均治愈 时间为10天。假设新患者出现的数量与现有患者的数量成正比，也与现有易感者的数 量成正比，即发病率是患者人数和易感者人数的双线性函数。则有：M (t 1) M(t) Z(t 1) tZ(t 1) M(t) Z(t ) ( ) t 0对其进行整理可得:M (tM(t)1)t(1 Z(t ) ( )模型VI其中，M为t时刻易感人群总数，Z（t）为t时刻新增病人数，（）为病人从患病起经过 时间仍为病人的概率（图中用p表示）。假设病人开始,患病记为第1天，最迟到第10天病愈。那么病人从患病起经过时间仍为病人呈逐步递减的概率参数，如下图：图10 概率一时

22、间图由图10可看出，如果病人发病后5天才开始隔离的话，病人仍患病的概率相当大(图10阴影区域D),即病人在社会上与易感人群的接触率也相对较大。由模型V可得：M(t 1)(1' Z(t ) ( ) 1M即M(t 1) M(t),说明易感人群总数将会以较大的数值递减，给疫情的控制带来更大的困难。所以，如果在病人发病前提前 5天隔离的话，新增病人数将变得很小。5.2 问题二的模型的建立与求解问题二要我们收集甲流对经济某个方面影响的数据并建立相应的数学模型并进行 预测，针对该问题二，我们充分利用附件二，建立甲流对旅游带来的经济影响，而旅游 经济与游客数目成正比例关系，故建立预测游客数目模型来预

23、测旅游经济。5.2.1 香港接待海外旅游人数折线图根据附件2,利用Excel2003作出2003年至2009年各个月份香港接待海外旅游人 数的折线图，如下：碧港揍得海外旅游人数(屋位万人) 国关区I图11香港接待海外旅游人数折线图从图11可看出，2003至2008年各整年的海外游客人数的增长率相对稳定，2009年前三月份海外游客人数稳定，从四月份至六月份是因受到H1N1影响而急促下降，从七月份至九月份海外游客又逐步的上升，十月份至十二月份就是要预测的。5.2.2 2009 年后三个月预测为了预测2009年后3个月的海外旅游人数，根据图11折线变化，利用2003年至 2008年后7至9月份各个月

24、份的平均值与 2009做差值，利用其差值进行拟合，利用 Mtlab7.1求得2003年至2008年与2009年后三个月的差值为 2.4468 ,-2.2407 ,0.6516.i7891011126(ai)/6i 122.466727.466727.150027.966724.550018.8500Y (2009年)2.618.816.2y1y2y36(ai)/6 Yi 119.856718.666710.950027.9667-y124.5500-y218.8500-y3即 27.9667-y1 2.4468, 24.5500-y2 2.2407, 18.8500- y3 0.6516.所以

25、 2009 年后三个月香港接待海外旅游人数分别为：y1 25.5199 ,y2 26.7907 , y3 18.1984 (单位：万人).5.2.3 灰色预测模型为了预测2010年香港接待海外旅游总人数，先分别计算出 2003至2009年每年的 总人数，得出如下表(单位：万人)：年份2003200420052006200720082009旅客人数229.2217.3250292.7297326.1196.69假设设 X0 k 229.2,217.3,250,292. 7,297,326.1,196.695.2.3.1 GM(1,1) 模型的建立为了使其成为有规律的时间序列数据，对其作一次累加生

26、成运算，即令tX(1)(t) X(0)(k) (t,k 1,2 7) n 1从而得到新 的生成数列 X1 k 229.2 ,446.5,696.5,989.2,1286.2 ,1612.3 ,1809 .对X1 k做紧邻均值生成.则数据阵B和数据向量Yn为121 x21 x31-337.8510 x02217.3113141-571.501x3250xx-842.8510292.72Ynx04114151-1137.71x5297xx2-1449.31x06326.111 x51 x61-1710.71x07196.691212211 x1 x1 x1对参数列a, bT进行最小乘估计，可得BT

27、B1BT0.00050.0003 0.0001-0.0001 -0.0003 -0.00050.65930.4876 0.28820.0716 -0.1574 -0.3494T 1 tBTBBTYn-0.0124250.7669(其中，a为发展系数，反映X的发展趋势;b为灰色作用量，反映数据间的变化关系.)从而可得出GM(1,1)模型：dx 1一 0.0124x1250.7669出模型即(0)0 、b、a b0.0124xk 1 (x )e 20468e20468a a其中，x k 1 (x(0)(1) b)ea b 20468e 0.0124 20468 为时间响应函数形式。 a a5.2.

28、3.2 GM(1,1)模型的残差检验残差大小检验，即对模型值和实际值的残差进行逐点检验.1(1)根据预测公式，计算X k ,得1X k 229.2,446.5,696.5,989.2,1286.2,1612.3,1809 k 0,1, 6 0(2)累减生成X k序列，k 1,2 7 0X k 229.2 ,255.1972 ,258.389,261.6207,264.8928,268.2059,271.5603原始序列：X 0 k 229.2,217.3,250,292.7,297,326.1,196.69(3)计算绝对残差和相对残差序列绝对残差序列：0,0.1744 ,0.0336,0.10

29、62,0.1081,0.1775,0.3807相对残差：0,0.080258%,0.013440%,0.036283%,0.036397%,0.054431%,0.19%GM(1,1)模型的残差检验结果：相对残差不超过0.19%,精确度高。5.2.3.3GM(1,1)模型关联度检验关联度是用来定量描述各变化过程之间的差别.关联系数越大，说明预测值和实际 值越接近.0(1)计算序列X0与X 的绝对残差序列 0 k00,0.1744 ,0.0336,0. 1062,0.108 1,0.1775,0 .3807min 0 k min 0,0.1744 ,0.0336,0.1062,0.1081,0.

30、1775,0.38070max 0 k max 0,0.1744 ,0.0336,0.1062,0.1081,0.1775,0.38070.3807(2)计算关联度(k)min 0 (k) P max 0 (k)0 (k) P max 0 (k)(k 1,.,7,P 0.5)1 nrii(k) 0.9547n k 1精度检验等级如下表:好（1级）0.90合格（2级）0.80勉强（3级）0.70不合格（4级）0.70年份20032004200520062007200820092010实际旅客229.2217.3250292.7297326.1196.69预测旅客229.2255.197258.3

31、89261.621264.893268.206271.560274.957GM（1,1）模型关联度检验2果：关联度为0.9547 0.90,精确度高。5.2.3.4GM（1,1）模型求解利用Mtlab7.1进行预测（程序参见附件8）,得到实际值与预测值如下表（单位：万人）：即2010年香港接待海外旅游总人数为274.957万人六、模型评价与推广模型评价：在建模前期，全面分析影响甲流疫情的各种因素，找出各因素之间的关系以及作用 的时间段和范围，收集比较完整而准确的前期数据。在模型建立中我们采用了各种软件 （如Mtlab, Excel等）进行求解，制图精确，计算结果较为准确。但在预测模型中，时 间

32、序列数据的时间间隔不是稳定的，这对模型的求解结果的准确性有一定的影响。本文 所建立的控制模型忽略了人口流动、变化给该地区甲型H1N1流感带来的影响，从而模型预测结果会与实际情况有一定差距。模型推广：通过模型的分析可知，如果全社会的努力和投入的程度继续增加，即隔离措施的提 早进行、隔离率增大、防疫药品的早日研发、公众的防御意识提高，可使得疫情周期缩 短、患者人数逐步减少。实时监控甲流疫情走势，采集更多的数据以验证模型和改进模 型，若有预料之外的干扰因素出现，应及时修正模型，重新预测其后期走势。七、参考文献1姜启源谢金星叶俊，数学模型M,北京：高等教育出版社，2003 （第三版）.2刘国卫MATL

33、AE®序设计与应用M,北京：高等教育出版社,2011.3 SARSt播的数学模型及应用， 151K2008-4-14.4马知恩周义仓王稳地等.传染病动力学的数学建模与研究M.北京，科学出版 社，2004.2.八、附录附件1:香港疫情的数据日期已确诊病例累计现有疑似病例死计治愈出院累计5月20日33940218335月21日48261025435月22日58866628465月23日69378235555月24日77486339645月25日87795442735月26日988109348765月27日1114125556785月28日1199127559785月29日13471358

34、66835月30日1440140875905月31日15531415821006月01日16361468911096月02日17411493961156月03日180315371001186月04日189715101031216月05日196015231071346月06日204915141101416月07日213614861121526月08日217714251141686月09日222713971161756月10日226514111201866月11日230413781292086月12日234713381342446月13日237013081392526月14日23881317140

35、2576月15日240512651412736月16日242012501453076月17日243412501473326月18日243712491503496月19日244412251543956月20日244412211564476月21日245612051585286月22日246511791605826月23日249011341636676月24日249911051677046月25日250410691687476月26日251210051728286月27日25149411758666月28日25178031769286月29日252076017710066月30日252174718

36、110877月01日252273918111247月02日252273418111577月03日252272418111897月04日252271818112637月05日252271618113217月06日252271318314037月07日252366818314467月08日252255018415437月09日252245118416537月10日252235118617477月11日252325718618217月12日252315518718767月13日25227118719447月14日2522418919947月15日2522318920157月16日2521319020

37、537月17日2521519021207月18日2521419121547月19日2521319121717月20日2521319121897月21日2521219122317月22日2521219122577月23日252121912277附件2:香港接待海外旅游人数（单位：万人）年1月2月3月4月5月6月7月8月9月10月11月12月20039.4 11.3 16.8 19.8 20.3 18.8 20.9 24.9 24.7 24.3 19.4200418.620059.6 11.7 15.8 19.9 19.5 17.8 17.8 23.3 21.4 24.5 20.1200615.9

38、200710.1 12.9 17.7 21.0 21.0 20.4 21.9 25.8 29.3 29.8 23.6200816.5200911.4 26.0 19.6 25.9 27.6 24.3 23.0 27.8 27.3 28.5 32.818.511.5 26.4 20.4 26.1 28.9 28.0 25.2 30.8 28.7 28.1 22.220.713.7 29.7 23.1 28.9 29.0 27.4 26.0 32.2 31.4 32.6 29.222.915.4 17.1 23.5 29.2 11.6 1.78 2.61 8.8 16.2注：以上数据为模拟数据。部

39、分参考文献：100432420 515523.html附件3:x1=1:65;y1=339 482 588 693 774 877 988 1114 1199 1347 1440 1553 1636 1741 1803 1897 1960 2049 2136 2177 2227 2265 2304 2347 2370 2388 2405 2420 2434 2437 2444 2444 2456 2465 2490 2499 2504 2512 2514 2517 2520 2521 2522 2522 2522 2522 2522 2522 2523 2522 2522 2522 2523

40、2523 2522 2522 2522 2521 2521 2521 2521 2521 2521 2521 2521;plot(x1,y1, '-*')附件4:x1=1:65;y1=339 482 588 693 774 877 988 1114 1199 1347 1440 1553 1636 1741 1803 1897 1960 2049 2136 2177 2227 2265 2304 2347 2370 2388 2405 2420 2434 2437 2444 2444 2456 2465 2490 2499 2504 2512 2514 2517 2520 25

41、21 2522 2522 2522 2522 2522 2522 2523 2522 2522 2522 2523 2523 2522 2522 2522 2521 2521 2521 2521 2521 2521 2521 2521;log_y1=log(y1);p=polyfit(x1,log_y1,1)x= exp(7.0101)x1=1:65;y1=339 482 588 693 774 877 988 1114 1199 1347 1440 1553 1636 1741 1803 1897 1960 2049 2136 2177 2227 2265 2304 2347 2370 23

42、88 2405 2420 2434 2437 2444 2444 2456 2465 2490 2499 2504 2512 2514 2517 2520 2521 2522 2522 2522 2522 2522 2522 2523 2522 2522 2522 2523 2523 2522 2522 2522 2521 2521 2521 2521 2521 2521 2521 2521;y2=1107.8*exp(0.0175*x1);plot(x1,y1, ':o' ,x1,y2, '-*');附件5:t=0:0.01:13;i=1./(1+(1./0.

43、1-1).*exp(-1.*t);plot(t,i)title ('SI模型的it曲线');xlabel( 't');ylabel( 'i');axis(0 13 0 1.1);x=0:0.01:1;y=x-x.*x;plot(x,y)title ('SI模型的di/dti曲线');xlabel( 'i');ylabel( 'di/dt');axis(0 1 0 0.3);附件6:x=0:0.01:1;y=0.7.*x-x.A2;piot(x,y);title ('SIS 模型的 di/dti 曲线')；xlabel( 'i'); ylabel( 'di/dt');function y=aini(t,x)a=2;b=0.5;y=a*x(1)*x(2