Advances in Applied Mathematics
Vol. 08  No. 10 ( 2019 ), Article ID: 32655 , 7 pages
10.12677/AAM.2019.810191

Classification and Stability Analysis of Equilibrium Points in a Leslie Predator-Prey System with Generalized Holling Type III Functional Response

Chenrong Pan, Songlin Chen

School of Mathematical Science and Engineering, Anhui University of Technology, Ma’anshan Anhui

Received: Oct. 1st, 2019; accepted: Oct. 17th, 2019; published: Oct. 24th, 2019

ABSTRACT

In this paper, we consider a Leslie predator-prey system with Holling type III functional response. We construct a suitable variable transformation for the original system, and then obtain an equivalent system. By investigating the classification and stability of the equilibrium points of the equivalent system, we derive the results on the classification and stability of the equilibrium points of the original system.

Keywords:Leslie Predator-Prey System, Variable Transformation, Equilibrium Point, Stability

具有广义Holling III型功能反应的 Leslie捕食–食饵系统的平衡点分类 和稳定性分析

潘陈蓉,陈松林

安徽工业大学,数理科学与工程学院,安徽 马鞍山

收稿日期:2019年10月1日;录用日期:2019年10月17日;发布日期:2019年10月24日

摘 要

在本文中,我们考虑一个具有Holling III型功能反应的Leslie型捕食–食饵系统。我们通过对原系统构造合适的变量变换,将其化为一个等价的系统。通过对等价系统的平衡点分类和稳定性进行分析,我们得出了原系统的平衡点分类和稳定性结果。

关键词 :Leslie型捕食者–食饵系统,变量变换,平衡点,稳定性

Copyright © 2019 by author(s) and Hans Publishers Inc.

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

1. 引言

考虑如下具有功能反应的Leslie捕食–食饵系统 [1] :

