一类具有病例失踪的结核病模型的全局稳定性..doc

具有病例失踪的结核病模型的全局稳定性 肖雪霞, 郭玉玲, 刘璐菊 （河南科技大学数学与统计学院 河南洛阳 471003） 摘 要：建立了一类数学模型定义了模型的基本再生数通过构造适当的Lyapunov函数模型解的渐近性态. 证明了当基本再生数小于时无病平衡点是全局渐近稳定的当基本再生数大于时唯一的地方病平衡点是全局渐近稳定的. 关键词：模型；平衡点；基本再生数； Lyapunov函数；全局渐近稳定 中图分类号：O175 文献标识码：A 基金项目：国家自然科学基金项目河南科技大学博士启动基金项目 作者简介：，女，河人，. 结核病是由结核分枝杆菌感染所引起的一 种古老的慢性传染病,它可以影响全身多个器官,尤其是肺部.目前世界上大约有三分之一的人口已经受结核杆菌感染,并在每一秒钟产生一个新的结核杆菌感染者,每年全球约有一亿个新的结核杆菌感染者,全球现有两千万肺结核病人,每年新发结核病人约一千万,与此同时,每年约有三百万人死于结核病[8]. 我国是全球22个结核病高负担国家之一, 据2000年全国流行病学调查统计[9]：约有5.5亿人感染结核杆菌；传染性结核病人150万，患病人数仅次于印度，居世界第二位；每年约有12万人死于结核病；约有43.6%的结核病人在治疗过程中失踪(中断治疗),结核病例失踪最直接的后果是产生耐药性,从而使得耐药结核病例数目剧增,而耐药结核病已经被认为是上世纪80年代末90年代初我国结核病死灰复燃的三大原因之一，并且我国的耐药肺结核病例数目以11.2万而高居世界第二位.现有的统计数据表明，结核病已成为全世界由单一致病菌引起的死亡最多的疾病[8]，结核病严重地威胁着人类的健康和社会的发展.耐药结核病的研究是医学、公共卫生和应用数学的一个热点问题. 目前医学上尚未研制出有效的疫苗或药物 来预防和控制结核病的感染，只能通过对结核病人进行合理治疗,降低结核病人的传染期性,缩短传染期，延长结核病人的寿命来控制结核病.由于结核病例失踪等原因直接导致耐药或更严重的广泛耐药结核病的出现,而且由于缺乏新药、敏感药,再加上耐药结核病的治疗周期长,从而使得耐药结核病的治疗极其困难、治愈率也非常低、对健康人群造成的威胁也大[10].非常有意义的问题就是如何降低结核病例的失踪率、提高结核病例的规则治愈率、减少耐药结核病例的数目.国外已经有一些关于病例失踪的结核病模型的研究,已经有少量关于病例失踪随着时间变化的流行病学模型[3,7],而在我国,有关结核病例失踪的研究几乎是一片空白,也没有大规模和长期完整的数据,在本文中通过我国结核病的实际情况来建立一个简单的具有病例失踪的结核病模型,通过构造Lyapunov函数,分析模型的动力学性态. 1 建立模型 时刻的易感者类(), 时刻的携带结核杆菌的潜伏者类(), 时刻的具有传染性的结核病人类(), 时刻的失踪者类(), 和时刻的恢复者类(). 具有传染性的结核病人在经过一段抗结核药物治疗以后病情会好转, 部分病人会由于经济困难, 症状好转自行停药或药物反应等原因而中断治疗, 把这部分病人称为失踪者类, 之后这部分病人会在人群中自由流动, 对易感者可能造成感染或者会发展成为耐药结核病人. 易感者会被结核病人或失踪者感染, 有一部分人会在短时间内发病而成为结核病人; 另一部分会携带结核杆菌而进入潜伏者类. 潜伏者可能由于免疫系统的减弱等原因其体内的结核杆菌会再次活跃, 而发展成为结核病人; 也可以对结核杆菌携带者进行治疗而使其成为恢复者. 部分结核病人由于有效规则治疗而成为恢复者; 而另一部分会由于各种原因未能够痊愈而进入失踪者类, 失踪者类不会再次受到抗结核药物治疗, 在人群中自由活动, 他们可能会由于身体免疫系统的增强使体内的结核杆菌达到稳定状态而不再排出结核杆菌, 从而进入潜伏者类. 假设恢复者不会再次被感染而永久免疫.这五类人之间相互转移的转移框图如图1. 表示时刻总人口数目, 从而有 . 图1 具有病例失踪的结核病模型各仓室的转移框图 参数表示人口的输入率, 表示人口的自然死亡率, 表示传染率系数, 表示失踪者与结核病人对健康人的传染性之比(假设结核病人对健康人的传染性为1), 表示易感者被结核病人或失踪者感染之后快速进展为结核病人的比例, 而表示转化为携带结核杆菌的潜伏者的比例, 表示潜伏者转化为结核病人的内源复燃率, 表示对潜伏者的治疗率, 表示对结核病人的治愈率, 表示结核病人的失踪率(治疗中断率), 分别表示结核病人和失踪者的因病死亡率, 表示失踪者类的自愈率, 这里假设每个参数都是正数. 本文采用双线性发生率, 根据图1中各仓室的转移情况，此模型可以由下面的常微分方程组来 (1) 由于恢复者