
Citation: Xin Long. Novel stability criteria on a patch structure Nicholson's blowflies model with multiple pairs of time-varying delays[J]. AIMS Mathematics, 2020, 5(6): 7387-7401. doi: 10.3934/math.2020473
[1] | Zhibo Cheng, Lisha Lv, Jie Liu . Positive periodic solution of first-order neutral differential equation with infinite distributed delay and applications. AIMS Mathematics, 2020, 5(6): 7372-7386. doi: 10.3934/math.2020472 |
[2] | Wentao Wang . Positive almost periodic solution for competitive and cooperative Nicholson's blowflies system. AIMS Mathematics, 2024, 9(5): 10638-10658. doi: 10.3934/math.2024519 |
[3] | Li Wan, Qinghua Zhou, Hongbo Fu, Qunjiao Zhang . Exponential stability of Hopfield neural networks of neutral type with multiple time-varying delays. AIMS Mathematics, 2021, 6(8): 8030-8043. doi: 10.3934/math.2021466 |
[4] | Jenjira Thipcha, Presarin Tangsiridamrong, Thongchai Botmart, Boonyachat Meesuptong, M. Syed Ali, Pantiwa Srisilp, Kanit Mukdasai . Robust stability and passivity analysis for discrete-time neural networks with mixed time-varying delays via a new summation inequality. AIMS Mathematics, 2023, 8(2): 4973-5006. doi: 10.3934/math.2023249 |
[5] | Nayika Samorn, Kanit Mukdasai, Issaraporn Khonchaiyaphum . Analysis of finite-time stability in genetic regulatory networks with interval time-varying delays and leakage delay effects. AIMS Mathematics, 2024, 9(9): 25028-25048. doi: 10.3934/math.20241220 |
[6] | Er-yong Cong, Li Zhu, Xian Zhang . Global exponential synchronization of discrete-time high-order BAM neural networks with multiple time-varying delays. AIMS Mathematics, 2024, 9(12): 33632-33648. doi: 10.3934/math.20241605 |
[7] | R. Sriraman, P. Vignesh, V. C. Amritha, G. Rachakit, Prasanalakshmi Balaji . Direct quaternion method-based stability criteria for quaternion-valued Takagi-Sugeno fuzzy BAM delayed neural networks using quaternion-valued Wirtinger-based integral inequality. AIMS Mathematics, 2023, 8(5): 10486-10512. doi: 10.3934/math.2023532 |
[8] | Xingyue Liu, Kaibo Shi . Further results on stability analysis of time-varying delay systems via novel integral inequalities and improved Lyapunov-Krasovskii functionals. AIMS Mathematics, 2022, 7(2): 1873-1895. doi: 10.3934/math.2022108 |
[9] | Wenlong Xue, Zhenghong Jin, Yufeng Tian . Stability analysis of delayed neural networks via compound-parameter -based integral inequality. AIMS Mathematics, 2024, 9(7): 19345-19360. doi: 10.3934/math.2024942 |
[10] | Junkang Tian, Zerong Ren, Yanmin Liu . New stability criteria for systems with an interval time-varying delay. AIMS Mathematics, 2023, 8(1): 1139-1153. doi: 10.3934/math.2023057 |
As we all know, dynamic model analysis has developed in many disciplines such as biology, economic sciences, and engineer, etc [1,2,3,4,5,6,7,8,9,10]. Especially, using mathematical biological model to study the population dynamic behaviors plays an important role in maintaining the function of ecosystem [11,12,13], guiding production practice [14,15,16], designing artificial ecosystem [17,18] and promoting the development of information technology [19,20]. In order to better describe the influence of delays in mathematical biological model, the Functional Differential Equations (FDEs) have been applied extensively to the study of population dynamics [21], which can be traced back to the 1920s, when Volttera [22] studied the predator-prey model. In the past 100 years, the theory of population dynamics has made remarkable progress with main results scattered in numerous research papers [23,24,25,26,27,28,29,30,31,32,33,34,35,36]. In particular, there have been plenty of papers written about the dynamical characteristic analysis of FDEs, which include global stability [37,38,39,40], periodicity [19,20,27,32], almost periodicity [41,42,43], Hopf bifurcation [14], boundedness [16] and synchronization [17]. Just as pointed out by Berezansky and Braverman [18], the following delayed differential equation
x′(t)=m∑j=1Fj(t,x(t−τ1(t)),⋯,x(t−τl(t)))−G(t,x(t)), t≥t0, | (1.1) |
can be used to characterise the dynamics of many population models. Here m and l are positive integers, G represents the instantaneous mortality rate, and each Fj describes the feedback control depending on the values of the stable variable with respective delays τ1(t),τ2(t),⋯,τl(t). Obviously, (1.1) includes the modified Nicholson's blowflies equation [23,24]
x′(t)=−α(t)x(t)+m∑j=1βj(t)x(t−gj(t))e−γj(t)x(t−hj(t)), t≥t0, | (1.2) |
which in the case gj(t)≡hj(t) coincides with the classical models [25,26,27,28,29,30,31,32,33]. In (1.2), α(t),βj(t), gj(t),hj(t) and γj(t) are all continuous and nonnegative functions, α(t) is the death rate of the population which depends on time t and the current population level x(t), βj(t)x(t−gj(t))e−γj(t)x(t−hj(t)) is the time-dependent birth function which involves maturation delay gj(t) and incubation delay hj(t), and gets the reproduce at its maximum rate 1γj(t), and j∈I:={1,2,⋯,m}.
It should be noted that, in most of the aforementioned works [23,24,25,26,27,28,29,30,31,32,33], the per capita daily adult mortality term α(t) is uniformly positive (that means, there is a positive constant α− such that α(t)≥α− for all t≥t0). However, throughout the time of some seasons, the death or harvest rate may be less or greater than the birth rate in nature [33,34,35], it is more reasonable to study Nicholson's blowflies models without restriction on uniform positiveness of the per capita death rate. On the other hand, it is of great practical significance to investigate the dynamic behaviors of the patch structure Nicholson's blowflies model [36,37,38]. In recent years, the dynamic behaviors of the following modified Nicholson's blowflies model with patch structure
x′i(t)=−αii(t)xi(t)+n∑j=1,j≠iαij(t)xj(t)+m∑j=1βij(t)xi(t−gij(t))e−γij(t)xi(t−hij(t)), t≥t0, i∈Q:={1,2,⋯,n}, | (1.3) |
have aroused the current research interest, and some useful results have been obtained in the existing papers [38,39,40,41,42]. In particular, without restriction on uniform positiveness of αii(t)(i∈Q), Xu, Cao and Guo [39] established some new criteria for the global attractivity of model (1.3) with the following admissible initial conditions:
xi(t0+θ)=φi(θ), θ∈[−ri,0], φ∈C0+={φ∈C+|φi(0)>0,i∈Q}, | (1.4) |
where ri=maxj∈I{supt∈[t0,+∞) gij(t),supt∈[t0,+∞)hij(t)}, C+=n∏i=1C([−ri,0],[0,+∞)). Moreover, the authors in [39] obtained the main result as follows.
Theorem 1.1. Suppose that
(W1) limt→+∞∫t0α(s)ds=+∞, where α(t)=mini∈Q{αii(t)},
and
(W2) 1>lim supt→+∞[n∑j=1,j≠iαij(t)αii(t)+m∑j=1βij(t)αii(t)], i∈Q,
are satisfied. Then the zero equilibrium point in system (1.3) is globally generalized exponentially stable, i.e., there exist two constants T∗>max{t0,0} and M(φ)>0 such that every solution x(t)=x(t;t0,φ) of (1.3) and (1.4) agrees with that
xi(t)≤M(φ)e−λ∫tT∗α(s)ds for all t>T∗ and i∈Q. |
Unfortunately, there exist two mistakes in the proof of the main results in [39]. The first mistake is on lines 1-2 of page 3. More precisely, let xi(t)=xi(t;t0,φ), yi(t)=maxt0−ri≤s≤txi(s) and y(t)=maxi∈Qyi(t), the inequality
xi(t)≤N(t,i):=∥φ∥+∫tt0[n∑j=1,j≠iαij(v)+m∑j=1βij(v)]y(v)dv, | (1.5) |
can not lead to
y(t)=maxi∈Q{maxt0−ri≤s≤txi(s)}≤N(t,i), | (1.6) |
where t∈[t0−ri, η(φ)), i∈Q and ∥φ∥=maxi∈Q{maxt∈[−ri,0]|φi(t)|}. The same mistake also appeared in [40] [see, p. 4], [41] [see, p. 5] and [43] [see, p. 4]. The logical error for the above mistake is that the left side of (1.6) is the whole maximum value for all i, but the right side of (1.6) is only for a fixed i. Fortunately, this error will be corrected in the Lemma 2.1 of this present paper. The second mistake in literature [39] is on line 26 of page 3 that the inequality
y(t)≤‖y(t∗0)‖∞e−∫tt∗0α(v)dv+(1−ξ)y(t), 0<ξ<1, |
which is absurd. Actually, by taking ∥φ∥=100 and letting t→+∞, one can find that
100≤y(t)≤(1−ξ)y(t)<y(t), |
is a contradiction.
Furthermore, we consider the following concrete example
{x′1(t)=−tx1(t)+t3x2(t) +t4x1(t−50|sint|)e−(1+cos2t)x1(t−46|sint|) +t5x1(t−39|sint|)e−(1+sin2t)x1(t−46|sint|), x′2(t)=−(t+1)x2(t)+(t4+14)x1(t) +(t5+15)x2(t−50|cost|)e−(1+sint)x2(t−49|cost|) +(t6+16)x2(t−38|cost|)e−(1+cost)x2(t−39|cost|), | (1.7) |
and its numerical simulations in Figure 1 to illustrate the above mistakes. Obviously, the assumptions (W1) and (W2) are satisfied in (1.7). Therefore, from Theorem 1.1, we have that the zero equilibrium point of system (1.7) is globally generalized exponentially stable. However, the numerical solutions with three different initial values shown in Figure 1 reveal that three numerical trajectories of (1.7) are convergent to 0 as t→+∞, but they are not generalized exponential convergence to 0. This makes us doubt whether the conclusions of Theorem 1.1 are correct.
Based on the above theoretical analysis and numerical simulations, two problems naturally arise. One is whether the assumptions (W1) and (W2) can guarantee that every solution of (1.3) and (1.4) is convergent to 0 as t→+∞. The other is what kind of conditions can ensure the global generalized exponential convergence of the zero equilibrium point of (1.3).
Regarding the above discussions, in this manuscript, we first establish the global convergence of model (1.3) under the original conditions (W1) and (W2). It is worth noting that, two or more delays appearing in the same time-dependent birth function will cause complex dynamic behaviors, even chaotic oscillations. It is impossible to establish the global exponential convergence of (1.1) without appropriate restrictions on the distinctive delays and some corresponding examples are also given in [18]. Therefore, we add a new delay-dependent assumption to obtain some new criteria for the generalized exponential convergence and global exponential asymptotic stability on the zero equilibrium point of (1.3). In a nutshell, our results not only correct the errors in the existing literature [39,40,41,43], but also improve and complement the existing conclusions in the recent publications [23,24,26,32,39], and the effectiveness is demonstrated by a numerical example.
Hereafter, for i∈Q,j∈I, we assume that αii,γij∈C([t0, +∞),(0,+∞)) and αij(i≠j),βij,gij,hij∈C([t0,+∞),[0,+∞)). Label x(t;t0,φ) for a solution of (1.3) associating to (1.4), and [t0, η(φ)) be the maximal existence right-interval.
Lemma 2.1. Every solution x(t)=x(t;t0,φ) is positive on [t0,+∞).
Proof. According to (1.3), (1.4) and Theorem 5.2.1 in [44] [see p. 81], one can find that xi(t)≥0 for all t∈[t0,η(φ)) and i∈Q. Now, we reveal that η(φ)=+∞. For all t∈[t0,η(φ)) and i∈Q, define yi(t)=maxt0−ri≤s≤txi(s) and y(t)=maxi∈Qyi(t), we gain
x′i(s)≤n∑j=1,j≠iαij(s)xj(s)+m∑j=1βij(s)xi(s−gij(s))≤[n∑j=1,j≠iαij(s)+m∑j=1βij(s)]y(s) |
and
xi(s)≤∥φ∥+∫st0[n∑j=1,j≠iαij(v)+m∑j=1βij(v)]y(v)dv≤∥φ∥+∫tt0[n∑j=1,j≠iαij(v)+m∑j=1βij(v)]y(v)dv ≤∥φ∥+∫tt0maxi∈Q[n∑j=1,j≠iαij(v)+m∑j=1βij(v)]y(v)dv, |
where s∈[t0, t] and i∈Q. On combining this with the fact that
xi(s)≤∥φ∥ for all s∈[t0−ri,t0] and i∈Q, |
we deduce that
y(t)≤∥φ∥+∫tt0maxi∈Q[n∑j=1,j≠iαij(v)+m∑j=1βij(v)]y(v)dv for all t∈[t0, η(φ)). |
Hence, by the Gronwall-Bellman inequality, we obtain
xi(t)≤yi(t)≤y(t)≤∥φ∥e∫tt0maxi∈Q[n∑j=1,j≠iαij(v)+m∑j=1βij(v)]dv for all t∈[t0, η(φ)) and i∈Q. |
This and [45] [see Theorem 2.3.1] suggest that η(φ)=+∞, and xi(t)≥0 for all t∈[t0,+∞) and i∈Q. Note that xi(t0)=φi(0)>0, we get
xi(t)=xi(t0)e−∫tt0αii(v)dv+∫tt0e−∫tsαii(v)dv[n∑j=1,j≠iαij(s)xj(s)+m∑j=1βij(s)xi(s−gij(s))e−γij(s)xi(s−hij(s))]ds>0 for all t∈[t0, +∞) and i∈Q, |
which evidences Lemma 2.1.
Theorem 3.1. Under (W1) and (W2), every solution x(t)=x(t;t0,φ) is convergent to 0 as t→+∞.
Proof. First, we prove that every solution x(t)=x(t;t0,φ) is bounded on [t0,+∞). By (W2), we have
σ:=lim supt→+∞(n∑j=1,j≠iαij(t)αii(t)+m∑j=1βij(t)αii(t))<1, |
and hence, 1−σ2>0. Therefore, there exists T∗>t0−ri such that
n∑j=1,j≠iαij(t)αii(t)+m∑j=1βij(t)αii(t)<σ+1−σ2=1+σ2<1, where t∈[T∗,+∞) and i∈Q. | (3.1) |
Clearly,
xi(t)=xi(t;t0,φ)<M(φ):=n∑i=1maxt∈[t0−ri,T∗]xi(t)+1 forall t∈[t0−ri,T∗] and i∈Q. | (3.2) |
Furthermore, we state that
xi(t)<M(φ) forall t∈(T∗,+∞) and i∈Q. |
Otherwise, there exist i∗∈Q and T∗∗>T∗ such that
xi∗(T∗∗)=M(φ), xj(t)<M(φ) forall t∈[t0−rj, T∗∗) and j∈Q. |
Then
0≤x′i∗(T∗∗)=−αi∗i∗(T∗∗)xi∗(T∗∗)+n∑j=1,j≠i∗αi∗j(T∗∗)xj(T∗∗) +m∑j=1βi∗j(T∗∗)xi∗(T∗∗−gi∗j(T∗∗))e−γi∗j(T∗∗)xi∗(T∗∗−hi∗j(T∗∗))≤M(φ)αi∗i∗(T∗∗)[−1+n∑j=1,j≠i∗αi∗j(T∗∗)αi∗i∗(T∗∗)+m∑j=1βi∗j(T∗∗)αi∗i∗(T∗∗)]<0, |
which is a contradiction and proves that x(t) is bounded on [t0,+∞).
Now, it is sufficient to show that u=maxi∈Qlim supt→+∞xi(t)=0.
For any ε>0, one can choose Λ>T∗ such that
xi(t)<u+ε forall t∈[Λ−ri, +∞) and i∈Q. |
Thus,
xi(t)=xi(t0)e−∫tt0αii(v)dv+∫tt0e−∫tsαii(v)dv[n∑j=1,j≠iαij(s)xj(s)+m∑j=1βij(s)xi(s−gij(s))e−γij(s)xi(s−hij(s))]ds<xi(t0)e−∫tt0αii(v)dv+(u+ε)∫tt0e−∫tsαii(v)dv[n∑j=1,j≠iαij(s)+m∑j=1βij(s)]ds<xi(t0)e−∫tt0αii(v)dv+(u+ε)1+σ2∫tt0e−∫tsαii(v)dvαii(s)ds<xi(t0)e−∫tt0αii(v)dv+(u+ε)1+σ2 forall t∈[Λ, +∞) and i∈Q. |
By taking the upper limits, (W1) leads to
0≤u≤(u+ε)1+σ2, |
which, together with the arbitrariness of ε, implies that u=0. This finishes the proof of Theorem 3.1.
We now make the following assumptions:
(W∗2) 1>supt≥t0[n∑j=1,j≠iαij(t)αii(t)+m∑j=1βij(t)αii(t)] forall i∈Q.
(W3) There exists a positive constant m[α]>0 such that
supt−gij(t)≥0∫tt−gij(t)α(s)ds≤m[α] forall i∈Q and j∈I. |
Theorem 3.2. Under assumptions (W1), (W2) and (W3), all positive solutions of system (1.3) are generalized exponentially convergent to the zero equilibrium point. More precisely, there exist positive constants κ,λ,M(φ) such that every solution x(t)=x(t;t0,φ) of (1.3) and (1.4) agrees with that
xi(t)≤κM(φ)e−λ∫tt0α(s)ds, ∀i∈Q. |
Proof. From (3.1) and (W3), one can see that
supt∈[T∗,+∞)[n∑j=1,j≠iαij(t)+m∑j=1βij(t)−αii(t)]<0 forall i∈Q, | (3.3) |
and there exists λ>0 such that
λαii(t)+n∑j=1,j≠iαij(t)+m∑j=1βij(t)eλ∫tt−gij(t)α(s)ds<αii(t), | (3.4) |
where t∈[T∗,+∞) and i∈Q.
Let us consider the following function:
v(t)=κM(φ)e−λ∫tt0α(s)ds, |
where t≥t0−ri, κ=eλ∫T∗t0α(s)ds, and M(φ) is defined in (3.2). In view of (3.4) and the fact that
v(t−gij(t))=eλ∫tt−gij(t)α(s)dsv(t), |
we obtain
v′(t)=κM(φ)e−λ∫tt0α(s)ds[−λα(t)]=[−λα(t)]v(t)≥[−λαii(t)]v(t)≥[n∑j=1,j≠iαij(t)+m∑j=1βij(t)eλ∫tt−gij(t)α(s)ds−αii(t)]v(t)=[n∑j=1,j≠iαij(t)+m∑j=1βij(t)v(t−gij(t))v(t)−αii(t)]v(t)=−αii(t)v(t)+n∑j=1,j≠iαij(t)v(t)+m∑j=1βij(t)v(t−gij(t)) | (3.5) |
where t∈[T∗,+∞) and i∈Q.
Hereafter, we show that
xi(t)<v(t)=κM(φ)e−λ∫tt0α(s)ds for all t≥t0−ri and i∈Q. | (3.6) |
It is easy to see that
v(t)=κM(φ)e−λ∫tt0α(s)ds=M(φ)eλ∫T∗t0α(s)dse−λ∫tt0α(s)ds=M(φ)eλ∫T∗tα(s)ds≥M(φ)>xi(t) for all t∈[t0−ri,T∗]. |
Next, we prove that
xi(t)<v(t), ∀t∈(T∗,+∞). | (3.7) |
Suppose, for the sake of contradiction, there exist ˜t>T∗ and i0∈Q satisfying that
xi0(˜t)=v(˜t),xj(t)<v(t) for all t∈[t0−rj, ˜t) and j∈Q. |
Then, xi0(t)≤v(t) for all t∈[t0−ri,˜t]. From (3.5), we have
v′(t)≥−αi0i0(t)v(t)+n∑j=1,j≠i0αi0j(t)v(t)+m∑j=1βi0j(t)v(t−gi0j(t)), |
which together with (1.3) and (3.3), we gain
0≤x′i0(˜t)−v′(˜t)=−αi0i0(˜t)xi0(˜t)+n∑j=1,j≠i0αi0j(˜t)xj(˜t)+m∑j=1βi0j(˜t)xi0(˜t−gi0j(˜t))e−γi0j(˜t)xi0(˜t−hi0j(˜t))−v′(˜t)≤−αi0i0(˜t)xi0(˜t)+m∑j=1βi0j(˜t)xi0(˜t−gi0j(˜t))+n∑j=1,j≠i0αi0j(˜t)xj(˜t)+αi0i0(˜t)v(˜t)−m∑j=1βi0j(˜t)v(˜t−gi0j(˜t))−n∑j=1,j≠i0αi0j(˜t)v(˜t)<−αi0i0(˜t)v(˜t)+m∑j=1βi0j(˜t)v(˜t−gi0j(˜t))+n∑j=1,j≠i0αi0j(˜t)v(˜t)+αi0i0(˜t)v(˜t)−m∑j=1βi0j(˜t)v(˜t−gi0j(˜t))−n∑j=1,j≠i0αi0j(˜t)v(˜t)=0, |
which yields a contradiction and shows (3.7) holds. Consequently, (3.6) holds, and the proof of Theorem 3.2 is now completed.
Theorem 3.3. Let (W1), (W∗2) and (W3) be satisfied, the zero equilibrium point in (1.3) is globally asymptotically stable.
Proof. According to Theorem 3.2, it follows from (W1), (W∗2) and (W3) that the zero equilibrium point in system (1.3) is globally attractive. Now, we only need to show the local stability on the zero equilibrium point in system (1.3). By (W∗2), we gain
−αii(t)+n∑j=1,j≠iαij(t)+m∑j=1βij(t)<0 forall t∈[t0, +∞) and i∈Q. | (3.8) |
For ε>0, we choose δ∈(0, ε). Next, we claim that, for ‖φ‖<δ,
xi(t)=xi(t;t0,φ)<ε forall t∈[t0,+∞) and i∈Q. | (3.9) |
In the contrary case, there exist t∗∈(t0,+∞) and i∗∈Q such that
xi∗(t∗)=ε, xj(t)<ε forall t∈[t0−rj, t∗) and j∈Q. | (3.10) |
Consequently, (1.3), (3.8) and (3.10) involve that
0≤x′i∗(t∗)≤[−αi∗i∗(t∗)+n∑j=1,j≠i∗αi∗j(t∗)+m∑j=1βi∗j(t∗)]ε<0, |
which is absurd and proves (3.9). Therefore, the zero equilibrium point is globally asymptotically stable.
Remark 3.1. One can easily see that all convergence results in references [23,24,26,32] are special ones of Theorems 3.1, 3.2 and 3.3 with n=1 and gij(t)≡hij(t). In addition, the errors in the existing literatures [39,40,41,43] have been corrected in the proof of Lemma 2.1 and Theorem 3.2. Specially, we add a new delay-dependent condition W3 to ensure the correctness of the main conclusions Theorems 3.1 and 3.2 in [39]. Therefore, the obtained results of this present paper improve and complement the above mentioned references [23,24,26,39,40,41,43].
In this section, we present an example to check the validity of the main results obtained in Section 3.
Example 4.1. Consider the following patch structure Nicholson's blowflies model with multiple pairs of time-varying delays:
{x′1(t)=−112t+3x1(t)+12t+3x2(t) +62t+3x1(t−|sint|)e−(1+cos2t)x1(t−12|sint|) +32t+3x1(t−23|sint|)e−(1+sin2t)x1(t−35|sint|), x′2(t)=−132t+2x2(t)+32t+2x1(t) +52t+2x2(t−34|cost|)e−(1+sint)x2(t−|cost|) +12t+2x2(t−16|cost|)e−(1+cost)x2(t−13|cost|). | (4.1) |
One can see that (4.1) satisfies all the conditions in Theorem 3.3. It follows that the zero equilibrium point of system (4.1) is globally asymptotically stable on C+={φ∈C([−1, 0],[0, +∞))×C([−1, 0],[0, +∞))}. The numerical trajectories with different initial values are shown in Figure 2, which reveal the above conclusions.
Remark 4.1. It should be pointed out that the per capita daily adult mortality terms of (4.1) are not uniformly positive mortality terms since
limt→+∞α11(t)=limt→+∞112t+3=limt→+∞α22(t)=limt→+∞132t+2=0. |
In addition, the global stability of the patch structure Nicholson's blowflies model with multiple pairs of time-varying delays has not been touched in [46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66]. This implies that all the results in [23,24,25,26,27,28,29,30,31,32,33] and [46,47,48,49,50,51,52,53,54,55,56,57,58,59,60,61,62,63,64,65,66] cannot be used to show the global convergence and stability on system (3.1) where τij(t)≢hij(t)(i∈Q,j∈I). Moreover, the proposed techniques could be taken into consideration in the dynamics research on other patch structure population models involving two or more delays in the same time-dependent birth function.
In the present manuscript, we investigate the asymptotic behavior for a patch structure Nicholson's blowflies model involving multiple pairs of different time-varying delays. Here, we obtain some novel results about the asymptotic behavior on the zero equilibrium point of the addressed model without assuming the uniform positiveness of the death rate and the boundedness of coefficients, which complement some earlier publications to some extent. In addition, the method used in this paper provides a possible method for studying the global asymptotic stability of other patch structure population dynamic models with multiple pairs of different time-varying delays.
I would like to thank the anonymous referees and the editor for very helpful suggestions and comments which led to improvements of my original paper. This work was supported by the Postgraduate Scientific Research Innovation Project of Hunan Province (No. CX20200892) and "Double first class" construction project of CSUST in 2020 ESI construction discipline Grant No. 23/03.
We confirm that we have no conflict of interest.
[1] | W. Li, L. Huang, J. Ji, Globally exponentially stable periodic solution in a general delayed predator-prey model under discontinuous prey control strategy, Discrete Contin. Dyn. Syst. Ser. B, 25 (2020), 2639-2664. |
[2] | T. Chen, L. Huang, P. Yu, et al. Bifurcation of limit cycles at infinity in piecewise polynomial systems, Nonlinear Anal. Real World Appl., 41 (2018), 82-106. |
[3] | C. Huang, S. Wen, M. Li, et al. An empirical evaluation of the influential nodes for stock market network: Chinese A shares case, Finance Research Letters, (2020), 101517. Available from: https://doi.org/10.1016/j.frl.2020.101517. |
[4] | F. Karim, S. Chauhan, J. Dhar, On the comparative analysis of linear and nonlinear business cycle model: Effect on system dynamics, economy and policy making in general, Quantit. Finance Econ., 4 (2020), 172-203. |
[5] | X. Yang, S. Wen, Z. Liu, et al. Dynamic properties of foreign exchange complex network, Mathematics, 7 (2019). Available from: https://doi.org/10.3390/math7090832. |
[6] | F. Wen, Y. Yuan, W. Zhou, Cross-shareholding networks and stock price synchronicity: Evidence from China, Int. J. Fin. Econ., (2020). Available from: https://doi.org/10.1002/ijfe.1828. |
[7] | C. Huang, Y. Tan, Global behavior of a reaction-diffusion model with time delay and Dirichlet condition, J. Diff. Equations, (2020). Available from: https://doi.org/10.1016/j.jde.2020.08.008. |
[8] | J. Cao, F. Wen, The impact of the Cross-Shareholding network on extreme price movements: Evidence from China, J. Risk, 22 (2019), 79-102. |
[9] | L. Huang, H. Ma, J. Wang, et al. Global dynamics of a Filippov plant disease model with an economic threshold of infected-susceptible ratio, J. Appl. Anal. Comput., 10 (2020), 1-15. |
[10] | Z. Ye, C. Hu, L. He, et al. The dynamic time-frequency relationship between international oil prices and investor sentiment in China: A wavelet coherence analysis, Energy J., 41 (2020). Available from: https://doi.org/10.5547/01956574.41.5.fwen. |
[11] | Y. Tan, Dynamics analysis of Mackey-Glass model with two variable delays, Math. Biosci. Eng., 17 (2020), 4513-4526. |
[12] | H. Hu, X. Yuan, L. Huang, et al. Global dynamics of an sirs model with demographics and transfer from infectious to susceptible on heterogeneous networks, Math. Biosci. Eng., 16 (2019), 5729-5749. |
[13] | C. Huang, X. Long, L. Huang, et al. Stability of almost periodic Nicholson's blowflies model involving patch structure and mortality terms, Canad. Math. Bull., 63 (2020), 405-422. |
[14] | Q. Cao, X. Guo, Anti-periodic dynamics on high-order inertial Hopfield neural networks involving time-varying delays, AIMS Mathematics, 5 (2020), 5402-5421. |
[15] | C. Huang, X. Zhao, J. Cao, et al. Global dynamics of neoclassical growth model with multiple pairs of variable delays, Nonlinearity, (2020). Available from: https: //doi.org/10.1088/1361-6544/abab4e. |
[16] | C. Huang, Y. Qiao, L. Huang, et al. Dynamical behaviors of a food-chain model with stage structure and time delays, Adv. Difference Equ., 2018 (2018), 1-26. |
[17] | C. Song, S. Fei, J. Cao, et al. Robust synchronization of fractional-order uncertain chaotic systems based on output feedback sliding mode control, Mathematics, 7 (2019). Available from: https://doi.org/10.3390/math7070599. |
[18] | L. Berezansky, E. Braverman, A note on stability of Mackey-Glass equations with two delays, J. Math. Anal. Appl., 450 (2017), 1208-1228. |
[19] | C. Huang, X. Long, J. Cao, Stability of anti-periodic recurrent neural networks with multiproportional delays, Math. Methods Appl. Sci., 43 (2020), 6093-6102. |
[20] | C. Huang, H. Yang, J. Cao, Weighted pseudo almost periodicity of multi-proportional delayed shunting inhibitory cellular neural networks with D operator, Discrete Contin. Dyn. Syst. Ser. S, (2020). Available from: https://doi.org/10.3934/dcdss.2020372. |
[21] | W. Tang, J. Zhang, Symmetric integrators based on continuous-stage Runge-Kutta-Nystrom methods for reversible systems, Appl. Math. Comput., 361 (2019), 1-12. |
[22] | V. Volttera, Variazioni e fluttuazioni del numerou d'individui in specie animali conviventi, R. Comitato Talassografico Italiano, Memoria, 131 (1927), 1-142. |
[23] | C. Huang, L. Yang, J. Cao, Asymptotic behavior for a class of population dynamics, AIMS Mathematics, 5 (2020), 3378-3390. |
[24] | X. Long, S. Gong, New results on stability of Nicholson's blowflies equation with multiple pairs of time-varying delays, Appl. Math. Lett., 100 (2020), 106027. |
[25] | L. Berezansky, E. Braverman, Nicholson's blowflies differential equations revisited: Main results and open problems, Appl. Math. Model., 34 (2010), 1405-1417. |
[26] | V. H. Le, Global asymptotic behaviour of positive solutions to a non-autonomous Nicholson's blowflies model with delays, J. Biol. Dyn., 8 (2014), 135-144. |
[27] | B. Liu, Global exponential stability of positive periodic solutions for a delayed Nicholson's blowflies model, J. Math. Anal. Appl., 412 (2014), 212-221. |
[28] | F. Liu, L. Feng, V. Anh, et al. Unstructured-mesh Galerkin finite element method for the twodimensional multi-term time-space fractional Bloch-Torrey equations on irregular convex domains, Comput. Math. Appl., 78 (2019), 1637-1650. |
[29] | X. Chen, C. Shi, Permanence of a Nicholson's blowflies model with feedback control and multiple time-varying delays, Chinese Quart. J. Math., 1 (2015), 153-158. |
[30] | B. Liu, New results on global exponential stability of almost periodic solutions for a delayed Nicholson blowflies model, Ann. Polon. Math., 113 (2015), 191-208. |
[31] | S. Zhou, Y. Jiang, Finite volume methods for N-dimensional time fractional Fokker-Planck equations, Bull. Malays. Math. Sci. Soc., 42 (2019), 3167-3186. |
[32] | T. S. Doan, V. H. Le, T. A. Trinh, Global attractivity of positive periodic solution of a delayed Nicholson model with nonlinear density-dependent mortality term, Electron. J. Qual. Theory Differ. Equ., 8 (2019), 1-21. |
[33] | Z. Long, Exponential convergence of a non-autonomous Nicholson's blowflies model with an oscillating death rate, Electron. J. Qual. Theory Differ. Equ., 2016 (2016), 1-7. |
[34] | L. Berezansky, E. Braverman, On exponential stability of a linear delay differential equation with an oscillating coefficient, Appl. Math. Lett., 22 (2009), 1833-1837. |
[35] | J. Shao, Pseudo almost periodic solutions for a Lasota-Wazewska model with an oscillating death rate, Appl. Math. Lett., 43 (2015), 90-95. |
[36] | Q. Cao, X. Long, New convergence on inertial neural networks with time-varying delays and continuously distributed delays, AIMS Mathematics, 5 (2020), 5955-5968. |
[37] | B. Liu, Global stability of a class of Nicholson's blowflies model with patch structure and multiple time-varying delays, Nonlinear Anal. Real World Appl., 11 (2010), 2557-2562. |
[38] | G. Yang, Dynamical behaviors on a delay differential neoclassical growth model with patch structure, Math. Methods Appl. Sci., 41 (2018), 3856-3867. |
[39] | Y. Xu, Q. Cao, X. Guo, Stability on a patch structure Nicholson's blowflies system involving distinctive delays, Appl. Math. Lett., 105 (2020). Available from: https://doi.org/10.1016/j.aml.2020.106340. |
[40] | Y. Xu, Q. Cao, Global asymptotic stability for a nonlinear density-dependent mortality Nicholson's blowflies system involving multiple pairs of time-varying delays, Adv. Difference Equ., 2020 (2020), 1-14. |
[41] | H. Zhang, Q. Cao, H. Yang, Asymptotically almost periodic dynamics on delayed Nicholson-type system involving patch structure, J. Inequal. Appl., 2020 (2020). Available from: https://doi.org/10.1186/s13660-020-02366-0. |
[42] | C. Qian, Y. Hu, Novel stability criteria on nonlinear density-dependent mortality Nicholson's blowflies systems in asymptotically almost periodic environments, J. Inequal. Appl., 2020 (2020), 1-18. |
[43] | Q. Cao, G. Wang, H. Zhang, et al. New results on global asymptotic stability for a nonlinear density-dependent mortality Nicholson's blowflies model with multiple pairs of time-varying delays, J. Inequal. Appl., 2020 (2020), 1-12. |
[44] | H. L. Smith, Monotone Dynamical Systems: An Introduction to the Theory of Competitive and Cooperative Systems, Mathematical Surveys and Monographs, 1995. |
[45] | J. K. Hale, S. M. Verduyn Lunel, Introduction to Functional Differential Equations, SpringerVerlag, New York, 1993. |
[46] | W. Wang, W. Chen, Stochastic Nicholson-type delay system with regime switching, Systems Control Lett., 136 (2020), 104603. |
[47] | W. Wang, Positive periodic solutions of delayed Nicholson's blowflies models with a nonlinear density-dependent mortality term, Appl. Math. Model., 36 (2012), 4708-4713. |
[48] | C. Wang, R. P. Agarwal, S. Rathinasamy, Almost periodic oscillations for delay impulsive stochastic Nicholson's blowflies timescale model, Comput. Appl. Math., 37 (2018), 3005-3026. |
[49] | C. Wang, Existence and exponential stability of piecewise mean-square almost periodic solutions for impulsive stochastic Nicholson's blowflies model on time scales, Appl. Math. Comput., 248 (2014), 101-112. |
[50] | Z. Cai, J. Huang, L. Huang, Periodic orbit analysis for the delayed Filippov system, Proc. Amer. Math. Soc., 146 (2018), 4667-4682. |
[51] | C. Huang, H. Zhang, L. Huang, Almost periodicity analysis for a delayed Nicholson's blowflies model with nonlinear density-dependent mortality term, Pure Appl. Anal., 18 (2019), 3337-3349. |
[52] | H. Hu, T. Yi, X. Zou, On spatial-temporal dynamics of a Fisher-KPP equation with a shifting environment, Proc. Amer. Math. Soc., 148 (2020), 213-221. |
[53] | H. Hu, X. Zou, Existence of an extinction wave in the Fisher equation with a shifting habitat, Proc. Amer. Math. Soc., 145 (2017), 4763-4771. |
[54] | J. Wang, X. Chen, L. Huang, The number and stability of limit cycles for planar piecewise linear systems of node-saddle type, J. Math. Anal. Appl., 469 (2019), 405-427. |
[55] | J. Wang, C. Huang, L. Huang, Discontinuity-induced limit cycles in a general planar piecewise linear system of saddle-focus type, Nonlinear Anal. Hybrid Syst., 33 (2019), 162-178. |
[56] | R. Wei, J. Cao, C. Huang, Lagrange exponential stability of quaternion-valued memristive neural networks with time delays, Math. Methods Appl. Sci., 43 (2020), 7269-7291. |
[57] | C. Huang, H. Zhang, J. Cao, et al. Stability and hopf bifurcation of a delayed prey-predator model with disease in the predator, Internat. J. Bifur. Chaos Appl. Sci. Engrg., 29 (2019), 1950091. |
[58] | C. Huang, J. Wang, L. Huang, New results on asymptotically almost periodicity of delayed Nicholson-type system involving patch structure, Electron. J. Differential Equations, 2020 (2020), 1-17. |
[59] | C. Huang, J. Cao, F. Wen, et al. Stability analysis of SIR model with distributed delay on complex networks, Plos One, 11 (2016), 1-22. |
[60] | C. Huang, L. Liu, Boundedness of multilinear singular integral operator with non-smooth kernels and mean oscillation, Quaest. Math., 40 (2017), 295-312. |
[61] | X. Zhang, H. Hu, Convergence in a system of critical neutral functional differential equations, App. Math. Lett., 107 (2020). Available from: https://doi.org/10.1016/j.aml.2020.106385. |
[62] | C. Huang, S. Wen, L. Huang, Dynamics of anti-periodic solutions on shunting inhibitory cellular neural networks with multi proportional delays, Neurocomputing, 357 (2019), 47-52. |
[63] | H. Hu, L. Liu, Weighted inequalities for a general commutator associated to a singular integral operator satisfying a variant of Hormander's condition, Math. Notes, 101 (2017), 830-840. |
[64] | J. Wang, S. He, L. Huang, Limit cycles induced by threshold nonlinearity in planar piecewise linear systems of node-focus or node-center type, Internat. J. Bifur. Chaos Appl. Sci. Engrg., (2020). Available from: https://doi.org/10.1142/S0218127420501606. |
[65] | C. Huang, Z. Yang, T. Yi, et al. On the basins of attraction for a class of delay differential equations with non-monotone bistable nonlinearities, J. Differential Equations, 256 (2014), 2101-2114. |
[66] | J. Zhang, C. Huang, Dynamics analysis on a class of delayed neural networks involving inertial terms, Adv. Difference Equ., 2020 (2020), 1-12. |
1. | Changjin Xu, Wei Zhang, Peiluan Li, Chuangxia Huang, Periodic Oscillating Dynamics for a Delayed Nicholson-Type Model with Harvesting Terms, 2021, 2021, 1563-5147, 1, 10.1155/2021/9734342 | |
2. | Xianhui Zhang, Convergence analysis of a patch structure Nicholson’s blowflies system involving an oscillating death rate, 2021, 0952-813X, 1, 10.1080/0952813X.2021.1908433 | |
3. | Jian Zhang, Ancheng Chang, Gang Yang, Periodicity on Neutral-Type Inertial Neural Networks Incorporating Multiple Delays, 2021, 13, 2073-8994, 2231, 10.3390/sym13112231 | |
4. | Qian Cao, Positive almost periodicity on a n-species food chain system incorporating time-varying delays and polluted environments, 2022, 0952-813X, 1, 10.1080/0952813X.2022.2080279 | |
5. | Xiaojin Guo, Chuangxia Huang, Jinde Cao, Nonnegative periodicity on high-order proportional delayed cellular neural networks involving $ D $ operator, 2020, 6, 2473-6988, 2228, 10.3934/math.2021135 | |
6. | Wenke Wang, Le Li, Xuejun Yi, Chuangxia Huang, Convergence on Population Dynamics and High-Dimensional Haddock Conjecture, 2021, 13, 2073-8994, 2252, 10.3390/sym13122252 | |
7. | Xiaodan Ding, Global attractivity of Nicholson’s blowflies system with patch structure and multiple pairs of distinct time-varying delays, 2023, 16, 1793-5245, 10.1142/S1793524522500814 | |
8. | Yang Gao, Global stability for a new predator–prey model with cross-dispersal among patches based on graph theory, 2021, 2021, 1687-1847, 10.1186/s13662-021-03645-w | |
9. | Chuangxia Huang, Xiaojin Guo, Jinde Cao, Ardak Kashkynbayev, Bistable dynamics on a tick population equation incorporating Allee effect and two different time-varying delays, 2022, 15, 1937-1632, 3173, 10.3934/dcdss.2022122 | |
10. | Hong Zhang, Qian Cao, Hedi Yang, Dynamics analysis of delayed Nicholson-type systems involving patch structure and asymptotically almost periodic environments, 2022, 34, 0952-813X, 725, 10.1080/0952813X.2021.1924869 | |
11. | Xiaoyuan Chang, Junping Shi, Bistable and oscillatory dynamics of Nicholson's blowflies equation with Allee effect, 2022, 27, 1531-3492, 4551, 10.3934/dcdsb.2021242 | |
12. | Yanli Xu, Stability on a patch-structured Nicholson's blowflies system incorporating mature delay and feedback delay, 2022, 0020-7179, 1, 10.1080/00207179.2022.2081260 | |
13. | Chuangxia Huang, Lihong Huang, Jianhong Wu, Global population dynamics of a single species structured with distinctive time-varying maturation and self-limitation delays, 2022, 27, 1531-3492, 2427, 10.3934/dcdsb.2021138 | |
14. | Rajendiran Eswari, Jehad Alzabut, Mohammad Esmael Samei, Hui Zhou, On periodic solutions of a discrete Nicholson’s dual system with density-dependent mortality and harvesting terms, 2021, 2021, 1687-1847, 10.1186/s13662-021-03521-7 | |
15. | Yang Liu, Zhen Wang, Qian Ma, Hao Shen, Multistability analysis of delayed recurrent neural networks with a class of piecewise nonlinear activation functions, 2022, 152, 08936080, 80, 10.1016/j.neunet.2022.04.015 | |
16. | Xian Zhao, Chuangxia Huang, Bingwen Liu, Jinde Cao, Stability analysis of delay patch-constructed Nicholson’s blowflies system, 2024, 222, 03784754, 379, 10.1016/j.matcom.2023.09.012 | |
17. | Qian Wang, Lihong Huang, Almost Periodic Dynamics of a Delayed Patch-Constructed Nicholson’s Blowflies System, 2024, 23, 1575-5460, 10.1007/s12346-024-01129-2 |