Next Article in Journal
Function Reconstruction from Reflection Symmetric Radon Data
Previous Article in Journal
Forecasting Model for Stock Market Based on Probabilistic Linguistic Logical Relationship and Distance Measurement
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hyers–Ulam Stability and Existence of Solutions to the Generalized Liouville–Caputo Fractional Differential Equations

1
College of Science, Guizhou Institute of Technology, Guiyang 550025, China
2
Department of Mathematics, Guizhou University, Guiyang 550025, China
3
Department of Mathematical Analysis and Numerical Mathematics, Faculty of Mathematics, Physics and Informatics, Comenius University in Bratislava, Mlynská Dolina, 842 48 Bratislava, Slovakia
4
Mathematical Institute, Slovak Academy of Sciences, Štefánikova 49, 814 73 Bratislava, Slovakia
5
School of Mathematical Sciences, Qufu Normal University, Qufu 273165, China
*
Author to whom correspondence should be addressed.
Submission received: 4 April 2020 / Revised: 20 April 2020 / Accepted: 6 May 2020 / Published: 4 June 2020

Abstract

:
The aim of this paper is to study the stability of generalized Liouville–Caputo fractional differential equations in Hyers–Ulam sense. We show that three types of the generalized linear Liouville–Caputo fractional differential equations are Hyers–Ulam stable by a ρ -Laplace transform method. We establish existence and uniqueness of solutions to the Cauchy problem for the corresponding nonlinear equations with the help of fixed point theorems.

1. Introduction

Because fractional calculus has a good global correlation performance to reflect the historical dependence process of the development of system functions, and can also describe the attributes of the dynamic system itself, it becomes a powerful mathematical tool to describe some complex movements, irregular phenomena, memory features, and other aspects. Fractional calculus theory was widely used by mathematicians as well as chemists, engineers, economists, biologists, and physicists (see [1,2,3,4,5]). In 1876, Riemann proposed the definition of the Riemann–Liouville derivative. Caputo first proposed another definition of fractional derivative via a modified Riemann–Liouville fractional integral at the beginning of the 20th century, namely a Caputo fractional derivative. Caputo and Fabrizio [6] introduced a new nonlocal derivative without a singular kernel and obtained the new Caputo–Fabrizio fractional derivative of order 0 < α < 1 . Theoretical research and application of Caputo–Fabrizio fractional can be referred to [7,8,9,10,11,12,13]. Butzer et al. [14,15,16,17,18] study properties of the Hadamard fractional integral and the derivative. In [19,20], Katugampola introduced a new fractional integral and fractional derivative, which generalizes the Riemann–Liouville and the Hadamard integrals and derivative into a single form, respectively.
Hyers–Ulam stability has been one of the most active research topics in differential equations, and obtained a series of results (see [21,22,23,24,25,26,27,28,29,30]). Recently, Alqifiary et al. [22] obtained generalized Hyers–Ulam stability of linear differential equations. Razaei et al. [31] proved that the Hyers–Ulam stability of linear differential equations. Wang et al. [32] proved that two types of fractional linear differential equations are Hyers–Ulam stable. Shen et al. [33] deal with the Ulam stability of linear fractional differential equations with constant coefficients. Liu et al. [34] proved the Hyers–Ulam stability of linear Caputo–Fabrizio fractional differential equations. Liu et al. [35] studied the Hyers–Ulam stability of linear Caputo–Fabrizio fractional differential equations with the Mittag–Leffler kernel. Laplace transform method is used to deal with linear equations and fixed point approach and Gronwall inequality are used to deal with nonlinear equations.
For some differential equations describing physical models and practical problems, it is very difficult to find their exact solutions and the method of finding its exact solution (if they exist) is also very complicated. In order to construct explicit solutions to differential equations with constant coefficients and in the frame of Riemann–Liouville, Caputo and Riez fractional derivatives, integral transforms including Laplace, Mellin, and Fourier were found to be strong tools. One of the main difficulty is to find some appropriate transformations in order to find analytic solutions to some classes of fractional differential equations. In order to extend the possibility of working in a large class of functions, Jarad et al. [36] present a modified Laplace transform that it call ρ -Laplace transform, study its properties, and prove its own convolution theorem.
Motivated by [36], we apply the ρ -Laplace transform method to study the Hyers–Ulam stability of the following linear differential equations:
( D c α , ρ f ) ( t ) = g ( t ) , 0 < α < 1 , ρ > 0 ,
and
( D c α , ρ f ) ( t ) λ f ( t ) = g ( t ) , 0 < α < 1 , ρ > 0 ,
and
( D c α , ρ f ) ( t ) λ ( D c β , ρ f ) ( t ) = g ( t ) , 0 < α < 1 , ρ > 0 ,
where D c α , ρ f denotes the left generalized α order Liouville–Caputo fractional derivative for f with the parameter ρ (see Definition 2).
Next, we study to Cauchy problem for nonlinear equations as follows:
( D c α , ρ f ) ( t ) = g ( t , f ( t ) ) , 0 < α < 1 , f ( 0 ) = f 0 ,
and show the existence and uniqueness of solutions via Banach fixed point theorem and Schaefer’s fixed point theorem and obtain the generalized Hyers–Ulam-Rassias stability via an extended Gronwall’s inequality.

2. Preliminaries

