申佩嫻，薛亞奎

(中北大學(xué) 理學(xué)院， 太原 030051)

根據(jù)實(shí)驗(yàn)結(jié)論[10]，只有當(dāng)食餌感知到捕食者散發(fā)出的氣味時(shí)才會(huì)減低其自身的繁殖率，造成種群整體數(shù)量的減少，進(jìn)而影響捕食者種群的數(shù)量。假設(shè)捕食者只能捕食避難所之外的食餌，并且捕食者散發(fā)出的氣味只對(duì)避難所之外的食餌的出生率造成影響，由此考慮避難所保護(hù)食餌的數(shù)量對(duì)模型的影響：當(dāng)避難所保護(hù)食餌的數(shù)量逐漸增多時(shí)，會(huì)導(dǎo)致捕食者可捕食的食餌數(shù)量減少。本文用r0[1-cy]+刻畫捕食者釋放的氣味對(duì)食餌出生率所造成的影響，其中把單個(gè)捕食者產(chǎn)生的氣味干擾看做常數(shù)c，捕食者數(shù)量用y(t)表示。當(dāng)氣味干擾足夠小或捕食者種群數(shù)量足夠少時(shí)，避難所外食餌的出生率將不受捕食者氣味影響；反之，食餌出生率將為0。由此建立如下模型：

(1)

其中：x(t)表示捕食種群數(shù)量；y(t)表示食餌種群數(shù)量；r0表示食餌的自然出生率；c表示氣味對(duì)食餌出生率的干擾系數(shù)；mx(t)表示避難所保護(hù)的食餌數(shù)量；(1-m)x(t)表示能被捕食者能探測(cè)到的食餌數(shù)量；a表示食餌種群的密度制約系數(shù)；p表示捕食者的捕獲率；u表示捕食者捕獲食餌的轉(zhuǎn)化率(u

假設(shè)m<1且c足夠小，也就是說(shuō)：避難所保護(hù)的食餌數(shù)量要少于食餌種群的總數(shù)量。由于氣味對(duì)食餌種群的干擾系數(shù)依賴于捕食者種群的數(shù)量，因此只有氣味干擾系數(shù)足夠小時(shí)，避難所外的食餌才可以繼續(xù)繁殖和生長(zhǎng)。如果不滿足以上條件,捕食者種群將會(huì)滅絕。因此，系統(tǒng)(1)變?yōu)椋?/p>

(2)

1 系統(tǒng)解的有界性

定理1 系統(tǒng)(2)的全部正解都是一致有界的。

證明定義一個(gè)關(guān)于解的和函數(shù)[11]z(t)=x(t)+y(t)，則

-ax2+r0x-[p(1-m)-u(1-m)]xy-d1x-d2y≤x(-ax+r0)-mind1,d2}z

(3)

所以總存在一個(gè)正數(shù)β1，使得式(3)變?yōu)?/p>

2 模型的分析

2.1 平衡點(diǎn)的存在性

證明顯然系統(tǒng)存在滅絕平衡點(diǎn)E0；若(T1)成立，即r0>d1，存在邊界平衡點(diǎn)E1(x1,0)。

令系統(tǒng)(2)右端為零，即

(4)

由方程組(4)的第2個(gè)方程可得

(5)

將式(5)代入方程組(4)的第1個(gè)方程,可得

(6)

2.2 平衡點(diǎn)的局部穩(wěn)定性

定理3 設(shè)：(T3)r0

1) 若(T3)成立，則E0局部漸近穩(wěn)定；若(T1)成立，則E0不穩(wěn)定。

2) 若(T1)、(T4)同時(shí)成立，則E1局部漸近穩(wěn)定；若(T1)、(T2)同時(shí)成立，則E1不穩(wěn)定；若(T5)成立，則E1是一個(gè)鞍結(jié)點(diǎn)。

3) 若(T1)、(T2)同時(shí)成立，則E2局部漸近穩(wěn)定。

證明系統(tǒng)(2)在滅絕平衡點(diǎn)E0處的Jacobian矩陣為:

系統(tǒng)在平衡點(diǎn)E0處的特征方程為：

λ2+(d2+d1-r0)λ+(d1-r0)d2=0

(7)

顯然，λ1=-d2是方程(7)的一個(gè)根，若(T3)成立，則方程(7)的所有根均為負(fù)，E0局部漸近穩(wěn)定；若(T1)成立，則方程(7)有一個(gè)正根，一個(gè)負(fù)根，故E0不穩(wěn)定。

隨著我國(guó)市場(chǎng)經(jīng)濟(jì)的發(fā)展，保險(xiǎn)商業(yè)化的趨勢(shì)不斷加強(qiáng)。很多地區(qū)也開始試水醫(yī)療保險(xiǎn)的商業(yè)化。如果醫(yī)療機(jī)構(gòu)違反了醫(yī)療服務(wù)合同，給患者造成了人生傷害，可以通過購(gòu)買商業(yè)保險(xiǎn)的方式分散風(fēng)險(xiǎn)?，F(xiàn)在的商業(yè)醫(yī)療保險(xiǎn)的狀況不容樂觀，因?yàn)榛颊咛岢龅母哳~損害賠償常常令保險(xiǎn)公司措手不及，使得很多商業(yè)保險(xiǎn)公司放棄了醫(yī)療保險(xiǎn)業(yè)務(wù)。但是如果我們可以把醫(yī)療服務(wù)合同的地位法定化，約定好雙方當(dāng)事人的權(quán)利和義務(wù)，這樣商業(yè)保險(xiǎn)公司理賠的時(shí)候才會(huì)有依據(jù)，商業(yè)醫(yī)療保險(xiǎn)也會(huì)慢慢發(fā)展起來(lái)。醫(yī)療保險(xiǎn)制度的構(gòu)建，是一個(gè)“三贏”的局面，可以分散醫(yī)療機(jī)構(gòu)的經(jīng)營(yíng)風(fēng)險(xiǎn)，擴(kuò)大商業(yè)保險(xiǎn)公司的業(yè)務(wù)范圍，更有利于保護(hù)患者的利益。

