双时滞四维捕食网络的分析【基于matlab的动力学模型学习笔记_6】

/*仅当作学习笔记，若有纰漏欢迎友好交流指正，此外若能提供一点帮助将会十分荣幸*/

本系列谈论过单时滞，但还没提及过双时滞，本文将着重介绍一种双时滞系统并对其进行简单处理分析。

摘 要：本文针对一个捕食网络，根据其特性建立起一个双时滞四维系统。首先分析其稳定性，找到其正平衡点，然后通过调试时滞量τ来模拟预测网络模型的发展前景。并根据所求得临界滞量τ0，带入实值仿真出其Hopf分岔图形。

0引言

自然界中，广泛存在着种群之间的捕食与被捕食关系，他们一起构成了生物链，为流传的“大鱼吃小鱼，小鱼吃虾米”就通俗易懂的道出了生物链之间的关系。而为了保护生态链的完整，甚至保证整个生态系统的有序发展，对其中捕食系统的研究就存在相当的意义。

在捕食系统中，最经典的就是Lotka-Volterra捕食模型(简称L-V模型)，其模型方程为：

-V模型中，x1(t)、x2(t)分别表示再t时刻食饵和捕食者的密度，而x1、x2则分别表示食饵和捕食者在t时刻的种群密度变化量。a10代表食饵种群内自然增长率、a12代表捕食者对食饵捕食的所造成的影响（或者说威胁）、a20代表捕食者种群内部因素带来的数量变化，比如生老病死亦或同族之间的相互争斗带来的影响、a21代表食饵对捕食者所提供的积极影响，即提供食物来源[2]。

显而易见，L-V模型相当的粗略，考虑的因素也很少。因此，在此基础上本文通过增加模型的维数，来加长食物链的长度；通过加入更多系数来模拟捕食者捕食时间、捕食者成长期等因素的影响。试图通过考虑到更多的因素来提高模型的真实性，然后可以人为的调整控制因素来使生物链达到一个动态平衡，即达到一个良性。

1.建立系统模型

自然界中，广泛存在着种群之间的捕食与被捕食关系，他们一起构成了生物链，为流传的“大鱼吃小鱼，小鱼吃虾米”就通俗易懂的道出了生物链之间的关系。而为了保护生态链的完整，甚至保证整个生态系统的有序发展，对其中捕食系统的研究就存在相当的意义。

根据实际情况，并借鉴文献[1]我们这里设置的捕食模型应该具备如下前提条件：

 该双时滞的四维捕食模型为：

在该模型中，x1(t)、x2(t)、x3(t)、x4(t)分别代表一二三级食饵，和食物链顶端的捕食者。其中，一级食饵为食物链最低端，会被二级食饵捕食，二三级食饵以此类推（即生物链的概念，但捕食者只捕食比自己低一级的食饵，这里不考虑跨级捕食）。

而dx1/dt、dx2/dt、dx3/dt、dx4/dt分别代表一二三级食饵，和食物链顶端的捕食者在时刻的数量变化情况。

1.1参数设置

①对一级食饵，即x1(t)而言：

 ②对二级食饵，即x2(t)而言：