基于水环境污染日益严重，本文建立了一类毒素影响下鱼类具有常数投放的捕食–食饵模型，研究了系统平衡点的存在性与稳定性，利用后继函数法证明了系统阶一周期解存在性，其次利用类Poincaré准则，得到了阶一周期解稳定的条件，最后通过数值模拟验证结论的重要性。 In this paper, a predator-prey model of fish with constant release under the influence of toxins is established based on the increasing pollution of water environment. The existence and stability of the equilibrium point of the system are studied, and the existence of the order-1 periodic solution is proved by using the method of successor function. Secondly, the condition of the stability of the or-der-1 periodic solution is obtained by using the analogue of Poincaré’s criterion. Finally, the im-portance of the conclusion is verified by numerical simulation.
基于水环境污染日益严重，本文建立了一类毒素影响下鱼类具有常数投放的捕食–食饵模型，研究了系统平衡点的存在性与稳定性，利用后继函数法证明了系统阶一周期解存在性，其次利用类Poincaré准则，得到了阶一周期解稳定的条件，最后通过数值模拟验证结论的重要性。
毒素，常数投放，阶一周期解，类Poincaré准则，数值模拟
Yan Zhang*, Meng Zhang
School of Science, Beijing University of Civil Engineering and Architecture, Beijing
Received: Jul. 24th, 2022; accepted: Aug. 17th, 2022; published: Aug. 29th, 2022
In this paper, a predator-prey model of fish with constant release under the influence of toxins is established based on the increasing pollution of water environment. The existence and stability of the equilibrium point of the system are studied, and the existence of the order-1 periodic solution is proved by using the method of successor function. Secondly, the condition of the stability of the order-1 periodic solution is obtained by using the analogue of Poincaré’s criterion. Finally, the importance of the conclusion is verified by numerical simulation.
Keywords:Toxin, Constant Release, Order-1 Periodic Solution, The Analogue of Poincaré’s Criterion, Numerical Simulation
Copyright © 2022 by author(s) and beplay安卓登录
This work is licensed under the Creative Commons Attribution International License (CC BY 4.0).
http://creativecommons.org/licenses/by/4.0/
随着现代化工业的迅速发展，在促进国民经济发展的同时，造成的环境污染问题也是不容忽视的，这严重影响了空气、土壤以及水源。环境污染会对生态系统造成直接的影响和破坏。同时环境污染问题已是人类面临的重大课题之一。
近些年来，许多学者发现脉冲微分方程的理论和方法在生命科学领域中有很好的应用前景。例如陈 [
本文考虑水环境污染导致鱼类的减少，对捕食者(鱼类)进行常量投放，可以防止浮游植物(食饵)占据大量生存空间，同时利于鱼类的可持续发展。因此建立了在毒素影响下捕食者具有常数投放和状态反馈脉冲控制模型。
本文模拟了在毒素影响下捕食者具有常数投放和HollingΙ功能反应的捕食系统模型，来反映鱼类的养殖与捕捞整个过程。模型如下
{ d x d t = x ( r − x ) − k x y d y d t = h x y − f y + u Δ x = − α x Δ y = − β y } y < m } y = m . (1)
其中x与y分别为浮游植物与鱼类在t时刻的数量，r表示浮游植物的生长率(出生率减去死亡率)，k为鱼类对浮游植物的消耗率，h为鱼类的捕食率，f表示为毒素影响下鱼类的死亡率，u表示鱼类的投放率。当鱼类的数量达到m时，收获鱼类，为了避免浮游植物大量繁殖，也要捕捞一定量的浮游植物。 0 ≤ α ≤ 1 为收获浮游植物的比率， 0 ≤ β ≤ 1 为收获鱼类的比例。考虑种群实际情况，以上参数均为正常数且两种群具有正初始值情况研究将在 R 2 + = { ( x , y ) | x ≥ 0 , y ≥ 0 } 内进行。
当不考虑脉冲时，系统(1)退化为
{ d x d t = x ( r − x ) − k x y d y d t = h x y − f y + u . (2)
求方程组
{ d x d t = 0 d y d t = 0 . (3)
当 y = 0 时，x无解。则在横坐标轴上不存在奇点。
当 x = 0 时， y = u f ，故存在奇点 E 1 ( 0, u f ) 。
当 x > 0 , y > 0 时，有 r − x − k y = 0 ，将其代入式子 h x y − f y + u = 0 中，得
− h k y 2 + ( h r − f ) y + u = 0 . (4)
由韦达定理可得
y 1 y * = u − h k < 0 ， (5)
所以必有一正根 y * ，由求根公式可得
y * = f − h r − ( h r − f ) 2 + 4 u k h − 2 h k ， (6)
x * = r − k y * = r + f − h r − ( h r − f ) 2 + 4 u k h 2 h . (7)
由于当 x * > 0 时，有
r + f − h r − ( h r − f ) 2 + 4 u k h 2 h > 0 ， (8)
2 h r + f − h r > ( h r − f ) 2 + 4 u k h ， (9)
( 2 h r + f − h r ) 2 > ( h r − f ) 2 + 4 u k h ， (10)
4 h 2 r 2 + 4 h r ( f − h r ) > 4 u k h ， (11)
可以得出
f r − u k > 0 . (12)
则在第一象限存在一个正平衡点 E 2 ( x * , y * ) 。根据实际情况，下面重点讨论系统存在正平衡点的情况。
系统(2)的Jocabi矩阵为
J = [ r − 2 x − k y − k x h y h x − f ] ， (13)
对于 E 1 ( 0 , u f ) 的特征方程
| r − k u f − λ ¯ 0 h u f − f − λ ¯ | = 0 . (14)
因为 f r − k u > 0 ，所以 λ ¯ 1 λ ¯ 2 < 0 ，则 E 1 ( 0 , u f ) 是鞍点。
通过计算，得到 E 2 ( x * , y * ) 的特征方程为
| − x * − λ ¯ − k x * h ( r − x * ) k h x * − f − λ ¯ | = 0 ， (15)
λ ¯ 1 λ ¯ 2 = det ( J ) = ( f − h x * ) x * + h ( r − x * ) x * = ( u y * + h k y * ) x * > 0 ， (16)
λ ¯ 1 + λ ¯ 2 = t r ( J ) = − x * + h x * − f ， (17)
因为 h x y − f y + u = 0 ，则 h x − f = − u y ，有 − x * + h x * − f < 0 ，即
λ ¯ 1 + λ ¯ 2 < 0 . (18)
所以 E 2 ( x * , y * ) 是稳定的结点或者焦点。
系统(2)的轨线图如图1所示。
图1. 系统(2)的轨线图
定理1：当 λ = k h 时，系统(2)在第一象限是全局渐近稳定。
证：建立一个Lyapunov函数
V ( x , y ) = [ ( x − x * ) − x * ln ( x x * ) ] + λ [ ( y − y * ) − y * ln ( y y * ) ] ， (19)
d V d t = x − x * x d x d t + y − y * y d y d t = ( x − x * ) ( r − x − k y ) + λ ( y − y * ) ( h x − f + u y ) ， (20)
考虑到平衡点方程
{ r − x * − k y * = 0 h x * y * − f y * + u = 0 ， (21)
d V d t = ( x − x * ) [ ( x * − x ) − k ( y − y * ) ] + λ ( y − y * ) [ h ( x − x * ) + u ( 1 y − 1 y * ) ] ， (22)
当 λ = k h 时，显然 d V d t < 0 ，所以系统(2)在第一象限是全局渐近稳定的。
{ d x d t = P ( x , y ) d y d t = Q ( x , y ) . (23)
令 B ( x , y ) = 1 x y ，则有
∂ ( B P ) ∂ x + ∂ ( B Q ) ∂ y = − 1 y − u x y 2 < 0 ， (24)
则无脉冲系统在第一象限不存在极限环。
根据
d x d t = 0 , d y d t = 0 ， (25)
得垂直等倾线
x = 0 , y = r − x k ， (26)
与水平等倾线
y = u f − h x . (27)
根据环域定理，构造外境界线，过点C作x轴平行线交水平等倾线于点A，过点A作AB垂直于x轴交x轴于点B。
对于OC，由于y轴是一条轨线，没有任何的轨线与它相交。
对于OB，当 y = 0 时，有 d y d t | y = 0 = u > 0 ，故当轨线与x轴相遇时，均从它的下面穿向上面。
对于AC，由于 0 < x < x A 时，有 d y d t | y = y A = h x y − f y + u < 0 ，故轨线与直线AC相遇时，从上面穿向下面。
可知 B ( f r − u k h r , r k ) ，当 f r − u k h r > r 时，可知B点在R点右侧。
对于AB，由于 d x d t | x = f r − u k h r = x ( r − x ) − k x y < 0 ，故当轨线与直线AB相遇时，均从它的右面穿向左面。
这样就围成了区域G的外境界线(如图2所示)。且G内除E2外无其他奇点，G的边界奇点E1是鞍点，所以从E2出发的解一致有界。
图2. 系统(2)的有界性
本节利用文献 [
定义1：系统为半连续动力系统，则
脉冲集为 M = { ( x , y ) ∈ R 2 + | y = m , x ≥ 0 } ，
脉冲映射为 φ ( x , y ) ∈ M → ( ( ( 1 − α ) x , ( 1 − β ) y ) ) ∈ R 2 + ，
相集为 N = φ ( M ) = { ( x , y ) ∈ R 2 + | y = ( 1 − β ) m , x ≥ 0 } 。
定义2：假设系统有一条从点 D ( x D , y D ) ∈ N 出发的轨线，到达 M 后与其相交于点 D ( x D 1 , h ) ∈ M ，经过脉冲作用后，跳跃至点 D 2 ( x D 2 , y D 2 ) ∈ N ，称 D 2 为 D 的后继点，且记 D 的后继函数为 f ( D ) = x D 2 − x D 。
定义3：对于系统来说，如果存在一点 D ∈ N 和 t 0 > 0 ，使得 f ( D , t 0 ) = D 1 ∈ M 和 φ ( f ( D , t 0 ) ) = D ∈ N 成立，那么 f ( D , t 0 ) 称为系统的阶一周期解。如图3所示。
图3. 阶一周期解
引理1：对于系统来说，如果存在两点 E ∈ N 和 F ∈ N ，使得 f ( E ) f ( F ) < 0 ，那么在点 E 和点 F 之间一定存在一点 D ∈ N ，使得 f ( D ) = 0 。
当 m = y * 时， Δ y = y + − y * = − β y * , y = ( 1 − β ) y * , 0 < β < 1 。相集N与垂直等倾线相交于点A，过A微分方程的轨线与脉冲集与一点 A 1 ，微分方程的轨线段记为 π ( A , A 1 ) ， A 1 经对应参数脉冲映射 φ 到相集上一点 A ′ ，映射线段记为 φ ( A 1 , A ′ ) ， π ( A , A 1 ) + φ ( A 1 , A ′ ) 称为脉冲参数为 − α 的阶一轨线段， A ′ 也称为A脉冲参数为 − α 的阶一后继点。设点A坐标为 α ， A ′ 点坐标为 a ′ ，则点A的后继函数 F − α 1 ( A ) = a ′ − a < 0 。
在相集N存在一点B，使B点坐标为 ( τ , ( 1 − β ) y * ) ， τ 为充分小的正数。过B微分方程的轨线与脉冲集与一点 B 1 ，微分方程的轨线段记为 π ( B , B 1 ) ， B 1 经对应参数脉冲映射 φ 到相集上一点 B ′ ，映射线段记为 φ ( B 1 , B ′ ) ， π ( B , B 1 ) + φ ( B 1 , B ′ ) 称为脉冲参数为 − α 的阶一轨线段， B ′ 也称为B脉冲参数为 − α 的阶一后继点。设点B坐标为b， B ′ 点坐标为 b ′ ，则点B的后继函数 F − α 1 ( B ) = b ′ − b > 0 (如图4所示)。
图4. 阶一周期解的存在性
则点A和点B之间必存在一点D，设D的后继点为 D ′ ，设D点坐标为d， D ′ 点坐标为 d ′ ，则点D的后继函数为 F ′ − α ( D ) = d ′ − d = 0 。这就证明了系统(1)存在阶1周期解 D D ′ (如图3所示)。
已经知道A2是A0的后继点，且A2在A0的右侧，那么从A2出发的轨线将会到达脉冲集，并形成交点A3，在脉冲作用后，轨线到达相集上的点A4，A4是A2的后继点，且A4位于A2的右侧。接着从A4出发的轨线也将与脉冲集相交于点A5，在脉冲作用后，轨线又达到相集上的点A6，显然A6是A4的后继点，且A6位于A4的右侧。以此类推，可以得到两个序列
x A 1 , x A 3 , x A 5 , x A 7 , ⋯ , x A 2 n − 1 , ⋯ ∈ A 1 B 1 ¯ ， (28)
x A 0 , x A 2 , x A 4 , x A 6 , ⋯ , x A 2 n , ⋯ ∈ A 0 B 0 ¯ . (29)
同理，还可以得到轨线从B0出发的两个序列
x B 1 , x B 3 , x B 5 , x B 7 , ⋯ , x B 2 n − 1 , ⋯ ∈ A 1 B 1 ¯ ， (30)
x B 0 , x B 2 , x B 4 , x B 6 , ⋯ , x B 2 n , ⋯ ∈ A 0 B 0 ¯ . (31)
根据这两列序列，可以看出线段 A 1 B 1 ¯ 通过脉冲作用被映射到线段 A 0 B 0 ¯ ，因此有 ( x A 1 , x B 1 ) ⊂ ( x A 0 , x B 0 ) ， d ( A 1 , B 1 ) < d ( A 0 , B 0 ) 。根据区间套定理，可以得到系统(1)在阶一周期解存在时是唯一的，且有 lim n → ∞ x A 2 n = x C ， lim n → ∞ x B 2 n = x C 成立(如图5所示)。
图5. 阶一周期解的唯一性
定理2：当 α h x 1 m > β u 1 − β 时，系统(1)过 A 点的阶1周期解轨道渐近稳定。
证：不妨设 A ( x 1 , m ) 点为脉冲点， A + ( ( 1 − α ) x 1 , ( 1 − β ) m ) 为A点的后继点。令
P ( x , y ) = x ( r − x ) − k x y ， Q ( x , y ) = h x y − f y + u ， (32)
α ( x , y ) = − α x ， β ( x , y ) = − β y ， Φ ( x , y ) = y − m ， (33)
计算得
∂ P ∂ x = r − 2 x − k y ， ∂ Q ∂ y = h x − f ， ∂ α ∂ x = − α ， ∂ β ∂ y = − β ， ∂ α ∂ y = ∂ β ∂ x = ∂ Φ ∂ x = 0 ， ∂ Φ ∂ y = 1 ， (34)
Δ 1 = P + ( ∂ β ∂ y ∂ Φ ∂ x − ∂ β ∂ x ∂ Φ ∂ y + ∂ Φ ∂ x ) + Q + ( ∂ α ∂ x ∂ Φ ∂ y − ∂ α ∂ y ∂ Φ ∂ x + ∂ Φ ∂ y ) P ( ∂ Φ ∂ x ) + Q ( ∂ Φ ∂ y ) = ( 1 − α ) Q ( A + ) Q ( A ) = ( 1 − α ) [ h ( 1 − α ) x 1 ( 1 − β ) m − f ( 1 − β ) m + u ] h x 1 m − f m + u . (35)
∫ 0 T ( ∂ P ∂ x + ∂ Q ∂ y ) d t = ∫ 0 T ( r − 2 x − k y + h x − f ) d t = ∫ ( 1 − α ) x 1 x 1 1 x d x + ∫ ( 1 − β ) m 1 m 1 1 y d y + ∫ 0 T ( − u y − x ) d t = ln 1 1 − α + ln 1 1 − β + ∫ 0 T ( − u y − x ) d t ， (36)
μ 2 = Δ 1 exp { ∫ 0 T ( ∂ P ∂ x + ∂ Q ∂ y ) d t } = h ( 1 − α ) x 1 m − f m + u 1 − β h x 1 m − f m + u exp ∫ 0 T ( − u y − x ) d t . (37)
当 α h x 1 m > β u 1 − β 时， | μ 2 | < 1 。由类Poincaré准则可知，系统(1)过A点的阶一周期解是轨道渐近稳定的。
基于前面对系统(1)的理论分析，通过数值模拟来验证污染环境下投放不同数量的捕食者(鱼类)对该系统的影响。当系统存在正平衡点时，取参数值为 r = 0.9 , k = 0.25 , h = 0.3 , f = 0.5 , u = 0.1 ，此时，系统(1)的捕食者与食饵的时间序列图如图6、图7所示。当我们保留其他参数不变，使 u = 0.5 ，系统(1)的捕食者与食饵的时间序列图如图8、图9所示。
可以说明当污染环境下，当我们投放少量的捕食者(鱼类)，食饵(浮游植物)占据了更多的生存空间，使鱼类的持续生存受到威胁，使收获量减小。但当我们加大投放量时，浮游动物适量减小，而鱼类的数量和收获量大幅度增加。我们通过人工对投放鱼类进行合理控制，这样既增加了资源，也能够维持渔业的正常生存。
图6. u = 0.1 ，x的时间序列图
图7. u = 0.1 ，y的时间序列图
图8. u = 0.5 ，x的时间序列图
图9. u = 0.5 ，y的时间序列图
由此可知，人工投放资源是维护生态平衡的重要手段。同时我们更应该为保护环境付诸行动，合理利用自然资源，防止环境的污染和破坏，这样才能使自然环境、人文环境、经济环境共同平衡可持续发展，扩大有用资源的再生产，保障社会的稳定可持续发展。
本研究由北京建筑大学研究生创新项目资助，项目号：PG2022141。
张 琰,张 蒙. 毒素影响下捕食种群具有常数投放率的捕食–食饵模型的动力学研究Dynamical Research of a Predator-Prey Model of Predator Population with Constant Release Rate under Toxin Influence[J]. 应用数学进展, 2022, 11(08): 6044-6054. https://doi.org/10.12677/AAM.2022.118636
https://doi.org/10.1142/S1793524512600066
https://doi.org/10.1016/j.apm.2008.06.008
https://doi.org/10.19943/j.2095-3070.jmmia.2021.01.01