Bifurcation analysis for a SIRS model with a nonlinear incidence rate

基本信息

• 论文ID: 2510.27292
• 标题: Bifurcation analysis for a SIRS model with a nonlinear incidence rate
• 作者: Xiaoling Wang, Kuilin Wu (贵州大学数学与统计学院)
• 分类: math.DS (动力系统)
• 发表时间: 2025年10月31日 (arXiv预印本)
• 论文链接: https://arxiv.org/abs/2510.27292

摘要

本文深入研究了具有一般非线性感染率 $f(I)S=\beta I(1+\upsilon I^{k-1})S$ ($k>0$) 的SIRS流行病模型。作者系统分析了平衡点的存在性和稳定性，应用局部分岔理论探索了模型丰富的动力学行为。文章推导了不同类型分岔的正规形式，包括Bogdanov-Takens分岔、幂零焦点分岔和Hopf分岔。通过计算前四个焦点值来确定Hopf分岔的余维数，揭示了系统可能产生多个极限环的复杂动力学现象。数值模拟结果验证了理论分析的正确性。

研究背景与动机

研究问题

传染病动力学中的分岔现象是指由环境、宿主或病原体因素引起的传播动态或流行行为的突变。这些分岔现象不仅决定传播的速度和范围，还可能改变病毒致病性，从而影响公共卫生干预措施的有效性。

问题重要性

1. 现实意义: 全球化加速和生态转型导致传染病传播模式日益复杂，分岔现象分析对于预测流行趋势、制定有效控制策略和高效配置公共卫生资源至关重要
2. 理论价值: 经典的双线性感染率模型 $\beta IS$ 只能描述理想和简单的传播过程，无法捕捉现实世界中的复杂动力学行为（如突然爆发、持续流行和周期性波动）

现有方法的局限性

1. 简单模型的限制: 经典双线性感染率模型通常最多只有一个地方病平衡点，无法表示双稳态或周期性等复杂现象
2. 特定参数研究: 以往对感染率 $f(I)S = \beta I(1+\upsilon I^{k-1})S$ 的研究主要集中在特定的 $k$ 值（如 $k=1$ 或 $k=2$）
3. 分岔分析不完整: 缺乏对一般参数 $k>0$ 情况下系统完整分岔结构的系统性研究

研究动机

通过引入饱和性和非线性特征，非线性感染率能够捕捉由心理行为、资源约束和超级传播等复杂因素引起的丰富动力学行为。本文旨在对一般参数 $k>0$ 的情况进行全面的分岔分析，揭示比以往研究更丰富的动力学现象。

核心贡献

1. 系统的平衡点分析: 对一般参数 $k>0$ 情况下，系统化地分析了平衡点的存在性、类型和稳定性，发现系统可能存在至多3个正平衡点
2. 高余维分岔分析:
   • 证明了系统在 $1<k\leq 2$ 时可发生余维3的Bogdanov-Takens分岔
   • 证明了系统在 $k>2$ 时可发生余维4的Bogdanov-Takens分岔
   • 分析了幂零焦点分岔（余维3和余维4）
3. Hopf分岔的精细分析:
   • 计算了前四个焦点值
   • 证明了系统可以产生至多4个极限环
   • 确定了不同参数区域下极限环的数量和稳定性
4. 完整的分岔图: 构建了详细的分岔图，展示了鞍结分岔、Hopf分岔、同宿轨分岔等多种分岔现象的共存和转换
5. 数值验证: 提供了大量数值模拟结果，直观展示了理论预测的各种动力学行为

方法详解

任务定义

研究具有非线性感染率的SIRS流行病模型：

其中：

• $S(t), I(t), R(t)$: 易感者、感染者、恢复者数量
• $\beta>0$: 传播系数
• $\upsilon>0, k>0$: 非线性感染率参数
• $b>0$: 人口补充率
• $d>0$: 自然死亡率
• $\mu>0$: 恢复率
• $\delta\geq 0$: 免疫丧失率

模型简化

通过总人口分析 $N(t)=S(t)+I(t)+R(t)$，得到 $\lim_{t\to\infty}N(t)=\frac{b}{d}\equiv\Lambda$。

