1. 引言
在海洋生态系统中,浮游生物处于所有水生生物组成的食物链中的最末端,且将矿物质营养转化为可供其它生物生存的原始生物原料。在实验性生物数学研究中,研究者已经开始对浮游植物和营养盐之间的关系和重要性进行了相关研究 [1] [2] [3] [4] [5] 。近年来,人们发现当大量的营养物质被投放到水体中,水体环境会发生变化,由于营养物质的大量供给,浮游植物会不断繁殖,导致了水体富营养化,使得水中大量鱼死亡、环境污染,甚至威胁到人类的健康 [6] - [15] 。因为,水体富营养化防治与预测问题已成为国内外重要的环境问题之一,迫切需要有效地解决方案。
许多研究者建立了数学模型来研究生物种群动力学问题,其目的是为了探讨如何通过数学方法精确地描述水体富营养化问题,进而提出相应的预测与防治策略 [16] 。论文 [17] 构建了一类水生生态系统中的氮-磷-浮游植物动力学新模型,通过Lyapunov function理论讨论系统平衡点的存在性和稳定性。论文 [18] 提出了一类具有self and cross扩散的有毒浮游植物-浮游动物动力学新模型,对其平衡点的局部存在性和全局稳定性进行系统研究。Smith等人 [19] 研究了一个具有状态脉冲反馈控制的半循环发酵模型,并利用Floquet理论得到了阶-1周期解的稳定判据。在论文 [20] 中,作者建立了一个具有脉冲反馈控制的浮游植物-浮游动物动力学新模型,利用后继函数研究了阶-1周期解和阶-2周期解的存在性。在论文 [21] 中,作者研究了系统在脉冲反馈控制下阶-1周期解的存在性,并通过微分方程几何理论和后继函数证明了阶-1周期解的存在性。在论文 [22] 中,作者研究了具有脉冲反馈控制的浮游植物-鱼类动力学新模型,并利用后继函数研究了阶-1周期解的存在性,并在极限环稳定性理论的基础上探讨了阶-1周期解的稳定性。
近年来,脉冲微分方程理论已得到了更深度地发展,许多学者在研究时都会考虑脉冲带来的影响作用 [23] [24] [25] 。实际上,状态脉冲微分方程在许多领域都有广泛的应用,且人类对自然资源的管理与开发都是离散的,而且还是瞬时完成的,从而在瞬时时刻也破坏了原有生态种群的固有状态,因此把这种瞬间扰动用脉冲微分方程来表示更加契合实际 [26] - [31] 。在论文 [32] 中,作者构建了一类恒化器下的状态脉冲反馈方程,并利用后继函数和庞加莱映射得到了系统阶-1周期解的存在性与稳定性。在论文 [33] 中,作者研究了一类乙醇抑制下的状态脉冲反馈控制的非线性模型,利用半连续动力系统证明了阶-1周期解和阶-2周期解的的存在性。
基于以上讨论结果,本论文考虑一类带有Michaelis-Menten功能反应项的两种群脉冲动力系统:
(1)
其中
和
分别代表营养盐的浓度和浮游植物在时刻
的种群密度。假设所有参数都为正数且有
。
,
。
表示流入到系统中的营养盐浓度;
代表营养盐的流失率;Michaelis-Menten功能反应函数可表示为
;其中
代表半饱和参数,
代表浮游植物对营养盐的摄入率;
代表浮游植物的死亡率。当浮游植物的数量在
时刻达到了阈值
,就会采取控制策略使得营养盐浓度和浮游植物种群密度分别变成
和
。
在系统(1)中,当不考虑脉冲效应时,计算得到系统的平衡点有:
;
,其中
和
。显然,
总是存在的;通过 [34] 可知,在满足
时,
是存在的。
2. 模型和准备
对于系统(1),脉冲集
,
脉冲映射:
,
相集:
.
定义2.1 [35] 若系统的轨迹始于点
交脉冲集
于点
,接着跳到点
,于是称点
为点
的后继点,表达式
称作点
的后继函数。
引理2.1 [35] 对于系统(1),假设
和
为相集
中的两点,它们的后继点分别为
和
。若满足条件
,则系统(1)存在一个阶1-周期解,其中
为其后继函数。
定义2.2 (轨道稳定性)
称作轨道稳定,如果满足:对于给定的
,存在
,使得对于系统(1)任意的解
都满足:
,于是当
为了讨论(1)的稳定性,定义向量场的两个横截面
。假设系统(1)的轨道
,以
为起点,在脉冲效应的作用下,轨道跳到
上的
,然后又回到
上的点
处。其中
由
和参数
决定。所以,可以定义
上庞加莱映射:
3. 正周期解的存在性与稳定性
情形
定理3.1 若满足条件
,系统(1)存在一个唯一的阶-1周期解且是轨道渐近稳定。
证明:若条件
成立,则可得到关系:
。从而得到直线
位于
和
。
等倾线
与两直线
和
分别交于
和
。从
作一条垂直于
的直线,并交于
。
为
的后继点,
为
的后继点。根据轨迹的趋势可知,
必定位于
的右边,
必定位于
的左边。固有
和
,后继函数表达式为:
,
,因此必定存在一点
满足
,由引理2.1知:系统(1)存在一个位于
和
之间的阶-1周期解。
接下来证明阶-1周期解的唯一性和轨道渐近稳定性。
根据上面分析知:轨迹从点
出发,交脉冲集
于点
,接着跳到后继点
。根据轨迹的趋势,点
必定会位于
的右方;因此有:
。
从
出发的轨迹交脉冲集
于点
,接着跳到后继点
,而
必定位于
的右边,因此有
。因此可得两个序列:
类似地,从
开始的轨迹,也可得到两个序列:
根据系统(1)的前两个方程知:线段
映射到线段
。由于脉冲效应,线段
映射到线段
。因而知道
根据以上分析过程可知:
利用数学分析中的闭区间套定理知,系统(1)存在唯一的阶-1周期解,且满足:
在
足够小的邻域内,选取任意一点
(不同于
),不失一般性,假设
(否则
,证明方法同理)。
根据序列
的单调性,必存在整数
使得:
。
开始于
的轨迹与脉冲集
相遇后跳到点
(相互独立的轨迹不会相交),因此
必定位于
与
之间,且有
。
的后继点表示为
,
必定位于
与
之间,且满足
。根据以上分析过程,得到序列:
满足以下条件:
1)
位于
与
之间;
2)
。
因此,
。
综上,由于连续脉冲效应,相点的后继点都被吸引到
周围,而由
的任意性可知,系统(1)的阶-1周期解为轨道渐近稳定。
完证。
情形2
考虑到系统相空间的几何构建,有一个轨道
交
于
,并与
交于
,其中
。所以对于系统(1)始于
的轨道不会与
相交。
假若系统(1)有一个轨道从
由于脉冲效应会与
相交无数次,其中
。假设
与
交于
。由于脉冲效应可知:
,
。重复以上的过程得到,
。另一方面,对于
上的任意两点
和
,
和
满足:
且
。由于脉冲效应,点
,
均在
的右侧。因此,可得
(2)
假设轨道始于
,与
第一次相交于
,
,若
和
,则系统(1)有正阶-1周期解。且对于任意从
出发的轨道都不会与
相交。因此,从
的庞加莱映射,得出
,
。重复上面的过程,得到:
。特别地,若
,则系统(1)有正阶-1周期解;若
,且
,则系统(1)有阶-2周期解。
(a) 若
,由(2)得
,在这种情形下,
的关系有如下两种情况:
(i)
若
,则有
。由(2)重复以上的过程得到:
.
(ii)
若
,类比(i),得到:
.
(b) 若
,由(2)可得
,这种情形下,
的关系有如下两种:
(i)
若
,则有
,由(2),并重复以上过程得到:
.
(ii)
若
,类比(i)得到:
.
因此,当出现情形(a)中的(i)时可知:
,
,且满足:
。因此,
可得
和
,也易得到,系统(1)有一个阶-2轨道渐近稳定的正周期解。类似地,当出现(a)中的(ii)和(b)中的(ii)这两种情形时,系统(1)有一个阶-1轨道渐近稳定周期解;当出现(b)中的(i)时,系统(1)有一个阶-2轨道渐近稳定周期解。
4. 结论
本文建立了一类状态脉冲反馈控制下的营养盐–浮游植物的动力系统,对该系统进行了相关定性分析,探讨了该系统阶-1正周期解的存在性、唯一性和轨道渐近稳定性。这些研究工作为全面研究状态脉冲反馈控制理论在水体富营养化防治中的应用提供理论支撑。
致 谢
首先,我向我的导师赵敏教授表达我最诚挚的谢意!感谢老师对我的论文细致而又严谨的指导,感谢老师对我生活及学习上无微不至的关怀。
其次,感谢我们实验室的于恒国师兄和戴传军师兄,感谢你们一直以来在实验室对我问题的细致引导、鼓励和解答,使我对我们专业有了更深更细致的理解。
最后,感谢温州大学数学与信息科学院的各位老师和领导,感谢他们的教导和帮助。同时感谢赵敏导师国家自然科学基金面上项目(KZ1111021)的资助。