一類帶有不育控制的單種群離散模型

李秋英， 張鳳琴， 劉漢武

(運城學院應用數學系 山西運城044000)

0 引言

鼠類在農、林、牧業等方面造成的損失遠大于病害和蟲害，它們還是多種病菌的寄主和傳播者，因此對害鼠的有效控制越來越重要．隨著對鼠類研究的深入，科學家們已經不再局限于研究鼠類對人類生產、生活的危害了，并展開鼠類對于生態系統的影響以及對鼠類的不育控制方面的研究．

不育控制是借助某種技術或方法使雄性或雌性絕育，或阻礙胚胎著床發育，甚至阻斷幼體生長發育，以降低生育率，從而控制種群數量的增長［1］．不育控制有許多優點:不會對環境造成污染，具有針對性，能大面積使用，在經濟上花費更少，效率更高等［2］．因此公眾和動物權益組織認為不育控制是人道的．目前，關于不育控制技術的研究主要分為2方面:一方面為藥效的研究即野外控制實驗［3］;另一方面為不育控制下害鼠種群動態的系統分析［4－5］．

種群動態的理論分析方面，一些學者采用微分方程模型［5－7］，也有學者采用差分模型研究種群動態．張知彬指出［4］，不育控制之前應從理論上分析，同時建立數學模型，并對其進行分析，預見不育控制下害鼠的種群動態，評價其效果，并提出指導意見．

1 模型建立

本文中，將害鼠種群分為可育者F和不育者S兩類．為建立模型，首先做以下假設:

1)生育函數為B(N)=b(1－rN)1r，b ＞0，r＜0．

2)不育控制方式是非自傳播方式．

3)單位時間內可育個體和不育個體的死亡率分別為d1，d2(d1，d2＜b)．

4)單位時間內可育個體到不育個體的轉化率為α．

5)單位時間不育個體恢復率(即不育的一部分會變為可育的)為δ．

基于以上假設，建立不育控制下的單種群模型為

由于一些害鼠物種之間沒有連續的重疊代，在(1)中令微分步長為1，利用歐拉法得到離散種群模型

根據實際意義假設1－d1－α＞0，1－d2－δ＞0，(2)的初始函數滿足 F(0)＞0，S(0)≥0．

2 平衡點的局部穩定性分析

定理1(i)當 R0=(d1(d2+δ)+d2α)/(b(d2+δ))＞1 時，(2)僅存在零平衡點 E0=(0，0);

(ii)當 R0=(d1(d2+δ)+d2α)/(b(d2+δ))＜1時，(2)除平衡點 E0外還存在正平衡點 E*(F*，S*)，其中，F*=((1－Δ)/r)·(d2+δ)/(d2+α+δ)，S*=((1－Δ)/r)·α/(d2+α+δ)，Δ =Rr0．

下面介紹差分方程平衡點穩定性的一般準則．

若某二維離散系統在其平衡點E處的線性化系統為X(m+1)=WX(m)，其中W為該差分系統在平衡點E處的Jacobi矩陣．若W特征值的絕對值都小于1，則平衡點E是局部漸近穩定的．即W同時滿足3個Jury條件:(i)1－tr W+det W＞0;(ii)1+tr W+det W＞0;(iii)1－det W＞0．

定理2 當R0＞1時，(2)的平衡點E0=(0，0)是局部漸近穩定的．

證明 令 ψ=1－d1－α，μ=1－d2－δ．則(2)在E0=(0，0)處的Jacobi矩陣

由矩陣W0得tr W0=b+ψ+μ＞0，det W0=(b+ψ)μ－δα＞0．于是

顯然，1+tr W+det W＞0，從而滿足3個Jury條件，因此平衡點E0是局部漸近穩定的．

定理3 當R0＜1時，(2)的正平衡點E*是局部漸近穩定的．

證明 (2)在E*處的Jacobi矩陣為

其中C=bF*=b(d2+δ)/(d2+δ+α)·R0·(1－R0－r)/(－r)．令 f(x)=(1－x－r)/(－r)，x＞d1/b．顯然有 f'(x)= －x－r－1＜0．故當 x＞d1/b時，f(x)≤f(d1/b)=(1 －(d1/b)－r)/(－r)．令 g(t)=1 －(d1/b)t－t(t＞0)．從而 t＞0 時，g'(t)=ln(d1/b)(d1/b)t－1 ＜0，且有 g(0)=0．因此 g(t)＜0，即 g(t)=(1－(d1/b)t)/t＜1(t＞0)．由于(1－R0－r)/(－r)＜1，從而 C=(b(d2+δ)/(d2+δ+α))·R0·(1－R0－r)/(－r)＜(d1(d2+δ)+d2α)/(d2+δ+α)．根據矩陣W*得tr W*=1－αδ/(d2+δ)－C+μ，det W*=(1－αδ/(d+δ)－C)μ－δα+αC．由Jury條件可知，當(2)在正平衡點 E*處的 Jacobi矩陣 W*滿足 tr W*＜1+

2det W*＜2時，正平衡點E*是局部漸近穩定的．而 tr W*=1－αδ/(d2+δ)－C+μ＞1－(d1(d2+δ)+d2α)/(d2+δ+α)－αδ/(d2+δ)＞0．因此只需證明 tr W*＜1+det W*＜2．因此1－tr W*+det W*=C(d2+δ)+αC ＞0，1－det W*=(d+δ)(1－C)+(1－α)C+αδ/(d+δ)＞0．故當R＜1時，tr W*＜1+det W*＜2．即

220

滿足Jury的3個條件，所以正平衡點E*(F*，S*)是局部漸近穩定的．

3 主要結果

