Advances in Applied Mathematics
Vol.
09
No.
11
(
2020
), Article ID:
38832
,
10
pages
10.12677/AAM.2020.911238
一类具有恐怖因子和羊群效应的捕食者–食饵模型的动力学行为
吴耀冲,温洁嫦
广东工业大学应用数学学院,广东 广州
收稿日期:2020年11月2日;录用日期:2020年11月19日;发布日期:2020年11月26日
摘要
在生物种群中,由于食饵会对天敌的捕食产生恐惧,所以他们会降低自身出生率来抵抗捕食者的捕抓,针对这一情况研究了一类具有恐怖因子和羊群效应的捕食者–食饵模型,首先证明了解的一致有界性;然后由Routh-Hurwitz准则讨论了系统平衡点的存在性和存在时平衡点的类型,以及稳定性;并且根据Poincare’-Andronov-Hopf分支定理得到了系统产生Hopf分支的条件。最后,通过Matlab软件仿生所得到的定理,并得出了只要选取恰当的参数便可使系统中的两类种群一直稳定存在下去的结论。
关键词
捕食者–食饵模型,恐怖因子,Hopf分支,平衡点
Dynamic Behavior of a Predator-Prey Model with Fear Factor and Herding Effect
Yaochong Wu, Jiechang Wen
School of Applied Mathematics, Guangdong University of Technology, Guangzhou Guangdong
Received: Nov. 2nd, 2020; accepted: Nov. 19th, 2020; published: Nov. 26th, 2020
ABSTRACT
In the biological population, the prey will be afraid of the predator, so they will reduce their own natality to prevent being caught. Considering this situation, a predator-prey model with fear factor and herd effect is studied. Firstly, we discuss uniform boundedness of this model. Then Routh-Hurwitz criterion is used to discuss the existence of the equilibrium point of the system, the type of the equilibrium point and its stability. Moreover, according to the Poincare’-Andronov-Hopf branch theorem, conditions for the generation of Hopf branch were obtained. Finally, the theoretical results are verified by numerical simulation, and a conclusion is drawn that the predator and the prey in the system can be kept stable as long as appropriate parameters are chosen.
Keywords:Predator-Prey Model, Fear Factor, Hopf Bifurcation, Equilibrium Point
Copyright © 2020 by author(s) and Hans Publishers Inc.
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
1. 引言
食饵–捕食者模型在动力学系统中的研究越来越受关注,并得到了很多相关的结论。比如文献 [1] 在确定性和随机环境下,对一种具有牧群效应的捕食者食饵模型进行了分析,得出在随机环境下,捕食者的死亡率和食饵的出生率均会受高斯白噪音的干扰的结论;文献 [2] 研究了具有时滞的捕食者–食饵模型,分析了系统平衡点稳定性,且讨论了单时滞和双时滞下系统产生Hopf分支的条件;文献 [3] [4] 讨论了各种捕食者–食饵模型的Hopf分支情况;后来很多学者考虑到人类在生态系统中的影响,文献 [5] 研究了食饵具有常数收获率的捕食者–食饵模型,而文献 [6] 则考虑了捕食者和食饵皆具有线性收获率的生态模型。在2011年,Zanette et al.通过实验发现,由于对天敌的恐惧,哥雀减少了约40%的出生率 [7],之后各界学者
开始越来越关注恐惧因子对生态系统的影响。文献 [8] 研究了带恐惧因子 的单物种模型、捕
食者–食饵模型和生态流行病学模型的动力学性质和生物学意义。生物学上的研究表明,捕食者的存在会减少猎物的能量储备,因为猎物会对捕食者的捕食行为产生应激激素的分泌 [9] [10]。这种应激反应可能因此“吓死猎物”,进一步增加猎物的死亡率 [11] [12] [13] [14] [15] 更多关于由于对捕食者的恐惧而导致的食饵死亡增加的讨论可以在 [14] [15] [16] 中找到。
受文献 [8] 的启发,对具有羊群效应的捕食者–食饵模型 [17],我们引入了恐惧因子 对
食饵出生率的影响,新的模型如下:
(1)
其中x为食饵,y为捕食者,r为食饵的最大出生率,k为环境容纳量,f为恐惧水平,a为捕食者对食饵的捕捉能力, 为捕食者将食饵转化的能力,m为捕食者的自然死亡率。
为了研究方便,对系统(1)进行了如下无量纲变换
此时系统(1)可化为
(2)
其中 。
考虑到生物学的意义,系统(2)只考虑 的情况。
2. 解的有界性
定理1 在系统(1)中,所有从 出发的解都是一致有界的。
证明
(3)
由积分不等式理论得
(4)
故 有界。
定义函数 ,那么
(5)
因此,
(6)
由微分不等式得
(7)
所以,系统(1)中,所有从 出发的解都是一致有界的。
3. 模型的平衡点分析
下面由常微分方程定性理论 [18] 讨论系统(2)的平衡点存在条件和稳定性。
系统(2)具有三个平衡点,分别是平凡平衡点 ,边界平衡点 和非负平衡点 当且仅当 。其中
(8)
(9)
定理2 平衡点 为鞍点,平衡点不稳定。
证明:系统(2)在 点的雅可比矩阵是:
(10)
那么矩阵 的特征方程为
(11)
因为 ,此时 为不稳定的鞍点。
定理3 当 时,平衡点 为稳定的结点。当 时,平衡点 为鞍点。
证明:系统(2)在 点的雅可比矩阵是:
(12)
的特征方程为
(13)
所以当
(14)
(15)
显然 ,所以当 时,平衡点 为结点。
进一步地,此时 ,此时平衡点为稳定结点。而当 ,此时 ,故平衡点为鞍点。
定理4 当 且 ,此时平衡点 是不稳定的焦点;当 且 ,此时平衡点 是稳定的焦点。当 且 ,此时平衡点 是不稳定的结点;当
且 ,此时平衡点 是稳定的结点。
证明:系统(2)在 点处的雅可比矩阵为
(16)
其特征方程为 ,其中 ,,所以当 时,
有以下结论:
H1:当 且 ,此时平衡点 是不稳定的焦点;当 且 ,
此时平衡点 是稳定的焦点。
H2:当 且 ,此时平衡点 是不稳定的结点;当 且 ,
此时平衡点 是稳定的结点。
4. Hopf分支分析
定理5 当 时,系统(2)在 附近产生超临界的Hopf分支。
证明:系统(2)在 处的特征矩阵为
(17)
特征方程为 ,当 时,特征方程存在一对纯虚数特征值 λ = ± ω 0 ,此时 。取 作为分支参数。
令 ,,则 ,当 时,有 ,,易证 知系统(2)产生Hopf分支的横截条件 [19] 满足,所以当 时,系统(2)在 附近产生Hopf分支。下面计算Hopf分支的方向。
令 ,,则系统变为
(18)
再做变换
(19)
系统(3)再次变为
(20)
引入复变量: ,则
(21)
其中:
(22)
将以上结果通过一阶Lyapunov系数计算公式 [19]
(23)
可得
(24)
所以 ,故系统(2)在 附近产生超临界的Hopf分支。
5. 数值模拟
我们选取合适的参数值之后,对定理2进行了数值模拟。取 ,此时平衡点 是鞍点,如图1所示。在生物学意义上,若 ,那么随着时间的推移,无论捕食者或食饵种群取何种初始条件,随着时间的流逝,捕食者和食饵最终都不能共存,最终捕食者种群灭绝而食饵种群生存下去。
对定理3进行了数值模拟。取 得到平衡点 为稳定结点的轨线走向(图2),取 时,得到平衡点 为鞍点的轨线走向(图3)。在生物背景下,当 时,任意取初始值,捕食者种群最终走向灭绝。
Figure 1. A phase diagram near a saddle point
图1. 鞍点附近的轨线图
Figure 2. A phase diagram near a stable crunode E1
图2. E1稳定结点附近的轨线图
Figure 3. A phase diagram near a saddle point E1
图3. E1鞍点附近的轨线图
选取合适的参数值之后,对定理4进行了数值模拟。取 得到不稳定的焦点(图4),当 时,得到稳定的焦点(图5)。此时捕食者和食饵最终共存。取 得到不稳定的结点(图6),当 时,得到稳定的结点(图7)。
Figure 4. A phase diagram near an unstable focus E2
图4. E2是不稳定焦点的轨线图
Figure 5. A phase diagram near a stable focus E2
图5. E2是稳定焦点的轨线图
Figure 6. A phase diagram near an unstable crunode E2
图6. E2是不稳定结点的轨线图
Hopf分支数值模拟:取 ,得到极限环,如图8所示。此时两物种呈周
期性变化而共存。
Figure 7. A phase diagram near a stable crunode E2
图7. E2是稳定结点的轨线图
Figure 8. gets the limit cycle
图8. 时得到极限环
6. 结论
文献 [9] 提出了具有羊群效应的捕食者–食饵模型,在此基础上,本文考虑到食饵具有反捕食的行为从而引入恐惧因子。文中的定理1证明了模型的一致有界性,定理2~4分别研究了系统平衡点的存在性和稳
定性,并进行了数值模拟证明;定理5得出了当 时,系统产生Hopf分支,为超临界Hopf分支的
结论,并且发现通过Hopf分支可产生稳定的极限环,表明恰当选取参数,便可使系统中的被捕食者和捕食者最终共同存在且呈周期性变化而稳定下来。
时得到极限
Fig.8 gets the limit cycle.
态系统加以控制,便可使系统中的食饵和捕食者最终共同存在且可使两物种呈周期性变化而稳定下来。
文章引用
吴耀冲,温洁嫦. 一类具有恐怖因子和羊群效应的捕食者–食饵模型的动力学行为
Dynamic Behavior of a Predator-Prey Model with Fear Factor and Herding Effect[J]. 应用数学进展, 2020, 09(11): 2053-2062. https://doi.org/10.12677/AAM.2020.911238
参考文献
- 1. Bera, S.P., Maiti, A. and Samanta, G. (2016) Stochastic Analysis of a Prey-Predator Model with Herd Behaviour of Prey. Nonlinear Analysis: Modelling and Control, 21, 345-361. https://doi.org/10.15388/NA.2016.3.4
- 2. 李玉叶, 石惠文, 李旭超. 具有时滞的食饵–捕食者模型Hopf分岔分析[J]. 赤峰学院学报(自然版), 2018, 34(3): 1-4.
- 3. 邓田. 一类生物动力系统的稳定性与Hopf分岔研究[D]: [硕士学位论文]. 兰州: 兰州交通大学, 2016.
- 4. 王新秀, 窦霁虹, 彭煜. 一类具有时滞的捕食–食饵系统发生Hopf分支的条件[J]. 西北大学学报(自然科学版), 2010, 40(6): 957-960.
- 5. 邹桂华, 刘宣亮. 一类具有常数收获率的捕食者–食饵模型的分支分析[J]. 动力系统与控制, 2019, 8(3): 181-190.
- 6. 吴凡, 焦玉娟. 一类带有线性收获率的捕食者–食饵模型的稳定性分析[J]. 湖北民族学院学报(自然科学版), 2019, 37(3): 282-286.
- 7. Zanette, L.Y., White, A.F., Allen, M.C. and Clinchy, M. (2011) Perceived Predation Risk Reduces the Number of Offspring Songbirds Produce per Year. Science, 334, 1398-1401. https://doi.org/10.1126/science.1210908
- 8. Sasmal, S.K. (2018) Population Dynamics with Multiple Allee Effects Induced by Fear Factors—A Mathematical Study on Prey-Predator Interactions. Applied Mathematical Modelling, 64, 1-14. https://doi.org/10.1016/j.apm.2018.07.021
- 9. Stoks, R., Govaert, L., Pauwels, K., et al. (2016) Resurrecting Complexity: The Interplay of Plasticity and Rapid Evolution in the Multiple Trait Response to Strong Changes in Predation Pressure in the Water Flea Daphnia Magna. Ecology Letters, 19, 180-190. https://doi.org/10.1111/ele.12551
- 10. Trussell, G.C., Ewanchuk, P.J. and Matassa, C.M. (2006) The Fear of Being Eaten Reduces Energy Transfer in a Simple Food Chain. Ecology, 87, 2979-2984. https://doi.org/10.1890/0012-9658(2006)87[2979:TFOBER]2.0.CO;2
- 11. Fraker, M.E. (2009) Predation Risk Assessment by Green Frog (Rana clamitans) Tadpoles through Chemical Cues Produced by Multiple Prey. Behavioral Ecology & Sociobiology, 63, 1397-1402. https://doi.org/10.1007/s00265-009-0822-6
- 12. McCauley, S.J., Rowe, L. and Fortin, M.-J. (2011) The Deadly Effects of “Nonlethal” Predators. Ecology, 92, No. 11. https://doi.org/10.1890/11-0455.1
- 13. Preisser, E.L. (2010) The Physiology of Predator Stress in Free-Ranging Prey. Journal of Animal Ecology, 78, 1103-1105. https://doi.org/10.1111/j.1365-2656.2009.01602.x
- 14. Zhu, Z., Wu, R., Lai, L., et al. (2020) The Influence of Fear Effect to the Lotka-Volterra Predator-Prey System with Predator Has Other Food Resource. Advances in Difference Equations, 2020, Article No. 237. https://doi.org/10.1186/s13662-020-02612-1
- 15. Vesely, L., Boukal, D.S., Buřič, M., et al. (2017) Effects of Prey Density, Temperature and Predator Diversity on Nonconsumptive Predator-Driven Mortality in a Freshwater Food Web. Scientific Reports, 7, Article No. 18075. https://doi.org/10.1038/s41598-017-17998-4
- 16. Mukherjee, D. (2020) Role of Fear in Predator-Prey System with Intraspecific Competition. Mathematics and Computers in Simulation, 177, 263-275. https://doi.org/10.1016/j.matcom.2020.04.025
- 17. Ajrald, V., Pittavino, M. and Venturino, E. (2011) Modeling Herd Behavior in Population Systems. Nonlinear Analysis: Real World Applications, 12, 2319-2338. https://doi.org/10.1016/j.nonrwa.2011.02.002
- 18. 马知恩, 周义仓. 常微分方程定性与稳定性方法[M]. 北京: 科学出版社, 2001.
- 19. Kuznetsov, Y.A. (2004) Elements of Applied Bifurcation Theory. Applied Mathematical Sciences, 288, 715-730. https://doi.org/10.1007/978-1-4757-3978-7