尚荣忠

温州大学，浙江 温州

收稿日期：2017年3月7日；录用日期：2017年3月25日；发布日期：2017年3月28日

本文主要在动态建模过程中考虑易感染浮游植物吸收营养盐转化为自身能量所造成的时间滞后和易感浮游植物变为感染浮游植物所造成的时间滞后，进而建立了一类具有时滞效应的营养盐–浮游植物动态系统。在此基础上，主要采用Routh-Hurwitz判据和Lyapunov-Krasovskii函数对该系统内平衡点的稳定性和Hopf分支进行了理论分析，获得了系统具有这些特定动力学行为的临界条件，这些研究结果对预防和控制浮游植物大面积爆发具有一定的指导意义。

关键词 :时滞，平衡点，稳定性，Hopf分支

Copyright © 2017 by author and beplay安卓登录

This work is licensed under the Creative Commons Attribution International License (CC BY).

http://creativecommons.org/licenses/by/4.0/

浮游生物在调节海洋和气候的变化中扮演着重要的角色，因为它们是食物链中最基础的生物部分，是全球碳循环不可缺少的一部分 [ 1 ] 。然而，随着人类社会经济的发展，化学物质已经被应用到日常生活的工业生产中，随着这些化学物质进入河流进入大海，从而造成浮游植物持续增长并导致水体富营养化，最终形成蓝藻水华和海洋赤潮，使得生态系统失稳，进而对人类造成极大损害 [ 2 ] [ 3 ] 。

为预防与控制有害藻类水华现象，并维持浮游动植物种群的长久共存，保护水生生态环境，使开发生态环境与人类生存相和谐 [ 4 ] ，有很多科学家都对浮游生物的生长繁殖进行了研究。数学生态学家对这一现象进行了研究，并建立了多种不同的数学模型来表述它们之间的关系。Chatterjee Anal. [ 5 ] 建立N-P-Z模型并且将Holling II模型作为功能反应项，文章中讨论了系统的边界性和Hopf分支，且采用数值模拟对结果进行了拟合。Mukhopadhyay B. [ 6 ] 研究的也是N-P-Z模型，他们将一种IV功能反应函数应用于系统当中，对ODE系统的平衡点稳定性及分支行为进行了研究，之后将模型空间扩展分析扩散对稳定性的影响。Chafferjee Anal and Yu Hengguo [ 7 ] [ 8 ] 分析这一系统由于妊辰引起的时间滞对系统稳定性和Hopf分支的影响。除了对N-P-Z模型研究，N-P模型和P-Z模型也吸引了很多的科学家。Chakraborty Sabhendu [ 9 ] 研究了关于营养盐和有毒浮游植物的空间动力学系统，考虑了系统的线性稳定性和图灵失稳的状态。Wang Yapei [ 10 ] 基于数值分析，发现浮游植物密度随着营养盐密度的增长而增长，并且由于扩散和浮游植物的下沉会影响系统的稳定。N-P系统中考虑了时滞的影响 [ 11 ] [ 12 ] ，对平衡点的存在性和稳定性、Hopf分支和方向以及分支周期解都进行了分析。Dai Chuanjun [ 1 ] 对影响N-P系统的反应扩散和时滞进行了系统的分析研究。总而言之，许多数学生态学家对P-Z模型做了大量的研究工作 [ 13 ] [ 14 ] [ 15 ] [ 16 ] ，他们考虑不同的影响因素包括扩散、脉冲、时滞等和不同的功能反应函数结合进行研究以达到更接近实际环境的目的。

在水生生态系统中，浮游植物种群不仅相互竞争，而且产生的毒素会抑制对方生长 [ 17 ] 。本文主要研究了带有时滞因子的N-P系统，考虑了营养盐、易受感染的浮游植物和已经被感染的浮游植物三者之间的关系并建立数学模型。事实上，已经有很多科学家对这类系统进行了研究，Chattopadhyay J. [ 18 ] 研究了共存和不共存的一般性条件，例如考虑营养吸收函数和Holling II功能反应函数。他们考虑脉冲对系统的影响 [ 19 ] ，得到了系统灭绝的周期解、全局渐近稳定性和持续性存在的条件，揭示了脉冲扩散在系统中扮演的重要角色。Li Jianjun [ 20 ] 研究了关于N-P的反应扩散系统，得到了系统不恒稳定和存在分支的临界条件。

在这些模型的基础上，考虑由于消化和感染引起的时间滞后现象，进而建立具有时滞效应的营养盐、易感染浮游植物和被感染浮游植物动力学模型如下：

这个模型的建立具有以下假设前提：N、P₁、P₂分别代表在时间t时刻营养盐浓度、易感染和被感染的浮游植物的浓度。此外，所有参数取值均为正数，其中 表示外部流入系统的营养盐总量，e代表营养盐从水面到水的深层造成的流失率，a₁和a₂是浮游植物吸收营养盐后的转化率，k是半饱和常数，b₁和b₂是浮游植物死亡后转化为营养盐的转化率， 和 分别表示浮游植物被感染率和易感染者的转化率，m₁和m₂表示易感染者和感染者的死亡率， 代表易感染者吸收营养盐消化的时间， 代表易感染的浮游植物从被感染到变成感染者的时间。

