一類捕食者染病的捕食者-食餌系統(tǒng)的隨機(jī)動(dòng)力學(xué)行為

馮 濤，孟新柱

(山東科技大學(xué) 數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院，山東 青島 266590)

一類捕食者染病的捕食者-食餌系統(tǒng)的隨機(jī)動(dòng)力學(xué)行為

馮 濤，孟新柱

(山東科技大學(xué) 數(shù)學(xué)與系統(tǒng)科學(xué)學(xué)院，山東 青島 266590)

考慮了一類具有Beddington-DeAngelis功能性反應(yīng)和Lévy跳的捕食者染病的捕食者-食餌系統(tǒng)的動(dòng)力學(xué)行為。利用Lyapunov方法和伊藤公式，本文討論了系統(tǒng)全局正解的存在唯一性；研究了隨機(jī)系統(tǒng)在其確定性模型的平衡點(diǎn)周圍的長(zhǎng)時(shí)間行為。研究結(jié)果表明，在一定條件下，隨機(jī)系統(tǒng)的解會(huì)在其確定性系統(tǒng)的平衡點(diǎn)周圍波動(dòng)，且波動(dòng)的幅度與隨機(jī)系統(tǒng)所受干擾的強(qiáng)度呈正相關(guān)。最后，本文運(yùn)用Matlab數(shù)值模擬對(duì)前述理論進(jìn)行了驗(yàn)證。

捕食者-食餌系統(tǒng)；Beddington-DeAngelis功能性反應(yīng)；漸近行為；隨機(jī)擾動(dòng)；Lévy跳

近年來，隨機(jī)干擾對(duì)生態(tài)系統(tǒng)的影響逐漸成為生物數(shù)學(xué)領(lǐng)域的研究熱點(diǎn)[1-6]。在生態(tài)系統(tǒng)中，捕食者和食餌之間的相互作用經(jīng)常受到一些環(huán)境因素的影響，例如海嘯、水源、地震、火山、疾病等。許多學(xué)者研究了傳染病對(duì)捕食者-食餌種群的影響[7-13]和環(huán)境干擾對(duì)捕食者-食餌種群的影響[14-17]。研究傳染病和環(huán)境噪聲對(duì)捕食者-食餌種群的影響具有重要的生物學(xué)意義。本文結(jié)合傳染病與環(huán)境噪聲對(duì)生物種群的影響，研究了一類帶有Beddington-DeAngelis功能性反應(yīng)和Lévy跳的捕食者染病的捕食者-食餌系統(tǒng)的動(dòng)力學(xué)行為。

1 數(shù)學(xué)模型

文獻(xiàn)[13]研究了如下具有Beddington-DeAngelis功能性反應(yīng)和捕食者染病的捕食者-食餌系統(tǒng)：

(1)

令

實(shí)際上，生態(tài)系統(tǒng)中的物種可能會(huì)遭受到一些較大的不確定因素的干擾，比如火山噴發(fā)、地震、海嘯和氣候驟變等，一般把這類較大的干擾因素稱為有色噪聲干擾，數(shù)學(xué)上使用Lévy跳表示。考慮這類有色噪聲干擾后，模型(1)變?yōu)椋?/p>

(2)

2 全局正解的存在唯一性

本文的研究對(duì)象為生物種群，由于負(fù)數(shù)解不存在生物學(xué)意義，因此本節(jié)首先證明模型(2)存在全局唯一正解。

假設(shè) 2.1 假定以下條件成立[15]：

(i) 1+γi(u)>0，

其中，i=1,2,3，這兩個(gè)假設(shè)意味著Lévy噪聲的強(qiáng)度不會(huì)無限大。

證明由引理2.1可知，只需要證明τe=∞。定義一個(gè)充分大的正常數(shù)k0，使得

定義停時(shí)

(3)

V(X,S,I)=X-1-lnX+S-1-lnS+I-1-lnI。

(4)

由伊藤公式得

(5)

其中

(6)

其中，H0>0是常數(shù)。

對(duì)方程(5)的兩端從0到τk∧N積分并取期望，得

(7)

(8)

結(jié)合方程(7)和方程(8)，得

