数学物理学报 ›› 2026, Vol. 46 ›› Issue (5): 1913-1931.

• • 上一篇    下一篇

具有CTL免疫调控的病毒感染模型的稳定性与分支分析

宋晨玮1(), 徐瑞2,*()   

  1. 1 晋中学院数学系 山西晋中 030619
    2 山西大学复杂系统研究所 太原 030006
  • 收稿日期:2025-09-15 修回日期:2026-04-30 出版日期:2026-10-26 发布日期:2026-09-03
  • 通讯作者: 徐瑞 E-mail:scwei1014@163.com;xurui@sxu.edu.cn
  • 作者简介:宋晨玮, E-mail:scwei1014@163.com
  • 基金资助:
    国家自然科学基金(11871316);国家自然科学基金(12501677);山西省自然科学基金(202503021212261);山西省高等学校科技创新项目(2025L117)

Stability and Bifurcation Analysis of a Virus Infection Model with CTL Immune Response

Chenwei Song1(), Rui Xu2,*()   

  1. 1 Department of Mathematics, Jinzhong University, Shanxi Jinzhong, 030619
    2 Complex Systems Research Center, Shanxi University, Taiyuan 030006
  • Received:2025-09-15 Revised:2026-04-30 Online:2026-10-26 Published:2026-09-03
  • Contact: Rui Xu E-mail:scwei1014@163.com;xurui@sxu.edu.cn
  • Supported by:
    NSFC(11871316);NSFC(12501677);Natural Science Foundation of Shanxi Province(202503021212261);Technological Innovation Programs of Higher Education Institutions in Shanxi Province(2025L117)

摘要:

该文基于细胞毒性T淋巴细胞(CTL)免疫应答机制和HIV 病毒感染的细胞内潜伏期, 建立了一个具有CTL 免疫反应和胞内时滞的HIV 感染动力学模型. 首先, 定义了两个关键阈值参数: 免疫失活再生数$\mathcal{R}_{0}$ 和免疫激活再生数$\mathcal{R}_{1}$. 其次, 在理论分析方面, 通过构造合适的Lyapunov 泛函, 证明了当$\mathcal{R}_{0}<1$ 时, 未感染平衡点是全局渐近稳定的. 进一步结合波动引理和Lyapunov 泛函方法, 建立了当$\mathcal{R}_{0}>1>\mathcal{R}_{1}$ 时, 免疫失活平衡点的全局渐近稳定性. 此外, 研究还发现, 当$\mathcal{R}_{1}>1$ 时, 系统在免疫激活平衡点附近展现出丰富的动力学特性: 无论是否引入胞内时滞, 模型均会出现Hopf 分支, 且在特定参数条件下可能出现更为复杂的双Hopf 分支. 最后, 在数值模拟部分, 针对无胞内时滞的情形, 通过规范性计算精确确定了分支方向、分支周期解的稳定性、振幅及周期等动力学性质; 对于存在胞内时滞的情形, 借助双参数分支分析验证了双Hopf 分支的存在性.

关键词: HIV 感染模型, CTL 免疫, 胞内时滞, 再生数, Hopf 分支

Abstract:

Based on the immune response mechanism of cytotoxic T lymphocytes (CTL) and the intracellular latency of HIV infection, a dynamic model of HIV infection with CTL immune response and intracellular delay is established. Firstly, two key threshold parameters are defined: the immune inactivation reproduction number $\mathcal{R}_{0}$ and the immune activation reproduction number $ \mathcal{R}_{1}$. Second, in theoretical analysis, by constructing a suitable Lyapunov functional, it is proved that the uninfected equilibrium is globally asymptotically stable when $\mathcal{R}_{0}<1$. Further, by combining the wave lemma and Lyapunov functional method, the global asymptotic stability of the immune inactivation equilibrium is established when $\mathcal{R}_{0}>1>\mathcal{R}_{1}$. In addition, when $\mathcal{R}_{1}>1$, it is found that the system exhibits rich dynamic behaviors near the immune activation equilibrium: Hopf bifurcation will occur in the model regardless of the intracellular time delay, and even more complicated double Hopf bifurcation may occur. Numerical simulations quantitatively characterize the dynamical properties, employing normalization methods to precisely determine the bifurcation direction, stability of the periodic solutions, and their amplitude and period. In the presence of intracellular delay, a two-parameter bifurcation analysis rigorously establishes the existence of a double Hopf bifurcation.

Key words: HIV infection model, CTL immunity, intracellular delay, reproduction number, Hopf bifurcation

中图分类号: 

  • O175.23