Next Article in Journal
Nearest Neighbor Decoding and Pilot-Aided Channel Estimation for Fading Channels
Next Article in Special Issue
Information Bottleneck Analysis by a Conditional Mutual Information Bound
Previous Article in Journal
On-The-Fly Syntheziser Programming with Fuzzy Rule Learning
Previous Article in Special Issue
Design and Practical Stability of a New Class of Impulsive Fractional-Like Neural Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Lotka-Volterra-Type Cooperation Models: Impulsive Control on Their Stability Behavior

1
Department of Biology, University of Texas at San Antonio, San Antonio, TX 78249, USA
2
Department of Mathematics, University of Texas at San Antonio, San Antonio, TX 78249, USA
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Submission received: 1 August 2020 / Revised: 26 August 2020 / Accepted: 27 August 2020 / Published: 31 August 2020
(This article belongs to the Special Issue Dynamics in Complex Neural Networks)

Abstract

:
We present a biological fractional n-species delayed cooperation model of Lotka-Volterra type. The considered fractional derivatives are in the Caputo sense. Impulsive control strategies are applied for several stability properties of the states, namely Mittag-Leffler stability, practical stability and stability with respect to sets. The proposed results extend the existing stability results for integer-order n species delayed Lotka-Volterra cooperation models to the fractional-order case under impulsive control.

1. Introduction

The dynamics of Lotka-Volterra and related systems are traditionally applied to problems in population dynamics [1,2,3,4]. For example, the work in Reference [5] studied a delayed cooperation ecosystem composed of a number (n) of species using the following Lotka-Volterra model:
x ˙ i ( t ) = r i ( t ) x i ( t ) 1 x i ( t τ i i ) a i ( t ) + j = 1 j i n b j ( t ) x j ( t τ i j ) c i ( t ) x i ( t τ i j ) ,
where i , j = 1 , , n ; t 0 , x i ( t ) represents the density of the species i at time t, r i ( t ) are the intrinsic growth rates at time t, a i ( t ) , b i ( t ) and c i ( t ) are the cooperative coefficients, τ i i denote the time delay due to gestation period of species i, τ i j , j i denote time delays related to the maturation of the species i and j, respectively, 0 τ i j τ , τ = c o n s t , ( i = 1 , 2 , , n ) . It is assumed that only adult individuals of species i can benefit the species j and vice versa.
The model (1) has been originally introduced in Reference [6] for n = 2 and τ i j = 0 , i , j = 1 , 2 , to describe the cooperative interaction between two species. Since then, it has been developed by several authors. See, for example, References [7,8,9].
An important objective in understanding biological systems is to determine their approach to equilibrium. Relative entropy, part of an information-theoretic approach, can be used to determine approach to equilibrium in biological systems, including the Lotka-Volterra type [10]. In addition, entropy methods are increasingly applied to biological models as robustness markers [11], as well as to determine information content in artificial neural networks [12].
Research on impulsive Lotka-Volterra and related systems has also increased because of their application to the study of population densities that are affected by some impulsive factors that cannot be ignored. Impulsive control strategies are also widely used in control theory of biological models [13,14,15,16,17,18,19,20]. In Reference [21] the author extended model (1) to the impulsive case and proved that the population dynamics of the system can be impulsively controlled. Furthermore, appropriate exterior or interior short-term changes in population densities may be used as controls. In general, impulsive controllers and impulsive control techniques demonstrate their efficiency in stability and synchronization of numerous applied models [22,23,24].
Fractional-order derivatives are increasingly being used to better reflect long-term dependence properties in biological systems [25,26,27,28]. Fractional Lotka-Voltera-type models have also been recently developed [29,30,31,32,33,34]. Indeed, the use of fractional derivatives and integrals has proven advantages in many applications across engineering and sciences [35,36,37]. In addition, the extended impulsive fractional-order case has also attracted the attention of many researchers [38,39,40]. Impulsive effects in fractional biological models are considered by several authors [38,39,41,42,43], including some fractional Lotka-Volterra models [39,44]. However, most of the studies considered competitive models, and to the best of our knowledge, there are no results on fractional cooperative models.
The asymptotic stability of equilibrium points is one of the most important tasks in the study of biological network systems. Asymptotic stability theory has been applied to a wide range of integer-order Lotka-Volterra systems to explain the behavior of the system’s solutions using stabilization control law [1,2,3,5,7,8,9,21,23]. For fractional-order models the concept of Mittag-Leffler stability, introduced in Reference [45], is adapted as a generalization of the exponential stability notion in integer-order systems [46,47,48]. Very recently, the Mittag-Leffler notion has been also applied to fractional-order control problems [49,50,51], including impulsive controls in neural network systems [47,52,53,54]. The Mittag-Leffler stability and impulsive control strategies have not been investigated for fractional generalizations of the cooperative Lotka-Volterra model (1) which are important and challenging problems and one of the main goals of the present research.
The concepts of practical stability [55,56] and stability of sets [57] have been applied to fractional-order systems [58,59,60,61]. Using a practical stability approaches is beneficial in the cases when the equilibrium of the considered system is unstable in the classical Lyapunov or Mittag-Leffler sense, but its performance may be sufficient for a particular application [56]. Stability of sets generalizes the concept of stability of a single solution and is suitable to numerous real-world applications [62,63]. However, due to the complexity of the problem, no results on practical stability and stability of sets properties have been reported so far for fractional Lotka-Volterra-type models.
Stimulated by the above discussion, in this paper, we study impulsive control strategies for the Mittag-Leffler stability, practical stability, and stability of sets behavior of the fractional-order modification of the n species delayed cooperation Lotka-Volterra type model (1). Caputo’s fractional-order approach is adopted in our research. In Section 2 we will first introduce the fractional generalization of the biological network model (1). Some definitions and lemmas are also presented. Section 3 is devoted to our Mittag-Leffler stability results via impulsive controllers. In Section 4 the practical stability results for the fractional impulsive control biological model are derived. Section 5 generalizes the results in Section 3 and Section 4, and offers stability of sets results. In Section 6 examples and simulations are given. Section 7 presents our conclusions and future directions.

2. Fractional Lotka-Volterra-Type Model and Preliminary Notes

