Next Article in Journal
A New Wavelet Tool to Quantify Non-Periodicity of Non-Stationary Economic Time Series
Next Article in Special Issue
On the Stability of Linear Incommensurate Fractional-Order Difference Systems
Previous Article in Journal
Fekete-Szegö Type Problems and Their Applications for a Subclass of q-Starlike Functions with Respect to Symmetrical Points
Previous Article in Special Issue
Asymptotic Stability of the Solutions of Neutral Linear Fractional System with Nonlinear Perturbation
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence and Stability Analysis for Fractional Impulsive Caputo Difference-Sum Equations with Periodic Boundary Condition

by
Rujira Ouncharoen
1,
Saowaluck Chasreechai
2,* and
Thanin Sitthiwirattham
3,*
1
Research Center in Mathematics and Applied Mathematics, Department of Mathematics, Faculty of Science, Chiang Mai University, Chiang Mai 50200, Thailand
2
Department of Mathematics, Faculty of Applied Science, King Mongkut’s University of Technology North Bangkok, Bangkok 10800, Thailand
3
Mathematics Department, Faculty of Science and Technology, Suan Dusit University, Bangkok 10300, Thailand
*
Authors to whom correspondence should be addressed.
Mathematics 2020, 8(5), 843; https://doi.org/10.3390/math8050843
Submission received: 6 May 2020 / Revised: 20 May 2020 / Accepted: 21 May 2020 / Published: 22 May 2020
(This article belongs to the Special Issue Stability Analysis of Fractional Systems)

Abstract

:
In this paper, by using the Banach contraction principle and the Schauder’s fixed point theorem, we investigate existence results for a fractional impulsive sum-difference equations with periodic boundary conditions. Moreover, we also establish different kinds of Ulam stability for this problem. An example is also constructed to demonstrate the importance of these results.

1. Introduction

In this paper, we study the following periodic boundary value problem for fractional impulsive difference-sum equations:
Δ C α u ( t ) = F t + α 1 , u ( t + α 1 ) , Ψ γ u ( t + α 1 ) , t N 0 , T , t + α 1 t k A u ( α 1 ) + B Δ β u ( α + β 1 ) = C u ( T + α ) + D Δ β u ( T + α + β )
where the impluse conditions subjected to I k , J k : R R are given by
Δ u ( t k ) = I k u ( t k 1 ) , k = 1 , 2 , , p Δ Δ β u ( t k + β ) = J k Δ β u ( t k + β 1 ) , k = 1 , 2 , , p
and α , β ( 0 , 1 ) ; N 0 , T : = { 0 , 1 , , T } ; t k N α , T + α 1 and t 0 = α 1 < t 1 < t 2 < < t p < T + α , t k + 1 t k 2 ; Δ u ( t k ) = u ( t k + 1 ) u ( t k ) ; Δ Δ β u ( t k + β ) = Δ β u ( t k + β + 1 ) Δ β u ( t k + β ) ; A , B , C , D R ; F : N α 1 , T + α × R × R R is a continuous function; and for φ : N α 1 , T + α × N α 1 , T + α [ 0 , ) ,
( Ψ γ u ) ( t ) : = [ Δ γ φ u ] ( t + γ ) = 1 Γ ( γ ) s = α γ 1 t γ ( t σ ( s ) ) γ 1 ̲ φ ( t , s + γ ) u ( s + γ ) .
Fractional calculus has been gaining more attention over the past decade as this calculus has been addressed to various problems used in science and engineering; see [1,2,3,4,5,6,7,8]. For fractional difference calculus theory, which is the discrete case of fractional calculus, there is still not much interest among mathematical researchers. Basic knowledge of fractional difference calculus can be found in [9]. Some interesting results on fractional difference calculus have been studied, which has helped to develop the basic theory of this calculus; see [10,11,12,13,14,15,16,17,18,19,20,21,22,23,24,25,26,27,28,29,30,31,32,33,34,35,36,37,38,39,40] and references cited therein. The extension of applications of fractional difference calculus; see [41,42,43] and references cited therein.
Currently, the studies of boundary value problems for fractional difference equations are shown both extensive and more complex of conditions. There are some recent papers to study the existence and stability results of fractional difference equation [44,45,46,47,48,49,50,51]. However, a few paper has been made to develop the theory of the existence and stability results of fractional difference equations with impulse [52,53].
These results are the motivation for this research. By using the Banach contraction principle and the Schauder’s fixed point theorem, we aim to prove the existence results for the problem (1) and (2). Moreover, we also provide a condition for the different kinds of Ulam stability for the problem (1) and (2). Then, we present an illustrative example.

2. Preliminaries