简化后的二维系统：

通过变量变换 $I=\frac{d+\delta}{\beta}x$, $R=\frac{d+\delta}{\beta}y$, $t=\frac{1}{d+\delta}\tau$，得到标准化系统：

其中 $p=(d+\delta)^{k-1}\upsilon/\beta^{k-1}$, $\Lambda_0=\beta\Lambda/(d+\delta)$, $\gamma=(d+\mu)/(d+\delta)$, $\eta=\mu/(d+\delta)$。

基本再生数: $R_0 := \Lambda_0/\gamma$

平衡点分析方法

无病平衡点 $E_0(0,0)$:

• $R_0<1$ 时为稳定结点
• $R_0>1$ 时为鞍点
• $R_0=1$ 时为鞍结点

地方病平衡点: 通过分析函数 $H(x) = (1+px^{k-1})\left(1-\frac{1+\eta}{\Lambda_0}x\right) - \frac{1}{R_0}$ 的零点来确定。根据 $k$ 的不同取值范围，系统可能有0-3个正平衡点。

分岔分析技术

1. 鞍结点分岔

当 $H(x)=0$ 且 $H'(x)=0$ 时，平衡点退化为鞍结点。通过中心流形定理和正规形理论，证明了系统在参数变化时会经历鞍结点分岔。

2. Bogdanov-Takens分岔

当平衡点的Jacobi矩阵有双零特征值时，系统可能发生Bogdanov-Takens分岔。文章通过一系列坐标变换，将系统化为标准正规形：

并验证了非退化条件 $\frac{\partial(\eta_1,\eta_2,\eta_3,\eta_4)}{\partial(\lambda_1,\lambda_2,\lambda_3,\lambda_4)}\neq 0$。

3. Hopf分岔分析

当 $\text{Tr}(J(E))=0$ 且 $\text{Det}(J(E))>0$ 时，平衡点可能是弱焦点。使用形式级数方法计算前四个焦点值：

等等。

通过分析焦点值的符号变化，确定弱焦点的阶数和极限环的数量。

技术创新点

1. 统一的参数分析框架: 不同于以往针对特定 $k$ 值的研究，本文建立了适用于一般 $k>0$ 的统一分析框架
2. 高余维分岔的完整刻画:
   • 首次证明了该模型在 $k>2$ 时可发生余维4的Bogdanov-Takens分岔
   • 发现了余维4的幂零焦点分岔
3. 极限环数量的精确确定: 通过计算四个焦点值，证明了系统最多可以有4个极限环，这比文献中已知结果（2个极限环）更精确
4. 系统的分类讨论: 根据 $k$ 的不同范围（$0<k<1$, $k=1$, $1<k<2$, $k=2$, $k>2$），给出了完整的动力学分类

实验设置

数值方法

• 使用MATLAB进行数值模拟
• 采用四阶Runge-Kutta方法求解微分方程
• 使用延拓方法追踪分岔曲线

参数选择

文章针对不同的分岔现象，精心选择了参数组合进行数值验证：

Bogdanov-Takens分岔 ($k=2$):

• $\Lambda_0=1.8993$, $p=5.7966$, $\eta=2.3072$
• 通过改变 $\gamma$ 观察分岔序列

高余维分岔 ($k=3$):

• $\Lambda_0=3.1832$, $p=0.9331$, $\eta=0.6355$
• 展示三极限环共存现象

Hopf分岔:

• 多组参数展示1-4个极限环的情况

可视化方法

1. 相图: 展示轨线在相空间中的行为
2. 分岔图: 在参数空间中标识不同分岔曲线
3. 极限环: 直观展示多个极限环的共存

实验结果

主要结果

1. 平衡点的存在性和稳定性（第2节）

定理 2.1: 无病平衡点 $E_0(0,0)$ 的性质完全由基本再生数 $R_0$ 决定

引理 2.2-2.4: 根据 $k$ 的不同范围，系统可能有：

• $k\leq 1$: 至多1个正平衡点
• $1<k\leq 2$: 至多2个正平衡点
• $k>2$: 至多3个正平衡点

