Next Article in Journal
Canonical Almost Geodesic Mappings of the First Type of Spaces with Affine Connections onto Generalized m-Ricci-Symmetric Spaces
Next Article in Special Issue
Integral Comparison Criteria for Half-Linear Differential Equations Seen as a Perturbation
Previous Article in Journal
Reliable Learning with PDE-Based CNNs and DenseNets for Detecting COVID-19, Pneumonia, and Tuberculosis from Chest X-Ray Images
Previous Article in Special Issue
Optimal Control for a Nonlocal Model of Non-Newtonian Fluid Flows
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Stability Concepts of Riemann-Liouville Fractional-Order Delay Nonlinear Systems

1
Department of Mathematics, Texas A&M University-Kingsville, Kingsville, TX 78363, USA
2
Mathematics Department, Florida Institute of Technology, Melbourne, FL 32901, USA
3
Faculty of Mathematics and Informatics, University of Plovdiv “Paisii Hilendarski”, 4000 Plovdiv, Bulgaria
4
School of Mathematics, Statistics and Applied Mathematics, National University of Ireland, H91 CF50 Galway, Ireland
*
Author to whom correspondence should be addressed.
Mathematics 2021, 9(4), 435; https://doi.org/10.3390/math9040435
Submission received: 4 January 2021 / Revised: 7 February 2021 / Accepted: 16 February 2021 / Published: 22 February 2021
(This article belongs to the Special Issue Nonlinear Equations: Theory, Methods, and Applications)

Abstract

:
First, we set up in an appropriate way the initial value problem for nonlinear delay differential equations with a Riemann-Liouville (RL) fractional derivative. We define stability in time and generalize Mittag-Leffler stability for RL fractional differential equations and we study stability properties by an appropriate modification of the Razumikhin method. Two different types of derivatives of Lyapunov functions are studied: the RL fractional derivative when the argument of the Lyapunov function is any solution of the studied problem and a special type of Dini fractional derivative among the studied problem.

1. Introduction

Various processes with anomalous dynamics in science and engineering can be formulated mathematically using fractional differential operators ([1,2,3]). When the Riemann-Liouville (RL) fractional derivative is applied in differential equations, the statement of initial conditions is important. It is worth mentioning that the physical and geometric interpretations of operations of fractional integration and differentiation were suggested by Podlubny [4]. Fractional differential equations in terms of the RL derivative require initial conditions expressed in terms of initial values of fractional derivatives of the unknown function ([5,6]). In [7], it was shown that the initial conditions for fractional differential equations with RL derivatives expressed in terms of fractional derivatives has physical meaning. In fact, it was shown that for any physically realistic model, zero initial conditions will be found for a continuous loading program or even in the case of a step discontinuity. Nonzero conditions will only be found in the case of an impulse, and this type of process can be found in physics, chemistry, engineering, biology and economics. In the case of zero initial conditions, the RL, Grünwald-Letnikov (GL) and Caputo fractional derivatives coincide ([4]). For this reason, some authors either study Caputo derivatives or use RL derivatives but avoid the problem of initial values of fractional derivatives by treating only the case of zero initial conditions. This leads to the consideration of mathematical correct problems, but without taking the physical nature of the described process into account. Sometimes, such as in the case of impulse response, nonzero initial conditions appear (see, for example, [7]).
In connection with the main idea of stability properties, we will consider in this paper nonzero initial conditions for RL fractional equations, and we will define in an appropriate way stability properties that are slightly different than those for Caputo fractional differential equations.
Note that stability properties of delay differential equations can be considered by an application of the Lyapunov-Krasovskii method by functionals or by the Razumikhin method by Lyapunov functions. It is worth mentioning that both mentioned methods are applied for stability study of Caputo fractional delay differential equations (see, for example, [8,9,10,11]).
In the case of delay fractional differential equations with the RL fractional derivative, following the idea of initial conditions in ordinary delay differential equations and the above-mentioned idea concerning the initial condition for RL fractional differential equations without any delay, we will set up initial conditions in an appropriate way. Note that any solution of the defined initial conditions with RL fractional derivatives is not continuous at zero (the initial point), which is the same as in the case without any delay. Delay RL fractional differential equations are set up and studied in [12], but the initial condition does not correspond to the idea of the case of delay differential equations with ordinary derivatives (the lower bound of the RL fractional derivative coincides with the left end side of the initial interval).
Asymptotic stability for RL fractional differential equations with delays is studied in [13,14,15], but only the autonomous case is considered. A Lyapunov functional and its integer order derivative is applied. This functional is similar to the one used in the theory of differential equations with ordinary derivative and delay. On one hand, the application of ordinary derivative of the Lyapunov functional is not similar to the used fractional derivatives in the equation; on the other hand, it leads to some restrictions on both the delay and the right side parts of the equation ([16]). Additionally, in [15], the initial condition is not adequately associated with the RL fractional derivative. RL fractional equations with delays were studied recently in [17,18], but there are unclear parts in the statement of the problem (the lower limit of the RL derivative is different than the initial time point) as well as in the initial condition (the RL fractional integral has no meaning, compared with [5] at the initial time). The Razumikhin method is applied to RL fractional differential equations in [11], but the initial condition is not connected with the RL fractional derivative.
In this paper, the initial value problem for nonlinear delay differential equations with the RL fractional derivative is studied. Based on the arguments in the books [5,6], we set up initial conditions expressed in terms of initial values of fractional derivatives of the unknown function. Any solution of the defined initial conditions with RL fractional derivatives is not continuous at zero (the initial point). We require a new definition for stability excluding a small interval around zero. We define stability in time and generalize Mittag-Leffler stability in time for RL fractional differential equations. The stability properties of the zero solution are studied by Lyapunov functions. An appropriate modification of the Razumikhin method is suggested. Two types of derivatives of Lyapunov functions are applied: the RL fractional derivative when the argument of the Lyapunov function is a solution of the studied problem and the Dini fractional derivative among the studied problem.
The main contribution in the paper could be summarized as follows:
-
the initial conditions connected with the RL fractional derivative are set up in an appropriate way;
-
new types of stability connected with the type of initial conditions are defined;
-
the RL fractional modification of the Razumikhin method is presented;
-
new sufficient conditions for the defined stability are obtained;
-
two types of fractional derivatives of the Lyapunov functions are used.

2. Preliminary Notes

