Exp-type Ulam-Hyers stability of fractional differential equations with

Gao et al. Advances in Difference Equations (2015) 2015:238 DOI 10.1186/s13662-015-0578-4

RESEARCH

Open Access

Exp-type Ulam-Hyers stability of fractional differential equations with positive constant coefficient Zhuoyan Gao1 , Xiulan Yu1 and JinRong Wang2* *

Correspondence: [email protected] 2 Department of Mathematics, Guizhou University, Guiyang, Guizhou 550025, P.R. China Full list of author information is available at the end of the article

Abstract In this paper, we apply asymptotic behavior on Mittag-Leffler functions Eα (z) and Eα ,α (z) for z > 0 to discuss exp-type Ulam-Hyers stability of c Dαt x(t) = λx(t) + f (t, x(t)) for the case λ > 0 on a finite time interval [0, 1] and an unbounded interval (1, ∞). MSC: exp-type Ulam-Hyers stability; fractional differential equations; asymptotic behavior Keywords: 26A33; 34D10; 45N05

1 Introduction Since SM Ulam raised the famous question on the stability of functional equations, more and more researchers generalize, improve, and extend Ulam’s type stability problems of differential equations. For more contributions on such area, one can refer to [–] and the references therein. Recently, the authors [] discussed Ulam’s type stability of fractional differential equations c

Dαt x(t) = λx(t) + f t, x(t) ,

α ∈ (, ), t ∈ J, λ < ,

()

on finite time interval J, where c Dαt is the Caputo fractional derivative of order α ∈ (, ) with the lower limit zero (see Definition .). The main tools and techniques are zk based on the estimation on Mittag-Leffler functions Eα (z) := ∞ k= (kα+) and Eα,α (z) := ∞ zk k= (αk+α) for z ≤  (see Lemma  in []) and a singular Gronwall type integral inequality (see Theorem  in []). The classical concept of Ulam-Hyers stability (see []) has applicable significance since it means that if we are dealing with Ulam-Hyers stable system then one does not seek the exact solution. All what is required is to find a function which satisfies a suitable approximation inequality. In other words, Ulam-Hyers stability can guarantee that there exists a close exact solution. This approach is quite useful in many applications such as numerical analysis and optimization, where seeking the exact solution is impossible. In this paper, we are interested in discussing exp-type Ulam-Hyers stability (see Definition .), a special Ulam-Hyers-Rassias stability, of equation () for the case λ > , on a © 2015 Gao et al. This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Gao et al. Advances in Difference Equations (2015) 2015:238

Page 2 of 10

finite time interval J = [, ] and an unbounded interval J = (, ∞), i.e., c

Dαt x(t) = λx(t) + f t, x(t) ,

α ∈ (, ), t ∈ J, λ > .

()

Let > , f : J × R → R be a continuous function. Consider equation () and the following inequality: c α D y(t) – λy(t) – f t, y(t) ≤ , t

α ∈ (, ), t ∈ J, λ > .

()

Definition . Let C(J, R) be the set of all continuous functions from J into R. Equation () is exp-type Ulam-Hyers stable if there exists c >  such that for each >  and for each solution y ∈ C(J, R) of inequality () there exists a solution x ∈ C(J, R) of equation () with  y(t) – x(t) ≤ ceλ α t ,

α ∈ (, ), t ∈ J, λ > .

Concerning the issues in [], the authors widely use certain estimation: Eα (z) ≤ ,  for z ≤ , which are reported in []. However, concerning our current isEα,α (z) ≤ (α) sues, we need more explicit asymptotic behavior on Mittag-Leffler functions Eα (z) and Eα,α (z) for z > . Following Theorem . in [], for z > , Eα,α (z) can be formulated by  –α  Eα,α (z) = z α ez α + α

∞

Kα,α (r, z) dr, 

where Kα,α (r, z) =

r sin(π( – α))  –α –r α r α e .  πα r – rz cos(πα) + z

It is remarkable that Cong et al. [] apply the above formula to obtain explicit asymptotic behavior on Mittag-Leffler functions Eα (z) and Eα,α (z) for z >  (see Lemma .). These explicit asymptotic properties are helpful to solve our problems. In Section , we introduce notations, definitions, and preliminary facts. In Section , we mainly prove exp-type Ulam-Hyers stability results for equation () on a compact interval [, ]. In Section , we mainly prove exp-type Ulam-Hyers stability result for equation () on an unbounded interval (, ∞). Finally, examples are given to illustrate our theoretical results.