2. 退化平衡点分析（第3节）

定理 3.1: 鞍结点的存在条件

• 当 $1<k<2$ 时，$E_2^*$ 是鞍结点
• 当 $k=2$ 时，$\bar{E}_0^*$ 是鞍结点
• 当 $k>2$ 时，多个退化平衡点可能是鞍结点

定理 3.2: 尖点分岔的存在性 当特定参数条件满足时，退化平衡点是余维2的尖点

定理 3.3: 高余维尖点

• $1<k\leq 2$: 余维3的尖点
• $k>2$ 且 $z=\tilde{z}$: 余维4的尖点

定理 3.4: 幂零焦点分岔 在特定参数条件下，系统可发生余维3或余维4的幂零焦点分岔

3. 分岔分析（第4节）

定理 4.1: 鞍结点分岔 系统在参数曲面 $SN_1, SN_2, SN_3$ 附近发生鞍结点分岔

定理 4.2: Bogdanov-Takens分岔（$k>2$） 系统在 $(\Lambda_0,\gamma,\eta,p)$ 接近 $(\tilde{\Lambda}_0,\tilde{\gamma},\tilde{\eta},\check{p})$ 时发生余维4的Bogdanov-Takens分岔，包含：

• 余维2的尖点分岔
• 余维3的尖点分岔
• 鞍结点分岔、Hopf分岔、同宿轨分岔等

定理 4.3: Bogdanov-Takens分岔（$1<k\leq 2$） 系统发生余维3的Bogdanov-Takens分岔

4. Hopf分岔（第5节）

定理 5.1: 当 $R_i=0$ ($i=0,1,2,3$) 时，平衡点是至多2阶的弱焦点

定理 5.2: 当 $l_1=0$, $\frac{3}{2}<k<2$ 且 $k\neq 1.9839$ 时，平衡点可以是3阶弱焦点

引理 5.1: 当 $k=2$ 时，系统可以有至多2个极限环

定理 5.3: 当 $k>2$ 且 $R_i=0$ 时，系统可以有至多3个极限环

定理 5.4: 当 $k>2$ 且 $l_2=0$ 或 $l_4=0$ 时，系统可以有至多3个极限环

注记 5.1: 数值结果显示系统可以有4个极限环（见图5.5）

分岔序列的数值验证

图4.2 ($k=2$): 展示了参数 $\gamma$ 减小时的分岔序列：

1. 鞍结点分岔 ($\gamma=2.6731$)
2. Hopf分岔 ($\gamma=2.6717$)
3. 同宿轨分岔 ($\gamma=2.6712$)
4. 双极限环 ($\gamma=2.6687$)

图4.3 ($k=3$): 展示了更复杂的分岔序列：

1. 鞍结点分岔
2. 双极限环
3. 三极限环
4. Hopf分岔

极限环的数值证据

图5.1: 展示了当 $\frac{3}{2}<k<2$ 且 $l_1=0$ 时，系统可以有3个极限环

图5.2: 展示了 $k=2$ 时的1-2个极限环

图5.3: 展示了 $k>2$ 时的3个极限环

图5.4: 展示了 $k>2$ 且 $l_2=0$ 时的3个极限环

图5.5: 展示了4个极限环共存的情况（参数: $k=3$）

实验发现

1. 参数 $k$ 的关键作用: $k$ 的取值范围决定了系统的动力学复杂度：
   • $k\leq 1$: 动力学相对简单
   • $1<k\leq 2$: 出现余维3分岔
   • $k>2$: 出现余维4分岔和更多极限环
2. 极限环的多样性: 系统可以表现出1-4个极限环，这在流行病模型中是罕见的
3. 分岔的丰富性: 系统展现了几乎所有经典的分岔类型：鞍结点、Hopf、同宿轨、Bogdanov-Takens等
4. 基本再生数的影响: $R_0$ 的变化会导致无病平衡点和地方病平衡点之间的稳定性转换

相关工作

经典SIRS模型

• Kermack-McKendrick (1927): 提出了经典的SIR模型框架
• Hethcote & van den Driessche (1991): 研究了非线性感染率模型