Let C ( I , R ) be the Banach space of all continuous functions from I into R with the norm y C : = sup { | y ( x ) | : x I } .
Definition 1. 
(see [35], Definition 2.1) Let α > 0 , t , β R . E α ( t ) : = n = 0 t n Γ ( α n + 1 ) is called the standard Mittag–Leffler function. E α , β ( t ) : = n = 0 t n Γ ( α n + β ) is called two-parameter Mittag–Leffler function.
Definition 2. 
(see [37], Definition 5) Let f : [ 0 , + ) R . The Liouville–Caputo generalized derivative of the function f is expressed in the form
( D c α , ρ f ) ( t ) = 1 Γ ( 1 α ) 0 t ( t ρ s ρ ρ ) α f ( s ) d s s 1 ρ , t > 0 ,
where the order α ( 0 , 1 ) , ρ > 0 , and Γ ( · ) is the Gamma function.
Definition 3. 
(see [36]) Let 0 < α < 1 , ρ > 0 . The generalized left fractional integrals of the function f is expressed in the form
a ( I α , ρ f ) ( t ) = 1 Γ ( α ) a t ( t ρ s ρ ρ ) α 1 f ( s ) ) d s s 1 ρ , t > a , ρ > 0 ,
Definition 4. 
(see [36], Corollary 3.3) Let α ( 0 , 1 ) , ρ > 0 . The ρ–Laplace transform of the function of the generalized fractional derivative in the Liouville–Caputo sense is expressed in the following form:
L ρ { ( D c α , ρ f ) ( t ) } ( s ) = s α L ρ { f ( t ) } ( s ) s α 1 f 0 , ρ > 0 , f ( 0 ) = f 0 .
The ρ–Laplace transform of the function f is given in the form
L ρ { f ( t ) } ( s ) = 0 e s t ρ ρ f ( t ) d t t 1 ρ , ρ > 0 ,
and
L ρ { f ( t ) } ( s ) = L { f ( ( ρ t ) 1 ρ ) } ( s ) , ρ > 0 ,
where L { f } is the usual Laplace transform of f.
Definition 5. 
(see [36], Definition 2.9) Let f and g be two functions which are piecewise continuous at each interval [ 0 , T ] and of exponential order. The ρ convolution of f and g is given by
( f ρ g ) ( t ) = 0 t f ( ( t ρ s ρ ) 1 ρ ) g ( s ) d s s 1 ρ , ρ > 0 .
Definition 6. 
(see [36], Definition 2.11) Let f and g be two functions which are piecewise continuous at each interval [ 0 , T ] and of exponential order e c t ρ ρ . Then,
f ρ g = g ρ f , ρ > 0 ,
and
L ρ { f ρ g } ( s ) = L ρ { f } L ρ { g } ( s ) .
Definition 7. 
Let R e ( α ) > 0 , ρ > 0 , and | λ s α | < 1 .
(i) 
L ρ { 1 } ( s ) = 1 s , s > 0 .
(ii) 
L ρ { e λ t ρ ρ } ( s ) = 1 s λ .
Proof. 
It is easy to check the following facts:
(i) 
L ρ { 1 } ( s ) = 0 e s t ρ ρ d t t 1 ρ = 1 s 0 e s t ρ ρ d ( s t ρ ρ ) = 1 s .
(ii) 
L ρ { e λ t ρ ρ } ( s ) = 0 e λ t ρ ρ e s t ρ ρ d t t 1 ρ = 1 s λ .
The proof is complete. □
Definition 8. 
(see [36], Lemma 3.4) Let R e ( α ) > 0 , ρ > 0 , and | λ s α | < 1 .
(iii) 
L ρ { E α ( λ ( t ρ ρ ) α ) } ( s ) = s α 1 s α λ .
(iv) 
L ρ { ( t ρ ρ ) α 1 ) E α , α ( λ ( t ρ ρ ) α ) } ( s ) = 1 s α λ .
From Definition 8, we derive the following result.
Definition 9. 
Let R e ( α ) > 0 , ρ > 0 , and | λ s α | < 1 . Then,
L ρ { ( t ρ ρ ) β 1 ) E α , β ( λ ( t ρ ρ ) α ) } ( s ) = s α β s α λ .
Proof. 
One can see
L ρ { ( t ρ ρ ) β 1 ) E α , β ( λ ( t ρ ρ ) α ) } ( s ) = L ρ { ( t ρ ρ ) β 1 k = 0 ( λ ( t ρ ρ ) α ) k Γ ( α k + β ) } ( s ) = L ρ { k = 0 λ k t α ρ k + ρ β ρ Γ ( α k + β ) ρ α k + β 1 } ( s ) = k = 0 λ k Γ ( α k + β ) ρ α k + β 1 L ρ { t α ρ k + ρ β ρ } ( s ) = k = 0 λ k Γ ( α k + β ) ρ α k + β 1 ρ α k + β 1 Γ ( 1 + k α + β 1 ) s 1 + k α + β 1 = k = 0 λ k Γ ( α k + β ) Γ ( k α + β ) s k α + β = k = 0 λ k s k α + β = 1 s β k = 0 ( λ s α ) k = s α β s α λ .
The proof is finished. □
Motivated by [34,38], we introduce the following definitions.
Definition 10.
The fractional differential equation Ψ ( g , f , ( D c α 1 , ρ f ) , , ( D c α n , ρ f ) ) = 0 has Hyers–Ulam stability if there exists a real number K > 0 , such that for a given ϵ > 0 and for each solution f C ( [ 0 , T ] , R ) of the inequality
Ψ ( g , f , ( D c α 1 , ρ f ) , , ( D c α n , ρ f ) ) ϵ
there exists a solution f a C ( [ 0 , T ] , R ) of the differential equation such that | f ( t ) f a ( t ) | K ϵ .
Definition 11.
The fractional differential equation Ψ ( g , f , ( D c α 1 , ρ f ) , , ( D c α n , ρ f ) ) = 0 has generalized Ulam–Hyers–Rassias stable with respect to G C ( [ 0 , T ] , R ) if there exists a real numbr c G > 0 , such that for each solution f C ( [ 0 , T ] , R ) of the inequality
Ψ ( g , f , ( D c α 1 , ρ f ) , , ( D c α n , ρ f ) ) G ( t ) , t [ 0 , T ]
there exists a solution f a C ( [ 0 , T ] , R ) of the differential equation such that | f ( t ) f a ( t ) | c G G ( t ) .
Definition 12. 
([39], Corollary 2.4) Let α , ρ > 0 , x ( t ) , a ( t ) be nonnegative functions and b ( t ) be nonnegative and nondecreasing function for t [ t 0 , T ) , T > 0 , b ( t ) M , where M is a constant. If
x ( t ) a ( t ) + b ( t ) t 0 t ( t ρ τ ρ ρ ) α 1 x ( τ ) d τ τ 1 ρ ,
then
x ( t ) a ( t ) E a b ( t ) Γ ( α ) ( t ρ t 0 ρ ρ ) α , t [ t 0 , T ) .

3. Hyers–Ulam Stability for Linear Problems

