Next Article in Journal
Nuclear Structure Evolution Reflected from Local Relations
Next Article in Special Issue
On the Crossing Numbers of the Join Products of Six Graphs of Order Six with Paths and Cycles
Previous Article in Journal
Optical Fiber Sensors for Monitoring Railway Infrastructures: A Review towards Smart Concept
Previous Article in Special Issue
Periodicity on Neutral-Type Inertial Neural Networks Incorporating Multiple Delays
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Convergence on Population Dynamics and High-Dimensional Haddock Conjecture

1
College of Meteorology and Oceanography, National University of Defense Technology, Changsha 410073, China
2
School of Mathematics, Hunan University, Changsha 410082, China
3
Hunan Provincial Key Laboratory of Mathematical Modeling and Analysis in Engineering, Changsha University of Science and Technology, Changsha 410114, China
*
Author to whom correspondence should be addressed.
Symmetry 2021, 13(12), 2252; https://doi.org/10.3390/sym13122252
Submission received: 10 October 2021 / Revised: 18 November 2021 / Accepted: 18 November 2021 / Published: 26 November 2021

Abstract

:
One fundamental step towards grasping the global dynamic structure of a population system involves characterizing the convergence behavior (specifically, how to characterize the convergence behavior). This paper focuses on the neutral functional differential equations arising from population dynamics. With the help of monotonicity techniques and functional methods, we analyze the subtle relations of both the ω -limited set and special point. Meanwhile, we prove that every bounded solution converges to a constant vector, as t tends to positive infinity. Our results correlate with the findings from earlier publications, and our proof yields an improved Haddock conjecture.

1. Introduction

It is well known that a delay differential equation model plays an important role in the process of population dynamics modeling, especially when the species clearly goes through the incubation period and evolves at different stages [1,2,3,4]. In addition to stability, limited cycles, bifurcation, and periodic solutions—how to characterize convergence behavior is a fundamental and important step toward grasping the global dynamic structure of a population system [5,6,7,8,9,10,11,12]. As early as 1976, at an international conference, Bernfeld and Haddock noted that every solution of the scalar population delay differential equation converges to a constant [13]
[ x ( t ) p x ( t r ) ] = x 1 3 ( t ) + x 1 3 ( t r ) , r > 0 .
Consequently, the above scalar neutral functional differential equation (NFDE) has been naturally generalized to the following autonomous NFDE [14]
[ x ( t ) p x ( t r ) ] = F ( x ( t ) ) + G ( x ( t r ) )
and non-autonomous NFDE [15]
[ x ( t ) p x ( t σ ( t ) ) ] = b ( t ) [ F ( x ( t ) ) + F ( x ( t σ ( t ) ) ) ] ,
respectively. Here, F C ( R 1 , R 1 ) is strictly increasing, G C ( R 1 , R 1 ) , and σ , b C ( R 1 , R 1 ) are bounded.
In fact, several important infinite dimensional dynamical systems arising from population growth are special cases of the following model [5,16,17]
( x 1 ( t ) c 2 x 2 ( t r 2 ) ) = F ( x 1 ( t ) ) + G ( x 2 ( t r 2 ) ) , ( x 2 ( t ) c 3 x 3 ( t r 3 ) ) = F ( x 2 ( t ) ) + G ( x 3 ( t r 3 ) ) , , ( x n 1 ( t ) c n x n ( t r n ) ) = F ( x n 1 ( t ) ) + G ( x n ( t r n ) ) , ( x n ( t ) c 1 x 1 ( t r 1 ) ) = F ( x n ( t ) ) + G ( x 1 ( t r 1 ) ) ,
where F describes the instantaneous mortality rate, and G C ( R 1 , R 1 ) is the feedback function incorporating with positive delays r 1 , , r n . Clearly, system (1) can be considered as a neutral delayed cellular neural network accompanying D operators [18]. When F is a strictly monotonically increasing function, c i ( 0 , 1 ) , r i > 0 for i Ω = { 1 , 2 , · · · , n } , and either G ( x ) F ( x ) or G ( x ) F ( x ) for all x, systems (1) has been used to confirm the high-dimensional Haddock conjecture [13,15,16].
Although the scalar Haddock conjecture has been successfully demonstrated by the authors in the literature [14,15], the high-dimensional structure brings practical difficulties to prove the high-dimensional Haddock conjecture based on the system (1). This makes the existing proofs of the high-dimensional Haddock conjecture to be done only with with additional conditions. In particular, when n = 2 , for the system (1), the authors in [16] proved that each solution converges to a constant vector, as long as the solution is bounded. Furthermore, the authors in [17] noted that it is necessary to add one of the following assumptions to ensure the correctness of the main results in [13].
  • (S + ) Assume G F , and if I is an any given bounded interval on R 1 , then there is a constant μ > 0 such that F ( x 2 ) F ( x 1 ) μ ( x 2 x 1 ) f o r a n y [ x 1 , x 2 ] I .
  • (S ) Assume G F , and if I is an any given bounded interval on R 1 , then there is a constant μ > 0 such that F ( x 2 ) F ( x 1 ) μ ( x 2 x 1 ) f o r a n y [ x 1 , x 2 ] I .
