β擾動的食餌有病的隨機食餌與捕食者系統的漸近行為

劉振文， 鄭凱鴻

(吉林建筑科技學院 基礎科學部， 吉林 長春 130022)

0 引 言

傳染病的傳播發生在人與人，人與動物或動物與動物之間。對于疾病流行規律的把握是對疾病防控的關鍵，傳染病動力學研究很好地處理了這類問題[1-2]。種群動力學揭示了種群的內在發展規律，學者們對于傳染病動力學和種群動力學的研究已有相當多的成果，對于這兩者相結合的生態流行病模型的研究也受到越來越多的關注[3-8]。Yanni Xiao等[9]建立了如下生態流行病模型

(1)

考慮到自然界中生態系統會受到環境噪聲的干擾，導致系統中某些參數發生變化，文中修正文獻[8]提出的系統(1)。首先，令

得到確定性系統

(2)

對系統(2)中β進行擾動，令系統中的β為

則系統(2)化為

(3)

式中：B(t)----布朗運動;

σ2----環境白噪聲的強度,σ2>0。

1 系統(3)的隨機漸近行為

1.1 系統(3)的正解存在唯一性

證明 對t≥0，考慮系統

(4)

(5)

則

S(t)+I(t)≤

(6)

式中：D3=S(0)+I(0)。

若證明這個解是幾乎必然全局的，就等價于證明τe=∞幾乎必然成立。選擇足夠大的k0≥0，使得S(0),I(0)和Y(0)全部位于區間[1/k0,k0]內，對每一個整數k≥k0，定義

或

max{S(t),I(t),Y(t)≥k}}。

(7)

P{τ∞≤T}>ε,

(8)

則存在一個整數k1≥k0,使得

P{τ∞≤T}≥ε,

(9)

(10)

所以上述函數的非負性是顯然的。使用Itǒ公式，有

(11)

其中

(12)

首先, 令

滿足

其次, 取a′>0充分小和合適的b′>0,滿足

和

和

則有

(13)

因此

(14)

則推得

E[V(S(τkΛT),I(τkΛT),Y(τkΛT))]≤

(15)

當k≥k1時，令Ωk={τk≤T} 且由式(9), 有P(Ωk) ≥ε，注意到對每一個ω∈Ωk, 在S(τk,ω),I(τk,ω)和Y(τk,ω)中至少有一個達到k或1/k,因此

V′(S(τkΛT),I(τkΛT),Y(τkΛT))≥

(16)

E[1Ωk(w)V′(S(τkΛT),I(τkΛT),Y(τkΛT))]≥

εh′(k),

(17)

從定理可以得到

I(t)>0,Y(t)>0,S(t)+I(t)≤K},

(18)

是系統(3)的正不變集，從現在開始總是假設初始值(S(0),I(0),Y(0))∈Г。

1.2 系統(3)的漸近性質

1.2.1 確定性系統的平衡點E1=(K，0,0)

顯然系統(2)一定有有界平衡點E1=(K，0,0)。但E1=(K，0,0)不是系統(3)的平衡點。此時我們給出系統(3),解得收斂率。

定理2令(S(t),I(t),Y(t))為系統(3)滿足初始條件 (S(0),I(0),Y(0))的解。若如下條件成立

1)Kβ>c;

則有

和

幾乎必然成立。

證明 由隨機比較定理有，I≤X，這里X是方程

(19)

即

(20)

由Itǒ公式，則

(21)

有

log(X(t))=log(X(0))+

KσB(t)。

(22)

令

F(t)=logX(0)+KσB(t),

(23)

則

幾乎必然成立。所以

(24)

其中，選取足夠大的噪聲強度σ2,使得

即I(t)以指數幾乎必然趨于0。換句話說，已感染者依概率1死亡。同樣有

dlogY=(kPI-d)dt,

(25)

則

(26)

所以

(27)

即Y(t)以指數幾乎必然趨于0。換句話說，捕食者依概率1死亡。

(28)

這里σ2充分小，且

(29)

定義

(30)

由Itǒ公式和式(27)可得

(31)

式中：Q=σSIdB(t)。

其中

(32)

(33)

其中

(34)

(35)

其中

kPR2Y+d(R-1)Y≤

kPR2Y,

(36)

則

(37)

定義

(38)

則

σ(F-1)Q。

(39)

令

(40)

則由Itǒ公式

(41)

其中

+k(σ)=

k(σ)=

(42)

其中

(43)

(44)

(45)

因此

(46)

其中

(47)

下面定義

(48)

其中V3=Y，從式(9)知道

LV3≤kPR2Y,

(49)

所以有

(50)

最后，定義

V=MV1+V4。

(51)

(52)

這樣得到

(53)

其中

(54)

σ2充分小，且

dV=LVdt+σ[MR1I+R1SI+FR2SI+

(55)

將上式兩邊從0到t積分得

(56)

這是一個實值連續局部鞅，且M(0)=0，而

(57)

則由強大數定律有

(58)

再由式(55)有

(59)

進一步有

(60)

定理得證。

從式(32)可以看出,系統(3)的解與E2在時間均值下的不同之處在于只與白噪聲的強度有關,噪聲強度越弱，系統(3)的解越穩定，所以當σ充分小時，可以認為有一種近似穩定性。