Advances in Applied Mathematics
Vol.06 No.07(2017), Article ID:22378,15 pages
10.12677/AAM.2017.67099

Dynamic Analysis of a Dissolved Oxygen-Plankton Model with Two Time Delays

Qing Guo

Wenzhou University, Wenzhou Zhejiang

Received: Sep. 30th, 2017; accepted: Oct. 15th, 2017; published: Oct. 19th, 2017

ABSTRACT

Based on the internal dynamic characteristics of marine ecosystem and the dynamic mechanism of marine biological population, this paper constructed a dissolved oxygen-plankton dynamic model involving two time delays, and made theoretical analysis of corresponding dynamics. We obtained the locally asymptotic stability of positive equilibrium and the threshold conditions of occurring Hopf bifurcation were gained. Therefore, we analyzed the dynamics of Hopf bifurcation in detail. This result provides a great help for the interaction between dissolved oxygen and marine plankton in dynamics, and is helpful to deeply understand how the delay affects the dynamic trend of the marine ecosystem, and furthermore provide certain theoretical support for the study of dynamic growth mechanism of mutual restriction and mutual coordination in marine plankton.

Keywords:Dissolved Oxygen, Plankton, Two Time Delays, Hopf Bifurcation, Stability

一类具有双时滞效应的溶解氧-浮游生物模型的动力学分析

郭庆

温州大学,浙江 温州

收稿日期:2017年9月30日;录用日期:2017年10月15日;发布日期:2017年10月19日

摘 要

基于海洋生态系统内在动态变化特性和海洋生物种群之间的动态作用机制,本论文构建了一类具有双时滞效应的溶解氧-浮游生物动态模型,并对其相关动力学性质进行理论分析,解析出模型内平衡点具有局部渐近稳定性和发生Hopf分支的阈值条件,并详细探析了Hopf分支的相关动力学性质。本研究成果有利于从动力学的角度揭示溶解氧与海洋浮游生物之间的相互作用机制,有助于深入理解时滞效应如何影响海洋生态系统的动态运行趋势,为进一步研究海洋浮游生物之间相互制约、相互协调的生长动态机制提供一定的理论支撑。

关键词 :溶解氧,浮游生物,时滞,Hopf分支,稳定性

Copyright © 2017 by author and Hans Publishers Inc.

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

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

1. 引言

海洋浮游生物是指悬浮在水层中常随水流移动的海洋生物,包括海洋浮游植物和海洋浮游动物,其特点是这类生物缺乏发达的运动器官,没有或仅有微弱的游动能力和多数分布于水体的上层或表层。海洋浮游植物是指在海水中营浮游生活的微小植物,主要包括单细胞藻类和细菌,其营养方式为自养方式,其中,藻类具有叶绿素或其它色素,能进行光合作用制造有机物,是水域生态系统中的生产者,细菌是生态系统的还原者(也可以是生产者),它们一般分布于海洋的真光层。海洋浮游动物一类经常在海水中浮游,本身不能制造有机物的异养型无脊椎动物和脊索动物幼体的总称,它们必须依赖已有的有机物为营养来源,多为滤食性。同时,浮游动物是海洋生态系统中的消费者,可分布于真光层,也可分布于较深的水层。

众所周知,海洋浮游植物是海洋生态系统中重要的生产者,是食物链的基础环节,在海洋生态系统的物质循环和能量转化过程中起着重要作用。为海洋中的生命活动提供了能量,有助于勘探海底石油,有助于研究海洋古地质和古环境,有助于研究、防治海洋环境污染。其生长繁殖除主要受到自身生物学特性影响外,还受到周围环境因素如海流、扰动、温度、盐度、营养盐和其他生物因素的影响。海洋浮游动物是海洋生态系统中非常重要的一大生态类群,种类组成繁殖、数量大、分布广,有着极其重要的生态学意义。浮游动物通过捕食作用控制浮游植物的数量,同时作为鱼类等高层营养者的饵料,其数量变化可以直接影响鱼类等的资源量,在海洋生态系统的结构和功能中起着重要的调控作用。同时,海洋浮游动物既是初级生产的消费者,也是次级生产者,在浮游植物和游泳动物食物链之间扮演重要角色,在海洋生态系统动态变化中起着重要的调控作用,在海洋环境保护研究中有重要意义。因此研究海洋浮游生物生长动态机制与相互制约、相互协调生长动态机制具有至关重要的现实意义。就目前来说,应用数学动态模型与数值模拟方法来研究此类问题是具有一定的优越性。因而数学模型已成为一些学者研究浮游生物动态生长及其变化趋势的常用工具之一 [1] - [7] 。例如,论文 [8] 建立了一类N-P-Z数学动态模型,讨论了系统的边界性和Hopf分支存在性等问题,给出了系统具有此类特性的阈值条件及其生物学意义。

溶解在海水中的氧(溶解氧)是海洋生命活动不可缺少的物质,它的含量在海洋中的分布,既受化学过程和生物过程的影响,还受物理过程的影响。最为重要的是海水中溶解氧是海洋水生动物、植物、微生物等进行新陈代谢的必要物质。随着研究的深入,一些学者发现溶解氧对藻类的爆发有着举足轻重的影响 [9] - [17] 。论文 [14] 建立了一类溶解氧-浮游生物动力学模型,研究了系统平衡点的存在性与稳定性,探讨了溶解氧如何影响浮游生物种群数量的动态变化趋势。论文 [15] 构建了一类数学动态模型去预测水体中溶解氧耗散动态变化情况,获得一些比较有意义的结果。论文 [18] 提出了一类数学动态模型来探析溶解氧如何影响浮游生物动力学性质,研究了模型平衡点的存在性与稳定性,探讨了系统分支的存在性及其相关动力学性质,其动态模型如下:

d P d t = r P ( ξ + D 0 D ) δ 1 P β P Z ( P + a ) , d D d t = η ( D 0 D ) α 1 P α 2 Z , d Z d t = β 1 P Z ( P + a ) δ 2 Z . (1)

考虑到生物学意义,所有的参数都为正, P , Z 分别代表浮游植物和浮游动物的种群数量, D 表示在任意时刻 t 的溶解氧的含量。此外 D 0 表示溶解氧在水体中的吸收量。浮游植物的内禀增长率随着 D 的增加而增加,且假设在 D = D 0 时有最大的增长率 r / ξ r / ( ξ + D 0 D ) 视为浮游植物吸收溶解氧的增长方程。 δ 1 δ 2 分别表示浮游植物和浮游动物的死亡率,其中包括水体的自然冲刷。 β 表示浮游植物被浮游动物觅食的比率, β 1 是浮游植物在浮游动物体内的转换率。 η 表示水体中溶解氧的补给率。 α 1 , α 2 分别表示每单位被浮游植物和浮游动物消耗掉的溶解氧。