In this paper, we will use the following definitions that are well known in the literature ([5,19]):
-
Riemann–Liouville fractional integral of order q > 0
0 I t q m ( t ) = 1 Γ q 0 t m ( s ) ( t s ) 1 q d s , t 0 ,
where Γ ( . ) is the gamma function.
Note that the notation 0 D t α m ( t ) = 0 I t α m ( t ) is sometimes used.
-
Riemann–Liouville fractional derivative of order q ( 0 , 1 )
0 R L D t q m ( t ) = 1 Γ 1 q d d t t 0 t t s q m ( s ) d s , t 0
-
Grünwald–Letnikov fractional derivative of order q ( 0 , 1 )
0 G L D t q m ( t ) = lim h 0 + 1 h q j = 0 [ t ] h ( 1 ) j q C j m ( t j h ) , t 0 ,
where
q C j = q ( q 1 ) ( q 2 ) ( q j 1 ) j ! = Γ ( q + 1 ) j ! Γ ( q j + 1 ) .
Remark 1.
If m C ( [ 0 , T ] , R ) , then (see, Theorem 2.25 [19]),
0 G L D t q m ( t ) = 0 R L D t q m ( t ) , t ( 0 , T ] .
The fractional derivatives for scalar functions could be easily generalized to the vector case, by taking fractional derivatives with the same fractional order for all components.
Define the class
C 1 q ( [ 0 , T ] , R ) = { m : [ 0 , T ] R : t 1 q m ( t ) C ( [ 0 , T ] , R ) } , with T .
We will provide some well-known results from the literature (see, for example, [5,19]):
Proposition 1.
For q ( 0 , 1 ) , β > 0 :
0 R L D t q t β 1 = Γ ( β ) Γ ( β q ) t β q 1 .
Proposition 2.
For q ( 0 , 1 ) :
0 R L D t q t q = Γ ( 1 q ) ,
0 R L D t q 1 = 1 Γ ( 1 q ) t q ,
0 R L D t q t q 1 = 0 .
Proposition 3
(Property 4 [20]). If the inequalities 0 R L D t q x ( t ) 0 R L D t q y ( t ) , t ( 0 , b ] and 0 R L D t q 1 x ( t ) | t = 0 0 R L D t q 1 y ( t ) | t = 0 hold, then x ( t ) y ( t ) , t ( 0 , b ] .
Remark 2.
As it is mentioned in [20] (see Example 1 [20]), a function might not be differentiable at one point in the classical sense, but it is RL differentiable. The positive RL fractional derivative 0 R L D t q m ( t ) > 0 of order q ( 0 , 1 ) only means that the RL fractional integral 0 R L I t 1 q m ( t ) is monotonously increasing with respect to t and it does not imply that the function m ( t ) is monotonously increasing.
So, we cannot regard RL and Caputo derivatives as the generalization of the ordinary derivative in a rigorous mathematical way.
Proposition 4
(Lemma 2.3 [21]). Let m C 1 q ( [ 0 , T ] , R ) . Suppose that for any t 1 ( 0 , T ] , we have m ( t 1 ) = 0 and m ( t ) < 0 for 0 t < t 1 . Then, it follows that 0 R L D t q m ( t ) | t = t 1 0 .
Proposition 5.
The initial value problem
a 0 R L D t q x ( t ) = A x ( t ) + h ( t ) for t > 0 , a 0 R L I t 1 q x ( t ) | t = 0 = x 0 ,
has a unique solution
x ( t ) = x 0 t q 1 E q , q ( A t q ) + 0 t ( t s ) q 1 E q , q ( A ( t s ) q ) h ( s ) d s ,
where E q , q ( z ) = k = 0 z k Γ ( q ( k + 1 ) ) is the two-parameter Mittag-Leffler function.
Now, based on Proposition 5, we will illustrate the importance of the initial condition when the RL fractional derivative is used in the equation. For simplicity, we will consider the case of equations without any delays.
Remark 3.
Consider the scalar linear RL fractional equation ( q ( 0 , 1 ) )
0 R L D t q x ( t ) = x ( t ) for t > 0 ,
It is well known that 0 R L D t q t q 1 E q , q ( t q ) = t q 1 E q , q ( t q ) , i.e., the solution of the above RL fractional differential equation, is x ( t ) = c t q 1 E q , q ( t q ) where c is a real constant.
Now, consider the initial condition x ( 0 ) = k where k is a real constant. However, t q 1 E q , q ( t q ) | t = 0 + = k . This illustrates that the initial condition of the type x ( 0 ) = k is not applicable for RL fractional equations (see, for example, [11]).
Now, consider the initial condition t 1 q x ( t ) | t = 0 + = k where k is a real constant. Then, t 1 q x ( t ) | t = 0 + = t 1 q c t q 1 E q , q ( t q ) | t = 0 + = c E q , q ( 0 ) = c Γ ( q ) , i.e., the initial condition t 1 q x ( t ) | t = 0 + = k has a meaning for the RL fractional derivative with k = c Γ ( q ) .
The practical definition of the initial condition of fractional differential equations with RL derivatives is based on the following result:
Lemma 1
([2]). Let q ( 0 , 1 ) and b > 0 , m : [ 0 , b ] R be a Lebesgue measurable function.
(a) 
If there exists a.e. a limit lim t 0 + [ t 1 q m ( t ) ] = c R , then there also exists a limit
0 I t 1 q m ( t ) | t = 0 = 0 D t q 1 m ( t ) | t = 0 : = lim t 0 + 1 Γ 1 q 0 t m ( s ) ( t s ) q d s = c Γ ( q ) .
(b) 
If there exists a.e. a limit lim t 0 + 0 I t 1 q m ( t ) = c R , and if there exists the limit lim t 0 + [ t 1 q m ( t ) ] , then
lim t 0 + [ t 1 q m ( t ) ] = c Γ ( q ) .

3. Statement of the Problem

Consider the following nonlinear Riemann–Liouville fractional delay differential equation (RLFrDDE) of fractional order q ( 0 , 1 ) :
0 R L D t q x ( t ) = f ( t , x t ) for t > 0 ,
with initial conditions
x ( t ) = ϕ ( t ) , for t [ τ , 0 ] , lim t 0 + [ t 1 q x ( t ) ] = ϕ ( 0 ) Γ ( q ) ,
where x R n , x t ( Θ ) = x ( t + Θ ) , Θ [ τ , 0 ] , ϕ : [ τ , 0 ] R and f : R + × R n R n .
Remark 4.
According to Lemma 1, the second equation in the initial conditions (3) could be replaced by the equality 0 I t 1 q x ( t ) | t = 0 = ϕ ( 0 ) .
We denote the solution of the initial value problem (IVP) for RLFrDDE (2) and (3) by x ( t ) = x ( t ; ϕ ) for t τ . In this paper, we will assume that the function f is such that for any continuous initial function ϕ the IVP for RLFrDDE (2) and (3) has a solution. Note that some existence and uniqueness results to RL fractional differential equations with delay were obtained in [22].
For any ϕ C ( [ τ , 0 ] , R n ) , we denote | | ϕ | | 0 = max t [ τ , 0 ] | | ϕ ( t ) | | where | | . | | is a norm in R n .
We will introduce the following conditions
Hypothesis 1
(H1). The function f C [ R + × R n , R n ] is such that for any initial function ϕ C ( [ τ , 0 ] , R n ) ) , the corresponding IVP for RLFrDDE (2) and (3) has a solution x ( t ; ϕ ) C 1 q ( [ 0 , ) , R n ) ;
Hypothesis 2
(H2). f ( t , 0 ) = 0 for all t R + .
Remark 5.
If ϕ ( t ) 0 and condition (H2) is satisfied, then, because of the equality 0 R L D t q 0 = 0 , the IVP for RLFrDDE (2) and (3) has the zero solution.
We will give the basic definitions for stability:
Definition 1.
The zero solution of RLFrDDE (2) and (3) (with the zero initial function) is said to be
-
stable in time if for any number ε > 0 there exist numbers δ > 0 and T ε > 0 depending on ε such that for any initial functions ϕ C 0 : | | ϕ | | 0 < δ , the corresponding solution x ( t ; ϕ ) of IVP (2) and (3) satisfies | | x ( t ; ϕ ) | | < ε for t T ε ;
-
asymptotically stable if it is stable in time and additionally x ( t ; ϕ ) 0 as t + ;
-
generalized Mittag-Leffler stable in time if there exist positive numbers λ , b and γ ( 0 , 1 ) and a locally Lipschitz function h C ( [ 0 , ) , [ 0 , ) ) : h ( 0 ) = 0 such that for any ε > 0 , there exists T ε > 0 the solution of IVP (2) and (3) satisfies
| | x ( t ; ϕ ) | | ε { h ( | | ϕ | | 0 ) t γ E q , 1 γ ( λ t q ) } b , t T ε .
As an example to discuss stability in time, we consider a scalar RL fractional equation without any delay, whose exact solution is known.
Example 1.
According to Proposition 5, the scalar initial value problem
0 R L D t q x ( t ) = a x ( t ) for t > 0 , 0 R L I t 1 q x ( t ) | t = 0 = x 0
has a unique solution
x ( t ) = x 0 t q 1 E q , q ( a t q ) .
Therefore, the zero solution is generalized Mittag-Leffler stable with γ = 1 q , m ( u ) = u and λ = a .
The zero solution is stable in time because for any ε > 0 there exist δ , T ε > 0 such that δ T ε 1 q E q , q ( a T ε q ) = ε . At the same time, the zero solution is not stable (in the regular sense) because lim t 0 t q 1 E q , q ( a t q ) = and | t q 1 E q , q ( a t q ) | < ε cannot be satisfied for values of t close to 0.
For example, if a = 1 , q = 0.5 , ε = 0.1 , δ = 0.2 > ε then T ε 0.257459 , if a = 1 , q = 0.5 , ε = 0.1 , δ = 0.01 < ε then T ε 0.00265746 .
Throughout the paper, we shall use the class
K = { w C ( R + , R + ) : w ( s ) is strictly increasing and w ( 0 ) = 0 }
and
B ¯ ( ρ ) = { x R n : | | x | | ρ } .

