葉 蔥，魏春金，張樹文

(集美大學(xué)理學(xué)院，福建 廈門 361021)

0 引言

長期以來，在種群動(dòng)力系統(tǒng)中，捕食與食餌模型是許多學(xué)者一個(gè)重要的研究方向，在此方面也取得了大量豐富的成果[1-5]。文獻(xiàn)[6]提出了模型

(1)

其中：x(t)、y(t)分別表示食餌和捕食者在t時(shí)刻的種群密度；a1>0、a2>0分別表示食餌和捕食者的內(nèi)稟增長率；b>0表示食餌種群的密度制約系數(shù)；f1y(t)/(c1+nx(t)+x2(t))為Holling IV功能反應(yīng)函數(shù)；f2y(t)/(c2+x(t))表示具有Leslies形式的捕食者數(shù)量反應(yīng)。近年來，學(xué)者們開始關(guān)注收獲對捕食-食餌系統(tǒng)動(dòng)力學(xué)的影響，并完成了大量的工作。在捕食-食餌系統(tǒng)的動(dòng)力學(xué)中，收獲函數(shù)起到了非常關(guān)鍵的作用。文獻(xiàn)[7]提出了常數(shù)收獲、比例收獲、非線性收獲3種類型的收獲功能，其中非線性收獲比常數(shù)收獲和比例收獲更具有現(xiàn)實(shí)性。然而，在現(xiàn)實(shí)世界中，隨機(jī)干擾是處處存在的，種群在環(huán)境中會受到各種隨機(jī)干擾的影響[8-10]。因此，與確定性模型相比，隨機(jī)模型更具有真實(shí)性。受文獻(xiàn)[11]的啟發(fā)，本文在模型(1)的基礎(chǔ)上獲得隨機(jī)模型

(2)

另一方面，周期現(xiàn)象在實(shí)際的生態(tài)系統(tǒng)中是普遍存在的，因?yàn)榉N群中的個(gè)體受生命周期、季節(jié)變化以及晝夜更替等影響，該物種的出生率、死亡率和其他參數(shù)都呈現(xiàn)出周期變化[12]。因此，考慮周期因素對相應(yīng)的動(dòng)力學(xué)系統(tǒng)行為的影響顯得尤為重要。所以，相應(yīng)系統(tǒng)(2)的非自治系統(tǒng)為

(3)

其中：a1(t)、b(t)、f1(t)、c1(t)、n(t)、p1(t)、m1(t)、σ1(t)、a2(t)、f2(t)、c2(t)、p2(t)、m2(t)、σ2(t)均為正的、有界的、連續(xù)的正w周期函數(shù)。

1 預(yù)備知識

隨機(jī)微積分方程

dx(t)=f(x(t)，t)dt+g(x(t)，t)dB(t)

(4)

的解記作x(t)=(x1(t)，x2(t)，…，xn(t))(t≥0)，其中f∈L1(Rn×R+，Rn)，g∈L2(Rn×R+，Rn×m)，B(t)是n維布朗運(yùn)動(dòng)。

定理1[13](存在唯一性定理) 假設(shè)f(x(t)，t)和g(x(t)，t)關(guān)于x(t)滿足下列條件：1)局部Lipschitz條件，存在ck>0(k=1，2，…)，使得?x，y∈Rn且|x|∨|y|≤k，有不等式|f(x，t)-f(y，t)|∨|g(x，t)-g(y，t)|≤ck|x-y|成立；2)線性增長條件，存在c>0，使得|f(x，t)|∨|g(x，t)|≤c|1+|x||，?(x，t)∈Rn×R+，則初始條件為x(0)=x0∈Rn的系統(tǒng)(4)存在唯一連續(xù)的局部解x(t)(t∈[0，τe))，τe是爆破時(shí)間。

定理2[13](It公式) 設(shè)x(t)(t≥0)是It過程，其隨機(jī)微分為dx(t)=f(t)dt+g(t)dB(t)，其中，f∈L1(R+，Rn)，g∈L2(R+，Rn×m)。若V(x(t)，t)∈C2，1(Rn×R+；R)，則V(x(t)，t)仍然是It過程，具有隨機(jī)微分dV(x(t)，t)=Vt(x(t)，t)dt+Vx(x(t)，t)dx(t)+dxT(t)Vxx(x(t)，t)dx(t)/2。

為了驗(yàn)證性質(zhì)2成立，只要證明存在非負(fù)的C2-函數(shù)及鄰域U，使得對于任意的X∈ElU，LV(X)是負(fù)的。

下面考慮方程

(5)

假設(shè)方程(5)的系數(shù)b(s，X(s))、σr(s，X(s))(r=1，2，…，k)滿足條件

(6)

(7)

其中：B是一個(gè)常數(shù)。

2 主要結(jié)果及其證明

證明首先考慮方程

(8)

令k0>0足夠大，使得(x(t)，y(t))∈[1/k0，k0]×[1/k0，k0]，對于任意的正數(shù)k≥k0，定義一個(gè)停時(shí)序列τk=inf{t∈[0，τe)：x(t)(1/k，k)或y(t)(1/k，k)}。