其中1nk 是Ωk的指標(biāo)函數(shù)。令k→+∞，得

+∞≥V(X(0),S(0),I(0))+H0N≥+∞。

顯然與假設(shè)矛盾。故τ∞=∞。

3 模型(2)解的漸近性質(zhì)

其中

證明定義

由伊藤公式得

(9)

其中

(10)

對(duì)方程(9)兩端由0到t積分并取期望，得

(11)

對(duì)方程(11)兩端同時(shí)除以t，令t→∞并取上界，得

證畢。

推論 3.1當(dāng)σ1,γ1=0時(shí)，模型(2)等同于模型(1)。由定理3.1中(10)可知，此時(shí)

所以模型(1)的平衡點(diǎn)E1(K,0,0)是全局漸近穩(wěn)定的。

注 3.1由定理3.1可知，如果干擾強(qiáng)度σ1,γ1足夠小，模型(2)的解會(huì)在模型(1)的邊界平衡點(diǎn)E1(K,0,0)周圍震蕩，且震蕩的幅度與干擾的強(qiáng)度呈正相關(guān)。

其中

定義

由伊藤公式得

(12)

其中

同理可得

(13)

其中

同理

(14)

結(jié)合方程(12)～(14)，得

(15)

其中

(16)

對(duì)方程(15)兩端由0到t積分并取期望，得

(17)

對(duì)方程(17)兩端同時(shí)除以t，令t→∞并取上界，得

證明完畢。

推論3.2當(dāng)σi,γi=0(i=1,2)時(shí)，模型(2)等同于模型(1)。由定理3.2中方程(16)知，此時(shí)

其中

證明由于(X*,S*,I*)是系統(tǒng)(1)的正平衡點(diǎn)，所以

定義

由伊藤公式得

(18)

其中

同理可得

(19)

其中

同理

(20)

其中

結(jié)合方程(18)～(20)，得

(21)

其中

對(duì)方程(21)兩端由0到t積分并取期望，得

(22)

對(duì)方程(22)兩端同時(shí)除以t，令t→∞并取上界，得

證明完畢。

推論 3.3當(dāng)σi,γi=0(i=1,2,3)時(shí)，模型(2)等同于模型(1)。由定理3.3知，此時(shí)

所以模型(1)的正平衡點(diǎn)E3(X*,S*,I*)是全局漸近穩(wěn)定的。

注 3.3由定理3.3知，如果干擾強(qiáng)度σi,γi=0(i=1,2,3)足夠小，模型(2)的解會(huì)在模型(1)的正平衡點(diǎn)E3(X*,S*,I*)周圍震蕩，且震蕩的幅度與干擾的強(qiáng)度呈正相關(guān)。

4 數(shù)值仿真

利用歐拉法和Matlab2014b[13]，本文進(jìn)行了數(shù)值模擬以支持獲得的結(jié)果。數(shù)值仿真中，使用以下參數(shù)：

初始值X(0)=3,S(0)=2,I(0)=2，p=q=1，Ζ=(0,+∞)，λ(Ζ)=1，步長(zhǎng)Δt=0.001。

在圖1中，b1=0.6，b2=0.4，b3=0.2，a11=0.5，a12=0.6，a21=0.3，a22=0.1，a33=0.1，β=0.2。圖1(a)為確定性模型(1)的時(shí)間序列圖，(b)-(d)為隨機(jī)模型(2)的時(shí)間序列圖，干擾強(qiáng)度分別取:(b)σi=0.05,γi=0.1，(c)σi=0.05,γi=0.15，(d)σi=0.05,γi=0.2，其中i=1,2,3。

(a) 確定性模型(1)；(b)-(d)為隨機(jī)模型(2)，干擾強(qiáng)度分別取值:(b)σi=0.05,ri=0.1，(c)σi=0.05,ri=0.15，(d)σi=0.05,ri=0.2，其中i=1,2,3。

在圖2中，b1=0.6，b2=0.12，b3=0.5，a11=0.5，a12=0.6，a21=0.8，a22=0.1，a33=0.4，β=0.1。圖2(a)為確定性模型(1)的時(shí)間序列圖，(b)～(d)為隨機(jī)模型(2)的時(shí)間序列圖，干擾強(qiáng)度分別取：(b)σi=0.05,γi=0.1，(c)σi=0.05,γi=0.15，(d)σi=0.05,γi=0.2，其中i=1,2,3。

