留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

具有分数阶扩散的捕食-食饵模型的共存性

上一篇

下一篇

柳文清, 陈清婉, 傅金波. 具有分数阶扩散的捕食-食饵模型的共存性[J]. 西南师范大学学报(自然科学版), 2020, 45(3): 16-20. doi: 10.13718/j.cnki.xsxb.2020.03.003
引用本文: 柳文清, 陈清婉, 傅金波. 具有分数阶扩散的捕食-食饵模型的共存性[J]. 西南师范大学学报(自然科学版), 2020, 45(3): 16-20. doi: 10.13718/j.cnki.xsxb.2020.03.003
Weng-qing LIU, Qing-wan CHEN, Jin-bo FU. On Coexistence of a Fractional Diffusion Predator-Prey System[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(3): 16-20. doi: 10.13718/j.cnki.xsxb.2020.03.003
Citation: Weng-qing LIU, Qing-wan CHEN, Jin-bo FU. On Coexistence of a Fractional Diffusion Predator-Prey System[J]. Journal of Southwest China Normal University(Natural Science Edition), 2020, 45(3): 16-20. doi: 10.13718/j.cnki.xsxb.2020.03.003

具有分数阶扩散的捕食-食饵模型的共存性

  • 基金项目: 福建省中青年教师教育科研项目(JAT191035,JAT191044);泉州科技高层次人才创新创业项目(2018C094R);福建省高校“新世纪优秀人才支持计划”(2018)
详细信息
    作者简介:

    柳文清(1984-), 男, 讲师, 主要从事偏微分方程研究 .

  • 中图分类号: O175.2

On Coexistence of a Fractional Diffusion Predator-Prey System

计量
  • 文章访问数:  1051
  • HTML全文浏览数:  1051
  • PDF下载数:  114
  • 施引文献:  0
出版历程
  • 收稿日期:  2018-07-02
  • 刊出日期:  2020-03-20

具有分数阶扩散的捕食-食饵模型的共存性

    作者简介: 柳文清(1984-), 男, 讲师, 主要从事偏微分方程研究
  • 闽南科技学院, 福建 泉州 362300
基金项目:  福建省中青年教师教育科研项目(JAT191035,JAT191044);泉州科技高层次人才创新创业项目(2018C094R);福建省高校“新世纪优秀人才支持计划”(2018)

摘要: 研究了一类分数阶扩散且具有B-D反应函数的捕食-食饵模型,通过构造Lyapunov函数,证明了平衡点的局部渐近稳定性和全局渐近稳定性;利用Leray-Schauder拓扑度方法,证明了满足一定条件时,非常数正平衡解存在.

English Abstract

  • 在种群动力学中,一般用算子Δ来表示种群从高密度处向低密度处扩散,这种扩散称之为正常的扩散,具有连续性,对于此类问题的研究已有丰富的成果[1-6].包括研究解的稳定性[1-2]、共存解问题[3-4]、行波解[5-6]等等.近来,有研究表明一些种群在空间扩散上有某种跳跃性[7-8],可以用算子(-Δ)αα∈(0,1)近似表示,为此考虑如下分数阶扩散捕食-食饵模型

    其中:uv分别表示捕食者和食饵在(xt)处的密度,边界条件表示系统是封闭的;r表示食饵内禀增长率;eβ分别表示捕食者死亡率和食饵被捕食后转化为食物的转化率.捕食关系利用B-D反应函数$\frac{uv}{1+au+bv}$来表示,模型中的参数均为正常数.

    算子(-Δ)αα∈(0,1)的特征问题为

    (-Δ)α定义域:

    这里μkα所对应的特征函数为ϕk,并且0=μ1α < μ2α < … < μkα +∞.

    由文献[9]可知,算子(-Δ)αu可定义为

    并且有如下计算结果[10]

    通过简单计算,系统存在唯一边界平衡点:E0=(k,0),并且当满足条件

    存在唯一的非常数正平衡点E*=(u*v*),其中

  • 定义基本再生数$R=\frac{r\beta }{kbe}$,易得R>1时,$\frac{r}{k}>\frac{b{{v}^{*}}}{1+a{{u}^{*}}+bv}$.

    定理1   当R>1时,平衡点E*是局部渐近稳定的.

      系统(1)在平衡点E*处的线性化系统为

    其中

    系统(2)所对应的特征方程为

    其中:

    由定理条件易知$\frac{r}{k}>\frac{b{{u}^{*}}}{{{\left(1+a{{u}^{*}}+b{{v}^{*}} \right)}^{2}}}$,则B1>0,B2>0,这说明特征方程的根均有负实部,所以平衡点E*是局部渐近稳定的.

    定理2   当R>1时,平衡点E*是全局渐近稳定的.

      定义Lyapunov函数

    其中$h=\frac{1+b{{v}^{*}}}{\beta +\beta au*}$,则有

    这里

    当$\frac{r}{k}>\frac{b{{u}^{*}}}{1+a{{u}^{*}}+b{{v}^{*}}}$时,$\frac{\text{d}V(t)}{\text{d}t}\le 0$,故平衡点E*是全局渐近稳定的.

  • 考虑系统(1)所对应的椭圆方程组

    引理1   存在常数CC,使系统(2)任一非负解(uv)满足

      根据比较原理[11],由方程组(3)中第二个方程容易得到,maxΩv(x)≤k.再将方程组(3)中两个方程相加,同样利用比较原理可得,${{\max }_{\Omega }}u(x)\le \beta k\left(1+\frac{r}{e} \right)$.取$\bar{C}=\max \left\{ k, \beta k\left(1+\frac{r}{e} \right) \right\}$即可证明uv的上界.

    下面证明uv的下界,将系统(3)中的第一个方程两端同乘以u并在Ω上积分,注意到

    则有

    由积分中值定理可知,存在x0Ω,使$\frac{\beta v\left({{x}_{0}} \right)}{1+au\left({{x}_{0}} \right)+bv\left({{x}_{0}} \right)}\ge e$,从而$v\left({{x}_{0}} \right)\ge \frac{e}{\beta }$.再由Hanack不等式[11]可得

    下面利用反证法证明u的下界.假设存在序列Wm=(umvm),使得m ∞时,min xum(x) 0.结合Hanack不等式,max xum(x) 0.则um(x) 0在Ω一致成立.另外,解序列Wm=(umvm)满足方程组(3),由方程组(3)中第二个方程积分

    m→∞,则vmk.同理将方程组(3)第一个式子积分两边同除以$\underset{x\in \bar{\Omega }}{\mathop{\max }}\, {{u}_{m}}(x)$积分可得

    由于$\frac{{{u}_{m}}}{\underset{x\in \bar{\Omega }}{\mathop{\max {{u}_{m}}}}\, (x)}>\frac{\underset{x\in \bar{\Omega }}{\mathop{\min {{u}_{m}}}}\, (x)}{\underset{x\in \bar{\Omega }}{\mathop{\max {{u}_{m}}}}\, (x)}\ge \frac{1}{{{A}_{1}}}$,再令m→∞,可得$\frac{\beta k}{1+bk}=e$,另一方面,由条件H1计算可得$\frac{\beta k}{1+bk}>\alpha $,矛盾.因此uvC立.引理证毕.

    记系统(3)为-(Δ)α W=F(W).令$B=\{x\in \Omega |u, v\in (\underset{\raise0.3em\hbox{$\smash{\scriptscriptstyle-}$}}{C}, \bar{C})\}$,W=(uv)且W满足方程

    其中(I-D)-1I-D的逆.由于算子G是单位算子I的一个等干扰算子.对∀W∈$\partial $BG(W)≠0,故可定义Leray-Schander拓扑度deg(G(·),0,B).直接计算得:

    假设$\ell _ { 1 } < 0, \ell _ { 1 } ^ { 2 } - 4 \ell _ { 2 } > 0$,此时H(μα)=0有两正根$\begin{array} { l } \wedge \quad \wedge \\ \mu _ { 1 }, \mu _ { 2 } \end{array}$,并且存在两正整数k1k2,使得${{\hat{\mu }}_{1}}\in \left[\mu _{{{k}_{1}}}^{\alpha } \right., \left. \mu _{{{k}_{1}}+1}^{\alpha } \right), {{\hat{\mu }}_{2}}\in \left(\mu _{{{k}_{2}}}^{\alpha }, \mu _{{{k}_{2}}+1}^{\alpha } \right).$.

    由Leray-Schander定理,有以下引理[9].

    引理2   对所有正整数i都有H(μiα)≠0,则

    定理3   若满足条件$\ell _ { 1 } < 0, \ell _ { 1 } ^ { 2 } - 4 \ell _ { 2 } > 0$,且$\sum\nolimits_{i={{k}_{1}}+1, H\left(\mu _{i}^{\alpha } \right) < 0}^{{{k}_{2}}}{\rm dim}E\left(\mu _{i}^{\alpha } \right)$为奇数,系统(3)至少有一个非常数正平衡解.

      利用反证法.设d为充分大的正数,定义

    考虑如下方程

    同理(5)式也可以记为

    t=1时,系统(5)即是系统(3),并且(5)只有唯一正常数解E*.由定理条件以及引理2可知

    另一方面,当t=0时,$H\left(0, {{\mu }^{\alpha }} \right)={{d}^{2}}{{\mu }^{2\alpha }}+d{{\ell }_{1}}{{\mu }^{\alpha }}+{{\ell }_{2}}$,由于${{\ell }_{2}}$,所以只要取足够大的d>0,满足$2{{d}^{2}}{{\mu }^{2\alpha }}+d{{\ell }_{1}}>0$,则对所有μ ≥0有H(0,μα)>0,故此时deg(F(D,·),E*)=(-1)0=1.由Leray-Schander度的同伦不变性[12]可得矛盾.这说明系统(3)中至少有一个不同于E*的解.定理证毕.

参考文献 (12)

目录

/

返回文章
返回