1. Introduction and Preliminaries
Fractional calculus, being considered as the generalization of the calculus dealing with integration and differentiation of functions with integer order, has fetched the interest of many researchers working on different scopes of science and engineering [
1,
2,
3,
4,
5].
Due to the spatial non-locality of the operators with fractional orders, these operators have become significant tools that enable researchers to bring new aspects to the dynamics of non-local complex systems [
1,
2,
3,
4,
5,
6,
7,
8]. In addition, there has been an intensive interest in dealing with differential equations embodying derivatives with fractional order from many point of views including the qualitative, theoretical and numerical aspects [
1,
2,
3] and studying the existence and uniqueness of solutions of differential equations in the frame of the traditional fractional derivative has been tackled in many works (see [
9] and the references therein).
The generalization of the traditional fractional operators has been an important issue and it has been tackled by many researchers. The author in [
10,
11] proposed what was called the generalized fractional operators. These operators involve two parameters and are reduced to well known operators such as the Riemann–Liouville and the Hadamard fractional operators [
3,
12] when one of these parameters approaches some values. In addition, the fractional derivatives considered there were also modified in the Caputo settings so that they merge the Caputo fractional derivative and the Caputo–Hadamard one [
3,
13,
14,
15,
16]. Nevertheless, the fractional operators considered [
10,
11] are still special cases of the so called the fractional operators of a function with respect to another function or fractional operators with function dependent kernels [
2,
3]. The Caputo modification of the fractional derivatives of functions with respect to another function was presented in [
17]. In that work, the authors considered a modified Laplace transform that enables the scientists to deal with such operators. In addition, in the literature one can find more generalized types of fractional derivatives, among which the so called tempered fractional derivatives attracted many researchers. It is worth mentioning that the tempered fractional derivatives utilize an exponential tempering to the power-law kernel. For more information on these derivatives we refer to [
18] and the references cited there.
In spite of the fact that there are many types of traditional fractional operators, these operators, having singularities in their kernels, might not be the most ideal operator a researcher could choose for the sake of better understanding and modeling the real world phenomena under investigation. For this reason, some researchers considered and explored new types of fractional operator embodying non-singular kernels that still possess the non-locality property, which is the most fundamental property of the fractional operators [
19,
20,
21,
22,
23,
24].
The main feature of the fractional derivatives is that the power law kernel codes the flexibility of an infinite number of relaxation times. This is one of the main reasons that fractional differential equations provide good experimental fits to a broad range of time scales [
25]. The singularities from the power-law kernels can be computationally dealt with using Lubich’s corrections and non-uniform grids.
There has been a large number of articles from those that discussed the extension of qualitative properties of differential equations with integer order derivatives, to those including fractional derivatives. Among these qualitative properties, stability plays a significant role [
26,
27,
28,
29,
30,
31,
32,
33,
34,
35,
36]. In this article, we study some qualitative properties of dynamic systems in the frame of fractional derivatives involving a function dependent kernel. We discuss the existence and uniqueness of a nonlinear Cauchy problem and obtain the solution in the linear case. Then, we investigate the Mittag-Leffler stability of these systems.
Before we start presenting the main results, we state some definitions, lemmas and properties from the literature.
The left Riemann–Liouville fractional integral of order , where has the following form
The left Riemann–Liouville fractional derivative of order , where, reads
The left Caputo fractional derivative of order takes the form
The left Hadamard fractional integral of order
is defined by
while the left Hadamard fractional derivative of order
is given as
In [
13,
14,
16], the left Caputo–Hadamard fractional derivatives of order
were defined as
and in the space
equivalently by
and for
and
, define the function space
For
,
. The generalized left fractional integrals in the sense of how they are defined in [
10] have the form
Remark 1. If the integrals in (8) become the Riemann–Liouville fractional integrals (1) and when (8) become the Hadamard fractional integrals (4). The left generalized fractional derivatives of order
are defined in [
11] as
where
and where
.
Remark 2. Once the derivatives in (9) reduce to the Riemann–Liouville fractional derivatives (2). Whereas, taking the limit of (9) as one obtains the Hadamard fractional derivatives (5). The Caputo modifications of the left generalized fractional derivatives were proposed in [
15] to have the form
Remark 3. Replacing ρ by 1 in (10) the Caputo fractional derivatives (3) are obtained and approaching ρ to 0 in (10) gives the Caputo–Hadamard fractional derivatives (6). Definition 1 ([
3]).
The left-fractional integral operator is defined bywhere Remark 4. We should remark that if , the left Riemann–Liouville fractional integrals (1) are obtained. If , the Hadamard integrals (4) are obtained and at last when the generalized fractional integrals in the sense of Katugampola (8) are obtained. Definition 2 ([
17]).
The left-fractional derivative of order in the Caputo setting is defined bywhere , and is the composition n times of the differential operator . Remark 5. One should notice that if one takes , the left Caputo fractional derivatives (3) are obtained. If , the Caputo–Hadamard integrals (6) are obtained and finally if the generalized Caputo fractional in the sense of Katugampola (10) are obtained. In particular, if , then In the Riemann–Liouville setting we have [3] Lemma 1 ([
3]).
Let and be an increasing and positive monotone function on which has a continuous derivative on . Then, Lemma 2 ([
3]).
For a continuous real-valued function f and an increasing and positive monotone function k on with continuous derivative on , we have Lemma 3 ([
17]).
For and f and k as desired we haveIn particular, if then This article is organized as follows: The second section discusses the existence uniqueness of solutions to a certain type of differential equations involving a fractional derivative with dependence on a kernel function and taking the advantage of successive approximation, we find the solution of a special case of the system under consideration. In third section, an upper bound of the solutions is obtained and the Mittag-Leffler stability of the considered system is debated. In the last section, conclusions and discussions are presented.
2. Existence and Uniqueness Theorems
Consider the following fractional nonlinear Cauchy problem
where
,
,
G an open subset of
or more generally of
, and
The space
is a Banach space when it is endowed by the norm
which is equivalent to the norm
. When
we accept that
the space of continuous functions on
and when
we accept
(see [
3], p. 4).
Definition 3. A function is said to be a solution of the initial value problem (18) if
- 1.
- 2.
satisfies (18).
Theorem 1. The fractional initial value problem (18) has a unique solution in the spacewith and , provided thatand f satisfies the Lipschitzian condition Proof. First let us prove the existence of a unique solution
in the space
. Define the operator
by
where the space
is endowed with the norm
, which is equivalent to the sup norm. For any
, by the help of the Lipschitz condition (21) and by writing
we have
Using the change of variable
, it follows that
and hence, by using the definition of Gamma function we have,
By assumption (20), the mapping
is a contraction and hence by Banach fixed point theorem it has a unique fixed point
. Moreover,
From the definition of
,
x has the form
From Lemma 3 with
, it is clear that if
x satisfies the initial value problem (18), then it has the representation (27). Conversely, if
x has the representation (27), then clearly
. By the help of the fact that for a continuous
, we have
x will satisfy the equation (18). Hence,
x has the representation (27) if and only if it satisfies the initial value problem (18). Finally, if
denotes the norm defined in (19) then we have
From (26), we conclude that . That is and . □
Consider the initial value problem
where
,
,
f real-valued function and
k is an increasing and positive monotone function on
which has a continuous derivative
on
. Throughout this article
k will be assumed to have the previously mentioned properties.
Theorem 2. The solution of the Caputo initial value problem (28) is given bywhere and are the Mittag-Leffler functions of 1 parameter and 2 parameters, respectively [3]. Proof. Consider the successive approximation
Then, by Lemma 1, by writing
we have
Proceeding inductively and by making use of Lemma 1, we have
Then, we reach our claim by expanding in the second summation, shifting the index k, interchanging the order of the integral and summation, and letting . □
Remark 6. The solution (29) coincides with the solution given in Theorem 5.2. in [17]. 3. Mittag-Leffler Bounds for Fractional Non-Autonomous Systems and -Mittag-Leffler Stability
Consider the system
where
, and
is piecewise continuous in
t and locally Lipchitz in
x on
with Lipchitz constant
and
is a domain that contains the origin
. Then we have
Definition 4. The constant is an equilibrium point of the system (33), if and when is replaced by if .
Definition 5. (k-Mittag-Leffler stability) The solution of the system (33) is called k-Mittag-Leffler stable ifwhere , , and is locally Lipchitz on on a ball B of with Lipchitz constant . Lemma 4. For the continuous function in (33), for any and a certain norm of we have Proof. The proof is direct and similar to the proof of Lemma 2 in [
27]. □
Theorem 3. The solution of system (33) satisfies Proof. By applying the fractional integral
to both sides of (33), and making use of Lemma 3 with
we have
If we apply the norm
to (37), use Lemma 4 and use the Lipchitz assumption we reach at
Find a nonnegative function
such that
By applying the following successive iteration in the frame of an existence uniqueness theorem similar to Theorem 1
and proceeding similar to what we have applied in Theorem 2, we conclude the solution representation
where the operator
. Since the second term in (41) is nonnegative, we conclude that
This completes the proof. □
Remark 7. In Theorem 3, if , it follows that Now, we extend the Lyapunov direct method to the case of fractional-order systems.
Theorem 4. Let be an equilibrium point for the system (33) and be a domain containing . Let be a continuously differentiable function and locally Lipschitz with respect to x such that, ,c and d are arbitrary positive constants. Then, is k-Mittag-Leffler stable. If the assumptions hold globally on , then is globally k-Mittag-Leffler stable. Proof. From (42) and (43), we obtain that
It follows that there exists a nonnegative function
such that
Now, applying Theorem 1 or Theorem 5.2. in [
17] (which was proved by the modified Laplace transforms), we get
Utilizing condition (42), we reach
Hence, the system is k-Mittag-Leffler stable. □
Using Remark 1 in [
17] with
or
we have the following identity
As a result of (44), For and we have and hence we can state the following stability theorem for the Riemann–Liouville generalized type fractional derivatives.
Theorem 5. Under the assumptions of Theorem 4, if we replace by , we have that is globally k-Mittag-Leffler stable.
Remark 8. If in the above, we take , we recover the Mittag-Leffler stability for the systems of Caputo and Riemann–Liouville fractional derivatives. If we take then the sufficient conditions for the Mittag-Leffler stability of Caputo–Hadamard and Hadamard fractional systems are determined by Theorems 4 and 6, respectively.
4. The Lyapunov Direct Method in the Term of Comparison Functions
In this section, we analyze the
k-Mittag-Leffler stability using the Lyapunov direct method. The Lyapunov direct method presented in this section is done in the term of the comparison functions. The comparison functions play an essential role in the stability analysis. Before starting this section, we introduce some definitions related to the comparison functions [
37,
38].
A continuous function is said to be positive definite when it satisfies the conditions , and for all . A class denotes a set of all increasing functions which are positive definite functions.
The class denotes a set of all unbounded k functions. The continuous function is said to be class if the function for any , is non-increasing and tends to zero as its arguments tend to infinity. For the characterization of the global k-Mittag-Leffler stability using the Lyapunov direct method, we introduce the following theorems which generalizes Theorem 4.
Theorem 6. Suppose there exists a positive continuous function , class functions and a class function β satisfying the following assumptions If in addition the following assumption holdswhere c is a positive constant. Then, the trivial solution of the k-fractional differential equation (33) is globally k-Mittag-Leffler stable. When the condition (45) holds the function V is said to be radially unbounded. In other words, the function V is said to be a Lyapunov candidate function when the assumption (45) holds. The assumption (47) represents the Mittag-Leffler estimate of the class function and is essential to study the k-Mittag-Leffler stability.
Proof. The first assumption (45) implies the function
V is a Lyapunov candidate function and combining it with the second assumption (46), we have the following relationship
From [
39], there exists a class
function
such that the following relationship is held
From which the following relationship follows
Using the
k-Mittag-Leffler estimate (47), we have the following identity
From which it follows the global k-Mittag-Leffler stability of the trivial solution of the fractional differential equation (33). □
When the condition (45) holds the function V is said to be radially unbounded. In other words, the function V is said to be a Lyapunov candidate function when the assumption (45) holds. The assumption (47) represents the Mittag-Leffler estimate of the class function and is essential to study the k-Mittag-Leffler stability.
Remark 9. The use of the Lyapunov direct method in the context of fractional order derivative is not trivial due to the fact that the determination of the generalized fractional derivative in the Caputo sense of a given function is not always possible to obtain. In many cases, we suppose the Lyapunov function for the fractional differential equation is quadratic. For these reasons, we propose the following assumption which is already proved in the literature in [40],where P represents a positive definite matrix. The dimension of the matrix P depends on the dimension of the fractional differential equation (33). 5. Examples
In this section, we provided some examples to illustrate Theorems 4 and 6.
Firstly, we illustrate Theorem 4 by the following example. Let the k-fractional differential equation described by the generalized fractional derivative in the Caputo sense given by
where
. Let the quadratic Lyapunov candidate function defined by
. Combining Equation (52) and the generalized fractional derivative in Caputo sense of the Lyapunov function along the trajectories, we obtain the following relationships
It follows from Theorem 4, the trivial solution of the fractional differential Equation (53) is global k-Mittag-Leffler stable.
Secondly, we illustrate Theorem 6 by the following example. Let the fractional differential equation described by the generalized fractional derivative in the Caputo sense be given by
where
. Let the quadratic Lyapunov candidate function defined by
. Combining Equation (52) and the generalized fractional derivative in the Caputo sense of the Lyapunov function along the trajectories, we obtain the following relationships:
Recalling Theorem 6 and the
-integration of Equation (55), we obtain the following inequality
Thus, the trivial solution of the k-fractional differential equation (54) is global k-Mittag-Leffler stable.