本文考虑一类具有时滞和收获项的离散Nicholson果蝇模型，通过构造适当的Lyapunov函数，得到了系统持久性和全局吸引性的一些充分条件，最后用Matlab进行仿真，验证本文主要结论的正确性。

Nicholson果蝇模型，时滞，持久性，全局吸引性

Zhifeng She, Xinyuan Liao^*, Wei Jin

School of Mathematics and Physics, University of South China, Hengyang Hunan

Received: Jul. 26^th, 2021; accepted: Aug. 23^rd, 2021; published: Aug. 30^th, 2021

In this paper, we consider a discrete Nicholson’s blowflies model with two delays. By constructing suitable Lyapunov functional, a sufficient condition for the permanence and global attractivity of the system is obtained. Numerical simulation shows the feasibility of our main results.

Keywords:Nicholson’s Blowflies Equation Model, Delay, Permanence, Global Attractivity

Copyright © 2021 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/

由于果蝇的生长周期很短，便于进行生物实验。因此，不论是在生物学方面，还是数学方面都相继出现了很多关于果蝇的讨论。20世纪80年代，Gurney [ 1 ] 首次提出了著名的Nicholson果蝇模型，此后，关于Nicholson果蝇模型的研究迅速发展，并已得到了一些很好的研究结果，如Chen [ 2 ] 考虑到外界环境周期性的变化，提出时滞周期型Nicholson苍蝇模型，得到了正周期解存在的一些充分条件；Yang [ 3 ] 研究了具有多时变时滞的广义Nicholson果蝇模型，利用一些微分不等式，证明了该模型的所有正解都指数收敛于零平衡点；Wang [ 4 ] 研究了一类具有多线性收获项的Nicholson果蝇模型正概周期解的存在性和指数收敛性；Berezansky [ 5 ] 等提出了一个带收获双时滞的微分方程模型，对于该模型，在文献 [ 6 ] [ 7 ] [ 8 ] 中已经得到了正平衡点的全局吸引性的相关结果，更多研究结果参见文献 [ 9 ] [ 10 ] [ 11 ] [ 12 ] [ 13 ]。

相对连续模型，离散模型的动力学性态更丰富也更复杂，许多无法求解或理论分析的连续模型往往需要离散模型进行数值模拟。关于离散的Nicholson模型正平衡点的振荡性、全局吸引性和概周期正解的存在性及指数收敛性等动力学性质得到了广泛地讨论 [ 14 ] [ 15 ] [ 16 ]。近年来，由于非自治离散模型能够更准确地刻画系统丰富且复杂的动力学性质，众多学者 [ 17 ] [ 18 ] [ 19 ] 把目光集中在非自治离散模型的研究上。本文基于文献 [ 5 ]，研究了一类具有收获项双时滞的离散非自治Nicholson果蝇模型的持久性和全局吸引性。

本文考虑下列具有双时滞的Nicholson果蝇离散模型

x ( n + 1 ) − x ( n ) = − α ( n ) x ( n ) + β ( n ) x ( n − τ ) e − γ ( n ) x ( n − τ ) − h ( n ) x ( n − σ ) , (2.1)

初值函数为 x ( k ) = φ ( k ) > 0 , k ∈ [ − τ * , 0 ] , τ * = max { τ , σ } (2.2)

其中， α ( n ) 为第n代成虫平均死亡率， τ 为成长期， β ( n ) 为第n代按种群数量计算平均产卵率的最大值， 1 / γ ( n ) 为第n代以最大产卵率繁殖的种群数量， h ( n ) 为第n代的收获率， σ 为成熟期。从可持续经济的角度出发，幼体应该不在捕捞范围内，模型(2.1)考虑到捕捞项，引入了一个具时滞的线性函数 h ( n )     x ( n − σ ) 。

首先，我们假设下列条件(H₁) (H₂)成立

(H₁) { α ( n ) } , { β ( n ) } , { γ ( n ) } , { h ( n ) } 都是有界非负序列，且

0 < α l < α u < 1 ,   0 < β l < β u ,   0 < γ l < γ u ,   0 < h l < h u < 1 ,     α u + h u < 1.

(H₂) 存在两个正常数 d 1 > d 2 > 0 ，使得