2 Preliminaries  Let C(J, R) be the Banach space of all continuous functions from J into R with the λ α norm y

 λα



= supt∈J e–λ α t |y(t)|.

Definition . (see []) The Caputo derivative of order γ for a function f : [, ∞) → R can be written as c

γ γ Dt f (t) = L Dt

f (t) –

n– k t k=

k!

f

(k)

() ,

t > , n –  < γ < n,


Page 3 of 10

γ

where L Dt f denotes the Riemann-Liouville derivative of order γ with the lower limit zero for a function f , which is given by L

γ Dt f (t) =

dn  (n – γ ) dt n

t 

f (s) ds, (t – s)γ +–n

t > , n –  < γ < n.

We need asymptotic behavior of Mittag-Leffler functions Eα (z) and Eα,α (z) for z > . Lemma . (see Lemma (i), (ii) and Lemma (ii); Lemma (ii) via their proofs in []) Let α ∈ (, ) and λ >  be arbitrary. 

(i) For all t > , |Eα (λt α ) – α eλ α t | ≤

m(α,λ) , tα

where

 ∞ 

∞ –r α dr  r α e–r α dr  e , . m(α, λ) = max sin(πα)παλ sin(πα)παλ

()

(ii) For all t > ,  α– m(α, λ) α t Eα,α λt α –  λ –α α eλ t ≤ . α t α+ (iii) For all t ∈ [, ],

t



α Eα,α (λ) Eα (λ) (t – s)α– Eα,α λ(t – s)α – λ –α α E λt α e–λ s ds ≤ + . α α λ

(iv) For all t ∈ (, ∞),

t

α (t – s)α– Eα,α λ(t – s)α – λ –α α E λt α e–λ s ds α



≤ m(α, λ)

 –α  m(α, λ)   + + Eα,α+ (λ) + λ α eλ α + := M(α, λ). α λ α λ

Remark . A function y ∈ C(J, R) is a solution of inequality () if and only if there exists a function g ∈ C(J, R) (which is not dependent on y) such that (i) |g(t)| ≤ , t ∈ J, (ii) c Dαt y(t) = λy(t) + f (t, y(t)) + g(t), t ∈ J. Lemma . Let y ∈ C([, ∞), R) be a solution of inequality (). Then y is a solution of the following integral inequality: t y(t) – Eα λt α y() – (t – s)α– Eα,α λ(t – s)α f s, y(s) ds 

t

≤ 



eλ α t m(α, λ) ds + , α+ (t – s) αλ

t ∈ [, ∞).

Proof Indeed, by Remark ., we have c

Dαt y(t) = λy(t) + f t, y(t) + g(t),

t ∈ [, ∞).


Page 4 of 10

From pp.-, (..)-(..) in [] we have y(t) = Eα λt α y() +

t

(t – s)α– Eα,α λ(t – s)α f s, y(s) ds



t

+

(t – s)α– Eα,α λ(t – s)α g(s) ds,

t ∈ [, ∞).



For t = , y() – Eα ()y() =  < , where we use the fact Eα () = . For all ∞ > t > , we have t y(t) – Eα λt α y() – (t – s)α– Eα,α λ(t – s)α f s, y(s) ds 

t = (t – s)α– Eα,α λ(t – s)α g(s) ds  t  –α λ α (t–s) α– α g(s) ds α ≤ (t – s) Eα,α λ(t – s) – α λ e  t  –α  + λ α eλ α (t–s) g(s) ds α  t t m(α, λ)  –α λ α (t–s) λ α e ds + ds ≤ α+  (t – s)  α



t

eλ α t m(α, λ) – ds + α+ (t – s) αλ αλ

t

eλ α t m(α, λ) , ds + (t – s)α+ αλ

= 

≤





where we use Lemma .(ii). The decried result is obtained.

Remark . Note that for some arbitrary fixed t >  and t > δ > ,

t



 ds = lim δ→ (t – s)α+



t –δ

  ds = lim δ –α – t–α = ∞. α+ δ→ (t – s) α