Definition 13. 
Let 0 < α < 1 , ρ > 0 and g ( t ) be a given real continuous function on [ 0 , ) . If a function f : [ 0 , ) R satisfies the following inequality
| ( D c α , ρ f ) ( t ) g ( t ) | ε , f ( 0 ) = f 0 R ,
for each t 0 and ε > 0 , then there exists a solution f a : [ 0 , ) R of (1) such that
| f ( t ) f a ( t ) | t ρ α Γ ( α + 1 ) ρ α ε .
Proof. 
Let
F ( t ) = ( D c α , ρ f ) ( t ) g ( t ) , t 0 .
Taking the ρ -Laplace transform of (9) via Definition 4, we have
L ρ { F ( t ) } ( s ) = L ρ { ( D c α , ρ f ) ( t ) g ( t ) } ( s ) = L ρ { ( D c α , ρ f ) ( t ) } ( s ) L ρ { g ( t ) } ( s ) = s α L ρ { f ( t ) } ( s ) s α 1 f 0 L ρ { g ( t ) } ( s ) ,
where L ρ { F ( · ) } denotes the ρ -Laplace transform of the function F. From (10), one has
L ρ { f ( t ) } ( s ) = 1 s f 0 + 1 s α L ρ { g ( t ) } ( s ) + 1 s α L ρ { F ( t ) } ( s ) = f 0 L ρ { 1 } ( s ) + 1 Γ ( α ) L ρ { ( t ρ ρ ) α 1 ρ g ( t ) } ( s ) + 1 Γ ( α ) L ρ { ( t ρ ρ ) α 1 ρ F ( t ) } ( s ) .
Set
f a ( t ) = f 0 + 1 Γ ( α ) ( t ρ ρ ) α 1 ρ g ( t ) = f 0 + 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 g ( s ) d s s 1 ρ .
Taking the ρ –Laplace transform of (11), one has
L ρ { f a ( t ) } ( s ) = 1 s f 0 + 1 Γ ( α ) L ρ 0 t ( t ρ s ρ ρ ) α 1 g ( s ) d s s 1 ρ ( s ) = 1 s f 0 + 1 Γ ( α ) L ρ ( t ρ ρ ) α 1 ρ g ( t ) ( s ) = 1 s f 0 + L ρ 1 Γ ( α ) ( t ρ ρ ) α 1 ( s ) · L ρ { g ( t ) } ( s ) = 1 s f 0 + 1 s α · L ρ { g ( t ) } ( s ) .
Note that
L ρ { ( D c α , ρ f a ) ( t ) } ( s ) = s α L ρ { f a ( t ) } ( s ) s α 1 f 0 = s α 1 f 0 + L ρ { g ( t ) } ( s ) s α 1 f 0 = L ρ { g ( t ) } ( s ) ,
which yields that f a ( · ) is a solution of Equation (1), since, according to the one-to-one transformation of L in (6), we can get that L ρ is the one-to-one transformation.
From (11) and (12), we have
L ρ { f ( t ) f a ( t ) } ( s ) = 1 s f 0 + 1 s α L ρ { g ( t ) } ( s ) + 1 s α L ρ { F ( t ) } ( s ) 1 s f 0 1 s α · L ρ { g ( t ) } ( s ) = 1 s α L ρ { F ( t ) } ( s ) = L ρ 1 Γ ( α ) ( t ρ ρ ) α 1 ( s ) · L ρ { F ( t ) } ( s ) = L ρ 1 Γ ( α ) ( t ρ ρ ) α 1 ρ F ( t ) ( s ) = L ρ 0 t 1 Γ ( α ) ( t ρ s ρ ρ ) α 1 F ( s ) d s s 1 ρ ( s ) .
This implies that
f ( t ) f a ( t ) = 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 F ( s ) d s s 1 ρ .
Thus,
| f ( t ) f a ( t ) | = 1 Γ ( α ) | 0 t ( t ρ s ρ ρ ) α 1 F ( s ) d s s 1 ρ | 1 Γ ( α ) 0 t | ( t ρ s ρ ρ ) α 1 1 s 1 ρ | | F ( s ) | d s ε Γ ( α ) 0 t | ( t ρ s ρ ρ ) α 1 1 s 1 ρ | d s = ε Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 1 s 1 ρ d s ε Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 d t ρ s ρ ρ = t ρ α Γ ( α + 1 ) ρ α ε .
The proof is complete. □
Definition 14. 
From Definition 10, (8) shows that (1) is Hyers–Ulam stable with the constant K = T ρ α Γ ( α + 1 ) ρ α provided that 0 t T . (1) is not Hyers–Ulam stable if t = .
Definition 15. 
Let 0 < α < 1 , ρ > 0 , λ R , and g ( t ) be a given real continuous function on [ 0 , ) . If a function f : [ 0 , ) R satisfies the following inequality:
| ( D c α , ρ f ) ( t ) λ f ( t ) g ( t ) | ε , f ( 0 ) = f 0 R ,
for each t 0 and ε > 0 , then there exists a solution f a : [ 0 , ) R of (2) such that
| f ( t ) f a ( t ) | ( t ρ ρ ) α E α , α + 1 ( | λ | ( t α ρ ) α ) ε .
Proof. 
Let
F 1 ( t ) = ( D c α , ρ f ) ( t ) λ f ( t ) g ( t ) , t 0 .
Taking the ρ -Laplace transform of (14) via Definition 4, we have
L ρ { F 1 ( t ) } ( s ) = L ρ { ( D c α , ρ f ) ( t ) λ f ( t ) g ( t ) } ( s ) = L ρ { ( D c α , ρ f ) ( t ) } ( s ) λ L ρ { f ( t ) } ( s ) L ρ { g ( t ) } ( s ) = s α L ρ { f ( t ) } ( s ) s α 1 f 0 λ L ρ { f ( t ) } ( s ) L ρ { g ( t ) } ( s ) ,
where L ρ { F 1 ( · ) } denotes the ρ -Laplace transform of the function F 1 .
From (15), one has
L ρ { f ( t ) } ( s ) = s α 1 s α λ f 0 + 1 s α λ L ρ { g ( t ) } ( s ) + 1 s α λ L ρ { F 1 ( t ) } ( s ) = f 0 L ρ { E α ( λ ( t ρ ρ ) α ) } ( s ) + L ρ { ( t ρ ρ ) α 1 E α , α ( λ ( t ρ ρ ) α ) } ( s ) · L ρ { g ( t ) } ( s ) + L ρ { ( t ρ ρ ) α 1 E α , α ( λ ( t ρ ρ ) α ) } ( s ) · L ρ { F 1 ( t ) } ( s ) = f 0 L ρ { E α ( λ ( t ρ ρ ) α ) } ( s ) + L ρ { ( t ρ ρ ) α 1 E α , α ( λ ( t ρ ρ ) α ) ρ g ( t ) } ( s ) + L ρ { ( t ρ ρ ) α 1 E α , α ( λ ( t ρ ρ ) α ) ρ F 1 ( t ) } ( s ) = f 0 L ρ { E α ( λ ( t ρ ρ ) α ) } ( s ) + L ρ { 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ t ρ ρ ) α ) · g ( t ) d s s 1 ρ } + L ρ { 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) · F 1 ( t ) d s s 1 ρ } .
Set
f a ( t ) = f 0 E α ( λ ( t ρ ρ ) α ) + 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ t ρ ρ ) α ) · g ( t ) d s s 1 ρ .
Taking the ρ -Laplace transform of (17), one has
L ρ { f a ( t ) } ( s ) = f 0 L ρ { E α ( λ ( t ρ ρ ) α ) } ( s ) + L ρ { 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ t ρ ρ ) α ) · g ( t ) d s s 1 ρ } = f 0 L ρ { E α ( λ ( t ρ ρ ) α ) } ( s ) + L ρ { ( t ρ ρ ) α 1 E α , α ( λ ( t ρ ρ ) α ) ρ g ( t ) } ( s ) = f 0 L ρ { E α ( λ ( t ρ ρ ) α ) } ( s ) + L ρ { ( t ρ ρ ) α 1 E α , α ( λ ( t ρ ρ ) α ) } ( s ) · L ρ { g ( t ) } ( s ) = s α 1 s α λ f 0 + 1 s α λ L ρ { g ( t ) } ( s ) .
By Definition 4 and (18), we obtain
L ρ { ( D c α , ρ f a ) ( t ) λ f a ( t ) } ( s ) = s α L ρ { f a ( t ) } ( s ) s α 1 f 0 λ L ρ { f a ( t ) } ( s ) = ( s α λ ) s α 1 s α λ f 0 + ( s α λ ) 1 s α λ L ρ { g ( t ) } ( s ) s α 1 f 0 = L ρ { g ( t ) } ( s ) ,
which yields that f a is a solution of Equation (2) , since L ρ is one-to-one.
From (16) and (18), we have
L ρ { f ( t ) f a ( t ) } ( s ) = f 0 L ρ { E α ( λ ( t ρ ρ ) α ) } ( s ) + L ρ { 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) · g ( t ) d s s 1 ρ } + L ρ { 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) · F 1 ( s ) d s s 1 ρ } f 0 L ρ { E α ( λ ( t ρ ρ ) α ) } ( s ) L ρ { 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) · g ( t ) d s s 1 ρ } = L ρ { 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) · F 1 ( s ) d s s 1 ρ } .
This implies that
f ( t ) f a ( t ) = 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) · F 1 ( s ) d s s 1 ρ .
Thus,
| f ( t ) f a ( t ) | = | 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) · F 1 ( s ) d s s 1 ρ | 0 t | ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) | | F 1 ( s ) | d s s 1 ρ ε 0 t | ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) | d s s 1 ρ ε 0 t | ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) | d ( t ρ s ρ ρ ) = ε 0 t | ( t ρ s ρ ρ ) α 1 k = 0 ( λ ( t ρ s ρ ρ ) α ) k Γ ( α k + α ) | d ( t ρ s ρ ρ ) = ε k = 0 | λ | k Γ ( α k + α ) 0 t ( t ρ s ρ ρ ) α k + α 1 d ( t ρ s ρ ρ ) = ε k = 0 | λ | k Γ ( α k + α + 1 ) ( t ρ ρ ) α k + α = ( t ρ ρ ) α E α , α + 1 ( | λ | ( t α ρ ) α ) ε .
The proof is complete. □
Definition 16. 
From Definition 10, (13) shows that (2) is Hyers–Ulam stable with the constant K = T α ρ ρ α ( k + 1 ) E α , α + 1 ( | λ | T α ρ ) provided that 0 t T . (2) is not Hyers–Ulam stable if t = .
Definition 17. 
Let 0 < α < 1 , ρ > 0 , λ R , and g ( t ) be a given real continuous function on [ 0 , ) . If a function f : [ 0 , ) R satisfies the following inequality
| ( D c α , ρ f ) ( t ) λ f ( t ) g ( t ) | G ( t ) , f ( 0 ) = f 0 R ,
then
| F 1 ( t ) | G ( t )
for each t 0 and some function G ( t ) > 0 , where F 1 is defined in (14).
From Definition 15, there exists a solution f a : [ 0 , ) R of (2) such that
| f ( t ) f a ( t ) | = | 0 t ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) · F 1 ( s ) d s s 1 ρ 0 t | ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) | | F 1 ( s ) | d s s 1 ρ G ( t ) 0 t | ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) | d s s 1 ρ G ( t ) 0 t | ( t ρ s ρ ρ ) α 1 E α , α ( λ ( t ρ s ρ ρ ) α ) | d ( t ρ s ρ ρ ) G ( t ) 0 t ρ ρ z α 1 E α , α ( λ z α ) d z = ( t ρ ρ ) α E α , α + 1 ( | λ | ( t ρ ρ ) α ) G ( t ) .
By Definition 11, (13) shows that (2) is generalized Hyers–Ulam-Rassias stable with the constant c G = ( T ρ ρ ) α E α , α + 1 ( | λ | ( T ρ ρ ) α ) for all t [ 0 , T ] .
Definition 18. 
Let 0 < β < α < 1 , ρ > 0 , λ R , and g ( t ) be a given real continuous function on [ 0 , ) . If a function f : [ 0 , ) R satisfies the following inequality
| ( D c α , ρ f ) ( t ) λ ( D c β , ρ f ) ( t ) g ( t ) | ε , f ( 0 ) = f 0 R ,
for each t 0 and ε > 0 , then there exists a solution f a : [ 0 , ) R of (3) such that
| f ( t ) f a ( t ) | ( t ρ ρ ) α E α β , α + 1 ( | λ | ( t ρ ρ ) α β ) ε .
Proof. 
Let
F 2 ( t ) = ( D c α , ρ f ) ( t ) λ ( D c β , ρ f ) ( t ) g ( t ) , t 0 .
Taking the ρ -Laplace transform of (20) via Definition 4, we have
L ρ { F 2 ( t ) } ( s ) = L ρ { ( D c α , ρ f ) ( t ) λ ( D c β , ρ f ) ( t ) g ( t ) } ( s ) = L ρ { ( D c α , ρ f ) ( t ) } ( s ) λ L ρ { ( D c β , ρ f ) ( t ) } ( s ) λ L ρ { f ( t ) } ( s ) L ρ { g ( t ) } ( s ) = s α L ρ { f ( t ) } ( s ) s α 1 f 0 λ s β L ρ { f ( t ) } ( s ) + λ s β 1 f 0 L ρ { g ( t ) } ( s ) = ( s α λ s β ) L ρ { f ( t ) } ( s ) ( s α 1 λ s β 1 ) f 0 L ρ { g ( t ) } ( s ) ,
where L ρ { F 2 ( · ) } denotes the ρ -Laplace transform of the function F 2 .
Note
s α 1 s α λ s β = s α β 1 s α β λ = L ρ { E α β ( λ ( t ρ ρ ) α β ) } ( s ) ,
s β 1 s α λ s β = s 1 s α β λ = L ρ { ( t ρ ρ ) α β E α β , α β + 1 ( λ ( t ρ ρ ) α β ) } ( s ) ,
and
1 s α λ s β = s β s α β λ = L ρ { ( t ρ ρ ) α 1 E α β , α ( λ ( t ρ ρ ) α β ) ( s ) .
From (21), one has
L ρ { f ( t ) } ( s ) = s α 1 λ s β 1 s α λ s β f 0 + 1 s α λ s β L ρ { g ( t ) } ( s ) + 1 s α λ s β L ρ { F 2 ( t ) } ( s ) = s α 1 s α λ s β f 0 λ s β 1 s α λ s β f 0 + 1 s α λ s β L ρ { g ( t ) } ( s ) + 1 s α λ s β L ρ { F 2 ( t ) } ( s ) = f 0 L ρ { E α β ( λ ( t ρ ρ ) α β ) } ( s ) λ f 0 L ρ { ( t ρ ρ ) α β E α β , α β + 1 ( λ ( t ρ ρ ) α β ) } ( s ) + L ρ { ( t ρ ρ ) α 1 E α β , α ( λ ( t ρ ρ ) α β ρ g ( t ) ) ( s ) + L ρ { ( t ρ ρ ) α 1 E α β , α ( λ ( t ρ ρ ) α β ρ F 2 ( t ) ) ( s ) = f 0 L ρ { E α β ( λ ( t ρ ρ ) α β ) } ( s ) λ f 0 L ρ { ( t ρ ρ ) α β E α β , α β + 1 ( λ ( t ρ ρ ) α β ) } ( s ) + L ρ 0 t ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β · g ( t ) ) d x x 1 ρ ( s ) + L ρ 0 t ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β · F 2 ( t ) ) d x x 1 ρ ( s ) .
Set
f a ( t ) = f 0 E α β ( λ ( t ρ ρ ) α β ) λ f 0 ( t ρ ρ ) α β E α β , α β + 1 ( λ ( t ρ ρ ) α β ) ( s ) + 0 t ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β · g ( t ) ) d x x 1 ρ .
Taking the ρ -Laplace transform of (23), one has
L ρ { f a ( t ) } ( s ) = f 0 L ρ { E α β ( λ ( t ρ ρ ) α β ) } ( s ) λ f 0 L ρ ( t ρ ρ ) α β E α β , α β + 1 ( λ ( t ρ ρ ) α β ) ( s ) ( s ) + L ρ 0 t ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β · g ( t ) ) d x x 1 ρ ( s ) = f 0 L ρ { E α β ( λ ( t ρ ρ ) α β ) } ( s ) λ f 0 L ρ { ( ( t ρ ρ ) α β E α β , α β + 1 ( λ ( t ρ ρ ) α β ) ( s ) + L ρ { ( t ρ ρ ) α 1 E α β , α ( λ ( t ρ ρ ) α β ρ g ( t ) ) ( s ) = s α 1 s α λ s β f 0 λ s β 1 s α λ s β f 0 + 1 s α λ s β L ρ { g ( t ) } ( s ) = s α 1 λ s β 1 s α λ s β f 0 + 1 s α λ s β L ρ { g ( t ) } ( s ) .
By Definition 4 and (24), we obtain
L ρ { ( D c α , ρ f a ) ( t ) λ ( D c β , ρ f a ) ( t ) } = s α L ρ { f a ( t ) } ( s ) s α 1 f 0 λ s β L { f a ( t ) } ( s ) + λ s β 1 f 0 = ( s α λ s β ) s α 1 λ s β 1 s α λ s β f 0 + ( s α λ s β ) 1 s α λ s β L ρ { g ( t ) } ( s ) ( s α 1 λ s β 1 ) f 0 = L ρ { g ( t ) } ( s ) ,
which yields that f a is a solution of Equation (3), since L ρ is one-to-one.
From (22) and (24), we have
L ρ { f ( t ) f a ( t ) } ( s ) = s α 1 λ s β 1 s α λ s β f 0 + 1 s α λ s β L ρ { g ( t ) } ( s ) + 1 s α λ s β L ρ { F 2 ( t ) } ( s ) s α 1 λ s β 1 s α λ s β f 0 1 s α λ s β L ρ { g ( t ) } ( s ) . = 1 s α λ s β L ρ { F 2 ( t ) } ( s ) = L ρ 0 t ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β · F 2 ( t ) ) d x x 1 ρ ( s ) .
This implies that
f ( t ) f a ( t ) = 0 t ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β · F 2 ( t ) ) d x x 1 ρ ,
so
| f ( t ) f a ( t ) | = | 0 t ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β ) · F 2 ( t ) d x x 1 ρ | 0 t | ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β ) | | F 2 ( s ) | d s s 1 ρ ε 0 t | ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β ) | d s s 1 ρ ε 0 t | ( t ρ x ρ ρ ) α 1 E α β , α ( λ ( t ρ x ρ ρ ) α β ) | d ( t ρ s ρ ρ ) = ε 0 t | ( t ρ s ρ ρ ) α 1 k = 0 ( λ ( t ρ s ρ ρ ) α β ) k Γ ( ( α β ) k + α ) | d ( t ρ s ρ ρ ) ( t ρ ρ ) α E α β , α + 1 ( | λ | ( t ρ ρ ) α β ) ε .
The proof is complete. □
Definition 19. 
From Definition 10, (19) shows that (3) is Hyers–Ulam stable with the constant K = ( T ρ ρ ) α E α β , α + 1 ( | λ | ( T ρ ρ ) α β ) provided that 0 t T .
Definition 20. 
If β = 0 , then ( D c α , ρ f ) ( t ) λ ( D c β , ρ f ) ( t ) = g ( t ) coincides with ( D c α , ρ f ) ( t ) λ f ( t ) = g ( t ) , and ( t ρ ρ ) α E α β , α + 1 ( | λ | ( t ρ ρ ) α β ) coincides with ( t ρ ρ ) α E α , α + 1 ( | λ | ( t ρ ρ ) α ) ε , so Definition 18 generalizes Definition 15.
Definition 21. 
Let 0 < β < α < 1 , ρ > 0 , λ R , and g ( t ) be a given real continuous function on [ 0 , ) . If a function f : [ 0 , ) R satisfies the following inequality
| ( D c α , ρ f ) ( t ) λ ( D c β , ρ f ) ( t ) g ( t ) | G ( t ) , f ( 0 ) = f 0 R ,
then
| F 2 ( t ) | G ( t )
for each t 0 and some function G ( t ) > 0 , where F 2 is defined in (14) .
From Definition 18, there exists a solution f a : [ 0 , ) R of (2) such that
| f ( t ) f a ( t ) | t α ρ ρ ( α β ) k + α E α β , α + 1 ( | λ | t ρ ( α β ) ) G ( t )
By Definition 11, (25) shows (3) is generalized Hyers–Ulam–Rassias stable with the constant c G = T α ρ ρ ( α β ) k + α E α β , α + 1 ( | λ | T ρ ( α β ) ) for all t [ 0 , T ] .

