國 會，胡志興

(北京科技大學(xué)數(shù)理學(xué)院，北京 100083)

艾滋病是由感染HIV 病毒引起的，HIV 病毒對人體的免疫系統(tǒng)有很強(qiáng)的攻擊性，至今仍沒有有效的藥物來治愈這類傳染病. 所以，近年來大量的數(shù)學(xué)模型來定量或定性的分析其傳染過程，這對弄清艾滋病的傳染機(jī)理以及對治療、控制艾滋病傳染具有重要的指導(dǎo)意義. 自Perelson 等人[1]提出經(jīng)典模型，之后很多學(xué)者對模型進(jìn)行了修改和完善，包括考慮潛伏期，CTL 免疫反應(yīng)的影響[2-5]. Wang Y 等人又在此基礎(chǔ)上考慮了時滯對HIV 病毒傳染的影響[6-9]，使模型更具有生物學(xué)意義. 很多HIV 數(shù)學(xué)模型都考慮的是飽和發(fā)生率，非線性發(fā)生率等等. Liao F 等人在此基礎(chǔ)上考慮了細(xì)胞免疫抑制的HIV 數(shù)學(xué)模型[10-12]， 但是卻忽略的細(xì)胞內(nèi)時滯以及CLT 免疫反應(yīng)時滯對動力學(xué)模型的影響.

本文研究了具有雙時滯且考慮CLT 免疫細(xì)胞抑制作用的HIV 病毒感染的動力學(xué)模型：

1 模型系統(tǒng)(1）的分析

系統(tǒng)(1)的初始條件滿足：

其中

下面討論解的非負(fù)性：由系統(tǒng)(1)我們可以得到

其中

2 平衡點的穩(wěn)定性

即

即

其中

的根.

因此，，

其中

根據(jù)文獻(xiàn)[13]的命題4.8.1，有如下引理

3 內(nèi)部平衡點的局部穩(wěn)定性和Hopf 分支

即

其中

由方程(3.2)得特征方程：

其中

因為

由Routh-hurwitz 判據(jù)可知，方程(3.3)僅有負(fù)實部的根當(dāng)且僅當(dāng)條件()成立. 所以()成立時，內(nèi)部平衡點是局部漸近穩(wěn)定的.

由方程(3.3)得特征方程：

其中

其中

定義

由文獻(xiàn)[14 -15]，得到下面引理.

引理3.1 如果方程(3.7)滿足如下條件.

若方程(3.7)滿足任一上述條件，則至少存在一個正實根.

則方程(3.7)無正實根.

不失一般性，我們不妨假設(shè)方程(3.7)有四個正實根，分別為然后方程(3.6)的四個正根分別為根據(jù)方程(3.5)得到

因此，

由方程(3.2)得特征方程：

其中

其中

由文獻(xiàn)[15]的定理2.1 得到如下引理.

引理3.2 對于方程(3.13)有如下條件：

若方程(3.13)滿足任一上述條件，則方程(3.13)至少存在一個正實根.

則方程(3.13)沒有正實根.

不失一般性，我們假設(shè)方程(3.13)有四個正實根，分別為而且方程(3.12)就有四個正實根

從方程(3.11)我們可以得到：

由方程(3.11)得：

從而

假設(shè)方程(3.15)有純虛根，帶入方程(3.15)，分離實部和虛部得：

解得：

其中

其中

其中

因此

其中

其中

其中

將方程(3.22)中的兩個方程平方后相加，得到：

其中

其中

其中

因此，

4 數(shù)值模擬

在這一部分中，通過MATLAB 對平衡點進(jìn)行了簡單的數(shù)值模擬.

對于情形1 選取參數(shù)的選取如下所示：

選取初值(9716，0.7961，0.0326，0.7822)，計算.當(dāng)時，所選參數(shù)滿足情形1，此時內(nèi)部平衡點是局部漸近穩(wěn)定的，如圖1 所示：

圖1 時，漸近穩(wěn)定Fig.1 is locally asymptotically stable when

在情形2 中選取如下參數(shù)：

圖2 時系統(tǒng)的相圖Fig.2 Phase diagram of the system for

在情形3 中選取參數(shù)如下：

圖3 時系統(tǒng)的相圖Fig.3 Phase diagram of the system for

在情形4 中選取參數(shù)如下：

圖4 時系統(tǒng)的相圖.Fig.4 Phase diagram of the system for .

圖5 時系統(tǒng)的相圖Fig.5 Phase diagram of the system for

5 結(jié)論

本文建立了具有免疫抑制及兩時滯的HIV 病毒感染模型. 主要由Routh -Hurwitz 判據(jù)等方法證明了兩邊部平衡點產(chǎn)生影響.將兩時滯作為作為分支參數(shù)分析了不同情形下內(nèi)部平衡點Hopf 分支的存在情況，結(jié)論表明時滯可能會使正平衡點擾動，進(jìn)而產(chǎn)生周期解.由于沒有臨床醫(yī)學(xué)數(shù)據(jù)，沒有具體研究病毒復(fù)制的抑制率對免疫細(xì)胞濃度的影響.事實上，未來還有更復(fù)雜的模型值得考慮，比如考慮細(xì)胞的Logistic 增長，空間擴(kuò)散等.