馬亞妮,袁海龍,2,王雅迪

(1.陜西科技大學(xué) 數(shù)學(xué)與數(shù)據(jù)科學(xué)學(xué)院,陜西 西安 710021;2.西安交通大學(xué) 數(shù)學(xué)與統(tǒng)計(jì)學(xué)院,陜西 西安 710049)

0 引 言

Sel’kov模型是一個(gè)典型的反應(yīng)擴(kuò)散系統(tǒng)[1],它最初是作為糖酵解模型提出,如下所示:

式中:Ω是RN(N>1)中具有光滑邊界?Ω的有界開區(qū)域,n是?Ω上的單位外法向量,u(x,t)和v(x,t)代表時(shí)間t>0和空間位置x∈Ω時(shí)的2種反應(yīng)物的濃度或2種物質(zhì)的密度;參數(shù)θ是物種的擴(kuò)散系數(shù),且θ、p和λ都是正常數(shù)。該反應(yīng)擴(kuò)散模型(1)已被用于描述各種形式的地貌形態(tài)[2]、種群動(dòng)力學(xué)[3]和自催化氧化反應(yīng)[4]的研究。更多模型的背景見文獻(xiàn)[5-8]。

近幾十年來,模型(1)的動(dòng)力學(xué)行為已經(jīng)被許多生物數(shù)學(xué)家所研究。特別地,WANG通過固定p和λ(或p和θ),以θ(或λ)為分歧參數(shù),研究了系統(tǒng)(1)的非常數(shù)正穩(wěn)態(tài),證明了其平衡態(tài)方程非常數(shù)正解的存在性和不存在性的條件[8]。此后,PENG討論了當(dāng)θ取很大值時(shí),若01,則會(huì)出現(xiàn)非常數(shù)的正穩(wěn)態(tài),表明參數(shù)p在兩反應(yīng)物的空間非均勻分布中起著關(guān)鍵作用[9]。同時(shí),PENG等也討論了在Dirichlet邊界條件下,利用不動(dòng)點(diǎn)指數(shù)理論,給出了非負(fù)穩(wěn)態(tài)解的存在性結(jié)果[10]。進(jìn)一步,HAN等以參數(shù)p為分歧參數(shù),證明在合適的條件下,該系統(tǒng)存在Hopf分歧解,同時(shí)建立了該周期解的穩(wěn)定性和分歧方向[11],詳見文獻(xiàn)[12-14]。

由于種群的成長(zhǎng)期,當(dāng)前時(shí)刻的種群增長(zhǎng)率總是與某個(gè)時(shí)刻之前的種群數(shù)量有關(guān)。這種由成長(zhǎng)期而導(dǎo)致的時(shí)滯現(xiàn)象在種群中普遍存在。因此,近年來,在種群生態(tài)學(xué)和生物學(xué)模型中提出了具有時(shí)滯的反應(yīng)擴(kuò)散系統(tǒng)。李冬梅等研究了一類含分布時(shí)滯和Michaelis-Menten型的擴(kuò)散模型,利用微分方程比較定理和Lyapunov函數(shù)方法,得到了一致持久和全局漸近穩(wěn)定的充分條件[15]。王長(zhǎng)有等利用上、下解方法及不動(dòng)點(diǎn)理論,研究了一類反應(yīng)項(xiàng)非單調(diào)的時(shí)滯反應(yīng)擴(kuò)散方程組,獲得了邊值問題周期解存在性的充分條件[16]。劉高楊等基于KADDAR和ALAOUT的模型,建立了一個(gè)帶有時(shí)滯的腫瘤免疫擴(kuò)散模型,通過分析免疫細(xì)胞識(shí)別腫瘤細(xì)胞時(shí)間的閾值,可判斷腫瘤的最佳治療時(shí)間,減少患者因腫瘤導(dǎo)致的死亡率[17]。李振振等研究一類具有 Dirichlet邊界條件的時(shí)滯合作-擴(kuò)散-平流系統(tǒng),證明了時(shí)滯的增加會(huì)使正穩(wěn)態(tài)解失穩(wěn),并且當(dāng)時(shí)滯值穿過臨界分支點(diǎn)時(shí)系統(tǒng)會(huì)存在Hopf分歧[18]?；谝陨夏Ｐ偷姆治?將時(shí)滯引入到該系統(tǒng)中。本文將研究時(shí)滯對(duì)u的影響,則系統(tǒng)(1)可寫為

式中:τ>0代表延遲效應(yīng)。

1 平衡點(diǎn)的穩(wěn)定性與分歧的存在性

在有界開區(qū)域Ω=(0,lπ),l∈R+上考慮系統(tǒng),

式中:u=u(x,t);uτ=u(x,t-τ);v=v(x,t)。定義X=C([0,lπ],R2),在抽象空間C([0,lπ],X)中,系統(tǒng)(4)可以寫成以下抽象的泛函微分方程

(5)

這里,dΔ=(d1Δ,d2Δ)

dom (dΔ)={(u,v)T:u,v∈C2([0,lπ],R),

ux,vx=0,x=0,lπ}

且L:C([-τ,0],X)→X,F:C([-τ,0],X)→X。

根據(jù)系統(tǒng)(4)有

L(φ)=[-λφ1(0)-λpφ2(0)

λφ1(-τ)+λ(p-1)φ2(0)]

F(φ)=F1(φ)F2(φ)

式中:φ=(φ1,φ2)T∈C([-τ,0],X);

F1(φ)=λ-λ[φ1(0)+1][φ2(0)+1]p+

λφ1(0)+λpφ2(0);

F2(φ)=λ[(φ1(-τ)+1)(φ2(0)+1)p]-

λ-λφ1(-τ)-λpφ2(0)。