Later on, under the technical assumptions, two-dimensional cases of systems (1) are investigated in [19,20,21,22] as well. Regrettably, if one chooses G ( x ) F ( x ) = x 1 3 designed in the Haddock conjecture, simple verification can guide that the function F and G do not meet the requirements of conditions (S + ) or (S ). Therefore, even involving the assumptions (S + ) and (S ), for high-dimensional neutral functional differential equations (1), the Haddock conjecture is still untenable. This motivates us to remove the assumptions (S + ) and (S ), and generalize the Haddock conjecture to a high-dimensional version (1) to reveal that every bounded solution converges to a constant.
Without adopting the above-mentioned assumptions (S + ) and (S ), in the following section, we develop the delicate monotonicity with systems (1) and analyze the subtle relation of both the ω -limited set and special point with the help of monotonicity techniques and functional methods. In Section 3, the main results are presented and demonstrated by applying the preliminary lemmas of Section 2. Moreover, our results in this paper correlate with earlier publications, and Corollary 1 yields an improved Haddock conjecture.

2. Preliminaries

For convenience, we present the following notations. We denote
C = i = 1 n C ( [ r i , 0 ] , R 1 ) , | r | = i = 1 n r i , Ω = { 1 , 2 , · · · , n } , τ = min j Ω r j , j ¯ = j for j Ω , 0 ¯ = n , n + 1 ¯ = 1 .
Equipping with the usual supremum norm, it is easy to see that C is a Banach space. We also define the order cones K in the following way:
K = { φ C + : φ i ( 0 ) c i + 1 ¯ φ i + 1 ¯ ( r i + 1 ¯ ) 0 f o r i Ω } ,
where
C + = i = 1 n C ( [ r i , 0 ] , R + 1 ) .
For A C , and ϕ 1 , ϕ 2 C , we designate
  • ϕ 1 ϕ 2 iff ϕ 2 ϕ 1 Int C + ; ϕ 1 K ϕ 2 iff ϕ 2 ϕ 1 Int K ; ϕ 1 K A iff ϕ 1 K ϕ 2 for all ϕ 2 A .
  • ϕ 1 ϕ 2 iff ϕ 2 ϕ 1 C + ; ϕ 1 K ϕ 2 iff ϕ 2 ϕ 1 K ; ϕ 1 K A iff ϕ 1 K ϕ 2 for all ϕ 2 A .
  • ϕ 1 < ϕ 2 iff ϕ 1 ϕ 2 and ϕ 1 ϕ 2 ; ϕ 1 < K ϕ 2 iff ϕ 1 K ϕ 2 and ϕ 1 ϕ 2 ; ϕ 1 < K A iff ϕ 1 < K ϕ 2 for all ϕ 2 A .