4. Stability of Nonlinear RL Fractional Differential Equations

4.1. Lyapunov Functions and Their Derivatives

One approach to study stability properties of nonlinear RL fractional differential equations is based on the application of Lyapunov functions and an appropriate modification of the Razumikhin method. The first step is to define a Lyapunov function. The second step is to define its derivative among the studied equation.
We will use the following class of functions called Lyapunov functions:
Definition 2
([8]). Let J = [ τ , T ) , T , be a given interval, and Δ R n , 0 Δ be a given set. We will say that the function V ( t , x ) : J × Δ R + belongs to the class Λ ( J , Δ ) if V ( t , x ) is continuous on J / { 0 } × Δ , and it is locally Lipschitzian with respect to its second argument.
In our study, we will use the Razumikhin condition for the Lyapunov function V Λ ( J , Δ ) and any ψ C ( [ τ , 0 ] , R n ) :
V ( t + Θ , ψ ( Θ ) ) V ( t , ψ ( 0 ) ) , Θ [ τ , 0 ] .
We will give a brief overview of the derivatives of Lyapunov functions among solutions of fractional differential equations in the literature. There are three main types of derivatives of Lyapunov functions from the class Λ ( J , Δ ) used in the literature to study stability properties of solutions of fractional differential Equation (2):
-
RL fractional derivative—Let x ( t ) : [ 0 , T ) Δ be a solution of the IVP for the RLFrDDE (2) and (3). Then, we consider
0 R L D q V ( t , x ( t ) ) = 1 Γ 1 q d d t 0 t t s q V ( s , x ( s ) ) d s , t [ 0 , T ) .
-
Dini fractional derivative—Let ϕ : [ τ , 0 ] Δ . Then, consider (see [8])
D ( 2 ) + V ( t , ϕ ( 0 ) , ϕ ) = lim   sup h 0 1 h q V ( t , ϕ ( 0 ) ) r = 1 [ t h ] ( 1 ) r + 1 q C r V ( t r h , ϕ ( 0 ) h q f ( t , ϕ 0 ) ) ,
where q C r is defined by (1) and ϕ 0 ( Θ ) = ϕ ( Θ ) , Θ [ τ , 0 ] .
The Dini fractional derivative is applicable for continuous Lyapunov functions.
We will provide an example concerning some Lyapunov functions and their Dini fractional derivatives. To simplify the calculations and to emphasize on the derivative, we will consider the scalar case, i.e., n = 1 .
Example 2.
Let V ( t , x ) = m ( t ) x 2 where m C ( R + , R + ) and ϕ C ( [ τ , 0 ] , R ) . Apply (5) and obtain
D ( 2 ) + V ( t , ϕ ( 0 ) , ϕ ) = lim   sup h 0 1 h q [ m ( t ) ( ϕ ( 0 ) ) 2 r = 1 [ t h ] ( 1 ) r + 1 q C r m ( t r h ) ( ϕ ( 0 ) h q f ( t , ϕ 0 ) ) 2 ] = lim   sup h 0 1 h q [ [ m ( t ) ( ϕ ( 0 ) ) 2 m ( t ) ( ϕ ( 0 ) ) 2 h q f ( t , ϕ 0 ) ) 2 ( ϕ ( 0 ) ) 2 h q f ( t , ϕ 0 ) ) 2 ( 1 ) 0 + 1 q C 0 m ( t 0 h ) ( ϕ ( 0 ) h q f ( t , ϕ ) ) 2 r = 1 [ t h ] ( 1 ) r + 1 q C r m ( t r h ) ] = lim   sup h 0 1 h q [ m ( t ) h q f ( t , ϕ 0 ) ( 2 ϕ ( 0 ) h q f ( t , ϕ 0 ) ) + ( ϕ ( 0 ) ) 2 h q f ( t , ϕ 0 ) ) 2 r = 0 [ t h ] ( 1 ) r q C r m ( t r h ) ] = 2 ϕ ( 0 ) m ( t ) f ( t , ϕ ) + ( ϕ ( 0 ) ) 2 0 R L D q m ( t ) .
Special case 1. Let V ( t , x ) = t 1 q x 2 for t > 0 , x R . According to Proposition 1 with β = 2 q , we get 0 R L D t q t 1 q = Γ ( 2 q ) Γ ( 2 2 q ) t 1 2 q and from (6), we get
D ( 2 ) + V ( t , ϕ ( 0 ) , ϕ ) = t 1 q ( 2 ϕ ( 0 ) f ( t , ϕ ) + ( ϕ ( 0 ) ) 2 Γ ( 2 q ) Γ ( 2 2 q ) t q ) .
Special case 2. Let V ( t , x ) = x 2 for x R . From Proposition 2, with m ( t ) 1 , we obtain
D ( 2 ) + V ( t , ϕ ( 0 ) , ϕ ) = 2 ϕ ( 0 ) f ( t , ϕ ) + ( ϕ ( 0 ) ) 2 t q Γ ( 1 q ) .
Special case 3. Let V ( t , x ) = t q 1 x 2 for t > 0 , x R . According to Proposition 2, we get from (6)
D ( 2 ) + V ( t , ϕ ( 0 ) , ϕ ) = 2 t q 1 ϕ ( 0 ) f ( t , ϕ ) .
Special case 4. Let V ( t , x ) = t 1 q m ( t ) x 2 where A t 1 q m ( t ) B for all t 0 . For example, if m ( t ) = 1 t 1 q + 1 1 + t 1 q , then t 1 q m ( t ) = 1 + t 1 q 1 + t 1 q ( 1 , 2 ) , 0 R L D t q 1 + t 1 q 1 + t 1 q exists and
D ( 2 ) + V ( t , ϕ ( 0 ) , ϕ ) = 2 ϕ ( 0 ) 1 + t 1 q 1 + t 1 q f ( t , ϕ ) + ( ϕ ( 0 ) ) 2 2 t q Γ ( 1 q ) 0 R L D t q 1 1 + t 1 q .
Remark 6.
Note that V ( t , x ( t ) ) = t 1 q m ( t ) i 1 n x i 2 ( t ) C 1 q ( [ 0 , ) , R + ) , m C ( [ 0 , ) , R + ) for any solution x ( t ) = ( x 1 ( t ) , x 2 ( t ) , , x n ( t ) ) of (2) and (3) because
t 1 q V ( t , x ( t ) ) | t = 0 = m ( t ) i = 1 n ( t 1 q x i ( t ) ) 2 | t = 0 = m ( 0 ) i = 1 n ϕ i 2 ( 0 ) Γ 2 ( q ) = m ( 0 ) Γ 2 ( q ) | | ϕ | | 0 2 .
If m ( t ) 1 , then the function V ( t , x ( t ) ) = t 1 q i 1 n x i 2 ( t ) C 1 q ( [ 0 , ) , R + ) for any solution x ( t ) of (2) and (3).
The quadratic Lyapunov function V ( t , x ( t ) ) = i = 1 n x i ( t ) 2 is not from the set C 1 q ( [ 0 , ) , R + ) for any solution x ( t ) of (2) and (3), because
t 1 q V ( t , x ( t ) ) | t = 0 = t 1 q i = 1 n x i ( t ) 2 | t = 0 = t q 1 i = 1 n ( t 1 q x i ( t ) ) 2 | t = 0 = .
The functions V ( t , x ( t ) ) = t q 1 i = 1 n x i 2 ( t ) is not from the set C 1 q ( [ 0 , ) , R + ) for any solution x ( t ) of (2) and (3).
We will study stability properties of the zero solution of RLFrDDE (2) by an application of both defined types of fractional derivatives of Lyapunov functions.

