发布时间 : 星期二 文章全国大学生数学建模竞赛全国一等奖论文资料设计更新完毕开始阅读b7193f42dc88d0d233d4b14e852458fb760b3899
^`
率，
1di表示平均传染期，建立模型 N??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.5 , 696.5 , 989.2 ,1286.2 ,1612.3 ,1809?，紧接着对X?1??k?作紧邻均值生成得出数据阵B和数据向量Yn，再对参数列??[a,b]T进行最
?小二乘估计最后建立出了灰色模型（GM(1,1)模型）。我们又经过对GM(1,1)模型的残差检验和关联度检验，最终得出了预测结果。
三、符号说明
符号 含 义 日接触率 日治愈率 疾病传播内所考察地区的总人数 整个传染病期内每个病人有效接触的平均人数 易感人群总数 隔离人数比例 未隔离人数比例 接触后没有及时隔离治疗的人数 新增病人数 病人仍患病的概率 单位 人 人 人 人 人 人 人 备注 常量 常量 常量 常量 常量 常量 ? ? N ? M g w p Z ?
^`
四、模型假设
1、假设已确诊人数作为主要的预测模型的指标，对于甲流感病情的预测没有影响。 2、假设所有的统计数据真实，没有遗漏现象。
3、假设与患者有效接触的易感染者（即未患过该病的健康者）均会被传染。
4、假设所考查人群的总数恒定，没有其他病源的输入和输出，不考虑总人口的出生率和自然死亡率。
五、模型的建立与求解
5.1 对问题一建立模型与求解 5.1.1 已确诊病例散点图
根据问题一，由附件1（香港疫情数据）中的已确诊病例数据,用Mtlab7.1作出如下散点图（程序参见附件3）：
图1 散点图
从图1可看出，前25天（即5月20日至6月15日），甲流的传播速度增长幅度较大，而后四十天，甲流的传播速度持续增长，但增长速度趋于平缓。 5.1.2 马尔萨斯模型（Malthusian 模型）
甲流传播预测模型类似于人口增长的预测模型，故首先采用马尔萨斯模型（Malthusian 模型）进行建模。设时刻t的病人人数x?t?是连续、可微函数，并且每天每个病人有效接触（足以使人致病的接触）的人数为常数?，考察t到t+?t病人人数的增加，则有
x?t??t??x?t???x?t??t 再设t?0时有x0个病人，即得微分方程
^`
解之可得：
?dx???x ?dt??x?0??x0 x?t??x0e?t 其中， x0,? 为常数。
根据香港疫情数据中的已确诊的病例数据散点图(图1)，考虑利用马尔萨斯模型
x?t??x0e?t来预测甲流的传播情况。用matlab7.1求得x0?1107.8 ??0.0175。即得
马尔萨斯模型如下（程序参见附件4）： x?t??1107.8e0.0175t 模型Ⅰ
图2 马尔萨斯拟合及预测图形
结果表明，随着t的增加，病人人数x?t?无限增长。即马尔萨斯拟合及预测图线与香港疫情中的已确诊病例数据图线拟合程度较差，且对未来预测情况跟实际显然是不太
相符合的，因此暂不考虑用该模型进行数据预测。
对模型Ⅰ的结果分析：马尔萨斯模型是关于人口或种群增长的模型，它发现人口或种群成指数增长。即在模型I中可引意为，患病人数随着时间得增长呈指数增长变化。但现实生活中，由于病人在有效接触的人群中，包含健康人和病人，而其中只有健康人才可以被传染为病人，因此在改进的模型中必须避免将健康人和病人混为一体这种情况，即要区别病人和健康人进行建模。 5.1.3 阻滞增长模型（Logistic模型）
在疾病传播期内所考察地区的总人数N不变，即不考虑生死，也不考虑迁移。人群分为易感染者和已感染者两类，以下简称健康者和病人。时刻t这两类人在总人数中占得比例分别记作s?t?和i?t?。
^`
假设病人每天的有效接触的平均人数是常数?，?成为日接触率。当病人与健康者有效接触时，使健康者受感染变为病人。
根据假设，每个病人每天可使?s?t?个健康者变为病人，因为病人数为Ni?t?，所以每天共有?Ns?t?i?t?个健康者被感染，于是?Nsi就是病人数Ni的增加率，即有
又因为
s?t??i?t??1
Ndi??Nsi dt再记初始时刻?t?0?病人的比例为i0，则
?di???i?1?i? ?dt??i?0??i0解之得：
i?t??
1?1???t1???i?1??e?0? 模型Ⅱ
用Mtlab7.1作出i?t?~t和
di： ~i的图形如下（程序参见附件5）
dt
图3 Logistic模型i?t?~t曲线 图4 Logistic模型模型Ⅱ结果分析：由图4可知，当i?0.5时
di~i曲线 dt
di?di?达到最大值??，这个时刻为 dt?dt?m
?1??tm???1ln??1?i? ?0?此时病人数增加得最快，预示着传染病的高潮的到来。tm与?成反比，由于日接触