{ d x d t = r x ( 1 x K ) y p ( x ) , d y d t = s y ( 1 y h x ) . (1)

其中, x , y 分别表示食饵种群和捕食者种群的密度; r , s > 0 ,分别表示食饵和捕食者的内在增长率; K > 0 ,表示食饵的环境承载力; h > 0 1 h 是支撑一个捕食者所需的食饵数量。根据生态系统的意义,我们在区域 G = { ( x , y ) | x > 0 , y 0 } 上讨论 [2]。

近几年来,Ali Atabaigi等研究了由Holling IV型或Monod-Haldane功能反应产生的Leslie型捕食者–食饵系统的平衡点的类型和稳定性 [3]。黄秀琴研究了一般形式的捕食–食饵模型的平衡点的类型和稳定性 [4]。本文受其启发,研究由广义Holling III型功能性反应 p ( x ) = m x 2 a x 2 + b x + 1 产生的Leslie型捕食者–食饵系统:

{ d x d t = r x ( 1 x K ) m x 2 y a x 2 + b x + 1 , d y d t = s y ( 1 y h x ) . (2)

p ( x ) 是捕食者的摄取率作为食饵密度的函数 [5]。其中, b > 2 a , a > 0 , m > 0 ,对所有的 x 0 ,使得 a x 2 + b x + 1 > 0 ,因此,对所有的 x > 0 p ( x ) > 0

对系统(2)作变换:

t ¯ = r t , x ¯ = x K , y ¯ = m K y r , a ¯ = a K 2 , b ¯ = b K , δ = s r , β = s h m K 2 . (3)

x ¯ , y ¯ , t ¯ , a ¯ , b ¯ 仍记为 x , y , t , a , b ,则系统(1)可以化为其等价系统:

{ d x d t = x ( 1 x ) x 2 y a x 2 + b x + 1 , d y d t = y ( δ β y x ) . (4)

2. 系统的平衡点分析

对于任意的参数,系统(4)总是存在一个边界平衡点 E 0 = ( 1 , 0 ) ,并且 E 0 是一个双曲鞍点。这个边界平衡点的生态学解释是:当缺乏捕食者时,食饵种群达到它的环境承载能力 [6]。平衡点 E 0 将x轴正轴分成两个部分,这两个部分是 E 0 的两个稳定流形,并且在G的内部,存在 E 0 的唯一不稳定流形。

不失一般性,假设点 E 1 = ( x ¯ , y ¯ ) 是系统(4)的任意正平衡点,则 x ¯ 是方程

a x ¯ 3 + ( δ β + b a ) x ¯ 2 + ( 1 b ) x ¯ 1 = 0 (5)

在区间 ( 0 , 1 ) 上的根。由于系统(4)的正平衡点个数由方程(5)在区间 ( 0 , 1 ) 上的根的个数决定,并且我们注意到,方程(5)在区间 ( 0 , 1 ) 上可能有一个、两个或三个正根,则系统(4)在区间 ( 0 , 1 ) 上可能有一个、两个或三个正平衡点。系统(4)在正平衡点 E 1 = ( x ¯ , y ¯ ) 处的雅可比矩阵为

J ( E 1 ) = [ 1 2 x ¯ ( b x ¯ 2 + 2 x ¯ ) y ¯ ( a x ¯ 2 + b x ¯ + 1 ) 2 x ¯ 2 a x ¯ 2 + b x ¯ + 1 β y ¯ 2 x ¯ 2 δ 2 β y ¯ x ¯ ] , (6)

并且雅可比矩阵的行列式和迹分别为

D e t ( J ( E 1 ) ) = ( δ 2 β y ¯ x ¯ ) ( 1 2 x ¯ δ ( b x ¯ 3 + 2 x ¯ 2 ) β ( a x ¯ 2 + b x ¯ + 1 ) 2 ) + β y ¯ 2 a x ¯ 2 + b x ¯ + 1 , (7)

T r ( J ( E 1 ) ) = 1 2 x ¯ ( b x ¯ 2 + 2 x ) y ¯ ( a x ¯ 2 + b x ¯ + 1 ) 2 + δ 2 β y ¯ x ¯ . (8)

D e t ( J ( E 1 ) ) = 0 时, E 1 ( x ¯ , y ¯ ) 是一个退化平衡点;当 D e t ( J ( E 1 ) ) 0 时,是一个基本的平衡点;当 D e t ( J ( E 1 ) ) < 0 时, E 1 ( x ¯ , y ¯ ) 是一个双曲鞍点。

3. 主要结论及证明

类似于参考文献 [7],考虑系统(4)的平衡点个数,我们能够得到下面的结论。

引理1:令 Δ 0 = ( δ β + b a ) 2 + 3 a ( b 1 ) Δ 1 = 4 Δ 0 3 + ( 27 a 2 + 9 a ( 1 b ) ( δ β + b a ) 2 ( δ β + b a ) 3 ) 2

1) 如果 Δ 1 > 0 ,则系统(4)有一个唯一的正平衡点 E = ( x , y ) ,它是一个基本的平衡点,并且是一个反鞍点;

2) 如果 Δ 1 = 0 ,并且

2.1) Δ 0 > 0 ,系统(4)有两个正平衡点:一个是基本的反鞍平衡点 E 2 ( x 2 , y 2 ) = ( β δ β δ + ( 1 δ ) 2 , δ 2 β δ + ( 1 δ ) 2 ) ,另一个是退化平衡点 E ( x , y ) = ( 1 δ , δ β ( 1 δ ) )

2.2) Δ 0 = 0 ,则系统(4)有唯一的正平衡点 E ( x , y ) = ( 3 1 b , δ β 3 1 b ) ,它是一个退化平衡点;

3) 如果 Δ 0 < 0 δ β a b 3 a ( 1 b ) 2 a < b < 1 ,则系统(4)有三个正平衡点 E 1 ( x 1 , y 1 ) E 2 ( x 2 , y 2 ) ,和 E 3 ( x 3 , y 3 ) ,它们都是基本的平衡点,并且 E 3 ( x 3 , y 3 ) 是一个鞍点。

引理1的证明:对于引理中的情形1),唯一的正平衡点的稳定性容易证明。本文只证明引理中的情形2.1)。

首先寻找一些参数,使得系统(4)有一个非双曲平衡点 E 2 ( x 2 , y 2 ) ,且 D e t ( J ( E 2 ) ) > 0 , T r ( J ( E 2 ) ) = 0 ,和一个退化平衡点 E ( x , y ) ,且 D e t ( J ( E ) ) = 0 , T r ( J ( E ) ) = 0

我们能够验证,如果 D e t ( J ( E ) ) = 0 , T r ( J ( E ) ) = 0 ,则

x = 1 δ , y = δ β ( 1 δ ) , a = ( δ 1 ) 2 + β δ β δ ( 1 δ ) 2 , b = ( δ 1 ) 3 2 β δ β δ ( 1 δ ) . (9)