4.2. Stability by the RL Fractional Derivative of Lyapunov Functions

We will obtain some sufficient conditions for stability with applications of the RL fractional derivative of Lyapunov functions.
Theorem 1.
Let conditions (H1) and (H2) be satisfied, and there exists a function V Λ ( R + , R n ) such that
(i) 
for any ε > 0 there exists T ε > 0 such that
ε a ( | | x | | ) V ( t , x ) for t > T ε , x R n ,
where a K ;
(ii) 
there exists an increasing function g C ( R + , R + ) such that for any function y C 1 q ( [ 0 , ) , R n ) : t 1 q y ( t ) | t = 0 + = y 0 the inequality
t 1 q V ( t , y ( t ) ) | t = 0 + = lim t 0 + t 1 q V ( t , y ( t ) ) g ( | | y 0 | | )
holds;
(iii) 
for any point t > 0 such that ( t + Θ ) 1 q V ( t + Θ , x ( t + Θ ) ) < t 1 q V ( t , x ( t ) ) for Θ ( min { t , τ } , 0 ) , the RL fractional derivative 0 R L D t q V ( t , x ( t ) ) exists and the inequality
0 R L D t q V ( t , x ( t ) ) = 1 Γ 1 q d d t 0 t t s q V ( s , x ( s ) ) d s < 0 ,
holds where x ( t ) is a solution of the IVP for RLFrDDE (2) and (3).
Then, the zero solution of (2) and (3) with the zero initial function is stable in time.
Proof. 
Let ϵ > 0 be an arbitrary number. According to condition (i), there exists T ˜ ε > 0 such that inequality (9) holds for t > T ˜ ε , x R n .
Let δ > 0 be such that for | | ϕ | | 0 < δ , the inequality g ( | | ϕ | | 0 Γ ( q ) ) < ε holds.
Consider the solution x ( t ) = x ( t ; ϕ ) of the IVP for RLFrDDE (2) and (3) with the initial function ϕ : | | ϕ | | 0 < δ .
From condition 2(ii) with y ( t ) x ( t ) and y 0 = ϕ ( 0 ) Γ ( q ) , we get
t 1 q V ( t , x ( t ) ) | t = 0 + = lim t 0 + t 1 q V ( t , x ( t ) ) g ( | | ϕ ( 0 ) | | Γ ( q ) ) g ( | | ϕ | | 0 Γ ( q ) ) < ε .
Therefore, there exists δ 1 = δ 1 ( ε ) > 0 such that V ( t , x ( t ) ) < ε t q 1 for t < δ 1 .
Consider the function H ( t ) = ε t q 1 C 1 q ( [ 0 , ) , R + ) and lim t H ( t ) = 0 . Therefore, there exists T ^ ε > δ 1 such that
H ( t ) < ε a ( ε ) for t > T ^ ε .
Define the function m ( t ) = V ( t , x ( t ) ) for t 0 . From (11), it follows that m C 1 q ( [ 0 , ) , R + ) .
We will prove that
m ( t ) < H ( t ) , t > 0 .
Note that the inequality (13) holds for t ( 0 , δ 1 ) . Assume inequality (13) is not true for all t > 0 . Therefore, there exists a point ξ δ 1 > 0 such that
m ( ξ ) = H ( ξ ) , and m ( t ) < H ( t ) , t ( 0 , ξ ) .
Therefore, m ( t ) H ( t ) C 1 q ( [ 0 , ξ ] , R ) . According to Proposition 4 with t 1 = ξ , the inequality 0 R L D t q m ( t ) H ( t ) | t = ξ 0 holds. From Proposition 2, we get 0 R L D t q t q 1 = 0 and therefore,
0 R L D t q m ( t ) | t = ξ = 0 R L D t q m ( t ) H ( t ) | t = ξ 0 .
Case 1. Let ξ > τ . Then, min { ξ , τ } = τ . From (14), it follows that ξ q 1 m ( ξ ) = ε > t q 1 m ( t ) , t ( 0 , ξ ) , or ( ξ + Θ ) q 1 m ( ξ + Θ ) = ( ξ + Θ ) q 1 V ( ξ + Θ , x ( ξ + Θ ) ) < ξ q 1 m ( ξ ) = ξ q 1 V ( ξ , x ( ξ ) ) for Θ ( τ , 0 ) . According to condition 2(iii)
0 R L D t q V ( t , x ( t ) ) | t = ξ < 0 .
The inequality (16) contradicts (15).
Case 2. Let ξ τ . Then, min { ξ , τ } = ξ . From (14), it follows that ξ q 1 m ( ξ ) = ε > t q 1 m ( t ) , t ( 0 , ξ ) , or ( ξ Θ ) q 1 m ( ξ + Θ ) = ( ξ Θ ) q 1 V ( ξ + Θ , x ( ξ + Θ ) ) < ξ q 1 m ( ξ ) = ξ q 1 V ( ξ , x ( ξ ) ) for Θ ( ξ , 0 ) and the proof is similar to the one of Case 1.
From inequality (13) and condition (i), it follows that
ε a ( | | x ( t ) | | ) V ( t , x ( t ) ) = m ( t ) < H ( t ) < ε a ( ε ) for t > T ε ,
where T ε = max { T ^ ε , T ˜ ε } .
This proves the stability in time of the zero solution. □
Remark 7.
The main condition in Theorem 1 is condition (iii), which is connected with any solution of the IVP (2) and (3).
Remark 8.
According to Remark 6, the Lyapunov function V ( t , y ) = t 1 q m ( t ) i = 1 n y i 2 with m C ( R + , R + ) satisfies condition (ii) with g ( s ) = m ( 0 ) s 2 . Condition (i) is satisfied with a ( u ) = u 2 if m ( t ) C = c o n s t > 0 for t 0 .
Remark 9.
Condition (i) is different than the corresponding condition for Lyapunov functions for differential equations with ordinary derivatives as well as with Caputo fractional differential equations. This is because of the type of initial condition (see Remark 3).
Theorem 2.
Let conditions (H1) and (H2) be satisfied and there exists a function V Λ ( [ τ , ) , R n ) such that conditions (i) and (ii) of Theorem 1 hold with a ( s ) = k s p , k , p > 0 , locally Lipschitz function g C ( [ 0 , ) , [ 0 , ) ) : g ( 0 ) = 0 and( i i i ) for any point t > 0 such that V ( t + Θ , x ( t + Θ ) ) < V ( t , x ( t ) ) for Θ ( min { t , τ } , 0 ) , the RL fractional derivative 0 R L D t q V ( t , x ( t ) ) exists and the inequality
0 R L D t q V ( t , x ( t ) ) < c V ( t , x ( t ) ) ,
holds where c > 0 , x ( t ) is a solution of the IVP for RLFrDDE (2) and (3).
Then, the zero solution of (2) and (3) with the zero initial function is generalized Mittag-Leffler stable.
Proof. 
Consider any solution x ( t ) = x ( t ; ϕ ) of the IVP for RLFrDDE (2) and (3) with the initial function ϕ C ( [ τ , 0 ] , R n ) .
From condition (ii) of Theorem 1, similar to inequality (11), we get
t 1 q V ( t , x ( t ) ) | t = 0 + g | | ϕ | | 0 Γ ( q ) .
Consider the function H ( t ) = Γ ( q ) g ( | | ϕ | | 0 Γ ( q ) ) t q 1 E q , q ( c t q ) C 1 q ( [ 0 , ) , ( 0 , ) ) .
Define the function m ( t ) = V ( t , x ( t ) ) for t 0 . From (19), it follows that m C 1 q ( [ 0 , ) , R + ) .
Let ε > 0 be an arbitrary number. We will prove that
m ( t ) < H ( t ) + ε t q 1 , t > 0 .
We have t 1 q V ( t , x ( t ) ) | t = 0 + g | | ϕ | | 0 Γ ( q ) < g | | ϕ | | 0 Γ ( q ) + ε = t 1 q H ( t ) | t = 0 + + ε . Therefore, there exists δ 2 > 0 such that for t ( 0 , δ 2 ) , the inequality V ( t , x ( t ) ) H ( t ) + ε t q 1 holds.
Assume inequality (20) is not true for all t > 0 . Therefore, there exists a point ξ δ 2 > 0 such that
m ( ξ ) = H ( ξ ) + ε t q 1 , and m ( t ) < H ( t ) + ε t q 1 , t ( 0 , ξ ) .
Therefore, m ( t ) H ( t ) ε t q 1 C 1 q ( [ 0 , ξ ] , R ) . According to 0 R L D t q t q 1 = 0 and Proposition 4 with t 1 = ξ , the inequality 0 R L D t q m ( t ) H ( t ) | t = ξ 0 holds. From 0 R L D t q t q 1 E q , q ( c t q ) = c t q 1 E q , q ( c t q ) , we have,
0 R L D t q m ( t ) | t = ξ + c m ( ξ ) = 0 R L D t q m ( t ) | t = ξ + c H ( ξ ) = 0 R L D t q m ( t ) | t = ξ + c g ( | | ϕ | | 0 Γ ( q ) ) ξ q 1 E q , q ( c ξ q ) = 0 R L D t q m ( t ) H ( t ) | t = ξ 0 .
Case 1. Let ξ > τ . Then, min { ξ , τ } = τ . Therefore, m ( ξ + Θ ) = V ( ξ + Θ , x ( ξ + Θ ) ) < m ( ξ ) = V ( ξ , x ( ξ ) ) for Θ ( τ , 0 ) . Let ψ ( Θ ) = x ( ξ Θ ) C ( [ τ , 0 ] , R n ) , Θ [ τ , 0 ] . Then, V ( ξ + Θ , ψ ( Θ ) ) < V ( ξ , ψ ( 0 ) ) for Θ [ τ , 0 ) and, according to condition (iii*),
0 R L D q V ( ξ , x ( ξ ) ) < c V ( t , x ( t ) ) .
The inequality (23) contradicts (22).
Case 2. Let ξ τ . Then, min { ξ , τ } = ξ and V ( ξ + Θ , x ( ξ + Θ ) ) < V ( ξ , x ( ξ ) ) for Θ ( ξ , 0 ) and the proof is similar to the one of Case 1.
Since ε > 0 is an arbitrary number from inequality (20), it follows that
m ( t ) H ( t ) , t > 0 .
Now, let ε > 0 . According to condition (i), there exists T ε > 0 such that ( 1 ε ) p k | | x ( t ) | | p V ( t , x ( t ) ) for t T ε . Then, from inequality (24), it follows that
( 1 ε ) p k | | x ( t ) | | p V ( t , x ( t ) ) = m ( t ) < Γ ( q ) g ( | | ϕ | | 0 Γ ( q ) ) t q 1 E q , q ( c t q ) for t T ε .
This proves the generalized Mittag-Leffler stability in time of the zero solution with b = 1 p , γ = 1 q , λ = c and h ( s ) = g ( s Γ ( q ) ) k (see Definition 1). □
Corollary 1.
If all the conditions of Theorem 2 are satisfied, then the zero solution of (2) and (3) is asymptotically stable.

