具有扩散的生态-流行病模型的稳定性分析
2024-01-20王玲书陈若彤
王玲书,陈若彤
(河北经贸大学数学与统计学学院,河北 石家庄 050061)
0 引 言
扩散是自然界中普遍存在的现象,扩散可以使种群的个体迁出和迁入,增加或降低当地种群的密度,这对理解种群的进化具有重要意义。在生态系统中,捕食关系是种群之间相互作用的最普遍、最重要的形式之一。捕食系统是描述自然环境中捕食者以食饵为食的物种间关系的模型,捕食者或食饵具有扩散的捕食模型已经被许多学者所研究(参见文献[1-3])。在文献[1]中,作者研究了食饵具有扩散的捕食者-食饵模型的稳定性,并讨论了由捕食者的妊娠期引起的时滞对模型稳定性的影响。
ANDERSON和MAY在文献[4]中首次将种群模型和传染病模型结合起来提出了一个食饵染病的生态-流行病模型以来,生态-流行病模型已经被许多学者所研究(如文献[5-8])。由于传染病可以在捕食者之间、食饵之间或捕食者和食饵之间传播,从而使得生态-流行病模型有着丰富的研究内容。文献[5]假设疾病仅在捕食者之间传播,研究了下列生态-流行病模型:
(1)
其中:x(t)、S(t)和I(t)分别表示食饵、易感捕食者和染病捕食者种群在时刻t的密度。文献[5]讨论了模型(1)的非负平衡点的局部稳定性,给出了正平点处存在Hopf分支的充分条件。
基于文献[1]和[5],本文研究一个食饵具有扩散和捕食者染病的捕食-被捕食病模型,并讨论由捕食者种群的妊娠期引起的时滞对种群动力学性态的影响。为此,我们研究下列生态-流行病模型
(2)
其中:x1(t)和x2(t)分别表示斑块1和斑块2中食饵在时刻t的密度;S(t)和I(t)分别表示斑块1中易感捕食者和染病捕食者在时刻t的密度.参数r1、a11、a12、a21、d12、d21、r2、d1、d2和β均为正常数,其中r1表示斑块1中食饵的内禀增长率,a11为斑块1中食饵的种内竞争率,a12表示捕食者种群的捕获率,a21/a12表示捕食者种群的生育转化率,d12和d21表示食饵种群在两个斑块之间的扩散率,r2是斑块2中食饵种群的死亡率,d1和d2分别表示易感捕食者和染病捕食者的死亡率,β为疾病的传染率;τ≥0表示捕食者种群的妊娠期所引起的时滞。模型(2)假设生态系统由两个孤立的斑块组成,食饵种群在两个斑块间扩散而捕食者种群仅在斑块1中捕食食饵。
本文将基于以下初始条件来讨论模型(2)的稳定性
(3)
其中φi(θ)和φi(θ)(i=1,2)均为[-τ,0]上的非负连续函数。由泛函微分方程的基本理论[9]可知,模型(2)存在满足初始条件(3)的唯一正解。
1 无病平衡点的稳定性
模型(2)总存在一个平凡平衡点EO(0,0,0,0)。当r1(r2+d21)>r2d12时,模型(2)有一个捕食者灭绝平衡点E1(x′1,x′2,0,0),其中
若下列条件成立:
进一步地,当下列条件成立时:
(λ+d2-βS+)[λ3+p2λ2+p1λ+p0+(q2λ2+q1λ+q0)e-λτ]=0 。
(4)
显然,方程(4)有一个实根λ=-(d2-βS+)=-d2(1-R1)。方程(4)的其他根由下列方程决定
λ3+p2λ2+p1λ+p0+(q2λ2+q1λ+q0)e-λτ=0 ,
(5)
其中
当τ=0时,方程(5)变为
λ3+(p2+q2)λ2+(p1+q1)λ+p0+q0=0 。
(6)
直接计算可得
假设方程(6)有纯虚根iv(v>0),则有
v6+h2v4+h1v2+h0=0 ,
(7)
其中
当h0>0且h1>0时,方程(7)的根均具有负实部,因此由Hurwitz判定定理可知,E2是局部渐近稳定的。当h1>0且h0<0时,方程(7)存在一个正实根v0,相应地,方程(6)存在一对纯虚根±iv0,直接计算可得
计算方程(5)关于τ的导数可得
直接计算有
综上分析,由文献[9]可以得到下面结论:
定理1对模型(2),设R0>1,则有下面结论:
(ii)如果R1>1,则E2是不稳定的。
证明 :设(x1(t),x2(t),S(t),I(t))为系统(2)的满足初始条件(3)的任一正解,构造函数
沿模型(2)的解计算V21(t)的导数可得
定义
直接计算可得
2 共存平衡点的稳定性
λ4+P3λ3+P2λ2+P1λ+P0+(Q3λ3+Q2λ2+Q1λ)e-λτ=0,
(8)
其中
当τ=0时,方程(8)变为
λ4+(P3+Q3)λ3+(P2+Q2)λ2+(P1+Q1)λ+P0=0 。
直接计算可得
Δ4=P0Δ3。
若Δ2>0,Δ3>0,则有Δ4>0.由Routh-Hurwitz判断准则可得,平衡点E*是局部渐近稳定的。
如果λ=±iω(ω>0)是方程(8)的一对纯虚根,将λ=iω代入式(8),可得下列等式
ω8+H3ω6+H2ω4+H1ω2+H0=0 ,
(9)
其中
设z=ω2,方程(9)变为
(10)
计算方程(8)关于τ的导数可得
直接计算有
综上分析,由文献[10]可以得到下面结论:
定理3对于模型(2),假设R1>1,Δ2>0和Δ3>0,我们可以得到下面结论:
定理4假设R1>1,则当下列条件成立时
证明:设(x1(t),x2(t),S(t),I(t))为模型(2)满足初始条件(3)的任一正解。构造函数
定义
直接计算可得
解限定在M上,由模型(2)的第3个方程可得
注1对于模型(2),用同样的方法可以证明:
(ⅰ)当r2d12>r1(r2+d21)时,模型(2)的平凡平衡点E0(0,0,0,0)是全局渐近稳定的,即捕食者和食饵种群都将灭绝;当r2d12 (ⅱ)设r2d12 注2本文研究了一个食饵种群具有扩散、捕食者种群染病的生态流行病模型的稳定性,并讨论了由捕食者种群的妊娠期引起的时滞对种群增长的影响。通过构造适当的Lyapunov泛函,我们分别建立了保证模型(2)的无病平衡点和共存平衡点全局吸引的充分条件。