We will use the following notations: R + = [ 0 , ) , R n is the n-dimensional Euclidean space and | | x | | = | x 1 | + + | x n | denote the norm of x R n .
Let u C 1 [ [ 0 , b ] , R ] , b > 0 . For any 0 < q < 1 , following Reference [36], we consider the Caputo fractional derivative of order q with the lower limit 0 for a function u, defined as
C D t q u ( t ) = 1 Γ ( 1 q ) 0 t u ( σ ) ( t σ ) q d σ , t 0 ,
where Γ is the Gamma function and u is the first-order ordinary derivative of the function u.
Since introducing fractional-order derivatives in population dynamics models has important practical significance and research values [25,26,27,28,29,30,31,32,33,34], we generalize the model (1), and introduce the following model with fractional-order derivatives
C D t q x i ( t ) = r i ( t ) x i ( t ) 1 x i ( t τ i i ) a i ( t ) + j = 1 j i n b j ( t ) x j ( t τ i j ) c i ( t ) x i ( t τ i j ) ,
where i = 1 , 2 , , n , n 2 , 0 < q < 1 , the functions r i ( t ) , a i ( t ) , b i ( t ) and c i ( t ) are continuous, positive and bounded on R + .
The fractional derivatives in the Lotka-Volterra-type cooperation model (2) will contribute for a better description of the population interactions characteristics and will improve the degree of freedom in comparison with the model (1) [35,36].
We will consider the following initial condition
x i ( s ; 0 , φ 0 ) = φ i 0 ( s ) , s [ τ , 0 ]
associated with the model (2), where the initial function φ 0 : [ τ , 0 ] R n , φ 0 = ( φ 10 , φ 20 , , φ n 0 ) T is continuous, positive and bounded on [ τ , 0 ] .
For a continuous function g ( t ) defined on [ 0 , ) , we denote
g L = inf t [ 0 , ) g ( t ) , g M = sup t [ 0 , ) g ( t ) .
The goal of our research is to investigate the effects of impulsive external perturbations on the stability behavior of the model (2) at some fixed points t 1 , t 2 , such that
0 < t 1 < t 2 < < t k <
and lim k t k = . To this end we will add impulsive controllers that can be used to synchronize the densities of all species onto that of system (2), and introduce the following impulsive control system
C D t q y i ( t ) = r i ( t ) y i ( t ) 1 y i ( t τ i i ) a i ( t ) + j = 1 j i n b j ( t ) y j ( t τ i j ) c i ( t ) y i ( t τ i j ) + η i ( y i ( t ) ) ,
where
η i ( y i ( t ) ) = k = 1 I i k ( y i ( t ) ) δ ( t t k )
are the control inputs, δ ( t ) is the Dirac impulsive function, y i : R + R + , I i k : R + R , i = 1 , , n , k = 1 , 2 , .
Due to the incorporation of the Delta impulsive function the controller η ( t ) has an effect on sudden changes in the states of the model (4) at the time instants t k , that is, η ( t ) is an impulsive control of (4). For t > 0 and i = 1 , 2 , , n , we denote
f i ( t , y i ( t ) ) = r i ( t ) y i ( t ) 1 y i ( t τ i i ) a i ( t ) + j = 1 j i n b j ( t ) y j ( t τ i j ) c i ( t ) y i ( t τ i j ) .
From (4), we have that [35,36,39,61], η i ( y i ( t ) ) = 0 for t t k , k = 1 , 2 , . , i = 1 , 2 , , n , and
y i ( t k + h ) y i ( t k ) = 1 Γ ( q ) 0 t k + h ( t k + h σ ) q 1 [ f i ( σ , y i ( σ ) ) + η i ( y i ( σ ) ) ] d σ
1 Γ ( q ) 0 t k ( t k σ ) q 1 [ f i ( σ , y i ( σ ) ) + η i ( y i ( σ ) ) ] d σ
= 1 Γ ( q ) t k t k + h ( t k + h σ ) q 1 [ f i ( σ , y i ( σ ) ) + η i ( y i ( σ ) ) ] d σ
+ 1 Γ ( q ) 0 t k ( t k + h σ ) q 1 ( t k σ ) q 1 [ f i ( σ , y i ( σ ) ) + η i ( σ ) ] d σ ,
where h > 0 is sufficiently small. As h 0 + , we obtain
Δ y i ( t k ) = y i ( t k + ) y i ( t k ) = I i k ( y i ( t k ) ) ,
where y i ( t k + ) = lim h 0 + y i ( t k + h ) .
Therefore, the correspondent n-dimensional fractional-order driven system is given by
C D t q y i ( t ) = y i ( t ) r i ( t ) 1 y i ( t τ i i ) a i ( t ) + j = 1 j i n b j ( t ) y j ( t τ i j ) c i ( t ) y i ( t τ i j ) , t t k , y i ( t k + ) = y i ( t k ) + I i k ( y i ( t k ) ) , i = 1 , , n , k = 1 , 2 , ,
where the quantities y i ( t k ) and y i ( t k + ) = lim h 0 + y i ( t k + h ) represent, respectively, the population densities of species i before and after an impulsive perturbation at the moment t k and the continuous in R + functions I i k characterize the impact of the impulse effect on the species i at the moments t k .
For more results on fractional-order models under impulsive control where the Dirac delta function is incorporated in a similar sense, we refer the reader to References [39,61] and the references therein.
Let ϕ 0 : [ τ , 0 ] R n , ϕ 0 = ( ϕ 10 , ϕ 20 , , ϕ n 0 ) T be a piecewise continuous function with points of discontinuity on [ τ , 0 ] at which it is continuous from the left. The class of all such functions will be denoted by PC , and by P C B [ [ τ , 0 ] , R n ] we will denote the class of all bounded functions ϕ 0 PC . We denote by y ( t ) = y ( t ; 0 , ϕ 0 ) = ( y 1 ( t ; 0 , ϕ 0 ) , y 2 ( t ; 0 , ϕ 0 ) , , y n ( t ; 0 , ϕ 0 ) ) T the solution of system (5), satisfying the initial conditions
y i ( s ; 0 , ϕ 0 ) = ϕ i 0 ( s ) , s [ τ , 0 ] , y i ( 0 + ; 0 , ϕ 0 ) = ϕ i 0 ( 0 ) , i = 1 , , n .
Furthermore, from the biological point of view, we will restrict our attention only to the nonnegative solutions of the control system (5), and we will assume that if ϕ i 0 ( s ) 0 , sup ϕ i 0 ( s ) < , s [ τ , 0 ] , ϕ i 0 ( 0 ) > 0 , 1 i n , then y i ( t ) 0 , 1 i n , t [ 0 , ) . We also assume that y i + I i k ( y i ) 0 for y i R + , i = 1 , 2 , , n , k = 1 , 2 , .
For some results on positive solutions of fractional Lotka-Volterra-type models we refer to References [32,33].
In addition, we assume that, if the functions I i k are such that y i I i k ( y i ) 0 for y i R + , i = 1 , 2 , , n , k = 1 , 2 , , then there exist positive constants m and M < such that
m y i ( t ) M , t [ 0 , ) .
For integer-order Lotka-Volterra models with and without impulses, similar results can be found in References [1,5,7,8,9,21,23]. Since one of the main advantages of Caputo’s fractional derivatives is that the initial conditions for fractional differential equations with such derivatives have a similar form as for integer-order differential equations, the validity of (7) can be verified using similar steps as in the integer-order models.
Define the norm | | ϕ | | τ = sup s [ τ , 0 ] | | ϕ ( s ) | | of the function ϕ P C B [ [ τ , 0 ] , R n ] .
It is well known [1,2,4,7,8,9,13,21,23] that Lyapunov functions can be a successful instrument in the investigation of the stability properties of a Lotka-Volterra type system. To apply the Lyapunov function method, we introduce the following notations:
G k = ( t k 1 , t k ) × R + n , k = 1 , 2 , ; G = k = 1 G k ;
V 0 = V : [ 0 , ) × R + n R + : V C [ G , R + ] , t [ 0 , ) , V is locally Lipschitzian in y R + n on each of the sets G k , V ( t k , y ) = V ( t k , y ) and V ( t k + , y ) = lim t t k t > t k V ( t , y ) exists}.
For V V 0 and t [ t k 1 , t k ) , we will use the following fractional derivative of V of order q, 0 < q < 1 with respect to system (1) defined [39] by
C D + q V ( t , ϕ ( 0 ) ) = lim h 0 + sup 1 h q [ V ( t , ϕ ( 0 ) ) V ( t h , ϕ ( 0 ) h q f ( t , ϕ ) ) ] ,
where f ( t , ϕ ) = ( f 1 ( t , ϕ ) , f 2 ( t , ϕ ) , , f n ( t , ϕ ) ) T denotes the right-hand side of (5), ϕ P C B [ [ τ , 0 ] , R n ] ,
f i ( t , ϕ ) = ϕ i ( t ) r i ( t ) 1 ϕ i ( t τ i i ) a i ( t ) + j = 1 j i n b j ( t ) ϕ j ( t τ i j ) c i ( t ) ϕ i ( t τ i j ) .
We also recall the following basic comparison result of fractional calculus from Reference [39].
Lemma 1.
Assume that the function V V 0 is such that for t R + , ϕ PC ,
V ( t + , ϕ ( 0 ) + Δ ϕ ) V ( t , ϕ ( 0 ) ) , t = t k ,
and for a constant μ > 0 the inequality
C D + q V ( t , ϕ ( 0 ) ) μ V ( t , ϕ ( 0 ) ) , t t k , k = 1 , 2 ,
is valid whenever V ( t + s , ϕ ( s ) ) V ( t , ϕ ( 0 ) ) for τ s 0 .
Then s u p τ s 0 V ( s , ϕ 0 ( s ) ) V ( t , ϕ ( 0 ) ) implies
V ( t , y ( t ; 0 , ϕ 0 ) ) sup τ s 0 V ( 0 + , ϕ 0 ( s ) ) E q ( μ t q ) , t [ 0 , ) ,
where E q is the Mittag-Leffler function defined as E q ( z ) = κ = 0 z κ Γ ( q κ + 1 ) , q > 0 .
For more comparison results and properties on Mittag-Leffler functions we refer to References [35,36,37,39].

