微分方程模型的甲型H1N1病毒传播分析与预测.docVIP

基于延迟微分方程模型的甲型H1N1病毒传播分析与预测 摘要 对于问题（1），本文以经典的SIR模型为基础建立了含有延迟参数的四房室微分方程模型，并对其进行了简化，利用线性拟合的方式确定了传染率的变化函数，然后将它代入微分方程，最后求出香港的确诊病人模拟变化情况，并与实际情况比，有较好的模拟效果，发现现在香港的发病情况正处于一种比较缓和的状态。接着，又建立一个灰色预测模型，对香港地区疫情情况进行拟合和预测，结果也能与实际情况很好的吻合。最后对、澳门、台湾等地疫情传播规律进行了分析，依此做出短期和长期预测。 然而，问题（2）中，由于省份与地区的不同，引入输入因子将模型进行了改进，通过对现有广东数据的处理计算，建立新的的含延迟参数的微分方程模型，同样利用线性拟合的方式确定传染率的变化函数，最终将广东的实际数据于拟合数据比较，效果较为满意。另外，还对北京、上海、浙江等地的传染病情况作出相应的分析与预测。 最后，对于问题（3），本文进行了敏感性分析，定量的说明了政府控制力度，自由带菌者进入速率，传染率的减慢速率等参数对疫情的扩散速度及最终规模产生的影响，并依据此分析给出了防止病毒大面积扩散的方案和措施。 关键词： 房室模型 微分方程 线性拟合 敏感性分析 灰色预测 一，问题重述： 甲型H1N1流感于2009年爆发，是猪的一种急性、传染性呼吸器官疾病，病毒侵入人体后使人致病。2009年5月初，甲型H1N1流感病毒感染给人导致人类致病。从墨西哥和美国开始，在人与人之间传播，极短时间内迅速扩散。世界卫生组织WHO告知，甲型H1N1病毒传染人后有1—7天的潜伏期，发病早期症状与普通人流感相似。部分患者病情可迅速进展，急剧恶化，最终导致死亡。甲型H1N1流感与普通流感和2003年的SARS有类似之处，也有不同之处，根据所给数据资料及相关信息以香港为例，建立特定区域的甲型H1N1病毒传播模型，定量分析特定区域内该传染病的传播规律，提出预测和控制方案，并用所建立的模型检验中国内地、中国台湾、中国澳门等三区域的发病情况；香港和内地在很大程度上有较大的区别，根据提高的各个省份信息，内陆可以省为单位进行细化。以广东为例，建立中国内地省份的甲型H1N1病毒传播模型，定量分析该传染病的传播规律，提出预测和控制方案，并用所建立的模型检验北京、上海、浙江、山东、江西等五省市的发病情况；利用所建立的模型对我国内地和港澳台地区未来甲型H1N1流感的发展趋势进行一周短期预测和一月长期预测，为疾病的预防和控制提出建议。 二，问题分析 本题是关于传染病的扩散与防治问题，需要通过数学模型来分析甲型H1N1传染病传播的特定规律，并且以此为基础作短期和长期预测。 分析一般传染病的传播规律，人群中的患者作为传染源可以使健康人被感染，直接转化为患者或是病毒携带者，接着在患者，携带者，易感人群中，上述过程将不断进行下去。由此可建立经典的传染病SIR模型[1]，其将人群分为三类，易感染者，患病者，移出者。然后通过三者之间的关系来确定传染病的扩散规律。 本题虽可利用经典SIR模型的思想，但许多地方与它有很大不同。首先，影响H1N1扩散的因素有很多，如政府的防控能力、医院的治愈能力、人群的防范意识、人体自身的自愈能力等因素都对病毒的蔓延趋势具有有效的限制能力。其次，SIR模型假设没有病源的输入和输出，这显然与中国的情况是不相符的，因为早期大部分病例是从国外输入的。 考虑到上述这些因素，我们决定改进 SIR模型，将人群做更细致的划分，考虑更多的因素，以做出更准确的分析与预测。 三，模型假设： 1时间起点为5月13日，记此时t=0. 2不考虑出生和死亡的影响。 3. 将H1N1 的传播途径视为与自由带菌者在传播范围内的有效接触； 4确诊病人被完全隔离，不再具有传染性 5.患者治愈后具有对病毒的免疫力 6.自由带菌者的活动是随机的 四，符号说明 符号 说明 单位 当前的天数，以5月13日为基准 天 时刻健康者人数 人 t时刻自由带菌者人数 人 t时刻累计确诊病人人数（含治愈者） 人 t时刻治愈患者人数 人 传染率，每个病人每天能够传染的有效人数 1/天 每日确诊率，每天每个病人转化为确诊的概率 无 每日治愈率 无 每日死亡率 无 五，模型建立与求解 5.1.针对问题（1）建立模型（1） 5.1.1 模型预处理 我们争对香港等特定区域建立模型，由于可利用的数据有限，我们决定建立四房室模型解决此问题： 将人群划分分为四类： 健康者(易受感染者)----用S表示健康者人数 自由带菌者(已受感染者) ----用I表示未被确诊具有传染性的病人的人数。 确诊者（已被隔离者）----用P表示，已经被隔离，退出传染系统，不会感染别人 治愈及死亡者(包括“被治愈者”和“死亡者”)----用R表示其人数，不再参与病