傅金波， 陳蘭蓀，2

(1. 福建師范大學 閩南科技學院， 福建 泉州 362332；2. 中國科學院 數(shù)學與系統(tǒng)科學研究院， 北京 100080)

食餌有病的生態(tài)-流行病模型的穩(wěn)定性分析

傅金波1， 陳蘭蓀1，2

(1. 福建師范大學 閩南科技學院， 福建 泉州 362332；2. 中國科學院 數(shù)學與系統(tǒng)科學研究院， 北京 100080)

研究一類具有雙線性發(fā)生率和功能反應(yīng)且食餌染病的生態(tài)-流行病模型的動力學行為.通過構(gòu)造適當?shù)腖yapunov函數(shù)，運用LaSalle不變集原理，獲得保證系統(tǒng)的無捕食者無病平衡點、疾病主導平衡點、捕食者主導平衡點和正平衡點全局漸近穩(wěn)定的閥值條件.通過疾病流行的閥值和捕食機制形成的閥值，以及疾病與捕食兩者競爭占優(yōu)的閥值，共同刻畫生態(tài)-流行病系統(tǒng)的演變規(guī)律性. 關(guān)鍵詞： 生態(tài)-流行病模型； Lyapunov函數(shù)； LaSalle不變集原理； 功能性反應(yīng)； 平衡點； 全局穩(wěn)定性

1 預(yù)備知識

近年來，許多學者利用種群動力學原理[1]和Kermack[2]的流行病學建模機理，分析了大量的食餌-捕食者系統(tǒng)[3-8]和多種形式的傳染病模型[9-14].然而，在已有的生物動力學模型建模中考慮疾病影響卻不多見.疾病對種群的影響是一個不可忽視的重要問題.將種群動力學與傳染病學結(jié)合起來所建立的生態(tài)-流行病模型，對于探討疾病流行對生態(tài)環(huán)境的影響，乃至改善生態(tài)環(huán)境以控制疾病的流行，無疑是一個新的研究方向，更富有實際應(yīng)用價值[15].據(jù)此，考慮疾病只在食餌之間傳播,且食餌種群分為易感者和染病者兩類，建立食餌有病的生態(tài)-流行病模型為

(1)

式(1)中：S，I分別為t時刻食餌種群中易感者和染病者兩類的數(shù)量；Y為t時刻捕食者種群的數(shù)量；A為易感者種群的增長率；β為接觸率；p為捕食率；di(i=1，2)為自然死亡率；g為捕食者種群的增長率；所有參數(shù)均為正的常數(shù)，且g

基于生物學意義，模型(1)的初值條件為S(0)>0，I(0)>0，Y(0)>0.定義3個閥值依次為

且

文中主要研究模型(1)的無捕食者無病平衡點、疾病主導平衡點、捕食者主導平衡點和正平衡點的全局漸近穩(wěn)定性，從理論上闡明捕食與疾病流行的演變規(guī)律性.

2 相關(guān)引理及其證明

引理1 模型(1)滿足初值條件的任意解(S，I，Y)皆為正解.

證明 模型(1)等價于

(2)

由此可見，在初值條件下，模型(1)的任意解(S，I，Y)皆為正解.證畢.

引理2 模型(1)滿足初值條件的任意正解(S，I，Y)均為最終有界的.

證明 令L=S+I+Y，由模型(1)得L′≤A-d1S-d1I-d2Y.取μ=min{d1，d2}，故有

