Next Article in Journal
Connectedness and Path Connectedness of Weak Efficient Solution Sets of Vector Optimization Problems via Nonlinear Scalarization Methods
Next Article in Special Issue
Impulsive Delayed Lasota–Wazewska Fractional Models: Global Stability of Integral Manifolds
Previous Article in Journal
Informal Norm in Hyperspace and Its Topological Structure
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On Dynamic Systems in the Frame of Singular Function Dependent Kernel Fractional Derivatives

1
Department of Mathematics and General Sciences, Prince Sultan University, P. O. Box 66833, Riyadh 11586, Saudi Arabia
2
Department of Medical Research, China Medical University, Taichung 40402, Taiwan
3
Department of Mathematics, University of Mohamed Boudiaf-PB 166, M’sila 28000, Algeria
4
Department of Mathematics, Çankaya University, 06790 Etimesgut, Ankara, Turkey
5
Laboratoire Lmdan, Dèpartement de Mathèmatiques de la Dècision, Universitè Cheikh Anta Diop de Dakar, Facultè des Sciences Economiques et Gestion, Dakar Fann BP 5683, Senegal
*
Authors to whom correspondence should be addressed.
Mathematics 2019, 7(10), 946; https://doi.org/10.3390/math7100946
Submission received: 15 August 2019 / Revised: 17 September 2019 / Accepted: 9 October 2019 / Published: 11 October 2019
(This article belongs to the Special Issue Stability Analysis of Fractional Systems)

Abstract