4.3. Stability by the Dini Fractional Derivative of Lyapunov Functions

We will study stability by the application of the defined above Dini fractional derivative of Lyapunov functions among the studied delay fractional differential equations.
Initially, we will prove a comparison result for Lyapunov functions.
Lemma 2.
Assume:
1. 
The function x ( t ) = x ( t ; ϕ ) C 1 q ( [ 0 , T ] , Δ ) is a solution of the IVP for RLFrDDE (2) and (3) with ϕ C ( [ τ , 0 ] , Δ ) where 0 < T .
2. 
The function V Λ ( [ τ , T ] , Δ ) , Δ R n , is such that:
(i) 
There exists an increasing function g C ( R + , R + ) such that the inequality t 1 q V ( t , x ( t ) ) | t = 0 + = lim t 0 + t 1 q V ( t , x ( t ) ) g ( | | ϕ ( 0 ) | | Γ ( q ) ) g ( | | ϕ | | 0 ) holds;
(ii) 
for any point t > 0 such that V ( t + Θ , x ( t + Θ ) ) < V ( t , x ( t ) ) for Θ ( min { t , τ } , 0 ) , the inequality
D ( 2 ) + V ( t , ψ ( 0 ) , ψ ) < 0
holds where D ( 2 ) + V ( t , ψ ( 0 ) , ψ ) is the Dini fractional derivative defined by (5) and ψ ( Θ ) = x ( t + Θ ) , Θ [ τ , 0 ] .
Then, V ( t , x ( t ; ϕ ) ) t q 1 g ( | | ϕ | | 0 ) for t [ 0 , T ] .
Remark 10.
Let us, for simplicity, again consider the RL fractional differential equation without any delay 0 R L D t q x ( t ) = x ( t ) for t > 0 , 0 R L I t 1 q x ( t ) | t = 0 + = c with a solution x ( t ) = c t q 1 E q , q ( t q ) (see Example 1). If we consider the quadratic Lyapunov function V ( x ) = x 2 , then D ( 2 ) + V ( t , ψ ( 0 ) , ψ ) 0 and ( c t q 1 E q , q ( t q ) ) 2 c 2 is not satisfied. However, if V ( x ) = t 1 q x 2 then D ( 2 ) + V ( t , ψ ( 0 ) , ψ ) 0 and t 1 q ( c t q 1 E q , q ( t q ) ) 2 c 2 t q 1 is satisfied. This example again illustrates the changes in the applied Lyapunov functions and their conditions in the application of RL fractional derivatives comparatively with the application in Caputo fractional derivatives.
Remark 11.
Let V ( t , x ) = t 1 q m ( t ) i = 1 n x i 2 , where m C ( [ 0 , T ] ) , T , and x ( t ) be a solution of (2) and (3). Then, the following
t 1 q V ( t , x ( t ) ) | t = 0 + = lim t 0 + m ( t ) t 2 2 q i = 1 n x i 2 ( t ) = m ( 0 ) i = 1 n ϕ i ( 0 ) 2 Γ 2 ( q )
holds, i.e., t 1 q V ( t , x ( t ) ) | t = 0 + m ( 0 ) Γ 2 ( q ) ( | | ϕ | | 0 ) 2 < and condition 2 (i) of Lemma 2 is satisfied with g ( s ) = m ( 0 ) Γ 2 ( q ) s 2 .
Proof. 
Define the function m ( t ) = V ( t , x ( t ) ) for t ( 0 , T ] . According to condition 2 (i), m C 1 q ( [ 0 , T ] , R ) .
Let H ( t ) = m ( t ) t q 1 ( B + ε ) where B = g ( | | ϕ | | 0 ) 0 and ε > 0 be an arbitrary number. We will prove
H ( t ) < 0 , t [ 0 , T ]
For t = 0 from condition 2(i), we get t 1 q H ( t ) | t = 0 = t 1 q V ( t , x ( t ) ) | t = 0 ( B + ε ) g ( | | ϕ | | 0 ) B ε < 0 , i.e., the inequality (27) is true.
Assume (27) is not true. Therefore, there exist t ( 0 , T ] such that
H ( t ) < 0 , t [ 0 , t ) , H ( t ) = 0 .
From Proposition 4, we have the inequality 0 R L D t q H ( t ) | t = t 0 . Then, applying Proposition 2 and 0 R L D t q t q 1 = 0 , we obtain
0 R L D t q H ( t ) | t = t = 0 G L D t q m ( t ) | t = t = 0 R L D t q m ( t ) | t = t 0 .
For any t ( 0 , t ] and h > 0 , we let
S ( x ( t ) , h ) = r = 1 [ t h ] ( 1 ) r + 1 q C r x ( t r h ) .
From Remark 1 and Equation (2), it follows that the function x ( t ) satisfies for t [ t 0 , t ] , the equalities 0 R L D t q x ( t ) = 0 G L D t q x ( t ) = lim   sup h 0 + 1 h q x ( t ) S ( x ( t ) , h ) = f ( t , x t ) and
lim   sup h 0 + 1 h q x ( t ) S ( x ( t ) , h ) = f ( t , x t ) .
Therefore,
S x ( t ) , h = x ( t ) h q f ( t , x t ) Λ ( h q )
or
x ( t ) h q f ( t , x t ) = S x ( t ) , h + Λ ( h q )
with | | Λ ( h q ) | | h q 0 as h 0 . Then, for any t [ 0 , t ] we obtain
m ( t ) r = 1 [ t h ] ( 1 ) r + 1 q C r m ( t r h ) = V ( t , x ( t ) ) r = 1 [ t h ] ( 1 ) r + 1 q C r V ( t r h , x ( t ) h q f ( t , x t ) + r = 1 [ t h ] ( 1 ) r + 1 q C r { V ( t r h , S x ( t ) , h + Λ ( h q ) ) V ( t r h , x ( t r h ) ) } .
Since V is locally Lipschitzian in its second argument with a Lipschitz constant L > 0 ,, we obtain
r = 1 [ t h ] ( 1 ) r + 1 q C r V ( t r h , S x ( t ) , h + Λ ( h q ) ) V ( t r h , x ( t r h ) ) L | | r = 1 [ t h ] ( 1 ) r + 1 q C r S x ( t ) , h + Λ ( h q ) ( x ( t r h ) ) | | L | | r = 1 [ t t 0 h ] ( 1 ) r + 1 q C r j = 1 [ t h ] ( 1 ) j + 1 q C j x ( t j h ) r = 1 [ t h ] ( 1 ) r + 1 q C r x ( t r h ) | | + L | | Λ ( h q ) | | | r = 1 [ t t 0 h ] ( 1 ) r + 1 q C r | = L | | r = 0 [ t h ] ( 1 ) r + 1 q C r j = 1 [ t h ] ( 1 ) j + 1 q C j x ( t j h ) | | + L | | Λ ( h q ) | | | r = 1 [ t h ] ( 1 ) r + 1 q C r | .
Substitute (32) in (31), divide both sides by h q , take the limit as h 0 + , use r = 0 q C r z r = ( 1 + z ) q if | z | 1 and we obtain for any t ( 0 , t ] the inequality
0 G L D t q m ( t ) lim h 0 + 1 h q { V ( t , x ( t ) ) r = 1 [ t h ] ( 1 ) r + 1 q C r V ( t r h , x ( t ) h q f ( t , x ( t ) ) } + L lim h 0 + | | Λ ( h q ) | | h q lim h 0 + | r = 1 [ t h ] ( 1 ) r + 1 q C r | + L lim h 0 + sup | r = 1 [ t h ] ( 1 ) r + 1 q C r | | | 1 h q j = 1 [ t h ] q C j x ( t j h ) | | lim h 0 + 1 h q V ( t , x ( t ) ) r = 1 [ t h ] ( 1 ) r + 1 q C r V ( t r h , x ( t ) h q f ( t , x ( t ) ) ) .
Let t = t . Define the function ψ ( Θ ) = x ( t + Θ ) , Θ [ τ , 0 ] . From the choice of the point t , it follows that V ( t + Θ , x ( t + Θ ) ) V ( t , x ( t ) ) , Θ [ τ , 0 ] and from inequalities (26) and (33) for t = t , we get
0 G L D t q m ( t ) | t = t lim h 0 + 1 h V ( t , ψ ( 0 ) ) r = 1 [ t h ] ( 1 ) r + 1 q C r V ( t r h , ψ ( 0 ) h q f ( t , ψ ) ) = D ( ) + V ( t , ψ ( 0 ) , ψ ) < 0 .
Now (34) contradicts (29). Therefore, inequality (27) holds for an arbitrary ε > 0 . Thus, the claim in our Lemma is true. □
Theorem 3.
Let conditions (H1) and (H2) be satisfied and there exists a function V Λ ( R + , R n ) such that V ( t , 0 ) = 0 , t 0 and
(i) 
for any ε > 0 , there exists T ε > 0 such that
ε a ( | | x | | ) V ( t , x ) for t > T ε , x R n ,
where a K ;
(ii) 
there exists an increasing function g C ( R + , R + ) such that for any function y C 1 q ( [ 0 , ) , R n ) : t 1 q y ( t ) | t = 0 + = y 0 the inequality
t 1 q V ( t , y ( t ) ) | t = 0 + = lim t 0 + t 1 q V ( t , y ( t ) ) g ( | | y 0 | | )
holds;
(iii) 
for any function ψ C ( [ τ , 0 ] , R n such that if for a point t we have V ( t + Θ , ψ ( Θ ) ) < V ( t , ψ ( 0 ) ) for Θ [ τ , 0 ) , then the inequality
D ( 2 ) + V ( t , ψ ( 0 ) , ψ ) < 0
holds.
Then, the zero solution of (2) with the zero initial function is stable.
Proof. 
Let ϵ > 0 be a positive number. According to condition (i), there exists T ˜ ε > 0 such that inequality (9) holds for t > T ˜ ε , x R n .
There exists a positive number δ = δ ( ε ) such that g ( u ) < a ( ε ) for u R + : u < δ . Choose the function ϕ C ( [ τ , 0 ] , R n ) such that | | ϕ | | 0 < δ . Consider the solution x ( t ) = x ( t ; ϕ ) of the IVP for RLFrDDE (2) and (3) with initial function ϕ and define the function m ( t ) = V ( t , x ( t ) ) for t 0 . From condition 2(ii), it follows that m C 1 q ( [ 0 , ) , R + ) .
Since lim t t q 1 = 0 , there exists T ^ ε > 0 such that t q 1 < ε a ( ε ) for t > T ^ ε . Therefore,
t q 1 g ( | | ϕ | | 0 ) < ε a ( ε ) for t > T ^ ε .
According to Lemma 2 with T = , ϕ ( 0 ) Γ ( q ) = y 0 , and ψ ( Θ ) x ( t + Θ ) , Θ [ τ , 0 ] , , we obtain the following inequality
V ( t , x ( t ) ) t q 1 g ( | | ϕ | | 0 ) , t 0 .
From inequalities (37), (38) and condition (i), it follows that ε a ( | | x ( t ) | | ) V ( t , x ( t ) ) t q 1 g ( | | ϕ | | 0 ) < ε a ( ε ) for t > T ε with T ε = max { T ^ ε , T ˜ ε } .
This proves the stability in time of the zero solution. □
Example 3.
Consider the scalar RL fractional differential equation
0 R L D t q x ( t ) = f ( t , x ( t ) , x ( t 1 ) ) for t > 0 ,
with initial conditions
x ( t ) = ϕ ( t ) , for t [ 1 , 0 ] , lim t 0 + [ t 1 q x ( t ) ] = ϕ ( 0 ) Γ ( q ) ,
where x R , ϕ C ( [ 1 , 0 ] , R ) and f : R + × R × R R is defined by
f ( t , x , y ) = 1.6 x Γ ( 2 q ) Γ ( 2 2 q ) t 1 2 q i f 0 < t 1 , x , y R 1.6 x Γ ( 2 q ) Γ ( 2 2 q ) t 1 2 q + t q ( t 1 ) 1 q y Γ ( 2 q ) Γ ( 2 2 q ) i f t > 1 , x , y R .
Consider the Lyapunov function V ( t , y ) = t 1 q y 2 . According to Remarks 8, this function satisfies condition (i) of Theorem 3 with a ( u ) = u 2 , u R + and T ε = ε 1 q . It also satisfies condition (ii) with g ( u ) = u 2 , u R + .
Now, let the function ψ C ( [ 1 , 0 ] , R n ) and the point t > 0 be such that V ( t + Θ , ψ ( Θ ) ) < V ( t , ψ ( 0 ) ) for Θ ( min { t , 1 } , 0 ) , i.e., ( t 1 ) 1 q ψ ( 1 ) 2 < t 1 q ψ ( 0 ) 2 .
For t > 1 , apply 2 ϕ ( 0 ) ( t 1 ) 1 q ϕ ( 1 ) ( t 1 ) 1 q ϕ 2 ( 0 ) + ( t 1 ) 1 q ϕ 2 ( 1 ) < 2 t 1 q ϕ 2 ( 0 ) and Example 2 (Special case 1) and obtain
D ( 2 ) + V ( t , ϕ ( 0 ) , ϕ ) = t 1 q 2 ϕ ( 0 ) f ( t , ϕ ( 0 ) , ϕ ( 1 ) ) + ( ϕ ( 0 ) ) 2 Γ ( 2 q ) Γ ( 2 2 q ) t q = t 1 q 2 ϕ ( 0 ) ( t 1 ) 1 q ϕ ( 1 ) Γ ( 2 q ) t q t 1 q Γ ( 2 2 q ) 3.2 ϕ 2 ( 0 ) Γ ( 2 q ) Γ ( 2 2 q ) t q + ( ϕ ( 0 ) ) 2 Γ ( 2 q ) Γ ( 2 2 q ) t q t 1 q 2 t 1 q ϕ 2 ( 0 ) Γ ( 2 q ) t q t 1 q Γ ( 2 2 q ) 3.2 ϕ 2 ( 0 ) Γ ( 2 q ) Γ ( 2 2 q ) t q ) + ϕ 2 ( 0 ) Γ ( 2 q ) Γ ( 2 2 q ) t q = 0.2 ϕ 2 ( 0 ) Γ ( 2 q ) Γ ( 2 2 q ) t q ) < 0 .
Let t 1 . Then, we get
D ( 2 ) + V ( t , ϕ ( 0 ) , ϕ ) = t 1 q 2 ϕ ( 0 ) f ( t , ϕ ( 0 ) , ϕ ( 1 ) ) + ( ϕ ( 0 ) ) 2 Γ ( 2 q ) Γ ( 2 2 q ) t q = t 1 q 3 ϕ 2 ( 0 ) Γ ( 2 q ) Γ ( 2 2 q ) t q + ( ϕ ( 0 ) ) 2 Γ ( 2 q ) Γ ( 2 2 q ) t q < 0 .
Therefore, all the conditions of Theorem 3 are satisfied, and thus the zero solution of (39) and (40) is stable in time.

5. Conclusions

The nonlinear RL fractional differential equation is studied. The initial value problem is a subject that remains quite up-to-date (see, for example, the books [5,6]). Note the initial condition imposed to study fractional kinetic equations with RL fractional derivative. This point is critical in many physical situations, especially in astrophysical problems and the problem of anomalous subdiffusion ([23]). A good overview of the physical interpretation of initial conditions for fractional differential equations with Riemann–Liouville fractional derivatives was most clearly formulated by Diethelm [19] and it is detailed discussed in [7], where it is shown that initial conditions for RL fractional differential equations have physical meaning, and that the corresponding quantities can be obtained from measurements.
In this paper, some new definition for stability excluding a small interval around zero is defined and studied. These types of stability are called stability in time and generalize Mittag-Leffler stability in time for RL fractional differential equations. The definitions are deeply connected with the singularity at the initial time point. The stability properties of the zero solution are studied by Lyapunov functions. Two types of derivatives of Lyapunov functions: the RL fractional derivative when the argument of the Lyapunov function is a solution of the studied problem and the Dini fractional derivative among the studied problem.

Author Contributions

Conceptualization, R.A., S.H. and D.O.; methodology, R.A., S.H. and D.O.; validation, R.A., S.H. and D.O.; formal analysis, R.A., S.H. and D.O.; writing—original draft preparation, R.A., S.H. and D.O.; writing—review and editing, R.A., S.H. and D.O.; funding acquisition, S.H. All authors have read and agreed to the published version of the manuscript.

Funding

The research is supported by the Bulgarian National Science Fund under Project KP-06-N32/7.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Hilfer, R. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  2. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier Science B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  3. Petras, I.; Vinagre, B.M. Practical application of digital fractional-order controller to temperature control. Acta Montan. Slovaca 2002, 7, 131–137. [Google Scholar]
  4. Podlubny, I. Geometric and physical interpretation of fractional integration and fractional differentiation. Fract. Calc. Appl. Anal. 2002, 5, 367–386. [Google Scholar]
  5. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  6. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: Amsterdam, The Netherlands, 1993. [Google Scholar]
  7. Heymans, N.; Podlubny, I. Physical interpretation of initial conditions for fractional differential equations with Riemann-Liouville fractional derivatives. Rheol Acta 2006, 46, 765–771. [Google Scholar] [CrossRef] [Green Version]
  8. Agarwal, R.; Hristova, S.; O’Regan, D. Lyapunov functions and stability of Caputo fractional differential equations with delays. Differ. Eq. Dyn. Syst. 2018. [Google Scholar] [CrossRef]
  9. Baleanu, D.; Ranjbar, A.; Sadati, S.J.; Delavari, H.; Maraaba, T.A.; Gejji, V. Lyapunov-Krasovskii stability theorem for fractional systems with delay. Rom. J. Phys. 2011, 56, 636–643. [Google Scholar]
  10. Baleanu, D.; Sadati, S.J.; Ghaderi, R.; Ranjbar, A.; Maraaba, T.A.; Jarad, F. Razumikhin stability theorem for fractional systems with delay. Abstr. Appl. Anal. 2010, 2010, 124812. [Google Scholar] [CrossRef] [Green Version]
  11. Chen, B.; Chen, J. Razumikhin-type stability theorems for functional fractional-order differential systems and applications. Appl. Math. Comput. 2015, 254, 63–69. [Google Scholar] [CrossRef]
  12. Li, M.; Wang, J.R. Representation of solution of a Riemann-Liouville fractional differential equation with pure delay. Appl. Math. Lett. 2018, 85, 118–124. [Google Scholar] [CrossRef]
  13. Liu, S.; Wu, X.; Zhou, X.F.; Jiang, W. Asymptotical stability of Riemann-Liouville fractional nonlinear systems. Nonlinear Dyn. 2016, 2016 86, 65–71. [Google Scholar] [CrossRef]
  14. Liu, S.; Zhou, X.-F.; Li, X.; Jiang, W. Asymptotical stability of Riemann-Liouville fractional singular systems with multiple time-varying delays. Appl. Math. Lett. 2017, 65, 32–39. [Google Scholar] [CrossRef]
  15. Chartbupapan, W.; Bagdasar, O.; Mukdasai, K. A novel delay-dependent asymptotic stability conditions for differential and Riemann-Liouville fractional differential neutral systems with constant delays and nonlinear perturbation. Mathematics 2020, 8, 82. [Google Scholar] [CrossRef] [Green Version]
  16. Liu, S.; Wu, X.; Zhang, Y.-J.; Yang, R. Asymptotical stability of Riemann-Liouville fractional neutral systems. Appl. Math. Lett. 2017, 69, 168–173. [Google Scholar] [CrossRef]
  17. Altun, Y. Further results on the asymptotic stability of Riemann-Liouville fractional neutral systems with variable delays. Adv. Differ. Eq. 2019, 2019, 437. [Google Scholar] [CrossRef]
  18. Altun, Y.; Tunc, C. On the Asymptotic Stability of a Nonlinear Fractional-order System with Multiple Variable Delays. Appl. Appl. Math. 2020, 15, 458–468. [Google Scholar]
  19. Diethelm, K. The Analysis of Fractional Differential Equations; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  20. Li, C.; Qian, D.; Chen, Y.Q. On Riemann-Liouville and Caputo Derivatives. Discret. Dyn. Nat. Soc. 2011, 2011, 562494. [Google Scholar] [CrossRef] [Green Version]
  21. Devi, J.V.; Rae, F.A.M.; Drici, Z. Variational Lyapunov method for fractional differential equations. Comput. Math. Appl. 2012, 64, 2982–2989. [Google Scholar] [CrossRef] [Green Version]
  22. Sin, C.-X. Existence and uniqueness of global solutions of fractional functional differential equations with bounded delay. arXiv 2017, arXiv:1705.06046v1. [Google Scholar]
  23. Hilfer, R. Fractional diffusion based on Riemann-Liouville fractional derivatives. J. Phys. Chem. B 2000, 104, 3914. [Google Scholar] [CrossRef] [Green Version]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Agarwal, R.; Hristova, S.; O’Regan, D. Stability Concepts of Riemann-Liouville Fractional-Order Delay Nonlinear Systems. Mathematics 2021, 9, 435. https://doi.org/10.3390/math9040435

AMA Style

Agarwal R, Hristova S, O’Regan D. Stability Concepts of Riemann-Liouville Fractional-Order Delay Nonlinear Systems. Mathematics. 2021; 9(4):435. https://doi.org/10.3390/math9040435

Chicago/Turabian Style

Agarwal, Ravi, Snezhana Hristova, and Donal O’Regan. 2021. "Stability Concepts of Riemann-Liouville Fractional-Order Delay Nonlinear Systems" Mathematics 9, no. 4: 435. https://doi.org/10.3390/math9040435

APA Style

Agarwal, R., Hristova, S., & O’Regan, D. (2021). Stability Concepts of Riemann-Liouville Fractional-Order Delay Nonlinear Systems. Mathematics, 9(4), 435. https://doi.org/10.3390/math9040435

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