在圖3中，b1=1，b2=0.1，b3=0.1，a11=0.2，a12=0.6，a21=0.1，a22=0.1，a33=0.2，β=0.5。圖3(a)為確定性模型(1)的時(shí)間序列圖，(b)-(d)為隨機(jī)模型(2)的時(shí)間序列圖，干擾強(qiáng)度分別取：(b)σi=0.05,γi=0.1，(c)σi=0.05,γi=0.15，(d)σi=0.05,γi=0.2，其中i=1,2,3。

在這種情況下，a11=0.2>p(b1-a11X*)=0.062 24，平衡E3(X*,S*,I*)=(4.69,0.66,1.15)R1=6.72>1。圖3顯示系統(tǒng)(2)的解在系統(tǒng)(1)的平衡點(diǎn)E3周圍浮動(dòng)，且浮動(dòng)的幅度與干擾的強(qiáng)度σi,γi(i=1，2，3)的取值呈正相關(guān)關(guān)系。這與定理3.3中的結(jié)論一致。

(a) 確定性模型(1)；(b)-(d)為隨機(jī)模型(2)，干擾強(qiáng)度分別取值:(b)σi=0.05,ri=0.1，(c)σi=0.05,ri=0.15，(d)σi=0.05,ri=0.2，其中i=1,2,3。

(a) 確定性模型(1)；(b)-(d)為隨機(jī)模型(2)，干擾強(qiáng)度分別取值:(b)σi=0.05,ri=0.1，(c)σi=0.05,ri=0.15，(d)σi=0.05,ri=0.2，其中i=1,2,3。

5 總結(jié)

本文研究了一類帶有Beddington-DeAngelis功能性反應(yīng)和Lévy跳的捕食者染病的捕食者-食餌系統(tǒng)。運(yùn)用李雅普諾夫方法和推廣的伊藤公式，本文首先證明了系統(tǒng)(2)全局正解的存在唯一性，然后討論了模型(2)的解在其確定性模型的平衡點(diǎn)周圍的漸近行為。Lévy跳在生物學(xué)上表示一些大的環(huán)境干擾，比如地震、火山、海嘯等。當(dāng)這類噪聲干擾發(fā)生的時(shí)候，會(huì)對(duì)生物種群的穩(wěn)定性產(chǎn)生一定的影響。文中干擾強(qiáng)度參數(shù)分別表示三個(gè)種群所受自然環(huán)境干擾的強(qiáng)度，比如地震、火山、暴雨等自然災(zāi)害的強(qiáng)度。定理3.1～3.3的理論結(jié)果表明，當(dāng)這類噪聲干擾不是特別大的時(shí)候，并不會(huì)使物種滅絕，但會(huì)使種群的數(shù)量在平衡狀態(tài)周邊浮動(dòng)，且這種浮動(dòng)的強(qiáng)度與噪聲干擾的強(qiáng)度呈正相關(guān)關(guān)系。也就是說，種群對(duì)環(huán)境的變化具有一定的適應(yīng)能力，當(dāng)環(huán)境變化的強(qiáng)度很小的時(shí)候，它們并不會(huì)導(dǎo)致物種的滅亡，而是使物種密度在一定范圍內(nèi)波動(dòng)。并且，當(dāng)這種干擾強(qiáng)度不斷增大時(shí)，它們對(duì)物種密度產(chǎn)生的影響也會(huì)越大。

[1]MENG X Z,ZHAO S N,FENG T,et al.Dynamics of a novel nonlinear stochastic SIS epidemic model with double epidemic hypothesis[J].Journal of Mathematical Analysis and Applications,2015,433:227-242.

[3]YAGI A,TON T V.Dynamic of a stochastic predator-prey population[J].Applied mathematics and Computation,2011,218(7):3100-3109.

[4]劉文昌,孟新柱.具有脈沖毒素輸入的隨機(jī)收獲模型最優(yōu)捕獲策略[J].山東科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2015,34(5):98-103. LIU Wenchang,MENG Xinzhu.Optimal harvesting strategies for stochastic harvest model with impulsive toxicant input[J].Journal of Shandong University of Science and Technology(Natural Science),2015,34(5):98-103.

