Next Article in Journal
Development of Preform for Simulation of Cold Forging Process of A V8 Engine Camshaft Free from Flash & Under-Filling
Next Article in Special Issue
Existence of Bounded Solutions to a Modified Version of the Bagley–Torvik Equation
Previous Article in Journal
Mathematical Analysis of an Autoimmune Diseases Model: Kinetic Approach
Previous Article in Special Issue
On Dynamic Systems in the Frame of Singular Function Dependent Kernel Fractional Derivatives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impulsive Delayed Lasota–Wazewska Fractional Models: Global Stability of Integral Manifolds

1
Department of Mathematics, Technical University of Sofia, 8800 Sliven, Bulgaria
2
Department of Mathematics, University of Texas at San Antonio, San Antonio, TX 78249, USA
*
Author to whom correspondence should be addressed.
Current address: Affiliation 2.
Mathematics 2019, 7(11), 1025; https://doi.org/10.3390/math7111025
Submission received: 9 October 2019 / Revised: 24 October 2019 / Accepted: 26 October 2019 / Published: 31 October 2019
(This article belongs to the Special Issue Stability Analysis of Fractional Systems)

Abstract

:
In this paper we deal with the problems of existence, boundedness and global stability of integral manifolds for impulsive Lasota–Wazewska equations of fractional order with time-varying delays and variable impulsive perturbations. The main results are obtained by employing the fractional Lyapunov method and comparison principle for impulsive fractional differential equations. With this research we generalize and improve some existing results on fractional-order models of cell production systems. These models and applied technique can be used in the investigation of integral manifolds in a wide range of biological and chemical processes.

1. Introduction

The intensive investigation of blood cell dynamics in numerous articles in the past 50 years leads to a huge forward movement in the growth of mathematical methods and models, numerical results, schemes to estimate parameters and prognosticate optimal treatments to particular diseases. The paper [1] offers a very completed overview of the main mathematical models related to blood formation, disorders and treatments.
A turning point in the establishment of new models and methods was the publication of the paper authored by Wazewska–Czyzewska and Lasota in 1976 [2]. In order to describe the survival of red blood cells in animals, Wazewska–Czyzewska and Lasota proposed in [2] the following delayed equation
x ˙ ( t ) = γ x ( t ) + β e α x ( t τ ) ,
where x ( t ) represents the number of red blood cells at time t, γ > 0 is the death probability for a red blood cell, α and β are positive constants related to the production of red blood cells per unit time and τ is the time delay between the production of immature red blood cells and their maturation for release in circulating blood stream.
The well known Lasota–Wazewska model seen in Equation (1) was extended and generalized by many authors, see, for example, refs. [3,4,5,6,7,8], including some recent publications on the topic [9,10,11,12,13]. This model can be also considered as one of the motivations to the development of the theory of delay differential equations, since delays are often considered in the hematopoiesis processes.
One of the directions in which the Lasota–Wazewska models have been extended and generalized is related to considering impulsive effects in the cells dynamic. Indeed, momentary (impulsive) changes at certain instants exists and can often affect the behavior of a real-world process. That is why the impulsive mathematical models has attracted a large amount of research interest and has become an emerging trend [14,15,16,17,18,19,20]. For some excellent results on impulsive Lasota–Wazewska models, we refer the reader to [21,22,23,24,25]. However, in all of the existing results, the authors considered fixed moments of impulsive effects. The aim of our research is to study some qualitative properties of a generalized impulsive Lasota–Wazewska model with variable impulsive perturbations. Indeed, considering impulsive effects at variable times is more general and close to reality [26,27,28].
On the other hand, fractional-order models are found to provide more advantages in describing memory effects and chaos [29,30]. The last two decades witnessed the rise of the development of the theory of fractional-order models and equations involved in fractional-order modeling. See, for example some recent publications [31,32,33,34,35] and the references therein. Some researchers used classical fractional-order models while new fractional differential systems, including impulsive fractional models started to be explored [36,37,38,39,40].
In addition, in relation to mathematical simulations in biology, chemistry and medicine, fractional calculus has been incorporated into some population dynamics models [41,42,43]. Since the survival of red blood cells is a deeply composite process, different methods have been developed throughout the years to get the most appropriate answers to some particular practice problems. To better reflect the dependence of Lasota–Wazewska-type models on their past history, an impulsive Lasota–Wazewska model of fractional order with time-varying delays has been introduced in [44] and some results related to existence and stability of almost periodic solutions of the model have been established. To the best of our knowledge no work on fractional Lasota–Wazewska systems was developed besides the one found in [44]. However, the paper [44] considered only fixed moments of impulsive perturbations.
Our aim in this paper is to extend the results in [44], and to consider variable impulsive perturbations. Furthermore, we will investigate the global stability behavior of an integral manifold related to the model. Indeed, the concept of integral manifolds is more general and includes, as particular cases, the global stability of zero solutions, equilibrium states, almost periodic solutions and etc. [45,46,47,48]. Despite the high importance of the concept of global stability of integral manifolds, the theory is not yet developed for fractional-order biological systems and this is another aim of the proposed research.
In the present paper, motivated by the above considerations, we focus on an impulsive Lasota–Wazewska model of fractional order with time-varying delays and variable impulsive perturbations. The remainder of this paper is organized as follows. In Section 2, we introduce the fractional Lasota–Wazewska model under consideration. Some notations and preliminaries are also given. In Section 3 first existence and boundedness results are proposed. Also, we obtain conditions for an integral manifold to be a global attractor of the survival model of red blood cells proposed by Wazewska and Lasota. The conditions are used in our global asymptotic stability analysis. In addition, the notion of global Mittag–Lerffler stability is defined and criteria for global Mittag–Leffler stability of integral manifolds with respect to the model are also proved. An example is presented in Section 4, to show the validity and efficiency of the obtained results. Finally, some ending remarks are stressed in Section 5.

2. Preliminaries