其中, 0 < δ < 1 ,如果 a = ( δ 1 ) 2 + β δ β δ ( 1 δ ) 2 b = ( δ 1 ) 3 2 β δ β δ ( 1 δ ) 0 < δ < 1 ,则系统(4)可以简化为其等价系统。

{ d x d t = x ( 1 x ) x 2 y ( δ 1 ) 2 + β δ β δ ( 1 δ ) 2 x 2 + ( δ 1 ) 3 2 β δ β δ ( 1 δ ) x + 1 , d y d t = y ( δ β y x ) . (10)

解得系统(10)有两个正平衡点 E 2 ( x 2 , y 2 ) = ( β δ β δ + ( 1 δ ) 2 , δ 2 β δ + ( 1 δ ) 2 ) E ( x , y ) = ( 1 δ , δ β ( 1 δ ) )

定理1:如果 ( a , b , β ) = ( 1 + δ 2 ( 1 δ ) 3 , δ 2 + δ + 2 1 δ , ( 1 δ ) 3 δ 2 ( 1 + δ ) ) ,且 0 < δ < 1 ,则系统(4)有两个正平衡点:

a) E 2 ( x 2 , y 2 ) = ( 1 δ 1 + δ 2 , δ 3 ( 1 + δ ) ( 1 δ ) 2 ( 1 + δ 2 ) ) 是一个多重数为1的不稳定多重焦点;

b)是一个余维为2的尖点。

定理1的证明:如果 ( a , b , β ) = ( 1 + δ 2 ( 1 δ ) 3 , δ 2 + δ + 2 1 δ , ( 1 δ ) 3 δ 2 ( 1 + δ ) ) ,且 0 < δ < 1 ,则系统(4)化为等价系统。

{ d x d t = x ( 1 x ) x 2 y 1 + δ 2 ( 1 δ ) 3 x 2 δ 2 + δ + 2 1 δ x + 1 , d y d t = y ( δ ( 1 δ ) 3 δ 2 ( 1 + δ ) y x ) . (11)

系统(10)有两个正平衡点 E 2 ( x 2 , y 2 ) = ( 1 δ 1 + δ 2 , δ 3 ( 1 + δ ) ( 1 δ ) 2 ( 1 + δ 2 ) ) E ( x , y ) = ( 1 δ , δ 3 ( 1 + δ ) ( 1 δ ) 2 )

a) 令 u = x 1 δ 1 + δ 2 , v = y δ 3 ( 1 + δ ) ( 1 δ ) 2 ( 1 + δ 2 ) ,将点 E 2 ( x 2 , y 2 ) 转化为原点,系统(11)围绕原点的泰勒展开形式为

d u d t = δ u ( 1 δ ) 3 δ 2 ( 1 + δ 2 ) v + ( 1 + δ + δ 3 ) u 2 + ( 1 δ ) 4 δ 3 u v + ( 1 + δ ) ( 1 + δ 2 ) ( 1 + δ δ 2 + 3 δ 3 δ 4 + δ 5 ) ( 1 + δ ) δ 2 u 3 + ( 1 δ ) 2 ( 2 δ + 3 δ 2 δ 3 + δ 4 ) δ 3 u 2 v + Ο ( | u , v | 4 ) (12)

d v d t = δ 4 ( 1 + δ ) ( 1 δ ) 3 u δ v δ 4 ( 1 + δ ) ( 1 + δ 2 ) ( 1 δ ) 4 u 2 + 2 δ ( 1 + δ 2 ) 1 δ u v ( 1 δ ) 2 ( 1 + δ 2 ) δ 2 ( 1 + δ ) v 2 + δ 4 ( 1 + δ ) ( 1 + δ 2 ) 2 ( 1 δ ) 5 u 3 2 δ ( 1 + δ 2 ) 2 ( 1 δ ) 2 u 2 v + ( 1 δ ) ( 1 + δ 2 ) 2 δ 2 ( 1 + δ ) u v 2 + Ο ( | u , v | 4 ) (13)

对(12)和(13)作变换:

u = ( 1 δ ) 7 δ 5 ( 1 + δ 2 ) 2 ( 1 + δ 2 ) x + ( 1 δ ) 3 δ 3 ( 1 + δ ) y , v = y . (14)

则系统(12)和(13)可以化为其等价系统