d 1 ≥ β u α l γ l e ,   1 γ u ≤ d 2 ≤ β l α u d 1 e − γ u d 1 − h u α u d 1 .

其中，对任意有界序列 { f ( n ) } ，定义 f l = inf k ∈ Z f ( k ) , f u = sup k ∈ Z f ( k ) 。

引理 2.1 [ 20 ] 若 A > 0 , y ( 0 ) > 0 ，且

1) 若 y ( n + 1 ) ≤ A y ( n ) + B ( n ) , n = 1 , 2 , ⋯ ，则对任意正整数 k ≤ n 有，

y ( n ) ≤ A k y ( n − k ) + ∑ i = 0 k − 1 A i B ( n − i − 1 ) ,       n = 1 , 2 , ⋯ .

特别地，当 A < 1 ，B有上确界M时，则

lim n → ∞ sup y ( n ) ≤ M 1 − A .

2) 若 y ( n + 1 ) ≥ A y ( n ) + B ( n ) , n = 1 , 2 , ⋯ ，则对任意正整数 k ≤ n 有，

y ( n ) ≥ A k y ( n − k ) + ∑ i = 0 k − 1 A i B ( n − i − 1 ) ,       n = 1 , 2 , ⋯ .

特别地，当 A < 1 ，B有下确界m时，则

lim n → ∞ sup y ( n ) ≥ m 1 − A .

引理2.2 若 x ( n ) 是方程(2.1)在初值条件(2.2)下的任意解，则

lim n → ∞ sup x ( n ) ≤ M ,     M = β u α l γ l e .

证明：由方程(2.1)可得

x ( n + 1 ) = ( 1 − α ( n ) ) x ( n ) + β ( n ) x ( n − τ ) e − γ ( n ) x ( n − τ ) − h ( n ) x ( n − σ )                         ≤ ( 1 − α ( n ) ) x ( n ) + β u x ( n − τ ) e − γ l x ( n − τ ) .

根据 sup u > 0 u e − γ   u = 1 γ e ，可以得到

x ( n + 1 ) ≤ ( 1 − ∂ ( n ) ) x ( n ) + β u γ l e . (2.3)

对于式(2.3)，由引理2.1可得

lim n → ∞ sup x ( n ) ≤ β u α l γ l e = d e f M .

证毕。

引理2.3 [ 16 ] 令 S 0 = { φ | d 2 < φ ( k ) < d 1 , ∀ k ∈ [ − τ * , 0 ] } ,   τ * = max { τ , σ } ，那么方程(2.1)在初值函数 φ ∈ S 0 下的任意解 x ( k ) 满足 d 2 < x ( k ) < d 1 ,   k ∈ Z + 。

定理3.1 若条件(H₁)，(H₂)成立，那么具有初值函数 φ ∈ S 0 的系统(2.1)是持久的，即存在常数 M , m > 0 ，使系统(2.1)的任意正解 x ( n ) 满足

m ≤ lim n → ∞ inf x ( n ) ≤ lim n → ∞ sup x ( n ) ≤ M .

证明：由引理2.2、引理2.3可以得到

lim n → ∞ sup x ( n ) ≤ M = β u α l γ l e ,   lim n → ∞ inf x ( n ) > 0.

接下来，只需证明存在正常数m使得 lim n → ∞ inf x ( n ) ≥ m ，

令

η = lim n → ∞ inf x ( n ) ,   G = min { g ( η ) , g ( M ) } ,   g ( x ) = x e − γ u x ,

由方程(2.1)和 0 < α u < 1 可以得到

η = lim n → ∞ inf x ( n ) ≥ [ ( x ( 0 ) + β l G 1 + h u ) ( 1 − α u ) n + β l G α u ( 1 + h u ) ] = β l G α u ( 1 + h u ) ,

如果 G = g ( η ) ，那么

η ≥ β l α u ( 1 + h u ) η e − γ u η ,

即

η ≥ 1 γ u ln β l α u ( 1 + h u ) ,

如果 G = g ( M ) ，那么

η ≥ β l α u ( 1 + h u ) M e − γ u M .

则

lim n → ∞ inf x ( n ) ≥ min { 1 γ u ln β l α u ( 1 + h u ) , β l α u ( 1 + h u ) M e − γ u M } = d e f m .

