李　婧, 竇家維

(陜西師范大學(xué) 數(shù)學(xué)與信息科學(xué)學(xué)院, 陜西 西安　710062)

?

具有擴(kuò)散及功能性反應(yīng)的周期脈沖捕食-被捕食系統(tǒng)研究

李婧, 竇家維*

(陜西師范大學(xué) 數(shù)學(xué)與信息科學(xué)學(xué)院, 陜西 西安710062)

摘要：主要研究了一類(lèi)具有擴(kuò)散及功能性反應(yīng)的周期脈沖捕食-被捕食系統(tǒng), 系統(tǒng)由一個(gè)具有脈沖的反應(yīng)擴(kuò)散方程組描述.首先通過(guò)拋物偏微分方程的不變矩形理論獲得了解的正性性質(zhì),然后利用比較原理研究系統(tǒng)解的有界性及持續(xù)生存性質(zhì);最后應(yīng)用緊性準(zhǔn)則討論了系統(tǒng)周期解的存在唯一性及全局漸近穩(wěn)定性.

關(guān)鍵詞：捕食-被捕食系統(tǒng); 功能性反應(yīng); 脈沖收獲; 持續(xù)生存性; 周期解

0引言

種群動(dòng)力學(xué)模型對(duì)種群生態(tài)學(xué)研究有著非常重要的意義,在種群動(dòng)力學(xué)中,捕食關(guān)系因其普遍存在性和重要性一直是數(shù)學(xué)與生態(tài)學(xué)界研究的主要課題[1-4].在種群動(dòng)力學(xué)模型的早期研究中,人們常假設(shè)種群在各處的密度分布是均勻的, 因此可以應(yīng)用常微分方程模型來(lái)描述種群的發(fā)展過(guò)程. 然而, 在種群的棲息地內(nèi), 如果各處種群密度分布不均勻, 那么種群密度函數(shù)不僅與時(shí)間有關(guān), 也與位置變量有關(guān), 并且種群中的個(gè)體將自發(fā)地從高密度地區(qū)向低密度地區(qū)擴(kuò)散, 這種情況下, 種群系統(tǒng)就需要用反應(yīng)擴(kuò)散方程來(lái)描述[5-8].

本文主要研究一類(lèi)具有擴(kuò)散及比率依賴(lài)功能性反應(yīng)的捕食系統(tǒng)的動(dòng)力學(xué)行為. 在捕食系統(tǒng)中考慮功能性反應(yīng)能更確切地描述食餌與捕食者的捕食關(guān)系. 功能性反應(yīng)函數(shù)有多種形式, 常見(jiàn)的有比率依賴(lài)型和食餌依賴(lài)型, 比率依賴(lài)意指功能性反應(yīng)函數(shù)是食餌與捕食者比率的函數(shù). 近年來(lái), 許多生物學(xué)家研究發(fā)現(xiàn)比率依賴(lài)型的功能性反應(yīng)更具有實(shí)用性, 因此得到了廣泛的應(yīng)用[9-15].

本文還假設(shè)由于季節(jié)變化等因素, 種群生長(zhǎng)環(huán)境是周期性變化的, 并且每隔一定的時(shí)間, 對(duì)種群進(jìn)行脈沖比例收獲[16-17].

本文所研究的系統(tǒng)將由下面具有周期系數(shù)的脈沖反應(yīng)擴(kuò)散方程組描述:

(1)

其中,u(t,x),v(t,x)分別表示食餌及捕食者的密度.系統(tǒng)的前兩個(gè)方程表示捕食系統(tǒng),Δu為L(zhǎng)aplace算子,μ1Δu,μ2Δv表示兩種群的擴(kuò)散,擴(kuò)散系數(shù)μ1,μ2均為正常數(shù).系統(tǒng)第三,四方程表示在時(shí)刻tk,k∈N+對(duì)兩種群進(jìn)行比例收獲,收獲系數(shù)分別為δ1k及δ2k;系統(tǒng)的最后兩個(gè)方程是邊界條件.

對(duì)于(1),本文假設(shè)下面條件(H)成立:

(H2)函數(shù)ai(t,x),ci(t,x)(i=1,2)以及b(t,x),r(t,x)關(guān)于t是ω-周期的;