饱和感染率模型

• Capasso & Serio (1978): 引入饱和感染率 $\frac{\beta IS}{1+\alpha I}$
• Xiao & Ruan (2007): 研究了非单调感染率 $\frac{\beta IS}{1+\alpha I^2}$
• Zhou et al. (2007): 研究了 $\frac{kIS}{1+\beta I+\alpha I^2}$

非线性感染率模型

• Ruan & Wang (2003): 研究了 $\frac{\beta I^2 S}{1+\alpha I^2}$，发现了余维2的Bogdanov-Takens分岔和2个极限环
• Zhang et al. (2022): 研究了一般饱和感染率 $\frac{\beta I^p S}{1+\alpha I^q}$
• Cui & Zhao (2024): 研究了一般参数 $p, q$ 的Bogdanov-Takens分岔

本文研究的感染率

• Lu et al. (2023): 研究了 $k=1$ 的情况，发现余维3的Bogdanov-Takens分岔
• Jin et al. (2007): 研究了 $k=2$ 的情况，发现余维2的Bogdanov-Takens分岔

本文的优势

1. 一般性: 对任意 $k>0$ 进行统一分析，而非特定值
2. 完整性: 发现了更高余维的分岔（余维4）
3. 精确性: 证明了最多4个极限环，比已知结果更精确
4. 系统性: 提供了完整的分岔图和动力学分类

结论与讨论

主要结论

1. 平衡点结构:
   • 系统总有唯一的无病平衡点 $E_0$
   • 当 $R_0>1$ 时，可能存在1-3个地方病平衡点
   • 平衡点的数量和类型由参数 $k$ 的范围决定
2. 分岔现象:
   • 系统可发生鞍结点分岔、Hopf分岔、同宿轨分岔
   • 存在余维3的Bogdanov-Takens分岔（$1<k\leq 2$）
   • 存在余维4的Bogdanov-Takens分岔（$k>2$）
   • 存在余维3和余维4的幂零焦点分岔
3. 极限环:
   • 系统最多可以有4个极限环
   • 极限环的数量由Hopf分岔的余维数决定
   • 不同参数区域可以观察到1-4个极限环共存
4. 流行病学意义:
   • 非线性感染率能够产生复杂的动力学行为
   • 系统可以表现出双稳态、周期振荡等现象
   • 参数 $\beta$（传播系数）在 $R_0$ 中起关键作用

局限性

1. 理论完整性:
   • 对于某些参数区域（如 $R_i\neq 0$ 且 $l_1\neq 0$），焦点值的完整分析尚未完成
   • 4个极限环的存在性仅有数值证据，缺乏严格的理论证明
2. 实际应用:
   • 模型假设人口总数趋于常数，可能不适用于人口快速变化的情况
   • 未考虑空间异质性和年龄结构
   • 参数估计的实际困难未讨论
3. 数值验证:
   • 某些高余维分岔的数值验证较为复杂
   • 未提供参数敏感性的系统分析
4. 模型假设:
   • 假设恢复者最终会失去免疫力，这可能不适用于所有疾病
   • 忽略了疾病引起的额外死亡率

未来方向

1. 理论完善:
   • 完成所有参数区域的焦点值分析
   • 严格证明4个极限环的存在性
   • 研究全局动力学行为，如全局稳定性和全局分岔
2. 模型扩展:
   • 考虑空间扩散效应（反应扩散系统）
   • 引入时滞效应
   • 考虑随机扰动（随机微分方程）
   • 加入年龄结构或异质性
3. 应用研究:
   • 将模型应用于具体疾病（如COVID-19、流感等）
   • 研究控制策略（疫苗接种、隔离等）的效果
   • 参数估计和模型验证
4. 计算方法:
   • 开发更高效的分岔分析算法
   • 构建参数空间的完整分岔图
   • 研究混沌和复杂动力学

深度评价

优点

1. 理论深度:
   • 运用了高级的动力系统理论（中心流形定理、正规形理论、分岔理论）
   • 计算了前四个焦点值，这在技术上非常困难
   • 发现了余维4的分岔，这在流行病模型中极为罕见