因此系統(tǒng)(4)在平衡點(diǎn)(1,1)處的線性化可寫為

(6)

其對(duì)應(yīng)的特征方程為

λy-dΔy-L[exp(λ)y]=0,

y∈dom(dΔ),y≠0

(7)

根據(jù)特征值問題

-ψ″=μψ,x∈(0,lπ),ψ′(0)=ψ″(lπ)=0

根據(jù)特征方程(7)知

λ2+Anλ+Bn+Cexp(-λτ)=0,n=0,1,2,…

(8)

C=λ2p。

若±iσ(σ>0)是式(8)的一對(duì)純虛根,則

由上述計(jì)算化簡(jiǎn)可得

式中:

對(duì)于0≤n≤N0(p),式(10)有一個(gè)正根σn滿足

那么當(dāng)

證明 將λn(τ)代入式(8)并對(duì)τ求導(dǎo)有

因此,根據(jù)以上引理1～3,可以得到如下結(jié)論。

2 Hopf分歧的穩(wěn)定性和方向

(18)

其中,對(duì)于φ∈C([-1,0],X)有

G(φ,μ)=μdΔφ(0)+μL0(φ)+

(μ+τ0)F0(φ)F0(φ)=

由第二部分可知,±iσ0τ0是線性系統(tǒng)的一對(duì)純虛特征值。因此有

(19)

其線性泛函微分方程為

(20)

應(yīng)用Riesz表示定理可知,存在一個(gè)有界變差2×2的矩陣函數(shù)η(θ,μ)(-1≤θ≤0),使得

事實(shí)上,

η(θ,μ)=(τ0+μ)Eδ(θ)-(τ0+μ)Fδ(θ+1)

其中,

?φ(θ)∈C1([-1,0],R2),定義A(0)為

以及?ψ=(ψ1,ψ2)∈C1([-1,0],(R2)*),定義

定義雙線性函數(shù)

其中,A(0)和A*是正規(guī)伴隨算子。

易證±iσ0τ0是A(0)和A*的特征值

q(θ)=(q1,q2)Texp(iσ0τ0θ)(θ∈[-1,0])

和

分別是A(0)和A*對(duì)應(yīng)于特征值iσ0τ0和-iσ0τ0的特征向量,其中

iσ0τ0exp(-iσ0τ0)

這里,u=(u1,u2),v=(v1,v2)∈X=C([0,lπ],R2)。并且對(duì)φ∈C([-1,0],X),有

線性方程(19)在μ=0時(shí)的中心子空間是PCNC,且

PCNφ=Φ(Ψ,〈φ,f0〉)0·f0,φ∈C,

對(duì)C進(jìn)行空間分解,C=PCNC⊕PSC,其中,PSC表示穩(wěn)定子空間。

由文獻(xiàn)[25]可知,線性系統(tǒng)(19)無(wú)窮小生成元AU滿足

由于只考慮在μ=0附近的分歧周期解的穩(wěn)定性和方向,因此令系統(tǒng)(18)中的μ=0,可得中心流形

則系統(tǒng)(18)的解可以表示為

式中:

從而

記

定義

f(u,v)=λ(u+1)(v+1)p。

根據(jù)泰勒展開式,得

其中,O(4)=O(‖(u,v)‖4),又根據(jù)

G(φ,0)=τ0(G1,G2)T。

那么

從而由式(25)、(27)和(28)可以得到:

如上表達(dá)式可以看出,為了計(jì)算g21,需要計(jì)算W20(θ)和W11(θ)。

根據(jù)

可得

當(dāng)-1≤θ<0時(shí),

于是結(jié)合式(25)和(31),當(dāng)-1≤θ<0時(shí),有

那么聯(lián)立式(30)的第一個(gè)方程和式(32)可得

E1exp(2iσ0τ0θ)

同理再聯(lián)立式(30)的第二個(gè)方程和式(33)可得

令式(30)中θ=0,再利用AU的定義和

分別得到

E1=E11E12和E2=E21E22

其中,

綜上,g21可由計(jì)算得到。基于以上分析可以看出,每個(gè)gij可以通過參數(shù)來確定,于是

定理2 根據(jù)系統(tǒng)(3)有如下結(jié)論。

(ⅱ)β2決定分歧周期解的穩(wěn)定性:若β2<0,(β2>0),則分歧周期解是漸近穩(wěn)定的(不穩(wěn)定的)。

(ⅲ) T2決定分歧周期解的周期變化:若T2>0,(T2<0),則周期增大(減小)。

3 數(shù)值模擬

選取時(shí)滯τ作為分歧參數(shù),其中τ=0.11,初值u(x,t)=1+0.1tcos x,v(x,t)=1-0.011cos x,t∈[-0.11,0],x∈[0,3π]。得到模型(3)平衡點(diǎn)(1,1)的穩(wěn)定性,數(shù)值模擬結(jié)果見圖1。

(a) u(x,t)數(shù)值模擬

由圖1知,當(dāng)τ∈(0,0.339 3)時(shí),系統(tǒng)(3)的解趨向于平衡點(diǎn)(1,1),此時(shí)正平衡點(diǎn)(1,1)是穩(wěn)定的。

選取時(shí)滯τ作為分歧參數(shù),其中τ=0.41,初值u(x,t)=1+0.1cos(4x),v(x,t)=1-0.011cos(x2),t∈[-0.41,0],x∈[0,3π]。那么周期解存在,數(shù)值模擬結(jié)果見圖2。

(a) u(x,t)數(shù)值模擬

由圖2知,當(dāng)τ在0.339 3的一個(gè)很小的右鄰域內(nèi)時(shí),系統(tǒng)(3)在正常數(shù)平衡點(diǎn)(1,1)處產(chǎn)生分歧周期解。