It is worth noting that ϕ 1 K ϕ 2 and ϕ 1 K ϕ 2 can be defined naturally. For i Ω , σ > 0 , x i C ( [ r i , σ ] , R 1 ) , it is obvious that x t C can be presented as x t = ( x 1 , t , x 2 , t , , x n , t ) and x i , t ( ξ ) = x i ( t + ξ ) for   all   ξ [ r i , 0 ] , t [ 0 , σ ] . Furthermore, we also define
β ^ = ( ( β ^ ) 1 , ( β ^ ) 2 , , ( β ^ ) n ) ,
in which ( β ^ ) i ( ξ ) = β , i Ω , ξ [ r i , 0 ] . For the sake of simplicity, if φ C is an initial value, we use the notation x t ( φ ) ( x ( t , φ ) ) to represent the solution of the initial value problem of systems (1).
Remark 1.
Assume there are α R 1 and φ C , such that φ K α ^ , then it is easy to establish two statements as below:
  • (i) If φ i ( 0 ) c i + 1 ¯ φ i + 1 ¯ ( r i + 1 ¯ ) > ( 1 c i + 1 ¯ ) α f o r   s o m e   i Ω , then φ i ( 0 ) > α .
  • (ii) If φ i ( 0 ) c i + 1 ¯ φ i + 1 ¯ ( r i + 1 ¯ ) = ( 1 c i + 1 ¯ ) α f o r   s o m e   i Ω , then either φ i ( 0 ) = φ i + 1 ¯ ( r i + 1 ¯ ) = α or φ i ( 0 ) < φ i + 1 ¯ ( r i + 1 ¯ ) .
Remark 2.
Assume that σ > 0 , x i C ( [ r i , σ ] , R 1 ) f o r i Ω and x 0 K α ^ , it is easy to obtain the following statements:
  • (i) For all t [ 0 , σ ] and i Ω , if x i ( t ) c i + 1 ¯ x i + 1 ¯ ( t r i + 1 ¯ ) ( 1 c i + 1 ¯ ) α , then, for all t [ 0 , σ ] , we have x t K α ^ .
  • (ii) If σ max i Ω r i , and x t K α ^ holds for any t [ 0 , σ ] , furthermore, x i ( t ) c i + 1 ¯ x i + 1 ¯ ( t r i + 1 ¯ ) > ( 1 c i + 1 ¯ ) α holds for all t [ σ r i , σ ] and i Ω , then we have x σ K α ^ .