3. Mittag-Leffler Stability Results

Let ϕ 0 * P C B [ [ τ , 0 ] , R + n ] , ϕ 0 * = ( ϕ 10 * , ϕ 20 * , , ϕ n 0 * ) T and y * ( t ) = y * ( t ; 0 , ϕ 0 * ) = ( y 1 * ( t ; 0 , ϕ 0 * ) , y 2 * ( t ; 0 , ϕ 0 * ) , , y n * ( t ; 0 , ϕ 0 * ) ) T be the solution of the impulsive control system (5) with initial function ϕ 0 * . In the next, we shall suppose that ϕ i 0 * ( s ) 0 , sup ϕ i 0 * ( s ) < , ϕ i 0 * ( 0 ) > 0 , i = 1 , 2 , , n .
We will use the following definition for Mittag-Leffler stability of the state y * ( t ) of the model (5) [45,46,47,48,49,50,51,52,53,54].
Definition 1.
The state y * ( t ) of system (5) is said to be globally Mittag-Leffler stable, if
| | y ( t ; 0 , ϕ 0 ) y * ( t ; 0 , ϕ 0 * ) | | { m [ | | ϕ 0 ϕ 0 * | | τ ] E q ( μ t q ) } d ,
where E q is the corresponding Mittag-Leffler function, q ( 0 , 1 ) , μ > 0 , d > 0 , m ( 0 ) = 0 , m ( ϕ ) 0 , and m ( ϕ ) is locally Lipschitz with respect to ϕ P C B [ [ τ , 0 ] , R n ] .
In some of the literature [50], the constant μ is also called the degree of Mittag-Leffler state estimator, which can be considered as an equivalence of the convergence rate as state estimator error tends to zero when time t goes to infinity.
Our main result in this section is the next theorem.
Theorem 1.
Assume that:
1. 
The model’s parameters satisfy
2 m min 1 i n r i L a i M + M j = 1 j i n b j M + min 1 i n r i L c i L > max 1 i n r i M + max 1 i n j = 1 j i n b j M r i M M 2 a i L + m κ = 1 κ j n b κ L 2 .
2. 
The impulsive functions I i k are such that
I i k ( y i ( t k ) ) = γ i k y i ( t k ) y i k * ( t k ) , 0 < γ i k < 2 , i = 1 , 2 , n , k = 1 , 2 , .
Then, the state y * ( t ) of system (5) is globally Mittag-Leffler stable.
Proof. 
The proof is based on the Lyapunov-Razumikhin technique [39]. We consider the following Lyapunov function candidate
V ( t , y ) = i = 1 n | y i y i * | .
For k = 1 , , and ϕ P C B , according to (8), we have
V ( t k + , ϕ ( 0 ) + I k ( ϕ ) ) = i = 1 n | ϕ i 0 ( 0 ) + I i k ( ϕ 0 ) ϕ i 0 * ( 0 ) I i k ( ϕ 0 * ) | i = 1 n | ( 1 γ i k ) ( ϕ i 0 ( 0 ) ϕ i 0 * ( 0 ) ) | V ( t k , ϕ ( 0 ) ) .
Also, for t 0 , t t k , k = 1 , 2 , , and ϕ P C B , we estimate the derivative C D + q V ( t , ϕ ( 0 ) ) along the system (5), and get that
C D + q V ( t , ϕ ( 0 ) ) max 1 i n r i M i = 1 n | ϕ i 0 ( 0 ) ϕ i 0 * ( 0 ) | min 1 i n 2 r i L c i L m i = 1 n | ϕ i 0 ( 0 ) ϕ i 0 * ( 0 ) | min 1 i n 2 r i L m a i M + M j = 1 j i n b j M i = 1 n | ϕ i 0 ( 0 ) ϕ i 0 * ( 0 ) | + max 1 i n j = 1 j i n b j M r i M M 2 a i L + m κ = 1 κ j n b κ L 2 i = 1 n sup s [ τ , 0 ] | ϕ i 0 ( s ) ϕ i 0 * ( s ) | ( 2 m min 1 i n r i L a i M + M j = 1 j i n b j M + min 1 i n r i L c i L + max 1 i n r i M + max 1 i n j = 1 j i n b j M r i M M 2 a i L + m κ = 1 κ j n b κ L 2 ) V ( t , φ ( 0 ) )
is valid whenever V ( t + s , φ ( s ) ) V ( t , φ ( 0 ) ) for τ s 0 .
In view of (10) and condition 1 of Theorem 1, there exists a positive μ such that
C D + q V ( t , ϕ ( 0 ) ) μ V ( t , ϕ ( 0 ) ) , t t k , k = 1 , 2 , ,
whenever V ( t + s , ϕ ( s ) ) V ( t , ϕ ( 0 ) ) for τ s 0 .
From (9) and (11), by Lemma 1, it follows that
V ( t , y ( t ; 0 , ϕ 0 ) ) sup τ s 0 V ( 0 + , ϕ 0 ( s ) ) E q ( μ t q ) , t [ 0 , ) .
So,
| | y ( t ) y * ( t ) | | | | ϕ 0 ϕ 0 * | | τ E q ( μ t q ) , t [ 0 , ) ,
and this completes the proof of the theorem. □
Remark 1.
The assumptions of Theorem 1 guarantee the Mittag-Leffler stability of the states of the impulsive control system (5). The impulsive control strategy is designed by (8). Therefore, the proposed structure of the control law can be applied to design impulsive based observers to estimate the states of (5) and globally Mittag-Leffler synchronize the states of the system (5) onto that of system (2). Note that, the synchronization scheme is independent on the lengths of the impulsive intervals.
Remark 2.
One of the states that is of a great interest to researchers of Lotka-Volterra and related models are the so called “steady" or equilibrium states [1,3,21,47]. Let
y * = ( y 1 * , y 2 * , , y n * ) T
be the equilibrium of the control system (5) such that
1 = y i * a i ( t ) + j = 1 j i n b j ( t ) y j * c i ( t ) y i * , t t k ,
I i k ( y i * ) = 0 , i = 1 , 2 , , n , k = 1 , 2 , .
Since I i k ( y i * ) = 0 , i = 1 , 2 , , n , k = 1 , 2 , , the state y * is also a steady state for the controlled model (2). Hence, in the case when y * ( t ) = y * = c o n s t , the providing results in Theorem 1 can be used so by means of convenient impulsive control strategy to keep the Mittag-Leffler stability properties of the equilibrium states of system (2). Note that the Mittag-Leffler stability for fractional-order models corresponds to the exponential stability for integer-order systems, which guarantees the fast convergence rate.

