李 良， 王會(huì)超

(1.四川大學(xué)數(shù)學(xué)學(xué)院， 成都 610064； 2.許昌學(xué)院數(shù)學(xué)與統(tǒng)計(jì)學(xué)院， 許昌 461000)

1 引 言

本文所研究的具有擴(kuò)散項(xiàng)抗病毒治療模型[1-7]是在文獻(xiàn)[1]中的基礎(chǔ)模型后添加擴(kuò)散項(xiàng)后得到的.原模型只考慮了均勻分布情形，添加了擴(kuò)散項(xiàng)的模型則更貼近實(shí)際.模型的定義如下：

(1)

其中u1表示未感染細(xì)胞數(shù)，u2表示感染細(xì)胞數(shù)，u3表示游離病毒數(shù)，未感染細(xì)胞以常熟率μ生成，以速率γu1死亡，βu1u3表示未感染細(xì)胞與游離病毒接觸之后變成感染細(xì)胞的速率，δu2表示感染細(xì)胞裂解死亡的速率，αu3表示游離病毒被清除的速率，Δu表示擴(kuò)散率，ν表示感染細(xì)胞回復(fù)至未感染細(xì)胞的比率，θu2表示感染細(xì)胞裂解后產(chǎn)生病毒速率.抗病毒藥物主要以干擾素和一拉米夫啶為代表，其中1-ε表示治療效果(0<ε<1)，1-λ表示核苷類藥物阻礙從DNA到CCCDNA循環(huán)轉(zhuǎn)移的效率.此外，Neumann邊界條件表示細(xì)胞和病毒與外界沒有進(jìn)行交換.這里我們只考慮x∈Ω=[0,l]和t>0的情形.

2 預(yù)備知識(shí)

針對(duì)非線性演化方程的抽象形式[2]

(2)

其中λ為參數(shù)，我們假設(shè)Lλ:H1→H為一個(gè)全連續(xù)算子，且為扇形算子，H1和H都為Hilbert空間(H1→H緊稠)， 并假設(shè)Lλ的特征值{βj(λ)∈C|j=1,2,…}滿足條件

(3)

由線性全連續(xù)場的譜理論可知,H1和H有分解

其中

算子Lλ有分解

另外，非線性項(xiàng)算子G(u,λ)在u=0處可展開成如下形式：

(4)

其中n≥2是一個(gè)整數(shù)，且Gn(u,λ)是一個(gè)n重線性算子.空間Hσ的定義參見文獻(xiàn)[3].

定理2.1[3]假設(shè)條件(3)成立，且G(u,λ)有形如(4)式的一階Taylor展開式.那么在λ0附近中心函數(shù)Φ(x,λ)可表示為如下形式：

O(|Reβ(λ)|‖x‖k)+o(‖x‖k),

另外，文中需要用到的Morse指數(shù)理論可參見文獻(xiàn)[3].

3 方程處理

(5)

我們將方程(5)化成形如方程(2)的形式.為此，首先建立空間如下：

H={(u1,u2,u3)∈L2(Ω,R3)}.

令Lλ=A+Bλ，其中

易知,Lλ:H1→H是全連續(xù)算子，并且為扇形算子，G(·,λ):H1→H滿足條件(4).程化為

現(xiàn)在我們來算方程(5)的特征根和特征向量.

首先，求拉普拉斯算子的特征值.由

可得

其中

易知

且在λ0附近有

再根據(jù)線性全連續(xù)場譜理論得到

其中

因此，任意u∈H1均可表達(dá)為下列形式：

4 結(jié) 論

其中

式中涉及的一些量的定義可參見文獻(xiàn)[3].

證明 令u∈H1.那么

P1ut=P1Lλu+P1G(u,λ),

P2ut=P2Lλu+P2G(u,λ)，

及

化簡可得分歧方程為

(6)

再令

則

其中

因此，

其中

這里0<|λ-λ0|<ε.因此方程(7)的穩(wěn)態(tài)解為

從而方程(5)的分歧解為

進(jìn)而方程(1)的分歧解為

o(|λ-λ0|).

結(jié)論(iii)得證.

當(dāng)λ0<λ<λ0+ε時(shí)，若x0(0)<0，則x0(t)<0，從而

這就證明了方程(1)在區(qū)域x0(0)<0是不連續(xù)躍遷的.

同樣方法可以證明

這就證明了方程(1)在區(qū)域x0(0)>0里是連續(xù)躍遷的.結(jié)論(i)得證.證畢.

注易求得方程(7)的解為

5 數(shù)值算例

當(dāng)λ=0.45時(shí)，

從圖1可見，Neumann邊界條件下的解(u1,u2,u3) →(0.4,0,0),與理論得出的結(jié)論一致.

再令λ=0.55，理論得出的奇點(diǎn)分歧解為(u1,u2,u3)→(0.33,0.44,0.41).

圖1 λ=0.45時(shí)模型(1)的平衡解Fig.1 The equilibrium solutions of Model (1) when λ=0.45

圖2表明，方程解(u1,u2,u3)→(0.36,0.17,0.38),此時(shí)理論解與數(shù)值解誤差較小.

圖2 λ=0.55時(shí)模型(1)的分歧解Fig.2 The birfurcation solutions of Model (1) when λ=0.55