In this section, for all x R 1 , we always assume that the inequality G ( x ) F ( x ) holds. Hereafter, we denote j = n + 1 n r j = 0 , j = n + 2 n r j = r 1 , and suppose that F is a strictly increasing function defined on R 1 .
Lemma 1.
The solution of the initial value problem of the system (1) exists and is unique on R + 1 .
Proof. 
To begin with, we show that x 1 ( t , φ ) exists and is unique on [ 0 , τ ] for any given initial value function φ C . For the convenience of expression, we denote a ( t ) = φ 2 ¯ ( t r 2 ¯ ) for t [ 0 , τ ] . Obviously, according to Lemma 2.2 in [16], one can obtain that the solution of the initial value problem of the below system exists and is unique on [ 0 , τ ] ,
y ( t ) = F ( y ( t ) + c 2 ¯ a ( t ) ) + G ( a ( t ) ) , y ( 0 ) = φ 1 ( 0 ) c 2 ¯ φ 2 ¯ ( r 2 ¯ ) .
Denote y ˜ ( t ) = x 1 ( t , φ ) c 2 ¯ φ 2 ¯ ( t r 2 ¯ ) , then y ˜ ( t ) satisfies (2); therefore, x 1 ( t , φ ) exists and is unique on [ 0 , τ ] . Using a similar method, one can also prove that x i ( t , φ ) exists and is unique on [ 0 , τ ] for i = 2 , 3 , , n . Therefore, the mathematical induction can be directly applied to know that Lemma 1 is indeed true.  □
As usual, we denote O ( φ ) = { x t ( φ ) φ C , t 0 } , O ( φ ) ¯ is the closure of O ( φ ) , and ω ( φ ) = t 0 O ( x t ( φ ) ) ¯ . If O ( φ ) is bounded, then O ( φ ) ¯ is compact in C, and ω ( x ) is nonempty, compact, invariant, and connected.
Next, we give the following lemmas, describing the monotonicity of the system (1).
Lemma 2.
If φ K α ^ holds for some φ C and α R 1 , then x t ( φ ) K α ^ holds for all t R + 1 .
Proof. 
For all t R + 1 and i Ω , we designate y i ( t ) = x i ( t , φ ) c i + 1 ¯ x i + 1 ¯ ( t r i + 1 ¯ , φ ) . Then for all t R + 1 and i Ω , we assert that y i ( t ) ( 1 c i + 1 ¯ ) α . In fact, assuming the opposite is true, one can choose i 0 Ω and t 0 ( 0 , τ ] , such that
y i 0 ( t 0 ) < ( 1 c i 0 + 1 ¯ ) α , a n d y i 0 ( t 0 ) < 0 .
Therefore, we have
x i 0 ( t 0 , φ ) < ( 1 c i 0 + 1 ¯ ) α + c i 0 + 1 ¯ x i 0 + 1 ¯ ( t 0 r i 0 + 1 ¯ , φ ) x i 0 + 1 ¯ ( t 0 r i 0 + 1 ¯ , φ ) .
By systems (1) and the monotonicity of function F, we have the following inequality
y i 0 ( t 0 ) = F ( x i 0 ( t 0 , φ ) ) + G ( x i 0 + 1 ¯ ( t 0 r i 0 + 1 ¯ , φ ) ) F ( x i 0 ( t 0 , φ ) ) + F ( x i 0 + 1 ¯ ( t 0 r i 0 + 1 ¯ , φ ) ) > 0 ,
which contradicts the assumption y i 0 ( t 0 ) < 0 , and proves the assertion.
Applying the assertion and Remark 2 (i), we have x t ( φ ) K α ^ for all t [ 0 , τ ] . Then Lemma 2 follows by the method of mathematical induction.  □
Lemma 3.
Under the hypotheses of Lemma 2, if φ i ( 0 ) c i + 1 ¯ φ i + 1 ¯ ( r i + 1 ¯ ) > ( 1 c i + 1 ¯ ) α for some i Ω , then x i 1 ¯ ( r i , φ ) c i x i ( 0 , φ ) > ( 1 c i ) α .
Proof. 
Otherwise, x i 1 ¯ ( r i , φ ) c i x i ( 0 , φ ) = ( 1 c i ) α . Applying Lemma 2, we get ( x i 1 ¯ ( t , φ ) c i x i ( t r i , φ ) ) | t = r i = 0 .
It follows from (1) that F ( x i 1 ¯ ( r i , φ ) ) = G ( x i ( 0 , φ ) ) F ( x i ( 0 , φ ) ) and hence x i 1 ¯ ( r i , φ ) x i ( 0 , φ ) = φ i ( 0 ) > α .
However, by Remark 1, we then obtain x i 1 ¯ ( r i , φ ) < x i ( 0 , φ ) , which yields a contradiction.  □
Lemma 4.
Assume O ( φ ) is bounded for some φ C . Define α * = sup { α R 1 , α ^ K ω ( φ ) } . Then we have the below statements:
  • (i) For i Ω , if ψ ω ( φ ) , then one can choose s * [ 0 , | r | ] , such that
    x i ( s * , ψ ) c i + 1 ¯ x i + 1 ¯ ( s * r i + 1 ¯ , ψ ) = ( 1 c i + 1 ¯ ) α * .
  • (ii) For i Ω , if ψ ω ( φ ) , then one can choose s * [ 0 , | r | ] such that
    x i ( s * + k | r | + j = i + 1 n r j , ψ ) c i + 1 ¯ x i + 1 ¯ ( s * + k | r | + j = i + 2 n r j , ψ ) = ( 1 c i + 1 ¯ ) α * ,
    where k 0 is an integer.
Proof. 
(i) By the symmetry of (1), we only need to consider the case i = n . Suppose the result is not true, one can choose ψ ω ( φ ) such that
x n ( t , ψ ) c 1 x 1 ( t r 1 , ψ ) > ( 1 c 1 ) α * , t [ 0 , | r | ] .
Iterating Lemma 3 n 1 times, we therefore have
x i ( t , ψ ) c i + 1 ¯ x i + 1 ¯ ( t r i + 1 ¯ , ψ ) > ( 1 c i + 1 ¯ ) α * ,
where t [ j = i + 1 n r j , | r | ] , i { 1 , 2 , , n 1 } .
It follows from Remark 2 that we have x | r | ( ψ ) K α * ^ . Applying Lemma 2 and combining the concept of ω ( φ ) , we can choose β * > α * such that ω ( φ ) K β * ^ K α * ^ . Note that this conflicts with the definition of α * . This completes the proof of the statement (i).
(ii) For all i Ω and k 0 , we denote
A i k = { s [ 0 , | r | ] : x i ( s + k | r | + j = i + 1 n r j , ψ ) c i + 1 ¯ x i + 1 ¯ ( s + k | r | + j = i + 2 n r j , ψ ) = ( 1 c i + 1 ¯ ) α * } .
Then A i k ϕ follows from (i).
For i = 1 , 2 , , n 1 , k 0 , applying Lemma 3, we obtain A i k A i + 1 k and A n k + 1 A 1 k . Consequently, k = 0 i = 1 n A i k ϕ . We choose s * k = 0 i = 1 n A i k as required.  □