(H3)對(duì)所有i=1,2.k∈N+,0<δik<1存在正整數(shù)p,使tk+p=tk+ω,δi(k+p)=δik.

在本文的研究中, 下面幾個(gè)基本結(jié)論將多次用到, 現(xiàn)敘述如下(參見(jiàn)[6]) 首先, 記

Qω=(0,ω]×Ω,Sω=(0,ω]×?Ω.則有

如果又有當(dāng)x∈Ω時(shí)u(0,x)不恒等于v(0,x),則u(t,x)

1系統(tǒng)的持續(xù)生存性

定理1假設(shè)條件(H)成立, 則對(duì)于系統(tǒng)(1),R2的非負(fù)象限和正象限是不變的.

(2)

由于當(dāng)(u,v)∈[0,M1]×[0,M2],且t≠tk(k∈N+)時(shí),

f1(t,x,u,v)|u=0=0,f1(t,x,u,v)|u=M1<0;

f2(t,x,u,v)|v=0=0,f2(t,x,u,v)|v=M2<0.

因此,根據(jù)文獻(xiàn)[6],并考慮到假設(shè)條件(H3),可知Π=[0,M1]×[0,M2]是系統(tǒng)(1)的不變矩形.所以,當(dāng)u(0,x)≥0,v(0,x)≥0時(shí),系統(tǒng)(1)的(u(t,x),v(t,x))滿(mǎn)足u(t,x)≥0,v(t,x)≥0,因此R2的非負(fù)象限是不變的. 又由(2),u(t,x)滿(mǎn)足

(3)

下面的結(jié)果給出了系統(tǒng)(1)解的最終有界性. 便于敘述, 記(u0,v0)=(u0(x),v0(x)), 其中u0(x)=u(0,x),v0(x)=v(0,x).

定理2假設(shè)條件(H)成立, 進(jìn)一步假設(shè)

(4)

(5)

證明: 假設(shè)(u(t,x,v(t,x))是系統(tǒng)(1)具有非負(fù)初值(u0,v0)的解.

(i)首先證明u(t,x)的最終有界性.

(6)

的解. 因?yàn)?/p>

(ii)證明v(t,x)的最終有界性.

(7)

定理3假設(shè)(H)成立,進(jìn)一步假設(shè)(5)式和下面兩個(gè)不等式成立:

(8)

(9)

證明:設(shè)(u(t,x),v(t,x))為系統(tǒng)(1)具有非負(fù)初值(u0,v0)的解. 由于(5)式和(8)式保證了定理2的條件成立, 無(wú)妨設(shè)當(dāng)t≥0時(shí),有0≤u(t,x)≤N1,0≤v(t,x)≤N2(N1,N2為兩個(gè)正常數(shù)).

(i)首先證明m1的存在性.

(10)

(ii)下面研究捕食者的持續(xù)生存問(wèn)題.

(11)

由引理1知

(12)

由條件(9), 可以選擇α>0充分小, 使得下式成立:

上面結(jié)果說(shuō)明:

2周期解的存在唯一性

為了討論系統(tǒng)(1)正周期解的存在唯一性及穩(wěn)定性, 需要下面的緊性準(zhǔn)則[16]引入記號(hào)Y=(u,v)∈Lp(Ω)×Lp(Ω),其中p>n為正整數(shù), 將系統(tǒng)(1)的前四式重寫(xiě)如下:

(13)

F(t,Y)=

引理3[16]假設(shè)函數(shù)Gk是連續(xù)可微的, 并且存在一個(gè)正值函數(shù)η(M)使

關(guān)于周期解的存在唯一性, 可以證明下面結(jié)論.

定理4假設(shè)條件(H)成立,而且系統(tǒng)是持續(xù)生存的,即存在正數(shù)β和N,當(dāng)t充分大時(shí),有β≤u(t,x)≤N,β≤v(t,x)≤N成立.若系統(tǒng)(1)還滿(mǎn)足下面不等式

(14)

其中λM為矩陣(aij)2×2的最大特征值,

則系統(tǒng)(1)有全局漸近穩(wěn)定的正ω-周期解.

證明:設(shè)(u(t,x),v(t,x)),(u1(t,x),v1(t,x))為系統(tǒng)(1)的任意兩個(gè)有界解,設(shè)其下界為β,上界為N.考慮函數(shù)