:
In this article, we discuss the existence and uniqueness theorem for differential equations in the frame of Caputo fractional derivatives with a singular function dependent kernel. We discuss the Mittag-Leffler bounds of these solutions. Using successive approximation, we find a formula for the solution of a special case. Then, using a modified Laplace transform and the Lyapunov direct method, we prove the Mittag-Leffler stability of the considered system.

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 ε C , R e ( ε ) > 0 has the following form
( a I ε f ) ( t ) = 1 Γ ( ε ) a t ( t - v ) ε - 1 f ( v ) d v .
The left Riemann–Liouville fractional derivative of order ε , where, ε C , R e ( ε ) 0 reads
( a D ε f ) ( t ) = ( d d t ) n ( a I n - ε f ) ( t ) , n = [ ε ] + 1 .
The left Caputo fractional derivative of order ε , R e ( ε ) 0 takes the form
( a D ε f ) ( t ) = ( a I n - ε f ( n ) ) ( t ) , n = [ ε ] + 1 .
The left Hadamard fractional integral of order ε C , R e ( ε ) > 0 is defined by
( a J ε f ) ( t ) = 1 Γ ( ε ) a t ( ln t - ln v ) ε - 1 f ( v ) v d v ,
while the left Hadamard fractional derivative of order ε C , R e ( ε ) 0 is given as
( a D ε f ) ( t ) = ( t d d t ) n ( a J n - ε f ) ( t ) , n = [ ε ] + 1 .
In [13,14,16], the left Caputo–Hadamard fractional derivatives of order ε C , R e ( ε ) 0 were defined as
( a C D ε f ) ( t ) = a D ε f ( u ) - i = 0 n - 1 δ i f ( a ) i ! ( ln v a ) i ( t ) , δ = t d d t
and in the space A C δ n [ a , b ] = { f : [ a , b ] C : δ n - 1 [ f ( x ) ] A C [ a , b ] } equivalently by
a C D ε f ( t ) = a J n - ε ( v d d v n f ) ( t ) , n = [ ε ] + 1 ,
and for a < b , c R and 1 p < , define the function space
X c p ( a , b ) = f : [ a , b ] R : f X c p = a b | t c f ( t ) | p d t t 1 / p < .
For p = , f X c p = ess sup a t b [ t c | f ( t ) | ] . The generalized left fractional integrals in the sense of how they are defined in [10] have the form
( a I ε , ρ f ) ( t ) = 1 Γ ( ε ) a t ( t ρ - v ρ ρ ) ε - 1 f ( v ) d v v 1 - ρ .
Remark 1.
If ρ = 1 the integrals in (8) become the Riemann–Liouville fractional integrals (1) and when ρ 0 (8) become the Hadamard fractional integrals (4).
The left generalized fractional derivatives of order ε > 0 are defined in [11] as
( a D ε , ρ f ) ( t ) = γ n ( a I n - ε , ρ f ) ( t ) = γ n Γ ( n - ε ) a t ( t ρ - v ρ ρ ) n - ε - 1 f ( v ) d u v 1 - ρ
where ρ > 0 and where γ = t 1 - ρ d d t .
Remark 2.
Once ρ = 1 the derivatives in (9) reduce to the Riemann–Liouville fractional derivatives (2). Whereas, taking the limit of (9) as ρ 0 one obtains the Hadamard fractional derivatives (5).
The Caputo modifications of the left generalized fractional derivatives were proposed in [15] to have the form
( a C D ε , ρ f ) ( t ) = ( a I n - ε , ρ γ n f ) ( t ) = 1 Γ ( n - ε ) a t ( t ρ - v ρ ρ ) n - ε - 1 γ n f ( v ) d v v 1 - ρ .
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 k - left-fractional integral operator is defined by
a I ε , k f ( x ) = 1 Γ ( ε ) a x k ( x ) - k ( t ) ε - 1 f ( t ) k ( t ) d t ,
where ε C , R e ( ε ) 0 .
Remark 4.
We should remark that if k ( x ) = x , the left Riemann–Liouville fractional integrals (1) are obtained. If k ( x ) = ln x , the Hadamard integrals (4) are obtained and at last when k ( x ) = x ρ ρ the generalized fractional integrals in the sense of Katugampola (8) are obtained.
Definition 2
([17]). The k - left-fractional derivative of order ε C , R e ( ε ) 0 in the Caputo setting is defined by
a C D ε , k f ( x ) = ( a I n - ε , k ) ( T k n f ) ( x ) = 1 Γ ( n - ε ) a x k ( x ) - k ( t ) n - ε - 1 T k n f ( t ) k ( t ) d t ,
where n = [ R e ( ε ) ] + 1 , and T g n f ( x ) = ( D k ( x ) ) n f ( x ) is the composition n times of the differential operator D k ( x ) .
Remark 5.
One should notice that if one takes k ( x ) = x , the left Caputo fractional derivatives (3) are obtained. If k ( x ) = ln x , the Caputo–Hadamard integrals (6) are obtained and finally if k ( x ) = x ρ ρ the generalized Caputo fractional in the sense of Katugampola (10) are obtained. In particular, if 0 < ε < 1 , then
a C D ε , k f ( x ) = ( a I 1 - ε , k ) T k f ( x ) = 1 Γ ( 1 - ε ) a x k ( x ) - k ( t ) - ε f ( t ) d t .
In the Riemann–Liouville setting we have [3]
a D ε , k f ( x ) = T k ( a I 1 - ε , k ) f ( x ) = T k Γ ( 1 - ε ) a x k ( x ) - k ( t ) - ε f ( t ) k ( t ) d t .
Lemma 1
([3]). Let R e ( β ) > 0 , R e ( ε ) > 0 and k ( t ) be an increasing and positive monotone function on ( a , b ] which has a continuous derivative k ( t ) on ( a , b ) . Then,
( a I ε , k ( k ( t ) - k ( a ) ) β - 1 ) ( x ) = Γ ( β ) Γ ( β + ε ) ( k ( x ) - k ( a ) ) β + ε - 1 .
Lemma 2
([3]). For a continuous real-valued function f and an increasing and positive monotone function k on ( a , b ] with continuous derivative k on ( a , b ) , we have
a I ε , k a I β , k f ( x ) = a I ε + β , k f ( x ) , ε > 0 , β > 0 .
Lemma 3
([17]). For ε > 0 and f and k as desired we have
a I ε , k a C D ε , k f ( x ) = f ( x ) - i = 0 n - 1 ( k ( x ) - k ( a ) ) i i ! T k i f ( a ) .
In particular, if 0 < ε < 1 then
a I ε , k a C D ε , k f ( x ) = f ( x ) - f ( a ) .
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
a C D ε , k x ( t ) = f ( t , x ( t ) ) , x ( a ) = x a , t ( a , b ] ,
where ε ( 0 , 1 ) , f : [ a , b ) × G , G an open subset of R or more generally of C , and
h ( t ) = f ( t , x ( t ) ) C γ , k [ a , b ] = { y : ( a , b ] R : ( k ( t ) - k ( a ) ) γ y ( t ) C [ a , b ] } , 0 γ < 1 .
The space C γ , k [ a , b ] is a Banach space when it is endowed by the norm
y = sup t | e - N ( k ( t ) - k ( a ) ) ( k ( t ) - k ( a ) ) γ y ( t ) | , N > 0 ,
which is equivalent to the norm y γ , k = sup t | ( k ( t ) - k ( a ) ) γ y ( t ) | . When γ = 0 we accept that C γ , ρ [ a , b ] = C [ a , b ] the space of continuous functions on [ a , b ] and when k ( t ) = t we accept C γ , k [ a , b ] = C γ [ a , b ] (see [3], p. 4).
Definition 3.
A function x ( t ) is said to be a solution of the initial value problem (18) if
1. 
( t , x ( t ) ) D , D = [ a , b ] × B , B = { x R : | x | L } G , L > 0
2. 
x ( t ) satisfies (18).
Theorem 1.
The fractional initial value problem (18) has a unique solution in the space
C γ , k ε , 0 [ a , b ] = { y ( t ) C [ a , b ] : a C D ε , k y ( t ) C γ , k [ a , b ] } ,
with 0 γ < 1 and γ ε , provided that
A N ε < 1 ,
and f satisfies the Lipschitzian condition
| f ( t , y 1 ) - f ( t , y 2 ) | A | y 1 - y 2 | , A > 0 .
Proof. 
First let us prove the existence of a unique solution y ( t ) in the space C [ a , b ] . Define the operator Ψ : C [ a , b ] C [ a , b ] by
( Ψ x ) ( t ) = x a + a I ε , k f ( t , x ( t ) ) .
where the space C [ a , b ] is endowed with the norm y C = sup t | e - N ( k ( t ) - k ( a ) ) y ( t ) | , which is equivalent to the sup norm. For any y 1 , y 2 B , by the help of the Lipschitz condition (21) and by writing k ( t ) - k ( a ) = k ( t ) - k ( s ) + k ( s ) - k ( a ) we have
| e - N ( k ( t ) - k ( a ) ) ( Ψ y 1 ( t ) - Ψ y 2 ( t ) ) |
A Γ ( ε ) a t e - N . K ( t , s ) K ( t , s ) ε - 1 k ( s ) d s y 1 - y 2 C .
Using the change of variable u = K ( t , s ) = ( k ( t ) - k ( s ) ) , it follows that
Ψ y 1 - Ψ y 2 A Γ ( ε ) 0 ( k ( t ) - k ( a ) ) u ε - 1 e - N u d u y 1 - y 2 C ,
and hence, by using the definition of Gamma function we have,
Ψ y 1 - Ψ y 2 C A N ε y 1 - y 2 C .
By assumption (20), the mapping Ψ is a contraction and hence by Banach fixed point theorem it has a unique fixed point x C [ t 0 , b ] . Moreover,
lim m T m x a - x C = 0 .
From the definition of Ψ , x has the form
x ( t ) = x a + a I ε , k f ( t , x ( t ) ) .
From Lemma 3 with n = 1 , 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 x ( t 0 ) = x 0 . By the help of the fact that for a continuous z ( t ) , we have
a C D ε , k z ( x ) = a D ε , k [ z ( t ) - z ( a ) ] ( x ) , a D ε , k a I ε , k z ( t ) = z ( t ) ,
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
a C D ε , k T m x a - a C D ε , k x A T m x a - x A ( b ρ - a ρ ) γ T m x a - x C .
From (26), we conclude that lim m a D ε , k T m x a - a D ε , k x = 0 . That is a D ε , k x C γ , k [ a , b ] and x C γ , k ε , 0 [ a , b ] . □
Consider the initial value problem
a C D ε , k x ( t ) = λ x ( t ) + f ( t ) , t > a , x ( a ) = x a ,
where ε ( 0 , 1 ) , λ R , f real-valued function and k is an increasing and positive monotone function on ( a , b ] which has a continuous derivative k on ( a , b ) . 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 by
x ( t ) = E ε λ ( k ( t ) - k ( a ) ) ε x 0 +
a t k ( t ) - k ( s ) ε - 1 E ε , ε λ k ( t ) - k ( s ) ε f ( s ) k ( s ) d s ,
where E ε ( . ) and E ε , ε ( . ) are the Mittag-Leffler functions of 1 parameter and 2 parameters, respectively [3].
Proof. 
Consider the successive approximation
4 x m ( t ) = x a + λ a I ε , k x m - 1 ( t ) + a I ε , k f ( t ) , m = 1 , 2 , . . . , x 0 ( t ) = x 0 .
Then, by Lemma 1, by writing x a = x a ( k ( t ) - k ( a ) ) 1 - 1 we have
x 1 ( t ) = x a + λ Γ ( 1 ) Γ ( ε + 1 ) ( k ( t ) - k ( a ) ) ε x a + a I ε , k f ( t ) .
Proceeding inductively and by making use of Lemma 1, we have
x m ( t ) = x a i = 0 m λ i ( k ( t ) - k ( a ) ) i ε Γ ( k ε + 1 ) + i = 1 m λ i - 1 a I i ε , k f ( t ) .
Then, we reach our claim by expanding t 0 I k ε , k in the second summation, shifting the index k, interchanging the order of the integral and summation, and letting m . □
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 k -Mittag-Leffler Stability