3. Main Results

Based on the relevant preparations in the previous section, we can now present and demonstrate the main results of this article in this section.
Theorem 1.
If O ( φ ) is bounded for some φ C , and G ( x ) F ( x ) for all x R 1 , then one can choose α * R 1 , such that ω ( φ ) = { α * ^ } .
Proof. 
Define α * = sup { α R 1 : α ^ K ω ( φ ) } . We will prove that ω ( φ ) = { α * ^ } ; otherwise, we have ω ( φ ) \ { α * ^ } ϕ .
By combining Lemma 4 with the invariance of ω ( φ ) , one can choose ψ ω ( φ ) \ { α * ^ } , such that
x i ( k | r | + j = i + 1 n r i , ψ ) c i + 1 ¯ x i + 1 ¯ ( k | r | + j = i + 2 n r j , ψ ) = ( 1 c i + 1 ¯ ) α *
for i Ω and k 0 .
Integrating (1), we have
k | r | + s + j = i + 1 n r j ( k + 1 ) | r | + j = i + 1 n r j ( x i ( t , ψ ) c i + 1 ¯ x i + 1 ¯ ( t r i + 1 ¯ , ψ ) ) d t = k | r | + s + j = i + 1 n r j ( k + 1 ) | r | + j = i + 1 n r j F ( x i ( t , ψ ) ) + G ( x i + 1 ¯ ( t r i + 1 ¯ , ψ ) ) d t .
Let
H i k ( s ) = k | r | + s + j = i + 1 n r j ( k + 1 ) | r | + j = i + 1 n r j G ( x i ( t , ψ ) ) d t , h i k ( s ) = k | r | + s + j = i + 1 n r j ( k + 1 ) | r | + j = i + 1 n r j F ( x i ( t , ψ ) ) d t ,
where s [ 0 , | r | ] , k 0 and i Ω . It follows from (3) and (4), we have
( 1 c i + 1 ) α ( x i ( k | r | + s + j = i + 1 n r j , ψ ) c i + 1 x i + 1 ( k | r | + s + j = i + 1 n r j , ψ ) ) = h i k ( s ) + H i + 1 k ( s ) ,
with k 0 , i = 1 , 2 , · · · , n 1 , and
( 1 c 1 ) α * ( x n ( k | r | + s , ψ ) c 1 x 1 ( k | r | r 1 , ψ ) ) = h n k ( s ) + H 1 k 1 ( s ) , k 1 .
Let ξ i k ( s ) = x i ( k | r | + s + j = i + 1 n r j , ψ ) α * for all i Ω , k 0 and s [ 0 , | r | ] . Then for k 1 and s [ 0 , | r | ] , we have
ξ 1 k ( s ) + c 2 ξ 2 k ( s ) h 1 k ( s ) + h 2 k ( s ) , ξ 2 k ( s ) + c 3 ξ 3 k ( s ) h 2 k ( s ) + h 3 k ( s ) , , ξ n 1 k ( s ) + c n ξ n k ( s ) h n 1 k ( s ) + h n k ( s ) , ξ n k ( s ) + c 1 ξ 1 k 1 ( s ) h n k ( s ) + h 1 k 1 ( s ) .
Applying (5), we obtain
ξ 1 k ( s ) c 1 ξ 1 k 1 ( s ) + i = 2 n ( 1 c i ) ξ i k ( s ) h 1 k ( s ) h 1 k 1 ( s )
Consequently, for any positive integers N, by (6), we have
k = 1 N i = 1 n ( 1 c i ) ξ i k ( s ) h 1 N ( s ) h 1 0 ( s ) + c 1 ξ 1 0 ( s ) c 1 ξ 1 N ( s ) .
Making full use of the compactness property of ω ( φ ) , one can choose M > 0 , such that
| h i k ( s ) | M and | ξ i k ( s ) | M .
Hence, k = 1 N ( 1 c i ) ξ i k ( s ) 4 M holds for N 1 , i Ω and s [ 0 , | r | ] . Thus, it is easy to see that for any i Ω , ξ i k ( · ) uniformly convergent to 0 in [ 0 , | r | ] as k . Therefore,
lim k h i k ( 0 ) = F ( α * ) | r | .
On the other hand, applying (5) with s = 0 , we obtain h i + 1 k ( 0 ) h i k ( 0 ) for all i = 1 , 2 , · · · , n 1 and h 1 k 1 ( 0 ) h n k ( 0 ) , so that h i k ( 0 ) h i 0 ( 0 ) for all i Ω and k 0 . Since ψ ω ( φ ) \ { α * ^ } and F is strictly increasing, there exists i 0 { 1 , 2 , · · · , n } , such that h i 0 0 ( 0 ) > F ( α * ) | r | .
Thus h i 0 k ( 0 ) h i 0 0 ( 0 ) > F ( α * ) | r | , which contradicts with lim k h i 0 k = F ( α * ) | r | . Therefore, ω ( φ ) = { α * ^ } . This finishes the proof of Theorem 1.  □
Theorem 2.
If O ( φ ) is bounded for some φ C , and G ( x ) F ( x ) for all x R 1 , then one can choose α * R 1 , such that ω ( φ ) = { α * ^ } .
Proof. 
To begin with, for all x R 1 , we define f ( x ) , y ( t ) , g ( x ) as below
f ( x ) = F ( x ) , y ( t ) = x ( t , φ ) , and g ( x ) = G ( x ) .
Then it is obvious that y ( t ) satisfies system (1) by replacing F , G with f , g , respectively. Consequently, directly using the method similar to Theorem 1, we can easily prove that Theorem 2 is true.  □
Remark 3.
For all x R 1 , if G ( x ) F ( x ) and φ K α ^ holds for some φ C and α R 1 , then directly employing the method similar to Lemma 2, one can easily establish the statement that x t ( φ ) K α ^ for all t R + 1 .
Thus, when G ( x ) F ( x ) , applying Lemma 2 and the above statement, it is easy to reveal that every solution of (1) is bounded. Furthermore, based on the above two theorems, we can directly establish Corollary 1 as below.
Corollary 1.
For system (1), if G ( x ) F ( x ) is always true, then every solution converges to a constant vector as t tends to infinity.
Remark 4.
Clearly, it follows Corollary 1 that the high-dimensional Haddock conjecture is true. Our main results in this paper improve and generalize the corresponding results obtained in [16,17,19,20,21,22]. In particular, our proofs are quite different from those of [13,14,15,16,17,19,20,21,22]. Furthermore, it should be noted that system (1) and its generalizations have been used by extensive recent literature studies [23,24,25,26,27,28,29,30] to simulate the dynamical behaviors of delayed neural network models. However, these literature studies treat the global Lipschz condition as the basic additional assumption to the functions F and G, which indicates that the results of this present paper are completely different from the corresponding conclusions in the above literature.
Remark 5.
It is worth noting that the critical case of system (1) with c i 1 for i Ω has been investigated in the literature [31] under the relatively strict additional conditions(S + )or(S ). The differences between our paper and the literature [31] mainly lie in the following two aspects. On the one hand, if one choose G ( x ) F ( x ) = x 1 3 designed in the Haddock conjecture, simple verification can guide that the function F and G do not meet the requirements of conditions(S + )or(S ). Obviously, the investigation in our paper removes the assumptions(S + )and(S )and overcomes its shortcomings. On the other hand, the main results are completely different; the ω limit set in literature [31] may be periodic solutions, which cannot ensure that each solution converges to a constant vector, but our conclusions clearly point out that each solution must converge to a constant vector. Therefore, on a broader level, our conclusion can be regarded as a high-dimensional version of the famous Haddock conjecture. Since the assumption of a monotonic function in biology is rather vague, and most stochastic dynamics do not always tend towards a strict monotonicity, it would be interesting to further study the deviations of the conjecture in such typical cases. This is our future research direction.