4. Existence and Stability Results for the Nonlinear Equation

We introduce the following conditions:
[ A 1 ] : g : [ 0 , T ] × R R is continuous.
[ A 2 ] : There exists a L > 0 such that
| g ( t , y ) g ( t , f ) | L | y f | , y , f R , t [ 0 , T ] .
[ A 3 ] : There exists a constant L g > 0 such that
| g ( t , f ) | L g ( 1 + | f | )
for each t [ 0 , T ] and all t R .
Definition 22. 
Let 0 < α < 1 , ρ > 0 . Assume that [ A 1 ] and [ A 2 ] hold. If L ( T ρ ρ ) α < 1 , then (4) has a unique solution on [ 0 , T ] .
Proof. 
Consider Λ : C ( [ 0 , T ] , R ) C ( [ 0 , T ] , R ) as follows:
( Λ f ) ( t ) = f 0 + 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 g ( s , f ( s ) ) d s s 1 ρ , t [ 0 , T ] .
Note that Λ is well defined because of [ A 1 ] .
For all f 1 , f 2 C ( [ 0 , T ] , R ) and all t [ 0 , T ] , using [ A 2 ] , we have
| ( Λ f 1 ) ( t ) ( Λ f 2 ) ( t ) | 0 t ( t ρ s ρ ρ ) α 1 | g ( s , f 1 ( s ) ) g ( s , f 2 ( s ) ) | d s s 1 ρ L 0 t ( t ρ s ρ ρ ) α 1 | f 1 ( s ) f 2 ( s ) | d s s 1 ρ = L f 1 f 2 C ( t ρ ρ ) α L ( T ρ ρ ) α f 1 f 2 C ,
which implies
Λ f 1 Λ f 2 C L ( T ρ ρ ) α f 1 f 2 C .
From the condition L ( T ρ ρ ) α < 1 , Λ is a contraction mapping, and, by applying the Banach contraction mapping principle, we know that the operator Λ has a unique fixed point on [ 0 , T ] . □
Next, we show that the existence of solutions for (4) via Schaefer’s fixed point theorem.
Definition 23. 
Assume that [ A 1 ] and [ A 3 ] hold. Then, (4) has at least one solution.
Proof. 
Consider Λ as in (26). We divide our proof into several steps.
Step 1. Λ is continuous.
Let { f n } be a sequence such that f n f in C ( [ 0 , T ] , R ) . For all t [ 0 , T ] , we get
| Λ f n ( t ) Λ f ( t ) | = 1 Γ ( α ) | 0 t ( t ρ s ρ ρ ) α 1 g ( s , f n ( s ) ) d s s 1 ρ 0 t ( t ρ s ρ ρ ) α 1 g ( s , f ( s ) ) d s s 1 ρ | 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 | g ( s , f n ( s ) ) g ( s , f ( s ) ) | d s s 1 ρ g ( · , f n ) g ( · , f ) C 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 d s s 1 ρ 1 Γ ( α ) ( T ρ ρ ) α g ( · , f n ) g ( · , f ) C .
This shows that Λ is continuous since g ( · , f n ) g ( · , f ) C 0 when n .
Step 2. Λ maps bounded sets into bounded sets of C ( [ 0 , T ] , R ) .
Indeed, we prove that for all r > 0 , there exists a k > 0 such that for every f B r = { f C ( [ 0 , T ] , R ) : f C r } , we have Λ f C k . In fact, for any t [ 0 , T ] , from [ A 3 ] , we have
| Λ f ( t ) | | f 0 | + 1 Γ ( α ) | 0 t ( t ρ s ρ ρ ) α 1 g ( s , f ( s ) ) d s s 1 ρ | | f 0 | + 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 | g ( s , f ( s ) ) | d s s 1 ρ | f 0 | + L g ( 1 + f C ) Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 d s s 1 ρ | f 0 | + L g ( 1 + r ) Γ ( α ) ( T ρ ρ ) α ,
which implies that
Λ f C | f 0 | + + L g ( 1 + r ) Γ ( α ) ( T ρ ρ ) α : = k .
Step 3. P maps bounded sets into equicontinuous sets in C ( [ 0 , T ] , R ) .
Let t 1 , t 2 [ 0 , T ] , with 0 t 1 < t 2 T , f B r . From [ A 3 ] , we have
| Λ f ( t 1 ) Λ f ( t 2 ) | = 1 Γ ( α ) | 0 t 1 ( t 1 ρ s ρ ρ ) α 1 g ( s , f ( s ) ) d s s 1 ρ 0 t 2 ( t 2 ρ s ρ ρ ) α 1 g ( s , f ( s ) ) d s s 1 ρ | L g ( 1 + r ) Γ ( α ) 0 t 1 ( t 1 ρ s ρ ρ ) α 1 ( t 2 ρ s ρ ρ ) α 1 d s s 1 ρ + t 1 t 2 ( t 2 ρ s ρ ρ ) α 1 d s s 1 ρ L g ( 1 + r ) Γ ( α ) α ρ α t 1 α ρ t 2 α ρ 2 ( t 2 ρ t 1 ρ ) α .
Then, the right-hand side of the above inequality tends to zero as t 1 t 2 . Thus, Λ is equicontinuous.
We can conclude that Λ is completely continuous from Steps 1–3 with the Arzela–Ascoli theorem.
Step 4. A priori bounds.
Now, we show that the set E ( Λ ) = { f C ( [ 0 , T ] , R ) : f = ι f for some ι ( 0 , 1 ) } is bounded.
Let f E ( Λ ) . Then, f = ι Λ f for some ι ( 0 , 1 ) } . For each t [ 0 , T ] , we have
| f ( t ) | f 0 + 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 | g ( s , f ( s ) ) | d s s 1 ρ f 0 + L g Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 | 1 + f ( s ) | d s s 1 ρ f 0 + L g Γ ( α + 1 ) ρ α t α ρ + L g Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 | f ( s ) | d s s 1 ρ .
By Definition 12, we obtain
| f ( t ) | ( f 0 + L g Γ ( α + 1 ) ρ α t α ρ ) E α ( L g ( t ρ ρ ) α ) ( f 0 + L g Γ ( α + 1 ) ρ α T α ρ ) E α ( L g ( T ρ ρ ) α ) < .
Then, the set E ( Λ ) is bounded.
Schaefer’s fixed point theorem guarantees that Λ has a fixed point, which is a solution of (4). The proof is finished. □
For the sake of discussion, the following inequality is given
| ( D c α , ρ f ) ( t ) g ( t , f ( t ) ) | G ( t ) , t [ 0 , T ] .
In the following, we consider (4) and (27) to discuss the generalized Ulam–Hyers–Rassias stability.
We need the following condition.
[ A 4 ] : Let G C ( [ 0 , T ] , R + ) be an increasing function and there exists λ G > 0 such that
0 t ( s ρ ) ( t ρ s ρ ) α 1 G ( s ) d s λ G G ( t ) , t [ 0 , T ] .
Definition 24. 
Assumptions [ A 1 ] , [ A 2 ] , and [ A 4 ] hold. If L ( T ρ ρ ) α < 1 , then (4) is generalized Ulam–Hyers–Rassias stable with respect to G on [ 0 , T ] .
Proof. 
Let f C ( [ 0 , T ] , R ) be a solution of (4). From Definition 22,
( D c α , ρ f ) ( t ) = g ( t , f ( t ) ) , 0 < α < 1 , t [ 0 , T ] , f ( 0 ) = f 0
has the unique solution
h ( t ) = f 0 + 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 g ( s , h ( s ) ) d s s 1 ρ .
Integrating the inequality (27) from 0 to t and using the condition [ A 4 ] , we have
| f ( t ) f 0 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 g ( s , f ( s ) ) d s s 1 ρ | 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 G ( s ) d s s 1 ρ = 1 Γ ( α ) ρ α 0 t ( s ρ ) ( t ρ s ρ ) α 1 G ( s ) d s 1 Γ ( α ) ρ α λ G G ( t ) .
Thus,
| f ( t ) h ( t ) | | f ( t ) f 0 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 g ( s , h ( s ) ) d s s 1 ρ | | f ( t ) f 0 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 g ( s , f ( s ) ) d s s 1 ρ + 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 g ( s , f ( s ) ) d s s 1 ρ 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 g ( s , h ( s ) ) d s s 1 ρ | | f ( t ) f 0 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 g ( s , f ( s ) ) d s s 1 ρ | + 1 Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 | g ( s , f ( s ) ) g ( s , h ( s ) ) | d s s 1 ρ 1 Γ ( α ) ρ α λ G G ( t ) + L Γ ( α ) 0 t ( t ρ s ρ ρ ) α 1 | f ( s ) h ( s ) | d s s 1 ρ .
From Definition 12, we obtain
| f ( t ) h ( t ) | ( λ G Γ ( α ) ρ α ) E α ( L Γ ( α ) ρ α Γ ( α ) t ρ α ) G ( t ) ( λ G Γ ( α ) ρ α ) E α ( L ( T ρ ρ ) α ) G ( t ) , t [ 0 , T ] .
Set K = ( λ G Γ ( α ) ρ α ) E α ( L ( T ρ ρ ) α ) . One has
| f ( t ) h ( t ) | K G ( t ) , t [ 0 , T ] .
From Definition 11, (4) is generalized Ulam–Hyers–Rassias stable with respect to G on [ 0 , T ] . The proof is complete. □