Let R = ( , ) , R + = [ 0 , ) , t 0 R + and let h = c o n s t > 0 . Here and in what follows we will use the Caputo fractional derivative of order q, 0 < q < 1 with the lower limit t 0 for a function l C 1 [ [ t 0 , b ] , R ] , b > t 0 ,
t 0 C D t q l ( t ) = 1 Γ ( 1 q ) t 0 t l ( σ ) ( t σ ) q d σ ,
where Γ as usually means the Gamma function.
Using the above fractional derivative, to consider long range memory in the survival of red blood cells, we introduce the following impulsive fractional-order Lasota–Wazewska model with time-varying delays and variable impulsive perturbations
t 0 C D t q x ( t ) = γ ( t ) x ( t ) + j = 1 m β j ( t ) e α j ( t ) x ( t s j ( t ) ) , t τ k ( x ( t ) ) , Δ x ( t ) = α k x ( t ) + ν k , t = τ k ( x ( t ) ) , k = 1 , 2 , ,
where:
(i) the model’s parameters α j ( t ) , β j ( t ) , γ ( t ) , s j ( t ) C [ R + , R + ] , 0 s j ( t ) h , j = 1 , 2 , , m , h = c o n s t > 0 , τ k : R + ( t 0 , ) , k = 1 , 2 , ;
(ii) Δ x ( t ) = x ( t + ) x ( t ) , α k , ν k R , k = 1 , 2 ,
The second equation of the impulsive system in Equation (2) is considered as a control or jump condition. The parameters α k and ν k , determine the controlled outputs x ( t + ) . In an impulsive control system of type seen in Equation (2), the functions Δ x ( t ) are considered as control forces at the variable times for t = τ k ( x ( t ) ) , k = 1 , 2 , For the basic concepts and theorems of such systems, we refer the reader to [14,20,21,22,23,24,25].
The model in Equation (2) is a generalization of many existing integer-order Lasota–Wazewska models with time-varying delays and impulsive perturbations. The consideration of time-varying delays s j ( t ) is motivated by the fact that “a long time delay will increase the length of time a population will spend in the neighborhood of an unstable steady state …” examined in [49]. The studies in [13,23,25] and some of the references therein also support the viewpoint that equations with time-varying delays provide a more realistic description for blood flow models.
Also, the fractional-order approach will better model the long-term memory phenomena and provides with a conceptually straightforward mathematical representation of rather complex processes. Indeed, a large number of empirical studies have investigated the long-term dependence of red blood cells and cell production systems in general. We refer the reader to [50,51,52] for some studies on long-range feedback effects in such systems. The very comprehensive analysis made in [44] shows that the population size of the red blood cells in an erythropoietic system has a very long memory. Therefore, the formulated fractional-order model in [44] of blood flow phenomenon is more realistic and takes into account the long-range hereditary properties of red blood cells populations. In addition, we consider variable impulsive perturbations in our research which are more useful and applicable in biological systems. Thus, our results will improve and generalize some known results obtained in [2,3,5,20,21,22,23,24,25,44]. For equations with non-fixed moments of impulsive perturbations, a number of difficulties related to the phenomena of ‘beating’ of the solutions, bifurcation, loss of the property of autonomy, etc. appeared. But the wide application of this type of equation requires their study.
The solutions x ( t ) of impulsive models with variable impulsive perturbations of the type in Equation (2) are piecewise continuous functions [14,20,26,28] that have first kind points of discontinuity at which they are continuous from the left, i.e., at the moments t l k when the integral curve of a solution x ( t ) meets the hypersurfaces
σ k = ( t , x ) [ t 0 , ) × Ω : t = τ k ( x ) ,
where it is continuous from the left, the following relations are satisfied:
x ( t l k ) = x ( t l k ) , x ( t l k + ) = x ( t l k ) + I l k ( x ( t l k ) ) .
The points t l 1 , t l 2 , ( t 0 < t l 1 < t l 2 < ) are the impulsive control instants. Let us note that, in general, k l k .
Let J R + be an interval. Define the following classes of piecewise continuous functions:
P C [ J , R n ] = σ : J R n : σ ( t ) is continuous everywhere except at some points t l k J at which σ ( t l k ) and σ ( t l k + ) exist and σ ( t l k ) = σ ( t l k ) } ;
P C = P C [ [ h , 0 ] , R + ] ;
P C B [ J , R + ] = { σ P C [ J , R + ] : σ ( t ) is bounded on J } .
Let φ 0 P C B [ [ h , 0 ] , R + ] . Denote by x ( t ) = x ( t ; t 0 , φ 0 ) , x R + the solution of the system in Equation (2) that satisfies the initial conditions:
x ( t ; t 0 , φ 0 ) = φ 0 ( t t 0 ) , t 0 h t t 0 , x ( t 0 + ; t 0 , φ 0 ) = φ 0 ( 0 ) .
In this paper, we will investigate such trajectories of solutions x ( t ) the motion along with which must be ensured by an appropriate choice of impulsive forces. To guarantee the existence, uniqueness and continuability of the solution x ( t ) = x ( t ; t 0 , φ 0 ) of the initial value problem (IVP) in Equation (2), Equation (3) on the interval [ t 0 , ) for φ 0 P C B [ [ h , 0 ] , R + ] and t 0 R + , as well as, the absence of the phenomenon ’beating’, we assume that:
1. τ 0 ( x ) t 0 for x R + , the functions τ k ( x ) are continuous and the following relations hold
t 0 < τ 1 ( x ) < τ 2 ( x ) < , τ k ( x ) as k
uniformly on x R + .
2. The functions α ( t ) , β j ( t ) , γ j ( t ) , s j ( t ) and τ k are continuous on R + .
For φ P C , we define | . | h as | φ | h = sup h ξ 0 | φ ( ξ ) | . In the case h = we have | φ | h = | φ | = sup ξ ( , 0 ] | φ ( ξ ) | .
We shall use the following definition for integral manifolds connected with Equation (2) [20,48].
Definition 1.
An arbitrary manifold M in the extended phase space [ t 0 h , ) × R + of (3) is called an integral manifold , if for any solution x ( t ) = x ( t ; t 0 , φ 0 ) , ( t , φ 0 ( t t 0 ) ) M , t [ t 0 h , t 0 ] implies ( t , x ( t ) ) M , t t 0 .
For a manifold M [ t 0 h , ) × R + we introduce the following sets and distances:
M ( t ) is the set of all x R + such that ( t , x ) M for t [ t 0 , ) ;
M 0 ( t ) is the set of all x R + such that ( t , x ) M for t [ t 0 h , t 0 ] ;
d ( x , M ( t ) ) = inf y M ( t ) | x y | is the distance between x R + and M ( t ) ;
M ( t , ε ) = { x R + : d ( x , M ( t ) ) < ε } ( ε > 0 ) is an ε -neighborhood of M ( t ) ;
d 0 ( φ , M 0 ( t ) ) = sup t [ t 0 h , t 0 ] d ( φ ( t t 0 ) , M 0 ( t ) ) is the distance between a function φ P C [ [ h , 0 ] , R + ] and M 0 ( t ) ;
M 0 ( t , ε ) = { φ P C [ [ h , 0 ] , R + ] : d 0 ( φ , M 0 ( t ) ) < ε } is an ε -neighborhood of M 0 ( t ) ;
S a ¯ ( P C 0 ) = { φ P C : | φ | h a , a = c o n s t > 0 } .
In order to realize our investigations we will need the following assumptions:
A1. The set M ( t ) is nonempty for t [ t 0 , ) .
A2. The set M 0 ( t ) is nonempty for t [ t 0 h , t 0 ] .
A3. The distance d ( x , M ( t ) ) is Lipschitz with respect to t on any compact subset F of [ t 0 , ) × R + .
We also will use the following definitions.
Definition 2.
The integral manifold M of system in Equation (2) is said to be:
(a) equi-bounded, if for any initial point t 0 R + and any positive constants η > 0 and a > 0 there exists a positive number b = b ( t 0 , η , a ) > 0 such that for any initial function φ S a ¯ ( P C 0 ) M 0 ( t , η ) the solution x ( t ; t 0 , φ ) M ( t , b ) for all t t 0 ;
(b) t (or a)- uniformly bounded, if the number b from (a) is independent of t 0 ( o r o f a ) ;
(c) uniformly bounded, if the number b from (a) depends only on η.
Definition 3.
The integral manifold M is said to be:
(a) stable with respect to the system in Equation (2), if for any initial point t 0 R + and any positive constants a > 0 and ε > 0 there exists a positive δ = δ ( t 0 , a , ε ) > 0 such that for any initial function such that φ S a ¯ ( P C 0 ) M 0 ( t , δ ) the solution x ( t ; t 0 , φ ) M ( t , ε ) for all t t 0 ;
(b) uniformly stable with respect to the system in Equation (2) if the number δ from (a) depends only on ε;
(c) uniformly globally attractive with respect to the system in Equation (2), if for any η > 0 and ε > 0 there exists a positive number T = T ( η , ε ) > 0 such that for any initial point t 0 R + , any constant a > 0 and any initial function such that φ S a ¯ ( P C 0 ) M 0 ( t , η ) the solution x ( t ; t 0 , φ ) M ( t , ε ) for all t t 0 + T ;
(d) uniformly globally asymptotically stable with respect to the system in Equation (2) if M is a uniformly stable, uniformly globally attractive and uniformly bounded with respect to the system in Equation (2).
Remark 1.
Point (d) of Definition 3 can be reduced to the following specific stability notions:
1. Lyapunov uniform global asymptotic stability of the zero solution of Equation (2), if
M = [ t 0 h , ) × { x 0 } .
2. Lyapunov uniform global asymptotic stability of a non-zero equilibrium state x * = x * ( t ) of Equation (2), if
M = [ t 0 h , ) × { x R + : x x * } .
3. Uniform global asymptotic stability of conditionally integral manifold B with respect to a integral manifold A, where A B R + , if M ( t ) = B for t t 0 and M 0 ( t ) = A for t [ t 0 h , t 0 ] .
Next, motivated by [38,53], for an integral manifold M with respect to Equation (2) we will define a generalization of the global exponential stability notion to the fractional-order case, called global Mittag–Leffler stability. To this end we will need the Mittag–Leffler function [29,30] defined as
E q ( z ) = κ = 0 z κ Γ ( q κ + 1 ) , q > 0 .
Definition 4.
The integral manifold M is said to be globally Mittag–Leffler stable with respect to Equation (2), if there exist constants μ > 0 and d ˜ > 0 such that
x ( t ; t 0 , φ 0 ) M t , { m [ φ 0 ] E q ( μ ( t t 0 ) q ) } d ˜ , t t 0 ,
where E q is the corresponding Mittag–Leffler function, φ 0 P C B [ [ h , 0 ] , R + ] , m ( 0 ) = 0 , m ( φ ) 0 and m ( φ ) is Lipschitz with respect to φ P C B [ [ h , 0 ] , R + ] .
Let τ 0 ( x ) = t 0 , x R + . Next, we need the following sets
G k = ( t , x ) [ t 0 , ) × R + : τ k 1 ( x ) < t < τ k ( x ) , k = 1 , 2 , , G = k = 1 G k
and in the future considerations, we will adopt the Lyapunov–Razumikhin approach. That is why in the proofs of our main results we will use the class of piecewise continuous auxiliary functions given as:
V M = { V : [ t 0 , ) × R + R + : V C [ G , R + ] , t [ t 0 , ) ,
V is locally Lipschitz with respect to its second argument on each of the sets G k ; V ( t , x ) = 0 for ( t , x ) M , t t 0 , V ( t , x ) > 0 for ( t , x ) [ t 0 , ) × R + \ M , for ( t 0 * , x 0 * ) σ k , V ( t 0 * , x 0 * ) and V ( t 0 * + , x 0 * ) exist and V ( t 0 * , x 0 * ) = V ( t 0 * , x 0 * ) } .
Lets, for simplicity, denote
X ( t , x ) = γ ( t ) x ( t ) + j = 1 m β j ( t ) e α j ( t ) x ( t s j ( t ) ) , t τ k ( x ) , x R + .
Using the above notation, for a function V V M we define the following fractional-order derivative.
Definition 5.
Let t ( τ k 1 ( x ) , τ k ( x ) ) , k = 1 , 2 , , x R + and φ PC . For V V M the upper right-hand derivative of V in Caputo’s sense of order q, 0 < q < 1 with respect to the system in Equation (2) is defined by
C D + q V ( t , φ ( 0 ) ) = lim χ 0 + sup 1 χ q V ( t , φ ( 0 ) ) V ( t χ , φ ( 0 ) χ q X ( t , φ ) ) .
The following class of weight functions will also be useful in the proofs of our main theorems:
K = { w C [ R + , R + ] : w ( r ) is strictly increasing and w ( 0 ) = 0 , w ( r ) , r } .
Let t 1 , t 2 , ( t 0 < t 1 < t 2 < ) be the impulsive control instants at which the integral curve ( t , x ( t ; t 0 , φ 0 ) ) of the IVP Equation (2), Equation (3) meets the hypersurfaces σ k , k = 1 , 2 , , i.e., each of the points t k is a solution of some of the equations t = τ k ( x ( t ) ) , k = 1 , 2 ,
In the next section we shall use the following Lemma from [38]. Similar comparison results can be found in [44] and the references therein.
Lemma 1.
Assume that the function V V M is such that for t [ t 0 , ) , φ PC ,
V ( t + , φ ( 0 ) + Δ φ ) V ( t , φ ( 0 ) ) , t = t k ,
and for a continuous function λ : [ t 0 , ) R the inequality
C D + q V ( t , φ ( 0 ) ) λ ( t ) V ( t , φ ( 0 ) ) , t t k , k = 1 , 2 ,
is valid whenever V ( t + ξ , φ ( ξ ) ) V ( t , φ ( 0 ) ) for h ξ 0 .
Then s u p h ξ 0 V ( t 0 + ξ , φ 0 ( ξ ) ) V ( t , φ ( 0 ) ) implies
V ( t , x ( t ; t 0 , φ 0 ) ) s u p h ξ 0 V ( t 0 + , φ 0 ( ξ ) ) E q ( λ ( t ) ( t t 0 ) q ) , t [ t 0 , ) .
In the case when λ ( t ) = 0 for t [ t 0 , ) the following corollary follows directly from Lemma 1.
Corollary 1.
Assume that the function V V M is such that for t [ t 0 , ) , φ PC ,
V ( t + , φ ( 0 ) + Δ φ ) V ( t , φ ( 0 ) ) , t = t k ,
and the inequality
C D + q V ( t , φ ( 0 ) ) 0 , t t k , k = 1 , 2 ,
is valid whenever V ( t + ξ , φ ( ξ ) ) V ( t , φ ( 0 ) ) for h ξ 0 .
Then
V ( t , x ( t ; t 0 , φ 0 ) ) s u p h ξ 0 V ( t 0 + , φ 0 ( ξ ) ) , t [ t 0 , ) .

