王藝霖,王 新,劉 兵,侯 祥,董秀輝

(1.遼寧師范大學(xué) 數(shù)學(xué)學(xué)院,遼寧 大連 116029;2.鞍山師范學(xué)院 數(shù)學(xué)與信息科學(xué)學(xué)院,遼寧 鞍山 114007;3.鞍山市氣象局,遼寧 鞍山 114004)

水母是水生環(huán)境中重要的浮游生物,有著復(fù)雜的結(jié)構(gòu)和生命史,是海洋生態(tài)系統(tǒng)中的一個(gè)重要物種.水母暴發(fā)是指水母在特定季節(jié)、特定海域內(nèi)數(shù)量激增的現(xiàn)象.以往大約40年才有一次水母暴發(fā),近年來越來越頻繁.研究發(fā)現(xiàn),水母數(shù)量超過一定水平,就會(huì)對旅游業(yè)、發(fā)電廠和核電站的正常工作造成影響[1]；水母與魚類爭奪餌料也會(huì)使各大漁場獲魚量持續(xù)下降[2].為此，要對水母進(jìn)行控制,把這種需要控制的危害水平稱為經(jīng)濟(jì)閾值.

目前，中外學(xué)者針對水母暴發(fā)問題已做了一些研究,包括缽水母生活史的長期動(dòng)力學(xué)行為、建立不同生命周期階段受溫度影響的過程模型，等等[3-4],但關(guān)于控制水母方面的研究成果卻很少.本文利用具有切換策略的動(dòng)力學(xué)模型模擬水母的控制過程,當(dāng)水母數(shù)量未超過經(jīng)濟(jì)閾值時(shí),不加干預(yù);超過經(jīng)濟(jì)閾值時(shí)，再對水母進(jìn)行捕撈,將水母數(shù)量控制在經(jīng)濟(jì)閾值內(nèi),為此建立Filippov水母控制模型.關(guān)于Filippov模型,在害蟲治理[5]、植物疾病研究[6]、傳染病研究[7]等多個(gè)方面已有廣泛應(yīng)用.本文得到的結(jié)果希望能為沿海生態(tài)和漁業(yè)、旅游業(yè)等發(fā)展提供相應(yīng)的理論依據(jù).

1 模型的建立和預(yù)備知識

1.1 模型的建立

假設(shè)水母種群是密度制約，而水螅體種群為非密度制約.水螅體通過自我復(fù)制和水母的有性繁殖補(bǔ)充,而水母只能通過水螅體的橫裂補(bǔ)充.當(dāng)水母種群數(shù)量小于經(jīng)濟(jì)閾值ET時(shí),建立如下方程：

(1)

其中，x(t),y(t)分別表示t時(shí)刻水螅體和水母的數(shù)量,a為實(shí)數(shù),b,c,d,b2為非負(fù)實(shí)數(shù).a表示水螅體的自我復(fù)制率與死亡率的差,假設(shè)a始終小于0;b表示水母有性繁殖成水螅體的存活率;c為水螅體無性繁殖成水母的存活率;d為水母的死亡率;b2為水母種內(nèi)競爭率.

當(dāng)水母種群數(shù)量大于ET時(shí),按比例δ(為非負(fù)實(shí)數(shù))打撈水母,此時(shí)方程為

(2)

系統(tǒng)(1)、(2)等價(jià)于

(3)

其中,

1.2 Filippov系統(tǒng)預(yù)備知識

令

H(Z)=y(t)-ET,

其中,

Z=(x,y)T,

且

FI1(Z)=(ax+by,cx-dy-b2y2)T,

FI2(Z)=(ax+by,cx-dy-b2y2-δy)T,

那么系統(tǒng)(1)和系統(tǒng)(2)合寫為Filippov系統(tǒng)

(4)

其中,

另外,記

為分割兩個(gè)區(qū)域I1和I2的分界線.

定義1令

Σs={Z∈Σ|〈HZ,FI1(Z)〉〈HZ,FI2(Z)〉<0},

稱Σs為滑線區(qū)域,其中〈·〉表示內(nèi)積,H(Z)是一個(gè)光滑的純量函數(shù),在Σ上關(guān)于H(Z)的梯度為

HZ=(0,1).

引理1(Filippov凸理論)如果滑線是光滑的,滑線系統(tǒng)可以表示為

其中,

(ii) 如果系統(tǒng)(4)的滑線區(qū)域Σs的平衡點(diǎn)Z1滿足

λFI1(Z1)+(1-λ)FI2(Z1)=0,0<λ<1,

則稱Z1為系統(tǒng)(4)的偽平衡點(diǎn).

2 子系統(tǒng)動(dòng)力學(xué)