L(t)=∫Ω[(u(t,x)-u1(t,x))2+

(v(t,x)-v1(t,x))2]dx.

對(duì)其求導(dǎo), 得到:

-2μ1∫Ω|(u-u1)|2dx-

2μ2∫Ω|(v-v1)|2dx+

2∫Ω(u-u1)2(a1(t,x)-b(t,x)(u+u1)-

2∫Ω(v-v1)2(-a2(t,x)+

2∫Ω(u-u1)(v-v1)

∫Ωa11(u-u1)2dx+∫Ωa22(v-v1)2dx+

2∫Ωa12(u-u1)(v-v1)dx≤

λM∫Ω[(u-u1)2+(v-v1)2]dx.

(15)

下面進(jìn)一步討論系統(tǒng)(1)正周期解的存在唯一性. 記Y=(u,v)首先有

因此,‖Gk(Y)‖≤δ‖Y‖,其中

δ=max{δik+1,i=1,2,k=1,2,…,p}.

sup‖Y‖a≤M‖Gk(Y)‖α≤sup‖Y‖a≤Mδ‖Y‖α≤δM=∶η(M).由于引理3的條件滿(mǎn)足,因此,(u(t,x),v(t,x)),(u1(t,x),v1(t,x))在C1+v上有界,結(jié)合(15)可知:當(dāng)t→∞時(shí)

(16)

由(16)式可知, 如果系統(tǒng)(1)的周期解存在, 則周期解唯一, 并且是全局漸近穩(wěn)定的.

‖Y(ω,Y(knω)-Y(knω))‖C+

3數(shù)值例子

本節(jié)通過(guò)一個(gè)具體例子驗(yàn)證前面所得到的理論結(jié)果. 對(duì)于系統(tǒng)(1), 令

a1(t,x)=0.2sint+0.25+0.000 1cosx,

a2(t,x)=0.4,b(t,x)=0.1cost+0.2,

c1(t,x)=0.05cost+0.08,

c2(t,x)=0.45-0.03sint,

r(t,x)=15-5cost,

δ1k=0.8,δ2k=0.89,ω=2π,p=1.

顯然滿(mǎn)足條件(H), 又經(jīng)計(jì)算可知:

所以式(5)、(8)、(9)成立. 進(jìn)一步驗(yàn)證(14)式.又令β=10,N=25,算得a11=-1.101,a12=a21=0.178,a22=-0.608.所以λM=-0.550 4,進(jìn)而ln(0.82+0.892)-0.550 4×2π=-3.099 1<0. 即(14)式成立. 故滿(mǎn)足定理4的條件, 所以由定理4可知, 系統(tǒng)有全局漸近穩(wěn)定的周期解.

4結(jié)論

本文研究了一類(lèi)由脈沖反應(yīng)擴(kuò)散方程描述的捕食系統(tǒng)的動(dòng)力學(xué)行為. 研究的基本思想是運(yùn)用比較原理, 即引理1, 應(yīng)用脈沖常微分系統(tǒng)的解對(duì)脈沖反應(yīng)擴(kuò)散方程的解進(jìn)行估計(jì), 主要應(yīng)用脈沖常微分方程的相關(guān)理論和分析技巧研究解的有界性和持續(xù)生存性, 并應(yīng)用一些分析技巧獲得了周期解的存在唯一性. 在文獻(xiàn)[15]中, 作者研究了類(lèi)似問(wèn)題, 但所用的比較系統(tǒng)大多是常系數(shù)系統(tǒng),本文應(yīng)用周期脈沖常微分系統(tǒng)作為比較系統(tǒng),主要結(jié)論推廣改進(jìn)了文獻(xiàn)[15]的結(jié)論. 本文的結(jié)果可為控制生物種群的共存及發(fā)展趨勢(shì)提供理論依據(jù), 研究結(jié)果有一定的現(xiàn)實(shí)意義.

參考文獻(xiàn)

[1]肖燕妮,周義倉(cāng),唐三一.生物數(shù)學(xué)原理[M].西安:西安交通大學(xué)出版社,2008:41-50.

[2]陳蘭蓀,孟新柱,焦建軍.生物動(dòng)力學(xué)[M].北京:科學(xué)出版社,2009:87-102.

[3]Z.Amine,R.Ortega.Aperiodicpredator-preysystem[J].JournalofMathematicalAnalysisandApplications, 1994(185):477-489.

[4]盧琨.一類(lèi)具有脈沖效應(yīng)的捕食系統(tǒng)分析[J].陜西科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2011,29(3)：155-158.

[5]LakshmiNarayanGuin,PrashantaKumarMandal.Spatiotemporaldynamicsofreaction-diffusionmodelsofinteract-ingpopulations[J].AppliedMathematicalModelling,2014,38(17-18):4 417-4 427.

[6]葉其孝,李正元.反應(yīng)擴(kuò)散方程引論[M].北京:科學(xué)出版社,1990:31-57.

[7]MurryJD.Mathematicalbiology[M].NewYork:Springer-Verlag,1989.

[8]岳宗敏,白云霄,曹慧,等.庇護(hù)效應(yīng)下一類(lèi)改進(jìn)的Holling-Tanner反應(yīng)擴(kuò)散捕食模型的定性分析[J].陜西科技大學(xué)學(xué)報(bào)(自然科學(xué)版),2014,32(2)：159-163.

[9]MFan,QWang,XZou.Dynamicsofanon-autono-mousratio-dependentpredator-preysystem[J].Proc.Roy.Soc.EdinburghSect.A,2003,133:97-118.

[10]TKKar,SMisra,BMukhopadhyay.Abio-economicmodelofaratio-dependentpredator-preysystemandoptimalharvesting[J].JournalofAppliedMathematicsandComputing,2006,9(11):387-401.

[11]RArditi,LRGinzburg.Couplinginpredator-preydynamics:Ratio-dependence[J].Theret.Biol,1989,139:1 123-1 128.

[12]HRAkcakaya,RArditi,LRGinzburg.Ratio-depen-dentpredation:Anabstractionthatworks[J].Ecology,1995,76:995-1 004.

[13]XiaoHongLi,ChunLu,XiuFengDu.Permanenceandglobalattractivityofadiscretesemi-ratio-dependentpredator-preysystemwithholling-IVtypefunctionalresponse[J].JournalofMathematicalResearch&Exposition,2010,30(3):442-450.

[14]PYHPang,MWang.Qualitativeanalysisofaratio-depe-ndentpredator-preysystemwithdiffusion[J].Proc.Roy.Soc.EdinburghSect.A,2003,133:919-942.

[15]MUAkhmet,MBeklioglu,TErgenc.Animpulsiveratio-dependentpredator-preysystemwithdiffusion[J].NonlinearAnalysis:RealWorldApplications,2006,7:1 255-1 267.

[16]DDBainov,PSSimeonov.Impulsivedifferentialequations:Periodicsolutionsandapplications[M].NewYork:LongmanScientific&TechnicalPress,1993:58-79.

[17]VLakshmikantham,DDBainov,PSSimeonov.Theoryofimpulsivedifferentialequations[M].Singapore:WorldScientific,1989.

Thestudyonaperiodicimpulsivepredator-preysystem

withdiffusionandfunctionalresponse

LIJing,DOUJia-wei*

(SchoolofMathematicsandInformationScience,ShaanxiNormalUniversity,Xi′an710062,China)

Abstract:In this paper,an impulsive predator-prey system with diffusion and functional response is investigated.The system is modeled by a reaction-diffusion equation with periodic coefficients and impulses.The positive property of the solution is obtained by invariant rectangle of parabolic partial differential equations.Further, the comparison principle is applied to study the boundedness of the solution and the persistence of the system.Finally,the condition to ensure the existence and stability of the periodic solution is obtained by a compactness criterion.

Key words:predator-prey system; functional response; impulsive harvest; persistence; periodic solution

中圖分類(lèi)號(hào)：O175.1

文獻(xiàn)標(biāo)志碼：A

文章編號(hào)：1000-5811(2015)05-0186-06

作者簡(jiǎn)介:李婧(1990-),女,河南三門(mén)峽人,在讀碩士研究生,研究方向:脈沖微分方程理論及應(yīng)用

通訊作者:竇家維(1963-),女,陜西西安人,副教授,博士,研究方向： 脈沖微分方程理論及應(yīng)用，jiawei@snnu.edu.cn

基金項(xiàng)目：國(guó)家自然科學(xué)基金項(xiàng)目(61272435)

收稿日期:*2015-04-13