3. Main Results

In this section, for a bounded continuous function f defined on R + , we set
f ¯ = sup t R + f ( t ) , f ̲ = inf t R + f ( t ) .

3.1. Existence and Boundedness Results

Theorem 1.
Assume that:
1. Assumptions A1–A3 are satisfied.
2. M is a manifold in the extended phase space of the system in Equation (2).
3. The functions α j ( t ) , β j ( t ) , γ ( t ) , s j ( t ) are bounded on R + , s j ( t ) < min { t , h } , j = 1 , 2 , , m , t [ t 0 , ) .
4. The sequences of constants { α k } are such that
1 < α k 0 , k = 1 , 2 , .
5. The model’s parameters are such that for t τ k ( x ) , k = 1 , 2 , , x R +
j = 1 m β j α j ¯ γ ̲ < 1 .
Then M is an integral manifold of Equation (2).
Proof. 
Let t 0 R + , φ 0 P C . Let x ( t ; t 0 , φ 0 ) be the solution of the IVP Equation (2), Equation (3) and ( t , φ 0 ( t t 0 ) ) M for t [ t 0 h , t 0 ] . We will prove that M is an integral manifold of Equation (2). If we suppose that this is not true, then there exists a t , t > t 0 such that ( t , x ( t ; t 0 , φ 0 ) ) M for t 0 < t t and ( t , x ( t ; t 0 , φ 0 ) ) M for t > t .
Consider the Lyapunov function V V M defined as
V ( t , x ) = d ( x , M ( t ) ) ,
where t [ t 0 h , ) .
It is easy to see that there exists t , t > t such that ( t , x ( t ; t 0 , φ 0 ) ) M and V ( t , x ( t ; t 0 , φ 0 ) ) > 0 .
In the case when t = τ k ( x ) , x R + , by condition 4, we get
V ( t + , φ ( 0 ) + Δ φ ) = d ( φ ( 0 ) + Δ φ , M ( t + ) )
inf φ 1 ( 0 ) M ( t ) | φ ( 0 ) + Δ φ φ 1 ( 0 ) Δ φ 1 ( 0 ) | = inf φ 1 ( 0 ) M ( t ) | φ ( 0 ) + α k φ ( 0 ) φ 1 ( 0 ) α k φ 1 ( 0 ) |
( 1 + α k ) inf φ 1 ( 0 ) M ( t ) | φ ( 0 ) φ 1 ( 0 ) | V ( t , φ ( 0 ) ) .
Let t t 0 , t τ k ( x ) , x R + . We can see that
t 0 C D t q | x ( t ) x * ( t ) | = s i g n ( x ( t ) x * ( t ) ) t 0 C D t q ( x ( t ) x * ( t ) ) .
Then for t t 0 and t ( τ k 1 ( x ) , τ k ( x ) ) and for the derivative C D + q V ( t , φ ( 0 ) ) along the solutions of the system in Equation (2) we have
C D + q V ( t , φ ( 0 ) ) γ ( t ) | φ ( 0 ) φ 1 ( 0 ) | + j = 1 m β j ( t ) | e α j ( t ) φ ( s j ( 0 ) ) e α j ( t ) φ 1 ( s j ( 0 ) ) |
γ ( t ) | φ ( 0 ) φ 1 ( 0 ) | + j = 1 m β j ( t ) α j ( t ) | φ ( s j ( 0 ) ) φ 1 ( s j ( 0 ) ) | , φ 1 ( 0 ) M ( t ) .
From the above estimate, we have
C D + q V ( t , φ ( 0 ) ) γ ̲ V ( t , φ ( 0 ) ) + j = 1 m β j α j ¯ sup h ξ 0 V ( t + ξ , φ ( ξ ) )
and for any solution x ( t ) that satisfies the Razumikhin condition
V ( t + ξ , φ ( ξ ) ) V ( t , φ ( 0 ) ) , h ξ 0
by virtue of condition 5 of Theorem 1, it follows that
C D + q V ( t , φ ( 0 ) ) 0 , t τ k ( x ) , t > t 0 .
Then using Equations (5) and (6) and Corollary 1, we get
V ( t , x ( t ) ) sup h ξ 0 V ( t 0 + , φ 0 ( ξ ) ) , t [ t 0 , ) .
From the last inequality we get V ( t , x ( t ; t 0 , φ 0 ) ) sup h ξ 0 V ( t 0 + , φ 0 ( ξ ) ) = 0 , t [ t 0 , ) which contradicts the fact that V ( t , x ( t ; t 0 , φ 0 ) ) > 0 and therefore, it shows that M is an integral manifold of Equation (2).
The proof of Theorem 1 is complete. □
Now, we will proof some boundedness results.
Theorem 2.
Let conditions of Theorem 1 are satisfied. Then the integral manifold M of Equation (2) is uniformly bounded.
Proof. 
Let a > 0 , η > 0 and t 0 R + . Consider again the Lyapunov function seen in Equation (4). It follows from the choice of the function V that there exist functions w 1 , w 2 K that satisfy
w 1 ( d ( x , M ( t ) ) ) V ( t , x ) w 2 ( d ( x , M ( t ) ) ) , ( t , x ) [ t 0 , ) × R + .
Now, we suppose that x ( t ) = x ( t ; t 0 , φ 0 ) is a solution of problem Equations (2) and (3). From the definition of the weight functions w 1 , w 2 K , it follows that we can choose the number b = b ( η ) > 0 so that w 2 ( η ) < w 1 ( b ) .
Let φ 0 S a ¯ ( P C 0 ) M 0 ( t , η ) . Since all conditions of Theorem 1 are satisfied, we get Equations (5) and (7). Then by Equations (5), (7) and (8), we have
w 1 ( d ( x ( t ) , M ( t ) ) ) V ( t , x ( t ) ) s u p h ξ 0 V ( t 0 + , φ 0 ( ξ ) )
w 2 ( d 0 ( φ 0 , M 0 ( t ) ) ) w 2 ( η ) < w 1 ( b ) , t t 0 .
Therefore, x ( t ; t 0 , φ 0 ) M ( t , b ) for t [ t 0 , ) and the theorem is proved. □
The proofs of the next two theorems is analogous to the proof of Theorem 2. They present sufficient conditions for t (respectively for a)-uniform boundedness of the integral manifold M with respect to Equation (2).
Theorem 3.
Let the conditions of Theorem 2 hold and for the Lyapunov function in Equation (4) there exists w 2 ( . , s ) K such that for each fixed s 0
V ( t , x ) w 2 ( d ( x , M ( t ) ) , x ) f o r ( t , x ) [ t 0 , ) × R + .
Then the integral manifold M of Equation (2) is t-uniformly bounded.
Theorem 4.
Let the conditions of Theorem 2 hold and for the Lyapunov function in Equation (4) there exists w 2 ( t , . ) K such that for each t [ t 0 , )
V ( t , x ) w 2 ( t , d ( x , M ( t ) ) ) f o r ( t , x ) [ t 0 , ) × R + .
Then the integral manifold M of Equation (2) is α- uniformly bounded.
Remark 2.
Theorems 1–4 offer efficient criteria for existence and boundedness of integral manifolds related to the fractional-order Lasota–Wazewska model in Equation (2). It is well known that boundedness is an important property that plays a significant role in the existence of permanent, periodic and almost-periodic solutions of different systems [20,38,54,55,56,57]. Since integral manifolds are sets of solutions of the system (see, [45,46,47,48]), our results generalize and complement the existing qualitative results for separate solutions of Equation (2), biological models [4,10,11,12,19,20,21,22,23,24,25,38,44].