4. Practical Stability Results

In this section, following Reference [56] we will provide practical stability results for the impulsive control model (5).
We introduce the following definition.
Definition 2.
The state y * ( t ) of the fractional impulsive control system (5) is said to be:
(a) practically uniformly globally stable with respect to ( λ , A ) , if given ( λ , A ) with 0 < λ < A , we have | | ϕ 0 ϕ 0 * | | τ < λ implies
| | y ( t ; 0 , ϕ 0 ) y * ( t ; 0 , ϕ 0 * ) | | A , t 0 ;
(b) practically uniformly globally Mittag-Leffler stable with respect to ( λ , A ) , if given ( λ , A ) with 0 < λ < A , we have | | ϕ 0 ϕ 0 * | | τ < λ implies
| | y ( t ; 0 , ϕ 0 ) y * ( t ; 0 , ϕ 0 * ) | | A + ν | | ϕ 0 ϕ 0 * | | τ E q [ μ t q ] , t 0 , ν , μ > 0 .
Condition 1 of Theorem 1 guarantees the existence of a positive constant μ such that (11) holds. The application of (11) together with the appropriate choice of the Lyapunov-candidate function V imply the global Mittag-Leffler stability of the state y * ( t ) of the control system (5) by means of the designed impulsive control law (8). This section will firstly address the case, when μ = 0 and the state y * ( t ) is not globally Mittag-Leffler stable. We will prove that in such case the impulsive controller (8) may assure its practical uniform global stability.
We will use again the Lyapunov function candidate
V ( t , y ) = i = 1 n | y i y i * |
and functions of the following class:
K = { a C [ R + , R + ] : a ( u ) is strictly increasing in u and a ( 0 ) = 0 }.
Theorem 2.
Assume that the impulsive controller satisfies (8) and:
1. 0 < λ < A are given.
2. The model’s parameters satisfy
2 m min 1 i n r i L a i M + M j = 1 j i n b j M + min 1 i n r i L c i L = max 1 i n r i M + max 1 i n j = 1 j i n b j M r i M M 2 a i L + m κ = 1 κ j n b κ L 2 .
3. There exists a function b K such that
V ( t , y ) b ( | | y y * | | ) , t R + , y , y * R n .
4. b ( λ ) < A holds.
Then, the state y * ( t ) of system (5) is practically uniformly globally stable with respect to ( λ , A ) .
Proof. 
Let 0 < λ < A . We will prove that | | ϕ 0 ϕ 0 * | | τ < λ implies | | y ( t ; 0 , ϕ 0 ) y * ( t ; 0 , ϕ 0 * ) | | A , t 0 .
If this is not true, there exists a corresponding solution y ( t ; 0 , ϕ 0 ) of (5) with | | ϕ 0 ϕ 0 * | | τ < λ , and 0 < t * < , t k < t * t k + 1 for some k, such that
| | y ( t ; 0 , ϕ 0 ) y * ( t ; 0 , ϕ 0 * ) | | < A , t [ 0 , t k ] and | | y ( t * ) y * ( t * ) | | A .
Note that the impulsive control low (8) guarantees that if
| | y ( t k ; 0 , ϕ 0 ) y * ( t k ; 0 , ϕ 0 * ) | | < A ,
then
| | y ( t k + ; 0 , ϕ 0 ) y * ( t k + ; 0 , ϕ 0 * ) | | = i = 1 n | y i ( t k ) + I i k ( y i ( t k ) ) y i * ( t k ) I i k ( y i * ( t k ) ) |
i = 1 n | ( 1 γ i k ) ( y i ( t k ) y i * ( t k ) ) | < A .
Hence t * t k for any k = 1 , 2 , .
From condition 2 of Theorem 2 by similar arguments as in the proof of Theorem 1, we have
C D + q V ( t , ϕ ( 0 ) ) 0 , t t k , k = 1 , 2 , ,
whenever V ( t + s , ϕ ( s ) ) V ( t , ϕ ( 0 ) ) for τ s 0 .
Also, using (12) and (8), we obtain by Lemma 1, for μ = 0
V ( t , y ( t ; 0 , ϕ 0 ) ) sup τ s 0 V ( 0 + , ϕ 0 ( s ) ) , t [ 0 , t * ] .
Hence, (13), (8) and condition 4 of Theorem 2 imply
A V ( t * , y ( t * ; 0 , ϕ 0 ) ) sup τ s 0 V ( 0 + , ϕ 0 ( s ) ) b ( λ ) < A .
The contradiction obtained proves our claim that system (5) is practically uniformly globally stable with respect to ( λ , A ) . □
Remark 3.
In Theorem 2 we proved that even in some concrete cases, when the system (5) is not globally Mittag-Leffler stable, by the impulsive controller (8) we can made the vicinity of a state y * arbitrarily small so that all other states to remain arbitrarily close to it provided that their initial functions were sufficiently near the initial data of the particular state y * . The property of practical stability is very important in many applications when the Mittag-Leffler stability is conservative [58,60,61]. For type (2) models from the population biology this can compensate a not stable initial population by keeping the population densities between particular bounds for t 0 .
Theorem 3.
Assume that:
1. There exists a function V V 0 for which conditions of Lemma 1 are met.
2. The function V ( t , y ) is such that
| | y y * | | A V ( t , y ) a ( u ) | | y y * | | , t R + , y , y * R n ,
where the continuous function a ( u ) 1 exists for any u > 0 .
3. Condition 1 of Theorem 2 holds.
Then, the state y * ( t ) of system (5) is practically uniformly globally Mittag-Leffler stable with respect to ( λ , A ) .
Proof. 
Let 0 < λ < A , and | | ϕ 0 ϕ 0 * | | τ < λ for an initial function ϕ 0 P C B [ [ τ , 0 ] , R n ] of an arbitrary solution y ( t ) = y ( t ; 0 , ϕ 0 ) of the initial value problem (5), (6).
From Lemma 1, we have that
V ( t , y ( t ; 0 , ϕ 0 ) ) sup τ s 0 V ( 0 + , ϕ 0 ( s ) ) E q ( μ t q ) , t [ 0 , ) ,
where μ > 0 .
From (14) and condition 2 of Theorem 3, we have
| | y ( t ; 0 , ϕ 0 ) y * ( t ; 0 , ϕ 0 * ) | | A V ( t , y ( t ; 0 , ϕ 0 ) )
sup τ s 0 V ( 0 + , ϕ 0 ( s ) ) E q ( μ t q ) a ( u ) | | ϕ 0 ϕ 0 * | | τ E q ( μ t q ) , t 0 .
Therefore, for any ν a ( u ) , u > 0 , we have
| | y ( t ; 0 , ϕ 0 ) y * ( t ; 0 , ϕ 0 * ) | | A + ν | | ϕ 0 ϕ 0 * | | τ E q ( μ t q ) , t 0 ,
which shows that system (5) is practically uniformly globally Mittag-Leffler stable with respect to ( λ , A ) . □
Remark 4.
For A = 0 Theorem 3 implies global Mittag-Leffler stability. However, Theorem 3 can be applied when this stability notion is pointless from the applied point of view. In such cases the designed impulsive controller allows the practical uniform global Mittag-Leffler stability.
Remark 5.
More detailed explanations about the importance of the practical stability notion for applied problems, as well as, for relations between stability and practical stability concepts can be found in References [56,61,64].

5. Stability of Sets

