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
where
is the gamma function.
Note that the notation is sometimes used.
- -
Riemann–Liouville fractional derivative of order
- -
Grünwald–Letnikov fractional derivative of order
where
Remark 1. If , then (see, Theorem 2.25 [19]), 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.
We will provide some well-known results from the literature (see, for example, [
5,
19]):
Proposition 1. For : Proposition 3 (Property 4 [
20]).
If the inequalities and hold, then . 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 of order only means that the RL fractional integral is monotonously increasing with respect to t and it does not imply that the function 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 . Suppose that for any , we have and for . Then, it follows that . Proposition 5. The initial value problemhas a unique solutionwhere 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 ()It is well known that , i.e., the solution of the above RL fractional differential equation, is where c is a real constant. Now, consider the initial condition where k is a real constant. However, . This illustrates that the initial condition of the type is not applicable for RL fractional equations (see, for example, [11]). Now, consider the initial condition where k is a real constant. Then, , i.e., the initial condition has a meaning for the RL fractional derivative with .
The practical definition of the initial condition of fractional differential equations with RL derivatives is based on the following result:
Lemma 1 ([
2]).
Let and , be a Lebesgue measurable function.- (a)
If there exists a.e. a limit , then there also exists a limit - (b)
If there exists a.e. a limit and if there exists the limit , then
3. Statement of the Problem
Consider the following nonlinear Riemann–Liouville fractional delay differential equation (RLFrDDE) of fractional order
:
with initial conditions
where
,
,
and
.
Remark 4. According to Lemma 1, the second equation in the initial conditions (3) could be replaced by the equality . We denote the solution of the initial value problem (IVP) for RLFrDDE (
2) and (
3) by
for
. 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 , we denote where is a norm in .
We will introduce the following conditions
Hypothesis 1 (H1).
The function is such that for any initial function , the corresponding IVP for RLFrDDE (2) and (3) has a solution ; Hypothesis 2 (H2). for all .
Remark 5. If and condition (H2) is satisfied, then, because of the equality , 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 there exist numbers and depending on ε such that for any initial functions , the corresponding solution of IVP (2) and (3) satisfies for ; - -
asymptotically stable if it is stable in time and additionally as ;
- -
generalized Mittag-Leffler stable in time if there exist positive numbers and and a locally Lipschitz function such that for any , there exists the solution of IVP (2) and (3) satisfies
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 problemhas a unique solution Therefore, the zero solution is generalized Mittag-Leffler stable with , and .
The zero solution is stable in time because for any there exist such that . At the same time, the zero solution is not stable (in the regular sense) because and cannot be satisfied for values of t close to 0.
For example, if then if then
Throughout the paper, we shall use the class
and
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 , , be a given interval, and be a given set. We will say that the function belongs to the class if is continuous on , and it is locally Lipschitzian with respect to its second argument. In our study, we will use the Razumikhin condition for the Lyapunov function
and any
:
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
used in the literature to study stability properties of solutions of fractional differential Equation (
2):
- -
RL fractional derivative—Let
be a solution of the IVP for the RLFrDDE (
2) and (
3). Then, we consider
- -
Dini fractional derivative—Let
. Then, consider (see [
8])
where
is defined by (
1) and
.
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., .
Example 2. Let where and . Apply (5) and obtain Special case 1. Let for . According to Proposition 1 with , we get and from (6), we get Special case 2. Let for . From Proposition 2, with , we obtain Special case 3. Let for . According to Proposition 2, we get from (6) Special case 4. Let where for all . For example, if , then , exists and Remark 6. Note that , for any solution of (2) and (3) because If , then the function for any solution of (2) and (3). The quadratic Lyapunov function is not from the set for any solution of (2) and (3), because The functions is not from the set for any solution 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 such that
- (i)
for any there exists such thatwhere ; - (ii)
there exists an increasing function such that for any function the inequality holds;
- (iii)
for any point such that for , the RL fractional derivative exists and the inequalityholds where 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
be an arbitrary number. According to condition (i), there exists
such that inequality (
9) holds for
.
Let be such that for , the inequality holds.
Consider the solution
of the IVP for RLFrDDE (
2) and (
3) with the initial function
.
From condition 2(ii) with
and
, we get
Therefore, there exists such that for .
Consider the function
and
. Therefore, there exists
such that
Define the function
for
. From (
11), it follows that
.
Note that the inequality (
13) holds for
. Assume inequality (
13) is not true for all
. Therefore, there exists a point
such that
Therefore,
. According to Proposition 4 with
, the inequality
holds. From Proposition 2, we get
and therefore,
Case 1. Let
. Then,
. From (
14), it follows that
, or
for
. According to condition 2(iii)
The inequality (
16) contradicts (
15).
Case 2. Let
. Then,
. From (
14), it follows that
, or
for
and the proof is similar to the one of Case 1.
From inequality (
13) and condition (i), it follows that
where
.
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 with satisfies condition (ii) with . Condition (i) is satisfied with if for .
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 such that conditions (i) and (ii) of Theorem 1 hold with , locally Lipschitz function and() for any point such that for , the RL fractional derivative exists and the inequalityholds where , 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
of the IVP for RLFrDDE (
2) and (
3) with the initial function
.
From condition (ii) of Theorem 1, similar to inequality (
11), we get
Consider the function .
Define the function
for
. From (
19), it follows that
.
Let
be an arbitrary number. We will prove that
We have . Therefore, there exists such that for , the inequality holds.
Assume inequality (
20) is not true for all
. Therefore, there exists a point
such that
Therefore,
. According to
and Proposition 4 with
, the inequality
holds. From
, we have,
Case 1. Let
. Then,
. Therefore,
for
. Let
. Then,
for
and, according to condition (iii*),
The inequality (
23) contradicts (
22).
Case 2. Let . Then, and for and the proof is similar to the one of Case 1.
Since
is an arbitrary number from inequality (
20), it follows that
Now, let
. According to condition (i), there exists
such that
for
. Then, from inequality (
24), it follows that
This proves the generalized Mittag-Leffler stability in time of the zero solution with , and (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 is a solution of the IVP for RLFrDDE (2) and (3) with where . - 2.
The function , , is such that:
- (i)
There exists an increasing function such that the inequality holds;
- (ii)
for any point such that for , the inequalityholds where is the Dini fractional derivative defined by (5) and .
Then, for
Remark 10. Let us, for simplicity, again consider the RL fractional differential equation without any delay with a solution (see Example 1). If we consider the quadratic Lyapunov function , then and is not satisfied. However, if then and 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 , where , , and be a solution of (2) and (3). Then, the followingholds, i.e., and condition 2 (i) of Lemma 2 is satisfied with . Proof. Define the function for . According to condition 2 (i), .
Let
where
and
be an arbitrary number. We will prove
For
from condition 2(i), we get
, i.e., the inequality (
27) is true.
Assume (
27) is not true. Therefore, there exist
such that
From Proposition 4, we have the inequality
. Then, applying Proposition 2 and
, we obtain
For any
and
, we let
From Remark 1 and Equation (
2), it follows that the function
satisfies for
, the equalities
and
Therefore,
or
with
as
. Then, for any
we obtain
Since
V is locally Lipschitzian in its second argument with a Lipschitz constant
,, we obtain
Substitute (
32) in (
31), divide both sides by
, take the limit as
, use
if
and we obtain for any
the inequality
Let
. Define the function
. From the choice of the point
, it follows that
,
and from inequalities (
26) and (
33) for
, we get
Now (
34) contradicts (
29). Therefore, inequality (
27) holds for an arbitrary
. Thus, the claim in our Lemma is true. □
Theorem 3. Let conditions (H1) and (H2) be satisfied and there exists a function such that and
- (i)
for any , there exists such thatwhere ; - (ii)
there exists an increasing function such that for any function the inequality holds;
- (iii)
for any function such that if for a point t we have for , then the inequality holds.
Then, the zero solution of (2) with the zero initial function is stable.
Proof. Let
be a positive number. According to condition (i), there exists
such that inequality (
9) holds for
.
There exists a positive number
such that
for
. Choose the function
such that
. Consider the solution
of the IVP for RLFrDDE (
2) and (
3) with initial function
and define the function
for
. From condition 2(ii), it follows that
.
Since
, there exists
such that
for
. Therefore,
According to Lemma 2 with
,
, and
, we obtain the following inequality
From inequalities (
37), (
38) and condition (i), it follows that
for
with
.
This proves the stability in time of the zero solution. □
Example 3. Consider the scalar RL fractional differential equationwith initial conditionswhere , and is defined by Consider the Lyapunov function . According to Remarks 8, this function satisfies condition (i) of Theorem 3 with and . It also satisfies condition (ii) with .
Now, let the function and the point be such that for , i.e., .
For , apply and Example 2 (Special case 1) and obtain 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.