1. 背景
对红松鼠的保护刻不容缓 [1],利用数学模型研究红松鼠保护引起了学者的兴趣。近年来,灰松鼠种群增长比较快,而红松鼠受灰松鼠身上的疱疹病毒传染,红松鼠死亡率增高,而病毒对灰松鼠伤害较小 [2],因此灰松鼠种群增长比较快。英国本土红松鼠因外来灰松鼠的入侵而面临灭绝危机 [3],为了保护红松鼠的生存量,英国政府提出捕杀灰松鼠的方法限制灰松鼠的增长,这样就产生了一个问题:我们每年要捕杀多少灰松鼠,才能既保护了英国本土红松鼠达到人们所希望的数量而灰松鼠也能保留人们所希望的种群规模 [4]。
本文考虑疱疹病毒建立状态脉冲模型 [5]。
2. 模型建立
由于疱疹病毒对红松鼠影响较大,对灰松鼠影响较小,我们首先建立非Lota-Volterra竞争模型 [6]:
(1)
其中,x和y分别表示红松鼠和灰松鼠在t时刻的密度,系统(2)中微分方程组即为系统(1)两种群相互作用模型,
和
分别表示红松鼠和灰松鼠的环境容纳量,
表示每个灰松鼠占用红松鼠环境容纳量的数量,
表示每个红松鼠占用灰松鼠环境容纳量的数量,r表示每个灰松鼠身上的病毒使红松鼠种群产生的死亡率。
表示允许灰松鼠保存的密度,也就是说在这个密度以内不会危及红松鼠的生存,
表示灰松鼠生存的密度达到
时,人们开始捕杀灰松鼠时的捕杀率,系统(1)中的:
,
,
表示在捕杀灰松鼠时不要伤害红松鼠 [7]。
3. 无脉冲模型分析
模型(1)对应的无脉冲模型为:
(2)
其中,x和y分别表示红松鼠和灰松鼠在t时刻的密度,
和
分别表示红松鼠和灰松鼠的环境容纳量,
表示每个灰松鼠占用红松鼠环境容纳量的数量,
表示每个红松鼠占用灰松鼠环境容纳量的数量,r表示每个灰松鼠身上的病毒使红松鼠种群产生的死亡率。
表示灰松鼠t时刻占用红松鼠环境容纳量的数量。
3.1. 模型的向量场与一致有界性
定理1:系统(2)所有界最终一致有界。
证明:根据
,
得到水平等倾线为:
,为直线。
得到垂直等倾线为:
,为抛物线。
设直线
与等倾线交于一点A,G与x轴交于
,过A作直线L,设L与y轴交于点B,直线L
记为 [8]:见图1。
直线段
,
,
和
,围成一个闭区域G,见图2。
这说明轨线与线段L相交 [9],必进入区域G。
因此,我们有:系统(2)所有界最终一致有界。
3.2. 系统的平衡点分析
系统(2)存在四个平衡点:
,
,
,
系统(2)对应的雅各比矩阵:
由于
得到方程:
其中:
1) 当
,
,
,
M1:
,
M2:
,
,
M3:
,
,
M4:
系统(2)有一个正平衡点
。
,
此时系统的稳定性:
结论一:
为不稳定的结点。
证明:
,
,
,
由
可知说明平衡点是结点。
所以
是不稳定的结点。
结论二:
为鞍点。
证明:平衡点,
考虑M3成立
,
,
所以
是鞍点。
结论三:
为鞍点。
证明:平衡点
,
,
,
所以
是鞍点。
结论四:
为稳定的结点。
证明:平衡点
,
,
考虑
成立
所以
是稳定的结点。
3.3. 系统极限环不存在性
定理2:系统(1)在第一象限不存在极限环。
证明:我们用DuLac极限环的存在性定理,记
(2)
因此,由DuLac定理得到结论,系统(1)在第一象限不存在极限环 [10]。
定理3:平衡点
为全局渐近稳定。
证明:由结论4可知
为稳定的焦点,而且
是区域G中唯一的平衡点。又由定理2:系统(2)在第一象限不存在极限环。再由定理1:系统(1)所有界最终一致有界 [11]。
由此,我们就可以得出结论:平衡点
为全局渐进稳定 [12]。
3.4. 全局稳定性
通过构造V函数证明无控制系统的全局稳定性 [13]。我们构造V函数:
考虑到平衡点方程:
当
时,
因为
在第一象限内是定负的,所以在
平衡点处是全局渐近稳定的 [14]。
4. 脉冲模型分析
4.1. 阶一周期解的存在性
定理4:系统(1)存在阶一周期解。
证明:脉冲集M:
,相集N:
,直线
与相集交于一点A,过A点的轨线与脉冲集M交于一点
,由向量场可以看出轨线
的形状如图3所示,脉冲到相集上一点
,这说明
为A的后继点 [15],同样可知从原点出来的一条轨线
,有
为B的后继点,在相集上建立坐标系,以N与Y轴的交点为坐标原点,分别设点A的坐标为a,点B的坐标为b,分别设点
的坐标为
,
点的坐标为
,后继函数 [16]:
,
这样,由文 [17] 的后继函数判别法证明了定理4,即脉冲状态反馈控制系统(1)存在阶一周期解
。如图4。
4.2. 阶一周期解的稳定性
定理5:系统(1)的阶一周期解是稳定的。
证明:由定理2的证明知系统(2)中DuLac函数
变换后得系统(3)。
(3)
由定理2的证明中的计算有:
又系统(1)的轨线与系统(3)的轨线是完全相同的,由文献 [17] 即证明了阶一周期解
是稳定的。