3.2. Global Asymptotic Stability

In this section, we shall use the measurable function λ : [ t 0 , ) R + . It is integrally positive if
J ^ λ ( t ) d t =
whenever J ^ = k = 1 [ a k , b k ] , a k < b k < a k + 1 and b k a k θ > 0 , k = 1 , 2 ,
Theorem 5.
Assume that:
1. Conditions 1–4 of Theorem 1 are satisfied.
2. There exists an integrally positive function λ = λ ( t ) : R + R + such that for t τ k ( x ) , k = 1 , 2 , , x R +
γ ̲ j = 1 m β j α j ¯ > λ ( t ) > 0 .
3. η Γ ( q ) 0 ( t s ) q 1 λ ( s ) d s = for sufficiently small values of η > 0 .
Then the integral manifold M is uniformly globally asymptotically stable with respect to Equation (2).
Proof. 
Let t 0 R + . The fact that M is a uniformly bounded integral manifold of Equation (2) follows from Theorem 2.
Consider the Lyapunov function V V M defined by Equation (4).
We will first prove the uniform stability of the integral manifold M with respect to Equation (2). Let ε > 0 and take δ = δ ( ε ) > 0 so that δ < ε .
Let a > 0 be arbitrary, φ 0 S a ¯ ( P C 0 ) M 0 ( t , δ ) and x ( t ) = x ( t ; t 0 , φ 0 ) be the solution of Equation (2) through ( t 0 , φ 0 ) .
By condition 2 of Theorem 5, using similar arguments as in the proof of Theorem 1, for the case t t 0 , t τ k ( x ) , x R + , we have
C D + q V ( t , φ ( 0 ) ) λ ( t ) V ( t , φ ( 0 ) ) , t τ k ( x ) , t > t 0 .
From Equations (5) and (9) and Corollary 1, we get
d ( x ( t ; t 0 , φ 0 ) , M ( t ) ) = V ( t , x ( t ) ) s u p h ξ 0 V ( t 0 + , φ 0 ( ξ ) )
= s u p h ξ 0 d ( φ 0 ( ξ ) , M 0 ( t 0 + ) ) d 0 ( φ 0 , M 0 ( t ) ) < δ < ε , t t 0 ,
hence the integral manifold M is uniformly stable with respect to Equation (2).
Next, for the given ε > 0 and η < ε , in view of condition 3 of Theorem 5, we can choose the number T = T ( η , ε ) > 0 so that
ε Γ ( q ) t 0 t 0 + T ( t 0 + T s ) q 1 λ ( s ) d s > η .
Let a > 0 be arbitrary, φ 0 S a ¯ ( P C 0 ) M 0 ( t , η ) and if we assume that for any t [ t 0 , t 0 + T ] we have
d ( x ( t ; t 0 , φ 0 ) , M ( t ) ) ε
then by Equations (5), (9) and (11) it follows that
V ( t , x ( t ; t 0 , φ 0 ) ) V ( t 0 + , φ 0 )
1 Γ ( q ) t 0 t 0 + T ( t 0 + T s ) q 1 λ ( s ) d ( x ( s ; t 0 , φ 0 ) , M ( s ) ) d s
V ( t 0 + , φ 0 ) ε Γ ( q ) t 0 t 0 + T ( t 0 + T s ) q 1 λ ( s ) d s , t [ t 0 , t 0 + μ ] .
From the above estimate, Equations (5) and (10) for t = t 0 + T we obtain
V ( t 0 + T , x ( t 0 + T ; t 0 , φ 0 ) ) η ε Γ ( q ) t 0 t 0 + T ( t 0 + T s ) q 1 λ ( s ) d s < 0 ,
which contradicts the fact that V V M . Hence, there exists a t * [ t 0 , t 0 + T ] , such that
d ( x ( t * ; t 0 , φ 0 ) , M ( t * ) ) < ε .
From the lack of increase of the function V along the solution x ( t ) = x ( t ; t 0 , φ ) it follows that for t t * (hence for any t t 0 + T as well) we have
d ( x ( t ; t 0 , φ 0 ) , M ( t ) ) ) = V ( t , x ( t ; t 0 , φ 0 ) ) s u p t * h t t * V ( t * + , x ( t t * ) )
= s u p t * h t t * d ( x ( t t * ) , M ( t t * ) ) < η < ε .
Therefore, x ( t ) M ( t , ϵ ) for t t 0 + T . Hence, the integral manifold M is uniformly globally attractive with respect to Equation (2).
The proof is complete. □
The proof of the next theorem is similar to the proof of Theorem 5 and we will omit it here.
Theorem 6.
Let the conditions 1 and 2 of Theorem 5 hold and for the Lyapunov function in Equation (4) there exists a function w K and integrally positive functions ζ , λ : [ t 0 , ) R + such that
ζ ( t ) w ( d ( x ( t ) , M ( t ) ) ) λ 1 ( t ) V ( t , x ) , t τ k ( x ) , k = 1 , 2 ,
for t [ t 0 , ) .
Then the integral manifold M is uniformly globally asymptotically stable with respect to Equation (2).
Remark 3.
Theorems 5 and 6 offer new sufficient conditions on the model’s parameters of the impulsive fractional Lasota–Wazewska Equation (2) that yield the global asymptotic stability of an integral manifold of states. Also, the obtained results extend the existing global attractivity and stability results for impulsive Lasota–Wazewska models [21,22,23,24,25] to the fractional-order case which better reflects the multiple history dependent phenomena in the red blood cells dynamics. In fact, stability properties are the most investigated qualitative properties of biological systems concerned with the study of their behavior over a finite or infinite interval of time. Stability is also a very essential and crucial issue in control systems (such as impulsive control systems [20]). Hence, formulating rigorous stability results is a research problem of theoretical and practical significance.
Remark 4.
With this research we generalize the results in [44], considering variable impulsive perturbations and integral manifolds. For example, if M = [ t 0 h , ) × { x R + : x x * } , then two solutions x ( t ) and x * ( t ) M ( t ) may have different impulsive moments t k and t k * which is more general than the case of fixed moments of impulsive perturbations, considered in [44]. Since taking into account variable impulsive perturbations in impulsive control systems is more natural and realistic, our results have great opportunities for applications. If the impulses are realized at fixed times, and the integral manifold M = [ t 0 h , ) × { x R + : x w } , where w is the unique almost periodic solution of the Lasota–Wazewska model, then the results in [44] follow as corollaries of our results.