In this Section, we extend the Lyapunov function method to stability of sets criteria for the impulsive control system (5). In fact, considering closed sets of states as attractors instead of equilibrium states is common in population dynamics [1].
Let B R n be an arbitrary set (compact or not compact). Define the norms:
| | y | | B = inf { | | y z | | : z B } is the distance from y R n to B;
| | ϕ | | τ B = sup s [ τ , 0 ] | | ϕ ( s ) | | B is the distance between a function ϕ P C B [ [ τ , 0 ] , R n ] and B.
Definition 3.
The solutions of system (5) are said to be uniformly bounded with respect to the set B, if for any positive constant η > 0 there exists a constant β > 0 such that ϕ 0 P C B [ [ τ , 0 ] , R n ] and | | ϕ 0 | | τ B η imply
| | y ( t ; 0 , ϕ 0 ) | | B β , t 0 .
Definition 4.
The set B is said to be globally Mittag-Leffler stable with respect to system (5), if the solutions of system (5) are uniformly bounded with respect to the set B, and
| | y ( t ; 0 , ϕ 0 ) | | B { m [ | | ϕ 0 | | τ B ] E q ( μ t q ) } d ,
where E q is the corresponding Mittag-Leffler function, q ( 0 , 1 ) , μ > 0 , d > 0 , m ( 0 ) = 0 , m ( ϕ ) 0 , and m ( ϕ ) is locally Lipschitz with respect to ϕ P C B [ [ τ , 0 ] , R n ] .
Remark 6.
The notion of stability of sets with respect to system (5) generalizes the stability of single states concepts of system (5). For example, when B = { y R n : y = y * ( t ) , t 0 } , Definition 3 reduces to Definition 1 for the global Mittag-Leffler stability of the state y * ( t ) , t 0 .
Remark 7.
The requirement for uniform boundedness of all solutions of (5) with respect to B is required by the fact that for noncompact sets B there is an opportunity that some solutions may tend to infinity in finite time within the set. If the set B is compact (for example, the equilibrium state), then the boundedness condition is surplus.
Theorem 4.
Assume that:
1. There exists a function V V 0 such that V ( t , y ) = 0 for y B and V ( t , y ) > 0 for y R + n \ B , t 0 .
2. There exists a function b K such that
| | y | | B V ( t , y ) b ( | | y | | B ) , ( t , y ) R + × R + n .
3. For the function V conditions of Lemma 1 are met.
Then, the set B is globally Mittag-Leffler stable with respect to system (5).
Proof. 
We will first prove the uniform boundedness of the solutions of (5) with respect to the set B. Let η > 0 be given and let ϕ 0 P C B [ [ τ , 0 ] , R n ] be such that | | ϕ 0 | | τ B η . We can choose β > 0 so that b ( η ) < β .
Let y ( t ; 0 , ϕ 0 ) be the solution of the system (5) corresponding to the initial function ϕ 0 .
Since the conditions of Lemma 1 are met, then for μ = 0 , we have
V ( t , y ( t ; 0 , ϕ 0 ) ) sup τ s 0 V ( 0 + , ϕ 0 ( s ) ) , t 0 .
Now, condition 2 of Theorem 4 and (15) imply the inequalities
| | y ( t ; 0 , ϕ 0 ) | | B V ( t , y ( t ; 0 , ϕ 0 ) ) sup τ s 0 V ( 0 + , ϕ 0 ( s ) )
b ( | | ϕ 0 | | τ B ) b ( η ) < β , t 0 .
Therefore,
| | y ( t ; 0 , ϕ 0 ) | | B β , t 0 .
Since the solution y ( t ; 0 , ϕ 0 ) is arbitrary the above proves the uniform boundedness of the solutions of (5) with respect to the set B.
In a similar way, as above, using Lemma 1 and condition 2 of Theorem 4, we obtain
| | y ( t ; 0 , ϕ 0 ) | | B V ( t , y ( t ; 0 , ϕ 0 ) ) sup τ s 0 V ( 0 + , ϕ 0 ( s ) ) E q ( μ t q )
b ( | | ϕ 0 | | τ B ) E q ( μ t q ) , t 0 ,
which according to Definition 4 means that the set B is globally Mittag-Leffler stable set with respect to (5). □
Remark 8.
Theorem 4 shows that the impulsive control law (8) is designed so that gets global Mittag-Leffler stability of a set B with respect to the fractional model (5). For the integer-order systems this is equivalent to the fact that all trajectories arrive at a given target set B. Hence, the proposed result is of a great importance in the cases when the information on the system’s parameters is incomplete and ones has to conduct a robust stability analysis.
Remark 9.
It is well know that if B is globally asymptotically stable with respect to (5), then the trivial solution of (5) is practically stable [65].
Remark 10.
We can also define the hybrid concept of practical stability of an arbitrary set B with respect to the impulsive control system (5). In fact, this concept is more natural than that of a neighbourhood of a state. In this case, one important question will be how well the stability behavior can be impulsively controlled, which is a subject for our future research.

6. Examples and Simulations

