1. 引言
众所周知，捕食者种群与食饵种群之间的捕食作用关系是生态种群之间最基本作用关系之一，作用过程具有复杂动态特性和随机不确定性 [1] 。此外，捕食者种群与食饵种群捕食作用关系一直并将继续在生态学和生物数学中占据主导地位 [2] [3] 。如何刻化这类捕食动态过程至关重要，目前大多数研究专家利用捕食者种群的功能反应函数来描述 [4] [5] [6] [7] 。1969年，Hassell和Varley [8] 首先发现捕食者种群之间的干扰行为可以影响捕食效率的实验证据，并提出了现在被称作Hassell-Varley型功能反应函数，这类功能反应函数可以表示如下：
.
其中
和
分别代表食饵种群和捕食者种群密度，
表示捕食者种群搜索食饵种群的效率，
表示食饵种群生物量转化为捕食者种群生物量的转化率，
称为Hassell-Varley常数 [8] ，用来表示捕食者种群之间的干扰。若
，即捕食者种群没有形成群体，则该模型就变为比率型捕食者–食饵系统；若
，即捕食者种群之间没有干扰，则该模型就变为Michaelis-Menten型(或是Holling typeⅡ型)捕食者–食饵系统。文献 [9] [10] [11] 已经研究了几类带有Hassell-Varley型功能反应的捕食-食饵系统，获得了一些比较好的结果。
自然生物种群受到周期性变化的自然环境的影响，就会本能地改变自己的生活习惯，例如季节迁徙现象和生物的不同季节的繁殖现象，因此研究受到周期性扰动的生态系统的动力学性态问题是十分重要的 [12] 。同时，在现实自然界中，生物种群自身的生长即不是连续的也不是离散的，而是两种并行存在的 [13] [14] [15] 。此外人类对自然资源的管理与开发利用都是离散的，而且还是瞬间完成的，从而在瞬间时刻也破坏了原有生态种群的固有状态，因此把这种瞬间扰动用脉冲微分方程来表示生物系统的数学模型更加符合实际 [16] - [21] 。文献 [22] [23] 详细介绍了脉冲微分方程理论以及它在现实其它学科中的运用，为后期做研究提供了宝贵的数学基础与依据。文献 [24] [25] [26] [27] [28] 详细研究了以时间作为实施脉冲控制取得的成果，证明了带有脉冲扰动的生物系统的半平凡周期解的存在性、局部渐近稳定性和全局渐近稳定性以及有界性和保持生态种群持久生存的关键条件；同时通过数值模拟显示系统所具有的特定动力学性态。
根据以上讨论，在这篇文章中考虑带有Hassell-Varley功能反应和季节效应的三种群脉冲动力学系统：
(1.1)
其中
和
分别表示食饵、中级捕食者和高级捕食者在时刻
的种群密度。
。假设所有参数都是正数。
表示捕食者搜索食饵的效率；
表示食饵生物量转化为捕食者生物量的转化率；
表示环境最大容纳量；
表示半饱和常数；
表示高级捕食者的死亡率；
表示脉冲效应周期；
表示在
时刻，高级捕食者的释放量或储存量；
表示在固定时刻
的收获率；
和
分别表示内禀增长率
与
的平均值；
表示季节扰动的角频率；参数
表示季节周期性振荡的强度。
和
是内禀增长函数的扰动最大量，参数
是平面三角度，其中
。当
时，表示食饵
与中级捕食者
的内禀增长率是同步受季节扰动的；当
时，表示食饵
与中级捕食者
的内禀增长率是异步受季节扰动的；当
时，表示食饵
与中级捕食者
的内禀增长率是反同步受季节扰动的。
首先假设
代入公式(1.1)，能够得到下面的辅助系统：
(1.2)
假设
代入公式(1.1)，能够得到下面的辅助系统：
(1.3)
这个生态系统表示食饵种群、中级捕食者种群和高级捕食者种群生活在一个生活格局中，高级捕食者种群主要依靠捕食食饵种群和中级捕食者种群而获得生存的物质能力，中级捕食者种群也捕食食饵种群获得生存的物质能量，但是中级捕食者种群也会从其他地方获得生存的物质能量。因为高级捕食者种群密度过高与过低都会影响其他两类种群密度的变化，为了维持三种群持久生存，在固定时间对高级捕食者种群进行一定量捕杀，减少成年高级捕食者种群的数量，以至于促进自然资源更好的繁衍生息。同时释放一定量的幼年高级捕食者种群，以防止食饵种群和中级捕食者种群繁衍太快，产生内部竞争过大、或者其他自然资源走向灭绝。
2. 基础知识与主要引理
定义2.1：设
和
。定义
是所有非负数的集合，
表示系统(1.1)的一个映射。设
，那么
属于
类，其中
满足以下条件：
1) 对
，
存在，那么
在
上是连续的。
2)
在
上满足局部Lipschitzian。
定义2.2 [19] [26] [30] [32] [33] 设
，那么
，定义系统(1.1)右导数
为
.
系统(1.1)的任意解
在
上是连续的，并且
是一定存在的，那么系统(1.1)的任意解是连续分段光滑的。根据文献 [29] - [31] ，映射
的光滑性可以保证系统(1.1)的解是唯一存在的。
定义2.3：如果存在一个紧集
，系统(1.1)的任意解
都最终进入或者驻留在紧集
内，那么称系统(1.1)是持久存在的。
引理2.1如果
是系统(1.1)的任意解，其初始值为
，那么
，就有
。进一步讲，如果
，就有
。
引理2.2(Lakshmikantham, 1989 [29] )。假设
并且
(2.1)
其中
在
上是连续的，且对
，
是存在的，
是单调不减函数。设
是脉冲微分方程
(2.2)
在
上的最大解，那么
，由此可得
。其中
是系统(1.1)的任意解。
如果引理中的不等式条件反过来，那么可得出相类似的结果。特别注意的是如果
的条件是光滑的，那么系统(2.2)的解是存在且唯一的，那么
是系统(2.2)的唯一解。
现在假设种群
和种群
是灭绝的，则系统(1.1)就变成子系统(2.3)。
(2.3)
可求得子系统(2.3)的周期解
，其中初始值为
。进一步可求得子系统(2.3)任意解与周期解之间的关系，即为
，其中初始值为
，
。
引理2.3 设
和
分别是子系统(2.3)的正周期解与任意解，其中初始值为
，则当
时，有
。
由引理2.3得知，当
充分大时，任意解
趋近于子系统(2.3)的正周期解
。因而系统(1.1)的半平凡周期解可写为
.
为了研究半平凡周期解的稳定性，给出了具有周期
线性脉冲方程的Floquet定理：
(2.4)
满足以下条件(H)：
(H1)
与
，其中
是在
时刻所有左连续的分段光滑连续矩阵函数的集合，
是所有
矩阵的集合；
(H2)
，
；
(H3) 存在某个数
使得
。
设
为方程(2.4)的基本矩阵，那么存在唯一一个非奇异矩阵
使得
. (2.5)
对于等式(2.5)，其相对应的基本矩阵
和常数矩阵或是系统(2.4)的单值矩阵
。
系统(2.4)的所有单值矩阵都是相似的并且有相同的特征值。单值矩阵的特征值
称之为系统(2.4) Floquet因子。
引理2.4 (Floquet定理 [26] [31] [32] [33] ) 设条件(H1)~(H3)成立。那么具有周期
的线性脉冲方程(2.4)是
1) 稳定的当且仅当系统(2.4)的所有因子
of (2.4)满足不等式
，而且，对于这些
，
，那么对应的是简单的初等因子；
2) 渐近稳定的当且仅当系统(2.4)的所有因子
满足不等式
；
3) 不稳定的，如果对于某些
，不等式
。
3. 数学分析
现在研究半平凡周期解
的局部渐近稳定性和全局渐近稳定性。
定理3.1 设
是系统(1.1)的任意解，如果
与
成立，那么，系统(1.1)的半平凡周期解
是局部渐近稳定的。
证明：首先，为了证明周期解
的局部渐近稳定性，考虑辅助系统(1.2)，那么，根据引理2.2，可得
和
，其中
与
分别是系统(1.1)和(1.2)的任意解。因此只需证明系统(1.2)的周期解
局部渐近稳定即可，其中
。考虑任意解的微扰为
. (2.6)
结合公式(2.6)和公式(1.2)，通过线性化表示，系统变为
(2.7)
因此
其中
满足
，
是单位矩阵并且
.
则半平凡周期解
的稳定性是由矩阵
的特征值决定。
如果所有
特征值的绝对值不超过1，那么半平凡周期解
是局部渐近稳定的。因此，对于
的所有特征值，即
根据引理2.4，如果
与
，
是局部渐近稳定的。即
与
.
定理证明完毕。
下面证明系统(1.1)任意解的一致有界性。
定理3.2 对于足够大的
，一定存在一个常数
使得对系统(1.1)的任意解
，都有
。
证明：设
，很显然
。计算
的右导数，可得
(2.8)
设
，那么
是有界的。选择适当的
与
使得
其中
是两个正数。
由引理2.2，易得
,
其中
。因此
所以
是一致有界的。这个结果也暗含了系统(1.1)的任意解也是有界的。
定理3.3 设
是系统(1.1)的任意解，如果
,
与
成立，那么半平凡周期解
是全局渐近稳定性的。
证明 根据定理3.1，易得
是局部渐近稳定的。下面只需证明
是全局吸引即可。设
，则有
根据定理3.2，对于足够大的时间
，存在一个常数
使得对系统(1.1)的任意解
都有
。因此
.
易得
(2.9)
由引理2.2与2.3知道，存在某一时刻
和充分小的
，当
，使得
。
设
，
那么
,
如果
与
，即
.
因此
，则有
.
因此当
时，
，
，
。因此半平凡周期解
是全局吸引子。从而半平凡周期解
是全局渐近稳定的。
推论1 设
是系统(1.2)的任意解，如果
,
与
成立，那么半平凡周期解
是全局渐近稳定的，其中
。
推论2 设
是系统(1.2)的任意解，如果
,
与
成立，半平凡周期解
是全局渐近稳定的，其中
。
现在研究系统(1.1)的持久生存性。
定理3.4 如果
,
与
成立，那么系统(1.1)是持久生存的。
证明：为了完成定理的证明，首先考虑辅助系统(1.3)，根据引理2.2，易得
,
.
.
其中
和
分别是系统(1.1)和(1.3)的任意解。
假设
的系统(1.3)的任意解，其初始值为
。由定理3.2，假设对
，有
.
由引理2.2，存在参数
和充分小的
使得
.
那么对于充分大的时间
，有
.
因此，对于充分大的时间
，只需找到
与
使得
与
即可。
选择
和充分小的
使得
，并记
,
其中
.
现在证明对任意的时间
，
不成立。否则，可以得到
.
根据引理2.2与2.1，当
时，易知
，
，其中
是
(2.10)
的解。因此，存在一个
使得
且
. (2.11)
设
及
并在
上积分(2.11)，易得
(2.12)
那么，当
时，得到
，这与
有界相矛盾。因此，存在一个
使得
。
第二，如果
，有
，那么就得到上面的结论。因此，只需考虑这些解不落入区域
。设
，易得
与
。由于
是连续的，很容易就能证明
。
可以断言一定存在一个
使得
，否则，
.
选择适当的
使得
及
，其中
.
考虑(2.10)及
，
，
，
那么
, 
,
意味着
。(2.11)是成立的。在
上积分(2.11)，可以得到
.
对于
，考虑以下两种情形：
情形Ⅰ:如果对
，
，那么
，则有
。
因此
. (2.13)
在
上积分公式(2.13)，就有
.
那么
,
这与已知相矛盾。
设
，则
，
及公式(2.13)成立。在
上积分公式(2.13)，则有
.
当
时，由于
，同理可证，
，
成立。
情形Ⅱ:存在一个
使得
。设
，则对于
，
及
。对于
，则公式(2.13)成立。在
上积分公式(2.13)，得到
。由于
，同理可证，可以得到当
时，
。
因此在这两种情形下，可以得出结论
，
。
类似的，可以证明
，
。
设集合
。
很显然，可以得出集合
是全局吸引子。系统(1.1)的任意解最终进入并保持在区域
内。因此，系统(1.1)是持久生存的。证明完毕。
4. 结论
本论文基于生态学理论与数学生物学知识，在动态建模过程中加入了Hassell-Varley功能反应函数，建立了一类具有季节效应的脉冲控制生物动力系统。基于脉冲微分方程理论、比较定理和生物学意义，研究了季节效应如何影响这类生态动力系统的特定动力学性态，分析了系统半平凡周期解的存在性、局部渐近稳定性和全局渐近稳定性，明确了该系统具有这些特定动力学性态的阈值条件，同时讨论了系统生物种群的灭绝性与持久生存性。这些研究结果为进一步研究如何运用脉冲控制策略维持生态种群持久生存提供了一定的理论支撑。
致谢
首先，我向我的导师赵敏教授表达我最诚挚的谢意！感谢老师对我的论文细致而又严谨的指导，感谢老师对我生活及学习上无微不至的关怀。
其次，感谢我们实验室的于恒国师兄和戴传军师兄，感谢你们一直以来在实验室对我问题的细致引导、鼓励和解答，使我对我们专业有了更深更细致的理解。
最后，感谢温州大学数学与信息科学院的各位老师和领导，感谢他们的教导和帮助。同时感谢赵敏导师国家自然科学基金面上项目(KZ1111021)的资助。
基金项目
国家自然科学基金面上项目(KZ1111021)。