3.3. Global Mittag–Leffler Stability

In this Section we will present our Mittag–Leffler stability results for the integral manifold M.
Theorem 7.
Assume that conditions 1 and 2 of Theorem 5 are met.
Then the integral manifold M is globally Mittag–Leffler stable with respect to Equation (2).
Proof. 
Let t 0 R + and φ 0 P C B [ [ h , 0 ] , R + ] and x ( t ) = x ( t ; t 0 , φ 0 ) be the solution of the IVP Equations (2) and (3). From conditions 1 and 2 of Theorem 5 we get Equations (5) and (9).
Consider again the Lyapunov function Equation (4). For it there exists a constant ν > 0 such that
V ( t , x ) < ν | x | , t [ t 0 , ) , x R + .
Then from the above estimate, Equations (5) and (9), using Lemma 1, we obtain
d ( x ( t ; t 0 , φ 0 ) , M ( t ) ) = V ( t , x ( t ) ) sup h ξ 0 V ( t 0 + , φ 0 ( ξ ) ) E q ( λ ( t ) ( t t 0 ) q ) < ν | φ 0 | h E q ( λ ( t ) ( t t 0 ) q ) , t t 0 .
Let m [ φ 0 ] = ν | φ 0 | h . Then
d ( x ( t ; t 0 , φ 0 ) , M ( t ) ) < m [ φ 0 ] E q ( λ ( t ) ( t t 0 ) q ) , t t 0 .
Since in the above estimate m 0 and m = 0 is true only if φ 0 = 0 , this proves that the integral manifold M is globally Mittag–Leffler stable with respect to Equation (2). □
Remark 5.
If in Theorem 7 we consider q = 1 , then we will have
d ( x ( t ; t 0 , φ 0 ) , M ( t ) ) ν | φ 0 | h exp ( λ ( t ) ( t t 0 ) ) , t t 0 ,
which implies the global exponential stability of the integral manifold M with respect to Equation (2). This shows that the notion of Mittag–Leffler stability for fractional-order differential equations is an extension of the notion of exponential stability for integer-order systems [20,53]. Therefore, with this research we generalize and improve the existing global exponential stability results for survival of red blood cells models [3,5,8,10,13,20,21,25] to the fractional-order case. It is also well known that the exponential stability guarantees the fast convergence rate which is preferable and desirable for mathematical models in biology, chemistry and medicine [20,24].