Example 1.
Consider the following two-species fractional cooperation Lotka-Volterra system
C D t q x 1 ( t ) = x 1 ( t ) 1 x 1 ( t τ 11 ) 1 + 5 x 2 ( t τ 12 ) 3 x 1 ( t τ 12 ) , C D t q x 2 ( t ) = x 2 ( t ) 1 x 2 ( t τ 22 ) 0.3 + 0.7 x 1 ( t τ 21 ) 3 x 2 ( t τ 21 ) ,
t 0 , 0 τ i j 1 , i , j = 1 , 2 .
We can easily check that the point ( 2 7 , 1 5 ) is an equilibrium for the model (16).
The corresponding fractional control Lotka-Volterra system is given by
C D t q y 1 ( t ) = y 1 ( t ) 1 y 1 ( t τ 11 ) 1 + 5 y 2 ( t τ 12 ) 3 y 1 ( t τ 12 ) , t t k , C D t q y 2 ( t ) = y 2 ( t ) 1 y 2 ( t τ 22 ) 0.3 + 0.7 y 1 ( t τ 21 ) 3 y 2 ( t τ 21 ) , t t k ,
under impulsive control law
Δ y 1 ( t k ) = y 1 ( t k + ) y 1 ( t k ) = 2 5 y 1 ( t k ) 2 7 , k = 1 , 2 , , Δ y 2 ( t k ) = y 2 ( t k + ) y 2 ( t k ) = 1 3 y 2 ( t k ) 1 5 , k = 1 , 2 , ,
where 0 < t 1 < t 2 < and lim k t k = .
We have that for y 1 * = 2 7 and y 2 * = 1 5 , C D t q y 1 * = 0 and C D t q y 2 * = 0 , that is, the state ( 2 7 , 1 5 ) is an equilibrium for the model (17), (18) too.
Also, obviously all conditions of Theorem 1 are satisfied for m = 1 5 , M = 2 7 , a 1 L = a 1 M = 1 , a 2 L = a 2 M = 0.3 , b 1 L = b 1 M = 0.7 , b 2 L = b 2 M = 5 , c i L = c i M = 3 , i = 1 , 2 , γ 1 k = 2 5 , γ 2 k = 1 3 , k = 1 , 2 , . In fact, we have that
1.364 2 m min 1 i 2 r i L a i M + M j = 1 j i n b j M + min 1 i 2 r i L c i L > max 1 i 2 r i M + max 1 i 2 j = 1 j i 2 b j M r i M M 2 a i L + m κ = 1 κ j 2 b κ L 2 1.314 .
Hence, the equilibrium state ( 2 7 , 1 5 ) of the impulsive system (17), (18) is globally Mittag-Leffler stable. This example shows the effectiveness of the obtained results under an impulsive controller (18) of type (8). It also, demonstrates that the proposed sufficient conditions are reasonably strong. The global Mittag-Leffler stable behavior of the state variables y 1 ( t ) and y 2 ( t ) is demonstrated on Figure 1 for q = 0.4 , x 1 ( s ) = y 1 ( s ) = 2 , x 2 ( s ) = y 2 ( s ) = 3 , s [ 0.01 , 0 ] , y 1 ( 0 + ) = 2 , y 2 ( 0 + ) = 3 and t k = [ 0.03 , 0.07 , 0.11 , 0.15 , 0.19 , 0.23 ] .
Example 2.
Consider again the fractional impulsive Lotka-Volterra system (17) and (18) with a globally Mittag-Leffler stable equilibrium state y * = ( 2 7 , 1 5 ) . Let the region of the Mittag-Leffler stability be determined by
| | ϕ y * | | τ < k .
If λ and A are such that
k < λ < A ,
then for an initial function ϕ 0 = ( ϕ 01 , ϕ 02 ) T such that
k | | ϕ y * | | τ = sup s [ τ , 0 ] | | ϕ ( s ) y * | | = sup s [ τ , 0 ] | ϕ 01 ( s ) 2 7 | + | ϕ 02 ( s ) 1 5 | < λ
the state ( 2 7 , 1 5 ) will not be practically uniformly globally Mittag-Leffler stable. Note that in the particular example, condition 2 of Theorem 2 is not satisfied.
The example again shows that Mittag-Leffler stability is not enough for practical Mittag-Leffler stability.
Example 3.
Consider the two-species fractional cooperation Lotka-Volterra model
C D t q x 1 ( t ) = x 1 ( t ) 1 x 1 ( t τ 11 ) 1 + 5 x 2 ( t τ 12 ) 13 4 x 1 ( t τ 12 ) , C D t q x 2 ( t ) = x 2 ( t ) 1 x 2 ( t τ 22 ) 1 + 7 x 1 ( t τ 21 ) 4 3 x 2 ( t τ 21 ) ,
t 0 , 0 τ i j 1 , i , j = 1 , 2 .
We can easily check that the point ( 2 7 , 3 5 ) is an equilibrium for the model (19). However, since
2 m min 1 i 2 r i L a i M + M j = 1 j i n b j M + min 1 i 2 r i L c i L < max 1 i 2 r i M + max 1 i 2 j = 1 j i 2 b j M r i M M 2 a i L + m κ = 1 κ j 2 b κ L 2
we cannot make any conclusion about its global Mittag-Leffler stability.
For the corresponding fractional impulsive control Lotka-Volterra system
C D t q y 1 ( t ) = y 1 ( t ) 1 y 1 ( t τ 11 ) 1 + 5 y 2 ( t τ 12 ) 13 4 y 1 ( t τ 12 ) , t t k , C D t q y 2 ( t ) = y 2 ( t ) 1 y 2 ( t τ 22 ) 1 + 7 y 1 ( t τ 21 ) 4 3 y 2 ( t τ 21 ) , t t k , Δ y 1 ( t k ) = 2 3 y 1 ( t k ) 2 7 , k = 1 , 2 , , Δ y 2 ( t k ) = 3 7 y 2 ( t k ) 3 5 , k = 1 , 2 ,
with 0 < t 1 < t 2 < , lim k t k = consider the Lyapunov function
V ( t , y ) = i = 1 n | y i y i * | .
Let λ and A be given, and a ( u ) = u + 1 , u R + . Then all conditions of Theorem 3 are satisfied, and hence, the equilibrium ( 2 7 , 3 5 ) of (20) is practically uniformly globally Mittag-Leffler stable, which means that the impulsive controllers contribute to keep the population densities of both spaces between given bounds. For example, in the case q = 0.4 , t k = [ 0.03 , 0.07 , 0.11 , 0.15 , 0.19 , 0.23 ] , τ i j = 0.01 , i , j = 1 , 2 , x 1 ( s ) = y 1 ( s ) = 2 , x 2 ( s ) = y 2 ( s ) = 3 , s [ 0.01 , 0 ] , y 1 ( 0 + ) = 2 , y 2 ( 0 + ) = 3 , the practical uniform global Mittag-Leffler stable behavior of the model (20) is shown in Figure 2.
Example 4.
For the fractional impulsive Lotka-Volterra model (17), (18) consider the set
B = { z R + 2 : 0 z y * } ,
where y * = ( y 1 * , y 2 * ) = ( 2 7 , 1 5 ) .
Following the steps in the proof of Theorem 1 using the Lyapunov function
V ( t , y ) = | | y | | B
it is trivial to show that the set B is globally Mittag-Leffler stable with respect to (17), (18).

7. Conclusions

In this paper, we introduce a fractional-order Lotka-Volterra type cooperation model. Based on Lyapunov function theory we present an impulsive control strategy that can be applied as a Mittag-Leffler stability and synchronization mechanism. In addition, more general practical stability results and stability of sets results with respect to the impulsive control system are established. Examples are presented to demonstrate the obtained criteria and the relations between the investigated stability behaviors. It is shown that the proposed impulsive control technique provides a useful tool for the design of appropriate population models. It can be also beneficial to models with reaction-diffusion terms which are the subject of our future research.

Author Contributions

Conceptualization, F.S. and I.S.; methodology, R.T., F.S. and I.S.; formal analysis, R.T., F.S. and I.S.; investigation, R.T., F.S. and I.S.; visualization, R.T.; writing—original draft preparation, I.S. All authors have read and agreed to the published version of the manuscript.

Funding