证毕。

定理3.2 若条件(H₁)，(H₂)成立，且 γ l m > 1 成立，那么具有初值函数 φ ∈ S 0 的系统(2.1)是全局吸引的，

即对系统(2.1)的任意解 x ( n ) , y ( n ) 有， lim n → ∞ ( x ( n ) − y ( n ) ) = 0 。

证明：由定理3.1可得

m ≤ lim n → ∞ inf x ( n ) ≤ lim n → ∞ sup x ( n ) ≤ M ,   m ≤ lim n → ∞ inf y ( n ) ≤ lim n → ∞ sup y ( n ) ≤ M .

对任意足够小的正常数 ε ，存在一个正整数 n 0 ，使得对任意 n ≥ n 0 ，

m ≤ x ( n ) ,   y ( n ) ≤ M , (3.1)

由均值定理可以得到

x ( n ) e − x ( n ) − y ( n ) e − y ( n ) = ( 1 − θ ( n ) ) e − θ ( n ) ( x ( n ) − y ( n ) ) ,   x ( n ) < θ ( n ) < y ( n ) . (3.2)

令

w ( n ) = x ( n ) − y ( n ) + ∑ s = n − τ n − 1 β ( s + τ ) ( x ( s ) e − γ ( s + τ ) x ( s ) − y ( s ) e − γ ( s + τ ) y ( s ) ) − ∑ s = n − σ n − 1 β ( s + σ ) ( x ( s ) − y ( s ) )   ,

则

w ( n + 1 ) = ( 1 − α ( n ) ) ( x ( n ) − y ( n ) ) + β ( n ) ( x ( n − τ ) e − γ ( n ) x ( n − τ ) − y ( n − τ ) e − γ ( n ) y ( n − τ ) )     − h ( n ) ( x ( n − σ ) − y ( n − σ ) ) + ∑ s = n − τ + 1 n β ( s + τ ) ( x ( s ) e − γ ( s + τ ) x ( s ) − y ( s ) e − γ ( s + τ ) y ( s ) )     − ∑ s = n − σ + 1 n β ( s + σ ) ( x ( s ) − y ( s ) ) ,

Δ w ( n ) = w ( n + 1 ) − w ( n ) = − α ( n ) ( x ( n ) − y ( n ) ) + β ( n + τ ) ( x ( n ) e − γ ( n + τ ) x ( n ) − y ( n ) e − γ ( n + τ ) y ( n ) )     − h ( n + σ ) ( x ( n ) − y ( n ) ) .

构造非负函数 V ( n ) = w 2 ( n ) ，则

由(3.2)、(3.3)可得

Δ V ( n ) = − ( α ( n ) + h ( n + σ ) ) ( 2 − α ( n ) − h ( n + σ ) ) ( x ( n ) − y ( n ) ) 2     + 2 ( 1 − α ( n ) − h ( n + σ ) ) β ( n + τ ) ( 1 − θ ( n ) γ ( n + τ ) ) e − θ ( n ) γ ( n + τ ) ( x ( n ) − y ( n ) ) 2     + 2 ( α ( n ) + h ( n + σ ) ) ( x ( n ) − y ( n ) ) ∑ s = n − τ n − 1 β ( s + τ ) ( 1 − θ ( s ) γ ( s + τ ) ) e − θ ( s ) γ ( s + τ ) ( x ( s ) − y ( s ) )     + ( β ( n + τ ) ( 1 − θ ( n ) γ ( n + τ ) ) e − θ ( n ) γ ( n + τ ) ( x ( n ) − y ( n ) ) ) 2       + 2 β ( n + τ ) ( 1 − θ ( n ) γ ( n + τ ) ) e − θ ( n ) γ ( n + τ ) ( x ( n ) − y ( n ) )       ⋅ ∑ s = n − τ n − 1 β ( s + τ ) ( 1 − θ ( s ) γ ( s + τ ) ) e − θ ( s ) γ ( s + τ ) ( x ( s ) − y ( s ) )       + 2 ( α ( n ) + h ( n + σ ) ) ( x ( n ) − y ( n ) ) ∑ s = n − σ n − 1 h ( s + σ ) ( x ( s ) − y ( s ) )       − 2 β ( n + τ ) ( 1 − θ ( n ) γ ( n + τ ) ) e − θ ( n ) γ ( n + τ ) ( x ( n ) − y ( n ) ) ⋅ ∑ s = n − σ n − 1 h ( s + σ ) ( x ( s ) − y ( s ) ) . (3.4)