{ d x d t = δ 3 ( 1 δ ) 1 + δ 2 y + f ( x , y ) , d y d t = δ 3 ( 1 δ ) 1 + δ 2 x + g ( x , y ) . (15)

其中,

f ( x , y ) = ( 1 δ ) 7 δ 5 ( 1 + δ ) 2 ( 1 + δ 2 ) ( 1 δ + δ 2 2 δ 3 + δ 4 ) δ 5 ( 1 δ + δ 4 + δ 5 ) x 2 + ( 1 δ ) 3 ( 3 + 2 δ δ 2 + 2 δ 3 ) δ 3 ( 1 + δ ) x y + δ 5 ( 1 + δ ) 2 ( 1 + δ 2 ) ( 1 δ ) 6 ( 2 + δ δ 2 + δ 3 ) ( 1 δ ) 7 δ 6 ( 1 + δ ) 2 y 2 + ( 1 δ ) 5 ( 1 δ + 2 δ 2 3 δ 3 + δ 5 δ 6 + δ 7 ) δ 7 ( 1 + δ ) 2 x 3 + δ 5 ( 1 + δ ) 2 ( 1 + δ 2 ) ( 1 δ ) ( 3 3 δ + 8 δ 2 8 δ 3 δ 4 + 6 δ 5 4 δ 6 + 3 δ 7 ) δ 10 ( 1 + δ ) 3 ( 1 δ ) 5 ( 3 + 3 δ 4 δ 2 + 14 δ 3 12 δ 4 + 14 δ 5 5 δ 6 + 3 δ 7 ) δ 8 ( 1 + δ ) x y 2 ( 1 δ ) 7 ( 1 δ ) ( 1 + δ 2 ) 2 ( 1 + δ 3 δ 2 + 4 δ 3 2 δ 4 + δ 5 ) δ 6 δ 5 ( 1 + δ ) 2 ( 1 + δ 2 ) y 3 + Ο ( | x , y | 4 ) (16)

g ( x , y ) = ( 1 δ ) 3 δ ( 1 + δ ) x 2 + ( 1 + δ ) 2 ( 1 δ ) 7 δ 5 ( 1 + δ ) 2 ( 1 + δ 2 ) δ 6 ( 1 + δ ) 3 x 3 + ( 1 δ ) 5 ( 1 + δ 2 ) δ 4 ( 1 + δ ) 2 x 2 y + Ο ( | x , y | 4 ) (17)

Liapunov常数 [8] [9] 为

Re c 1 = 1 16 { ( f x x x + f x y y + g x x y + g y y y ) + 1 ω [ f x y ( f x x + f y y ) g x y ( g x x + g y y ) f x x g x x + f y y g y y ] } | x = y = 0 = ( 1 + δ ) 2 ( 3 + 2 δ 8 δ 2 + δ 3 δ 5 + δ 6 ) 8 δ 8 ( 1 + δ ) > 0 (18)

其中, ω = δ 3 ( 1 δ ) 14 π ( 1 + δ 2 ) ,由于 0 < δ < 1 ,则 E 2 ( x 2 , y 2 ) = ( 1 δ 1 + δ 2 , δ 3 ( 1 + δ ) ( 1 δ ) 2 ( 1 + δ 2 ) ) 是一个多重数为1的不稳定多重焦点。

b) 作变换:

x 1 = x ( 1 δ ) x 2 = y δ 3 ( 1 + δ ) ( 1 δ ) 2 . (19)

则系统(15)简化为其等价系统