5. An Example

In this section, an example is given to illustrate our main results.
Definition 25. 
We consider the following fractional problem
( D 1 2 , 2 f ) ( t ) = e 2 t 5 | f ( t ) | 1 + | f ( t ) | , t [ 0 , 3 ] ,
and the inequality
| ( D 1 2 , 2 f ) ( t ) e 2 t 5 | f ( t ) | 1 + | f ( t ) | | G ( t ) , t [ 0 , 3 ] .
Set α = 1 2 , ρ = 2 , T = 3 and g ( t , f ) = e 2 t 5 | f | 1 + | f | , ( t , f ) [ 0 , 3 ] × R . For all t [ 0 , 3 ] and f , q R ,
| g ( t , f ( t ) ) g ( t , q ( t ) ) | e 2 t 5 | | f ( t ) | | q ( t ) | | 1 5 | f ( t ) q ( t ) | .
Set L = 1 5 . Then, L ( T ρ ρ ) α = 1 5 3 2 2 = 3 2 10 < 1 .
Let G ( t ) = e t , t [ 0 , 3 ] and λ G = 3 2 > 0 . Note
0 t ( s ρ ) ( t ρ s ρ ) α 1 G ( s ) d s = 0 t 2 s ( t 2 s 2 ) 1 2 e s d s 1 2 t e t λ G G ( t ) , t [ 0 , 3 ] .
Thus, all the assumptions in Definition 22 and Definition 24 being satisfied, our results can be applied to the problem (29).