So it is not possible to obtain some explicit estimation in this case. Next, we divide our time interval [, ∞) into two subintervals [, ] and (, ∞). Remark . Let y ∈ C([, ], R) be a solution of inequality (). Then y is a solution of the following integral inequality: t y(t) – Eα λt α y() – (t – s)α– Eα,α λ(t – s)α f s, y(s) ds 

t

α (t – s)α– Eα,α λ(t – s)α – λ –α α E λt α e–λ s g(s) ds α

≤ 

+

≤



t

α λ –α α E λt α e–λ s g(s) ds α

t  Eα,α (λ) Eα (λ) –α + + λ α Eα λt α e–λ α s ds α λ 


Page 5 of 10

 Eα,α (λ) Eα (λ) Eα (λ) ≤ + +  – e–λ α t α λ λ

Eα,α (λ) Eα (λ) ≤ + , t ∈ [, ], α λ

()

where we use Remark ., Lemma .(iii) and the fact that Eα is an increasing function. Remark . Let y ∈ C((, ∞), R) be a solution of inequality (). Then y is a solution of the following integral inequality: t y(t) – Eα λt α y() – (t – s)α– Eα,α λ(t – s)α f s, y(s) ds 

t

α (t – s)α– Eα,α λ(t – s)α – λ –α α E λt α e–λ s g(s) ds α

≤ 

+

t

α λ –α α E λt α e–λ s g(s) ds α



t

≤ M(α, λ) +

λ

–α α

 Eα λt α e–λ α s ds



≤ M(α, λ) +

 λ α t m(α, λ) , e + λα λt α

t ∈ (, ∞),

()

where we use Remark ., Lemma .(iv) and apply Lemma .(i) to derive the following fact:

t

λ

–α α

  Eα λt α e–λ α s ds = λ– Eα λt α  – e–λ α t



≤ λ– ≤

  λ α t m(α, λ)  – e–λ α t e + α α t

 λ α t m(α, λ) e + . λα λt α

3 Exp-type Ulam-Hyers stability results on J = [0, 1] Next, we introduce the following assumptions: (H ) f : J × R → R is jointly continuous. (H ) There exists L >  such that f (t, x) – f (t, y) ≤ L|x – y| (H ) ω :=  – L( Eα,αα (λ) +

for each t ∈ J and all x, y ∈ R.

–α Eα (λ) ) – Lλ α Eα (λ) > . λ

Theorem . Assume that (H ), (H ), and (H ) are satisfied. Then equation () is exp-type Ulam-Hyers stable on J = [, ]. Proof Let y ∈ C(J, R) be a solution of inequality (). Denote by x the unique solution of the Cauchy problem

c

Dαt x(t) = λx(t) + f (t, x(t)), x() = y(),

t ∈ J,


Page 6 of 10

that is, x(t) = Eα λt α y() +

t

(t – s)α– Eα,α λ(t – s)α f s, x(s) ds.



By using Lemma .(iii) and (), we have y(t) – x(t) t ≤ y(t) – Eα λt α y() – (t – s)α– Eα,α λ(t – s)α f s, x(s) ds 

t ≤ y(t) – Eα λt α y() – (t – s)α– Eα,α λ(t – s)α f s, y(s) ds  t + (t – s)α– Eα,α λ(t – s)α f s, y(s) ds 

t

(t – s)

–

α–



α Eα,α λ(t – s) f s, x(s) ds

Eα,α (λ) Eα (λ) + α λ t  –α + L (t – s)α– Eα,α λ(t – s)α – λ α Eα λt α e–λ α s y(s) – x(s) ds

≤



t

α λ –α α E λt α e–λ s y(s) – x(s) ds α

+L 

Eα,α (λ) Eα (λ) Eα,α (λ) Eα (λ) λ α t ≤ + +L + e x – y α λ α λ α λ –α + Ltλ α Eα λt α x – y α λ

 Eα,α (λ) Eα (λ) Eα,α (λ) Eα (λ) + + Leλ α t + x – y α ≤ λ α λ α λ + Lλ

–α α

Eα (λ)x – y



λα

,

which yields that

≤

x – y

 λα

Eα,α (λ) Eα (λ) Eα,α (λ) Eα (λ) + +L + x – y α λ α λ α λ

+ Lλ

–α α

Eα (λ)x – y



λα

.

Thus,

ω x – y



λα

Eα,α (λ) Eα (λ) ≤ + . α λ

