[论文解读] A note on Exact solution of SIR and SIS epidemic models

Ghulam Shabbir, Haroon Ahmed Khan+1 more|arXiv (Cornell University)|Dec 22, 2010

Mathematical and Theoretical Epidemiology and Ecology Models参考文献 12被引用 26

一句话总结

本文通过直接积分技术，为恒定人口下SIR与SIS流行病模型的特定情形提供了精确的解析解。作者在无需限制参数条件的情况下，推导出两种模型的闭式解，显著推进了数学流行病学的发展，使在一般初始条件下对疾病动态的精确分析成为可能。

ABSTRACT

In this article we have successfully obtained an exact solution of a particular case of SIR and SIS epidemic models given by Kermack and Mckendrick [1] for constant population, which are described by coupled nonlinear differential equations. Our result has no limiting conditions for any parameter involved in the given models. In epidemiology many researchers believe that it is very hard to get an exact solution for such models. We hope this solution will be an opening window and good addition in the area of epidemiology.

研究动机与目标

• 为恒定人口下的SIR与SIS流行病模型获得精确的解析解，解决数学流行病学中长期存在的挑战。
• 克服描述疾病传播的非线性耦合微分方程常见的不可积性难题。
• 提供一种不受模型参数限制条件影响的解框架，提升模型的适用性与可解释性。
• 通过精确的函数形式为疾病动态的更深层次定性分析提供基础。

提出的方法

• 对于SIS模型，作者利用守恒律 s(t) + i(t) = k 将系统约化为仅含 i(t) 的单个非线性常微分方程，随后通过代换 y = i⁻¹ 求解。
• 在将 y = i⁻¹ 代入约化后的常微分方程后，通过直接积分得到SIS解，其形式为含积分常数C的指数函数。
• 对于SIR模型，通过将方程相加导出守恒律 s(t) + i(t) = 1 + C e⁻ᵘᵗ，将系统约化为仅含 i(t) 的单个常微分方程。
• SIR解采用代换 z = i⁻¹，并对 e⁻ᵘᵗ 进行级数展开，忽略高阶项以简化积分过程。
• 两种模型均利用初始条件求解积分常数，最终得到 s(t) 与 i(t) 的闭式表达式。
• 该方法依赖于直接积分与代换技术，而非数值或微扰方法，确保了解的精确性。

实验结果

研究问题

• RQ1 在恒定人口条件下，能否为SIR与SIS模型推导出无参数限制的精确解析解？
• RQ2 如何通过解析积分方法精确求解这些经典流行病模型的非线性耦合微分方程？
• RQ3 守恒律 s(t) + i(t) = k（或 s(t) + i(t) = 1 + C e⁻ᵘᵗ）在实现精确解中起到何种作用？
• RQ4 精确解在多大程度上能提升对疾病动态的理解，并促进流行病学中的模型验证？

主要发现

• SIS模型的精确解为 i(t) = β / (r + βC e⁻ᵝᵗ)，其中 β = rk - α，C 由初始条件 i(0) = i₀ 确定。
• SIS解精确满足 s(t) + i(t) = k，验证了总人口的守恒性。
• 对于SIR模型，解表示为 i(t) = λ / (β + λD e⁻ᵘᵗ e^(βC/μ))，其中 λ = β - μ + βC，且 C = s₀ + i₀ - 1。
• SIR解通过初始条件 s(0) = s₀ 与 i(0) = i₀ 完全确定积分常数 D 与 C。
• SIS情形中未使用近似，而SIR情形中为简化积分采用级数展开（e⁻ᵘᵗ ≈ 1 - μt），忽略高阶项。
• 所推导的解对所有正参数 r, α, β, μ 及初始条件 s₀, i₀ > 0 均有效，且无任何限制条件。