Consider the system
a C D ε , k x ( t ) = f ( t , x ) , t > a , x ( a ) = x a ,
where ε ( 0 , 1 ) , , and f : [ a , ) × G R n is piecewise continuous in t and locally Lipchitz in x on [ a , ) × G with Lipchitz constant A > 0 and G R n is a domain that contains the origin x = 0 . Then we have
Definition 4.
The constant x 0 is an equilibrium point of the system (33), if f ( t , x 0 ) = 0 and when a C D ε , k is replaced by a D ε , k if a D ε , k x 0 = f ( t , x 0 ) .
Definition 5.
(k-Mittag-Leffler stability) The solution x ( t ) of the system (33) is called k-Mittag-Leffler stable if
x ( t ) h ( x 0 ) E ε ( λ ( k ( t ) - k ( a ) ) ε ) c ,
where λ > 0 , c > 0 , h ( 0 ) = 0 , h ( x ) 0 and h ( x ) is locally Lipchitz on on a ball B of R n with Lipchitz constant h 0 .
Lemma 4.
For the continuous function f ( t , x ) in (33), for any ε > 0 and a certain norm of R n we have
a I ε , k f ( t , x ( t ) ) a I ε , k f ( t , x ( t ) ) .
Proof. 
The proof is direct and similar to the proof of Lemma 2 in [27]. □
Theorem 3.
The solution of system (33) satisfies
x ( t ) x a E ε ( A ( k ( t ) - k ( a ) ) ) ε ) .
Proof. 
By applying the fractional integral a I ε , k to both sides of (33), and making use of Lemma 3 with n = 1 we have
x ( t ) - x ( a ) = a I ε , k f ( t , x ( t ) ) .
If we apply the norm . to (37), use Lemma 4 and use the Lipchitz assumption we reach at
| x ( t ) - x ( a ) | x ( t ) - x ( a ) = a I ε , k f ( t , x ( t ) ) A a I ε , k x ( t ) .
Find a nonnegative function L ( t ) such that
x ( t ) - x ( a ) = A a I ε , k x ( t ) - L ( t ) .
By applying the following successive iteration in the frame of an existence uniqueness theorem similar to Theorem 1
y m ( t ) = y 0 + A a I ε , k y m - 1 ( t ) - L ( t ) , y 0 ( t ) = x a , y ( t ) = x ( t ) ,
and proceeding similar to what we have applied in Theorem 2, we conclude the solution representation
x ( t ) = x a E ε ( A ( k ( t ) - k ( a ) ) ε ) - a t i = 0 1 Γ ( ε i ) A k ( k ( t ) - k ( s ) ) i ε - 1 L ( s ) k ( s ) d s = x a E ε ( A ( k ( t ) - k ( a ) ) ε ) - a t T k E ε ( A ( k ( t ) - k ( s ) ) ε ) L ( s ) k ( s ) d s ,
where the operator T k = 1 k ( t ) d d t . Since the second term in (41) is nonnegative, we conclude that
x ( t ) x a E ε ( A ( k ( t ) - k ( a ) ) ) ε ) .
This completes the proof. □
Remark 7.
In Theorem 3, if ε = 1 , it follows that
x ( t ) x a e A ( k ( t ) - k ( a ) ) .
Now, we extend the Lyapunov direct method to the case of k - fractional-order systems.
Theorem 4.
Let x 0 be an equilibrium point for the system (33) and D R n be a domain containing x 0 . Let V ( t , x ( t ) ) : [ a , ) × D R be a continuously differentiable function and locally Lipschitz with respect to x such that
β 1 x ( t ) 1 c V ( t , x ( t ) ) β 2 x ( t ) d
a C D ε , k V ( t , x ( t ) ) - β 3 x ( t ) d ,
x D , ε ( 0 , 1 ) , β 1 ; β 2 , β 3 ,c and d are arbitrary positive constants. Then, x 0 is k-Mittag-Leffler stable. If the assumptions hold globally on R n , then x 0 is globally k-Mittag-Leffler stable.
Proof. 
From (42) and (43), we obtain that
a C D ε , k V ( t , x ( t ) ) - β 3 β 2 V ( t , x ( t ) ) .
It follows that there exists a nonnegative function L ( t ) such that
a C D ε , k V ( t , x ( t ) ) = - β 3 β 2 V ( t , x ( t ) ) - L ( t ) .
Now, applying Theorem 1 or Theorem 5.2. in [17] (which was proved by the modified Laplace transforms), we get
V ( t , x ( t ) ) = V ( a , x a ) E ε - β 3 β 2 ( k ( t ) - k ( a ) ) ε - a t k ( t ) - k ( s ) ε - 1 E ε , ε - β 3 β 2 k ( t ) - k ( s ) ε L ( s ) k ( s ) d s .
Thus, we have
V ( t , x ( t ) ) V ( a , x a ) E ε - β 3 β 2 ( k ( t ) - k ( a ) ) ε .
Utilizing condition (42), we reach
β 1 x ( t ) 1 c V ( a , x a ) E ε - β 3 β 2 ( k ( t ) - k ( a ) ) ε .
Therefore,
x ( t ) V ( a , x a ) β 1 E ε - β 3 β 2 ( k ( t ) - k ( a ) ) ε c .
Hence, the system is k-Mittag-Leffler stable. □
Using Remark 1 in [17] with n = 1 or ε ( 0 , 1 ) we have the following identity
a C D ε , k x ( t ) = a D ε , k x ( t ) - x ( a ) Γ ( 1 - ε ) ( k ( t ) - k ( a ) ) - ε .
As a result of (44), For ε ( 0 , 1 ) and x ( a ) 0 we have a C D ε , k x ( t ) a D ε , k x ( t ) 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 a C D ε , k x ( t ) by a D ε , k x ( t ) , we have that x 0 is globally k-Mittag-Leffler stable.
Remark 8.
If in the above, we take k ( t ) = t , we recover the Mittag-Leffler stability for the systems of Caputo and Riemann–Liouville fractional derivatives. If we take k ( t ) = ln t 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 α : R 0 R 0 is said to be positive definite when it satisfies the conditions α ( 0 ) = 0 , and α ( s ) > 0 for all s > 0 . A class K denotes a set of all increasing functions which are positive definite functions.
The class K denotes a set of all unbounded k functions. The continuous function β : R 0 × R 0 R 0 is said to be class KL if the function β ( . , t ) K for any t 0 , 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 V : R + × R n R , class K functions α 1 , α 2 , α and a class KL function β satisfying the following assumptions
α 1 x V ( t , x ) α 2 x .
a C D ε , k V ( t , x ) - α x .
If in addition the following assumption holds
α 1 - 1 2 β x a , k ( t ) - k ( a ) ε x a E ε - c k ( t ) - k ( a ) ε
where 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 KL 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
a C D ε , k V ( t , x ) - α α 2 - 1 V ( t , x ) .
From [39], there exists a class KL function β such that the following relationship is held
α 1 x V ( t , x ) β x a , k ( t ) - k ( a ) ε .
From which the following relationship follows
x α 1 - 1 2 β x a , k ( t ) - k ( a ) ε .
Using the k-Mittag-Leffler estimate (47), we have the following identity
x α 1 - 1 2 β x a , k ( t ) - k ( a ) ε x a E ε - k k ( t ) - k ( a ) ε .
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 KL 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],
a C D ε , k x T P x 2 x T P a C D ε , k x ,
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
a C D ε , k x 1 = - x 1 + x 2 a C D ε , k x 2 = - x 2
where x T = ( x 1 , x 2 ) R 2 . Let the quadratic Lyapunov candidate function defined by V ( t , x ) = 1 2 x 1 2 + x 2 2 . Combining Equation (52) and the generalized fractional derivative in Caputo sense of the Lyapunov function along the trajectories, we obtain the following relationships
a C D ε , k V ( t , x ) x T a C D ε , k x = - x 1 2 + x 1 x 2 - x 2 2 = - x 1 2 + 1 2 x 1 2 + 1 2 x 2 2 - x 2 2 = - 1 2 x 1 2 - 1 2 x 2 2 = - V ( t , x ) .
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
a C D ε , k x 1 = - x 1 + x 2 2 a C D ε , k x 2 = - x 1 x 2 - x 2
where x T = ( x 1 , x 2 ) R 2 . Let the quadratic Lyapunov candidate function defined by V ( t , x ) = 1 2 x 1 2 + x 2 2 = 1 2 x 2 . Combining Equation (52) and the generalized fractional derivative in the Caputo sense of the Lyapunov function along the trajectories, we obtain the following relationships:
a C D ε , k V ( t , x ) x T a C D ε , k x = - x 1 2 + x 1 x 2 2 - x 1 x 2 2 - x 2 2 = - x 1 2 - x 2 2 = - 2 V ( t , x ) = - 2 V ( t , x ) .
Recalling Theorem 6 and the ε -integration of Equation (55), we obtain the following inequality
x V ( x a ) E ε - 2 k ( t ) - k ( a ) ε .
Thus, the trivial solution of the k-fractional differential equation (54) is global k-Mittag-Leffler stable.

