艾滋病传播模型的研究摘 要:艾滋病是人体的免疫系统被艾滋病病毒破坏,使人体对威胁生命的各种病原体丧失了抵抗能力,从而发生多种感染或肿瘤,最后导致死亡的一种严重传染病。国际医学界至今尚无防治艾滋病的有效药物和疗法。因此,做好艾滋病的有效预防和控制应是我们抗击艾滋病最有效的手段。本文通过建立参数规划数学模型以Matlab 软件包为工具平台研究艾滋病的传播过程及流行趁向,期望能为政府做好新时期艾滋病预防控制工作提供理论参考。关键词:艾滋病;传播模型;参数规划;MATLAB中图分类号: 文献标识码: A0. 引言艾滋病在世界范围内的传播越来越迅猛,严重威胁着人类的健康和社会的发展,已成为威胁人们健康的第四大杀手。联合国艾滋病规划署2006 年5 月30 日宣布自1981 年6 月首次确认艾滋病以来,25 年间全球累计有6500 万人感染艾滋病毒,其中250 万人死亡。尤其忧虑的是,全世界约95%的艾滋病患者来自防治能力薄弱的发展中国家,如非洲、南亚、东南亚、中美洲等地。在我国,《中国艾滋病防治联合评估报告(2007)》显示,截至2007年底,现存艾滋病毒感染者和病人约70 万。疫情处于总体低流行、特定人群和局部地区高流行态势,性传播逐渐成为HIV 主要传播途径,意味着艾滋病防控形势更加严峻。未来我国艾滋病的流行是趁向平稳还是快速增长,取决于能否大面积地积极开展艾滋病预防活动以及提供有效的治疗。建立数学模型研究流行病的发展机理和传播过程,已有一个多世纪的历史,艾滋病出现以后,更引起了生物数学家们的注意,且在这方面做了较多的研究。本文就是在前人的研究基础上,通过建立参数规划数学模型借助Matlab 软件包为工具对一个艾滋病传播模型的探讨,期望能为政府做好我国新时期艾滋病预防和控制工作提供理论参考。1. 艾滋病简介艾滋病的医学中文全名为“获得性免疫缺损综合症”,英文全名为“Acquired ImmunoDeficieney Syndrome”,简称AIDS,它是由艾滋病(中英文全名为“人体免疫缺损病毒”、“Human Immunodeficiency Virus”,简称HIV)引起的。这种病毒终身传染,破坏人的免疫系统,使人体丧失抵抗各种疾病的能力,从而严重危害人的生命[1]。HIV 进入人体后,它就把人体免疫系统中最重要的CD4+T 淋巴细胞作为攻击目标,大量吞噬、破坏CD4+T 淋巴细胞,从而破坏人的免疫系统,最终使人体免疫系统崩溃,使人体因丧失对各种疾病的抵抗能力而发病并死亡。AIDS 从感染到发作的进程大致可分3 个阶段:感染初期、潜伏期、发作期。在感染初期,HIV 进入人体的感染巨噬细胞,将病毒带到局部淋巴结,引起各种急性症状,接着,CD8+T 淋巴细胞、抗体做出反应,从而控制疾病的发展,血液中的HIV 持续在一个较稳定的水平,疾病进入潜伏期。随着HIV 对巨噬细胞、CD4+T 细胞等的感染,免疫系统逐步被破坏,被感染的CD4+T 细胞由裂解而大量产生HIV。1本课题获广西教育科研立项项目“离散空间的模糊多属性决策理论与方法研究”()的资助。- 2 -当正常的CD4+T 细胞急剧减少、HIV 迅速增加时,病情发作,随时出现的各种病原体都可能引起感染,病人最后因各种机能衰竭而死亡。一般在未治疗情况下,AIDS 从感染到发作的平均时间约为10 年,但是不同的病人的差异较大,临床上可观察到2 到18 年的潜伏期,这主要取决于CD4+T 细胞浓度下降和HIV浓度上升的速度。通常健康人每1mm3 血液中平均有1000 个CD4+T 细胞,当HIV 的携带者的CD4+T 细胞降至200 个/mm3 时,疾病发作。2. 艾滋病传播模型艾滋病传播途径主要有性传播、血液传播、共用针具的传播和母婴传播等四种,其中性传播已成为当今艾滋病传播的主要途径。故下面的模型主要研究通过性活动引起AISD 的传播 ,通过其他因素引起的传播可以建立类似的模型。 模型建立将目标人群(具有性活动者)分为3 类,x(t)为t 年易感染的人数,y(t)为被HIV 感染的人数,z(t)为已患AIDS 的人数,在没有特定目标的情况下,假定x,y,z 的初始值分别为15×106,3×106,×106。参照其它传染病的传播模型及参数规划模型案例得到模型为[2],[3]( ) , 1 s c xdtdx = − λ + μ (1)( ) , 2 c x ydtdy = λ − γ + μ (2)( ) , 3 y zdtdz = γ − μ +δ (3)其中各个参数的定义及其定值如下:s ~易感染者加入目标人群的速率(106/年);c ~获得新的性伴侣的平均速率(2/年);λ ~性伴侣被HIV 感染者的概率(); 1μ ~易感染者退出目标人群的比例(年);γ ~HIV 感染者进入AIDS 的比例(年); 2 μ ~HIV感染者退出目标人群的比例(年); 3 μ ~AIDS 退出目标人群的比例(年);δ ~AIDS 死亡的比例(年)。虽然线性常系数微分方程(1)~(3)有解析解,可是我们只想了解数值结果和观察直观的变化趁势,于是在上面的参数和初值下运用 软件包[4]求数值解便得到了图1。- 3 -可以看出,人们最关心的被HIV 感染的人数y 在约5 年时达到最大值约1200 万,20年后趁向稳定值约750 万。 模型分析由劳斯-霍尔维茨(Routh-Hurwitz)判据,容易得到方程(1)~(3)的唯一平衡点是,1*λ + μ=cx s * ,2y* c xγ μλ+= *3z* yμ δγ+= 。 (4)因为方程(1)~(3)的3 个特征根均为负值,所以平衡点是全局稳定的,与初始值无关。稳态下HIV 感染者在目标人群中的比例为λγ μμ δγβcx y zy23* * **11++++=+ += (5)当一个参数值增加时引起平衡点3 个坐标及β 值的变化见下表1(如γ 增加时,x* 不变,y*减少, z* 增加,β 减少)。表1 一个参数值增加时引起3 个坐标及β 值的变化情况参数 x* y* z* βs ↑ ↑ ↑  ̄λ ↓ ↑ ↑ ↑c ↓ ↑ ↑ ↑γ  ̄ ↓ ↑ ↓图1 方程(1)~(3)的数值解- 4 -1 μ↓ ↓ ↓  ̄2 μ  ̄ ↓ ↓ ↓3 μ  ̄  ̄ ↓ ↑δ  ̄  ̄ ↓ ↑由于平衡点的全局稳定性,所以只要采取适当措施使各个参数向正确的方向(增加或减少)改变,长期说来就可以使HIV 感染者和AISD 的人数减少,而不管目前情况如何。 接种疫苗模型用弱化的HIV 作疫苗帮助人体建立对病毒的免疫性,是一种人为的干预措施,为了建立这种情况模型,需要增加两个函数:目标人群中接种疫苗的人数( ) 1 y t ,接种疫苗后又被病毒感染的人数( ) 2 y t ,假定1 y , 2 y 的初始值分别为1000 和0。模型为(1 ) ( ) , 1 1 p s c c xdtdx = − − λ + λ + μ (6)( ) , 2 c x ydtdy = λ − γ + μ (7)(1 ) ( ) , 1 1 1 4 11 ps c x c y ydtdy = + λ − −ϕ λ − γ + μ (8)(1 ) ( ) , 1 2 5 22 c y ydtdy = −ϕ λ − γ + μ (9)( ) , 1 1 2 2 3 y y y zdtdz = γ +γ +γ − μ +δ (10)其中新增加的参数及其设定值如下: p ~目标人群接种疫苗的比例(); 1λ ~性伴侣接种疫苗的概率();ϕ ~接种疫苗后起到预防作用的概率(); 1 γ ~接种疫苗者进入AIDS 的比例(年); 4 μ ~接种疫苗人群退出目标人群的比例(年); 2 γ ~接种疫苗者被病毒感染进入AIDS 的比例(年); 5 μ ~接种者被病毒感染退出目标人群的比例(年)。虽然接种疫苗还不能预防HIV 的初期感染,但是我们可以作简单的计算来预测如果疫苗取得进展后的效果[4]。 比如在p =,ϕ = 的情况下利用 软件包求(6)~(10)的数值解即得到图2,从图2 可以看出,被HIV 感染的人数y 在约3 年时达到最大值约600 万,比模型(1)~(3)的结果减少了一半,20 年后趁向稳定值约200 万,比模型(1)~(3)的结果减少得更多。- 5 结论分析上面这个艾滋病传播模型概念上虽然较简单,但涉及到的参数很多,对于这些模型来说,关键在于如何确定其中的参数。应用它的主要困难是确定用于特定国家或地区的一组参数,虽然乌干达等非洲国家已经在这方面做了大量的统计研究,但是目前还不能得到确定这些参数的较有效方法,另外,在临床上较精确地得到被HIV 感染的和已患AIDS 的人数也是困难的。3. 结束语尽管目前我国艾滋病的疫情上升速度有所减缓,还没有出现艾滋病大规模流行之势,但是我们要清醒地看到,疫情存在潜在的流行趁势,HIV 的传播途径已演变成以性传播途径为主,已经与国际上的流行趁势一样,艾滋病疫情地区分布差异大,艾滋病流行因素广泛存在,局势越来越严峻,一触即发,并可能出现灾难性后果。因此,从现在到本世纪末将是我国预防控制艾滋病的关键时期,如果我们现在不积极采取预防控制措施,将错失良机。当务之急是要全面了解我国艾滋病传播途径的变化、流行趁势、受影响人群有关的危险行为等情况,建立一个有效的监测系统,为决策者提供有关艾滋病传播的准确数据,预测艾滋病流行疫情和趁势,为全国艾滋病规划策略的制定提供依据。随着艾滋病在我国不同地区不断蔓延扩大,其流行模式将变得越来越复杂。因此,我们的监测系统不仅仅是数据的收集,而应当注重数据的分析以帮助对策的制定。图2 方程(6)~(10)的数值解- 6 -参考文献[1]曹毅.警惕艾滋病[M] .北京:清华大学出版社,2005.[2]谭永基,蔡志杰. 数学模型[M].上海:复旦大学出版社,2005. 310~320.[3] Christelle,Christian Prins, Marc Sevaux.运筹学案例[M] .北京:北京林森科技发展有限公司,2007.[4]赵东方. 数学模型与计算[M]. 北京:科学出版社, of the spread of AIDS on the basis of estimationprogrammingZhu JiarongNanning Teacher’s College,Department of Mathematics & Computer Science, P. Province, CongZuo(532200)AbstractAIDS is a serious infectious disease,it’s caused by HIV infection, which damage the body’s immunesystem and make the body losing their resistance to various life-threatening pathogens. Internationalmedical profession has no effective drugs and treatments of preventing or healing AIDS. Therefore, doa good job in AIDS prevention and control should the most effective means to fight AIDS. Our paperdiscussed the process of AIDS’ spreading by establishing math model about estimation programmingbased on Matlab, and except the government can take it as a reference for AIDS prevention in the newera.