This yields that

 y(t) – x(t) ≤ Eα,α (λ) + Eα (λ) eλ α t , ω α λ The proof is completed.

t ∈ J.


Page 7 of 10

Remark . Under the assumptions of Theorem ., consider

t ∈ J, λ > ,

c

Dαt x(t) = λx(t) + f (t, x(t)), x() = x .

()

Define P : C(J, R) → C(J, R) as follows: (Px)(t) = Eα λt α x +

t

(t – s)α– Eα,α λ(t – s)α f s, x(s) ds.



Then take t ∈ J and x, y ∈ C(J, R), one can obtain Px – Py



λα

≤  – ω x – y



λα

< x – y



λα

.

In view of the contraction mapping principle, that P has a unique fixed point. Thus, problem () has a unique solution.

4 Exp-type Ulam-Hyers stability result on J = (1, ∞) Denote L ([, ∞), R+ ) by all integrable functions from [, ∞) into R+ . Denote C ∞ ([, ∞), R+ ) by all continuous functions from [, ∞) into R+ with y∞ = supt∈[,∞) |y(t)|. Next, we introduce the following assumptions: (H ) There exists L(·) ∈ L ([, ∞), R+ ) ∩ C ∞ ([, ∞), R+ ) such that f (t, x) – f (t, y) ≤ L(t)|x – y|

for each t ∈ J and all x, y ∈ R.

(H ) Set L∞ = supt∈[,∞) |L(t)| and

t

k = sup

L(s) ds =

t∈[,∞) 

∞

L(s) ds < ∞



such that ω =  – L∞ M(α, λ) – k λ

–α  α

α

> .

Theorem . Assume that (H ), (H ), and (H ) are satisfied. Then equation () is exp-type Ulam-Hyers stable on J = (, ∞). Proof To achieve our aim, we apply Lemma .(i), (vi), and () to obtain y(t) – x(t)

≤ M(α, λ) + + L∞ +

 λ α t m(α, λ) e + λα λt α

t

α (t – s)α– Eα,α λ(t – s)α – λ –α α E λt α e–λ s y(s) – x(s) ds α

 t

–α  L(s)λ α Eα λt α e–λ α s y(s) – x(s) ds



≤ M(α, λ) +



Page 8 of 10



+ L∞ e

λα t

M(α, λ)x – y



λα

≤ M(α, λ) +

∞

+

L(s) ds λ 




+ L∞ eλ α t M(α, λ)x – y

 λα

–α α

+ k λ

–α α

 λ α t m(α, λ) + e x – y α λ α tα

 λ α t m(α, λ) x – y α , e + λ α tα

which yields that x – y

 –α  α + L∞ M(α, λ) + k λ x – y α . ≤ λ λα α



λα

This implies that ω x – y



λα

≤

. αλ

Finally, y(t) – x(t) ≤

  eλ α t , αλω

t ∈ J.

Thus, the proof is finished.

Remark . Under the assumptions of Theorem ., problem () also has a unique solution. The uniqueness result can be shown in a similar way as in Remark .. So we omit the details here.

5 Examples Now we give some examples to illustrate our theoretical results. Example . Let α = c



Dt x(t) =

 

and λ =  . We consider the fractional order differential equation

  x(t) + sin x(t),  l

t ∈ [, ], l > ,

()

and the inequality c  Dt y(t) –  y(t) –  sin y(t) ≤ ,  l

t ∈ [, ].

()

Let y ∈ C([, ], R) be a solution of inequality (). Then there exists a function g(t) = 

et– ∈ C([, ], R) such that |g(t)| ≤ , t ∈ [, ], and c Dt y(t) =  y(t) + l sin y(t) + g(t), t ∈ [, ]. Define f (t, x(t)) = l sin x(t) for t ∈ [, ] and set L = l . Then (H ) and (H ) hold. Moreover, we choose l = [(E  ,  (  ) + E  (  )) +  E  (  )], then  

ω :=  –

 l





     E  ,  + E  + E = > ,         

which implies that (H ) holds.


Page 9 of 10

By Theorem ., we have

 y(t) – x(t) ≤  E    + E   e  t , ,   

t ∈ [, ].

Thus, equation () is exp-type Ulam-Hyers stable on [, ] with c = (E  ,  (  ) + E  (  )).  



Example . Let α =  , λ = . We consider the fractional order differential equation c



Dt x(t) = x(t) +

 |x(t)| , let  + |x(t)|

t ∈ (, ∞), l > ,

()

and the inequality c  Dt y(t) – y(t) –  |y(t)| ≤ et , let  + |y(t)|

t ∈ (, ∞).

()

Let y be a solution of (). Then there exists a function g(t) = e–t ∈ C((, ∞), R) such that 

 |y(t)| + g(t), t ∈ (, ∞). let +|y(t)| |x(t)|  Define f (t, x(t)) = let +|x(t)| for t ∈ (, ∞) and set L(t) = let , t ∈ (, ∞). Then (H ) and (H ) ∞ hold. Moreover, L∞ = l and k =  L(s) ds = l < ∞. Now, choose l = (M(  , ) + ), M(  ,) where M(  , ) = √π + E  ,  () + e. Then ω =  – l – l =  , which implies that (H )  

|g(t)| ≤ , t ∈ [, ], and c Dt y(t) = y(t) +

holds. By Theorem ., we have y(t) – x(t) ≤ et ,

t ∈ (, ∞).

Thus, equation () is exp-type Ulam-Hyers stable on (, ∞) with c = . Competing interests The authors declare that they have no competing interests. Authors’ contributions This work was carried out in collaboration between all authors. ZG, XY, and JRW proved the theorems, interpreted the results, and wrote the article. All authors defined the research theme, read and approved the manuscript. Author details 1 College of Applied Mathematics, Shanxi University of Finance and Economics, Taiyuan, Shanxi 030006, P.R. China. 2 Department of Mathematics, Guizhou University, Guiyang, Guizhou 550025, P.R. China. Acknowledgements The authors are grateful to the referees for their careful reading of the manuscript and valuable comments. The authors thank for the help from the editor too. This work is partially supported by the Youth Science Foundation of Shanxi University of Finance and Economic (2014003) and National Natural Science Foundation of China (11201091). Received: 2 March 2015 Accepted: 18 July 2015 References 1. Jung, S-M: Hyers-Ulam-Rassias Stability of Functional Equations in Mathematical Analysis. Hadronic Press, Palm Harbor (2001) 2. Jung, S-M: Hyers-Ulam-Rassias Stability of Functional Equations in Nonlinear Analysis. Springer, New York (2011) ´ 3. Brzde¸k, J, Brillouët-Belluot, N, Cieplinski, K, Xu, B: Ulam’s type stability. Abstr. Appl. Anal. 2012, Article ID 329702 (2012) 4. Rus, IA: Ulam stability of ordinary differential equations. Stud. Univ. Babe¸s-Bolyai, Math. 54, 125-133 (2009) 5. Popa, D, Ra¸sa, I: On the Hyers-Ulam stability of the linear differential equation. J. Math. Anal. Appl. 381, 530-537 (2011) 6. Hegyi, B, Jung, S-M: On the stability of Laplace’s equation. Appl. Math. Lett. 26, 549-552 (2013) 7. Wang, J, Li, X: Eα -Ulam type stability of fractional order ordinary differential equations. J. Appl. Math. Comput. 45, 449-459 (2014)


8. Wang, J, Fe˘ckan, M, Zhou, Y: Presentation of solutions of impulsive fractional Langevin equations and existence results. Eur. Phys. J. Spec. Top. 222, 1855-1872 (2013) 9. Ye, H, Gao, J, Ding, Y: A generalized Gronwall inequality and its application to a fractional differential equation. J. Math. Anal. Appl. 328, 1075-1081 (2007) 10. Hyers, DH: On the stability of the linear functional equation. Proc. Natl. Acad. Sci. USA 27, 222-224 (1941) 11. Gorenflo, R, Loutchko, J, Luchko, Y: Computation of the Mittag-Leffler function Eα,β (z) and its derivative. Fract. Calc. Appl. Anal. 5, 491-518 (2002). Correction: Fract. Calc. Appl. Anal. 6, 111-112 (2003) 12. Cong, ND, Doan, TS, Siegmund, S, Tuan, HT: On stable manifolds for planar fractional differential equations. Appl. Math. Comput. 226, 157-168 (2014) 13. Kilbas, AA, Srivastava, HM, Trujillo, JJ: Theory and Applications of Fractional Differential Equations. Elsevier, Amsterdam (2006)

Page 10 of 10