2.1 系統(tǒng)(1)的動(dòng)力學(xué)性質(zhì)

系統(tǒng)(1)始終存在滅絕平衡點(diǎn)A0(0,0),當(dāng)ad+bc>0時(shí),存在正平衡點(diǎn)A1(x1,y1),其中，

定理1對于系統(tǒng)(1),當(dāng)ad+bc<0時(shí),正平衡點(diǎn)不存在,A0為全局漸近穩(wěn)定的;當(dāng)ad+bc>0時(shí),A0為鞍點(diǎn),系統(tǒng)(1)存在唯一的正平衡點(diǎn)A1,且它是全局漸近穩(wěn)定的.

證明系統(tǒng)(1)在平衡點(diǎn)A(x,y)處的Jacobian行列式為

Tr(J(A0))=a-d,Det(J(A0))=-(ad+bc).

當(dāng)

ad+bc<0

時(shí),容易看出

Tr(J(A0))<0,Det(J(A0))>0,

此時(shí)A0為局部漸近穩(wěn)定的.

當(dāng)

ad+bc>0

時(shí),顯然此時(shí)A0為鞍點(diǎn);此時(shí)存在唯一正平衡點(diǎn)A1,

Det(J(A1))=ad+bc.

可以得出

Tr(J(A1))<0,Det(J(A1))>0,

此時(shí)A1為局部漸近穩(wěn)定的.又知不存在極限環(huán),因此A1為全局漸近穩(wěn)定的.

2.2 系統(tǒng)(2)的動(dòng)力學(xué)性質(zhì)

類似系統(tǒng)(1)的分析,顯然系統(tǒng)(2)始終存在滅絕平衡點(diǎn)B0(0,0),當(dāng)

ad+bc+aδ>0

時(shí),存在正平衡點(diǎn)B1(x2,y2),其中，

類似于定理1的證明，可得如下結(jié)論：

定理2對于系統(tǒng)(2),當(dāng)ad+bc+aδ<0時(shí),正平衡點(diǎn)不存在,B0為全局漸近穩(wěn)定的;當(dāng)ad+bc+aδ>0時(shí),系統(tǒng)(2)存在唯一的正平衡點(diǎn)B1,且它是全局漸近穩(wěn)定的.

3 滑動(dòng)模態(tài)及平衡點(diǎn)

由定義1知,滑面Σs存在條件為F2>0,F4<0,得

Σs={(x,y)|y=ET,x′≤x≤x″},

其中,

由引理1,滑面上的動(dòng)力系統(tǒng)為

當(dāng)

時(shí),滑面存在唯一偽平衡點(diǎn)

容易看出該條件等價(jià)于

y2

令

F(x)=ax+bET,F′(xEp)=a<0,

故偽平衡點(diǎn)存在時(shí),為局部漸近穩(wěn)定的.

4 全局動(dòng)力學(xué)

類似于文獻(xiàn)[7],有如下結(jié)論：

定理3Filippov系統(tǒng)(4)中不存在以下幾種極限環(huán)：

(1)包含于區(qū)域Ii(i=1,2)內(nèi)部的極限環(huán)(如圖1(a)所示);

(2)環(huán)繞滑線的穿越極限環(huán)(如圖1(b)所示);

(3)不環(huán)繞滑線的穿越極限環(huán)(如圖1(c)所示);

(4)鴨型環(huán)(如圖1(d)所示).

圖1為Filippov系統(tǒng)(4)用以表示吸引域內(nèi)各種可能的不同類型極限環(huán)的相圖,x′x″表示滑線,細(xì)虛線表示橫截區(qū)域,Li(i=1,2)表示可能的極限環(huán).

令

A1=(x1,y1),B1=(x2,y2),

根據(jù)正平衡點(diǎn)A1、B1,顯然有y2

由以上分析,可以得到：

圖1 各極限環(huán)相圖

5 結(jié)論

本文基于水螅體和水母的關(guān)系建立了一個(gè)Filippov水母控制模型,通過理論推理發(fā)現(xiàn)ET取值不同,對應(yīng)的動(dòng)力學(xué)性質(zhì)也不同.由定理4可知:

(1)當(dāng)y2

(2)當(dāng)y2

(3)當(dāng)y1>y2>ET時(shí),B1為真平衡點(diǎn)且全局漸近穩(wěn)定,此時(shí)水母數(shù)量超過經(jīng)濟(jì)閾值,對生態(tài)造成危害,不是想要的結(jié)果,在實(shí)際捕撈水母中,要盡量避免這種情況發(fā)生(見圖2(c)).

圖2 平衡點(diǎn)穩(wěn)定性相圖

a=-0.298,b=0.4,c=0.02,d=0.01,b2=0.05,δ=0.01.