Author Contributions

M.F. and J.W. contributed to the supervision and project administration, K.L., M.F., and J.W. contributed to the conceptualization and methodology. All authors have read and approved the final manuscript.

Funding

This work is partially supported by Guizhou Provincial Science and Technology Foundation ([2020]1Y002), the Slovak Research and Development Agency under the Contract No. APVV-18-0308, and the Slovak Grant Agency VEGA No. 1/0358/20 and No. 2/0127/20.

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.; Kilbas, A.; Marichev, O. Fractional Integrals and Derivatives, Theory and Applications; Gordon and Breach Science Publishers: Amsterdam, The Netherlands, 1993. [Google Scholar]
  3. Kilbas, A.; Srivastava, H.; Trujillo, J. Theory and Applications of Fractional Differential Equations; Elsevier B.V.: Amsterdam, The Netherlands, 2006. [Google Scholar]
  4. Martínez-García, M.; Zhang, Y.; Gordon, T. Memory pattern identification for feedback tracking control in human-machine systems. Hum. Factors 2019. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  5. Martínez-García, M.; Gordon, T. Human control of systems with fractional order dynamics. In Proceedings of the IEEE International Conference on Systems, Man, and Cybernetics, Budapest, Hungary, 9–12 October 2016. [Google Scholar]
  6. Caputo, M.; Fabrizio, M. A new definition of fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 73–85. [Google Scholar]
  7. Liu, K.; Wang, J.; O’Regan, D. Ulam-Hyers-Mittag–Leffler stability for ψ-Hilfer fractional-order delay differential equations. Adv. Differ. Equ. 2019, 50, 1–12. [Google Scholar] [CrossRef]
  8. Losada, J.; Nieto, J. Properties of a new fractional derivative without singular kernel. Prog. Fract. Differ. Appl. 2015, 1, 87–92. [Google Scholar]
  9. Baleanu, D.; Mousalou, A.; Rezapour, S. On the existence of solutions for some infinite coefficient-symmetric Caputo–Fabrizio fractional integro-differential equations. Bound. Value Prob. 2017, 145, 1–9. [Google Scholar] [CrossRef]
  10. Goufo, F.E.D. Application of the Caputo–Fabrizio fractional derivative without singular kernel to Korteweg-de Vries-Burgers equations. Math. Model. Anal. 2016, 21, 188–198. [Google Scholar] [CrossRef]
  11. Moore, E.; Sirisubtawee, S.; Koonparasert, S. A Caputo–Fabrizio fractional differential equation model for HIV/AIDS with treatment compartment. Adv. Differ. Equ. 2019, 2019, 200. [Google Scholar] [CrossRef] [Green Version]
  12. Dokuyucu, M.; Celik, E.; Bulut, H.; Baskonus, H. Cancer treatment model with the Caputo–Fabrizio fractional derivative. Eur. Phys. J. Plus 2018, 133, 92. [Google Scholar] [CrossRef]
  13. Başcı, Y.; Öğrekçi, S.; Mısır, A. On Hyers–Ulam stability for fractional differential equations including the new Caputo–Fabrizio fractional derivative. Mediterr. J. Math. 2019, 16, 131. [Google Scholar] [CrossRef]
  14. Butzer, P.; Kilbas, A.; Trujillo, J. Compositions of Hadamard-type fractional integration operators and the semigroup property. J. Math. Anal. Appl. 2002, 269, 387–400. [Google Scholar] [CrossRef] [Green Version]
  15. Butzer, P.; Kilbas, A.; Trujillo, J. Fractional calculus in the mellin setting and Hadamard-type fractional integrals. J. Math. Anal. Appl. 2002, 269, 1–27. [Google Scholar] [CrossRef] [Green Version]
  16. Butzer, P.; Kilbas, A.; Trujillo, J. Mellin transform analysis and integration by parts for hadamard-type fractional integrals. J. Math. Anal. Appl. 2002, 270, 1–15. [Google Scholar] [CrossRef] [Green Version]
  17. Kilbas, A. Hadamard-type fractional calculus. J. Korean Math. Soc. 2001, 38, 1191–1204. [Google Scholar]
  18. Kilbas, A.; Trujillo, J. Hadamard-type integrals as g-transforms. Integr. Transf. Spec. F. 2003, 14, 413–427. [Google Scholar] [CrossRef]
  19. Katugampola, U. New approach to a genaralized fractional integral. Appl. Math. Comput. 2011, 218, 860–865. [Google Scholar]
  20. Katugampola, U. A new approach to generalized fractional derivatives. Bull. Math. Anal. Appl. 2014, 6, 1–15. [Google Scholar]
  21. Başcı, Y.; Mısır, A.; Öğrekçi, S. On the stability problem of differential equations in the sense of Ulam. Results Math. 2020, 75, 6. [Google Scholar] [CrossRef]
  22. Alqifiary, Q.; Jung, S. Laplace transform and generalized Hyers–Ulam stability of linear differential equations. Electron. J. Differ. Equ. 2014, 2014, 1–11. [Google Scholar]
  23. Huang, J.; Li, Y. Hyers–Ulam stability of delay differential equations of first order. Math. Nachr. 2016, 289, 60–66. [Google Scholar] [CrossRef]
  24. da C. Sousa, J.; Capelas de Oliveira, E. Ulam-Hyers stability of a nonlinear fractional Volterra integro-differential equation. Appl. Math. Lett. 2018, 81, 50–56. [Google Scholar]
  25. Wang, J.; Zhou, Y.; Fečkan, M. Nonlinear impulsive problems for fractional differential equations and Ulam stability. Comput. Math. Appl. 2012, 64, 3389–3405. [Google Scholar] [CrossRef] [Green Version]
  26. Capelas de Oliveira, E.; da C. Sousa, J. Ulam-Hyers-Rassias stability for a class of fractional integro-differential equations. Results Math. 2018, 73, 111. [Google Scholar] [CrossRef] [Green Version]
  27. da C. Sousa, J.; Kucche, K.; Capelas de Oliveira, E. Stability of ψ-Hilfer impulsive fractional differential equations. Appl. Math. Lett. 2018, 88, 73–80. [Google Scholar]
  28. da C. Sousa, J.; Capelas de Oliveira, E. On the Ulam-Hyers-Rassias stability for nonlinear fractional differential equations using the ψ-Hilfer operator. J. Fixed Point Theory Appl. 2018, 20, 5–21. [Google Scholar]
  29. Shah, K.; Ali, A.; Bushnaq, S. Hyers–Ulam stability analysis to implicit Cauchy problem of fractional differential equations with impulsive conditions. Math. Meth. Appl. Sci. 2018, 41, 8329–8343. [Google Scholar] [CrossRef]
  30. Satco B., R. Ulam-type stability for differential equations driven by measures. Math. Nachr. 2020, 293, 147–157. [Google Scholar] [CrossRef]
  31. Rezaei, H.; Jung, S.; Rassias, T. Laplace transform and Hyers–Ulam stability of linear differential equations. J. Math. Anal. Appl. 2013, 403, 244–251. [Google Scholar] [CrossRef]
  32. Wang, C.; Xu, T. Hyers–Ulam stability of fractional linear differential equations involving Caputo fractional derivatives. Appl. Math. 2015, 60, 383–393. [Google Scholar] [CrossRef] [Green Version]
  33. Shen, Y.; Chen, W. Laplace transform mathod for the ulam stability of linear fractional differential equations with constant coefficients. Mediterr. J. Math. 2017, 14, 25. [Google Scholar] [CrossRef]
  34. Liu, K.; Fečkan, M.; O’Regan, D.; Wang, J. Hyers–Ulam stability and existence of solutions for differential equations with Caputo–Fabrizio fractional derivative. Mathematics 2019, 7, 333. [Google Scholar] [CrossRef] [Green Version]
  35. Liu, K.; Wang, J.; Zhou, Y.; O’Regan, D. Hyers–Ulam stability and existence of solutions for fractional differential equations with Mittag–Leffler kernel, Chaos. Solitons Fractals 2020, 132, 109534. [Google Scholar] [CrossRef]
  36. Fahd, J.; Thabet, A. A modified Laplace transform for certain generalized fractional operators. Results Nonlinear Anal. 2018, 1, 88–98. [Google Scholar]
  37. Ndolane, S.; Gautam, S. Generalized Mittag–Leffler input stability of the fractional differential equations. Symmetry 2019, 11, 608. [Google Scholar]
  38. Wang, J.; Lv, L.; Zhou, Y. Ulam stability and data depenaence for fractional differential equations with Caputo derivative. Elect. J. Qual. Theory Diff. Equ. 2011, 63, 1–10. [Google Scholar]
  39. Adjabi, Y.; Jarad, F.; Abdeljawad, T. On generalized fractional operators and a Gronwall type inequality with applications. Filomat 2017, 31, 5457–5473. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Liu, K.; Fečkan, M.; Wang, J. Hyers–Ulam Stability and Existence of Solutions to the Generalized Liouville–Caputo Fractional Differential Equations. Symmetry 2020, 12, 955. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12060955

AMA Style

Liu K, Fečkan M, Wang J. Hyers–Ulam Stability and Existence of Solutions to the Generalized Liouville–Caputo Fractional Differential Equations. Symmetry. 2020; 12(6):955. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12060955

Chicago/Turabian Style

Liu, Kui, Michal Fečkan, and Jinrong Wang. 2020. "Hyers–Ulam Stability and Existence of Solutions to the Generalized Liouville–Caputo Fractional Differential Equations" Symmetry 12, no. 6: 955. https://0-doi-org.brum.beds.ac.uk/10.3390/sym12060955

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