同时，系统(1)的初始函数应该满足下面的条件：

其中：

这部分主要讨论系统(1.1)内平衡点的稳定性。内平衡点 代表三个种群共存，其中：

存在需满足下面的条件：

第一种情况：

定理2.1：当 时，满足条件(H2)时，平衡点 是局部渐近稳定的。

证明：在假设(H1)条件下，对系统(1.1)进行换元变形：

将内平衡点 转化到原点，系统(1.1)可以重新写为：

其中：

系统(2.1)对应的特征方程为：

其中：

根据Routh-Hurwitz准则，满足条件(H2)时， 是局部渐近稳定的。

证明完毕。

第二种情况：

引理2.1：当 ，满足条件(H3)，系统内平衡点是局部渐近稳定的。

证明：当 时，系统(1.1)可以写成：

我们选择Lyapunov-Krasovskii函数做变换得到如下：

对 进行求导并且等于系统(2.3)的迹，得到：

因此，得到：

,

其中：

所以我们得到结果：如果 ，系统(2.3)是局部渐近稳定的，所以系统内平衡点 是局部渐近稳定的。证毕

系统(2.3)的特征方程可以写成如下：

其中：

假设方程(2.4)的根 是纯虚根，将其带入方程(2.4)并且分离实部和虚部，得到下面的方程组(2.5)：

平方并且相加两个方程得到方程(2.6)：

其中： 。定义函数 。很明显如果 那么 。当 ,方程 有两个实根：

，

注意到 和

,

通过方程(2.6)讨论 的稳定性。

引理 2.2

i) 如果(2.6)没有正根，满足条件(H1)、(H2)， ， 是局部渐近稳定的。

ii) 如果(2.6)存在几个正根 ，可以从(2.5)得到 ，令 。当 时，特征方程存在一对纯虚根 ；

当 时，特征方程没有纯虚根。

在第二种情况下讨论Hopf分支。

分别对方程(2.4)两边关于 进行微分求导，令 为方程(2.4)的根，并且 ，将 代入(2.4)求导，可得：

所以

其中

定理 2.2当 并且满足条件(H1)和(H2)

i) 如果(2.6)没有正根，当 ，系统(1.1)的内平衡点 是局部渐近稳定的。

ii) 如果(2.6)有个不同的正根 ， ，系统是局部渐近稳定的；当 且 ，系统(1.1)在 产生Hopf分支。

第三种情况：

当 时，系统(1.1)可以写成如下形式：

我们将 放在其在第二种情况的稳定区间，将 看作变量进行研究

系统(1.1)的特征方程变为

其中

并且

假设 是方程(2.7)的纯虚根且将 代入方程(2.7)，两边同时乘以 ，可以得到：

其中

分离实部和虚部得到

对(2.8)进行处理得到

且

所以代入(2.9)并进行简化，得到关于 的函数

1) 如果(2.10)不存在正根，(2.8)没有关于 的正根。

2) 如果(2.10)存在几个不同的正根 ，我们就可以从(2.5)得到 ，取 ，当 ，特征方程存在一对纯虚根 ；当 ，特征方程没有纯虚根。

通过(2.8)，可以得到：

,

令 ，

定理2.3i) 当 且 时， 是局部渐近稳定的。

ii) 当 且 时，方程有一对纯虚根 。 ，

系统(1.1)在 处出现Hopf 分支。

本文建立一类具有时滞效应的营养盐-浮游植物动态系统，基础Routh-Hurwitz判据和Lyapunov-Krasovskii函数，详细讨论了系统正平衡点的局部渐近稳定性和Hopf分岔动力学行为，获得了系统在正平衡点发生Hopf分岔行为的临界条件，这些研究结果对预防和控制浮游植物大面积爆发具有一定的指导意义。

首先，我向我的导师赵敏教授表达我最诚挚的谢意！感谢老师对我的论文细致而又严谨的指导，感谢老师对我生活及学习上无微不至的关怀。

其次，感谢我们实验室的于恒国师兄和戴传军师兄，感谢你们一直以来在实验室对我问题的细致引导、鼓励和解答，使我对我们专业有了更深更细致的理解。

最后，感谢温州大学数学与信息科学院的各位老师和领导，感谢他们的教导和帮助。同时感谢赵敏导师国家自然科学基金面上项目(KZ1111021)的资助。

国家自然科学基金面上项目(KZ1111021)。

尚荣忠. 一类具有时滞效应的营养盐–浮游植物系统的动力学研究Dynamic Analysis of a Time Delayed Nutrient-Phytoplankton System[J]. 应用数学进展, 2017, 06(02): 165-173. http://dx.doi.org/10.12677/AAM.2017.62019