一类具免疫控制的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 传染病 dt dt 3 模型，将人群分为易感染者(Susceptible)、感染者(Infective)、移出者(Removed)三类， 则所研究的数学模型如下： 〈 = - u2 I - rI dR = rI - u R + u 〈 = - u2 I - rI dR = rI - u R + uS (1) 其中 S (t ) ， I (t ) ， I (t )分别表示 t 时刻易感染者、感染者、移出者的数量， u ，u ， u 分 1 2 3 别表示各阶段的死亡率， u 表示接种率， b 表示传染率， r 表示移出率，系统中的所有参数 均为正值。 系统(1)的前两个方程不依赖于第三个方程，因此本文中仅考虑由系统(1)的前 两个方程构成的系统为： (| dS = B - u S - b SI - uS 〈|dtdI b SI 1 1 + a I dt 1 + aI 2l = dt 1 + aI 2 对系统(2)作变换dt = (1 + aI )dt ，仍记 dt 为dt ，则(2)化为： (| dS = B(1 + aI ) - u S (1 + aI ) - b SI - uS (1 + aI ) = P (S , I ) 〈 dt 1 l = l = b SI - u I (1 + aI ) - rI (1 + aI ) = Q(S , I )? dt 2 (2) (3) 考虑到系