2015年南京理工大学研究生数模大赛A题论文解析.doc

小反刍兽疫疫情研究 摘要 本文通过对小反刍兽疫疫情统计数据的分析，建立了经典的SIR模型，并通过运用数学软件mathematica9.0对该模型进行数值求解，并将数值结果与实际进行对比分析。通过对经典的SIR模型进行修改，得到较为合理的捕杀模型，并从验证捕杀模型经典的SIR模型相比，捕杀模型能更好的描述小反刍兽疫疫情传播的规律，小反刍兽疫疫情的预防与控制。关键词：小反刍，SIR模型，捕杀模型，mathematica 一 问题重述 小反刍兽疫（peste des petits ruminants，PPR）俗称羊瘟，又名小反刍兽假性牛瘟（pseudorinderpest）、肺肠炎（pneumoenteritis）、口炎肺肠炎复合症（stomatitis-pneumoenteritis complex），是由小反刍兽疫病毒引起的一种急性病毒性传染病，主要感染小反刍动物，以发热、口炎、腹泻、肺炎为特征。我国部分地区的羊发生小反刍兽疫疫情，对当地的畜牧业造成了严重的损失。所以要切实做好小反刍兽疫的综合防控工作，严防本病传播，此病引起场区广大养殖户高度重视。小反刍兽疫i(t)和r(t)。假设疫区羊群总数为常数N，N足够大，可以把变量s(t)、i(t)、r(t)视为连续变量，并进一步假定为连续可微变量； 3.每只感染的羊单位时间有效接触的平均羊的数量是λ，λ称为日接触率且始终为一定值。当感染的羊与健康羊有效接触时，使健康羊受感染变为病羊； 4.捕杀前，疫点地区只有感染羊群； 5.除少数县外，大部分地区羊群的数量基本相等，超过平均数量的发病数据不具有代表性； 6.捕杀的羊群包括感染却未死亡的羊群和与感染羊群接触的健康羊群，忽略健康羊群。 SIR模型示意图如图1所示。 图1 SIR类模型示意图 2.2模型的符号定义 :t时刻易感类羊群占疫区羊群总数的比例； :t时刻感染羊群占疫区羊群的总数的比例； :t时刻移除羊群占疫区羊群总数的比例； ：单位时间内每只感染的羊有效接触的平均羊的数量； :单位时间内感染羊群因病死亡的数量与t时刻感染羊群数量的比值。 2.3机理分析 为了简化模型，这里不考虑疫区羊群在疫情持续期间出生率和死亡率带来的影响，以全国羊群为研究对象。羊群可以分为3类： 类：称为未感染羊群，该类羊群没有感染，但对疾病没有免疫能力，与患病羊群接触，容易感染该疾病； 类：称为已感染羊群，该类羊群已感染，而且可以传染给类羊群； 类：称为移除类羊群，包括类被捕杀和类因捕杀死亡的羊群。 、、分别表示t时刻上述3类羊群占全国羊群总数的比例。 羊群中3类羊群均匀分布，传播方式为接触性传播。单位时间内一个感染羊与其它羊的接触率为，则一个感染羊在单位时间内与类羊群的接触率为，所以单位时间内S类羊被感染的数量为。单位时间内感染羊群死亡的个数为。于是时单位时间内类羊群增加的数量。这就是时单位时间内类羊群增加的数量，称为发病率，它是和的双线性函数。 单位时间内共有只羊被感染为类，于是就是单位时间内类羊群的增加率。 单位时间内移除的羊群数量包括，死亡与捕杀的羊群。单位时间内死亡的羊群数量为。单位时间内类羊群的增加率为。 感染羊群的死亡可以分为两类：(1)因病自然死亡；(2)感染后虽未因病死亡，却被捕杀。因此，疫情初始阶段，未进行捕杀时，死亡率为()；进入捕杀阶段以后，死亡率。 2.4SIR模型建立 基于以上机理分析，针对小反刍兽疫 （1） （2） （3） （4） 2.5相轨线分析 s~i平面称为相平面，相轨线在相平面上的定义域为 在方程（1）（2）消去，可得 （5） 令，得出 时，，i随s的增加单调减小； 时，，i随s的增加单调增加； 时，，i在处取得极大值。 容易求出方程（5）的解为 （6） 在定义域D内，(6)式表示的曲线即为相轨线，如图2所示，其中箭头表示了随时间t的增加和的变化趋势。 下面根据（1），（2），（6）式和图2分析，和的变化情况（时它们的极限值分别记作，和）。 1.不论初始条件，如何，感染羊群终将消失，即 证明如下： 首先，由（1）式得，而，故存在；由（3）式得，而，故存在，在由（4）式知存在。 其次，若，则由（3）式，对于充分大的t有，这将导致，与存在相矛盾。 从图形上看，不论相轨线从或从点出发，它终将与s轴相交（t充分大）。 2.最终未被感染的羊群的比例是，在（6）式中令得到，是方程