2. 系统性:
   • 对一般参数 $k>0$ 进行了统一分析
   • 提供了完整的平衡点分类和稳定性分析
   • 构建了详细的分岔图
3. 创新性:
   • 首次证明了该模型可以有4个极限环
   • 首次发现了余维4的Bogdanov-Takens分岔
   • 发现了多种高余维分岔的共存
4. 数值验证:
   • 提供了大量高质量的相图和分岔图
   • 数值结果与理论预测高度一致
   • 可视化清晰，易于理解
5. 写作质量:
   • 逻辑清晰，结构合理
   • 数学推导详细（虽然某些长表达式被省略）
   • 图表丰富，说明性强

不足

1. 理论缺口:
   • 某些参数区域的焦点值分析未完成（如第5.1.2节）
   • 4个极限环的严格存在性证明缺失
   • 全局动力学分析不足
2. 实际应用:
   • 缺乏与真实流行病数据的对比
   • 未讨论参数的生物学意义和实际可测量性
   • 控制策略的应用研究缺失
3. 技术细节:
   • 许多中间计算结果被省略（"omitted for brevity"）
   • 某些符号定义不够清晰
   • 数值方法的具体实现细节不足
4. 可读性:
   • 对于非专家读者，某些部分过于技术化
   • 缺乏对生物学意义的直观解释
   • 某些定理的条件过于复杂
5. 实验设计:
   • 缺乏参数敏感性分析
   • 未提供代码或数据的可复现性材料
   • 某些数值结果的参数选择依据不明确

影响力

1. 学术贡献:
   • 为流行病模型的分岔理论提供了重要补充
   • 发现的高余维分岔和多极限环现象具有理论价值
   • 为相关研究提供了新的分析框架
2. 实用价值:
   • 揭示了非线性感染率可能导致的复杂流行病动态
   • 为公共卫生决策提供了理论支持
   • 表明简单的控制策略可能无效（由于多稳态和周期振荡）
3. 可复现性:
   • 理论推导详细，可以验证
   • 数值模拟参数明确
   • 但缺乏代码和详细的数值方法说明
4. 引用价值:
   • 对研究非线性流行病模型的学者有重要参考价值
   • 对研究高余维分岔的数学家有方法论价值
   • 可能成为该领域的重要文献

适用场景

1. 理论研究:
   • 动力系统理论的应用
   • 分岔理论的发展
   • 流行病模型的数学分析
2. 流行病学应用:
   • 具有复杂传播机制的疾病（如超级传播事件）
   • 需要考虑心理因素和行为变化的疾病
   • 表现出周期性爆发的疾病
3. 政策制定:
   • 理解疾病传播的复杂动态
   • 评估干预措施的潜在效果
   • 识别可能的突变和临界点
4. 教学用途:
   • 动力系统课程的案例研究
   • 分岔理论的应用示例
   • 数学建模课程的教材

参考文献（精选）

本文引用了37篇重要文献，主要包括：

经典流行病模型:

• 28 Kermack & McKendrick (1927): 经典SIR模型
• 18 Levin et al. (1989): 应用数学生态学

非线性感染率:

• 14 Capasso & Serio (1978): 饱和感染率
• 15 Xiao & Ruan (2007): 非单调感染率
• 9 Ruan & Wang (2003): 非线性感染率的分岔

分岔理论:

• 32 Zhang et al. (1992): 微分方程定性理论
• 33 Khibnik et al. (1998): 三次Liénard方程的全局研究
• 34 Dumortier et al. (2001): 余维4幂零鞍点的展开

相关工作:

• 30 Lu et al. (2023): $k=1$ 情况的研究
• 31 Jin et al. (2007): $k=2$ 情况的研究
• 20-22 Zhang, Cui, Hu等: 一般饱和感染率的研究

总体评价: 这是一篇高质量的动力系统理论论文，在流行病模型的分岔分析方面取得了重要进展。理论深度和系统性是其最大优势，但在实际应用和可复现性方面还有改进空间。对于研究非线性动力系统和流行病数学模型的学者，这是一篇值得深入研究的文献。