将 2 a b ≤ a 2 + b 2 运用到式(3.4)得

Δ V ( n ) ≤ − ( α ( n ) + h ( n + σ ) ) ( 2 − α ( n ) − h ( n + σ ) ) ( x ( n ) − y ( n ) ) 2                         + 2 ( 1 − α ( n ) − h ( n + σ ) ) β ( n + τ ) ( 1 − θ ( n ) γ ( n + τ ) ) e − θ ( n ) γ ( n + τ ) ( x ( n ) − y ( n ) ) 2                         + ( α ( n ) + h ( n + σ ) ) ∑ s = n − τ n − 1 β ( s + τ ) ( 1 − θ ( s ) γ ( s + τ ) ) e − θ ( s ) γ ( s + τ ) ( x ( n ) − y ( n ) ) 2                         + ( α ( n ) + h ( n + σ ) ) ∑ s = n − τ n − 1 β ( s + τ ) ( 1 − θ ( s ) γ ( s + τ ) ) e − θ ( s ) γ ( s + τ ) ( x ( s ) − y ( s ) ) 2                         + ( β ( n + τ ) ( 1 − θ ( n ) γ ( n + τ ) ) e − θ ( n ) γ ( n + τ ) ( x ( n ) − y ( n ) ) ) 2                           + β ( n + τ ) ( 1 − θ ( n ) γ ( n + τ ) ) e − θ ( n ) γ ( n + τ )                           ⋅ ( ∑ s = n − τ n − 1 β ( s + τ ) ( 1 − θ ( s ) γ ( s + τ ) ) e − θ ( s ) γ ( s + τ ) ( x ( n ) − y ( n ) ) 2                           + ∑ s = n − τ n − 1 β ( s + τ ) ( 1 − θ ( s ) γ ( s + τ ) ) e − θ ( s ) γ ( s + τ ) ( x ( s ) − y ( s ) ) 2 )                           + ( α ( n ) + h ( n + σ ) ) ∑ s = n − σ n − 1 h ( s + σ ) ( x ( n ) − y ( n ) ) 2                           + ( α ( n ) + h ( n + σ ) ) ∑ s = n − σ n − 1 h ( s + σ ) ( x ( s ) − y ( s ) ) 2                           + β ( n + τ ) ( 1 − θ ( n ) γ ( n + τ ) ) e − θ ( n ) γ ( n + τ ) ∑ s = n − σ n − 1 h ( s + σ ) ( x ( n ) − y ( n ) ) 2                           + β ( n + τ ) ( 1 − θ ( n ) γ ( n + τ ) ) e − θ ( n ) γ ( n + τ ) ∑ s = n − σ n − 1 h ( s + σ ) ( x ( s ) − y ( s ) ) 2 . (3.5)

依据 γ l m > 1 ,   α u + h u < 1 ,   max x ∈ [ 1 , + ∞ ) ( 1 − x ) e − x < 0 ，不等式(3.5)可写成

Δ V ( n ) ≤ − ( α ( n ) + h ( n + σ ) ) ( 2 − α ( n ) − h ( n + σ ) ) ( x ( n ) − y ( n ) ) 2                     ≤ − ( α u + h u ) ( 2 − α u − h u ) ( x ( n ) − y ( n ) ) 2                     ≤ − δ ( x ( n ) − y ( n ) ) 2 .

取一个足够小的正常数 δ > 0 使得

( α u + h u ) ( 2 − α u − h u ) ≤ δ ,

则

Δ V ( n ) ≤ − δ ( x ( n ) − y ( n ) ) 2 . (3.6)

将不等式(3.6)左右两边分别从 n 0 + τ 累加到n可得

∑ s = n 0 + τ n ( V ( s + 1 ) − V ( s ) ) ≤ − δ ∑ s = n 0 + τ n ( x ( s ) − x ( s ) ) 2 ,

即

