一类分数阶HIV模型的稳定性分析.pdfVIP

第31卷第 1期 佳 木 斯 大 学 学 报 ( 自然 科 学 版 ) Vo1．31 No．1 2013 年 01月 JournalofJiamusiUniversity(NaturalScienceEdition) Jan． 2013 文章编号：1008—1402(2013)01-0129—02 一 类分数阶 HIV模型的稳定性分析① 庄科俊 (安徽财经大学统计与应用数学学院应用数学研究所。安徽 蚌埠 233030) 摘 要： 对一类分数阶的HIV病毒的CIM T细胞模型，利用分数阶的Routh—Hurwitz准则，研 究了平衡点的稳定性，并给出了相应的数值例子． 关键词： 分数阶；稳定性；平衡点 中图分类号： 0175．1 文献标识码：A 其中 ∈(0，1]，D 是 阶Caputo微分算子 ． 0 引 言 Arafa等人利用广义Taylor公式与同伦分析法，研 传染病学中数学模型的引入，为理解疾病传播 究了系统(2)的简化模型的解㈨．而本文的主要 机制及其控制策略提供 了一个全新的视角…． 目的在于，研究分数阶系统 (2)的平衡点的稳定 Tuck与Wan建立了如下的HIV型种群动力学模 性，并通过数值模拟验证理论分析的结果． 型 ： 1 稳定性分析 = s 一 一 由于分数阶导数不影响平衡点的位置，只需令 dy(t)= 一 口)， (1) 系统(2)的右端函数为0，即可得到系统的两个平 a 衡点：E1=(s／g，0，0)，E2=(2，Y2，7,2)，其中 2= = 一 一 = 音一 = 一 其中 与Y分别是t时刻未感染和已感染的CD4 T 细胞的数量，表示 t时刻血浆中病毒粒子的数 目． 系统中所有系数均为正数，s是CD4T细胞的生成 率，是人均死亡率，后是CD4 T细胞被病毒的感 染率，a是被感染细胞的人均消失率，C为感染细胞 与病毒的转化率，而 则是病毒粒子的死亡率．对 系统(1)，已有学者研究其简化系统的动力学性 0 质[3】、同伦摄动解[引． 最近，分数阶微积分引起了众多学者的广泛关 注 由于分数阶模型的非局部性、记忆性等性质， 图 1 n 0．95时系统 (2)的平衡点 渐近稳定 分数阶微分方程在生物系统的建模 中更加贴切实 系统(2)在平衡点 处的Jacobi矩阵为 际，并且出现了一些描述传染病动力学、种群动力 r一 0 一ks ] 学的分数阶模型-5]．因此，这里将考虑系统 (1) = 10 一o ks／g I 的分数阶情形： L 0 c — —ks J (￡)=s一 一kx= 相应的特征方程为 {D【Y(t)=kxz—ay (2) (A+ )[(A+0)(A+ +堕)一坐 ]=0 Da()： 一 一 ① 收稿 日期：2012—12—14 基金项 目：安徽省自然科学基金青年项 目(1208085Q