通過不等式y(tǒng)≤2(y-1-lny)+ln 4≤2V(x，y)+ln 4，可得LV(x，y)≤N1+(a2+f2/c2+f1/c1)ln 4+2(a2+f2/c2+f1/c1)V(x，y)≤N2(1+V(x，y))，其中，N2=max{N1+(a2+f2/c2+f1/c1)ln 4，2(a2+f2/c2+f1/c1)}，所以有

dV(x，y)≤N2(1+V(x，y))dt+(x-1)σ1x/(1+m1x)dB1+(y-1)σ2x/(1+m2x)dB2。

(9)

對式(9)從0到τk∧T1進(jìn)行積分可得

通過Gronwall不等式，得EV[x(τk∧T1)，y(τk∧T1)]≤(V[x(0)，y(0)]+N2T)eN2T。

后續(xù)證明與文獻(xiàn)[19]類似，在此省略。

證明首先考慮方程

接下來，考慮食餌種群x(t)。應(yīng)用It公式可得

定義一個(gè)Lyapunov函數(shù)lnx，通過It公式得，又因?yàn)閒1y/(c1+nx+x2)=y/(c1/f1+n/f1+x2/f1)。比較c1/f1+n/f1+x2/f1與c2+x的大小，c1/f1+n/f1+x2/f1-(c2+x)=(x2+(n-f1)x+c1-f1c2)/f1，所以Δ=(n-f1)2-4(c1-f1c2)。

(10)

證明i)定義Lyapunov函數(shù)lnx，由It公式得，對上述不等式兩邊從0到t進(jìn)行積分，并同時(shí)除以t，得即因此，對于任意小的ε>0，存在t0和一個(gè)集合Ωε，使得對t>t0和ω∈Ωε，都有P(Ωε)≥1-ε和x(t)/(c2+x(t))≤ε成立。由等式dy=y[a2-f2y/(c2+x)-p2/(1+m2y)]dt-σ2y/(1+m2y)dB2(t)=y[a2-f2y/c2+f2yx/((c2+x)c2)-p2/(1+m2y)]dt-σ2/(1+m2y)dB2(t)，可得y[a2-f2y/(c2+x)-p2/(1+m2y)]dt-σ2y/(1+m2y)dB2(t)≤dy≤y[a2-f2y/c2+εf2y/c2-p2/(1+m2y)]dt-σ2y/(1+m2y)dB2(t)。

對任意的ε>0，有

dy=y[a2-f2y/(c2+x)-p2/(1+m2y)]dt-σ2y/(1+m2y)dB2(t)。

(11)

(12)

(13)

(14)

(15)

3 數(shù)值模擬結(jié)果

對于隨機(jī)模型(2)，首先取初值x=0.1，y=0.2。選擇適當(dāng)參數(shù)：a1=0.4，b=0.1，m1=1/6，m2=1/6，f1=0.1，a2=0.4，f2=1.1，c2=5，c1=4，n=1，p1=0.6，p2=0.6，σ1=0.2，σ2=0.2。通過計(jì)算，顯然滿足定理5的條件ii)。由圖3可以看出，食餌x與捕食者y均滅絕，說明過度捕獲將導(dǎo)致種群減少。

對于隨機(jī)模型(3)，首先取初值x=0.1，y=0.2。選擇適當(dāng)參數(shù)：a1=0.4+0.02sin(πt/5)，b=0.1+0.02sin(πt/5)，m1=1+0.02sin(πt/5)，m2=1+0.02sin(πt/5)，f1=0.1+0.02sin(πt/5)，a2=0.4+0.02sin(πt/5)，f2=1.1+0.02sin(πt/5)，c2=5+0.02sin(πt/5)，c1=4+0.02sin(πt/5)，n=1+0.02sin(πt/5)，p1=0.06+0.02sin(πt/5)，p2=0.06+0.02sin(πt/5)，σ1=0.01+0.001sin(πt/5)，σ2=0.01+0.001sin(πt/5)，顯然符合定理7的ξ1>0、ξ2>0條件。由定理7可以看出，系統(tǒng)(3)至少存在一個(gè)周期解，這意味著捕食者和食餌種群將長期共存并表現(xiàn)出周期性。由圖4和圖5可以看出，對于任何正的初始值，確定性模型的解將在一段時(shí)間后進(jìn)入周期軌道，而當(dāng)噪聲強(qiáng)度較小時(shí)，隨機(jī)模型(3)的解會在周期軌道的小范圍內(nèi)振動(dòng)。

模擬結(jié)果說明，環(huán)境噪聲強(qiáng)度和可捕獲性對隨機(jī)捕食-食餌模型有重要的影響。

4 結(jié)論

本文研究了具有Holling IV功能性反應(yīng)和非線性收獲的隨機(jī)捕食-食餌系統(tǒng)。證明了對于任意給定的初始值，系統(tǒng)(2)都存在唯一的全局正解；得到了系統(tǒng)(2)的平均持續(xù)生存與滅絕的充分條件。此外，在滿足一定條件下，使用Has’minskii的平穩(wěn)分布理論及周期性理論，得到了系統(tǒng)(2)存在唯一的平穩(wěn)分布且具有遍歷性，并獲得了系統(tǒng)(3)存在非平凡的正周期解。