F.S. and R.T. partially funded by National Institutes of Health BRAIN Theories NIMH-NIBIB R01EB026939.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Ahmad, S.; Stamova, I.M. (Eds.) Lotka-Volterra and Related Systems: Recent Developments in Population Dynamics, 1st ed.; Walter de Gruyter: Berlin, Germany, 2013; ISBN 978-3-11-026951-2. [Google Scholar]
  2. Gopalsamy, K. Stability and Oscillations in Delay Differential Equations of Population Dynamics, 1st ed.; Springer: Boston, MA, USA, 1992; ISBN 978-90-481-4119-7. [Google Scholar]
  3. Takeuchi, Y. Global Dynamical Properties of Lotka-Volterra Systems, 1st ed.; World Scientific: Singapore, 1996; ISBN 9810224710/9789810224714. [Google Scholar]
  4. Zhang, T.; Yang, L.; Xu, L. Stage-structured control on a class of predator-prey system in almost periodic environment. Intern. J. Control 2020, 93, 1442–1460. [Google Scholar] [CrossRef]
  5. Wei, F.; Wang, K. Asymptotically periodic solution of n-species cooperation system with time delay. Nonlinear Anal. Real World Appl. 2006, 7, 591–596. [Google Scholar] [CrossRef]
  6. May, R.M. Theoretical Ecology. Principles and Applications, 1st ed.; Saunders: Philadelphia, PA, USA, 1976; ISBN 0721662056/978-0721662053. [Google Scholar]
  7. Chen, F.; Liao, X.; Huang, Z. The dynamic behavior of N-species cooperation system with continuous time delays and feedback controls. Appl. Math. Comput. 2006, 181, 803–815. [Google Scholar] [CrossRef]
  8. Xu, R.; Chaplain, M.A.J.; Davidson, F.A. Global asymptotic stability of periodic solution for a cooperative system with time delays. Indian J. Pure Appl. Math. 2004, 35, 915–936. [Google Scholar]
  9. Yang, P.; Xu, R. Global asymptotic stability of periodic solution in n-species cooperative system with time delays. J. Biomath. 1998, 13, 841–846. [Google Scholar]
  10. Baez, J.C.; Pollard, B.S. Relative entropy in biological systems. Entropy 2016, 18, 46. [Google Scholar] [CrossRef]
  11. Rachdi, M.; Waku, D.; Hazgui, H.; Demongeot, J. Entropy as a robustness marker in genetic regulatory networks. Entropy 2020, 22, 260. [Google Scholar] [CrossRef] [Green Version]
  12. finnegan, A.; Song, J.S. Maximum entropy methods for extracting the learned features of deep neural networks. PLoS Comput. Biol. 2017, 13, e1005836. [Google Scholar] [CrossRef]
  13. Ahmad, S.; Stamova, I.M. Asymptotic stability of an n-dimensional impulsive competitive system. Nonlinear Anal. 2007, 8, 654–663. [Google Scholar] [CrossRef]
  14. Ballinger, G.; Liu, X. Permanence of population growth models with impulsive effects. Math. Comput. Model. 1997, 26, 59–72. [Google Scholar] [CrossRef]
  15. Dai, C.; Zhao, M.; Chen, L. Complex dynamic behavior of three-species ecological model with impulse perturbations and seasonal disturbances. Math. Comput. Simulation 2012, 84, 83–97. [Google Scholar] [CrossRef]
  16. Dong, L.; Chen, L.; Sun, L. Extinction and permanence of the predator-prey system with stocking of prey and harvesting of predator impulsively. Math. Methods Appl. Sci. 2006, 29, 415–425. [Google Scholar] [CrossRef]
  17. Dong, L.; Takeuchi, Y. Impulsive control of multiple Lotka-Volterra systems. Nonlinear Anal. Real World Appl. 2013, 14, 1144–1154. [Google Scholar] [CrossRef]
  18. Struk, O.O.; Tkachenko, V. On impulsive Lotka-Volterra systems with diffusion. Ukr. Math. J. 2002, 54, 629–646. [Google Scholar] [CrossRef]
  19. Tang, S.; Chen, L. The periodic predator-prey Lotka-Volterra model with impulsive effects. J. Mech. Med. Biol. 2002, 2, 267–296. [Google Scholar] [CrossRef]
  20. Yu, X.; Wang, Q.; Bai, Y. Permanence and almost periodic solutions for N-species nonautonomous Lotka-Volterra competitive systems with delays and impulsive perturbations on time scales. Complexity 2018, 2018, 2658745. [Google Scholar] [CrossRef] [Green Version]
  21. Stamova, I.M. Impulsive control for stability of n-species Lotka-Volterra cooperation models with finite delays. Appl. Math. Lett. 2010, 23, 1003–1007. [Google Scholar] [CrossRef]
  22. Haddad, W.M.; Chellaboina, V.S.; Nersesov, S.G. Impulsive and Hybrid Dynamical Systems, Stability, Dissipativity, and Control, 1st ed.; Princeton University Press: Princeton, NJ, USA, 2006; ISBN 9780691127156. [Google Scholar]
  23. Stamova, I.M.; Stamov, G.T. Applied Impulsive Mathematical Models, 1st ed.; Springer: Cham, Switzerland, 2016; ISBN 978-3-319-28060-8/978-3-319-28061-5. [Google Scholar]
  24. Yang, X.; Peng, D.; Lv, X.; Li, X. Recent progress in impulsive control systems. Math. Comput. Simul. 2019, 155, 244–268. [Google Scholar] [CrossRef]
  25. Abbas, S.; Banerjee, M.; Momani, S. Dynamical analysis of fractional-order modified logistic model. Comput. Math. Appl. 2011, 62, 1098–1104. [Google Scholar] [CrossRef] [Green Version]
  26. Mohyud-Din, S.T.; Ali, A.; Bin-Mohsin, B. On biological population model of fractional order. Int. J. Biomath. 2016, 9, 1650070. [Google Scholar] [CrossRef]
  27. Teka, W.; Marinov, T.; Santamaria, F. Neuronal spike timing adaptation described with a fractional leaky integrate-and-fire model. PLoS Comput. Biol. 2014, 10, e1003526. [Google Scholar] [CrossRef] [PubMed]
  28. Vázquez-Guerrero, P.; Gómez-Aguilar, J.F.; Santamaria, F.; Escobar-Jiméneza, R.F. Design of a high-gain observer for the synchronization of chimera states in neurons coupled with fractional dynamics. Physical A 2020, 539, 122896. [Google Scholar] [CrossRef]
  29. Agrawal, S.K.; Srivastava, M.; Das, S. Synchronization between fractional-order Ravinovich-Fabrikant and Lotka-Volterra systems. Nonlinear Dyn. 2012, 69, 2277–2288. [Google Scholar] [CrossRef]
  30. Das, S.K.; Gupta, P.K. A mathematical model on fractional Lotka-Volterra equations. J. Theor. Biol. 2011, 277, 1–6. [Google Scholar] [CrossRef]
  31. Gatabazi, P.; Mba, J.C.; Pindza, E. Fractional gray Lotka-Volterra models with application to cryptocurrencies adoption. Chaos 2019, 29, 073116. [Google Scholar] [CrossRef]
  32. Jun, Z.; Kim, C.-G. Positive solutions for a Lotka-Volterra prey-predator model with cross-diffusion of fractional type. Results Math. 2014, 65, 293–320. [Google Scholar] [CrossRef]
  33. Matlob, M.A.; Towers, I.N.; Jovanoski, Z.; Irwin, A.J. Memory and mutualism in species sustainability: A time-fractional Lotka-Volterra model with harvesting. arXiv 2020, arXiv:1904.12340v2. [Google Scholar]
  34. Wang, Y.; Liu, S. Fractal analysis and control of the fractional Lotka-Volterra model. Nonlinear Dyn. 2019, 95, 1457–1470. [Google Scholar] [CrossRef]
  35. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, 1st ed.; Elsevier: Amsterdam, The Netherlands, 2006; ISBN 0444518320/9780444518323. [Google Scholar]
  36. Podlubny, I. Fractional Differential Equations, 1st ed.; Academic Press: San Diego, CA, USA, 1999; ISBN 558840-2. [Google Scholar]
  37. Cattani, C.; Srivastava, H.M.; Yang, X.-J. (Eds.) Fractional Dynamics, 1st ed.; Springer: Berlin, Germany, 2015; ISBN 978-3-11-047209-7. [Google Scholar]
  38. Baleanu, D.; Diethelm, K.; Scalas, E.; Trujillo, J.J. Fractional Calculus. Models and Numerical Methods, 2nd ed.; World Scientific: Hackensack, NJ, USA, 2016; ISBN 9813140038/978-9813140035. [Google Scholar]
  39. Stamova, I.M.; Stamov, G.T. Functional and Impulsive Differential Equations of Fractional Order: Qualitative Analysis and Applications, 1st ed.; CRC Press; Taylor and Francis Group: Boca Raton, MA, USA, 2017; ISBN 9781498764834. [Google Scholar]
  40. Zhang, T.; Xiong, L. Periodic motion for impulsive fractional functional differential equations with piecewise Caputo derivative. Appl. Math. Lett. 2020, 101, 106072. [Google Scholar] [CrossRef]
  41. Debbouche, A.; Baleanu, D. Controllability of fractional evolution nonlocal impulsive quasilinear delay integro-differential systems. Comput. Math. Appl. 2011, 62, 1442–1450. [Google Scholar] [CrossRef] [Green Version]
  42. Stamov, G.; Stamova, I. Modelling and almost periodic processes in impulsive Lasota-Wazewska equations of fractional order with time-varying delays. Quaest. Math. 2017, 40, 1041–1057. [Google Scholar] [CrossRef]
  43. Stamov, G.; Stamova, I. Impulsive delayed Lasota-Wazewska fractional models: Global stability of integral manifolds. Mathematics 2019, 7, 1025. [Google Scholar] [CrossRef] [Green Version]
  44. Stamov, G.; Stamova, I. On almost periodic processes in impulsive fractional-order competitive systems. J. Math. Chem. 2018, 56, 583–596. [Google Scholar] [CrossRef]
  45. Li, Y.; Chen, Y.; Podlubny, I. Mittag-Leffler stability of fractional order nonlinear dynamic systems. Automatica 2009, 45, 1965–1969. [Google Scholar] [CrossRef]
  46. Li, H.; Kao, Y.G. Mittag-Leffler stability for a new coupled system of fractional-order differential equations with impulses. Appl. Math. Comput. 2019, 361, 22–31. [Google Scholar] [CrossRef]
  47. Stamova, I.M. Global Mittag-Leffler stability and synchronization of impulsive fractional-order neural networks with time-varying delays. Nonlinear Dyn. 2014, 77, 1251–1260. [Google Scholar] [CrossRef]
  48. Zhang, X.; Niu, P.; Ma, Y.; Wei, Y.; Li, G. Global Mittag-Leffler stability analysis of fractional-order impulsive neural networks with one-side Lipschitz condition. Neural Netw. 2017, 94, 67–75. [Google Scholar] [CrossRef]
  49. Li, H.-L.; Zhang, L.; Hu, C.; Jiang, H.; Cao, J. Global Mittag-Leffler synchronization of fractional-order delayed quaternion-valued neural networks: Direct quaternion approach. Appl. Math. Comput. 2020, 373, 125020. [Google Scholar] [CrossRef]
  50. Pratap, A.; Raja, R.; Rajchakit, G.; Cao, J.; Bagdasar, O. Mittag-Leffler state estimator design and synchronization analysis for fractional-order BAM neural networks with time delays. Intern. J. Adapt. Control Signal Process. 2019, 33, 855–874. [Google Scholar] [CrossRef]
  51. Zhou, H.-C.; Lv, C.; Guo, B.-Z.; Chen, Y.Q. Mittag-Leffler stabilization for an unstable time-fractional anomalous diffusion equation with boundary control matched disturbance. Intern. J. Robust Nonlinear Control 2019, 29, 4384–4401. [Google Scholar] [CrossRef]
  52. Ding, X.; Cao, J.; Zhao, X.; Alsaadi, F.E. Mittag-Leffler synchronization of delayed fractional-order bidirectional associative memory neural networks with discontinuous activations: State feedback control and impulsive control schemes. Proc. R. Soc. A Math. Phys. Eng. Sci. 2017, 473, 2204. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Pratap, A.; Raja, R.; Sowmiya, C.; Bagdasar, O.; Cao, J.; Rajchakit, G. Robust generalized Mittag-Leffler synchronization of fractional order neural networks with discontinuous activation and impulses. Neural Netw. 2018, 103, 128–141. [Google Scholar] [CrossRef]
  54. Stamova, I.; Stamov, G. Mittag-Leffler synchronization of fractional neural networks with time-varying delays and reaction-diffusion terms using impulsive and linear controllers. Neural Netw. 2017, 96, 22–32. [Google Scholar] [CrossRef] [PubMed]
  55. Ballinger, G.; Liu, X. Practical stability of impulsive delay differential equations and applications to control problems. In Optimization Methods and Applications. Applied Optimization; Yang, X., Teo, K.L., Caccetta, L., Eds.; Kluwer: Dordrecht, The Netherlands, 2001; Volume 52, pp. 3–21. [Google Scholar]
  56. Lakshmikantham, V.; Leela, S.; Martynyuk, A.A. Practical Stability of Nonlinear Systems; World Scientific: Teaneck, NJ, USA, 1990; ISBN 981-02-0351-9/981-02-0356-X. [Google Scholar]
  57. Yoshizawa, T. Asymptotic behavior of solutions of non-autonomous system near sets. J. Math. Kyoto Univ. 1962, 1, 303–323. [Google Scholar] [CrossRef]
  58. Çicek, M.; Yaker, C.; Gücen, M.B. Practical stability in terms of two measures for fractional order systems in Caputo’s sense with initial time difference. J. Frankl. Inst. 2014, 351, 732–742. [Google Scholar] [CrossRef]
  59. Bohner, M.; Stamova, I.; Stamov, G. Impulsive control functional differential systems of fractional order: Stability with respect to manifolds. Eur. Phys. J. Spec. Top. 2017, 226, 3591–3607. [Google Scholar] [CrossRef]
  60. Stamov, G.; Stamova, I.M.; Li, X.; Gospodinova, E. Practical stability with respect to h-manifolds for impulsive control functional differential equations with variable impulsive perturbations. Mathematics 2019, 7, 656. [Google Scholar] [CrossRef] [Green Version]
  61. Stamova, I.M.; Henderson, J. Practical stability analysis of fractional-order impulsive control systems. ISA Trans. 2016, 64, 77–85. [Google Scholar] [CrossRef]
  62. Hingu, D. Asymptotic stability of strongly uninvadable sets. Ann. Oper. Res. 2020, 287, 737–749. [Google Scholar] [CrossRef]
  63. Li, Y.; Sanfelice, R.G. finite time stability of sets for hybrid dynamical systems. Automatica 2019, 100, 200–211. [Google Scholar] [CrossRef] [Green Version]
  64. Mironchenko, A. Uniform weak attractivity and criteria for practical global asymptotic stability. Syst. Control Lett. 2017, 105, 92–99. [Google Scholar] [CrossRef] [Green Version]
  65. Stocker, C. Event-Based State-Feedback Control of Physically Interconnected Systems, 1st ed.; Logos: Berlin, Germany, 2014; ISBN 3832536825. [Google Scholar]