6. Conclusions

The fractional derivatives, in either Riemann–Liouville or Caputo settings, of function with respect to function, unify and extend Riemman–Liouville, Caputo, Hadamard, Caputo–Hadamard and other fractional derivatives. In this article, we proved the existence and uniqueness theorem for the generalized fractional derivative systems in newly defined function spaces. Then, we applied the successive approximation to obtain the solution representation for nonhomogeneous fractional differential systems with a constant coefficient. We employed these results to obtain Mittag-Leffler bounds for the solutions under Lipschitzity. The proved tools were used to prove k-Mittag-Leffler stability theorems of the Caputo and Riemann–Liouville types using a Lyapunov direct method. Alternatively, it is worth mentioning that the Laplace transforms introduced recently in [17] can be used to proceed in the proof of the main results instead of applying an operator method through successive approximations. Finally, we solved and discussed some examples. We should note that in the classical stability analysis, studying the global k-Mittag-Leffler stability using the analytical solution is not trivial and, in fact, almost impossible. An alternative issue is in using the Lyapunov direct method. The inconvenience of this proposed method is in finding the exact Lyapunov function for the given model. Another inconvenience of this method is in the calculation of the fractional derivative. To avoid these difficulties we propose using (52) which is fundamental in our studies.

Author Contributions