Author Contributions

Formal analysis, C.H. and W.W.; methodology, L.L. and X.Y.; writing—original draft preparation, C.H. and W.W.; writing—review and editing, C.H., W.W., L.L. and X.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the National Natural Science Foundation of China (No. 11971076).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The authors would like to thank the reviewers and the editor for the helpful suggestions and comments, which led to improvements of our original paper.

Conflicts of Interest

We confirm that we have no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
NFDENeutral functional differential equation
IffIf and only if

References

  1. Wu, J.; Zhang, X. Transmission Dynamics of Tick-Borne Diseases with Co-Feeding, Developmental and Behavioural Diapause; Lecture Notes on Mathematical Modelling in the Life Sciences; Springer: Cham, Switzerland, 2020. [Google Scholar]
  2. Zhu, Q. Stabilization of stochastic nonlinear delay systems with exogenous disturbances and the event-triggered feedback control. IEEE Trans. Autom. Control 2019, 64, 3764–3771. [Google Scholar] [CrossRef]
  3. Zhu, Q. Stability analysis for a class of stochastic delay nonlinear systems driven by G-Brownian motion. Syst. Control Lett. 2020, 140, 104699. [Google Scholar] [CrossRef]
  4. Ding, K.; Zhu, Q. Extended dissipative anti-disturbance control for delayed switched singular semi-Markovian jump systems with multi-disturbance via disturbance observer. Automatica 2021, 18, 109556. [Google Scholar] [CrossRef]
  5. Hale, J.K.; Lunel, S.M.V. Introduction to Functional Differential Equations; Springer: New York, NY, USA, 1993. [Google Scholar]
  6. Tan, Y. Dynamics analysis of mackey-glass model with two variable delays. Math. Biosci. Eng. 2020, 17, 4513–4526. [Google Scholar] [CrossRef] [PubMed]
  7. Huang, C.; Tan, Y. Global behavior of a reaction-diffusion model with time delay and Dirichlet condition. J. Differ. Equ. 2021, 271, 186–215. [Google Scholar] [CrossRef]
  8. Long, X. Novel stability criteria on a patch structure Nicholson’s blowflies model with multiple pairs of time-varying delays. AIMS Math. 2020, 5, 7387–7401. [Google Scholar] [CrossRef]
  9. Tan, Y.; Huang, C.; Sun, B.; Wang, T. Dynamics of a class of delayed reaction-diffusion systems with Neumann boundary condition. J. Math. Anal. Appl. 2018, 458, 1115–1130. [Google Scholar] [CrossRef]
  10. Hu, H.; Yi, T.; Zou, X. On spatial-temporal dynamics of a fisher-kpp equation with a shifting environment. Proc. Am. Math. Soc. 2020, 148, 213–221. [Google Scholar] [CrossRef]
  11. Hu, H.; Zou, X. Existence of an extinction wave in the fisher equation with a shifting habitat. Proc. Am. Math. Soc. 2017, 145, 4763–4771. [Google Scholar] [CrossRef]
  12. Long, X.; Gong, S. New results on stability of Nicholson’s blowflies equation with multiple pairs of time-varying delays. Appl. Math. Lett. 2020, 2020, 106027. [Google Scholar] [CrossRef]
  13. Haddock, J. Functional differential equations for which each constant function is a solution: A narrative. In Proceedings International Conference on Nonlinear Oscillations; Farkas, M., Kertesz, V., Stepan, G., Eds.; Janos Bolyai Math. Soc.: Budapest, Hungary, 1987; pp. 86–93. [Google Scholar]
  14. Yi, T. A generalization of the Haddock conjecture and its proof. Nonlinear Anal. Real World Appl. 2008, 9, 1112–1118. [Google Scholar] [CrossRef]
  15. Xiao, S. A New Generalization of the Haddock Conjecture. Math. Model. Anal. 2017, 22, 634–642. [Google Scholar] [CrossRef]
  16. Xiong, W.; Liu, B. Asymptotic behavior of bounded solutions for a class of systems of neutral functional differential equations. J. Math. Anal. Appl. 2006, 313, 754–760. [Google Scholar] [CrossRef] [Green Version]
  17. Xiong, W.; Wang, L. Corrigendum to “Asymptotic behavior of bounded solutions for a system of neutral functional differential equations” [J. Math. Anal. Appl. 313 (2006) 754–760]. J. Math. Anal. Appl. 2009, 350, 423. [Google Scholar] [CrossRef] [Green Version]
  18. Yang, H. Weighted pseudo almost periodicity on neutral type CNNs involving multi-proportional delays and D operator. AIMS Math. 2021, 6, 1865–1879. [Google Scholar] [CrossRef]
  19. Meng, J. Asymptotic behavior of bounded solutions for a system of neutral functional differential equations. Nonlinear Anal. Real World Appl. 2010, 11, 1634–1639. [Google Scholar] [CrossRef]
  20. Gao, H.; Long, H. Dynamic behavior for a system of neutral functional differential equations. Comput. Math. Appl. 2009, 58, 2046–2050. [Google Scholar] [CrossRef] [Green Version]
  21. Xu, M.; Yuan, Z.; Wang, W.; Long, H. ω limit sets of solutions for a class of neutral functional differential equations. Nonlinear Anal. Real World Appl. 2010, 11, 2345–2349. [Google Scholar] [CrossRef]
  22. Kong, F.; Luo, Z. Asymptotic behavior of bounded solutions to a system of neutral functional differential equations in critical case. Appl. Math. Lett. 2018, 81, 44–49. [Google Scholar] [CrossRef]
  23. Xiao, S. Global exponential convergence of HCNNs with neutral type proportional delays and D operator. Neural Process. Lett. 2019, 49, 347–356. [Google Scholar] [CrossRef]
  24. Xu, C.; Li, P.; Xiao, Q.; Yuan, S. New results on competition and cooperation model of two enterprises with multiple delays and feedback controls. Bound. Value Probl. 2019, 2019, 36. [Google Scholar] [CrossRef]
  25. Wang, J.; Chen, X.; Huang, L. The number and stability of limit cycles for planar piecewise linear systems of node-saddle type. J. Math. Anal. Appl. 2019, 469, 405–427. [Google Scholar] [CrossRef]
  26. Duan, L.; Huang, L.; Guo, Z.; Fang, X. Periodic attractor for reaction-diffusion high-order Hopfield neural networks with time-varying delays. Comput. Math. Appl. 2017, 73, 233–245. [Google Scholar] [CrossRef]
  27. Cao, Q.; Long, X. New convergence on inertial neural networks with time-varying delays and continuously distributed delays. AIMS Math. 2020, 5, 5955–5968. [Google Scholar] [CrossRef]
  28. Yang, C.; Huang, L.; Li, F. Exponential synchronization control of discontinuous nonautonomous networks and autonomous coupled networks. Complexity 2018, 2018, 6164786. [Google Scholar] [CrossRef]
  29. Wang, H.; Liu, L.; Tan, Y. Dynamics of Third-Order Nonlinear Neutral Equations. Discrete Dyn. Nat. Soc. 2014, 2014, 1–4. [Google Scholar] [CrossRef]
  30. Wang, J.; Huang, L. Limit cycles bifurcated from a focus-fold singularity in general piecewise smooth planar systems. J. Differ. Equ. 2021, 304, 491–519. [Google Scholar] [CrossRef]
  31. Zhang, X.; Hu, H. Convergence in a system of critical neutral functional differential equations. Appl. Math. Lett. 2020, 107, 106385. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Wang, W.; Li, L.; Yi, X.; Huang, C. Convergence on Population Dynamics and High-Dimensional Haddock Conjecture. Symmetry 2021, 13, 2252. https://doi.org/10.3390/sym13122252

AMA Style

Wang W, Li L, Yi X, Huang C. Convergence on Population Dynamics and High-Dimensional Haddock Conjecture. Symmetry. 2021; 13(12):2252. https://doi.org/10.3390/sym13122252

Chicago/Turabian Style

Wang, Wenke, Le Li, Xuejun Yi, and Chuangxia Huang. 2021. "Convergence on Population Dynamics and High-Dimensional Haddock Conjecture" Symmetry 13, no. 12: 2252. https://doi.org/10.3390/sym13122252

APA Style

Wang, W., Li, L., Yi, X., & Huang, C. (2021). Convergence on Population Dynamics and High-Dimensional Haddock Conjecture. Symmetry, 13(12), 2252. https://doi.org/10.3390/sym13122252

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop