Nothing Special   »   [go: up one dir, main page]

Next Article in Journal
Vibration Systems with Fractional-Order and Distributed-Order Derivatives Characterizing Viscoinertia
Next Article in Special Issue
Fractional Integral Inequalities for Exponentially Nonconvex Functions and Their Applications
Previous Article in Journal
Adsorption on Fractal Surfaces: A Non Linear Modeling Approach of a Fractional Behavior
Previous Article in Special Issue
Fractals Parrondo’s Paradox in Alternated Superior Complex System
You seem to have javascript disabled. Please note that many of the page functionalities won't work as expected without javascript enabled.
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence, Uniqueness, and Eq–Ulam-Type Stability of Fuzzy Fractional Differential Equation

1
Faculty of Mathematics and Computational Science, Xiangtan University, Xiangtan 411105, China
2
Department of Mathematics and Statistics, University of Lahore, Sargodha 40100, Pakistan
3
College of Mathematics and Information Science, Guangxi University, Nanning 530004, China
*
Author to whom correspondence should be addressed.
Fractal Fract. 2021, 5(3), 66; https://doi.org/10.3390/fractalfract5030066
Submission received: 26 April 2021 / Revised: 4 July 2021 / Accepted: 5 July 2021 / Published: 11 July 2021
(This article belongs to the Special Issue Fractional Order Systems and Their Applications)

Abstract

:
This paper concerns with the existence and uniqueness of the Cauchy problem for a system of fuzzy fractional differential equation with Caputo derivative of order q ( 1 , 2 ] , 0 c D 0 + q u ( t ) = λ u ( t ) f ( t , u ( t ) ) B ( t ) C ( t ) , t [ 0 , T ] with initial conditions u ( 0 ) = u 0 , u ( 0 ) = u 1 . Moreover, by using direct analytic methods, the E q –Ulam-type results are also presented. In addition, several examples are given which show the applicability of fuzzy fractional differential equations.

1. Introduction

In real-life phenomena, numerous physical processes are used to present fractional-order sets that may change with space and time. The operations of differentiation and integration of fractional order are authorized by fractional calculus. The fractional order may be taken on imaginary and real values [1,2,3]. The theory of fuzzy sets is continuously drawing the attention of researchers. This is mainly due to its extended adaptability in various fields including mechanics, engineering, electrical, processing signals, thermal system, robotics, control, signal processing, and in several other areas [4,5,6,7,8,9,10]. Therefore, it has been a topic of increasing concern for researchers during the past few years.
Fuzzy fractional differential equations appeared for the first time in 2010 when an idea of the solution was initially proposed by Agarwal et al. [11]. However, the Riemann–Liouville H derivative based on the strongly generalizing Hukuhara differentiability [12,13] was defined by Allahviranloo and Salahshour [14,15]. They worked on solutions to Cauchy problems under this kind of derivative.
0 R L D a + q u ( t ) = λ u ( t ) + f ( t ) , t [ a , b ] , 0 R L D a + q 1 u ( t ) = u 0 E 1 .
In the above, q ( 0 , 1 ] , through using Laplace transforms [13] and Mittag–Leffler functions [12]. By using fractional hyperbolic functions and the properties of these functions, Chehlabi et al. obtain some new results [16]. More latest studies on fuzzy fractional differential equations can be found through references [17,18,19,20,21,22].
In 1940, Ulam promoted the Ulam stability. Lately, Hyers and Rassias used this concept of stability. Since then, in mathematical analysis and differential equations, the Ulam-type stability has had great significance. In fractional differential equations, E α –Ulam-type stabilities were promoted by Wang in 2014 [23].
0 c D t α u ( t ) + λ u ( t ) = f ( t , u ( t ) ) , t [ 0 , T ] u ( 0 ) = u 0 R .
In the above equation, α ( 0 , 1 ] and λ > 0 . Shen studied the Ulam stability under the generalization of Hukuhara differentiability of a first-order linear fuzzy differential equation in 2015 [24]. Later, Shen et al. investigated the Ulam stability of a nonlinear fuzzy fractional equation with the help of fixed-point techniques in 2016 [25],
0 R L D 0 + q u ( t ) = λ u ( t ) f ( t , u ( t ) ) , t [ 0 , T ] ,
by focusing on the initial condition
0 R L D 0 + q 1 u ( 0 ) = u 0 E 1 ,
where 0 R L D 0 + q denoted Riemann–Liouville H derivative with respect to order q ( 0 , 1 ] , f : ( 0 , T ] × E 1 E 1 , T R + , and λ R .
More results can be observed that are related to Ulam-type stability in [26,27,28]. Motivated by the above-cited papers, we aim to deal with fuzzy fractional differential equations of the form,
0 c D 0 + q u ( t ) = λ u ( t ) f ( t , u ( t ) ) B ( t ) C ( t ) , t [ 0 , T ] ,
with initial conditions
u ( 0 ) = u 0 , u ( 0 ) = u 1 ,
Here, 0 c D 0 + q denotes the Caputo derivative of order q ( 1 , 2 ] , f : ( 0 , T ] × E 1 E 1 , T R + and λ R .
This paper focuses on facilitating, with as few conditions as possible, to assure the uniqueness and existence of a solution to Cauchy problems (1) and (2). It establishes a link between fuzzy fractional differential equations and the Ulam-type stability, which enhances and generalizes some familiar outputs in the existing literature.

2. Basic Concepts

Assume that P k ( R ) denotes the collection of all nonempty convex and compact subsets of R and define sums and scalar products in P k ( R ) in the usual manner. Let A and B be two nonempty bounded subsets in R . The distance between A and B is defined through the Hausdorff metric,
D ( A , B ) = max { sup a A inf b B a b , sup b B inf a A a b } .
In the above equality, x stands for the usual Euclidean norm in R . Now it is well known that the metric D turns the space ( P k ( R ) , D ) into a complete and separable metric space [26].
Denote
E 1 = { u : R [ 0 , 1 ] |   u   s a t i s f i e s ( 1 ) ( 4 ) }
where (1)–(4) stands for the following properties of the function u:
(1)
u is normal in the sense that there exists an s 0 R such that u ( s 0 ) = 2 ;
(2)
u is fuzzy convex, that is u ( q s + ( 1 q ) y ) m i n { u ( s ) , u ( y ) } for any s , y R and q ( 1 , 2 ] ;
(3)
u is an upper semicontinuous function on R ;
(4)
The set [ u ] 1 defined by [ u ] 1 = { t R | u ( t ) > 1 } ¯ is compact.
For 1 < q 2 , denote [ u ] q = { t R | u ( t ) q } ¯ . Now, from (1)–(4), it follows that the q-level set [ u ] q P k ( R ) 1 q 2 .
Define u ̲ as the lower branch and u ¯ as the upper branch of the fuzzy number u E 1 . The set [ u ] q = { t R | u ( t ) q } : = [ u ̲ q , u ¯ q ] is known as the q-level set of fuzzy number u, where q ( 1 , 2 ] . The length of q-level set is calculated as diam [ u ] q = u ¯ q u ̲ q .
Lemma 1
([29,30]). If u , v , s , y E 1 , then
(i)
( E 1 , D ) is a complete metric space;
(ii)
D ( u s , v s ) = D ( u , v ) ;
(iii)
D ( λ u , λ v ) = | λ | D ( u , v ) λ R ;
(iv)
D ( u s , v y ) D ( u , v ) + D ( s , y ) ;
(v)
D ( λ u , μ u ) = | λ μ | D ( u , 0 ^ ) , λ , μ 0 .
Let C E [ a , b ] and L E [ a , b ] be spaces for all continuous and Lebesgue integrable fuzzy-valued functions on [ a , b ] , respectively. Moreover, ( C E [ a , b ] , D ) stands for the complete metric space, where
D ( u , v ) = sup t [ a , b ] d ( u ( t ) , v ( t ) ) .
Remark 1.
On E 1 , we can define the subtraction , called the H difference as follows: u v makes sense if there exists ω E 1 such that u = v ω . Then, by definition, ω = u v .
Let u , v E 1 be such that u v is well defined. Then, its q- l e v e l is determined by
[ u v ] q = [ u ̲ q v ̲ q , u ¯ q v ¯ q ] .
Through a generalization of the Hausdorff–Pompeiu metric on convex and compact sets, the metric D on E 1 can be defined by
D ( u , v ) = sup 1 q 2 max { | u ̲ q v ̲ q | , | u ¯ q u ¯ q | } .
Definition 1
([13]). Assume that F C E ( a , b ] L E ( a , b ] . The fuzzy Riemann–Liouville integral for a fuzzy-valued function F is defined by
T q F ( t ) = 1 Γ ( q ) a t F ( x ) ( t x ) 1 q d x , t ( a , b ) ,
q ( 1 , 2 ] . For q = 1 we obtain T 1 F ( t ) = a t F ( x ) d s , which is the classical fuzzy integral operator.
Definition 2
([13]). Assume that F C E ( a , b ] L E ( a , b ] , t 1 ( a , b ) and
ϕ ( t ) = 1 Γ ( 1 q ) a t ( t x ) q F ( x ) d x .
It is said that F is Caputo H-differentiable of order 1 < q 2 at t 1 , if there exists an element 0 c D t q F ( t 1 ) E 1 such that the following fuzzy equalities are valid:
(i)
0 c D t q F ( t 1 ) = lim h 0 + ϕ ( t 1 + h ) ϕ ( t 1 ) h
(ii)
0 c D t q F ( t 1 ) = lim h 0 + ϕ ( t 1 ) ϕ ( t 1 h ) h
(iii)
0 c D t q F ( t 1 ) = lim h 0 + ϕ ( t 1 ) ϕ ( t 1 + h ) h
(iv)
0 c D t q F ( t 1 ) = lim h 0 + ϕ ( t 1 h ) ϕ ( t 1 ) h
Here, we use only the first two cases [23]. These derivatives are trivial because they reduce to crisp elements. Regarding other fuzzy cases, the reader is referred to [23]. Furthermore, regarding this simplicity, a fuzzy-valued function F is called c [(i)-GH]-differentiable or c [(ii)-GH]-differentiable if it is differentiable according to concept (i) or to (ii) of Definition 2, respectively.
The Mittag–Leffler and fractional hyperbolic functions frequently occur in solutions to fractional systems; see, e.g., [16,23]. The Mittag–Leffler functions in the form of a single and a double parameter are defined by, respectively,
E α ( x ) = k = 1 x k Γ ( α k + 1 ) E α , β ( x ) = k = 1 x k α k + β .
Some properties of these functions can be found in [31,32,33].
Lemma 2.
Let δ > 0 . Some properties of the functions E α ( . ) and E α , β ( . ) are listed below:
(i)
Let 1 < α < 2 . Then E α ( δ t α ) 2 and E α , α ( δ t α ) 1 Γ ( α ) ;
(ii)
Let 1 < α 2 and β < α + 1 . Then E α ( . ) and E α , β ( . ) are positive. If, moreover, 0 t 2 t 3 , then E α ( δ t 2 α ) E α ( δ t 3 α ) and E α , β ( δ t 2 α ) E α , β ( δ t 3 α ) ;
(iii)
0 z E α , β ( t α ) t β 1 d t = z β E α , β + 1 ( z α ) , α > 1 .
Remark 2.
According to the lemma given above, it can be observed that E α , α ( s ) 1 Γ ( α ) E α , α ( s ) for α ( 1 , 2 ] and s R + . Fractional hyperbolic functions that are generalizations of standard hyperbolic functions can be defined through Mittag–Leffler functions (see, e.g., [16]) as follows:
c o s h α , β ( s ) = k = 0 s 2 k Γ ( 2 α k + β ) = E 2 α , β ( s 2 ) ,
s i n h α , β ( s ) = k = 0 s 2 k + 1 Γ ( 2 α k + α + β ) = s E 2 α , α + β ( s 2 ) ,
for α , β > 1 . It is noticed that cosh α , β ( s ) is an even function and that sinh α , β ( s ) , s R , is an odd function. For α = β , we write C h α ( s ) and S h α ( s ) instead of cosh α , α ( s ) and sinh α , α ( s ) respectively. It is not difficult to observe that C h α ( s ) + S h α ( s ) = E α , α ( s ) and C h α ( s ) S h α ( s ) = E α , α ( s ) , s R (see, e.g., [16]).
Remark 3.
According to the above arguments and Remark 2, we have | C h α ( s ) ± S h β ( s ) | E α , α ( | s | ) for any s R .
Lemma 3.
(Gronwall lemma) [34] Let μ, v C ( [ 0 , 1 ] , R + ) . Suppose μ is increasing. If s C ( [ 0 , 1 ] , R + ) obeys the inequality
s ( t ) μ ( t ) + 0 t v ( x ) s ( x ) d x , t [ 0 , 1 ] ,
then
s ( t ) μ ( t ) e x p 0 t v ( x ) s ( x ) d x , t [ 0 , 1 ] .

3. Existence and Uniqueness Results

In this part, existence and uniqueness of solutions to the Cauchy problem in (1) and (2) are discussed. We can start with the lemma given below.
Lemma 4
([16]). When λ > 1 , the c [(i)-GH]-differentiable solution to problem (1) is given by
u ( t ) = E q , 1 ( λ t q ) u 0 t E q , 2 ( λ t q ) u 1 0 t E q , q ( λ ( t x ) q ) f ( x ) ( t x ) 1 q d x ;
when λ < 1 , the c [(ii)-GH]-differentiable solution to problem (1) is given by
u ( t ) = E q , 1 ( λ t q ) u 0 ( 1 ) t E q , 2 ( λ t q ) u 1 ( 1 ) 0 t E q , q ( λ ( t x ) q ) f ( x ) ( t x ) 1 q d x ;
when λ < 1 , the c [(i)-GH]-differentiable solution to problem (1) is given by
u ( t ) = [ C h q , 1 ( λ t q ) u 0 S h q , 1 ( λ t q ) u 0 ] [ t C h q , 2 ( λ t q ) u 1 t S h q , 2 ( λ t q ) u 1 ] 0 t C h q λ f ( x ) S h q λ f ( x ) ( t x ) 1 q d x ;
when λ > 1 , the c [(ii)-GH]-differentiable solution to problem (1) is given by
u ( t ) = [ C h q , 1 ( λ t q ) u 0 ( 1 ) S h q , 1 ( λ t q ) u 0 ] ( 1 ) [ t C h q , 2 ( λ t q ) u 1 ( 1 ) t S h q , 2 ( λ t q ) u 1 ] ( 1 ) 0 t C h q λ f ( x ) ( 1 ) S h q λ f ( x ) ( t x ) 1 q d x ;
when λ = 1 , the c [(i)-GH]-differentiable solution to problem (1) is given by
u ( t ) = u 0 t u 1 0 t f ( x ) ( t x ) 1 q d x ;
when λ = 1 , the c [(ii)-GH]-differentiable solution to problem (1) is given by
u ( t ) = u 0 ( 1 ) t u 1 ( 1 ) 0 t f ( x ) ( t x ) 1 q d x ;
Remark 4.
If λ = 0 , then problem (1) reduces to
0 c D t q u ( t ) = f ( t ) , 0 c D 0 + q 1 u ( 0 ) = u 0 E 1 , 0 c D 0 + q 1 u ( 0 ) = u 1
By applying Lemma 4 and Remark 4 with f ( t , u ( t ) ) B ( t ) C ( t ) instead of f ( t ) , it follows that the Cauchy problem in (1) and (2) possesses an integral version. In case λ 1 and the function t u ( t ) , t [ 0 , T ] is assumed to be c [(i)-GH]-differentiable, then the function u satisfies
u ( t ) = E q , 1 ( λ t q ) u 0 t E q , 2 ( λ t q ) u 1 ( 1 ) 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x .
In case λ 1 and the function t u ( t ) is supposed to be c [(ii)-GH]-differentiable, then the function u satisfies
u ( t ) = E q , 1 ( λ t q ) u 0 ( 1 ) t E q , 2 ( λ t q ) u 1 ( 1 ) 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x .
In case λ < 1 and the function t u ( t ) is c [(i)-GH]-differentiable, then the function u satisfies
u ( t ) = [ C h q , 1 ( λ t q ) u 0 S h q , 1 ( λ t q ) u 0 ] [ t C h q , 2 ( λ t q ) u 1 t S h q , 2 ( λ t q ) u 1 ] 0 t C h q λ [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] S h q λ [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x .
In case λ > 1 and the function t u ( t ) is c [(ii)-GH]-differentiable, then the function u satisfies
u ( t ) = [ C h q , 1 ( λ t q ) u 0 ( 1 ) S h q , 1 ( λ t q ) u 0 ] ( 1 ) [ t C h q , 2 ( λ t q ) u 1 ( 1 ) t S h q , 2 ( λ t q ) u 1 ] 0 t C h q λ [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( 1 ) S h q λ [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x .
We should formulate the basic assumptions before initiating our main work:
( H 1 )
The function f : [ 0 , T ] × E 1 E 1 is continuous;
( H 2 )
There exists a finite constant L > 0 such that for all t [ 0 , T ] and for all u , v E 1 the inequality
D ( f ( t , u ) , f ( t , v ) ) L D ( u , v )
is valid and such that λ R is satisfied;
( H 3 )
L T q E ( q , q + 1 ) ( | λ | T q ) < 1 .
Theorem 1.
Let λ 1 and suppose that the conditions ( H 1 ) ( H 3 ) are satisfied. Then, the Cauchy problem (1) and (2) has a unique c [(i)-GH]-differentiable solution u in C E [ 0 , T ] .
Proof. 
Let the operator P 1 : C E [ 0 , T ] C E [ 0 , T ] be defined as
P 1 u ( t ) = E q , 1 ( λ t q ) u 0 t E q , 2 ( λ t q ) u 1 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x
It is not difficult to see that u is a c [(i)-GH]-differentiable solution for Cauchy problem (1) and (2) if and only if u = P 1 u . Let u and v belong to E 1 . From the above Lemmas 1 and 2 we infer
D ( P 1 u ( t ) , P 1 v ( t ) ) = D [ 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x , 0 t E q , q ( λ ( t x ) q ) [ f ( x , v ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d s ] L D ( u , v ) 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) , f ( x , v ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x L D ( u , v ) 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) , f ( x , v ( x ) ) ] ( t x ) 1 q d x L D ( u , v ) L 0 t E q , q ( λ ( t x ) q ) f ( u ( x ) , v ( x ) ) ( t x ) 1 q d x L D ( u , v ) L D ( u , v ) 0 t E q , q ( λ ( t x ) q ) ( t x ) 1 q d x = L D ( u , v ) t q E q , q + 1 ( λ t q ) L D ( u , v )
for u , v E 1 , and for all t [ 0 , T ] , which means that
D ( P 1 u , P 1 v ) L [ 1 t q E q , q + 1 ( | λ | T q ) ] D ( u , v )
Thus, the Banach contraction mapping (BCM) principle shows the operator P 1 has a unique fixed point u * C E [ 0 , T ] . It represents the unique c [(i)-GH]-differentiable solution to the Cauchy problem (1) and (2). □
Theorem 2.
Let λ 1 and suppose the conditions ( H 1 ) ( H 3 ) are satisfied. Assume that ( H 4 ) for any t ( 0 , T ] ,
E q , 1 ( λ t q ) u 0 α ̲ + t E q , 2 ( λ t α ) u 1 α ̲ + 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] α ¯ ( t x ) 1 q d x
is non-decreasing in α,
E q , 1 ( λ t q ) u 0 α ¯ + t E q , 2 ( λ t q ) u 1 α ¯ + 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] α ̲ ( t x ) 1 q d x
is non-increasing in α, and for any q [ 1 , 2 ] and t ( 0 , T ]
t E q ( λ t q ) u 1 α + 0 t E q , q ( λ ( t x ) q ) d i a m [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] α ( t x ) 1 q d x t q E q , 2 ( λ t q ) d i a m [ u 1 ] α + t q 1 E q , 1 ( λ t q ) d i a m [ u 0 ] α .
Then, the Cauchy problem (1) and (2) has a unique c [(ii)-GH]-differentiable solution in C E [ 0 , T ] .
Proof. 
Let the operator P 2 : C E [ 0 , T ] C E [ 0 , T ] be defined by
P 2 u ( t ) = E q , 1 ( λ t q ) u 0 E q , 2 ( λ t q ) u 1 ( 1 ) 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x
For condition ( H 4 ) and [35], we know that P 2 is well defined on C E [ 0 , T ] . Moreover, it is not difficult to see that u is a c [(ii)-GH]-differentiable solution for Cauchy problem (1) and (2) if and only if u = P 2 u . Let u and v belong to C E [ 0 , T ] . From the above Lemmas 1 and 2 and Remark 2 we infer
D ( P 2 u ( t ) , P 2 v ( t ) ) L D ( u , v ) ( 1 ) L D ( u , v ) 0 t E q , q ( λ ( t x ) q ) ( t x ) 1 q d x L D ( u , v ) ( 1 ) L D ( u , v ) 0 t E q , q ( | λ | ( t x ) q ) ( t x ) 1 q d x = L D ( u , v ) ( 1 ) t q E q , q + 1 ( | λ | t q ) L D ( u , v )
for u , v E 1 and for all t [ 0 , T ] , which means that
D ( P 2 u ( t ) , P 2 v ( t ) ) L D ( u , v ) ( 1 ) L t q E q , q + 1 ( | λ | t q ) D ( u , v )
Thus, the Banach contraction mapping (BCM) principle shows the operator P 2 has a unique fixed point u * C E [ 0 , T ] . It represents the unique c [(ii)-GH]-differentiable solution to the Cauchy problem (1) and (2). Now, the proof is completed. □
Theorem 3.
Let λ < 1 , and suppose that the conditions ( H 1 ) ( H 3 ) are satisfied. Then, the Cauchy problem (1) and (2) has a c [(i)-GH]-differentiable solution u in C E [ 0 , T ] .
Proof. 
Let the operator P 3 : C E [ 0 , T ] C E [ 0 , T ] be defined as
P 3 u ( t ) = [ C h q , 1 ( λ t q ) u 0 S h q , 1 ( λ t q ) u 0 ] t [ C h q , 2 ( λ t q ) u 1 S h q , 2 ( λ t q ) u 1 ] 0 t C q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] S q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x ,
t [ 0 , T ] . It is not difficult to see that u is a c [(i)-GH]-differentiable solution for Cauchy problem (1) and (2) if and only if u = P 3 u . Let u and v belong to C E [ 0 , T ] . From the above Lemmas 1 and 2 and Remarks 2 and 3 we deduce
D ( P 3 u ( t ) , P 3 v ( t ) ) L D ( u , v ) L D ( u , v ) 0 t E q , q ( λ ( t x ) q ) ( t x ) 1 q d x L D ( u , v ) L D ( u , v ) 0 t E q , q ( | λ | ( t x ) q ) ( t x ) 1 q d x = L D ( u , v ) t q E q , q + 1 ( | λ | t q ) L D ( u , v )
For u , v E 1 and for all t [ 0 , T ] , which signifies as
D ( P 3 u ( t ) , P 3 v ( t ) ) L D ( u , v ) L t q E q , q + 1 ( | λ | t q ) D ( u , v )
Thus, the Banach contraction mapping (BCM) principle shows the operator P 3 has a unique fixed point u * C E [ 0 , T ] . It represents the unique c [(i)-GH]-differentiable solution to the Cauchy problem (1) and (2). Now the proof is done. □
Theorem 4.
Let λ > 0 and suppose that the conditions ( H 1 ) ( H 3 ) are satisfied. Assume that ( H 5 ) for all t ( 0 , T ] the functions
ξ 1 ( t , α ) = C h q , 1 ( λ t q ) u 0 ̲ α + S h q , 1 ( λ t q ) u 0 ¯ α μ 1 ( t , α ) = t C h q , 2 ( λ t q ) u 1 ̲ α + t S h q , 2 ( λ t q ) u 1 ¯ α
is non-decreasing in α. in addition, the function
ξ 2 ( t , α ) = C h q , 1 ( λ t q ) u 0 ¯ q + S h q , 1 ( λ t q ) u 0 ̲ α μ 2 ( t , α ) = t C h q , 2 ( λ t q ) u 1 ¯ q + t S h q , 2 ( λ t q ) u 1 ̲ α
are non-increasing in α. Furthermore, assume ( H 6 ) for all t ( 0 , T ] , the function
ψ 1 ( t , x , α ) = C h q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] α ̲ + S h q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] α ¯
is non-decreasing in α. In addition, the function
ψ 2 ( t , x , α ) = C h q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] α ¯ + S h q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] α ̲
is non-increasing in α. In addition, the function ( H 7 ) for all t ( 0 , T ]
ξ 1 ( t , α ) + μ 1 ( t , α ) + 0 t ψ 1 ( t , x , α ) ( t x ) 1 q d x
is non-decreasing in α, the expression
ξ 2 ( t , α ) + μ 2 ( t , α ) + 0 t ψ 2 ( t , x , α ) ( t x ) 1 q d x
is non-increasing in α, and for all q ( 1 , 2 ] and t ( 0 , T ] ,
0 t E q , q ( λ ( t x ) q ) d i a m [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] α ( t x ) 1 q d x t q 1 E q , q ( λ t q ) d i a m [ u 0 ] q
Then, the Cauchy problem (1) and (2) has a unique c [(ii)-GH]-differentiable solution for C E [ 0 , T ] .
Proof. 
Let the operator P 4 : C E [ 0 , T ] C E [ 0 , T ] be defined as
P 4 u ( t ) = [ C h q , 1 ( λ t q ) u 0 ( 1 ) S h q , 1 ( λ t q ) u 0 ] ( 1 ) t [ C h q , 2 ( λ t q ) u 1 ( 1 ) S h q , 2 ( λ t q ) u 1 ] ( 1 ) 0 t C q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( 1 ) S q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x
According to conditions ( H 5 ) ( H 7 ) and [21], it is known that P 4 is well illustrated on C E [ 0 , T ] . From the above Lemmas 1, 2, and Remark 2,
D ( P 4 u ( t ) , P 4 v ( t ) ) L D ( u , v ) 0 t E q , q ( λ ( t x ) q ) ( t x ) 1 q d x D ( P 4 u ( t ) , P 4 v ( t ) ) = t q E q , q + 1 ( λ t q ) L D ( u , v )
for u , v E 1 and for all t [ 0 , T ] , which means that
D ( P 4 u ( t ) , P 4 v ( t ) ) L t q E q , q + 1 ( | λ | t q ) D ( u , v )
Thus, the Banach contraction mapping (BCM) principle shows the operator P 4 has a unique fixed point u * C E [ 0 , T ] . It represents the unique c [(ii)-GH]-differentiable solution to the Cauchy problem (1) and (2). □

4. Stability Results

In various studies, E α –Ulam-type stability approaches regarding fractional differential equations [23] and Ulam-type stability approaches regarding fuzzy differential equations [24,25] were established. Afterward, Yupin Wang and Shurong Sun worked on E q –Ulam-type stability concepts regarding fuzzy fractional differential equation where q ( 0 , 1 ] . We offer some new E q –Ulam-type stability concepts regarding fuzzy fractional differential equation where q ( 1 , 2 ] .
Assume that γ > 0 is a constant and that t ζ ( t ) , t [ 0 , T ] is a positive continuous function. In addition, suppose that t u ( t ) , t [ 0 , T ] is a continuous function that solves the equation in (1) and consider the following related inequalities:
D ( 0 c D t q u ( t ) , λ u ( t ) ( f ( t , u ( t ) ) B ( t ) C ( t ) ) γ ,
D ( 0 c D t q u ( t ) , λ u ( t ) ( f ( t , u ( t ) ) B ( t ) C ( t ) ) ζ ( t ) ,
D ( 0 c D t q u ( t ) , λ u ( t ) ( f ( t , u ( t ) ) B ( t ) C ( t ) ) γ ζ ( t ) ,
where t [ 0 , T ] .
Definition 3.
Equation (1) is called E q –Ulam–Hyers stable in case there exist a finite constant c > 1 and a function v C E [ 0 , T ] that satisfies the equation in (1) such that for all γ > 1 and for all solutions u C E [ 0 , T ] of Equation (1) that satisfy the inequality in (3), the following inequality is valid:
D ( u ( t ) , v ( t ) ) c E q ( ξ f t q ) γ , ξ f 1 , t [ 0 , T ] .
Definition 4.
Equation (1) is called E q –Ulam–Hyers stable in case there exist a continuous function ϑ : R + R + with ϑ ( 1 ) = 1 and a function v C E [ 0 , T ] that satisfies the equation in (1) and for all solutions u C E [ 0 , T ] of Equation (1) that satisfy the inequality in (3), the following inequality is valid:
D ( u ( t ) , v ( t ) ) ϑ ( γ ) E q ( ξ f t q ) , ξ f 1 , t [ 0 , T ] .
Definition 5.
Equation (1) is called E q –Ulam–Hyers–Rassias stable in case with respect to ζ, when there exist c ζ > 1 and a function v C E [ 0 , T ] that satisfies the equation in (1) such that for all γ > 1 and for all solutions u C E [ 0 , T ] of the equation in (1) that satisfy the inequality in (5), the following inequality is valid:
D ( u ( t ) , v ( t ) ) c ζ γ ζ ( t ) E q ( ξ f t q ) , ξ f 1 , t [ 0 , T ] .
Definition 6.
Equation (1) is called E q –Ulam–Hyers–Rassias stable in case with respect to ζ if there exist c ζ > 1 and a function v C E [ 0 , T ] that satisfies the equation in (1) such that for all γ > 1 and for all solutions u C E [ 0 , T ] of the equation in (1) that satisfy the inequality in (4), the following inequality is valid:
D ( u ( t ) , v ( t ) ) c ζ ζ ( t ) E q ( ξ f t q ) , ξ f 1 , t [ 0 , T ] .
Lemma 5.
The function u C E [ 0 , T ] with the property that ( H 8 ) 0 c D t q u ( t ) [ λ u ( t ) ( f ( t , u ( t ) ) B ( t ) C ( t ) ) ] exists in E 1 for all t [ 0 , T ] satisfies the inequality (3) if and only if there exists a function h C E [ 0 , T ] such that
(i) D h ( t ) , 0 ^ γ , for all t ( 0 , T ] ,
and the function u C E [ 0 , T ] itself satisfies
(ii) 0 c D t q u ( t ) = λ u ( t ) ( f ( t , u ( t ) ) B ( t ) C ( t ) ) , for all t ( 0 , T ] .
Proof. 
The sufficiency begins obviously, and we will only prove the necessity. From condition ( H 8 ) , we observe that the function t h ( t ) , t [ 0 , T ] , defined by
h ( t ) = 0 c D t q u ( t ) [ λ u ( t ) ( f ( t , u ( t ) ) B ( t ) C ( t ) ) ]
belongs to C E [ 0 , T ] and that h(t) belongs to E 1 for all t ( 0 , T ] . Therefore, it follows that the equation in (ii) is satisfied. Additionally, we have
D ( 0 c D t q u ( t ) , λ u ( t ) ( f ( t , u ( t ) ) B ( t ) C ( t ) ) ) = D 0 c D t q u ( t ) [ λ u ( t ) f ( t , u ( t ) ) B ( t ) C ( t ) ] , 0 ^ = D h ( t ) , 0 ^ .
From the inequality (3), it then follows that D h ( t ) , 0 ^ γ , and therefore, (i) is satisfied. This completes the proof of Lemma 5. □
Remark 5.
Similar results as in Lemma 5 can be obtained by using the inequalities in (4) and (5).
Lemma 6.
Let u ( t ) be a c [ ( i ) -GH]-differentiable function that solves the Cauchy equation in (1) and (2) and satisfies the inequality in (3) and is such that 0 c D t q u ( 0 ) = u 0 . Let the condition in ( H 8 ) be satisfied. Then, for every t [ 0 , T ] , the function u ( t ) satisfies the inequality
D ( u ( t ) , G 1 ( f , t ) ) γ E q , q ( | λ | t q ) 0 t ζ ( x ) ( t x ) 1 q d x
when λ > 1 , and
D ( u ( t ) , G 2 ( f , t ) ) γ E q , q ( | λ | t q ) 0 t ζ ( x ) ( t x ) 1 q d x
when λ < 1 , and t [ 0 , T ] . Here, the functions G 1 ( f , t ) and G 2 ( f , t ) are defined by
G 1 ( f , t ) = E q , 1 ( λ t q ) u 0 t E q , 2 ( λ t q ) u 1 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x G 2 ( f , t ) = [ C h q , 1 ( λ t q ) u 0 S h q , 1 ( λ t q ) u 0 ] [ t C h q , 2 ( λ t q ) u 1 t S h q , 2 ( λ t q ) u 1 ] 0 t C q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] S q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x
Proof. 
Since the function u C E [ 0 , T ] is a solution to the Cauchy problem (1) and (2), we infer
0 c D t α u ( t ) = λ u ( t ) [ f ( t , u ( t ) ) + B ( t ) C ( t ) ] , u ( 0 ) = u 0 u ( 0 ) = u 1 .
Now, regarding clarity, the proof can be divided into two cases.
Case 1.
Suppose λ > 1 . Then, we write
C 1 ( B , C , t ) = 0 t E q , q ( λ ( t x ) q ) B ( x ) C ( x ) ( t x ) 1 q d x
Observing that u is a c [(i)-GH]-differentiable solution of Equation (6), then Lemma 4 with f ( t , u ( t ) ) + B ( t ) C ( t ) instead of f ( t ) shows the equality
u ( t ) = E q , 1 ( λ t q ) u 0 t E q , 2 ( λ t q ) u 1 0 t E q , q ( λ ( t x ) q ) ( f ( x , u ( x ) ) + B ( x ) C ( x ) ) ( t x ) 1 q d x = E q , 1 ( λ t q ) u 0 t E q , 2 ( λ t q ) u 1 0 t E q , q ( λ ( t x ) q ) f ( x , u ( x ) ) ( t x ) 1 q d x E q , q ( λ ( t x ) q ) B ( x ) C ( x ) ( t x ) 1 q d x = G 1 ( f , t ) C 1 ( B , C , t )
Now, it follows that
D ( u ( t ) , G 1 ( f , t ) ) = D ( u ( t ) C 1 ( B , C , t ) , G 1 ( f , t ) C 1 ( B , C , t ) ) = D ( u ( t ) C 1 ( B , C , t ) , u ( t ) ) = D ( C 1 ( B , C , t ) , 0 ^ ) = D 0 t E q , q ( λ ( t x ) q ) B ( x ) C ( x ) , 0 ^ ( t x ) 1 q d x = 0 t E q , q ( λ ( t x ) q ) D ( B ( x ) C ( x ) , 0 ^ ) ( t x ) 1 q d x γ E q , q ( | λ | t q ) 0 t ζ ( x ) ( t x ) 1 q d x
Case 2.
When λ < 1 , we denote
C 2 ( B , C , t ) = 0 t C q λ ( t , x ) B ( x ) C ( x ) S q λ ( t , x ) B ( x ) C ( x ) ( t x ) 1 q d x
It should be observed that u ( t ) is a c [(ii)-GH]-differentiable solution of Equation (6) that obeys the inequality in (3). An application of Lemma 5 then yields
u ( t ) = [ C h q , 1 ( λ t q ) u 0 S h q , 1 ( λ t q ) u 0 ] [ t C h q , 2 ( λ t q ) u 1 t S h q , 2 ( λ t q ) u 1 ] 0 t C q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] S q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x = G 2 ( f , t ) C 2 ( B , C , t ) .
Now, it follows that
D ( u ( t ) , G 2 ( f , t ) ) = D ( C 2 ( B , C , t ) , 0 ^ ) 0 t D ( C q λ ( t , x ) [ B ( x ) C ( x ) ] S q λ ( t , x ) [ B ( x ) C ( x ) ] , 0 ^ ) ( t x ) 1 q d x 0 t E q , q ( | λ | t x q ) D ( B ( x ) C ( x ) ) , 0 ^ ) ( t x ) 1 q d x γ E q , q ( | λ | t q ) 0 t ζ ( x ) ( t x ) 1 q d x
Now, the proof is completed. □
Lemma 7.
Let u ( t ) be a c [(ii)-GH]-differentiable function that solves the Cauchy equation in (1) and (2) and satisfies the inequality (3) and is such that 0 c D t q u ( 0 ) = u 0 . Let the condition in ( H 8 ) be satisfied. Then, for every t [ 0 , T ] the function u ( t ) satisfies the integral inequality
D ( u ( t ) , G 3 ( f , t ) ) γ E q , q ( | λ | t q ) 0 t ζ ( x ) ( t x ) 1 q d x
when λ < 1 , and
D ( u ( t ) , G 4 ( f , t ) ) γ E q , q ( | λ | t q ) 0 t ζ ( x ) ( t x ) 1 q d x
when λ > 1 and t [ 0 , T ] . Here, the functions G 3 ( f , t ) and G 4 ( f , t ) are defined by
G 3 ( f , t ) = E q , 1 ( λ t q ) u 0 ( 1 ) t E q , 2 ( λ t q ) u 1 ( 1 ) 0 t E q , q ( λ ( t x ) q ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x G 4 ( f , t ) = [ C h q , 1 ( λ t q ) u 0 ( 1 ) S h q , 1 ( λ t q ) u 0 ] ( 1 ) [ t C h q , 2 ( λ t q ) u 1 ( 1 ) t S h q , 2 ( λ t q ) u 1 ] ( 1 ) 0 t C q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( 1 ) S q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x
Proof. 
Now, regarding clarity, the proof can be divided into two cases.
Case 1.
When λ < 1 , observe that u is a c [(ii)-GH]-differentiable solution of Equation (5), then the Lemma 5 with f ( x , u ( x ) ) + B ( x ) C ( x ) instead of f ( t ) shows the equality
u ( t ) = G 3 ( f , t ) ( 1 ) C 1 ( B , C , t ) .
Now, it follows that
D ( u ( t ) , G 3 ( f , t ) ) = D ( 0 ^ , C 1 ( g , t ) ) 0 t E q , q ( | λ | ( t x ) q ) D ( B ( x ) C ( x ) , 0 ^ ) ( t x ) 1 q d x γ E q , q ( | λ | t q ) 0 t ζ ( t ) ( t x ) 1 q d x
Case 2.
Suppose λ > 1 . Then, we denote
C 3 ( B , C , t ) = 0 t C q λ ( t , x ) B ( x ) C ( x ) ( 1 ) S q λ ( t , x ) B ( x ) C ( x ) ( t x ) 1 q d x
Observing that u ( t ) is a c [(ii)-GH]-differentiable solution of Equation (5), then Lemma 5 with f ( t , u ( t ) ) + B ( t ) C ( t ) instead of f ( t ) shows the equality
u ( t ) = [ C h q , 1 ( λ t q ) u 0 ( 1 ) S h q , 1 ( λ t q ) u 0 ] ( 1 ) t [ C h q , 2 ( λ t q ) u 1 S h q , 2 ( λ t q ) u 1 ] ( 1 ) 0 t C q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] S q λ ( t , x ) [ f ( x , u ( x ) ) + B ( x ) C ( x ) ] ( t x ) 1 q d x = G 4 ( f , t ) ( 1 ) C 3 ( B , C , t ) .
Now, it follows that
D ( u ( t ) , G 4 ( f , t ) ) = D ( C 3 ( B , C , t ) , 0 ^ ) 0 t D ( C q λ ( t , x ) [ B ( x ) C ( x ) ] ( 1 ) S q λ ( t , x ) [ B ( x ) C ( x ) ] , 0 ^ ) ( t x ) 1 q d x 0 t E q , q ( | λ | t x q ) D ( B ( x ) C ( x ) ) , 0 ^ ) ( t x ) 1 q d x γ E q , q ( | λ | t q ) 0 t ζ ( x ) ( t x ) 1 q d x
Now, the proof is completed. □
Remark 6.
We can obtain similar results to those in Lemmas 6 and 7 for inequalities (3) and (4).
Theorem 5.
Suppose λ 1 , condition ( H 1 ) ( H 3 ) are satisfied, and the following condition holds ( H 9 ) ; there exists a positive, increasing, and continuous function ζ such that
E q , q ( | λ | t q ) 0 t ζ ( x ) ( t x ) 1 q d x c ζ ζ ( t ) , t [ 0 , T ] .
Assume further that c [ ( i ) -GH]-differentiable function u satisfied the inequality (5) with the function ζ in ( H 9 ) and that u satisfies condition ( H 8 ) . Then, Equation (1) is E q –Ulam–Hyers–Rassias stable.
Proof. 
According to Theorem 1, u is a c [(i)-GH]-differentiable solution to Cauchy problem (1) and (2). Let u be a c [(i)-GH]-differentiable solution to Equation (1), which satisfies inequality (5) with u ( 0 ) = u 0 . From Lemma 6, we obtain
D ( u ( t ) , G 1 ( f , t ) ) γ E q , q ( | λ | t q ) 0 t ζ ( x ) ( t x ) 1 q d x c ζ γ ζ ( t ) ,
t [ 0 , T ] . According to condition ( H 9 ) , it follows that
D ( u ( t ) , v ( t ) ) D ( u ( t ) , G 1 ( f , t ) ) + D ( G 1 ( f , t ) , v ( t ) ) c ζ γ ζ ( t ) + 0 t E q , q ( λ ( t x ) q ) D [ ( f ( x , u ( x ) ) + B ( x ) C ( x ) ) , ( f ( x , v ( x ) ) + B ( x ) C ( x ) ) ] ( t x ) 1 q d x c ζ γ ζ ( t ) + L 0 t E q , q ( λ ( t s ) q ) D [ u ( x ) , v ( x ) ] ( t x ) 1 q d x c ζ γ ζ ( t ) + L E q , q ( λ ( t x ) q ) 0 t ) D [ u ( x ) , v ( x ) ] ( t x ) 1 q d x
By the generalized Gronwall inequality [36], we obtain
D ( u ( t ) , v ( t ) ) c ζ γ ζ ( t ) E q ( L E q , q ( | λ | T q ) Γ ( q ) t q ) .
Thus, Equation (1) is E q –Ulam–Hyers–Rassias stable in view of Definition 5. □
Theorem 6.
Let λ 1 and let the condition ( H 1 ) ( H 4 ) , ( H 8 ) and ( H 9 ) hold for a c [(ii)-GH]-differentiable function u satisfy inequality (5). Then, Equation (1) is E q –Ulam–Hyers–Rassias stable.
Proof. 
According to Theorem 2, u is a c [ ( i i ) -GH]-differentiable solution to Cauchy problem (1) and (2). Let u be a c [(ii)-GH]-differentiable solution to Equation (1), which satisfies inequality (5) with u ( 0 ) = u 0 . From Lemma 7, we obtain
D ( u ( t ) , G 3 ( f , t ) ) c ζ γ ζ ( t ) ,
t [ 0 , T ] . According to condition ( H 9 ) it follows that
D ( u ( t ) , v ( t ) ) D ( u ( t ) , G 3 ( f , t ) ) + D ( G 3 ( f , t ) , v ( t ) ) c ζ γ ζ ( t ) + 0 t E q , q ( | λ | ( t x ) q ) D [ ( f ( x , u ( x ) ) + B ( x ) C ( x ) ) , ( f ( x , v ( x ) ) + B ( x ) C ( x ) ) ] ( t x ) 1 q d x c ζ γ ζ ( t ) + L 0 t E q , q ( | λ | ( t x ) q ) D [ u ( x ) , v ( x ) ] ( t x ) 1 q d x c ζ γ ζ ( t ) + L E q , q ( | λ | ( t x ) q ) 0 t ) D [ u ( x ) , v ( x ) ] ( t x ) 1 q d x
By the generalized Gronwall inequality, we obtain
D ( u ( t ) , v ( t ) ) c ζ γ ζ ( t ) E q ( L E q , q ( | λ | T q ) Γ ( q ) t q ) .
Thus, Equation (1) is E q –Ulam–Hyers–Rassias stable in view of Definition 5. □
Theorem 7.
Let λ < 1 , and let the condition ( H 1 ) ( H 3 ) , ( H 8 ) and ( H 9 ) hold for a c [(i)-GH]-differentiable function u satisfies inequality (5). Then Equation (1) is E q –Ulam–Hyers–Rassias stable.
Proof. 
According to Theorem 3, u is a c [(i)-GH]-differentiable solution to Cauchy problem (1) and (2). Let u be a c [(i)-GH]-differentiable solution to Equation (1), which satisfies inequality (5) with u ( 0 ) = u 0 . From Lemma 6, we obtain
D ( u ( t ) , G 2 ( f , t ) ) c ζ γ ζ ( t ) ,
t [ 0 , T ] . According to condition ( H 9 ) , it follows that
D ( u ( t ) , v ( t ) ) D ( u ( t ) , G 2 ( f , t ) ) + D ( G 2 ( f , t ) , v ( t ) ) c ζ γ ζ ( t ) + 0 t E q , q ( | λ | ( t x ) q ) D [ ( f ( x , u ( x ) ) + B ( x ) C ( x ) ) , ( f ( x , v ( x ) ) + B ( x ) C ( x ) ) ] ( t x ) 1 q d x c ζ γ ζ ( t ) + L 0 t E q , q ( | λ | ( t x ) q ) D [ u ( x ) , v ( x ) ] ( t x ) 1 q d x c ζ γ ζ ( t ) + L E q , q ( | λ | ( t x ) q ) 0 t ) D [ u ( x ) , v ( x ) ] ( t x ) 1 q d x
By the generalized Gronwall inequality, we obtain
D ( u ( t ) , v ( t ) ) c ζ γ ζ ( t ) E q ( L E q , q ( | λ | T q ) Γ ( q ) t q ) .
Thus, Equation (1) is E q –Ulam–Hyers–Rassias stable in view of Definition 5. □
Theorem 8.
Let λ > 1 , let the condition ( H 1 ) ( H 3 ) as well as ( H 5 ) ( H 7 ) , ( H 8 ) ( H 9 ) hold for a c [(ii)-GH]-differentiable function u, which satisfies inequality (5). Then, Equation (1) is E q –Ulam–Hyers–Rassias stable.
Proof. 
According to Theorem 4, set u is a c [(ii)-GH]-differentiable solution to Cauchy problem (1) and (2). let u be a c [(ii)-GH]-differentiable solution to Equation (1), which satisfies the inequality (5) with u ( 0 ) = u 0 . from Lemma 7, we obtain
D ( u ( t ) , G 4 ( f , t ) ) c ζ γ ζ ( t ) ,
t [ 0 , T ] . According to condition ( H 9 ) , it follows that
D ( u ( t ) , v ( t ) ) D ( u ( t ) , G 4 ( f , t ) ) + D ( G 4 ( f , t ) , v ( t ) ) c ζ γ ζ ( t ) + 0 t E q , q ( | λ | ( t x ) q ) D [ ( f ( x , u ( x ) ) + B ( x ) C ( x ) ) , ( f ( x , v ( x ) ) + B ( x ) C ( x ) ) ] ( t x ) 1 q d x c ζ γ ζ ( t ) + L 0 t E q , q ( | λ | ( t x ) q ) D [ u ( x ) , v ( x ) ] ( t x ) 1 q d x c ζ γ ζ ( t ) + L E q , q ( | λ | ( t x ) q ) 0 t ) D [ u ( x ) , v ( x ) ] ( t x ) 1 q d x
By the generalized Gronwall inequality, we obtain
D ( u ( t ) , v ( t ) ) c ζ γ ζ ( t ) E q ( L E q , q ( | λ | T q ) Γ ( q ) t q ) .
Thus, Equation (1) is E q –Ulam–Hyers–Rassias stable in view of Definition 5. Now, the proof is completed. □
Remark 7.
In view of Definition 6 can be verified as according to the assumption in Theorems 5–8, we assume Equation (1) and inequality (4). It can be verified that Equation (1) is generalized E q –Ulam–Hyers–Rassias stable with respect to Definition 6.
Remark 8.
Condition ( H 9 ) weakens 0 t ζ ( x ) d x c ζ ζ ( t ) E 2 ( L E 2 , 2 ( | 1 | T 2 ) Γ ( 2 ) t 2 t [ 0 , T ] when we assume q = 1 . This means that certain theorems in [25] are special cases of Theorem 5 and 6 in the present paper.
Remark 9.
According to the assumptions excluding ( H 9 ) in Theorems 5–8, we consider the equation in (1) and inequality in (3). It can be proved that in terms of Definitions 3 and 4, Equation (1) is E q –Ulam–Hyers.

5. Examples

In this part, we will show four examples to explain our main results.
Example 1.
Consider the following Cauchy problem in terms of a Fuzzy fractional differential equation
0 c D t 1.5 u ( t ) = u ( t ) 1.2 u ( t ) c o s ( t ) t 2 e t F
on ( 0 , 2 π ] , with initial conditions
u ( 0 ) = 0 ^ u ( 0 ) = 1 ^ .
Compared to Equation (1), in the above equations, q = 1.5 , λ = 2 , T = 2 π , f ( t , u ( t ) ) = 1.2 u ( t ) c o s ( t ) t 2 e t F , and F = ( 0 , 1 , 2 ) E 1 is a symmetric triangular fuzzy number. Hence, with L = 1.3 , the condition ( H 1 ) and ( H 2 ) are satisfied. It is not difficult to prove that condition ( H 3 ) is satisfied. Hence, as a consequence of Theorem 1, the Cauchy problem (7) and (8) has a c [(i)-GH]-differentiable solution. The numerical solutions with respect to the q = 1.5 level are provided by utilizing the Adams–Moultan predictor–corrector method.
Furthermore, for ε > 1 , assume that the c [(i)-GH]-differentiable fuzzy-valued function u : ( 0 , 2 π ] E 1 satisfies condition ( H 8 ) and
D ( 0 c D t 1.5 u ( t ) = u ( t ) 1.2 u ( t ) c o s ( t ) t 2 e t F ) ε t
Assuming ζ ( t ) = 1 , t [ 0 , 2 π ] and c ζ = 4 2 π 3 E 1.5 , 1.5 ( 2 π ) , this means that condition ( H 9 ) is satisfied. Hence, Equation (7) is E q –Ulam–Hyers–Rassias stable with respect to Theorem 5.
Example 2.
Consider the following Cauchy problem in terms of a Fuzzy fractional differential equation
0 c D t 1.5 u ( t ) = u ( t ) + t 2 + t + 4
with initial condition
u ( 0 ) = u 0 u ( 0 ) = u 1
Compared to Equation (1), in the above equation, q = 1 , f ( t , u ( t ) ) = t 2 + t + 4 , and u 0 = ( 0 , 1 , 2 ) E 1 is a symmetric triangular fuzzy number.
Hence, with L = 1.3 , the condition ( H 1 ) ( H 2 ) are satisfied. It is not difficult to prove that condition ( H 4 ) is satisfied. Hence, by employing Theorem 2, the Cauchy problem (9)–(10) has a different c [(ii)-GH]-differentiable solution. The numerical solution provides with respect to q-level by utilizing the Adams–Moultan predictor–corrector method.
Furthermore, for ε > 1 , assumes that the c [(ii)-GH]-differentiable fuzzy-valued function u : ( 0 , 2 π ] E 1 satisfies condition ( H 8 ) and
D ( 0 c D t 1.5 u ( t ) , u ( t ) + t 2 + t + 4 ) ε t , t ( 0 , 2 π ]
Assuming ζ ( t ) = t , t [ 0 , 2 π ] and c ζ = 4 2 π 3 E 1.5 , 1.5 ( 2 π ) , this means that condition ( H 9 ) is satisfied. Hence, Equation (9) is E q –Ulam–Hyers–Rassias stable with respect to Theorem 6.
Example 3.
Consider the following Cauchy problem in terms of a Fuzzy fractional differential equation
0 c D t 1.5 u ( t ) = u ( t ) 1.2 u ( t ) c o s ( t ) t 2 e t F
on ( 0 , 2 π ] , with initial conditions
u ( 0 ) = 0 ^ u ( 0 ) = 1 ^
Compared to equations (1), in the above equation, q = 1.5 , λ = 2 , T = 2 π , f ( t , u ( t ) ) = u ( t ) 1.2 u ( t ) c o s ( t ) t 2 e t F and F = ( 0 , 1 , 2 ) E 1 is a symmetric triangular fuzzy number. Hence, with L = 1.3 , it is not difficult to prove that condition ( H 1 ) and ( H 3 ) are satisfied. Hence, as a consequence of Theorem 3, the Cauchy problem (11) and (12) has a c [(i)-GH]-differentiable solution. The numerical solutions with respect to q = 1.5 level are provided by utilizing the Adams–Moultan predictor–corrector method.
Furthermore, for ε > 1 , assume that the c [(i)-GH]-differentiable fuzzy-valued function u : ( 0 , π ] E 1 satisfies condition ( H 8 ) and
D ( 0 c D t 1.5 u ( t ) , u ( t ) 1.2 u ( t ) c o s ( t ) t 2 e t F ε t
Assuming ζ ( t ) = 1 , t [ 0 , 2 π ] and c ζ = 4 2 π 3 E 1.5 , 1.5 ( 2 π ) , this means that condition ( H 9 ) satisfied. Hence, Equation (11) is E q –Ulam–Hyers–Rassias stable with respect to Theorem 7.
Example 4.
Consider the following Cauchy problem in terms of a Fuzzy fractional differential equation
0 c D t 1.5 u ( t ) = u ( t ) + t 2 + t + 4 , t ( 0 , 2 π ]
with initial condition
u ( 0 ) = u 0 u ( 0 ) = u 1
Compared to Equation (1), in the above equations, q = 1.5 , f ( t , u ( t ) ) = t 2 + t + 4 , and u 0 = ( 0 , 1 , 2 ) E 1 is a symmetric triangular fuzzy number.
Hence, with L = 1.3 , the condition ( H 1 ) ( H 3 ) are satisfied. Notice C h 1.5 ( x ) S h 1.5 ( x ) > 1 and C h 1.5 ( x ) S h 1.5 ( x ) < 1 for u ( 0 , 2 π ] . It is not difficult to prove that condition ( H 5 ) ( H 7 ) are satisfied. Hence, as a consequence of Theorem 4, the Cauchy problem (13) and (14) has a unique c [(ii)-GH]-differentiable solution. The numerical solutions with respect to the q = 1.5 level are provided by utilizing the Adams–Moultan predictor–corrector method.
Furthermore, for ε > 1 , assume that the c [(ii)-GH]-differentiable fuzzy-valued function u : ( 0 , 2 π ] E 1 satisfies condition ( H 8 ) and
D ( 0 c D t 1.5 u ( t ) , u ( t ) + t 2 + t + 4 ) ε t , t ( 0 , 2 π ]
Assuming ζ ( t ) = t , t [ 0 , 2 π ] and c ζ = 4 2 π 3 E 1.5 , 1.5 ( 2 π ) , this means that condition ( H 9 ) satisfied. Hence, Equation (13) is E q –Ulam–Hyers–Rassias stable with respect to Theorem 8.

6. Graphical Presentation

We used the Adams–Bashforth–Moulton technique to acquire the numerical solution for this fractional differential equation for graphical representation of the solution of the problem presented in Equations (7), (9), (11) and (13). For simulation, the modified predictor–corrector scheme is used to examine the effect and contribution of the time-delayed factor. A graphical representation of the solution with different variations of the time delay factor, as well as other parameters, is made to check and demonstrate the stability of the model under consideration. We are able to see the Ulam–Hyers stability of varied accuracies and delays from the numerical data. The system will attain Ulam–Hyers stability more quickly with greater accuracy. This is also true when the number of delays increases. Figure 1, Figure 2, Figure 3 and Figure 4 show the stability of the system (7), (9), (11) and (13) for various time delays and fractional derivatives.

7. Conclusions

This paper aims to define the uniqueness and existence of a group of nonlinear fuzzy fractional differential equation of solutions to the Cauchy problem. Moreover, E q –Ulam-type stability of Equation (1) is observed by applying the inequality technique. We obtain uniqueness and existence results with the help of nonlocal conditions of the Caputo derivative. Moreover, future work may include broadening the idea indicated in this task and familiarizing observability, and generalize other tasks. Ulam-type stability of fuzzy fractional differential equations, similar to crisp situations for approximate solutions, provides a reliable theoretical basis. This a fruitful area with wide research projects, and it can bring about countless applications and theories. We have decided to devote much attention to this area. Furthermore, it is fruitful to investigate stability problems in a classical sense for the fuzzy fractional differential equation.

Author Contributions

A.U.K.N., J.H., R.S. and B.A. contributed equally to the writing of this paper. All authors studied and validated the final document. All authors have read and agreed to the published version of the manuscript.

Funding

The work was of Azmat Ullah Khan Niazi, and supported by post doctoral funding of Xiangtan University, Hunan, China.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All the data is present within the manuscript.

Acknowledgments

The authors would like to express their sincere gratitude to the reviewers and the editors for their careful reviews and constructive recommendations. Jiawei He was supported by the Hunan key Laboratory for Computation and Simulation in Science and Engineering Xiangtan University.

Conflicts of Interest

All authors have no conflict of interest.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley: New York, NY, USA, 1993. [Google Scholar]
  2. Oldham, K.; Spanier, J. The Fractional Calculus Theory and Applications of Differentiation and Integration to Arbitrary Order; Elsevier: Amsterdam, The Netherlands, 1974. [Google Scholar]
  3. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives (Theory and Applications); Gordon and Breach Science Publishers: Amsterdam, The Netherlands, 1993. [Google Scholar]
  4. Ahmad, B.; Nieto, J.J. Existence results for a coupled system of nonlinear fractional differential equations with three-point boundary conditions. Comput. Math. Appl. 2009, 58, 1838–1843. [Google Scholar] [CrossRef] [Green Version]
  5. Ahmad, B.; Ntouyas, S.K.; Agarwal, R.P.; Alsaedi, A. On fractional differential equations and inclusions with nonlocal and average-valued (integral) boundary conditions. Adv. Differ. Equ. 2016, 2016, 80. [Google Scholar] [CrossRef] [Green Version]
  6. Ding, Z.; Ma, M.; Kandel, A. Existence of the solutions of fuzzy differential equations with parameters. Inf. Sci. 1997, 99, 205–217. [Google Scholar] [CrossRef]
  7. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
  8. Lakshmikantham, V.; Vatsala, A.S. Basic theory of fractional differential equations. Nonlinear Anal. Theory Methods Appl. 2008, 69, 2677–2682. [Google Scholar] [CrossRef]
  9. Podlubny, I. Fractional Differential Equations: An Introduction to Fractional Derivatives, Fractional Differential Equations, to Methods of Their Solution and Some of Their Applications; Elsevier: Amsterdam, The Netherlands, 1998. [Google Scholar]
  10. Sakulrang, S.; Moore, E.J.; Sungnul, S.; de Gaetano, A. A fractional differential equation model for continuous glucose monitoring data. Adv. Differ. Equ. 2017, 2017, 150. [Google Scholar] [CrossRef]
  11. Agarwal, R.; Lakshmikantham, V.; Nieto, J. On the concept of solution for fractional differential equations with uncertainty. Nonlinear Anal. 2010, 72, 2859–2862. [Google Scholar] [CrossRef]
  12. Allahviranloo, T.; Salahshour, S.; Abbasbandy, S. Explicit solutions of fractional differential equations with uncertainty. Soft Comput. 2012, 16, 297–302. [Google Scholar] [CrossRef]
  13. Salahshour, S.; Allahviranloo, T.; Abbasbandy, S. Solving fuzzy fractional differential equations by fuzzy Laplace transforms. Commun. Nonlinear Sci. Numer. Simul. 2012, 17, 1372–1381. [Google Scholar] [CrossRef]
  14. Bede, B.; Gal, S.G. Almost periodic fuzzy-number-valued functions. Fuzzy Sets Syst. 2004, 147, 385–403. [Google Scholar] [CrossRef]
  15. Bede, B.; Gal, S.G. Generalizations of the differentiability of fuzzy-number-valued functions with applications to fuzzy differential equations. Fuzzy Sets Syst. 2005, 151, 581–599. [Google Scholar] [CrossRef]
  16. Chehlabi, M.; Allahviranloo, T. Concreted solutions to fuzzy linear fractional differential equations. Appl. Soft Comput. 2016, 44, 108–116. [Google Scholar] [CrossRef]
  17. Allahviranloo, T.; Armand, A.; Gouyandeh, Z. Fuzzy fractional differential equations under generalized fuzzy Caputo derivative. J. Intell. Fuzzy Syst. 2014, 26, 1481–1490. [Google Scholar] [CrossRef]
  18. Armand, A.; Allahviranloo, T.; Abbasbandy, S.; Gouyandeh, Z. Fractional relaxation-oscillation differential equations via fuzzy variational iteration method. J. Intell. Fuzzy Syst. 2017, 32, 363–371. [Google Scholar] [CrossRef]
  19. Wang, Y.; Sun, S.; Han, Z. Existence of solutions to periodic boundary value problems for fuzzy fractional differential equations. Int. J. Dyn. Syst. Differ. Equ. 2017, 7, 195–216. [Google Scholar] [CrossRef]
  20. Agarwal, R.P.; Baleanu, D.; Nieto, J.J.; Torres, D.F.; Zhou, Y. A survey on fuzzy fractional differential and optimal control nonlocal evolution equations. J. Comput. Appl. Math. 2018, 339, 3–29. [Google Scholar] [CrossRef] [Green Version]
  21. Huang, L.L.; Baleanu, D.; Mo, Z.W.; Wu, G.C. Fractional discrete-time diffusion equation with uncertainty: Applications of fuzzy discrete fractional calculus. Phys. A Stat. Mech. Its Appl. 2018, 508, 166–175. [Google Scholar] [CrossRef]
  22. Wang, Y.; Sun, S.; Han, Z. On fuzzy fractional Schrödinger equations under Caputo’s H-differentiability. J. Intell. Fuzzy Syst. 2018, 34, 3929–3940. [Google Scholar] [CrossRef]
  23. Wang, J.; Li, X. Eα-Ulam type stability of fractional order ordinary differential equations. J. Appl. Math. Comput. 2014, 45, 449–459. [Google Scholar] [CrossRef]
  24. Shen, Y. On the Ulam stability of first order linear fuzzy differential equations under generalized differentiability. Fuzzy Sets Syst. 2015, 280, 27–57. [Google Scholar] [CrossRef]
  25. Wang, Y.; Sun, S. Existence, uniqueness and Eq-Ulam type stability of fuzzy fractional differential equations with parameters. J. Intell. Fuzzy Syst. 2019, 36, 5533–5545. [Google Scholar] [CrossRef]
  26. Shen, Y. Hyers-Ulam-Rassias stability of first order linear partial fuzzy differential equations under generalized differentiability. Adv. Differ. Equ. 2015, 2015, 351. [Google Scholar] [CrossRef] [Green Version]
  27. Wang, C.; Xu, T.Z. Hyers-Ulam stability of fractional linear differential equations involving Caputo fractional derivatives. Appl. Math. 2015, 60, 383–393. [Google Scholar] [CrossRef] [Green Version]
  28. Brzdęk, J.; Eghbali, N. On approximate solutions of some delayed fractional differential equations. Appl. Math. Lett. 2016, 54, 31–35. [Google Scholar] [CrossRef]
  29. Lakshmikantham, V.; Mohapatra, R.N. Theory of Fuzzy Differential Equations and Inclusions; CRC Press: Boca Raton, FL, USA, 2004. [Google Scholar]
  30. Lakshmikantham, V.; Bhaskar, T.G.; Devi, J.V. Theory of Set Differential Equations in Metric Spaces; Cambridge Scientific Publishers: Cambridge, UK, 2006. [Google Scholar]
  31. Gorenflo, R.; Kilbas, A.A.; Mainardi, F.; Rogosin, S.V. Mittag-Leffler Functions, Related Topics and Applications; Springer: Berlin, Germany, 2014; Volume 2. [Google Scholar]
  32. Wang, J.; Feckan, M.; Zhou, Y. Presentation of solutions of impulsive fractional Langevin equations and existence results. Eur. Phys. J. Spec. Top. 2013, 222, 1857–1874. [Google Scholar] [CrossRef]
  33. Peng, S.; Wang, J. Cauchy problem for nonlinear fractional differential equations with positive constant coefficient. J. Appl. Math. Comput. 2016, 51, 341–351. [Google Scholar] [CrossRef]
  34. Melliani, S.; El Allaoui, A.; Chadli, L.S. Relation between fuzzy semigroups and fuzzy dynamical systems. Nonlinear Dyn. Syst. Theory 2017, 17, 60–69. [Google Scholar]
  35. Bede, B.; Stefanini, L. Generalized differentiability of fuzzy-valued functions. Fuzzy Sets Syst. 2013, 230, 119–141. [Google Scholar] [CrossRef] [Green Version]
  36. Ye, H.; Gao, J.; Ding, Y. A generalized Gronwall inequality and its application to a fractional differential equation. J. Math. Anal. Appl. 2007, 328, 1075–1081. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Solution of Problems (7) and (8).
Figure 1. Solution of Problems (7) and (8).
Fractalfract 05 00066 g001
Figure 2. Solution of Problems (9) and (10).
Figure 2. Solution of Problems (9) and (10).
Fractalfract 05 00066 g002
Figure 3. Solution of Problems (11) and (12).
Figure 3. Solution of Problems (11) and (12).
Fractalfract 05 00066 g003
Figure 4. Solution of Problems (13) and (14).
Figure 4. Solution of Problems (13) and (14).
Fractalfract 05 00066 g004
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Niazi, A.U.K.; He, J.; Shafqat, R.; Ahmed, B. Existence, Uniqueness, and Eq–Ulam-Type Stability of Fuzzy Fractional Differential Equation. Fractal Fract. 2021, 5, 66. https://doi.org/10.3390/fractalfract5030066

AMA Style

Niazi AUK, He J, Shafqat R, Ahmed B. Existence, Uniqueness, and Eq–Ulam-Type Stability of Fuzzy Fractional Differential Equation. Fractal and Fractional. 2021; 5(3):66. https://doi.org/10.3390/fractalfract5030066

Chicago/Turabian Style

Niazi, Azmat Ullah Khan, Jiawei He, Ramsha Shafqat, and Bilal Ahmed. 2021. "Existence, Uniqueness, and Eq–Ulam-Type Stability of Fuzzy Fractional Differential Equation" Fractal and Fractional 5, no. 3: 66. https://doi.org/10.3390/fractalfract5030066

APA Style

Niazi, A. U. K., He, J., Shafqat, R., & Ahmed, B. (2021). Existence, Uniqueness, and Eq–Ulam-Type Stability of Fuzzy Fractional Differential Equation. Fractal and Fractional, 5(3), 66. https://doi.org/10.3390/fractalfract5030066

Article Metrics

Back to TopTop