[5]ZHANG Q M,JIANG D Q,LIU Z W,et al.Asymptotic behavior of a three species eco-epidemiological model perturbed by white noise[J].Journal of Mathematical Analysis and Applications,2016,433:121-148.

[6]LIU M,WANG K.Persistence,extinction and global asymptotical stability of a non-autonomous predator-prey model with random perturbation[J].Applied Mathematical Modelling,2012,36:5344-5353.

[7]HAQUE M,JIN Z,EZIO V.An ecoepidemiological predator-prey model with standard disease incidence [J].Mathematical Methods in the Applied Sciences,2009,32:875-898.

[8]CHATTOPADHYAY J,ARINO O.A predator-prey model with disease in the prey[J].Nonlinear Analysis-Series A Theory and Methods and Series B Real World Applications,1999,36:747-766.

[9]XIAO Y N,CHEN L S.Modeling and analysis of a predator-prey model with disease in the prey[J].Mathematical Biosciences,2001,171:59-82.

[10]HADELER K P,FREEDMAN H I.Predator-prey populations with parasitic infection[J].Journal of mathematical biology,1989,27:609-631.

[11]HAN L T,MA Z E,HETHCOTE H W.Four predator prey models with infectious diseases[J].Mathematical and Computer Modelling,2001,34:849-858.

[12]HETHCOTE H W,WANG W D,HAN L T,et al.A predator-prey model with infected prey[J].Theoretical Population Biology,2004,66:259-268.

[13]LI S,WANG X P.Analysis of stochastic predator-prey models with disease in the predator and Beddington-DeAngelis functional response[J].Advances in Difference Equations,2015,(1):1-21.

[14]ZU L,JIANG D Q,O'REGAN D.Conditions for persistence and ergodicity of a stochastic Lotka-Volterra predator-prey model with regime switching[J].Communications Nonlinear Science and Numerical Simulation,2015,29:1-11.

[15]LI D,CUI J A,SONG G H.Permanence and extinction for a single-species system with jump-diffusion[J].Journal of Mathematical Analysis and Applications,2015,430:438-464.

[16]LIU M,WANG K.Stochastic Lotka-Volterra systems with Levy noise[J].Journal of Mathematical Analysis and Applications,2014,410:750-763.

[17]ZOU X L,WANG K,Numerical simulations and modeling for stochastic biological systems with jumps[J].Communications Nonlinear Science and Numerical Simulation,2014,19:1557-1568.

(責(zé)任編輯：傅 游)

Stochastic Dynamics of a Predator-prey System with Disease in Predator

FENG Tao，MENG Xinzhu

(College of Mathematics and Systems Science，Shandong University of Science and Technology，Qingdao，Shandong 266590，China)

This paper investigated the stochastic dynamics of an infected predator-prey model with Beddington-DeAngelis functional response and Lévy jump. By using Lyapunov methods and It’s formula,this paper first discussed the existence and uniqueness of the global positive solution of the stochastic system,and then studied the asymptotic behaviors around the equilibrium points of its deterministic model. Results show that the solutions of the stochastic system fluctuate around the equilibrium points of its deterministic model under certain conditions,and the fluctuation intensity is positively correlated with the intensity of interference. Finally,numerical simulations were carried out to verify the theoretical findings.

predator-prey system; Beddington-DeAngelis functional response; asymptotic behaviors; random disturbance; Lévy jump

2016-05-16

國(guó)家自然科學(xué)基金項(xiàng)目(11371230,11501331)；山東省自然科學(xué)基金項(xiàng)目(ZR2015AQ001,BS2015SF002)；山東科技大學(xué)科研創(chuàng)新團(tuán)隊(duì)項(xiàng)目(2014TDJH102)

孟新柱(1972—)，男，山東菏澤人，教授，博士生導(dǎo)師，主要從事生物數(shù)學(xué)方面的研究，本文通信作者. E-mail:mxz721106@sdust.edu.cn

O175

A

1672-3767(2017)01-0099-12