Figure 1. The global Mittag-Leffler stable behavior of the state variables y 1 ( t ) and y 2 ( t ) of the model (17) and (18).
Figure 1. The global Mittag-Leffler stable behavior of the state variables y 1 ( t ) and y 2 ( t ) of the model (17) and (18).
Entropy 22 00970 g001
Figure 2. The practical uniform global Mittag-Leffler stable behavior of the model (20). The state variables x 1 ( t ) and x 2 ( t ) of (19) and y 1 ( t ) and y 2 ( t ) of (20).
Figure 2. The practical uniform global Mittag-Leffler stable behavior of the model (20). The state variables x 1 ( t ) and x 2 ( t ) of (19) and y 1 ( t ) and y 2 ( t ) of (20).
Entropy 22 00970 g002

Share and Cite

MDPI and ACS Style

Tuladhar, R.; Santamaria, F.; Stamova, I. Fractional Lotka-Volterra-Type Cooperation Models: Impulsive Control on Their Stability Behavior. Entropy 2020, 22, 970. https://0-doi-org.brum.beds.ac.uk/10.3390/e22090970

AMA Style

Tuladhar R, Santamaria F, Stamova I. Fractional Lotka-Volterra-Type Cooperation Models: Impulsive Control on Their Stability Behavior. Entropy. 2020; 22(9):970. https://0-doi-org.brum.beds.ac.uk/10.3390/e22090970

Chicago/Turabian Style

Tuladhar, Rohisha, Fidel Santamaria, and Ivanka Stamova. 2020. "Fractional Lotka-Volterra-Type Cooperation Models: Impulsive Control on Their Stability Behavior" Entropy 22, no. 9: 970. https://0-doi-org.brum.beds.ac.uk/10.3390/e22090970

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