然而在现实海洋生态系统中,浮游生物之间的相互作用并不是瞬时的,而是存在一些时滞效应,这些时滞可能是由于浮游动物的妊娠期或者浮游生物对食物的消化等其它因素引起的,因此将时滞引入到生态系统中可以更加合理地描述与揭示现实生态系统动态运行规律。大量的研究表明时滞系统比常规系统有更加丰富的动力学性质 [19] - [25] 。论文 [26] 研究了一个包含由妊娠期引起的时滞生态动力学模型,分别在六种情况下讨论了时滞参数如何影响系统的动力学性质和解析出相关的临界条件。论文 [27] 构建了一类具有双时滞效应的浮游生态动力学模型,其时滞效应是由浮游植物产生毒素和浮游动物消化食物所引起的。同时,作者将随机因素引入该系统中,探讨了随机稳定性和随机Hopf分支的存在性,并利用李雅普诺夫函数证明了系统含有双时滞情况下内平衡点的渐近稳定性。

基于论文 [18] 中的相关结果,本论文将浮游动物消化浮游植物以及浮游生物吸收溶解氧所引起的时滞效应考虑在内,建立了一类具有双时滞效应的海洋生态动力学新模型,其可以表示如下:

d D d t = η ( D 0 D ) α 1 P α 2 Z , d P d t = r P ( ξ + D 0 D ( t τ 1 ) ) δ 1 P β P Z ( P + a ) , d Z d t = β 1 P ( t τ 2 ) Z ( P ( t τ 2 ) + a ) δ 2 Z . (2)

其中 τ 1 是浮游植物吸收溶解氧引起的时滞项, τ 2 是由于浮游动物消化食物所引起的时滞项。本论文的主要工作是解析时滞参数如何影响该生态系统的动力学行为及其给出相关临界条件,为后续研究提供理论基础。

2. 局部稳定性和Hopf分支

现在主要讨论系统(2)内平衡点的稳定性和Hopf分支的存在性。从论文 [18] 中得知:当 β 1 > δ 2 时系统(2)存在一个内平衡点,记为 E * = ( D * , P * , Z * ) ,其中:

D * = D 0 a α 1 δ 2 η ( β 1 δ 2 ) α 2 Z * η P * = a δ 2 β 1 δ 2 Z * = a β 1 β ( β 1 δ 2 ) ( r ξ + D 0 D * δ 1 ) . (3)

首先对系统(2)进行换元变形: u 1 ( t ) = D ( t ) D * u 2 ( t ) = P ( t ) P * u 3 ( t ) = Z ( t ) Z * 。将内平衡点 E * 转化到原点,系统(2)可以重新写为

u 1 = a 11 u 1 + a 12 u 2 + a 13 u 3 + F 1 , u 2 = a 21 u 1 ( t τ 1 ) + a 22 u 2 + a 23 u 3 + F 2 , u 3 = a 32 u 2 ( t τ 2 ) + F 3 , (4)

其中:

a 11 = η < 0 ; a 12 = α 1 < 0 ; a 13 = α 2 < 0 ; a 21 = r P * ( ξ + D 0 D * ) 2 > 0 ; a 22 = r ξ + D 0 D * δ 1 β a Z * ( a + P * ) 2 > 0 ; a 23 = β P * a + P * < 0 ; a 32 = β 1 Z * a ( a + P * ) 2 .

f 1 = 0 , f 2 = ( a β ( a + P * ) 2 ) y ( t ) z ( t ) + a β Z * ( a + P * ) 3 y 2 ( t ) a β Z * ( a + P * ) 4 y 3 ( t ) + , f 3 = a β 1 ( a + P * ) 2 y ( t τ 2 ) z ( t ) a β 1 Z * ( a + P * ) 3 y 2 ( t τ 2 ) + a β 1 Z * ( a + P * ) 4 y 3 ( t τ 2 ) + .

则系统(4)对应的线性部分为:

u 1 = a 11 u 1 + a 12 u 2 + a 13 , u 2 = a 21 u 1 ( t τ 1 ) + a 22 u 2 + a 23 u 3 , u 3 = a 32 u 2 ( t τ 2 ) , (5)

系统(5)对应的特征方程如下:

λ 3 + A λ 2 + B λ + C λ e λ τ 1 + ( D λ + E ) e λ τ 2 + F e λ ( τ 1 + τ 2 ) = 0. (6)

其中:

A = ( a 11 + a 22 ) , B = a 11 a 22 , C = a 12 a 21 , D = a 32 a 23 , E = a 11 a 23 a 32 , F = a 13 a 21 a 32 .

第一种情况: τ 1 = τ 2 = 0

定理2.1: τ 1 = τ 2 = 0 时,满足条件 ( H 1 ) 时,内平衡点 E * 是局部稳定的。

证明:当 τ 1 = τ 2 = 0 时,方程(6)变为

λ 3 + A 11 λ 2 + A 12 λ + A 13 = 0 , (7)

其中:

A 11 = A , A 12 = B + C + D , A 13 = E + F .

根据Routh-Hurwitz判据,满足条件 ( H 1 ) A 11 > 0 A 11 A 12 > A 13 时, E * 是局部渐近稳定的。

第二种情况: τ 1 > 0 , τ 2 = 0

定理2.2: τ 1 > 0 , τ 2 = 0 时,满足条件 ( H 21 , H 22 ) ,系统内平衡点 E * τ 1 [ 0 , τ 10 ) 是局部渐近稳定的,并且当 τ 1 = τ 10 时出现Hopf分支。

证明:当 τ 1 > 0 , τ 2 = 0 时,方程(6)化为

λ 3 + A λ 2 + ( B + D ) λ + ( C λ + F ) e λ τ 1 + E = 0. (8)

假设 λ = i ω 1 ( ω 1 > 0 ) 是方程(8)的纯虚根,将其代入方程(8),得到如下方程组

ω 1 3 + ( B + D ) ω 1 F sin ( ω 1 τ 1 ) + C ω 1 cos ( ω 1 τ 1 ) = 0 , A ω 1 2 + E + C ω 1 sin ( ω 1 τ 1 ) + F cos ( ω 1 τ 1 ) = 0. (9)

将两个方程平方并且相加得到方程(10):

ω 1 6 + [ A 2 2 ( B + D ) ] ω 1 4 + [ ( B + D ) 2 2 A E C 2 ] ω 1 2 + E 2 F 2 = 0 , (10)

假设 v 1 = ω 1 2 ,然后方程(10)变为

v 1 3 + [ A 2 2 ( B + D ) ] v 1 2 + [ ( B + D ) 2 2 A E C 2 ] v 1 + E 2 F 2 = 0 , (11)

接下来定义一个方程:

f 1 ( v 1 ) = v 1 3 + [ A 2 2 ( B + D ) ] v 1 2 + [ ( B + D ) 2 2 A E C 2 ] v 1 + E 2 F 2 (12)

假设 ( H 21 ) :文献 [28] 中定理(2.1)情况(i),(iii)成立,则 f 1 ( v 1 ) 至少有一个正根,不失一般性,假设方程(11)的根为 v 11 , v 12 , v 13 ,则有 ω 1 k = v 1 k , ( k = 1 , 2 , 3 ) 根据方程组(9)得到:

τ 1 k j = 1 ω 1 k ( arccos C ω 1 k 4 + [ A F ( B + D ) C ] ω 1 k 2 E F C 2 ω 1 k 2 + F 2 + 2 j π ) , ( j = 0 , 1 , 2 , k = 1 , 2, 3 )

假设 τ 10 = min k = 1 , 2 , 3 τ 1 k ( 0 ) ,由Hopf分支定理 [29] 得知,还需要证明横街性条件,对方程(8)两边关于 τ 1 微分求导得

( d λ d τ 1 ) 1 = 3 λ 2 + 2 A λ + B + D + C e λ τ 1 λ ( C λ + F ) e λ τ 1 τ 1 λ , (13)

λ = i ω 10 代入(13),通过简单计算得到

Re ( d λ d τ 1 ) λ = i ω 10 1 = ω 10 2 Δ f 1 ( ω 10 2 ) , (14)

其中:

Δ = ( F ω 10 sin ( ω 10 τ 1 ) C ω 10 2 cos ( ω 10 τ 1 ) ) 2 + ( F ω 10 cos ( ω 10 τ 1 ) + C ω 10 2 sin ( ω 10 τ 1 ) ) 2 ,

得知当 ( H 22 ) : f 1 ( ω 10 2 ) 0 成立时,则 Re ( d λ d τ 1 ) λ = i ω 10 1 0 。证毕。

第三种情况: τ 2 > 0 , τ 1 = 0

定理2.3: τ 2 > 0 , τ 1 = 0 时,满足 ( H 21 , H 32 ) ,系统内平衡点 E * τ 2 [ 0 , τ 20 ) 是局部渐近稳定的,并且当 τ 2 = τ 20 时出现Hopf分支。

证明:当 τ 2 > 0 , τ 1 = 0 时,方程(6)化为

λ 3 + A λ 2 + ( B + C ) λ + D λ e λ τ 2 + ( E + F ) e λ τ 2 = 0. (15)

与第二种情况相似,假设 λ = i ω 2 ( ω 2 > 0 ) 是方程(15)的根,将其代入得到:

ω 2 3 + ( B + C ) ω 2 ( E + F ) sin ( ω 2 τ 2 ) + D ω 2 cos ( ω 2 τ 2 ) = 0 , A ω 2 2 + ( E + F ) ω 2 cos ( ω 2 τ 2 ) + D sin ( ω 2 τ 2 ) = 0. (16)

对方程组平方并且相加得到

ω 2 6 + ( A 2 2 ( B + C ) ) ω 2 4 + [ ( B + C ) 2 D 2 ] ω 2 2 ( E + F ) 2 = 0 , (17)

假设 v 2 = ω 2 2 ,则方程(17)化为

v 2 3 + ( A 2 2 ( B + C ) ) v 2 2 + [ ( B + D ) 2 D 2 ] v 2 ( E + F ) 2 = 0 , (18)

然后定义

f 2 ( v 2 ) = v 2 3 + ( A 2 2 ( B + C ) ) v 2 2 + [ ( B + D ) 2 D 2 ] v 2 ( E + F ) 2 (19)

假设 ( H 21 ) 成立,则方程(18)至少有一个正根,不失一般性,将方程(18)的正根记为 v 21 , v 22 , v 23 ,然后有 ω 2 k = v 2 k , ( k = 1 , 2 , 3 ) 。根据方程(16)得到:

τ 2 k j = 1 ω 2 k ( arccos D ω 2 k 4 + [ A ( E + F ) ( B + C ) D ] ω 2 k 2 D 2 ω 2 k 2 + ( E + F ) 2 + 2 j τ ) , ( j = 0 , 1 , 2 , k = 1 , 2, 3 )

假设 τ 20 = min k = 1 , 2 , 3 τ 2 k ( 0 ) 与第二种情况相似,还需要证明横街性条件,对方程(15)两边关于 τ 2 微分求导得

( d λ d τ 2 ) 1 = 3 λ 2 + 2 A λ + B + C + D e λ τ 2 λ ( D λ + E + F ) e λ τ 2 τ 2 λ , (20)

λ = i ω 20 代入(20),通过简单计算得到

Re ( d λ d τ 2 ) λ = i ω 20 1 = ω 20 2 Δ f 2 ( ω 20 2 ) (21)

其中:

Δ = ( D ω 20 2 cos ( ω 20 τ 2 ) + ( E + F ) ω 20 sin ( ω 20 τ 2 ) ) 2 + ( ( E + F ) ω 20 cos ( ω 20 τ 2 ) + D ω 20 2 ( sin ω 20 τ 2 ) ) 2 ,

得知当 ( H 32 ) : f 2 ( ω 20 2 ) 0 ,则 Re ( d λ d τ 2 ) λ = i ω 20 1 0 。证毕。

第四种情况: τ 1 = τ 2 = τ

定理2.4: τ 1 = τ 2 = τ 时,满足 ( H 41 , H 42 ) ,系统内平衡点 E * τ [ 0 , τ 30 ) 是局部渐近稳定的,并且当 τ = τ 30 时出现Hopf分支。

证明:当 τ 1 = τ 2 = τ 时,方程(6)化为

λ 3 + A λ 2 + B λ + [ ( C + D ) λ + E ] e λ τ + F e 2 λ τ = 0. (22)

方程两边同时乘以 e λ τ ,得到

( λ 3 + A λ 2 + B λ ) e λ τ + ( C + D ) λ + E + F e λ τ = 0. (23)

假设 λ = i ω 3 ( ω 3 > 0 ) 是方程(23)的根,将其代入得到

sin ( ω 3 τ ) ( B ω 3 ω 3 3 ) + cos ( ω 3 τ ) ( A ω 3 2 F ) = E , sin ( ω 3 τ ) ( A ω 3 2 + F ) + cos ( ω 3 τ ) ( ω 3 3 B ω 3 ) = ( C + D ) ω 3 . (24)

根据方程(24),得到如下方程组:

sin ( ω 3 τ ) = E ( ω 3 3 B ω 3 ) ( C + D ) ( A ω 3 2 F ) ω 3 ( ω 3 3 B ω 3 ) 2 + ( F A 2 ω 3 4 ) , cos ( ω 3 τ ) = ( C + D ) ( B ω 3 ω 3 3 ) ω 3 E ( A ω 3 2 + F ) ( ω 3 3 B ω 3 ) 2 + ( F A 2 ω 3 4 ) . (25)

利用 sin 2 ( ω 3 τ ) + cos 2 ( ω 3 τ ) = 1 ,就有

ω 3 12 + E 41 ω 3 10 + E 42 ω 3 8 + E 43 ω 3 6 + E 44 ω 3 4 + E 45 ω 3 2 + E 46 = 0 , (26)

其中:

E 41 = 2 A 2 4 B ; E 42 = 6 B 2 + A 4 4 A 2 B ( C + D ) 2 ; E 43 = 4 B 3 2 F 2 + 2 A 2 B 2 ( C + D ) 2 ( A 2 2 B ) E 2 ; E 44 = B 4 2 A F 2 + 4 B F 2 4 E F ( C + D ) ( C + D ) 2 ( B 2 2 A F ) + E 2 ( A 2 + 2 B ) ; E 45 = 2 B 2 F 2 + 4 B E F ( C + D ) ( C + D ) 2 F 2 E 2 ( B 2 + 2 A F ) ; E 46 = F 4 + E 2 F 2 ;

假设 v 3 = ω 3 2 ,则方程(26)变为

v 3 6 + E 41 v 3 5 + E 42 v 3 4 + E 43 v 3 3 + E 44 v 3 2 + E 45 v 3 + E 46 = 0 , (27)

接下来定义

f 3 ( v 3 ) = v 3 6 + E 41 v 3 5 + E 42 v 3 4 + E 43 v 3 3 + E 44 v 3 2 + E 45 v 3 + E 46 , (28)

假设 ( H 41 ) :方程(28)至少有一个正根成立。不失一般性,记方程(27)的正根为 v 31 , v 32 , , v 36 ,则有 ω 3 k = v 3 k , ( k = 1 , 2 , , 6 ) ,由方程(24)可以得到

τ 3 k j = 1 ω 3 k ( arccos ( C + D ) ( B ω 3 ω 3 3 ) ω 3 E ( A ω 3 2 + F ) ( ω 3 3 B ω 3 ) 2 + ( F A 2 ω 3 4 ) + 2 j π ) , ( j = 0 , 1 , 2 , k = 1 , 2, , 6 )

假设 τ 30 = min k = 1 , 2 , , 6 τ 3 k ( 0 ) ,然后对方程(23)两边关于 τ 微分求导得

( d λ d τ ) 1 = ( 3 λ 2 + 2 A λ + B ) e λ τ + C + D F λ e λ τ ( λ 3 + A λ 2 + B λ ) λ e λ τ τ λ , (29)

λ = i ω 30 代入(29),通过简单计算得到

Re ( d λ d τ ) λ = i ω 30 1 = F 41 F 43 + F 42 F 44 F 43 2 + F 44 2 ,

其中:

F 41 = 3 ω 30 2 cos ( ω 30 τ ) 2 A ω 30 sin ( ω 30 τ ) + B cos ( ω 30 τ ) + C + D ; F 42 = 3 ω 30 2 sin ( ω 30 τ ) + 2 A ω 30 cos ( ω 30 τ ) + B sin ( ω 30 τ ) ; F 43 = F ω 30 sin ( ω 30 τ ) ω 30 4 cos ( ω 30 τ ) + B ω 30 2 cos ( ω 30 τ ) A ω 30 3 sin ( ω 30 τ ) ; F 44 = F ω 30 cos ( ω 30 τ ) ω 30 4 sin ( ω 30 τ ) + B ω 30 2 sin ( ω 30 τ ) + A ω 30 3 cos ( ω 30 τ ) .

得知当 ( H 42 ) : F 41 F 43 + F 42 F 44 0 ,则 Re ( d λ d τ ) λ = i ω 30 1 0 。证毕。

第五种情况: τ 1 ( 0 , τ 10 ) , τ 2 > 0

定理2.5: τ 1 ( 0 , τ 10 ) , τ 2 > 0 时,满足 ( H 51 , H 52 ) ,系统内平衡点 E * τ 2 [ 0 , τ 20 * ) 是局部渐近稳定的,并且当 τ = τ 20 * 时出现Hopf分支。

证明:当 τ 1 ( 0 , τ 10 ) , τ 2 > 0 时,将 τ 2 作为参数,假设 λ = i ω 2 * ( ω 2 * > 0 ) 是方程(6)的根,将其代入得到

R 51 sin ( ω 2 * τ 2 ) + R 52 cos ( ω 2 * τ 2 ) = R 53 , R 51 cos ( ω 2 * τ 2 ) R 52 sin ( ω 2 * τ 2 ) = R 54 . (30)

其中:

R 51 = D ω 2 * F sin ( ω 2 * τ 1 ) ; R 52 = E + F cos ( ω 2 * τ 1 ) ; R 53 = A ω 2 * 2 C ω 2 * sin ( ω 2 * τ 1 ) ; R 54 = ω 2 * 3 B ω 2 * C ω 2 * cos ( ω 2 * τ 1 ) ;

对方程组(30)平方相加得

F 1 ( ω 2 * ) + F 2 ( ω 2 * ) sin ( ω 2 * τ 1 ) + F 3 ( ω 2 * ) cos ( ω 2 * τ 1 ) = 0 , (31)

其中:

F 1 ( ω 2 * ) = ω 2 * 6 + ( A 2 2 B ) ω 2 * 4 + ( C 2 + B 2 D 2 ) ω 2 * 2 ( E 2 + F 2 ) ; F 2 ( ω 2 * ) = 2 D F ω 2 * 2 A C ω 2 * 3 ; F 3 ( ω 2 * ) = 2 B C ω 2 * 2 2 C ω 2 * 4 2 E F ;

假设 ( H 51 ) :(31)存在有限个正根成立,且记为 ω 2 k * , ( k = 1 , 2 , , l 1 ) 。由方程(30)有

τ 2 k * ( j ) = 1 ω 2 k * ( arccos R 52 R 53 + R 51 R 54 R 51 2 + R 52 2 + 2 j π ) , ( j = 0 , 1 , 2 , k = 1 , 2, , l 1 )

假设 τ 20 * = min τ 2 k * ( 0 ) ( k = 1 , 2 , , l 1 ) ,然后对方程(6)两边关于 τ 2 微分求导得

( d λ d τ 2 ) 1 = f 20 + f 21 e λ τ 1 + f 22 e λ τ 2 + f 23 e λ ( τ 1 + τ 2 ) f 24 e λ τ 2 + f 25 e λ ( τ 1 + τ 2 ) τ 2 λ , (32)

其中:

f 20 = 3 λ 2 + 2 A λ + B ; f 21 = C C λ τ 1 ; f 22 = D ; f 23 = τ 1 F ; f 24 = ( D λ + E ) λ ; f 25 = F λ ;

λ = i ω 20 * , τ 2 = τ 20 * 代入(32)得

Re ( d λ d τ 2 ) λ = i ω 20 * 1 = F 51 F 53 + F 52 F 54 F 53 2 + F 54 2

其中:

F 51 = 3 ω 20 * 2 + B + C cos ( ω 20 * τ 1 ) C ω 20 * τ 1 sin ( ω 20 * τ 1 ) + D cos ( ω 20 * τ 20 * ) τ 1 F cos ( ω 20 * τ 1 ) ) cos ( ω 20 * τ 20 * ) + τ 1 F sin ( ω 20 * τ 1 ) sin ( ω 20 * τ 20 * ) ; F 52 = 2 A ω 20 * C sin ( ω 20 * τ 1 ) C ω 20 * τ 1 cos ( ω 20 * τ 1 ) D sin ( ω 20 * τ 20 * ) + τ 1 F sin ( ω 20 * τ 1 ) cos ( ω 20 * τ 20 * ) + τ 1 F cos ( ω 20 * τ 1 ) sin ( ω 20 * τ 20 * ) ; F 53 = D ω 20 * 2 cos ( ω 20 * τ 20 * ) + E ω 20 * sin ( ω 20 * τ 20 * ) + F ω 20 * sin ( ω 20 * τ 1 ) cos ( ω 20 * τ 20 * ) + F ω 20 * cos ( ω 20 * τ 1 ) sin ( ω 20 * τ 20 * ) ; F 54 = D ω 20 * 2 sin ( ω 20 * τ 20 * ) + E ω 20 * cos ( ω 20 * τ 20 * ) + F ω 20 * cos ( ω 20 * τ 1 ) cos ( ω 20 * τ 20 * ) F ω 20 * sin ( ω 20 * τ 1 ) sin ( ω 20 * τ 20 * ) .

假设 ( H 52 ) : F 51 F 53 + F 52 F 54 0 成立,则有 Re ( d λ d τ 2 ) λ = i ω 20 * 1 0 。证毕。

第六种情况: τ 2 ( 0 , τ 20 ) , τ 1 > 0

定理2.6: τ 2 ( 0 , τ 20 ) , τ 1 > 0 时,满足 ( H 61 , H 62 ) ,系统内平衡点 E * τ 1 [ 0 , τ 10 * ) 是局部渐近稳定的,并且当 τ = τ 10 * 时出现Hopf分支。

证明:当 τ 2 ( 0 , τ 20 ) , τ 1 > 0 时,将 τ 1 作为参数,假设 λ = i ω 1 * ( ω 1 * > 0 ) 是方程(6)的根,代入可得

R 61 sin ( ω 1 * τ 1 ) + R 62 cos ( ω 1 * τ 1 ) = R 63 , R 61 cos ( ω 1 * τ 1 ) R 62 sin ( ω 1 * τ 1 ) = R 64 . (33)

其中:

R 61 = C ω 1 * F sin ( ω 1 * τ 2 ) ; R 62 = F cos ( ω 1 * τ 2 ) ; R 63 = A ω 1 * 2 E cos ( ω 1 * τ 2 ) D ω 1 * sin ( ω 1 * τ 2 ) ; R 64 = ω 1 * 3 B ω 1 * + E sin ( ω 1 * τ 2 ) D ω 1 * cos ( ω 1 * τ 2 ) ;

对方称组(33)平方相加得

G 1 ( ω 1 * ) + G 2 ( ω 1 * ) sin ( ω 1 * τ 2 ) + G 3 ( ω 1 * ) cos ( ω 1 * τ 2 ) = 0. (34)

其中:

G 1 ( ω 1 * ) = ω 1 * 6 + ( A 2 2 B ) ω 1 * 4 + ( B 2 + D 2 C 2 ) ω 1 * 2 + E 2 F 2 ; G 2 ( ω 1 * ) = ( 2 E 2 A D ) ω 1 * 3 + ( 2 C F 2 B E ) ω 1 * ; G 3 ( ω 1 * ) = 2 D ω 1 * 4 + ( 2 B D 2 A E ) ω 1 * 2 ;

假设 ( H 61 ) :(33)存在有限个正根成立,且记为 ω 1 k * , ( k = 1 , 2 , , l 2 ) 。由方程(33)得到

τ 1 k * ( j ) = 1 ω 1 k * arccos ( R 62 R 63 + R 61 R 64 R 61 2 + R 62 2 + 2 j τ ) , ( j = 0 , 1 , 2 , k = 1 , 2, , l 2 )

假设 τ 10 * = min τ 1 k * ( 0 ) ( k = 1 , 2 , , l 2 ) ,对方称(6)关于 τ 1 微分求导,并将 λ = i ω 10 * τ 1 = τ 10 * 代入得

Re ( d λ d τ 1 ) λ = i ω 10 * 1 = F 61 F 63 + F 62 F 64 F 63 2 + F 64 2 (35)

其中:

F 61 = 3 ω 10 * 2 + B + C cos ( ω 10 * τ 10 * ) + D cos ( ω 10 * τ 2 ) τ 2 E cos ( ω 10 * τ 2 ) τ 2 D ω 10 * sin ( ω 10 * τ 2 ) τ 2 F cos ( ω 10 * τ 10 * ) cos ( ω 10 * τ 2 ) + τ 2 F sin ( ω 10 * τ 10 * ) sin ( ω 10 * τ 2 ) ; F 62 = 2 A ω 10 * C sin ( ω 10 * τ 10 * ) D sin ( ω 10 * τ 2 ) τ 2 D ω 10 * cos ( ω 10 * τ 2 ) + τ 2 E sin ( ω 10 * τ 2 ) + τ 2 F sin ( ω 10 * τ 10 * ) cos ( ω 10 * τ 2 ) + τ 2 F sin ( ω 10 * τ 2 ) cos ( ω 10 * τ 10 * ) ; F 63 = F ω 10 * sin ( ω 10 * τ 10 * ) cos ( ω 10 * τ 2 ) + F ω 10 * cos ( ω 10 * τ 10 * ) sin ( ω 10 * τ 2 ) C ω 10 * 2 cos ( ω 10 * τ 10 * ) ; F 64 = F ω 10 * cos ( ω 10 * τ 10 * ) cos ( ω 10 * τ 2 ) F ω 10 * sin ( ω 10 * τ 10 * ) sin ( ω 10 * τ 2 ) + C ω 10 * 2 sin ( ω 10 * τ 10 * ) ;

假设 ( H 62 ) : F 61 F 63 + F 62 F 64 0 成立,则有 Re ( d λ d τ 1 ) λ = i ω 10 * 1 0 。证毕。

3. 分支性质

本部分利用范数定理和中心流行定理 [29] ,讨论关于 τ 1 的Hopf分支性质,其中 τ 2 [ 0 , τ 20 ) 。假设在 τ 2 [ 0 , τ 20 ) 条件下有 τ 2 < τ 10 ,由前面的章节得知当 τ = τ 10 时系统(2)在正平衡点 E * = ( D * , P * , Z * ) 处存在一个Hopf分支。

假设 τ 1 = μ + τ 10 ,则 μ = 0 是系统(2)的Hopf分支值。令 u 1 ( t ) = D ( t ) D * u 2 ( t ) = P ( t ) P * u 3 ( t ) = Z ( t ) Z * ,并重新设置时滞为 t ( t / τ 1 ) ,将系统(2)重新写为

u ˙ ( t ) = L μ u t + F ( μ , u t ) , (36)

其中:

L μ ϕ = ( τ 10 + μ ) ( A ϕ ( 0 ) + B ϕ ( 1 ) + C ϕ ( τ 2 τ 10 ) ) ,

F ( μ , ϕ ) = ( τ 10 + μ ) ( F 1 , F 2 , F 3 ) T .

ϕ ( θ ) = ( ϕ 1 ( θ ) , ϕ 2 ( θ ) , ϕ 3 ( θ ) ) T C ( [ 1 , 0 ] , R 3 ) ,

A = ( a 11 a 12 a 13 0 a 22 a 23 0 0 0 ) , B = ( 0 0 0 a 21 0 0 0 0 0 ) , C = ( 0 0 0 0 0 0 0 a 32 0 ) ,

F 1 = 0 , F 2 = ( a β ( a + P * ) 2 ) ϕ 2 ( 0 ) ϕ 3 ( 0 ) + a β Z * ( a + P * ) 3 ϕ 2 2 ( 0 ) a β Z * ( a + P * ) 4 ϕ 2 3 ( 0 ) + , F 3 = a β 1 ( a + P * ) 2 ϕ 2 ( τ 2 τ 10 ) ϕ 3 ( 0 ) a β 1 Z * ( a + P * ) 3 ϕ 2 2 ( τ 2 τ 10 ) + a β 1 Z * ( a + P * ) 4 ϕ 2 3 ( τ 2 τ 10 ) + .

根据Riesz表示定理,对一个有界变量存在一个 3 × 3 的矩阵方程 η ( θ , μ ) : [ 1 , 0 ] R 3 2 ,其中 θ [ 1 , 0 ] ,则有

L μ ϕ = 1 0 d η ( θ , μ ) ϕ ( θ ) , ϕ C ( [ 1 , 0 ] , R 3 ) .

选择

ν ( θ ) = { ( τ 10 * + μ ) ( A + B + C ) , θ = 0 , ( τ 10 * + μ ) ( B + C ) , θ [ τ 2 τ 10 * , 0 ) , ( τ 10 * + μ ) B , θ ( 1 , τ 2 τ 10 * ) , 0 , θ = 1.

ϕ C ( [ 1 , 0 ] , R 3 ) ,定义

Γ ( μ ) ϕ = { d ϕ ( θ ) d θ , θ [ 1 , 0 ) , 1 0 d η ( θ , μ ) ϕ ( θ ) , θ = 0 , R ( μ ) ϕ = { 0 , θ [ 1 , 0 ) , F ( μ , ϕ ) , θ = 0.

因此,系统(36)可以重新写为

u ˙ t = Γ ( μ ) u t + R ( μ ) u t , (37)

其中 u t = u ( t + θ ) = ( u 1 ( t + θ ) , u 2 ( t + θ ) , u 3 ( t + θ ) ) θ [ 1 , 0 ]

φ C 1 ( [ 0 , 1 ] , ( R 3 ) * ) ,定义 Γ 的耦合算子 Γ * 如下:

Γ * ( φ ) = { d φ ( s ) d s , s ( 0 , 1 ] , 1 0 d η T ( s , 0 ) φ ( s ) , s = 0 ,

一个线性内积定义如下:

φ ( s ) , ϕ ( θ ) = φ ¯ ( 0 ) ϕ ( 0 ) θ = 1 0 ξ = 0 θ φ ¯ ( ξ θ ) d η ( θ ) ϕ ( ξ ) d ξ (38)

其中: η ( θ ) = η ( θ , 0 )

由前面的研究得知 ± i ω 1 * τ 10 * Γ ( 0 ) 的特征值,则 ± i ω 1 * τ 10 * 同样是 Γ * ( 0 ) 的特征值。

假设 q ( θ ) = ( 1 , q 2 , q 3 ) T e i ω 1 * τ 10 * θ Γ ( 0 ) 关于特征值 i ω 1 * τ 10 * 的特征向量,同样假设 q * ( s ) = D ( 1 , q 2 * , q 3 * ) e i ω 1 * τ 10 * s Γ * ( 0 ) 关于特征值 i ω 1 * τ 10 * 的特征向量。通过计算可以得到:

q 2 = ω 1 * 2 a 11 i ω 1 * a 12 i ω 1 * + a 13 a 32 e i ω 1 * τ 2 , q 3 = ( i ω 1 * a 11 ) ( a 32 e i ω 1 * τ 2 ) a 13 a 32 e i ω 1 * τ 2 a 12 i ω 1 * ; q 2 * = ( i ω 1 * a 11 ) a 21 e i ω 1 * τ 10 * , q 3 * = ( i ω 1 * + a 11 ) a 23 a 13 a 21 e i ω 1 * τ 10 * i ω 1 * a 21 e i ω 1 * τ 10 * .

从方程(38),得到

D = [ 1 + q 2 * q ¯ 2 + q 3 * q ¯ 3 + τ 2 * e i ω 1 * τ 10 * ( q 2 * τ 10 * a 21 + q ¯ 2 q 3 * τ 2 a 32 ) ] 1

满足 q * ( s ) , q ( θ ) = 1 q * ( s ) , q ¯ ( θ ) = 0

引用 [30] 给出的算法可以得到用来决定分支方向和分支周期解稳定性的参数如下:

g 20 = 2 D ¯ τ 10 * ( a q 2 q 3 ( a + P * ) 2 ( β 1 e i ω 1 * τ 2 * q 3 * ¯ β q 2 * ¯ ) + a Z * ( a + P * ) 3 ( β q 2 2 q 2 * ¯ β 1 q 2 2 e 2 i ω 1 * τ 2 * q 3 * ¯ ) ) ,

g 11 = D ¯ τ 10 * ( a q 3 q 2 ¯ ( a + P * ) 2 ( β 1 e i ω 1 * τ 2 * q 3 * ¯ β q 2 * ¯ ) + a q 3 ¯ q 2 ( a + P * ) 2 ( β 1 e i ω 1 * τ 2 * q 3 * ¯ β q 2 * ¯ ) + 2 a Z * q 2 q 2 ¯ ( a + P * ) 3 ( β q 2 * ¯ β 1 q 3 * ¯ ) ) ,

g 02 = 2 D ¯ τ 10 * ( a q 3 ¯ q 2 ¯ ( a + P * ) 2 ( β 1 e i ω 1 * τ 2 * q 3 * ¯ β q 2 * ¯ ) + a Z * q 2 ¯ 2 ( a + P * ) 3 ( β q 2 * ¯ β 1 e 2 i ω 1 * τ 2 * q 3 * ¯ ) ) ,

g 21 = 2 D ¯ τ 10 * { q 2 * ¯ [ ( a β ( a + P * ) 2 ) ( W 20 ( 2 ) ( 0 ) 1 2 q 3 ¯ + W 11 ( 2 ) ( 0 ) q 3 + q 2 ¯ 1 2 W 20 ( 3 ) ( 0 ) + W 11 ( 3 ) ( 0 ) q 2 ) ] + q 2 * ¯ a Z * β ( a + P * ) 3 ( 2 W 11 ( 2 ) ( 0 ) q 2 + W 20 ( 2 ) ( 0 ) q 2 ¯ ) + q 3 * ¯ a β 1 ( a + P * ) 2 ( W 11 ( 3 ) ( 0 ) e i ω 1 * τ 2 * q 2 + 1 2 e i ω 1 * τ 2 * q 2 ¯ W 20 ( 3 ) ( 0 ) + 1 2 W 20 ( 2 ) ( τ 2 * τ 10 * ) q 3 ¯ + W 11 ( 2 ) ( τ 2 * τ 10 * ) q 3 ) q 3 * ¯ a β 1 Z * ( a + P * ) 3 ( 2 W 11 ( 2 ) ( τ 2 * τ 10 * ) e i ω 1 * τ 2 * q 2 + W 20 ( 2 ) ( τ 2 * τ 10 * ) e i ω 1 * τ 2 * q 2 ¯ ) } .

并且可以得到:

W 20 ( θ ) = i g 20 q ( 0 ) ω 1 * τ 10 * e i ω 1 * τ 10 * θ + i g ¯ 02 q ¯ ( 0 ) 3 ω 1 * τ 10 * e i ω 1 * τ 10 * θ + E 1 e 2 i ω 1 * τ 10 * θ . W 11 ( θ ) = i g 11 q ( 0 ) ω 1 * τ 10 * e i ω 1 * τ 10 * θ + i g ¯ 11 q ¯ ( 0 ) ω 1 * τ 10 * e i ω 1 * τ 10 * θ + E 2 .

其中 E 1 E 2 可以由如下方程分别决定

( 2 i ω 1 * a 11 a 12 a 13 a 21 e 2 i ω 1 * τ 10 * 2 i ω 1 * a 22 a 23 0 a 32 e 2 i ω 1 * τ 2 * 2 i ω 1 * ) E 1 = 2 ( K 11 K 21 K 31 ) , ( a 11 a 12 a 13 a 21 a 22 a 23 0 a 32 0 ) E 2 = ( K 12 K 22 K 32 ) .

其中:

K 11 = 0 , K 21 = a β ( a + P * ) 2 q 2 q 3 + a β Z * ( a + P * ) 3 q 2 2 , K 31 = a β 1 ( a + P * ) 2 q 2 q 3 e i ω 1 * τ 2 * a β 1 Z * ( a + P * ) 3 q 2 2 e 2 i ω 1 * τ 2 * ;

K 12 = 0 , K 22 = a β ( a + P * ) 2 ( q 2 q 3 ¯ + q 2 ¯ q 3 ) + a β Z * ( a + P * ) 3 2 q 2 q 2 ¯ , K 31 = a β 1 ( a + P * ) 2 ( q 2 ¯ q 3 e i ω 1 * τ 2 * + q 2 q 3 ¯ e i ω 1 * τ 2 * ) a β 1 Z * ( a + P * ) 3 2 q 2 q 2 ¯ .

最后基于 g 20 , g 11 , g 02 , g 21 的表达式可以计算下面的值

c 1 ( 0 ) = i 2 ω 1 * τ 10 * ( g 20 g 11 2 | g 11 | 2 | g 02 | 2 3 ) + g 21 2 , μ 2 = Re { c 1 ( 0 ) } Re { λ ( τ 10 * ) } , β 2 = 2 Re { c 1 ( 0 ) } , T 2 = Im { c 1 ( 0 ) } + μ 2 Im { λ ( τ 10 * ) } ω 1 * τ 10 * .

根据之前的结论,可以得到关于Hopf分支性质的定理如下:

定理3.1:系统(2)在 τ 1 = τ 10 * 处的Hopf分支性质可由 c 1 ( 0 ) , μ 2 , β 2 , T 2 决定,且结论如下:

(i) 如果 μ 2 > 0 ( < 0 ) ,Hopf分支是超临界的(亚临界的),

(ii) 如果 β 2 < 0 ( > 0 ) ,分支周期解是稳定的(不稳定的),

(iii) 如果 T 2 > 0 ( < 0 ) ,系统(2)的分支周期解是增加的(减少的)。

4. 结论

本文建立了一类具有时滞效应的溶解氧-浮游生物动力学新模型,以此为基础对模型的相关动力学性质进行了数学分析与理论推导,解析出模型内平衡点具有局部渐近稳定性和发生Hopf分支的阈值条件,并进一步讨论了Hopf分支的相关性质,给出一些理论结果。本论文的研究结果有利于促进海洋生态动力学模型的长远发展,有助于剖析海洋生态系统中各类生态环境因素的影响驱动性。

致谢

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

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

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

文章引用

郭庆. 一类具有双时滞效应的溶解氧-浮游生物模型的动力学分析
Dynamic Analysis of a Dissolved Oxygen-Plankton Model with Two Time Delays[J]. 应用数学进展, 2017, 06(07): 816-830. http://dx.doi.org/10.12677/AAM.2017.67099

参考文献 (References)

  1. 1. Rana, S., Samanta, S., Bhattacharya, S., et al. (2015) The Effect of Nanoparticles on Plankton Dynamics: A Mathematical Model. BioSystems, 127, 28-41. https://doi.org/10.1016/j.biosystems.2014.11.003

  2. 2. Dai, C., Zhao, M. and Yu, H. (2016) Dynamics Induced by Delay in a Nutrient-Phytoplankton Model with Diffusion. Ecological Complexity, 26, 29-36. https://doi.org/10.1016/j.ecocom.2016.03.001

  3. 3. Edwards, A.M. and Brindley, J. (1999) Zooplankton Mortality and the Dynamical Behaviour of Plankton Population Models. Bulletin of Mathematical Biology, 61, 303-339. https://doi.org/10.1006/bulm.1998.0082

  4. 4. Chattopadhyay, J. and Pal, S. (2002) Viral Infection on Phytoplankton-Zooplankton System—A Mathematical Model. Ecological Modelling, 151, 15-28. https://doi.org/10.1016/S0304-3800(01)00415-X

  5. 5. Mei, L. and Zhang, X. (2012) Existence and Nonexistence of Positive Steady States in Multi-Species Phytoplankton Dynamics. Journal of Differential Equations, 253, 2025-2063. https://doi.org/10.1016/j.jde.2012.06.011

  6. 6. Ruan, S. (1995) The Effect of Delays on Stability and Persistence in Plankton Models. Nonlinear Analysis, 24, 575-585. https://doi.org/10.1016/0362-546X(95)93092-I

  7. 7. Zhao, J., Tian, J.P. and Wei, J. (2016) Minimal Model of Plankton Systems Revisited with Spatial Diffusion and Maturation Delay. Bulletin of Mathematical Biology, 78, 381. https://doi.org/10.1007/s11538-016-0147-3

  8. 8. Chatterjee, A. and Pal, S. (2016) Plankton Nutrient Interaction Model with Effect of Toxin in Presence of Modified Traditional Holling Type II Functional Response. Systems Science & Control Engineering, 4, 20-30. https://doi.org/10.1080/21642583.2015.1136801

  9. 9. Misra, A.K. (2011) Modeling the Depletion of Dissolved Oxygen due to Algal Bloom in a Lake by Taking Holling Type-III Interaction. Applied Mathematics & Computation, 217, 8367-8376. https://doi.org/10.1016/j.amc.2011.03.034

  10. 10. Smith, D.W and Piedrahita, R.H. (1988) The Relation Between Phytoplankton and Dissolved Oxygen in Fish Ponds. Aquaculture, 68, 249-265. https://doi.org/10.1016/0044-8486(88)90357-2

  11. 11. Sekerci, Y. and Petrovskii, S. (2015) Mathematical Modelling of Plankton-Oxygen Dynamics under the Climate Change. Bulletin of Mathematical Biology, 77, 2325-2353. https://doi.org/10.1007/s11538-015-0126-0

  12. 12. Sekerci, Y. and Petrovskii, S. (2015) Mathematical Modelling of Spatiotemporal Dynamics of Oxygen in a Plankton System. Mathematical Modelling of Natural Phenomena, 10, 96-114. https://doi.org/10.1051/mmnp/201510207

  13. 13. Murillo, A.A., Delong, E.F. and Ulloa, O. (2014) Enhanced Metabolic Versatility of Planktonic Sulfur-Oxidizing γ-Proteobacteria in an Oxygen-Deficient Coastal Ecosystem. Frontiers in Marine Science, 1, 18. https://doi.org/10.3389/fmars.2014.00018

  14. 14. Kharel, S., Kumar, S. and Singh, C. (2013) Modelling Effect of the Depleting Dissolved Oxygen on the Existence of Interacting Planktonic Population. Elixir. Applied Mathematics, 55, 12739-12742.

  15. 15. Naik, V.K. and Manjapp, S. (2010) Prediction of Dissolved Oxygen through Mathematical Modeling. International Journal of Environmental Research, 4, 153-160.

  16. 16. Misra, A.K. (2010) Modeling the Dep-letion of Dissolved Oxygen in a Lake Due to Submerged Macrophytes. Nonlinear Analysis: Modelling and Control, 15, 185-198.

  17. 17. Shukla, J.B., Misra, A.K. and Chandra, P. (2008) Mathematical Modeling and Analysis of the Depletion of Dissolved Oxygen in Eutrophied Water Bodies Affected by Organic Pollutants. Nonlinear Analysis Real World Applications, 9, 1851-1865.

  18. 18. Dhar, J. and Baghel, R.S. (2016) Role of Dissolved Oxygen on the Plankton Dynamics in Spatio-Temporal Domain. Modeling Earth Systems & Environment, 2, 6. https://doi.org/10.1007/s40808-015-0061-y

  19. 19. Sharma, A., Sharma, A.K. and Agnihotri, K. (2014) The Dynamic of Plankton-Nutrient Interaction with Delay. Applied Mathematics & Computation, 231, 503-515.

  20. 20. Lian, F.Y. and Xu, Y.T. (2009) Hopf Bifurcation Analysis of a Predator-Prey System with Holling Type IV Functional Response and Time Delay. Applied Mathematics & Computation, 215, 1484-1495.

  21. 21. Zhao, H.Y., Huang, X.X. and Zhang, X.B. (2015) Hopf Bifurcation and Harvesting Control of a Bioeconomic Plankton Model with Delay and Diffusion Terms. Physica A, 421, 300-315.

  22. 22. Das, K. and Ray, S. (2008) Effect of Delay on Nutrient Cycling in Phytop-lankton-Zooplankton Interactions in Estuarine System. Ecological Modelling, 215, 69-76.

  23. 23. Rehim, M., Zhang, Z.Z. and Muhammadhaji, A. (2016) Mathematical Analysis of a Nutrient-Plankton System with Delay. SpringerPlus, 5, 1-22. https://doi.org/10.1186/s40064-016-2435-7

  24. 24. Chen, S.S., Shi, J.P. and Wei, J.J. (2013) Time Delay-Induced Instabilities and Hopf Bifurcations in General Reaction-Diffusion Systems. Journal of Nonlinear Science, 23, 1-38. https://doi.org/10.1007/s00332-012-9138-1

  25. 25. Song, Y.L. and Wei, J.J. (2005) Local Hopf Bifurcation and Global Periodic Solutions in a Delayed Predator-Prey System. Journal of Mathematical Analysis & Applications, 301, 1-21.

  26. 26. Shi, R.X. and Yu, J. (2017) Hopf Bifurcation Analysis of Two Zooplankton-Phytoplankton Model with Two Delays. Chaos Solitons & Fractals, 100, 62-73.

  27. 27. Liu, C., Wang, L.P., Zhang, Q.L. and Yan, Y. (2017) Dynamical Analysis in a Bioeconomic Phytoplankton Zooplankton System with Double Time Delays and Environmental Stochasticity. Physica A Statistical Mechanics & Its Applications, 82, 682-698.

  28. 28. Liu, J. (2015) Hopf Bifurcation Analysis for an SIRS Epidemic Model with Logistic Growth and Delays. Journal of Applied Mathematics & Computing, 50, 1-20.

  29. 29. Hassard, B.D., Kazarinoff, N.D. and Wan, Y. (1981) Theory and Applications of Hopf Bifurcation. Cambridge University Press, Cambridge.

  30. 30. Song, Y.L. and Wei, J.J. (2004) Bifurcation Analysis for Chen’s System with Delayed Feedback and Its Application to Control of Chaos. Chaos Solitons Fractals, 22, 75-91.

期刊菜单