1、 1 毕业设计文献综述 信息与计算科学 常微分方程在数学建模中的应用 人们将数学方法应用到有关传染病方面的研究可追溯到 1760 年 , Bernoulli 在其论文中用数学模型评价天花对期望寿命的影响 . 上世纪初 , Kermark 和 Mckendrick 首先利用动力学方法建立了传染病的数学模型 . 1928 年 Reed 及 Frost 共同提出 Reed-Frost 模型 , 它的基本公式是确定性的 , 之后又获得了该模型的随机过程 . 确定性 Reed-Frost 模型 : 1 (1 )t t tC S qC . 指下一代 将发生的病例数为 t 代的易感者人群与当时有效接触率 1
2、 tqC 的乘积 . 其中 1tC 是在第 1t 代的病例者 , tS 是在第 t 代的易感者 , q 为单位时间内未发生有效接触的概率 . Reed-Frost 模型的假设条件有 5 个 . 条件一 : 研究的群体与其他群体隔绝 ; 条件二 : 在疾病流行期间 , 人群中任何个体间互相接触的机会均等 ; 条件三 : 易感者与感染者充分接触后 , 按一定概率变成新的感染者 ; 条件四 : 感染者在传染期间具有传染性 , 其后成为完全的免疫者 ; 条件五 : 以上条件在流行期间不变 . 该模型的缺陷是 : 研究结果常常与实际情况有一定程度差距 , 这是因为在模型中假设有效接触率传染力是不变的 .
3、 在国内 , 该模型的改进模式 , 认为易感者在感染过程中不仅会变成显性感染者（患者） , 也 可能会变成隐性感染者 , 而且患者和隐性感染者均具有传染能力 . 故对 Reed-Frost 的假设条件三和条件四做如下假设 . 条件三 : 易感者和感染者充分接触后 , 按一定的概率变成新感染者 , 新感染者中患者和隐性感染者成一定比例 ; 条件四 : 患者和隐性感染者在传染期间均具有传染力 , 传染期后均变成完全的免疫者 . 设单位时间内一名患者与易感者的有效接触率为 p , 一名隐性感染者与易感者的有效接触率为 1p , 则该期间 避免一名患者或一名隐性感染者的感染率为 1qp , 111qp
4、 . 设 b 为感染过程中隐性感染者与患者的比例常数 , 则单位时间内新患者数为 tC , 新隐性感染者数为 tbC , 易感者有 tC 名患者和 tbC 名 . 隐性感染者方面受到感染的概率为11 ttqCqbC , 该期间被感染者在下一期间均成为传染源 . 因此 , 1t 时点的新感染者的期2 待值 (E 表示期待 )为 1 1 1(1 )t t t t tb E C E C S q C q b C , 111 (1 )1t t t tE C S q C q b Cb ; 11( 1)t t tS S b C . 以 R 为非显性感染者 , 在 Reed-Frost 模型中 R 与易感者
5、S 完全相同 . 而在改进模型中 , R 包括易感者和隐性易感者 , 即 t t tS R b tC . 则 1t 时点的免疫者I 数为 1 ( 1)t t tI I b C . 随机性 Reed-Frost模型表示当时的易感者全部未受到感染至全部受到感染的各种情况概率 , 每种情况发生的概率可有不同 , r 是下一代病例数 . 该模型体现了传染力的可变性 , 适合广泛的情况 , 是模型发展的新趋向 . 而传染病动力学的常微分方程 SIR 模型便是由 May4 等在 1979 年提出的 , 该模型考虑了 3 类个体 : 正常可被感染者 ; 患病者 ; 已恢复且具有 免疫力者 . 该模型的假设条
6、件分别为 :条件一 : 人群分为易感者、患者和病愈免疫移出者 . 条件二 : 个体获得免疫是永久的 , 这意味着假若某个个体获得免疫 , 他们将永远不会再感染 . 这种模型适合于滤过性霉菌引起的流行病 , 如麻疹、天花、腮腺炎等 . 条件三 : 易感人群的减少速度与易感人群和被感染者数量的乘积呈正比 . 条件四 : 恢复者的增长速度与被感染者的数量成正比 . 该模型对乙型肝炎病毒在人群中的感染和传播有较好的应用 . 通过该模型的参数可以描述疫苗接种前人群 HBV 的动态传播过程 , 也可以预测不同疫苗接种覆盖率时免 疫后人群 HBV 的变化趋势 , 从而评价疫苗的远期效果 . 后来在 SIR
7、模型 10 考虑 3 类个体的基础上 , 增加了 1 类个体 : 已感染但处于潜伏期未发病者 . 上述 4 类个体及描述其相互关系的常微分方程组构成新的传染病动力学模型 : SEIR 模型 . 很多学者对这类模型进行了深入研究 , Michael8 研究了 SEIR 模型的全局动力学性质 . Langlais 9 等对治愈后不具有免疫能力的 SEIRS 模型进行的研究中 , 给出了最大静止状态的稳定性判据 , 他们还将该模型应用于动物传染病的研究 . 近几年 , 人们用数学方法来研究传染病的发病机理、动态过程和发展趋势 , 已逐步成为一个活跃的研究领域 . 在国外 , 数学预测模型已经能够成功
8、地应用于生物分子水平 , 模拟体内病毒的复制及半衰期 , 让我们更加全面地认识并了解了传染病的感染机制 . 而我们的 国内学者吴开琛等也成功的把该模型应用于非典型肺炎 (SARS)的研究 , 并在此基础上提出 5分室模型 , 即 : SEIDR, 其中的 D(death)为人群中感染发病者不治死亡的 . 本文是利用 SIR模型来研究传染病问题的 , 由于传染病流行过程的研究与其他学科有所不同 , 不能通过在人群中实验的方式来获得数据 , 所以有关传染病的数据、资料只能从已有3 的传染病流行的报告中获取 , 这些数据往往不够全面 , 难以根据这些数据来准确地确定参数 , 只能大概估计其范围 .
9、这次论文 主要是通过全面调查、收集相关的数据资料 , 有效应用常微分方程和数学建模的相关知识 , 并充分利用图书馆和互联网上的 丰富的资源来建立 SIR 模型 , 在 对建立好的数学模型进行定量和定性的分析 与探究的过程中 , 观察和研究实际对象的固有特征和内在规律 , 抓住问题的主要矛盾 , 对当今社会中经常爆发的传染病建立常微分方程模型 并利用常微分方程和数学建模的相关知识对它分别进行分析和研究 , 探讨了它的传播规律以及影响它们流行的因素、预测可能发生的后果及如何抑制其流行或恶化 . 这个模型的建立及探究说明了在反映客观现实世界运动过程的量与量之间的关系中 , 大量存在了满足常微分方程关
10、系式的模型 , 需要我们通过求解常微分方程来了解未知函数的性质 , 常微分方程是解决实际问题的重要工具 .所建立的模型 , 在常微分方程的观点剖析下 , 充分展现现代社会生活中常微分方程应用 . 4 参考文献 1 May RM et al, Nature J. Nature Publishing Group, 1979, 180: 455461. 2 Langlais M et al, Math Comp Model J. Elsevier Science, 2000, 31: 117124. 3 陈文江 , 吴开琛等 . 运用数学模型探讨 SARS 聚集性传播的机制 J. 中国热带医学 ,
11、2004, 4(1): 221228. 4 王高雄 , 周之铭 等 . 常微分方程 M. 北京 : 高等教育出版社 , 2006, 01: 131135. 5 丁慧 , 王亚男 . 从实践教学中谈常微分方程的发展及其应用 J. 科学时代 , 2010, 4(1): 121123. 6 赵静 , 但琦等 . 数学建模与数学实验 M. 北京 : 高等教育与出版社 . 2008, 01: 2631. 7 查淑玲 . 传染病的 SIR 模型 J. 山西中医学院学报 , 2003, 4(2): 5258. 8 黄其春 . 亚健康的产生及解决对策 J. 广西中医学院学报 , 2002, 03: 3238. 9 王育学 . 亚健康问题纵横谈 J. 解放军健康 , 2005, 01: 5561. 10 阳凌云 , 符云锦 . 一阶线性微分方程组的解法新探 J. 湖南工业大学学报 , 2010, 1(1): 6872