定理4 當R0＞1時，(2)的平衡點E0是全局漸近穩定的．

證明 因為 R0＞1，所以 b－d1＜αd2/(d2+δ)＜1．定義 Lypunov函數 V(F，S)=F(n)+S(n)．由(2)得

其中Θ=min{b－d1，d2}，顯然0＜Θ＜1．從而平衡點E0是全局漸近穩定的．為證明(2)的解是持續生存的，引進引理1．考慮

其中 ω≥0，0＜d＜1 是正常數．則對(3)有引理1［8］．

引理1 令u(n)是方程(3)滿足u(0)＞0的解．則

(a)nl→im∞u(n)=ω/d;(b)任給常數 ε ＞0，M ＞0，存在正常數δ(ε)(ε，M)使得當

定理5 當δ=0，R0＜1時，(2)的解是一致持續生存的．

證明 易證明(2)是耗散的．即存在常數M＞0，使得集合Ω={(F，S )0≤Fn≤M，0≤Sn≤M}為(2)的正向不變集．

考慮方程

令 v(n，n0，v0)為(4)的解．根據引理 1，對于 n0∈Z+，0≤v0≤M，給定 ε1＞0 和正常數 M，存在常數 δ0=δ0(ε1)＞0和 n*=n*(ε1，M)＞0 使得 e(n)＜ δ0，n≥n0+n*時，0≤v(n，n0，v0)＜ ε1．由本定理條件知，存在常數 η 和 ε1滿足(1－((d1+α)/b)1r)/r＞ε1+η，αη ＜ δ0．即

下面證明存在常數β使得對于(2)的任意解有

首先證明(2)的任意解(F(n)，S(n))有

由解的有界性知對每個 m，存在 K使得當 n＞K時，F(n，θ)＜M，S(n，θ)＜M．又因為當 q→+∞時，s(qm)→+∞．所以存在K1(m)，使得q＞K1(m)時，s(qm)＞K(m)．令q≥K1(m)，故n∈［s(qm)+1，t(qm)－1］時，，從而可得

故當 q≥K1(m)，m=1，2，…時，t(qm)－s(qm)≥(ln(m+1))/(－ ln γ)．進而，可取，使得 m≥．q≥K1(m)時，有t(qm)－s(qm)≥n*．對任意的 m≥，q≥K1(m)及 n∈［s(qm)，t(qm)］，Sn+1≤α·η/(m+1)+(1 －d2－δ)Sn成立．設v(n)是(4)滿足 v(s(qm))=S(s(qm))的解，則當 n∈［s(qm)，t(qm)］，m≥，q≥K1(m)時，有 S(n，θ(m))≤v(n)．取n0=s(qm)，v0=S(s(qm))，由于 0＜v0＜M，αη ＜ δ1，因此可得 n∈［s(qm)+n*，t(qm)］時，方程(4)的解 v(n，s(qm)，u(s(qm)))滿足 v(n)=v(n，s(qm)，u(s(qm)))＜ ε1．故 S(n，θ(m))＜ ε1，n∈［s(qm)+n*，t(qm)］，q≥K1(m)，m≥．根據(2)的第一個方程和比較定理，可得當m≥，q≥K1(m)，n∈［s(qm)+n*，t(qm)］時，

成立．因此η/(m+1)2≥F(s(qm)，θ(m))≥F(s(

qm)－1，θ(m))(b(1 －rη －rε1)1r+1－d1－ α)＞ α0/(m+1)2．

這導致矛盾．因此，存在η＞0使得對于(2)任意的解(F(n)，S(n))，有 lnim inf F(n)＞η．類似可得，存在 η'，

→∞使得．從而定理得證．

4 結論

本文主要討論了不育控制的作用．從定理4～5可以看出，當不育轉化率α比較小，使得R0＞1時，(2)的解是持續生存的，且可以借助Matlab軟件可得正平衡點是全局漸近穩定的;種群規模趨向于(1－((d1(d2+δ)+d2α)/(b(d2+δ)))r)/r．當α比較大，使得R0＜1時，種群逐漸消失．由定理4～5得，若不育者的恢復率δ比較小，使得R0＜1時，種群是滅絕的．若不育劑失效了，導致恢復率變大且使得R0＜1，從而種群將會持續生存．所以在對害鼠進行控制時，應使α比較大，而δ比較?。?/p>

［1］ 陳東平，王曉．鼠類不育技術控制鼠害的理論與實踐［J］．預防醫學情報雜志，2005，21(2):163－165．

［2］ 張濤，劉銘泉．不育劑在滅鼠中的應用研究［J］．中國媒介生物學及控制雜志，2005，16(4):327－328．

［3］ 宛新榮，石巖生，寶祥，等．不育劑對黑線毛足鼠種群繁殖的影響［J］．獸類學報，2006，26(4):392－397．

［4］ 張知彬．鼠類不育控制的生態學基礎［J］．獸類學報，1995，15(3):229－234．

［5］ 劉漢武，周立，劉偉，等．利用不育技術防治高原鼠兔的理論模型［J］．生態學雜志，2008，27(7):1238－1243．

［6］ 郝麗杰，蔣貴榮，鹿鵬．具有垂直傳染的SIRS傳染病模型分岔分析［J］．鄭州大學學報:理學版，2013，45(2):31－36．

［7］ 張美明，張鳳琴，呂江，等．具有性別結構的免疫不育控制模型［J］．信陽師范學院學報:自然科學版，2011，24(4):425－429．

［8］ Xu Jiabo，Teng Zhidong．Permanence for nonautonomous discrete single-species system with delays and feedback control［J］．Applied Mathematics Letters，20140，23(9):949 －954．