Conceptualization, T.A. and F.J.; methodology, F.J. and N.S.; formal analysis, F.M.; investigation, F.M.; resources, F.M.; writing—original draft preparation, T.A.; writing—review and editing, F.J. and N.S.; visualization, T.A.; supervision, N.S.

Funding

The first author would like to thank Prince Sultan University for funding this work through research group Nonlinear Analysis Methods in Applied Mathematics (NAMAM) group number RG-DES2017-01-17.

Acknowledgments

All authors would like to thank their institutions for the technical supports.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999. [Google Scholar]
  2. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach: Yverdon, Switzerland, 1993. [Google Scholar]
  3. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Application of Fractional Differential Equations; North Holland Mathematics Studies 204; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  4. Magin, R.L. Fractional Calculus in Bioengineering; Begell House Publishers: Redding, CT, USA, 2006. [Google Scholar]
  5. Hilfer, R. Applications of Fractional Calculus in Physics; Word Scientific: Singapore, 2000. [Google Scholar]
  6. Luchko, Y. Fractiona wave equations and damped waves. J. Math. Phys. 2013, 54, 031505. [Google Scholar] [CrossRef]
  7. Suzuki, J.L.; Zayernouri, M.; Bittencourt, M.L.; Karniadakis, G.E. Fractional-order uniaxial visco-elasto-plastic models for structural analysis. Comput. Meth. Appl. Mech. Eng. 2016, 308, 443–467. [Google Scholar] [CrossRef] [Green Version]
  8. Rakhshan, M.; Gupta, V.; Goodwine, B. On passivity of fractional order systems. SIAM J. Control Optim. 2019, 57, 1378–1389. [Google Scholar] [CrossRef]
  9. Ardjouni, A.; Djoudi, A. Existence and uniqueness of solutions for nonlinear implicit Caputo–Hadamard fractional differential equations with nonlocal conditions. Adv. Theor. Nonlinear Anal. Appl. 2019, 3, 46–52. [Google Scholar] [CrossRef]
  10. Katugampola, U.N. New approach to a generalized fractional integral. Appl. Math. Comput. 2011, 218, 860–865. [Google Scholar] [CrossRef] [Green Version]
  11. Katugampola, U.N. A new approach to generalized fractional derivatives. Bull. Math. Anal. Appl. 2014, 6, 1–15. [Google Scholar]
  12. Kilbas, A.A. Hadamard type fractional calculus. J. Korean Math. Soc. 2001, 38, 1191–1204. [Google Scholar]
  13. Jarad, F.; Abdeljawad, T.; Baleanu, D. Caputo-type modification of the Hadamard fractional derivatives. Adv. Differ. Equ. 2012, 2012, 142. [Google Scholar] [CrossRef] [Green Version]
  14. Gambo, Y.Y.; Jarad, F.; Baleanu, D.; Abdeljawad, T. On Caputo modification of the Hadamard fractional derivatives. Adv. Differ. Equ. 2014, 2014, 10. [Google Scholar] [CrossRef] [Green Version]
  15. Jarad, F.; Abdeljawad, T.; Baleanu, D. On the generalized fractional derivatives and their Caputo modification. J. Nonlinear Sci. Appl. 2017, 10, 2607–2619. [Google Scholar] [CrossRef] [Green Version]
  16. Adjabi, Y.; Jarad, F.; Baleanu, D.; Abdeljawad, T. On Cauchy problems with Caputo Hadamard fractional derivatives. J. Comput. Anal. Appl. 2016, 21, 661–681. [Google Scholar]
  17. Jarad, F.; Abdeljawad, T. Generalized fractional derivatives and Laplace transform. Discret. Contin. Dyn. Ser. S 2019, 709–722. [Google Scholar] [CrossRef]
  18. Meerschaert, M.M.; Zhang, Y.; Baeumer, B. Tempered anomlous diffusion in hetrogeneous systems. Geophys. Res. Lett. 2008, 35, L17403. [Google Scholar] [CrossRef]
  19. Atangana, A.; Baleanu, D. New fractional derivatives with non-local and non-singular kernel: Theory and application to heat transfer model. Therm. Sci. 2016, 20, 753–763. [Google Scholar] [CrossRef]
  20. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  21. Losada, J.; Nieto, J.J. Properties of a new fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 87–92. [Google Scholar]
  22. Abdeljawad, T.; Baleanu, D. Integration by parts and its applications of a new nonlocal fractional derivative with Mittag-Leffler nonsingular kernel. J. Nonlinear Sci. Appl. 2017, 10, 1098–1107. [Google Scholar] [CrossRef] [Green Version]
  23. Abdeljawad, T.; Baleanu, D. Monotonicity results for fractional difference operators with discrete exponential kernels. Adv. Differ. Equ. 2017, 2017, 78. [Google Scholar] [CrossRef]
  24. Abdeljawad, T.; Baleanu, D. On fractional derivatives with exponential kernel and their discrete versions. Rep. Math. Phys. 2017, 80, 11–27. [Google Scholar] [CrossRef]
  25. Jaishankar, A.; McKinley, G.H. Power-Law rheology in the bulk and at the interface: Quasi-properties and fractional constitutive equations. Proc. R. Soc. A 2013. [Google Scholar] [CrossRef]
  26. Li, Y.; Chen, Y.Q.; Podlubny, I. Stability of fractional-order nonlinear fractional systems: Lyapunov direct method and generalized Mittag-Leffler stability. Comput. Math. Appl. 2010, 59, 1810–1821. [Google Scholar] [CrossRef]
  27. Li, Y.; Chen, Y.Q.; Podlubny, I. Mittag-Leffler stability of fractional order nonlinear dynamic systems. Automatica 2009, 45, 1965–1969. [Google Scholar] [CrossRef]
  28. Jarad, F.; Abdeljawad, T.; Baleanu, D.; Biçen, K. On the Stability of Some Discrete Fractional Nonautonomous Systems. Abstr. Appl. Anal. 2012. [Google Scholar] [CrossRef]
  29. Jarad, F.; Abdeljawad, T.; Baleanu, D. Stability of q-fractional non-autonomous systems. Nonlinear Anal. Real World Appl. 2013, 14, 780–784. [Google Scholar] [CrossRef]
  30. Jarad, F.; Abdeljawad, T.; Baleanu, D.; Gündogdu, E. On the Mittag-Leffler stability of q-fractional nonlinear dynamical systems. Proc. Rom. Acad. Ser. A 2011, 12, 309–314. [Google Scholar]
  31. Baleanu, D.; Sadati, S.J.; Ghaderi, R.; Ranjbar, A.; Abdeljawad, T.; Jarad, F. Razumikhin stability theorem for fractional systems with delay. Abstr. Appl. Anal. 2010. [Google Scholar] [CrossRef]
  32. Sadati, S.J.; Baleanu, D.; Ranjbar, A.; Ghaderi, R.; Abdeljawad, T. Mittag-Leffler stability theorem for fractional nonlinear systems with delay. Abstr. Appl. Anal. 2010. [Google Scholar] [CrossRef]
  33. Baleanu, D.; Ranjbar, A.; Sadati, S.J.; Delvari, R.H.; Abdeljawad, T.; Gejji, V. Lyapunov-Krasovskii stability theorem for fractional systems with delay. Rom. J. Phys. 2011, 56, 636–643. [Google Scholar]
  34. Zhang, X. Some results of linear fractional order time delay systems. Appl. Math. Comput. 2008, 197, 407–411. [Google Scholar] [CrossRef]
  35. Deng, W.; Li, C.; Lü, J. Stability analysis of linear fractional differential system with multiple time delays. Nonlinear Dyn. 2007, 48, 409–416. [Google Scholar] [CrossRef]
  36. Momani, S.; Hadid, S. Lyapunov Stability solutions of fractional integrodifferential equations. Int. J. Math. Math. Sci. 2004, 47, 2503–2507. [Google Scholar] [CrossRef]
  37. Sene, N.; Srivastava, G. Generalized Mittag-Leffler Input Stability of the Fractional Differential Equations. Symmetry 2019, 11, 608. [Google Scholar] [CrossRef]
  38. Sene, N. Stability analysis of the generalized fractional differential equations with and without exogenous inputs. J. Nonlinear Sci. Appl. 2019, 12, 562–572. [Google Scholar] [CrossRef]
  39. Sene, N. Exponential form for Lyapunov function and stability analysis of the fractional differential equations. J. Math. Comput. Sci. 2018, 18, 388–397. [Google Scholar] [CrossRef] [Green Version]
  40. Aguila-Camacho, N.; Duarte-Mermoud, M.A.; Gallegos, J.A. Lyapunov functions for fractional order systems. Commun. Nonlinear Sci. Numer. Simul. 2014, 19, 2951–2957. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Abdeljawad, T.; Madjidi, F.; Jarad, F.; Sene, N. On Dynamic Systems in the Frame of Singular Function Dependent Kernel Fractional Derivatives. Mathematics 2019, 7, 946. https://doi.org/10.3390/math7100946

AMA Style

Abdeljawad T, Madjidi F, Jarad F, Sene N. On Dynamic Systems in the Frame of Singular Function Dependent Kernel Fractional Derivatives. Mathematics. 2019; 7(10):946. https://doi.org/10.3390/math7100946

Chicago/Turabian Style

Abdeljawad, Thabet, Fadila Madjidi, Fahd Jarad, and Ndolane Sene. 2019. "On Dynamic Systems in the Frame of Singular Function Dependent Kernel Fractional Derivatives" Mathematics 7, no. 10: 946. https://doi.org/10.3390/math7100946

APA Style

Abdeljawad, T., Madjidi, F., Jarad, F., & Sene, N. (2019). On Dynamic Systems in the Frame of Singular Function Dependent Kernel Fractional Derivatives. Mathematics, 7(10), 946. https://doi.org/10.3390/math7100946

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