由微分方程比較定理可知，存在t0>0，M≥M0，當t≥t0時，恒有L≤M.進而，集合Ω={(S，I，Y)∈R3∶0

模型(1)的非負平衡點滿足

(3)

根據(jù)方程組(3)，易獲得如下結(jié)論.

引理3 模型(1)總存在著無捕食者無病平衡點P0=(S0，0，0).當R0>1時，存在疾病主導平衡點P1=(S1，I2，0)；當R1>1時，存在捕食者主導平衡點P2=(S2，0，Y2)；當R1>1，且R2>1時，還存在正平衡點P3=(S3，I3，Y3).

3 主要結(jié)果與證明

定理1 當R0≤1，R1≤1時，模型(1)的無捕食無病平衡點P0在域Ω上是全局漸近穩(wěn)定.

證明 將模型(1)改為如下等價系統(tǒng)，即

(4)

設(shè)(S，I，Y)是系統(tǒng)(4)的任意正解，利用函數(shù)F(ξ)=ξ-1-lnξ在(0，+∞)存在唯一最小值點ξ=1且F(ξ)≥F(1)=0的性質(zhì)，構(gòu)造Lyapunov泛函為

易見正定函數(shù)V0(t)在無捕食無病平衡點P0處取得唯一最小值為零.當R0≤1，R1≤1時，直接計算V0(t)沿著系統(tǒng)(4)軌線的全導數(shù)，可得

根據(jù)LaSalle不變性原理[16]可知：模型(1)的無捕食無病平衡點P0在域Ω上是全局漸近穩(wěn)定的.證畢.

定理2 當R0>1≥R1時，模型(1)的疾病主導平衡點P1在域Ω上是全局漸近穩(wěn)定.

證明 將模型(1)改為如下等價系統(tǒng)，有

(5)

設(shè)(S，I，Y)是系統(tǒng)(5)的任意正解，利用函數(shù)F(ξ)=ξ-1-lnξ在(0，+∞)上的非負性質(zhì)，構(gòu)造Lyapunov泛函為

由此可見，正定函數(shù)V1(t)在疾病主導平衡點P1處取得唯一最小值為零.

當R0>1≥R1時，注意到βS1=d1，直接計算V1(t)沿著系統(tǒng)(5)軌線的全導數(shù)，可得

當R0>1，R1>1時，特征方程存在一個正特征值，故P1是不穩(wěn)定的；當R0>1≥R1時，特征方程的非零特征值均具有負實部，而且零特征值為特征單根，故P1是局部漸近穩(wěn)定.

由LaSalle不變性原理[16]可以得到，模型(1)的疾病主導平衡點P1在域Ω上是全局漸近穩(wěn)定的.證畢.

定理3 當R1>1≥R2時，模型(1)的捕食者主導平衡點P2在域Ω上是全局漸近穩(wěn)定.

證明 將模型(1)改為如下等價系統(tǒng)，有

(6)

設(shè)(S，I，Y)是系統(tǒng)(6)的任意正解，利用函數(shù)F(ξ)=ξ-1-lnξ在(0，+∞)上的非負性質(zhì)，構(gòu)造Lyapunov泛函為

由此可知，正定函數(shù)V2(t)在捕食者主導平衡點P2處取得唯一最小值為零.

當R1>1≥R2時，注意到gS2=d2，沿著系統(tǒng)(6)軌線計算V2(t)的全導數(shù)，可得

當R1>1，R2>1時，特征方程存在一個正特征值，故P2是不穩(wěn)定的；當R1>1≥R2時，特征方程的非零特征值均具有負實部，而且零特征值作為特征單根，故P2是局部漸近穩(wěn)定.

根據(jù)LaSalle不變性原理[16]可知：模型(1)的捕食者主導平衡點P2在域Ω上是全局漸近穩(wěn)定的.證畢.

定理4 當R1>1，R2>1時，模型(1)的正平衡點P3在域Ω內(nèi)是全局漸近穩(wěn)定.

證明 將模型(1)改為如下等價系統(tǒng)，有

(7)

設(shè)(S，I，Y)是系統(tǒng)(7)的任意正解，利用函數(shù)F(ξ)=ξ-1-lnξ在(0，+∞)上的非負性質(zhì)，構(gòu)造Lyapunov泛函為

易于驗證正定函數(shù)V3(t)在正平衡點P3處取得唯一最小值為零.

同時，模型(1)在P3處Jacobian矩陣的特征方程為

其中：a1=d1+βI3+pY3;a2=β2S3I3+pgS3Y3+pgI3Y3;a3=pgI3Y3(d1+βI3+pY3).

又因為Δ1=a1，Δ2=a1a2-a3=(d1+βI3+pY3)(β2S3I3+pgS3Y3)，Δ3=a3Δ2皆為正數(shù)，故由Routh-Hurwitz判據(jù)[16]可知：當R1>1，R2>1時，P3總是局部漸近穩(wěn)定.

由LaSalle不變性原理[16]可知：模型(1)的正平衡點P3在域Ω內(nèi)是全局漸近穩(wěn)定的.證畢.

4 結(jié)論

由定理1可知：當食餌種群的易感者染病數(shù)量不多于1個，且捕食者種群的捕食數(shù)量也不多于1個時，該生態(tài)-流行病系統(tǒng)中的疾病尚未流行，且捕食機制也尚未形成，疾病和捕食者種群將在該系統(tǒng)中很快被消除，食餌種群的易感者數(shù)量將全局漸近穩(wěn)定在一個正常數(shù)上.

由定理2可知：當食餌種群的易感者染病數(shù)量大于1個，且捕食者種群的捕食數(shù)量不多于1個時，該生態(tài)-流行病系統(tǒng)中疾病流行，且捕食機制尚未形成，捕食者種群在系統(tǒng)中趨于滅絕，食餌種群的易感者和染病者數(shù)量將全局漸近穩(wěn)定在一組正常數(shù)上.

由定理3可知：當捕食者種群的捕食數(shù)量大于1個，且疾病低于捕食對食餌種群的侵害，該生態(tài)-流行病系統(tǒng)中疾病尚未流行，且捕食機制已形成，疾病在系統(tǒng)中很快消除，食餌種群的易感者和捕食者數(shù)量將全局漸近穩(wěn)定在一組正常數(shù)上.

由定理4可知：當捕食者種群的捕食數(shù)量大于1個，且疾病還高于捕食對食餌種群的侵害，該生態(tài)-流行病系統(tǒng)中捕食機制已形成且疾病流行，食餌種群的易感者、染病者和捕食者種群共存并將全局漸近穩(wěn)定在一組正常數(shù)上.

綜上討論，R0是疾病是否流行的閥值，R1是捕食機制是否形成的閥值，R2是疾病與捕食兩者占優(yōu)比較的閥值，三者共同刻畫了生態(tài)-流行病系統(tǒng)(1)的演變規(guī)律性.

[1] 陳蘭蓀，陳鍵.非線性生物動力系統(tǒng)[M].北京：科學出版社，1993：102-130.

[2]KERMACKWO，MCKENDRICKAG.Contributionstothemathematicaltheoryofepidemics[J].BulletinofMathematicalBiology，1991，53(1):33-55.

[3] 李林.一個環(huán)境數(shù)學模型的一致持久性與穩(wěn)定性[J].應(yīng)用數(shù)學學報，2003，26(1)：149-157.

[4] 程榮福，李輝.一個具功能性反應(yīng)的微分生態(tài)系統(tǒng)的定性分析[J].東北師大學報(自然科學版)，2005，37(1):11-15.

[5]CHENFengde，SHIJinlin.Periodicityinalogistictypesystemwithseveraldelays[J].ComputersandMathematicswithApplications，2004，48(1/2):35-44.

[6] 張正球，王志成.基于比率的三種群捕食者-食餌擴散系統(tǒng)的周期解[J].數(shù)學學報，2004，47(3):531-540.

[7] 范猛，王克.一類具有HollingⅡ型功能性反應(yīng)的捕食者-食餌系統(tǒng)全局周期解的存在性[J].數(shù)學物理學報，2001，21(4):83-91.

[8] 徐瑞，陳蘭蓀.具有時滯和基于比率的三種群捕食系統(tǒng)的持久性與全局漸近穩(wěn)定性[J].系統(tǒng)科學與數(shù)學，2001，21(2):204-212.

[9]YOSHIDAN，HARAT.GlobalstabilityofadelayedSIRepidemicmodelwithdensitydependentbirthanddeathrates[J].ComputationalandAppliedMathematics，2007，201(2):339-347.

[10] 侯強，靳禎.一個基于修正的Leslie-Gower 生態(tài)傳染病模型的研究[J].中北大學學報(自然科學版)，2012，33(2):98-101.

[11] 張少輝，靳禎.具有非線性發(fā)生率的傳染病模型性態(tài)分析[J].中北大學學報(自然科學版)，2012，33(4):353-357.

[12] 劉玉英，肖燕妮.一類受媒體影響的傳染病模型的研究[J].應(yīng)用數(shù)學和力學，2013，34(4):399-407.

[13] 傅金波，陳蘭蓀，程榮福.具有潛伏期和免疫應(yīng)答的時滯病毒感染模型的全局穩(wěn)定性[J].高校應(yīng)用數(shù)學學報(A輯)，2015，30(4):379-388.

[14] 傅金波，陳蘭蓀，程榮福.具有Logistic增長和治療的SIRS傳染病模型的后向分支[J].吉林大學學報(理學版)，2015，53(6):1166-1170.

[15] LI Jianquan，ZHANG Juan，MA Zhien.Global analysis of some epidemic models with general contact rate and constant immigration[J].Applied Math Mech，2004，25(4):396-404.

[16] 馬知恩，周義倉.傳染病動力學的數(shù)學建模與研究[M].北京:科學出版社，2004:178-192.

(責任編輯： 錢筠 英文審校： 黃心中)

Stability Analysis of Eco-Epidemiologic Model With Disease in Prey

FU Jinbo1， CHEN Lansun1，2

(1. Minnan Science and Technology Institute, Fujian Normal University， Quanzhou 362332， China；2. Academy of Mathematics and Systems Science， Chinese Academy of Sciences， Beijing 100080， China)

In this paper， the dynamical behaviors of the eco-epidemiologic model with double linearity incidence rate， functional response and disease in the prey are studied. By constructing suitable Lyapunov function and using LaSalle invariance principle， the global asymptotic stable threshold conditions of non-predator disease-free equilibrium， disease dominant equilibrium， predator dominant equilibrium and positive equilibrium in the system are obtained. The threshold of disease popularity， the threshold of formation of predation mechanism and the threshold of dominance of disease or predator in their competition depict the evolvement law of the eco-epidemiologic system. Keywords： eco-epidemiologic model； Lyapunov function； LaSalle invariance principle； functional response； equilibrium point； global stability

10.11830/ISSN.1000-5013.201702025

2016-09-09

傅金波(1978-)，男，副教授，主要從事生物數(shù)學的研究.E-mail:fujinbomnkjxy@sina.com.

國家自然科學基金資助項目(11371306)； 福建省教育廳自然科學基金資助項目(JA13370)； 福建師范大學閩南科技學院青年骨干教師重點項目(MKQ201006)

O 175.13

A

1000-5013(2017)02-0266-05