In this section, we recall some notations, definitions and lemmas which will be used in the main results.
Definition 1.
The generalized falling function is defined as
t α ̲ : = Γ ( t + 1 ) Γ ( t + 1 α ) ,
when t + 1 α is not a pole of the Gamma function. If t + 1 is not a pole and when t + 1 α is a pole, then t α ̲ = 0 .
Definition 2.
For α > 0 and f defined on N a : = { a , a + 1 , } , the fractional sum of order α of f is defined as
Δ α f ( t ) : = 1 Γ ( α ) s = a t α ( t σ ( s ) ) α 1 ̲ f ( s ) ,
where t N a + α and σ ( s ) = s + 1 .
Definition 3.
For α > 0 , N N is satisfied with α ( N 1 , N ] and f defined on N a , the Riemann-Liouville fractional delta difference of order α of f is defined as
Δ α f ( t ) : = Δ N Δ ( N α ) f ( t ) = 1 Γ ( α ) s = a t + α ( t σ ( s ) ) α 1 ̲ f ( s ) ,
where t N a + N α . The Caputo fractional difference of order α of f is defined as
Δ C α f ( t ) : = Δ ( N α ) Δ N f ( t ) = 1 Γ ( N α ) s = a t ( N α ) ( t σ ( s ) ) N α 1 ̲ Δ N f ( s ) ,
where t N a + N α . If α = N , then Δ α f ( t ) = Δ C α f ( t ) = Δ N f ( t ) .
Lemma 1
([10]). Let α > 0 and α ( N 1 , N ] . Then
Δ α Δ C α y ( t ) = y ( t ) + C 0 + C 1 t 1 ̲ + C 2 t 2 ̲ + + C N 1 t N 1 ̲ ,
for some C i R , 0 i N 1 .
Now, we aim to study the following linear variant of the boundary value problem (1) and (2).
Lemma 2.
Let Λ 0 ; α , β ( 0 , 1 ) ; A , B , C , D R ; t 0 = α 1 < t 1 < t 2 < < t p < T + α ; for k = 1 , 2 , , p , t k + 1 t k 2 and t k N α , T + α 1 ; I k , J k : R R and h : N α 1 , T + α R be continuous. Then the following problem
Δ C α u ( t ) = h ( t + α 1 ) , t N 0 , T , t + α 1 t k A u ( α 1 ) + B Δ β u ( α + β 1 ) = C u ( T + α ) + D Δ β u ( T + α + β )
where the impluse conditions subjected to I k , J k are given by
Δ u ( t k ) = I k u ( t k 1 ) , k = 1 , 2 , , p Δ Δ β u ( t k + β ) = J k Δ β u ( t k + β 1 ) , k = 1 , 2 , , p
has the unique solution which is in a form
u ( t ) = 1 Λ [ C i = 1 p I i u ( t i 1 ) + D i = 1 p J i Δ β u ( t i + β 1 ) + C Φ 1 [ h ] + D Φ 2 [ h ] ] + 1 Γ ( α ) s = t 0 t 1 ( t s + α 2 ) α 1 ̲ h ( s ) , t N t 0 , t 1 1 Λ [ C i = 1 p I i u ( t i 1 ) + D i = 1 p J i Δ β u ( t i + β 1 ) + C Φ 1 [ h ] + D Φ 2 [ h ] ] + i = 1 k I i u t i 1 + 1 Γ ( α ) i = 1 k s = t i 1 t i 1 t i s + α 2 α 1 ̲ h ( s ) + 1 Γ ( α ) s = t k t 1 t s + α 2 α 1 ̲ h ( s ) , t N t k + 1 , t k + 1
where Δ u ( t k ) = u ( t k + 1 ) u ( t k ) ; Δ Δ β u ( t k + β ) = Δ β u ( t k + β + 1 ) Δ β u ( t k + β ) ; the functionals Φ 1 [ h ] , Φ 2 [ h ] and the constant Λ are defined by
Φ 1 [ h ] : = 1 Γ ( α ) i = 1 p s = t i 1 t i 1 t i s + α 2 α 1 ̲ h ( s ) + s = t p T + α 1 T + 2 α s 2 α 1 ̲ h ( s )
Φ 2 [ h ] : = 1 Γ ( α ) Γ ( β ) [ i = 1 p r = t i 1 + 1 t i s = t i 1 r 1 t i + β s 1 β 1 ̲ r s + α 2 α 1 ̲ h ( s ) + s = t p + 1 T + α s = t k r 1 T + α + β s 1 β 1 ̲ r s + α 2 α 1 ̲ h ( s ) ]
Λ : = A + B C D Γ ( β ) i = 1 p s = t i 1 t i t i s + β 1 β 1 ̲ + s = t p T + α T + α + β s 1 β 1 ̲ .
Proof. 
For t N t 0 , t 1 , by Lemma 1 and taking the fractional sum of order α for (3) and (4), a general solution can be written as
u ( t ) = C 0 + 1 Γ ( α ) s = 0 t α ( t σ ( s ) ) α 1 ̲ h ( s + α 1 ) = C 0 + 1 Γ ( α ) s = t 0 t 1 ( t s + α 2 ) α 1 ̲ h ( s )
By substituting t = t 1 into (9), we have
u ( t 1 ) = C 0 + 1 Γ ( α ) s = t 0 t 1 1 ( t 1 s + α 2 ) α 1 ̲ h ( s ) .
If t N t 1 + 1 , t 2 , then we have
u ( t ) = u ( t 1 + 1 ) + 1 Γ ( α ) s = t 1 t 1 ( t s + α 2 ) α 1 ̲ h ( s ) = = Δ u ( t 1 ) + u ( t 1 ) + 1 Γ ( α ) s = t 1 t 1 ( t s + α 2 ) α 1 ̲ h ( s ) = I 1 u t 1 1 + C 0 + 1 Γ ( α ) s = t 0 t 1 1 ( t 1 s + α 2 ) α 1 ̲ h ( s ) + 1 Γ ( α ) s = t 1 t 1 ( t s + α 2 ) α 1 ̲ h ( s ) .
If t N t 2 + 1 , t 3 , then we have
u ( t ) = u ( t 2 + 1 ) + 1 Γ ( α ) s = t 2 t 1 ( t s + α 2 ) α 1 ̲ h ( s ) = Δ u ( t 2 ) + u ( t 2 ) + 1 Γ ( α ) s = t 2 t 1 ( t s + α 2 ) α 1 ̲ h ( s ) = I 2 u t 2 1 + [ I 1 u t 1 1 + C 0 + 1 Γ ( α ) s = t 0 t 1 1 ( t 1 s + α 2 ) α 1 ̲ h ( s ) + 1 Γ ( α ) s = t 1 t 2 1 ( t 2 s + α 2 ) α 1 ̲ h ( s ) ] + 1 Γ ( α ) s = t 2 t 1 ( t s + α 2 ) α 1 ̲ h ( s ) .
Repeating the process, the solution u ( t ) for t N t k + 1 , t k + 1 ( k = 1 , 2 , , p ) can be written as
u ( t ) = C 0 + i = 1 k I i u t i 1 + 1 Γ ( α ) i = 1 k s = t i 1 t i 1 t i s + α 2 α 1 ̲ h ( s ) + 1 Γ ( α ) s = t k t 1 t s + α 2 α 1 ̲ h ( s ) .
Next, for t N t 0 + β , t 1 + β , taking the fractional sum of order β for (9), we have
Δ β u ( t ) = C 0 Γ ( β ) s = α 1 t β ( t σ ( s ) ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = α 1 t β s = 0 r α ( t σ ( r ) ) β 1 ̲ ( r σ ( s ) ) α 1 ̲ h ( s ) .
Thus, for t N t 0 , t 1 , we have
Δ β u ( t + β ) = C 0 Γ ( β ) s = t 0 t ( t + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = t 0 + 1 t s = t 0 r 1 ( t + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) .
By substituting t = t 1 into (15), we have
Δ β u ( t 1 + β ) = C 0 Γ ( β ) s = t 0 t 1 ( t 1 + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = t 0 + 1 t 1 s = t 0 r 1 ( t 1 + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) .
If t N t 1 + 1 , t 2 , then we have
Δ β u ( t + β ) = Δ β u ( t 1 + β + 1 ) + C 0 Γ ( β ) s = t 1 t ( t + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = t 1 + 1 t s = t 1 r 1 ( t + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) = Δ Δ β u ( t 1 + β ) + Δ β u ( t 1 + β ) + C 0 Γ ( β ) s = t 1 t ( t + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = t 1 + 1 t s = t 1 r 1 ( t + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) = J 1 Δ β u ( t 1 + β 1 ) + [ C 0 Γ ( β ) s = t 0 t 1 ( t 1 + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = t 0 + 1 t 1 s = t 0 r 1 ( t 1 + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) ] + C 0 Γ ( β ) s = t 1 t ( t + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = t 1 + 1 t s = t 1 r 1 ( t + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) .
If t N t 2 + 1 , t 3 , then we have
Δ β u ( t + β ) = Δ β u ( t 2 + β + 1 ) + C 0 Γ ( β ) s = t 2 t ( t + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = t 2 + 1 t s = t 2 r 1 ( t + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) = Δ Δ β u ( t 2 + β ) + Δ β u ( t 2 + β ) + C 0 Γ ( β ) s = t 2 t ( t + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = t 2 + 1 t s = t 2 r 1 ( t + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) = J 1 Δ β u ( t 1 + β 1 ) + J 2 Δ β u ( t 2 + β 1 ) + C 0 Γ ( β ) s = t 0 t 1 ( t 1 + β s 1 ) β 1 ̲ + s = t 1 t 2 ( t 2 + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) [ s = t 0 + 1 t 1 s = t 0 r 1 ( t 1 + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) + s = t 1 + 1 t 2 s = t 1 r 1 ( t 2 + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) ] + C 0 Γ ( β ) s = t 2 t ( t + β s 1 ) β 1 ̲ + 1 Γ ( α ) Γ ( β ) s = t 2 + 1 t s = t 2 r 1 ( t + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) .
Repeating the process, the solution u ( t ) for t N t k + 1 , t k + 1 ( k = 1 , 2 , , p ) can be written as
Δ β u ( t ) = i = 1 k J i Δ β u t i + β 1 + C 0 Γ ( β ) i = 1 k s = t i 1 t i 1 t i s + α 2 β 1 ̲ h ( s ) + s = t k t t s + β 1 β 1 ̲ + 1 Γ ( α ) Γ ( β ) [ i = 1 k s = t i 1 + 1 t i s = t i 1 r 1 ( t i + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) + s = t k + 1 t s = t k r 1 ( t + β r 1 ) β 1 ̲ ( r s + α 2 ) α 1 ̲ h ( s ) ] .
By substituting t = α 1 into (9), (16); t = T + α into (13), (19); and using the condition A u ( α 1 ) + B Δ β u ( α + β 1 ) = C u ( T + α ) + D Δ β u ( T + α + β ) , we have
C 0 = 1 Λ C i = 1 p I i u ( t i 1 ) + D i = 1 p J i Δ β u ( t i + β 1 ) + C Φ 1 [ h ] + D Φ 2 [ h ]
where Φ 1 [ h ] , Φ 2 [ h ] and Λ are defined in (6)–(8), respectively. Substituting C 0 into (9) and (13), we have (5). □

3. Main Results

3.1. Existence and Uniqueness Solution

Let C = C ( N α 1 , T + α , R ) be the Banach space equipped with the norm u = max t N α 1 , T + α | u ( t ) | . Defined the operator T : C C by
( T u ) ( t ) : = 1 Λ C i = 1 p I i u ( t i 1 ) + D i = 1 p J i Δ β u ( t i + β 1 ) + C Φ 1 * [ F u ] + D Φ 2 * [ F u ] + 1 Γ ( α ) s = t 0 t 1 ( t s + α 2 ) α 1 ̲ F s , u ( s ) , Ψ γ u ( s ) , t N t 0 , t 1 1 Λ [ C i = 1 p I i u ( t i 1 ) + D i = 1 p J i Δ β u ( t i + β 1 ) + C Φ 1 * [ F u ] + D Φ 2 * [ F u ] ] + i = 1 k I i u t i 1 + 1 Γ ( α ) i = 1 k s = t i 1 t i 1 t i s + α 2 α 1 ̲ F s , u ( s ) , Ψ γ u ( s ) + 1 Γ ( α ) s = t k t 1 t s + α 2 α 1 ̲ F s , u ( s ) , Ψ γ u ( s ) , t N t k + 1 , t k + 1
where the functionals Φ 1 * [ F u ] and Φ 2 * [ F u ] are defined as
Φ 1 * [ F u ] = 1 Γ ( α ) [ i = 1 p s = t i 1 t i 1 t i s + α 2 α 1 ̲ F s , u ( s ) , Ψ γ u ( s ) + s = t p T + α 1 T + 2 α s 2 α 1 ̲ F s , u ( s ) , Ψ γ u ( s ) ]
Φ 2 * [ F u ] = 1 Γ ( α ) Γ ( β ) [ i = 1 p r = t i 1 + 1 t i s = t i 1 r 1 t i + β s 1 β 1 ̲ r s + α 2 α 1 ̲ × F s , u ( s ) , Ψ γ u ( s ) + s = t p + 1 T + α s = t k r 1 T + α + β s 1 β 1 ̲ r s + α 2 α 1 ̲ × F s , u ( s ) , Ψ γ u ( s ) ]
and constant Λ is given in (8). Observe that the operator in (21) has the fixed points which are the solutions of the problem (1) and (2).
Theorem 1.
Let F : N α 1 , T + α × R × R R , I k , J k : R R , k = 1 , 2 , , p be continuous; φ : N α 1 , T + α × N α 1 , T + α [ 0 , ) be continuous with φ 0 = max { φ ( t 1 , s ) : ( t , s ) N α 1 , T + α × N α 1 , T + α } . In addition, suppose that
( H 1 )
There exist constants 1 , 2 > 0 such that
| F [ t , u ( t ) , ( Ψ γ u ) ( t ) ] F [ t , v ( t ) , ( Ψ γ v ) ( t ) ] | 1 | u v | + 2 | ( Ψ γ u ) ( Ψ γ v ) |
for each t N α 1 , α + T and u , v R .
( H 2 )
There exist constants λ 1 , λ 2 > 0 such that
| I k ( u ) I k ( v ) | λ 1 | u v | | J k Δ β u J k Δ β v | λ 2 | Δ β u Δ β v |
for each u , v C and k = 1 , 2 , , p .
( H 3 )
L P + λ 1 Q + λ 2 R < 1 ,
where
L = 1 + 2 φ 0 ( T + α + γ ) γ ̲ Γ ( γ + 1 )
P = C Λ + 1 Ω 1 + D Λ Ω 2
Q = p ( p + 1 ) 2 C Λ + 1
R = D Λ ( T + α + β ) β ̲ Γ ( β + 1 )
Ω 1 = p ( p + 1 ) 2 + 1 ( T + α ) α ̲ Γ ( α + 1 )
Ω 2 = p ( p + 1 ) 2 + 1 ( T + α ) α ̲ ( T + α + β ) β ̲ Γ ( α + 1 ) Γ ( β + 1 ) .
Therefore, the problem (1) and (2) has a unique solution on N α 1 , α + T .
Proof. 
We will show that T is a contraction. Let
H | u v | ( t ) : = | F [ t , u ( t ) , ( Ψ γ u ) ( t ) ] F [ t , v ( t ) , ( Ψ γ v ) ( t ) ] | .
For any u , v C , we have
| Φ 1 * [ F u ] Φ 1 * [ F v ] | = 1 Γ ( α ) | i = 1 p s = t i 1 t i 1 t i s + α 2 α 1 ̲ H | u v | ( s ) + s = t p T + α 1 T + 2 α s 2 α 1 ̲ H | u v | ( s ) | 1 | u v | + 2 | ( Ψ γ u ) ( Ψ γ v ) | Γ ( α ) | i = 1 p s = t i 1 t i 1 t i s + α 2 α 1 ̲ + s = t p T + α 1 T + 2 α s 2 α 1 ̲ | 1 + 2 φ 0 ( T + α + γ ) γ ̲ Γ ( γ + 1 ) | u v | p ( p + 1 ) 2 + 1 ( T + α ) α ̲ Γ ( α + 1 ) = L Ω 1 | u v | .
Similary, we get
Φ 2 * [ F u ] Φ 2 * [ F v ] L Ω 2 | u v | .
Next, for each t N t k + 1 , t k + 1 , k = 1 , 2 , , p , we obtain
| ( T u ) ( t ) ( T v ) ( t ) | 1 Λ | C i = 1 p | I i ( u t i 1 ) I k ( v t i 1 ) | + D i = 1 p | J i Δ β u ( t i + β 1 ) J i Δ β v ( t i + β 1 ) | + C Φ 1 * [ F u ] Φ 1 * [ F v ] + D Φ 2 * [ F u ] Φ 2 * [ F v ] | + i = 1 k | I i ( u t i 1 ) I k ( v t i 1 ) | + 1 Γ ( α ) i = 1 k s = t i 1 t i 1 t i s + α 2 α 1 ̲ H | u v | ( s ) + 1 Γ ( α ) s = t k t 1 t s + α 2 α 1 ̲ H | u v | ( s ) 1 | Λ | [ | C | p ( p + 1 ) 2 λ 1 | u v | + | D | p ( p + 1 ) 2 λ 2 Δ β u Δ β v + | C | L Ω 1 | u v | + | D | L Ω 2 u v ] + p ( p + 1 ) 2 λ 1 | u v | + p ( p + 1 ) 2 + 1 ( T + α ) α ̲ Γ ( α + 1 ) L | u v | { 1 | Λ | | C | p ( p + 1 ) 2 λ 1 + | C | L Ω 1 + | D | L Ω 2 + p ( p + 1 ) 2 λ 1 + p ( p + 1 ) 2 + 1 ( T + α ) α ̲ Γ ( α + 1 ) L } | u v | + D Λ p ( p + 1 ) 2 λ 2 ( T + α + β ) β ̲ Γ ( β + 1 ) u v L P + λ 1 Q + λ 2 R | u v | .
Obviously, for each t N t 0 , t 1 , we obtain | ( T u ) ( t ) ( T v ) ( t ) | < | u v | .
Thus, for each t N α 1 , T + α , we have ( T u ) ( t ) ( T v ) ( t ) u v .
So, T is a contraction. By the Banach contraction principle, we can conclude that T has a fixed point. Hence, the problem (1) and (2) has a unique solution on t N α 1 , α + T . □

3.2. Existence of at Least One Solution

In the next result, we use of the Schauder’s fixed point theorem to discuss the existence of at least one solution of (1) and (2).
Theorem 2.
Suppose that ( H 1 ) and ( H 3 ) hold. Then, there exists at least one solution for the problem (1) and (2).
Proof. 
We organize the proof as follows:
Step I. Examine that T map bounded sets into bounded sets in B R = { u C : u R } . Let max | F ( t , 0 , 0 , 0 ) | = K , max | I k ( u ) | = M , max | J k Δ β u | = M , and choose a constant
R P K 1 L P + M Q + N R .
Let
| S ( t , u , 0 ) | = | F t , u , Δ β u F t , 0 , 0 | + | F t , 0 , 0 | .
For any u , v C , we obtain
| Φ 1 * [ F u ] | = 1 Γ ( α ) | i = 1 p s = t i 1 t i 1 t i s + α 2 α 1 ̲ | S ( s , u , 0 ) | + s = t p T + α 1 T + 2 α s 2 α 1 ̲ | S ( s , u , 0 ) | | 1 + 2 φ 0 ( T + α + γ ) γ ̲ Γ ( γ + 1 ) | u | + K Ω 1 = L Ω 1 R + K Ω 1 .
Similarly, we obtain
Φ 2 * [ F u ] L Ω 2 R + K Ω 2 .
So, for each t N t k + 1 , t k + 1 , k = 1 , 2 , , p , we get
| ( T u ) ( t ) | L P + M Q + N R R + K P R .
Obviously, for each t N t 0 , t 1 , we have | T u ( t ) | < R .
Therefore, T u R for each t N α 1 , T + α , which implies that F is uniformly bounded.
Step II. It is obvious that the operator T is continuous on B R since the continuity of F.
Step III. Examine that T is equicontinuous on B R . For any ϵ > 0 , there exists δ > 0 such that for τ 1 , τ 2 N α 1 , T + α with τ 1 < τ 2
| τ 2 α ̲ τ 1 α ̲ | < ϵ Γ ( α + 1 ) F   whenever   | τ 2 τ 1 | < δ .
Then, we obtain
| ( T u ) ( τ 2 ) ( T u ) ( τ 1 ) | F ( T ω 0 ) α Γ q ( α + 1 ) Γ q ( β + 1 ) | 1 Γ ( α ) s = t k τ 2 1 τ 2 s + α 2 α 1 ̲ 1 Γ ( α ) s = t k τ 1 1 τ 1 s + α 2 α 1 ̲ | = F Γ ( α + 1 ) | τ 2 α ̲ τ 1 α ̲ | < ϵ .
This implies that the set T ( B R ) is an equicontinuous set.
Hence, by the Arzelá–Ascoli theorem, we deduce that T : C C is completely continuous. Therefore, it follows from the Schauder’s fixed point theorem that problem (1)-(2) has at least one solution. □

4. Ulam Stability Analysis Results

Based on the concept in Wang et al. [54], we provide Ulam’s type stability concepts for the problem (1) and (2). Consider the inequalities:
Δ C α z ( t ) F t + α 1 , z ( t + α 1 ) , Ψ γ z ( t + α 1 ) ϵ Δ u ( t k ) I k u ( t k 1 ) ϵ Δ Δ β u ( t k + β ) J k Δ β u ( t k + β 1 ) ϵ ,
Δ C α z ( t ) F t + α 1 , z ( t + α 1 ) , Ψ γ z ( t + α 1 ) ϵ ρ ( t ) Δ u ( t k ) I k u ( t k 1 ) ϵ ψ 1 Δ Δ β u ( t k + β ) J k Δ β u ( t k + β 1 ) ϵ ψ 2 ,
for t N 0 , T , t + α 1 t k , k = 1 , 2 , , p .
Definition 4.
The problem (1) and (2) is the Ulam–Hyers stable if there exists a real number c F , p > 0 such that for each ϵ > 0 and for each solution z C of the inequality (38), there exists a solution u C of problem (1) and (2) with
| z ( t ) u ( t ) | c F , p ϵ .
The problem (1) and (2) is the generalized Ulam–Hyers stable, if we substitute the function θ F , p ( ϵ ) C ( R + , R + ) , θ F , p ( 0 ) = 0 for the constant c F , p ϵ on inequality (40).
Definition 5.
The problem (1) and (2) is the Ulam–Hyers–Rassias stable with respect to ( ρ , ψ 1 , ψ 2 ) if there exists a real number c F , p , ρ > 0 such that for each ϵ , ψ 1 , ψ 2 > 0 and for each solution z C of the inequality (40), there exists a solution u C of problem (1) and (2) with
| z ( t ) u ( t ) | c F , p , ρ ϵ ρ ( t ) + ψ 1 + ψ 2 .
The problem (1) and (2) is the generalized Ulam–Hyers–Rassias stable, if we substitute the function ρ * ( t ) for the function ϵ ρ ( t ) and the constants ψ 1 * , ψ 2 * for on the constants ϵ ψ 1 , ϵ ψ 2 on inequalities (39) and (41).
Remark 1.
A function z C is a solution of the inequality (38) if and only if there exist ϕ C (depend on z) and sequence ω k , κ k , k = 1 , 2 , , p with ω = max { ω k } , κ = max { κ k } , such that
( i )
| ϕ ( t ) | ϵ for t N t k + 1 , t k + 1 , and | ω | , | κ | ϵ ;
( i i )
Δ C α z ( t ) = F t + α 1 , z ( t + α 1 ) , Ψ γ z ( t + α 1 ) + ϕ ( t + α 1 )
for t N 0 , T , t + α 1 t k ;
( i i i )
Δ u ( t k ) = I k u ( t k 1 ) + ω k ;
( i v )
Δ Δ β u ( t k + β ) = J k Δ β u ( t k + β 1 ) + κ k .
Remark 2.
A function z C is a solution of the inequality (39) if and only if there exist ϕ C (depend on z) and sequence ω k , κ k , k = 1 , 2 , , p with ω = max { ω k } , κ = max { κ k } ,
( i )
| ϕ ( t ) | ϵ ρ ( t ) for t N t k + 1 , t k + 1 , | ω | ϵ ψ 1 and | κ | ϵ ψ 2 ;
( i i )
Δ C α z ( t ) = F t + α 1 , z ( t + α 1 ) , Ψ γ z ( t + α 1 ) + ϕ ( t + α 1 )
for t N 0 , T , t + α 1 t k ;
( i i i )
Δ u ( t k ) = I k u ( t k 1 ) + ω k ;
( i v )
Δ Δ β u ( t k + β ) = J k Δ β u ( t k + β 1 ) + κ k .
Lemma 3.
If z C is a solution of the inequality (38), then for ϵ > 0 , z is solution of the inequality
| z ( t ) T z ( t ) | ϵ P + Q + R ,
where P , Q and R are defined in (25)–(27).
Proof. 
From Remark 1 and Lemma 2, we have
z ( t ) = 1 Λ [ C i = 1 p I i u ( t i 1 ) + ω i + D i = 1 p J i Δ β u ( t i + β 1 ) + κ i + C Φ 1 * [ F u + ϕ ] + D Φ 2 * [ F u + ϕ ] ] + i = 1 k I i u t i 1 + ω i + 1 Γ ( α ) i = 1 k s = t i 1 t i 1 t i s + α 2 α 1 ̲ F s , u ( s ) , Ψ γ u ( s ) + ϕ ( s ) + 1 Γ ( α ) s = t k t 1 t s + α 2 α 1 ̲ F s , u ( s ) , Ψ γ u ( s ) + ϕ ( s ) .
Hence, we obtain
| z ( t ) ( T z ) ( t ) | 1 | Λ | | C | i = 1 p ω i + | D | i = 1 p κ i + | C | Φ 1 * [ | ϕ | ] + D Φ 2 * [ | ϕ | ] + i = 1 k ω i + 1 Γ ( α ) i = 1 k s = t i 1 t i 1 t i s + α 2 α 1 ̲ ϕ ( s ) + 1 Γ ( α ) s = t k t 1 t s + α 2 α 1 ̲ ϕ ( s ) p ( p + 1 ) 2 C Λ + 1 | ω | + D Λ | κ | + ϕ | Λ | | C | Ω 1 + | D | Ω 2 + ϕ Γ ( α + 1 ) p ( p + 1 ) 2 ( T + α ) α ̲ + t α ̲ ϕ P + | ω | Q + | κ | R ϵ P + Q + R
where P , Q , R , Ω 1 , Ω 2 are defined in (25)–(29). This completes the proof. □
Theorem 3.
Suppose that (H1)–(H3) hold with 1 L P + λ 1 Q + λ 2 R 0 . Then, problem (1) and (2) is Ulam–Hyers stable.
Proof. 
Suppose z C is a solution of (38) and assume that x ( t ) is the unique solution of problem (1) and (2). Consider
| z ( t ) x ( t ) | = | z ( t ) ( T z ) ( t ) + ( T z ) ( t ) x ( t ) | | z ( t ) ( T z ) ( t ) | + | ( T z ) ( t ) x ( t ) | ϵ P + Q + R + L P + λ 1 Q + λ 2 R | z ( t ) x ( t ) | ,
from which we obtain
| z ( t ) x ( t ) | ϵ P + Q + R | 1 L P + λ 1 Q + λ 2 R |
By setting c F , p = P + Q + R | 1 L P + λ 1 Q + λ 2 R | , we have
| z ( t ) x ( t ) | c F , p ϵ .
Hence, the problem (1) and (2) is Ulam–Hyers stable. Next, by setting θ F , p ( ϵ ) = c F , p ϵ with θ F , p ( 0 ) = 0 , the problem (1) and (2) is generalized Ulam–Hyers stable. □
Lemma 4.
If z C is a solution of the inequality (39) and assume that
( H 4 )
p ( p + 1 ) 2 C Λ + 1 | ω | + D Λ | κ | + ϕ | Λ | | C | Ω 1 + | D | Ω 2 + ϕ Γ ( α + 1 ) p ( p + 1 ) 2 ( T + α ) α ̲ + t α ̲
ϵ O ρ ( t ) + ψ 1 + ψ 2 ,
then for ϵ > 0 , z is solution of the following inequality
| z ( t ) T z ( t ) | ϵ O ρ ( t ) + ψ 1 + ψ 2 ,
where O = max { P , Q , R } and P , Q , R , Ω 1 , Ω 2 are defined in (25)–(29).
Proof. 
By the same argument of the proof Lemma 3, we have
| z ( t ) ( T z ) ( t ) | p ( p + 1 ) 2 C Λ + 1 | ω | + D Λ | κ | + ϕ | Λ | | C | Ω 1 + | D | Ω 2 + ϕ Γ ( α + 1 ) p ( p + 1 ) 2 ( T + α ) α ̲ + t α ̲
ϵ ρ ( t ) P + ψ 1 Q + ψ 2 R ϵ O ρ ( t ) + ψ 1 + ψ 2 .
This completes the proof. □
Theorem 4.
Suppose that ( H 1 ) ( H 3 ) hold with 1 L P + λ 1 Q + λ 2 R 0 . Then, problem (1) and (2) is Ulam–Hyers–Rassias stable.
Proof. 
Suppose z C is a solution of (39) and assume that x ( t ) is the unique solution of problem (1) and (2). Consider
| z ( t ) x ( t ) | = | z ( t ) ( T z ) ( t ) + ( T z ) ( t ) x ( t ) | | z ( t ) ( T z ) ( t ) | + | ( T z ) ( t ) x ( t ) | ϵ O ρ ( t ) + ψ 1 + ψ 2 + L P + λ 1 Q + λ 2 R | z ( t ) x ( t ) | ,
from which we obtain
| z ( t ) x ( t ) | ϵ O ρ ( t ) + ψ 1 + ψ 2 | 1 L P + λ 1 Q + λ 2 R |
By setting c F , p = O | 1 L P + λ 1 Q + λ 2 R | , we have
| z ( t ) x ( t ) | c F , p ϵ ρ ( t ) + ψ 1 + ψ 2 .
Hence, the problem (1) and (2) is Ulam–Hyers–Rassias stable. Next, by setting ρ * ( t ) + ψ 1 * + ψ 2 * = ϵ ρ ( t ) + ψ 1 + ψ 2 , the problem (1) and (2) is generalized Ulam–Hyers–Rassias stable. □

5. An Example

To show the application of our theorems, we provide the fractional impulsive difference-sum equations with periodic boundary conditions of the form
Δ C 1 2 u ( t ) = cos 2 t 1 2 π 2 t 199 2 2 u 2 + | u | | u | + 1 + e sin 2 t 1 2 π Ψ 1 3 u t 1 2 , t N 0 , T , t 1 2 t k 10 u 1 2 + 20 Δ 2 3 u 1 6 = 5 u 21 2 + 15 Δ 2 3 u 67 6
where the impluse conditions subjected to I k , J k : R R are given by
Δ u ( t k ) = 1 k + 100 cos u ( t k 1 ) , k = 1 , 2 , 3 Δ Δ 2 3 u t k + 2 3 = 1 ( k + 50 ) 2 tan 1 Δ 2 3 u t k 1 3 , k = 1 , 2 , 3
and φ ( t , s ) = e | s t | ( t + 10 ) 3 and t k = 1 2 + 2 k t 1 = 5 2 , t 2 = 9 2 , t 3 = 13 2 .
Here α = 1 2 , β = 2 3 , γ = 1 3 , T = 10 , p = 3 , A = 10 , B = 20 , C = 5 , D = 15 , and F t , u ( t ) , Ψ γ u ( t ) = cos 2 π t 2 t + 100 2 u 2 + | u | | u | + 1 + e sin 2 π t Ψ 1 3 u t .
We can find that
Θ = 6.74071 , | Λ | = 75.57065 , Ω 1 = 25.90098 , Ω 2 = 144.78000 ,
P = 46.7733 , Q = 7.46312 and R = 0.73969 .
Let t N 1 2 , 21 2 and u , v R , we have
F t , u , Ψ γ u F t , v , Ψ γ v 1 2 1 2 + 100 2 | u v | + 1 2 e π 1 2 + 100 2 Ψ γ u Ψ γ v .
So, ( H 1 ) holds with 1 = 0.000495 , 2 = 5.6095 × 10 6 , and we have φ 0 = 0.00907 .
For all u , v C and k = 1 , 2 , 3
I k ( u ) I k ( v ) 1 101 | u v | and J k Δ β u I k Δ β v 1 51 2 Δ β u Δ β v
So, ( H 2 ) holds with λ 1 = 1 101 , λ 2 = 1 51 2 .
Finally, ( H 3 ) holds with
L P + λ 1 Q + λ 2 R 0.07418 < 1 .
Hence, by Theorem 1, the boundary value problem (52) and (53) has a unique solution.
In view of Theorem 3, we have
c F , p = P + Q + R | 1 L P + λ 1 Q + λ 2 R | 2128.38328 .
Therefore, problem (52) and (53) is Ulam–Hyers stable and hence generalized Ulam–Hyers stable.
By setting
c F , p , ρ = O | 1 L P + λ 1 Q + λ 2 R | 44.03976 , and
ρ ( t ) = 0.02412 t α ̲ + 1.12568 , ψ 1 = 0.13677 , ψ 2 = 0.01698 .
Thus, by Theorem 4, problem (52) and (53) is Ulam–Hyers–Rassias stable and also generalized Ulam–Hyers–Rassias stable.

Author Contributions

Conceptualization, R.O. and S.C.; Formal analysis, R.O. and S.C.; Funding acquisition, S.C.; Investigation, R.O.; Methodology, R.O., S.C. and T.S.; Writing—original draft, R.O., S.C. and T.S.; Writing—review and editing, R.O., S.C. and T.S. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by King Mongkut’s University of Technology North Bangkok. Contract no.KMUTNB-61-GOV-D-65.

Acknowledgments

This research was supported by Chiang Mai University.

Conflicts of Interest

The authors declare no conflicts of interest regarding the publication of this paper.

References

  1. Zhang, Y.; Liu, X.; Belić, M.R.; Zhong, W.; Zhang, Y.; Xiao, M. Propagation Dynamics of a Light Beam in a Fractional Schrödinger Equation. Phys. Rev. Lett. 2015, 115, 180403. [Google Scholar] [CrossRef] [PubMed]
  2. Zingales, M.; Failla, G. The finite element method for fractional non-local thermal energy transfer in non-homogeneous rigid conductors. Commun. Nonlinear Sci. Numer. Simul. 2015, 29, 116–127. [Google Scholar] [CrossRef] [Green Version]
  3. Lazopoulos, K.A.; Lazopoulos, A.K. On fractional bending of beams. Arch. Appl. Mech. 2016, 86, 1133–1145. [Google Scholar] [CrossRef]
  4. Sumelka, W.; Voyiadjis, G.Z. A hyperelastic fractional damage material model with memory. Int. J. Solids Struct. 2017, 124, 151–160. [Google Scholar] [CrossRef]
  5. Voyiadjis, G.Z.; Sumelka, W. Brain modelling in the framework of anisotropic hyperelasticity with time fractional damage evolution governed by the Caputo-Almeida fractional derivative. J. Mech. Behav. Biomed. 2019, 89, 209–216. [Google Scholar] [CrossRef] [PubMed]
  6. Caputo, M.; Ciarletta, M.; Fabrizio, M.; Tibullo, V. Melting and solidification of pure metals by a phase-field model. Rend Lincei-Mat. Appl. 2017, 28, 463–478. [Google Scholar] [CrossRef]
  7. Gómez-Aguilar, J.F. Fractional Meissner–Ochsenfeld effect in superconductors. Phys. Lett. B 2019, 33, 1950316. [Google Scholar] [CrossRef]
  8. Rahimi, Z.; Rezazadeh, G.; Sumelka, W.; Yang, X.-J. A study of critical point instability of micro and nano beams under a distributed variable-pressure force in the framework of the inhomogeneous non-linear nonlocal theory. Arch. Mech. 2017, 69, 413–433. [Google Scholar]
  9. Goodrich, C.S.; Peterson, A.C. Discrete Fractional Calculus; Springer: New York, NY, USA, 2015. [Google Scholar]
  10. Abdeljawad, T. On Riemann and Caputo fractional differences. Comput. Math. Appl. 2011, 62, 1602–1611. [Google Scholar] [CrossRef] [Green Version]
  11. Abdeljawad, T. Dual identities in fractional difference calculus within Riemann. Adv. Differ. Equ. 2013, 2013, 36. [Google Scholar] [CrossRef] [Green Version]
  12. Abdeljawad, T. On delta and nabla Caputo fractional differences and dual identities. Discrete Dyn. Nat. Soc. 2013, 2013, 406910. [Google Scholar] [CrossRef]
  13. Atici, F.M.; Eloe, P.W. Two-point boundary value problems for finite fractional difference equations. J. Differ. Equ. Appl. 2011, 17, 445–456. [Google Scholar] [CrossRef]
  14. Atici, F.M.; Eloe, P.W. A transform method in discrete fractional calculus. Int. J. Differ. Equ. 2007, 2, 165–176. [Google Scholar]
  15. Atici, F.M.; Eloe, P.W. Initial value problems in discrete fractional calculus. Proc. Amer. Math. Soc. 2009, 137, 981–989. [Google Scholar] [CrossRef]
  16. Agarwal, R.P.; Leanu, D.; Rezapour, S.; Salehi, S. The existence of solutions for some fractional finite difference equations via sum boundary conditions. Adv. Differ. Equ. 2014, 2014, 282. [Google Scholar] [CrossRef] [Green Version]
  17. Goodrich, C.S. On discrete sequential fractional boundary value problems. J. Math. Anal. Appl. 2012, 385, 111–124. [Google Scholar] [CrossRef] [Green Version]
  18. Goodrich, C.S. On a discrete fractional three-point boundary value problem. J. Differ. Equ. Appl. 2012, 18, 397–415. [Google Scholar] [CrossRef]
  19. Erbe, L.; Goodrich, C.S.; Jia, B.; Peterson, A. Survey of the qualitative properties of fractional difference operators: Monotonicity, convexity, and asymptotic behavior of solutions. Adv. Differ. Equ. 2016, 2016, 43. [Google Scholar] [CrossRef] [Green Version]
  20. Chen, Y.; Tang, X. Three difference between a class of discrete fractional and integer order boundary value problems. Commun. Nonlinear Sci. 2014, 19, 4057–4067. [Google Scholar] [CrossRef]
  21. Lv, W.; Feng, J. Nonlinear discrete fractional mixed type sum-difference equation boundary value problems in Banach spaces. Adv. Differ. Equ. 2014, 2014, 184. [Google Scholar] [CrossRef] [Green Version]
  22. Lv, W. Existence of solutions for discrete fractional boundary value problems witha p-Laplacian operator. Adv. Differ. Equ. 2012, 2012, 163. [Google Scholar] [CrossRef] [Green Version]
  23. Jia, B.; Erbe, L.; Peterson, A. Two monotonicity results for nabla and delta fractional differences. Arch. Math. 2015, 104, 589–597. [Google Scholar] [CrossRef]
  24. Jia, B.; Erbe, L.; Peterson, A. Convexity for nabla and delta fractional differences. J. Differ. Equ. Appl. 2015, 21, 360–373. [Google Scholar]
  25. Ferreira, R.A.C. Existence and uniqueness of solution to some discrete fractional boundary value problems of order less than one. J. Differ. Equ. Appl. 2013, 19, 712–718. [Google Scholar] [CrossRef]
  26. Ferreira, R.A.C.; Goodrich, C.S. Positive solution for a discrete fractional periodic boundary value problem. Dyn. Contin. Discrete Impuls. Syst. Ser. A Math. Anal. 2012, 19, 545–557. [Google Scholar]
  27. Kang, S.G.; Li, Y.; Chen, H.Q. Positive solutions to boundary value problems of fractional difference equations with nonlocal conditions. Adv. Differ. Equ. 2014, 2014, 7. [Google Scholar] [CrossRef] [Green Version]
  28. Dong, W.; Xu, J.; Regan, D.O. Solutions for a fractional difference boundary value problem. Adv. Differ. Equ. 2013, 2013, 319. [Google Scholar] [CrossRef] [Green Version]
  29. Holm, M. Sum and difference compositions in discrete fractional calculus. Cubo 2011, 13, 153–184. [Google Scholar] [CrossRef] [Green Version]
  30. Sitthiwirattham, T.; Tariboon, J.; Ntouyas, S.K. Existence results for fractional difference equations with three-point fractional sum boundary conditions. Discrete Dyn. Nat. Soc. 2013, 2013, 104276. [Google Scholar] [CrossRef] [Green Version]
  31. Sitthiwirattham, T.; Tariboon, J.; Ntouyas, S.K. Boundary value problems for fractional difference equations with three-point fractional sum boundary conditions. Adv. Differ. Equ. 2013, 2013, 296. [Google Scholar] [CrossRef] [Green Version]
  32. Sitthiwirattham, T. Existence and uniqueness of solutions of sequential nonlinear fractional difference equations with three-point fractional sum boundary conditions. Math. Method. Appl. Sci. 2015, 38, 2809–2815. [Google Scholar] [CrossRef]
  33. Sitthiwirattham, T. Boundary value problem for p-Laplacian Caputo fractional difference equations with fractional sum boundary conditions. Math. Method. Appl. Sci. 2016, 39, 1522–1534. [Google Scholar] [CrossRef]
  34. Chasreechai, S.; Kiataramkul, C.; Sitthiwirattham, T. On nonlinear fractional sum-difference equations via fractional sum boundary conditions involving different orders. Math. Probl. Eng. 2015, 2015, 519072. [Google Scholar] [CrossRef] [Green Version]
  35. Reunsumrit, J.; Sitthiwirattham, T. Positive solutions of three-point fractional sum boundary value problem for Caputo fractional difference equations via an argument with a shift. Positivity 2016, 20, 861–876. [Google Scholar] [CrossRef]
  36. Reunsumrit, J.; Sitthiwirattham, T. On positive solutions to fractional sum boundary value problems for nonlinear fractional difference equations. Math. Method. Appl. Sci. 2016, 39, 2737–2751. [Google Scholar] [CrossRef]
  37. Kaewwisetkul, B.; Sitthiwirattham, T. On nonlocal fractional sum-difference boundary value problems for Caputo fractional functional difference equations with delay. Adv. Differ. Equ. 2017, 2017, 219. [Google Scholar] [CrossRef] [Green Version]
  38. Reunsumrit, R.; Sitthiwirattham, T. A New class of four-point fractional sum boundary value problems for nonlinear sequential fractional difference equations involving shift operators. Kragujevac J. Math. 2018, 42, 371–387. [Google Scholar] [CrossRef] [Green Version]
  39. Chasreechai, S.; Sitthiwirattham, T. Existence results of initial value problems for hybrid fractional sum-difference equations. Discrete Dyn. Nat. Soc. 2018, 2018, 5268528. [Google Scholar] [CrossRef] [Green Version]
  40. Chasreechai, S.; Sitthiwirattham, T. On separate fractional sum-difference boundary value problems with n-point fractional sum-difference boundary conditions via arbitrary different fractional orders. Mathematics. 2019, 2019, 471. [Google Scholar] [CrossRef] [Green Version]
  41. Wu, G.C.; Baleanu, D. Discrete fractional logistic map and its chaos. Nonlinear Dyn. 2014, 75, 283–287. [Google Scholar] [CrossRef]
  42. Wu, G.C.; Baleanu, D. Chaos synchronization of the discrete fractional logistic map. Signal Process. 2014, 102, 96–99. [Google Scholar] [CrossRef]
  43. Wu, G.C.; Baleanu, D.; Xie, H.P.; Chen, F.L. Chaos synchronization of fractional chaotic maps based on stability results. Phys. A 2016, 460, 374–383. [Google Scholar] [CrossRef]
  44. Chen, C.; Bohner, M.; Jia, B. Ulam-Hyers stability of Caputo fractional difference equations. Math. Meth. Appl. Sci. 2019, 42, 7461–7470. [Google Scholar] [CrossRef]
  45. Wu, G.C.; Baleanu, D.; Zeng, S.D. Finite-time stability of discrete fractional delay systems: Gronwall inequality and stability criterion. Commun. Nonlinear Sci. Numer. Simulat. 2018, 57, 299–308. [Google Scholar] [CrossRef]
  46. Baleanu, D.; Wu, G.C.; Bai, Y.R.; Chen, F.L. Stability analysis of Caputo–like discrete fractional systems. Commun. Nonlinear Sci. Numer. Simulat. 2017, 48, 520–530. [Google Scholar] [CrossRef]
  47. Chen, F.; Zhou, Y. Existence and Ulam Stability of Solutions for Discrete Fractional Boundary Value Problem. Discrete Dyn. Nat. Soc. 2013, 2013, 459161. [Google Scholar] [CrossRef]
  48. Cermák, J.; Gőrigyori, I.; Nechvátal, L. On explicit stability conditions for a linear fractional difference system. Fract. Calc. Appl. Anal. 2015, 18, 651–672. [Google Scholar] [CrossRef]
  49. Lizama, C. The Poisson distribution, abstract fractional difference equations, and stability. Proc. Am. Math. Soc. 2017, 145, 3809–3827. [Google Scholar] [CrossRef]
  50. Cermák, J.; Kisela, T.; Nechvátal, L. Stability and asymptotic properties of a linear fractional difference equation. Adv. Differ. Equ. 2012, 2012, 122. [Google Scholar] [CrossRef] [Green Version]
  51. Dassios, I.K. Stability and Robustness of Singular Systems of Fractional Nabla Difference Equations. Circuits Syst. Signal Process 2017, 36, 49–64. [Google Scholar] [CrossRef] [Green Version]
  52. Wu, G.C.; Baleanu, D. Stability analysis of impulsive fractional difference equations. Fract. Calc. Appl. Anal. 2018, 21, 354–375. [Google Scholar] [CrossRef]
  53. Wu, G.C.; Baleanu, D.; Huang, L.L. Novel Mittag-Leffler stability of linear fractional delay difference equations with impulse. Appl. Math. Lett. 2018, 82, 71–78. [Google Scholar] [CrossRef]
  54. Wang, J.; Lv, L.; Zhou, Y. Ulam stability and data dependence for fractional differential equations with Caputo derivative. Electron. J. Qual. Theory Differ. Equ. 2011, 63, 1–10. [Google Scholar] [CrossRef]

Share and Cite

MDPI and ACS Style

Ouncharoen, R.; Chasreechai, S.; Sitthiwirattham, T. Existence and Stability Analysis for Fractional Impulsive Caputo Difference-Sum Equations with Periodic Boundary Condition. Mathematics 2020, 8, 843. https://doi.org/10.3390/math8050843

AMA Style

Ouncharoen R, Chasreechai S, Sitthiwirattham T. Existence and Stability Analysis for Fractional Impulsive Caputo Difference-Sum Equations with Periodic Boundary Condition. Mathematics. 2020; 8(5):843. https://doi.org/10.3390/math8050843

Chicago/Turabian Style

Ouncharoen, Rujira, Saowaluck Chasreechai, and Thanin Sitthiwirattham. 2020. "Existence and Stability Analysis for Fractional Impulsive Caputo Difference-Sum Equations with Periodic Boundary Condition" Mathematics 8, no. 5: 843. https://doi.org/10.3390/math8050843

APA Style

Ouncharoen, R., Chasreechai, S., & Sitthiwirattham, T. (2020). Existence and Stability Analysis for Fractional Impulsive Caputo Difference-Sum Equations with Periodic Boundary Condition. Mathematics, 8(5), 843. https://doi.org/10.3390/math8050843

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