4. An Example

In this section, we present an example to illustrate the obtained results.
Consider the following impulsive Lasota–Wazewska fractional-order model with time-varying delays and variable impulsive perturbations:
0 C D t q x ( t ) = ( 11 + sin t ) x ( t ) + ( 3.5 + 0.4 sin t ) e ( 0.6 + 0.4 cos t ) x ( t 2 cos t ) + ( 3.1 + 0.2 cos t ) e ( 0.3 + 0.2 sin t ) x ( t 2 sin t ) , t τ k ( x ) , x ( t + ) = 1 k x ( t ) + ν k , t = τ k ( x ) , k = 1 , 2 , ,
where m = 2 , t 0 , 0 < q < 1 , γ ( t ) = 11 + sin t , β 1 ( t ) = 3.5 + 0.4 sin t , β 2 ( t ) = 3.1 + 0.2 cos t , α 1 ( t ) = 0.6 + 0.4 cos t , α 2 ( t ) = 0.3 + 0.2 sin t , s 1 ( t ) = 2 + cos t , s 2 ( t ) = 2 + sin t , ν k R , τ k ( x ) = | x | + k , k = 1 , 2 ,
First, we have that the functions τ k ( x ) are continuous on R + and satisfy
τ 1 ( x ) < τ 2 ( x ) < , τ k ( x ) as k .
Consider the manifold
M = [ 3 , ) × { x R + : M ̲ x M ¯ } ,
where M ̲ , M ¯ R + are two constant solutions of Equation (12).
It is easy to check that for the manifold Equation (13) assumptions A1–A3 are satisfied.
Now, we have that α ¯ 1 = 1 , α ¯ 2 = 0.5 , β ¯ 1 = 3.9 , β ¯ 2 = 3.3 , γ ̲ = 10 , and
j = 1 m β j α j ¯ γ ̲ = 5.55 10 < 1 .
i.e., condition 5 of Theorem 1 is satisfied.
Moreover,
1 < α k = 1 k k 0 ,
for k = 1 , 2 ,
Therefore, according to Theorem 1, M is an integral manifold for Equation (12), and according to Theorem 2, M is uniformly bounded.
Also, if there exists an integrally positive function λ = λ ( t ) : R + R + such that for t τ k ( x ) , k = 1 , 2 , , x R +
γ ̲ j = 1 m β j α j ¯ = 4 . 45 > λ ( t ) > 0 ,
then according to Theorem 7 the integral manifold M is globally Mittag–Leffler stable with respect to Equation (12).
In addition, if the function λ is such that η Γ ( q ) 0 ( t s ) q 1 λ ( s ) d s = for each sufficiently small value of η > 0 , the according to Theorem 5 the integral manifold M is uniformly globally asymptotically stable with respect to Equation (12).
Remark 6.
With our example we illustrated the established theoretical results. Since the notion of stability of manifolds includes as particular cases stability of zero solutions, equilibrium states, almost periodic solutions, etc., our results have universal applicability and can be easily expanded in the study of many other fractional biochemical reactions processes.

5. Conclusions