V ( n + 1 ) + δ ∑ s = n 0 + τ n ( x ( s ) − x ( s ) ) 2 ≤ V ( n 0 + τ ) ,

则

∑ s = n 0 + τ n ( x ( s ) − x ( s ) ) 2 ≤ V ( n 0 + τ ) δ .

由(3.1)可知 V ( n 0 + τ ) 是有界的，因此

∑ s = n 0 + τ n ( x ( s ) − x ( s ) ) 2 ≤ V ( n 0 + τ ) δ < + ∞ ,

也就是说

∑ s = n 0 + τ + ∞ ( x ( s ) − x ( s ) ) 2 ≤ V ( n 0 + τ ) δ < + ∞ .

显而易见， lim n → ∞ ( x ( n ) − y ( n ) ) 2 = 0 ，即 lim n → ∞ ( x ( n ) − y ( n ) ) = 0 。

证毕。

考虑下列方程

x ( n + 1 ) − x ( n ) = − ( 0.5 + 0.01 sin ( 2 n ) ) x ( n ) + ( 2.79 + 0.01 sin ( 3 n ) ) x ( n − 1 )                                               ⋅ e − ( 0.82 + 0.01 sin n ) x ( n − 1 ) ) − ( 0.2 + 0.01 sin 2 n ) x ( n − 1 ) (4.1)

初始条件为

x ( − 1 ) = 1.4 ,     x ( 0 ) = 1.6

其中，

α u = 0.51 ,   α l = 0.49 ,       β u = 2.8 ,       β l = 2.78 ,       γ u = 0.83 ,       γ l = 0.81 ,       h u = 0.21 , h l = 0.19 ,     M = 2.25 ,   m = 1.56 ,   α u + h u < 1 ,   γ l m ≈ 1.26 > 1.

x ( n + 1 ) − x ( n ) = − ( 0.5 + 0.02 sin ( 2 n ) ) x ( n ) + ( 2.79 + 0.02 sin ( 3 n ) ) x ( n − 1 )                                               ⋅ e − ( 0.82 + 0.02 sin n ) x ( n − 1 ) ) − ( 0.2 + 0.02 sin 2 n ) x ( n − 1 ) (4.2)

初始条件为

x ( − 1 ) = 1.3 ,     x ( 0 ) = 1.4

x ( n + 1 ) − x ( n ) = − ( 0.51 + 0.02 sin ( 2 n ) ) x ( n ) + ( 2.80 + 0.02 sin ( 3 n ) ) x ( n − 1 )                                               ⋅ e − ( 0.80 + 0.02 sin n ) x ( n − 1 ) ) − ( 0.21 + 0.02 sin 2 n ) x ( n − 1 ) (4.3)

初始条件为

x ( − 1 ) = 1.5 ,     x ( 0 ) = 1.4

x ( n + 1 ) − x ( n ) = − ( 0.51 + 0.01 sin ( 2 n ) ) x ( n ) + ( 2.80 + 0.02 sin ( 3 n ) ) x ( n − 1 )                                               ⋅ e − ( 0.82 + 0.01 sin n ) x ( n − 1 ) ) − ( 0.21 + 0.02 sin 2 n ) x ( n − 1 ) (4.4)

初始条件为

x ( − 1 ) = 1.3 ,     x ( 0 ) = 1.5

显然，上述系统均满足条件(H₁)，(H₂) (H₃)，依据定理3.1和定理3.2，系统(4.1)、(4.2)、(4.3)、(4.4)的解 x ( n ) 是持久且全局吸引的的(如图1)。

图1. 系统(4.1)、(4.2)、(4.3)、(4.4)解 x ( n ) 的动态行为

本文考虑了具有时滞和收获项的离散非自治Nicholson果蝇模型。对于该模型，通过构造适当的Lyapunov函数，我们得到了它的持久性和全局吸引性存在的充分条件，并通过Matlab数值模拟充分验证了结果。

佘智凤,廖新元,金 薇. 一类具双时滞的离散Nicholson果蝇模型的持久性和全局吸引性Permanence and Global Attractivity of a Discrete Nicholson’s Blowflies Model with Two Delays[J]. 应用数学进展, 2021, 10(08): 2888-2896. https://doi.org/10.12677/AAM.2021.108302