一类具免疫控制的SIR传染病模型的稳定.docx

一类具免疫控制的SIR传染病模型的稳定性
信息与计算科学专业学生：肖宪伟指导教师：宫兆刚
摘要：利用微分方程理论研究了具有免疫控制的数学模型，考虑总人口数是常数输入的影响，讨论了模型无病平衡点和地方病平衡点的存在性，利用特征值方法和Jacobi矩阵得到了无病平衡点和地方病平衡点的局部稳定性。构造Dulac函数的方法，得到了无病平衡点和地方病平衡点的全局稳定性充分条件，利用Matlab软件进行了数值模拟。
关键词：免疫控制;Jacobi；Dulac；平衡点；全局稳定性
1引言
面对传染病长期严峻的威胁和日益出现的新的疫情，其严重的危害着人类健康与社会经济的发展。又由于人们不能在人群中进行传染病的试验，因此，对各类传染病的流行趋势、发病规律的预测以及防治策略的重要性日益突出。根据疾病的发生、发展以及与之有关的阐述流行过程的特征，利用动力学的方法来研究传染病模型是十分重要的，目前对传染病的研究方法主要有描述性方面的研究、理论性方面的研究、分析性方面的研究和实验性方面的研究。传染病动力学[1]是对传染病进行理论性定量分析的一种非常重要的方法，通过对动力学性态的定性分析和模拟实验[2]，来显示疾病的发展过程，揭示起流行规律，分析疾病流行的原因和关键因素，预测其变化发展趋势，为预防和控制的最优策略提供了有力的理论依据。
在早期的传染病动力学中大多数传染病模型都是假设种群的总是常数状态而保持人口数不变，而没有考虑到其它方面的因素，但这种假设仅存在于一些环境状态封闭，人口的生育率和自然死亡率相平衡，且不考虑其它各方面等因素的理想状态下成立。随着传染病模型的不断发展和研究的不断深入，对各方面因素做了大量的研究，极大地丰富了传染病动力学理论。程晓云，胡志兴等在2007年考虑了具有阶段结构因素研究了一类具有阶段结构的自治传染病模型的稳定性[3]；徐为坚研究了一类具有种群Logistic增长饱和传染率的SIS模型的稳定性和Hopf；杜艳可，徐瑞，段立江在经典的传染病模型上考虑了标准发生率[4-5]的因素，研究了一类具有标准发生率的传染病模型的全局稳定性；李健全，马知恩研究了一类带有一般接触率和常数输入的流行病模型的全局分析；付景超等在2008年研究了一类具有垂直传染和连续预防接种的SIRS传染病模型[6-8]，得出了垂直传染和连续预防接种的稳定性分析；徐文雄，张仲华等研究了一类具有预防接种免疫力的双线性传染率SIR流行病模型全局稳定性；高淑京，滕志懂在2008年研究了一类具有饱和传染力和常数输入的SIRS脉冲接种模型研究[9-11]。
2具有免疫控制的SIR模型
2.1模型的建立
本文将基于经典的具有常数输入率的SIR模型，建立一类具有免疫控制的SIR传染病模型，将人群分为易感染者（Susceptible）、感染者（Infective）、移出者（Removed）三类，则所研究的数学模型如下：
（1）
其中，，分别表示t时刻易感染者、感染者、移出者的数量，，，分别表示各阶段的死亡率，表示接种率，表示传染率，r表示移出率，系统中的所有参数均为正值。
系统（1）的前两个方程不依赖于第三个方程，因此本文中仅考虑由系统（1）的前两个方程构成的系统为：
（2）
对系统(2)作变换，仍记为，则(2)化为：
（3）
考虑到系统（1），（2），（3）的实际上的生物意义，其S，I只能为非负数，因此本文只在区域中讨论问题。
2.2平衡点的存在性
对系统（3）我们令：
，当时，得出无病平衡点；
令时，当时，得出地方病平衡点，其中：，，，，。
2.3平衡点的稳定性分析
定理2当时，无病平衡点是局部渐近稳定的。
证明：系统（3）在处的Jacobi矩阵为：
，
其中：，，
，，
则系统（3）所对应的特征方程为：，
其中：

当时，有，，所以其特征根为负实根，从而无病平衡点是局部渐近稳定的。
定理3令，当，时，地方病平衡点是局部渐近稳定的。
证明：系统（3）在地方病平衡点处的Jacobi矩阵为：
，
其中：
，
，
其中：，则系统（3）所对应的特征方程为：，
其中：

，

，
令，当，时，有，，所以其特征根为负实根，从而地方病平衡点是局部渐近稳定的。
2.4平衡点的全局渐近稳定性分析
定理4当时，无病平衡点在区域内是全局渐近稳定的。
证明：当时，系统（3）在区域
内仅存在唯一的一个无病平衡点，并且可以得出平衡点在其边界上，所以在区域W内不存在有闭轨线，而且系统（3）从区域W内出发的轨线都不会超出W，又考虑到区域W的有界性，则对任给区域W内的一个初始值，系统（3）的满足初始值的解（S,I）最终都将趋向于平衡点，又因为平衡点的局部渐近稳定性，可得出是全局渐近稳定的。这表明，在所给群体中无论初始值的染病者会有多少，传染病都将不会流行且会逐渐消失。
定理5当且成立时，地方病平衡点
在区