系統(tǒng)(2)在邊界平衡點(diǎn)E1處的Jacobian矩陣為：

則

故若(T1)、(T4)同時(shí)成立時(shí)，E1局部漸近穩(wěn)定；若(T1)、(T2)同時(shí)成立時(shí)，E1不穩(wěn)定。

(8)

系統(tǒng)(2)在正平衡點(diǎn)E2處的Jacobian矩陣為：

其中：

a12=-r0cx*+r0cmx*-p(1-m)x*

a21=u(1-m)y*

a22=u(1-m)x*-d2

則：

其中：

2.3 平衡點(diǎn)的全局穩(wěn)定性

[r0-d1+u(1-m)]x-r0cyx(1-m)-ax2-p(1-m)xy-d2

證明令

F1(x,y)=r0mx+r0(1-cy)(1-m)x-p(1-m)xy

-ax2-d1x

F2(x,y)=u(1-m)xy-d2y

3 討論

圖1 系統(tǒng)(2)隨避難所系數(shù)m及氣味干擾系數(shù)c變化的軌線

從圖1(d)和(e)可知：隨著避難所系數(shù)m的增大，系統(tǒng)趨于穩(wěn)定的時(shí)間延長(zhǎng)，食餌種群和捕食種群的數(shù)量明顯增加。從圖1(f)可以看出：當(dāng)避難所系數(shù)m足夠大時(shí)，食餌數(shù)量將會(huì)大量增加，而捕食種群將趨于滅絕。因此，在一定的捕食者氣味干擾下對(duì)食餌種群加入適當(dāng)?shù)谋茈y所，不僅可以增加食餌種群的數(shù)量，同時(shí)也有利于捕食種群數(shù)量的增加?？梢?，氣味干擾下對(duì)食餌種群加入避難所對(duì)捕食系統(tǒng)起著關(guān)鍵性的作用。

[1] CRESSMAN R.,GARAY J.A predator-prey refuge system:evolutionary stability in ecological systems[J].Theor.Popul.Biol,2009,76:248-257.

[2] HUANG Y,CHEN F,ZHONG L.Stability analysis of prey-predator model with Holling type-III response function incorporating a prey refuge[J].Appl Math.Comput.2006,182:672-683.

[3] MA Z H，WANG S F，LI W D,et al.The effect of prey refuge in a patchy predator-prey system[J].Mathematical Biosciences 2013,243:126-130.

[4] MUKHERJEE D.The effect of prey refuges on a three species food chain model[J].Differential Equations Dynam.Systems,2014,22:413-426.

[5] LV Y,YUAN R,PEI Y,et al.Two types of predator-prey models with harvesting:non-smooth and non-continuous[J].Comput Appl Math,2013,250(10):122-142.

[6] SURACI J P,CLINCHY M,DILL L M,et al.Fear of large carnivores causes a trophic cascade[J].Nature communications,2016,73 (5)：1-26.

[7] ZHU Y,WANG K.Existence and global attractivity of positive periodic solutions for a predator-prey model with modified Leslie-Gower Holling-type II schemes[J].Math.Anal.Appl.2011,384:400-408.

[8] PAL P J,SARWARDI S,SAHA T,et al.Mean square stability in a modified Leslie-Gower and holling-typeII predator-prey model[J].Appl.Math.Inform.,2011,29 (3/4):781-802.

[9] LAI X,LIU S,LIN R.Rich dynamical behaviours for predator-prey model with weak Allee effect[J].Appl.Anal.,2010,8(89):1271-1292.

[10] 王振龍,劉季科.艾鼬氣味對(duì)根田鼠社會(huì)行為和繁殖的影響 [J].動(dòng)物學(xué)報(bào),2002,48(1):20-26.

[11] GHOSH J,SAHOO B,PORIA S.Prey-predator dynamics with prey refuge providing additional food to predator[J].Chaos,Solitons and Fractals,2017,96:110-119.

[12] SHAN C,ZHU H.Bifurcations and complex dynamics of an SIR model with the impact of the number of hospital beds [J].Journal of Differential Equations,2014,257(5):1662-1688.

[13] 馬知恩,周義倉(cāng).常微分方程定性與穩(wěn)定性方法 [M].北京：科學(xué)出版社,2001.

[14] SAHABUDDIN S,PRASHANTA K M,SANTANU R.Analysis of a competitive prey-predator system with a prey refuge[J].BioSystems,2012,110：133- 148.

[15] DEBASIS M.The effect of refuge and immigration in a predator-prey system in the presence of a competitor for the prey[J].Nonlinear Analysis:Real World Applications 2016,31:277-287.

[16] WANG X,ZANETTE L,ZOU X.Modelling the fear effect in predator-prey interactions[J].Journal of mathematical biology,2016,73(5):1179-1204.