{ d x 1 d t = δ x 1 + ( 1 + δ ) 3 δ 2 ( 1 + δ ) x 2 2 + δ 1 + δ x 1 2 + ( 1 δ ) 3 δ 3 ( 1 + δ ) x 1 x 2 + Ο ( | x 1 , x 2 | 3 ) , d x 2 d t = δ 4 ( 1 + δ ) ( 1 δ ) 3 x 1 δ x 2 δ 4 ( 1 + δ ) ( 1 δ ) 4 x 1 2 + 2 δ 1 δ x 1 x 2 ( 1 δ ) 2 δ 2 ( 1 + δ ) x 2 2 + Ο ( | x 1 , x 2 | 3 ) . (20)

作变换:

y 1 = x 1 , y 2 = δ x 1 + ( 1 + δ ) 3 δ 2 ( 1 + δ ) x 2 . (21)

则系统(20)可以化为等价系统

{ d y 1 d t = y 2 1 1 + δ y 1 2 1 δ y 1 y 2 + Ο ( | y 1 , y 2 | 3 ) , d y 2 d t = δ 1 + δ y 1 2 y 1 y 2 + 1 1 δ y 2 2 + Ο ( | y 1 , y 2 | 3 ) . (22)

在原点 ( 0 , 0 ) 的一个小邻域内作变换:

z 1 = y 1 2 δ 1 2 δ ( 1 δ ) y 1 2 , z 2 = y 2 1 1 + δ y 1 2 1 1 δ y 1 y 2 . (23)

则系统(22)可以化为等价系统

{ d z 1 d t = z 2 + Ο ( | z 1 , z 2 | 3 ) , d z 2 d t = μ 1 z 1 2 + μ 2 z 1 z 2 + Ο ( | z 1 , z 2 | 3 ) . (24)

其中, μ 1 = δ 1 + δ , μ 2 = δ + 3 1 + δ ,因为 μ 1 μ 2 = δ ( δ + 3 ) ( 1 + δ ) 2 0 ,所以平衡点 E ( x , y ) 是一个余维为2的尖点。定理得证,从而引理中的2.1)得证,即 Δ 0 > 0 时,系统(4)有两个不同的正平衡点:一个是基本的反鞍平衡点 E 2 ( x 2 , y 2 ) = ( β δ β δ + ( 1 δ ) 2 , δ 2 β δ + ( 1 δ ) 2 ) ,另一个是退化平衡点 E ( x , y ) = ( 1 δ , δ β ( 1 δ ) )

致谢

感谢审稿老师及编辑老师提出的宝贵意见。

文章引用

潘陈蓉,陈松林. 具有广义Holling III型功能反应的Leslie捕食–食饵系统的平衡点分类和稳定性分析
Classification and Stability Analysis of Equilibrium Points in a Leslie Predator-Prey System with Generalized Holling Type III Functional Response[J]. 应用数学进展, 2019, 08(10): 1625-1631. https://doi.org/10.12677/AAM.2019.810191

参考文献

  1. 1. Huang, J.C., Ruan, S.G. and Song, J. (2014) Bifurcations in a Predator-Prey System of Leslie Type with Generalized Holling Type III Functional Response. Journal of Differential Equations, 257, 1721-1752.
    https://doi.org/10.1016/j.jde.2014.04.024

  2. 2. 卜令杰, 窦霁虹, 刘萌萌, 邢伟. 一类三次系统极限环的存在唯一性[J]. 延安大学学报(自然科学版), 2014, 33(2): 1-5.

  3. 3. Atabaigi, A. and Barati, A. (2017) Relaxation Oscillations and Canard Explosion in a Predator-Prey System of Holling and Leslie Types. Nonlinear Analysis: Real World Applications, 36, 139-153.
    https://doi.org/10.1016/j.nonrwa.2017.01.006

  4. 4. 黄秀琴. 一类Kolmogorov系统的极限环存在性与唯一性[J]. 南京晓庄学院学报, 2003(4): 63-66.

  5. 5. Dai, Y.F., Zhao, Y.L. and Sang, B. (2019) Four Limit Cycles in a Predator-Prey System of Leslie Type with Generalized Holling Type III Functional Response. Nonlinear Analysis: Real World Applications, 50, 218-239.
    https://doi.org/10.1016/j.nonrwa.2019.04.003

  6. 6. 林园, 高瑾. Lotka-Volterra竞争扩散系统连接边界平衡点和正平衡点行波解的存在性[J]. 教育教学论坛, 2019(27): 95-98.

  7. 7. Li, Y.L. and Xiao, D.M. (2007) Bifurcations of a Predator-Prey System of Holling and Leslie Types. Chaos, Solitons and Fractals: Applications in Science and Engineering: An Interdisciplinary Journal of Nonlinear Science, 34, 606-620.
    https://doi.org/10.1016/j.chaos.2006.03.068

  8. 8. Perko, L. (1996) Differential Equations and Dynamical Systems. Springer, New York.
    https://doi.org/10.1007/978-1-4684-0249-0

  9. 9. Guckenheimer, J. and Holmes, P.J. (1983) Nonlinear Oscillations, Dynamical Systems and Bifurcations of Vector Fields. Springer, New York.
    https://doi.org/10.1007/978-1-4612-1140-2

期刊菜单