A large number of empirical studies have investigated the long-term dependence of red blood cells and cell production systems in general. To better model the long-term memory phenomena in the survival of red blood cells models, in this paper we extend the existing impulsive Lasota–Wazewska models with time-varying delays to the fractional-order case. In addition, we consider variable impulsive perturbations and integral manifolds in our analysis. Existence, boundedness, uniform global asymptotic stability and global Mittag–Leffler stability results are established. Since the notion of stability of integral manifolds is much more general than that of trivial solutions, equilibrium points, periodic and almost-periodic solutions, etc., our results generalized many existing boundedness and stability criteria. The generalized concept and the results obtained can be applied to study other types of impulsive control fractional biochemical systems.

Author Contributions

Conceptualization, G.S. and I.S.; methodology, G.S. and I.S.; formal analysis, G.S. and I.S.; investigation, G.S. and I.S.; writing—original draft preparation, I.S.

Funding

This research was funded in part by the European Regional Development Fund through the Operational Programme “Science and Education for Smart Growth” under contract UNITe No. BG05M2OP001–1.001–0004 (2018–2023).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Pujo-Menjouet, L. Blood cell dynamics: Half of a century of modelling. Math. Model. Nat. Phenom. 2016, 11, 92–115. [Google Scholar] [CrossRef]
  2. Wazewska–Czyzewska, M.; Lasota, A. Mathematical problems of the dynamics of red blood cell population. Math. Stos. 1976, 6, 23–40. (In Polish) [Google Scholar]
  3. Gyori, I.; Ladas, G. Oscillation Theory of Delay Differential Equations. With Applications, 1st ed.; Oxford University Press: London, UK, 1991; ISBN 0750301422, 9780750301428. [Google Scholar]
  4. Kiskinov, H.; Zahariev, A.; Zlatev, S. Permanence of the positive solutions of the generalised haematopoiesis Mackey-Glass model. C. R. Acad. Bulgare Sci. 2014, 67, 745–752. [Google Scholar]
  5. Kulenovic, M.R.S.; Ladas, G.; Sficas, Y.G. Global attractivity in population dynamics. Comput. Math. Appl. 1989, 18, 925–928. [Google Scholar] [CrossRef] [Green Version]
  6. Lasota, A.; Loskot, K.; Mackey, M.C. Stability properties of proliferatively coupled cell replication models. Acta Biotheor. 1991, 39, 1–14. [Google Scholar] [CrossRef]
  7. Lasota, A.; Mackey, M.C. Cell division and the stability of cellular replication. J. Math. Biol. 1999, 38, 241–261. [Google Scholar] [CrossRef]
  8. Liu, G.; Zhao, A.; Yan, J. Existence and global attractivity of unique positive periodic solution for a Lasota–Wazewska model. Nonlinear Anal. 2006, 64, 1737–1746. [Google Scholar] [CrossRef]
  9. Chérif, F.; Miraoui, M. New results for a Lasota–Wazewska model. Int. J. Biomath. 2019, 12, 1950019. [Google Scholar] [CrossRef]
  10. Duan, L.; Huang, L.; Chen, Y. Global exponential stability of periodic solutions to a delay Lasota–Wazewska model with discontinuous harvesting. Proc. Am. Math. Soc. 2016, 144, 561–573. [Google Scholar] [CrossRef]
  11. Rihani, S.; Kessab, A.; Chérif, F. Pseudo-almost periodic solutions for a Lasota–Wazewska model. Electron. J. Differ. Equ. 2016, 2016, 1–17. [Google Scholar]
  12. Shao, J. Pseudo-almost periodic solutions for a Lasota–Wazewska model with anoscillating death rate. Appl. Math. Lett. 2015, 43, 90–95. [Google Scholar] [CrossRef]
  13. Xiao, S. Delay effect in the Lasota–Wazewska modelwith multiple time-varying delays. Int. J. Biomath. 2018, 11, 1850013. [Google Scholar] [CrossRef]
  14. Benchohra, M.; Henderson, J.; Ntouyas, S. Impulsive Differential Equations and Inclusions; Hindawi Publishing Corporation: New York, NY, USA, 2006. [Google Scholar]
  15. Kiskinov, H.; Zahariev, A. Nonlinear impulsive differential equations with weighted exponential or ordinary dichotomous linear part in a Banach space. Int. J. Differ. Equ. 2015, 748607. [Google Scholar] [CrossRef]
  16. Li, X.; Bohner, M. An impulsive delay differential inequality and applications. Comput. Math. Appl. 2012, 64, 1875–1881. [Google Scholar] [CrossRef] [Green Version]
  17. Li, X.; Bohner, M.; Wang, C.K. Impulsive differential equations: Periodic solutions and applications. Autom. J. IFAC 2015, 52, 173–178. [Google Scholar] [CrossRef]
  18. Li, X.; Ho, D.W.C.; Cao, J. Finite-time stability and settling-time estimation of nonlinear impulsive systems. Automatica 2019, 99, 361–368. [Google Scholar] [CrossRef]
  19. Li, X.; Yang, X.; Huang, T. Persistence of delayed cooperative models: Impulsive control method. Appl. Math. Comput. 2019, 342, 130–146. [Google Scholar] [CrossRef]
  20. 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]
  21. Liu, X.; Takeuchi, Y. Periodicity and global dynamics of an impulsive delay Lasota–Wazewska model. J. Math. Anal. Appl. 2007, 327, 326–341. [Google Scholar] [CrossRef]
  22. Stamov, G.T. On the existence of almost periodic solutions for the impulsive Lasota–Wazewska model. Appl. Math. Lett. 2009, 22, 516–520. [Google Scholar] [CrossRef]
  23. Wang, L.; Yu, M.; Niu, P. Periodic solution and almost periodic solution of impulsive Lasota–Wazewska model with multiple time-varying delays. Comput. Math. Appl. 2012, 64, 2383–2394. [Google Scholar] [CrossRef]
  24. Yan, J. Existence and global attractivity of positive periodic solution for an impulsive Lasota–Wazewska model. J. Math. Anal. Appl. 2003, 279, 111–120. [Google Scholar] [CrossRef]
  25. Yao, Z. Uniqueness and exponential stability of almost periodic positive solution for Lasota–Wazewska model with impulse and infinite delay. Math. Methods Appl. Sci. 2015, 38, 677–684. [Google Scholar] [CrossRef]
  26. Benchohra, M.; Henderson, J.; Ntouyas, S.K.; Ouahab, A. Impulsive functional differential equations with variable times. Comput. Math. Appl. 2004, 47, 1659–1665. [Google Scholar] [CrossRef] [Green Version]
  27. Graef, J.R.; Ouahab, A. Global existence and uniqueness results for impulsive functional differential equations with variable times and multiple delays. Dyn. Contin. Discret. Impuls. Syst. Ser. A Math. Anal. 2009, 16, 27–40. [Google Scholar]
  28. Song, Q.; Yang, X.; Li, C.; Huang, T.; Chen, X. Stability analysis of nonlinear fractional-order systems with variable-time impulses. J. Franklin Inst. 2017, 354, 2959–2978. [Google Scholar] [CrossRef]
  29. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations, 1st ed.; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006; ISBN 0444518320, 9780444518323. [Google Scholar]
  30. Podlubny, I. Fractional Differential Equations, Mathematics in Science and Engineering, 1st ed.; Academic Press: San Diego, CA, USA, 1999; ISBN 9780125588409, 9780080531984. [Google Scholar]
  31. 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]
  32. Bohner, M.; Rahman, G.; Mubeen, S.; Nisar, K.S. A further extension of the extended Riemann-Liouville fractional derivative operator. Turk. J. Math. 2018, 42, 2631–2642. [Google Scholar] [CrossRef]
  33. Boyadzhiev, D.; Kiskinov, H.; Zahariev, A. Integral representation of solutions of fractional system with distributed delays. Integral Transforms Spec. Funct. 2018, 29, 725–744. [Google Scholar] [CrossRef]
  34. Skovranek, T.; Podlubny, I.; Petráš, I. Modeling of the national economics instate-space: A fractional calculus approach. Econ. Model. 2012, 29, 1322–1327. [Google Scholar] [CrossRef]
  35. Veselinova, M.; Kiskinov, H.; Zahariev, A. Stability analysis of neutral linear fractional system with distributed delays. Filomat 2016, 30, 841–851. [Google Scholar] [CrossRef] [Green Version]
  36. Stamov, G.; Stamova, I. Impulsive fractional-order neural networks with time-varying delays: Almost periodic solutions. Neural Comput. Appl. 2017, 28, 3307–3316. [Google Scholar] [CrossRef]
  37. 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]
  38. Stamova, I.M.; Stamov, G.T. Functional and Impulsive Differential Equations of Fractional Order: Qualitative Analysis and Applications, 1st ed.; CRC Press: Boca Raton, FL, USA, 2017; ISBN 9781498764834. [Google Scholar]
  39. Wang, J.R.; Fečkan, M.; Zhou, Y. On the new concept of solutions and existence results for impulsive fractional evolution equations. Dyn. Partial Differ. Equ. 2011, 8, 345–361. [Google Scholar] [Green Version]
  40. 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]
  41. 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]
  42. Das, S.; Gupta, P.K. A mathematical model on fractional Lotka-Volterra equations. J. Theor. Biol. 2011, 277, 1–6. [Google Scholar] [CrossRef]
  43. El-Sayed, A.M.A.; El-Mesiry, A.E.M.; El-Saka, H.A.A. On the fractional-order logistic equation. Appl. Math. Lett. 2007, 20, 817–823. [Google Scholar] [CrossRef] [Green Version]
  44. 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]
  45. Akhmet, M.U. Integral manifolds of differential equations with piecewise constant argument of generalized type. Nonlinear Anal. 2007, 66, 367–383. [Google Scholar] [CrossRef] [Green Version]
  46. Bernfeld, S.R.; Corduneanu, C.; Ignatyev, A.O. On the stability of invariant sets of functional differential equations. Nonlinear Anal. 2003, 55, 641–656. [Google Scholar] [CrossRef]
  47. Smale, S. Stable manifolds for differential equations and diffeomorphisms. Ann. Scuola Norm. Super. Pisa 1963, 3, 97–116. [Google Scholar]
  48. Stamov, G. Lyapunov’s functions and existence of integral manifolds for impulsive differential systems with time-varying delay. Methods Appl. Anal. 2009, 16, 291–298. [Google Scholar]
  49. Blythe, S.P.; Nisbet, R.M.; Gurney, W.S.C. Instability and complex dynamic behaviour in population models with long time delays. Theor. Popul. Biol. 1982, 22, 147–176. [Google Scholar] [CrossRef]
  50. Arino, O.; Kimmel, M. Stability analysis of models of cell production systems. Mathematical models in medicine: Diseases and epidemics. Part 2. Math. Model. 1986, 7, 1269–1300. [Google Scholar] [CrossRef]
  51. Johnston, L. Non-Hodgkin’s Lymphomas: Making Sense of Diagnosis, Treatment, and Options, 1st ed.; O’Reilly & Associates, Inc.: Boston, MA, USA, 1999; ISBN 9781565924444. [Google Scholar]
  52. Pujo-Menjouet, L.; Bernard, S.; Mackey, M.C. Long period oscillations in a G0 model of hematopoietic stem cells. SIAM J. Appl. Dyn. Syst. 2005, 4, 312–332. [Google Scholar] [CrossRef]
  53. Li, Y.; Chen, Y.Q.; Podlubny, I. Stability of fractional-order nonlinear dynamic systems: Lyapunov direct method and generalized Mittag–Leffler stability. Comput. Math. Appl. 2010, 59, 1810–1821. [Google Scholar] [CrossRef]
  54. Liu, X.; Ballinger, G. Boundedness for impulsive delay differential equations and applications to population growth models. Nonlinear Anal. 2003, 53, 1041–1062. [Google Scholar] [CrossRef]
  55. Luo, Z.; Shen, J. Stability and boundedness for impulsive functional differential equations with infinite delays. Nonlinear Anal. 2001, 46, 475–493. [Google Scholar] [CrossRef]
  56. Stamova, I.M. Boundedness of impulsive functional differential equations with variable impulsive perturbations. Bull. Aust. Math. Soc. 2008, 77, 331–345. [Google Scholar] [CrossRef]
  57. Zhang, Y.; Sun, J. Boundedness of the solutions of impulsive differential systems with time-varying delay. Appl. Math. Comput. 2004, 154, 279–288. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Stamov, G.; Stamova, I. Impulsive Delayed Lasota–Wazewska Fractional Models: Global Stability of Integral Manifolds. Mathematics 2019, 7, 1025. https://doi.org/10.3390/math7111025

AMA Style

Stamov G, Stamova I. Impulsive Delayed Lasota–Wazewska Fractional Models: Global Stability of Integral Manifolds. Mathematics. 2019; 7(11):1025. https://doi.org/10.3390/math7111025

Chicago/Turabian Style

Stamov, Gani, and Ivanka Stamova. 2019. "Impulsive Delayed Lasota–Wazewska Fractional Models: Global Stability of Integral Manifolds" Mathematics 7, no. 11: 1025. https://doi.org/10.3390/math7111025

APA Style

Stamov, G., & Stamova, I. (2019). Impulsive Delayed Lasota–Wazewska Fractional Models: Global Stability of Integral Manifolds. Mathematics, 7(11), 1025. https://doi.org/10.3390/math7111025

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