Next Article in Journal
Potential Subgraphs of the Missing Moore Graph
Previous Article in Journal
Three-Dimensional Unmanned Aerial Vehicle Trajectory Planning Based on the Improved Whale Optimization Algorithm
Previous Article in Special Issue
Closed-Form Exact Solution for Free Vibration Analysis of Symmetric Functionally Graded Beams
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improved Fractional Differences with Kernels of Delta Mittag–Leffler and Exponential Functions

by
Miguel Vivas-Cortez
1,
Pshtiwan Othman Mohammed
2,*,
Juan L. G. Guirao
3,
Majeed A. Yousif
4,
Ibrahim S. Ibrahim
4 and
Nejmeddine Chorfi
5
1
Faculty of Exact and Natural Sciences, School of Physical Sciences and Mathematics, Pontifical Catholic University of Ecuador, Av. 12 de Octubre 1076 y Roca, Quito 17-01-2184, Ecuador
2
Department of Mathematics, College of Education, University of Sulaimani, Sulaimani 46001, Iraq
3
Department of Applied Mathematics and Statistics, Technical University of Cartagena, Hospital de Marina, 30203 Cartagena, Spain
4
Department of Mathematics, College of Education, University of Zakho, Zakho 42002, Iraq
5
Department of Mathematics, College of Science, King Saud University, P.O. Box 2455, Riyadh 11451, Saudi Arabia
*
Author to whom correspondence should be addressed.
Symmetry 2024, 16(12), 1562; https://doi.org/10.3390/sym16121562
Submission received: 25 October 2024 / Revised: 8 November 2024 / Accepted: 13 November 2024 / Published: 21 November 2024
(This article belongs to the Special Issue Symmetry and Asymmetry in Nonlinear Partial Differential Equations)

Abstract

:
Special functions have been widely used in fractional calculus, particularly for addressing the symmetric behavior of the function. This paper provides improved delta Mittag–Leffler and exponential functions to establish new types of fractional difference operators in the setting of Riemann–Liouville and Liouville–Caputo. We give some properties of these discrete functions and use them as the kernel of the new fractional operators. In detail, we propose the construction of the new fractional sums and differences. We also find the Laplace transform of them. Finally, the relationship between the Riemann–Liouville and Liouville–Caputo operators are examined to verify the feasibility and effectiveness of the new fractional operators.

1. Introduction

The discrete fractional operator theory in fractional calculus settings plays an important role in solving difference equations. Fractional difference models are a core issue in mathematical and numerical analysis. It is also widely applied in many fields, including dynamical systems, artificial intelligence, computer communication, compressed sensing and so on, e.g., see [1,2,3].
In particular, mathematical and physical models of continuous and discrete fractional calculus have become important tools for analyzing differential and difference equations [4,5,6,7]. However, many researchers have proposed and studied the existence and uniqueness of discrete fractional models in different systems [8,9,10,11,12].
Various fractional operators have been used worldwide to study mathematical models in the presence of continuous fractional differential equations, for example [13,14,15,16], and discrete fractional difference equations [17,18,19,20] for nabla and delta types in the sense of Riemann–Liouville (RL), Liouville–Caputo (LC), Caputo–Fabrizio (CF) and Attangana–Baleanu (AB) operators.
Discrete fractional calculus has several correlations with special functions. An important application of discrete fractional calculus in applied mathematics is proving new identities and relationships between special functions. Specifically, enormous efforts have been made to accurately calculate special functions, which can be challenging for researchers who model real-world problems in discrete times. A function with particularly strong ties to discrete fractional calculus is the discrete Mittag–Leffler (ML) function—see, e.g., [21,22,23,24,25]. Furthermore, different authors have developed the problem of discrete fractional calculus in the time scale Z to the time scale h Z , see [26,27,28].
This study proposes discrete Δ h -CF and Δ h -AB fractional operators, including delta exponential and ML functions, in their kernels. The key contributions of the proposed research article are summarized as follows:
  • A new delta exponential and ML function with some special cases are proposed in the context of h -discrete fractional calculus.
  • Some properties of the proposed delta ML functions have been addressed by utilizing Laplace transformations and sum/difference rules.
  • The concept of the Δ h -CF and Δ h -AB fractional sums and differences are introduced based on the new delta exponential and Mittag–Leffler functions.
  • In conclusion, some properties of the Δ h -CF and Δ h -AB fractional difference operators are provided based on the Laplace transformation on the left and right sides.
In addition, the article has the following sections: in Section 2, essential definitions in discrete fractional calculus and Laplace transformation with some of its properties are recalled. In Section 3, we formulate the new h -ML functions; we divide this section into two subsections: Δ h -CF fractional operators in Section 3.1 and Δ h -AB fractional operators in Section 3.2. Finally, some discussions finish this article in Section 4.

2. Preliminaries

Let α , h > 0 . Then, we consider the notations N ξ 0 , h : = ξ 0 + h N and N ξ , h : = ξ h N , for ξ 0 , ξ R . Furthermore, let T h : = N ξ 0 ; ξ h = { ξ 0 , ξ 0 + h , , ξ } such that ξ = ξ 0 + k h , for some k N . From Definition 2.25 in [1], we recall the definition of Δ h fractional sums as follows:
Δ ξ 0 h α v ( τ ) = = ξ 0 h τ h α α 1 , h ( τ , σ h ( ) ) v ( ) , τ N ξ 0 + α h , h ,
where v : N ξ 0 , h R , σ h ( τ ) : = τ + h and
α , h ( τ , ) = ( τ ) h ( α ) Γ ( α + 1 ) with ( τ ) h ( α ) = Γ τ h + 1 Γ τ h + 1 α ,
whenever Γ τ h + 1 Γ τ h + 1 α is well-defined for τ , α R .
Lemma 1 
(Theorem 2.50 in [1]). If v is defined on N ξ 0 , h , then for α > 0 , it can be expressed that
Δ ξ 0 h α Δ h v ( τ ) = Δ h Δ h α ξ 0 v ( τ ) ( τ ξ 0 ) h ( α 1 ) Γ ( α ) v ( ξ 0 ) ,
for τ N ξ 0 + α h , h .
Remark 1. 
For each α R and h > 0 , we have
(i) 
Δ h τ h ( α ) = α τ h ( α 1 ) .
(ii) 
Δ ξ 0 + μ h h α ( τ ξ 0 ) h ( μ ) = Γ ( μ + 1 ) Γ ( μ + 1 + α ) ( τ ξ 0 ) h ( α + μ ) , for α > 0 , μ 0 and τ N ξ 0 + ( α + μ ) h , h .
Lemma 2 
([26,27]). Let α > 0 . Then, for τ T h , we have
lim α 0 Δ h α ξ 0 v ( τ + α h ) = v ( τ ) , lim α 0 Δ ξ α h v ( τ α h ) = v ( τ ) .
The following definitions are given in [26].
Definition 1. 
For the function v defined on N ξ 0 , h and N ξ , h , the h-RL fractional difference on the left side is defined by
Δ h α ξ 0 v ( τ ) = Δ h ı Δ h ( ı α ) ξ 0 v ( τ ) , τ N ξ 0 + ( ı α ) h , h ,
and on the right side by
Δ ξ α h v ( τ ) = ( 1 ) ı Δ h ı Δ ξ ( ı α ) h v ( τ ) , τ N ξ ( ı α ) h , h ,
where ı = [ α ] + 1 .
Definition 2. 
For a function v defined on N ξ 0 , h , the h-LC fractional difference on the left side is defined by
Δ h α ξ 0 L C v ( τ ) = Δ h ( ı α ) ξ 0 Δ h ı v ( τ ) , τ N ξ 0 + ( ı α ) h , h ,
and on the right side by
Δ ξ α h L C v ( τ ) = Δ ξ ( ı α ) h ( 1 ) ı Δ h ı v ( τ ) , τ N h ξ ( ı α ) h , ,
where ı = [ α ] + 1 .
The operator of Q is defined by
Q v ( τ ) = v ( ξ 0 + ξ τ ) ,
and was used in [29,30] to obtain the right sides of fractional sum and difference operators by knowing the left sides of these operators without conducting their proofs.
Lemma 3 
([26]). Let α , h > 0 and v be defined on T h . Then,
(i) 
Δ h α ξ 0 Q v ( τ ) = Q Δ ξ α h v ( τ ) = Δ ξ α h v v ( ξ 0 + ξ τ ) .
(ii) 
Δ h α ξ 0 Q v ( τ ) = Q Δ ξ α h v ( τ ) = Δ ξ α h v v ( ξ 0 + ξ τ ) .
(iii) 
Δ h α ξ 0 L C Q v ( τ ) = Q Δ ξ α h L C v ( τ ) = Δ ξ α h L C v v ( ξ 0 + ξ τ ) .
According to Theorem 2.2 in the monograph [1] and Definition 14 in [26], we can introduce the following definition of the Δ h -Laplace transformation.
Definition 3. 
For the function v defined on N ξ 0 , h , the Δ h -Laplace transformation of v is defined by
L ξ 0 , h { v ( τ ) } ( s ) = ξ 0 e ˜ s h ( τ , ξ 0 ) h s + 1 v ( τ ) Δ h τ = h k = 0 v ( ξ 0 + k h ) ( h s + 1 ) k + 1 ,
where e ˜ s h is defined later in Section 3.1.
Definition 4. 
The Δ h -convolution of two functions v and g defined on N ξ 0 , h is given by
v g ( τ ) = ξ 0 τ g ( τ σ ( s ) + ξ 0 ) ) v ( s ) Δ h s = h k = 0 v ( k h ) g ( τ σ ( k h ) + ξ 0 ) ) .
Furthermore, the Laplace convolution of these functions can be expressed by
L ξ 0 , h ( v g ) ( τ ) ( s ) = L { v ( τ ) } ( s ) L { g ( τ ) } ( s ) ,
for s R and τ N ξ 0 , h .
Lemma 4. 
For v on N ξ 0 , h , we have
L ξ 0 , h { Δ h v ( τ ) } ( s ) = s L { v ( τ ) } ( s ) v ( ξ 0 ) .
Lemma 5 
([1]). Let ξ 0 , h 0 v and α R { 1 , 2 , } . Then, for h s + 1 < 1 , it can be expressed that
L ξ 0 + h α , h { H ˜ α h ( τ , ξ 0 ) } ( s ) = ( h s + 1 ) α s α + 1 ,
where H ˜ α h ( τ , ξ 0 ) = ( τ ξ 0 ) h ( α ) Γ ( α + 1 ) .
Lemma 6 
([1]). Let M N 0 and v : N ξ 0 M h , h R and g : N ξ 0 , h R be two functions of exponential order > 0 . Then, for h s + 1 > , it can be expressed that
L ξ 0 M h α , h { v ( τ ) } ( s ) = 1 ( h s + 1 ) M L ξ 0 , h { v ( τ ) } ( s ) + h j = 0 M 1 v ( ξ 0 + ( j M ) h ) ( h s + 1 ) j + 1 ,
and
L ξ 0 + M h α , h { g ( τ ) } ( s ) = ( h s + 1 ) M L ξ 0 , h { g ( τ ) } ( s ) h j = 0 M 1 ( h s + 1 ) M j 1 g ( ξ 0 + j h ) .

3. Delta ML Functions and Related Operators

In this part of our article, we state the delta h -ML function and derive the CF and AB fractional differences and sums including delta h -ML function in their kernel.
Definition 5. 
For any complex numbers α , μ γ , λ with ( α ) > 0 and λ h α < 1 , we express the delta h -ML function as follows:
E ( α , μ ) γ h ( λ , τ ) = ı = 0 λ ı ( τ + ı ( α 1 ) + ( μ 1 ) h ) h ( ı α + μ 1 ) ( γ ) ı Γ ( α ı + μ ) .
From (7), the following special cases can be observed:
  • If γ = 1 , then (7) becomes
    E ( α , μ ) h ( λ , τ ) = ı = 0 λ ı ( τ + ı ( α 1 ) h + ( μ 1 ) h ) h ( ı α + μ 1 ) Γ ( α τ + μ ) .
  • If γ = μ = 1 , then (7) becomes
    E ( α ) h ( λ , τ ) = ı = 0 λ ı ( τ + ı ( α 1 ) h ) h ( ı α ) Γ ( α ı + 1 ) .
For the main concepts about the discrete Mittag–Liffler functions and their properties, we refer the reader to [21,22,23,24,25,26,27].
Remark 2. 
One of the motivations make us to work on h Z is that a smaller h ( 0 , 1 ) allows us to use a larger interval of λ. For example, if h = 1 , we need 0 < α < 1 2 to guarantee convergence as λ = 1 1 α . However, as h 0 , we can attain that the series will be convergent for 0 < α < 1 .

3.1. Delta Fractional Differences with Exponential Kernels

Let us denote the Δ h exponential kernel by
e ˜ λ h ( τ , σ ( ) ) = ( 1 + λ h ) τ σ ( ) h = 1 α α h 1 α τ σ ( r ) h ,
where λ = α 1 α with λ h α < 1 . In addition, for s = 1 + h s s , it follows that
e ˜ s h ( τ , ξ 0 ) = 1 1 + h s τ ξ 0 h , τ N ξ 0 , h .
Definition 6. 
For 0 < α < 1 and with v defined on N ξ 0 , h in the left case and N ξ , h in the right case. Let D ( α ) be the normalization, as defined in [26]. Then, we define the following:
(a) 
The Δ h -CF of Liouville–Caputo type on the left side by
Δ h α ξ 0 C F L C v ( τ ) = D ( α ) 1 α = ξ 0 h τ h 1 h Δ h v ( h ) ( 1 + λ h ) τ σ ( r ) h = D ( α ) 1 α = ξ 0 h τ h 1 h Δ h v ( h ) 1 α α h 1 α τ σ ( r ) h ,
for τ N ξ 0 + h , h .
(b) 
The Δ h -CF of Liouville–Caputo type on the right side by
Δ ξ α h C F L C v ( τ ) = D ( α ) 1 α = τ h + 1 ξ h h h v ( h ) ( 1 + λ h ) h σ ( τ ) h = D ( α ) 1 α = τ h + 1 ξ h h h v ( h ) 1 α α h 1 α h σ ( τ ) h ,
for τ N ξ h , h .
(c) 
The Δ h -CF of Riemann–Liouville type on the left side by
Δ h α ξ 0 C F R L v ( τ ) = D ( α ) 1 α Δ h = ξ 0 h τ h 1 h v ( h ) ( 1 + λ h ) h σ ( τ ) h = D ( α ) 1 α Δ h = ξ 0 h τ h 1 h v ( h ) 1 α α h 1 α h σ ( τ ) h ,
for τ N ξ 0 + h , h .
(d) 
The Δ h -CF of Riemann–Liouville type on the right side as follows:
Δ ξ α h C F R L v ( τ ) = D ( α ) 1 α h = τ h + 1 ξ h h v ( h ) ( 1 + λ h ) h σ ( τ ) h = D ( α ) 1 α h = τ h + 1 ξ h h v ( h ) 1 α α h 1 α h σ ( τ ) h ,
for τ N ξ h , h .
Remark 3. 
Considering the above definition, we can note the following:
(i) 
As α 0 , we have
lim α 0 Δ h α ξ 0 C F L C v ( τ ) v ( τ ) v ( ξ 0 ) , and lim α 0 Δ ξ α h C F L C v ( τ ) v ( τ ) v ( ξ ) , lim α 0 Δ h α ξ 0 C F R L v ( τ ) v ( τ ) , and lim α 0 Δ ξ α h C F R L v ( τ ) v ( τ ) .
(ii) 
As α 1 , we have
lim α 1 Δ h α ξ 0 C F L C v ( τ ) Δ v ( τ ) , and lim α 1 Δ ξ α h C F L C v ( τ ) v ( τ ) , lim α 1 Δ h α ξ 0 C F R L v ( τ ) Δ v ( τ ) , and lim α 1 Δ ξ α h C F R L v ( τ ) v ( τ ) .
Remark 4. 
The expression of Q can be applied to the above operators to have
(i) 
Q ξ 0 C F L C Δ h α v ( τ ) = ( h C F L C Δ ξ α Q v ) ( τ ) ,
(ii) 
Q ξ 0 C F R L Δ h α v ( τ ) = ( h C F R L Δ ξ α Q v ) ( τ ) .
Consider the difference equation
( ξ 0 C F R L Δ h α v ) ( τ ) = g ( τ ) ,
or, equivalently,
D ( α ) 1 α Δ h v e ˜ λ h ( τ , ξ 0 ) = g ( τ ) .
By applying L ξ 0 , h on both sides of (15) and considering (4), Lemma 4 and
L ξ 0 , h { h e ˜ λ ( τ , ξ 0 ) } ( s ) = 1 s λ = ( 1 α ) ( 1 α ) s + α ,
we obtain
s D ( α ) ( 1 α ) s + α V ( s ) = G ( s ) ,
where G ( s ) = L ξ 0 , h { g ( τ ) } ( s ) and V ( s ) = L ξ 0 , h { v ( τ ) } . This leads to
V ( s ) = 1 α D ( α ) G ( s ) + α D ( α ) 1 s G ( s ) .
Apply L ξ 0 , h 1 on both sides to reach
v ( τ ) = 1 α D ( α ) g ( τ ) + α D ( α ) k = ξ 0 h τ h 1 h g ( k h ) .
This helps us introduce the following CF fractional definitions as follows:
Definition 7. 
For α ( 0 , 1 ) and v defined on T h , we express the Δ h -CF fractional sum in the left side case by
Δ h α ξ 0 C F v ( τ ) = 1 α D ( α ) v ( τ ) + α D ( α ) k = ξ 0 h τ h 1 h v ( k h ) ,
and in the right side case by
Δ ξ α h C F v ( τ ) = 1 α D ( α ) v ( τ ) + α D ( α ) k = τ h + 1 ξ h h v ( k h ) .
It can be easily shown that Δ h α ξ 0 C F Q v ( τ ) = Q Δ ξ α h C F v ( τ ) . Furthermore, we have Δ ξ α h C F h C F R L Δ ξ α v ( τ ) = v ( τ ) .
Conversely, consider the difference equation
Δ h α ξ 0 C F v ( τ ) = 1 α D ( α ) v ( τ ) + α D ( α ) k = ξ 0 h τ h 1 h v ( k h ) = g ( τ ) .
Applying L ξ 0 , h on both sides to obtain
1 α D ( α ) V ( s ) + α s D ( α ) V ( s ) = G ( s ) .
leads to
V ( s ) = s D ( α ) ( 1 α ) s + α G ( s ) = ( 16 ) b y L ξ 0 , h Δ h α ξ 0 C F R L g ( τ ) ( s ) = L ξ 0 , h Δ h α ξ 0 C F R L ξ 0 C F Δ h α v ( τ ) ( s ) .
It follows that
Δ h α ξ 0 C F R L ξ 0 C F Δ h α v ( τ ) = v ( τ ) .
By the same method, or by applying the operator of Q , we can show that
Δ ξ α h C F h C F R L Δ ξ α v ( τ ) = v ( τ ) .
It is always of interest to make a connection between Δ h -CF fractional differences of RL and LC types. For this reason, the following theorem will be raised.
Theorem 1. 
For the function v on N ξ 0 , h , we have the following relationships:
(i) 
( ξ 0 C F L C Δ h α v ) ( τ ) = ( ξ 0 C F R L Δ h α v ) ( τ ) D ( α ) 1 α v ( ξ 0 ) 1 α α h 1 α τ ξ 0 h ,
(ii) 
Δ ξ α h C F L C v ( τ ) = Δ ξ α h C F R L v ( τ ) D ( α ) 1 α v ( ξ ) 1 α α h 1 α ξ τ h .
Proof. 
Having (16), we see that
L ξ 0 , h { ξ 0 C F R L Δ h α v ( τ ) } ( s ) = s D ( α ) ( 1 α ) s + α V ( s ) .
By applying L ξ 0 , h on (11) with the use of (19), we obtain
L ξ 0 , h { ξ 0 C F L C Δ h α v ( τ ) } ( s ) = D ( α ) 1 α L ξ 0 , h ( Δ v ( τ ) ) e ˜ λ h ( τ , ξ 0 ) ( s ) = D ( α ) 1 α L ξ 0 , h Δ v ( τ ) ( s ) · L ξ 0 , h e ˜ λ h ( τ , ξ 0 ) ( s ) = D ( α ) 1 α s V ( s ) s α D ( α ) 1 α ( 1 α ) ( 1 α ) s + α v ( ξ 0 ) = L ξ 0 , h ( ξ 0 C F R L Δ h α v ) ( τ ) ( s ) D ( α ) 1 α ( 1 α ) ( 1 α ) s + α v ( ξ 0 ) .
Applying L ξ 0 , h 1 on both sides of (20), we obtain (i) as desired. The proof of (ii) will be obtained by applying the operator of Q on the first item (i). □

3.2. Delta Fractional Differences with ML Kernels

This subsection is dedicated to introduce the Δ h -AB fractional differences and sums, and some of their properties. We start by introducing the Δ h -AB fractional difference operators.
Definition 8. 
For the function v on T h and 0 < α 1 with λ h α < 1 , we express the Δ h -AB fractional difference of Liouville–Caputo type on the left side by
Δ h α ξ 0 A B L C v ( τ ) = D ( α ) 1 α = ξ 0 h τ h 1 h Δ h v ( h ) E ( α ) h ( λ , τ σ ( h ) ) = D ( α ) 1 α Δ h v ( τ ) E ( α ) h ( λ , τ ξ 0 ) ,
and of Riemann–Liouville type on the left side by
Δ h α ξ 0 A B R L v ( τ ) = D ( α ) 1 α Δ h = ξ 0 h τ h 1 h v ( h ) E ( α ) h ( λ , τ σ ( h ) ) = D ( α ) 1 α Δ h v ( τ ) E ( α ) h ( λ , τ ξ 0 ) .
Definition 9. 
For the function v on T h and 0 < α 1 with λ h α < 1 , we express the Δ h -AB fractional difference of Liouville–Caputo type on the right side by
Δ ξ α h A B L C v ( τ ) = D ( α ) 1 α = τ h + 1 ξ h h h v ( h ) E ( α ) h ( λ , σ ( h ) τ ) ,
and of Riemann–Liouville type on the right side by
Δ ξ α h A B R L v ( τ ) = D ( α ) 1 α h = τ h + 1 ξ h h v ( h ) E ( α ) h ( λ , σ ( h ) τ ) .
To find the Laplace of Δ h α ξ 0 v ( τ + ( α 1 ) h ) , we need the following lemma on h -ML functions.
Lemma 7. 
For 0 < α 1 , we have
L ξ 0 , h E ( α ) h ( λ , τ ξ 0 ) ( s ) = s α 1 ( h s + 1 ) 1 α s α ( h s + 1 ) 1 α λ ,
provided that λ ( h s + 1 ) α 1 < s α .
Proof. 
From Definitions 3 and 5, and making use of Lemmas 6 and 7, we see that
L ξ 0 , h E ( α ) h ( λ , τ ξ 0 ) ( s ) = ı = 0 λ ı L ξ 0 , h ( τ ξ 0 + ı ( α 1 ) h ) h ( α ı ) Γ ( α ı + 1 ) ( s ) = ı = 0 λ ı [ 1 ( h s + 1 ) ı L ξ 0 + ı h , h ( τ ξ 0 + ı ( α 1 ) h ) h ( α ı ) Γ ( α ı + 1 ) ( s ) + h ( h s + 1 ) ı Γ ( α ı + 1 ) j = 0 ı 1 ( h s + 1 ) ı j 1 ( ( j + ı ( α 1 ) ) h ) h ( α ı ) ] = ı = 0 λ ı [ 1 ( h s + 1 ) ı ( h s + 1 ) α ı s α ı + 1 + h α ı + 1 j = 0 ı 1 ( h s + 1 ) j 1 Γ ( j + ı ( α 1 ) + 1 ) Γ ( j ı + 1 ) ] = ı = 0 λ ı 1 ( h s + 1 ) ı ( h s + 1 ) α ı s α ı + 1 + 0 = 1 s ı = 0 λ ( h s + 1 ) α 1 s α ı .
This will be convergent if λ ( h s + 1 ) α 1 < s α . Therefore, by the test of geometric series, we obtain
L ξ 0 , h E ( α ) h ( λ , τ ξ 0 ) ( s ) = 1 s 1 1 λ ( h s + 1 ) α 1 s α = s α 1 ( h s + 1 ) 1 α s α ( h s + 1 ) 1 α λ ,
which completes the proof. □
Lemma 8. 
For each α > 0 , it can be
L ξ 0 , h Δ h α ξ 0 v ( τ + ( α 1 ) h ) ( s ) = ( h s + 1 ) α 1 s α V ( s ) ,
where V ( s ) = L ξ 0 , h v ( s ) .
Proof. 
In view of (1) and (3), we have
Δ h α ξ 0 v ( τ + ( α 1 ) h ) = 1 Γ ( α ) k = ξ 0 h τ h 1 ( τ σ ( k h ) ) h ( α 1 ) v ( k h ) h = k = ξ 0 h τ h 1 H ˜ α 1 h ( τ + ( α 1 ) h , σ ( k h ) ) v ( k h ) h = h k = ξ 0 h τ h 1 v ( k h ) H ˜ α 1 h ( τ + ( α 1 ) h σ ( k h ) + ξ 0 , ξ 0 ) = v H ˜ α 1 h ( τ , ξ 0 ) ( τ + ( α 1 ) h ) ,
where H ˜ α 1 h ( τ , ξ 0 ) is as defined in Lemma 5. Therefore, by considering Lemma 5, we obtain
L ξ 0 , h Δ h α ξ 0 v ( τ + ( α 1 ) h ) ( s ) = L ξ 0 , h v ( τ ) ( s ) · L ξ 0 , h H ˜ α 1 h ( τ + ( α 1 ) h , ξ 0 ) ( s ) = V ( s ) ( h s + 1 ) α 1 s α ,
which ends the proof. □
In association to Δ h α ξ 0 A B R L u ( τ ) = v ( τ ) , we can formulate the Δ h -AB fractional sum by considering the difference equation
Δ h α ξ 0 A B R L g ( τ ) = v ( τ ) .
To solve (25), we can apply L ξ 0 , h on both sides and apply (22) and Lemma 7 to obtain
D ( α ) 1 α s α 1 ( h s + 1 ) 1 α s α ( h s + 1 ) 1 α λ G ( s ) = V ( s ) .
From this, we conclude that
G ( s ) = 1 α D ( α ) V ( s ) + α D ( α ) 1 s α ( h s + 1 ) 1 α V ( s ) .
Then, we apply L ξ 0 , h 1 and use Lemma 8 to deduce that
g ( τ ) = 1 α D ( α ) v ( τ ) + α D ( α ) Δ h α ξ 0 v ( τ + ( α 1 ) h ) .
Conversely, if we solve Δ h α ξ 0 A B g ( τ ) = v ( τ ) , we find out that g ( τ ) = Δ h α ξ 0 A B R L v ( τ ) .
Depending on (27), we can define the Δ h -AB fractional sums on the left and right sides as follows.
Definition 10. 
Let h > 0 and 0 < α < 1 . We express the Δ h -AB fractional sum on the left side by
Δ h α ξ 0 A B v ( τ ) = 1 α D ( α ) v ( τ ) + Δ h α ξ 0 v ( τ + ( α 1 ) h ) ,
and on the right side by
Δ ξ α h A B v ( τ ) = 1 α D ( α ) v ( τ ) + Δ ξ α h v ( τ ( α 1 ) h ) .
It is worth noting that Δ h α ξ 0 A B Q v ( τ ) = Q Δ ξ α h A B v ( τ ) . In addition, by considering the operator of Q , it can be proved that ( Δ ξ α h A B h A B R L Δ ξ α v ) ( τ ) = v ( τ ) .
In the last stage of this section, it is worth making a connection between Δ h -AB fractional differences of RL and LC types. For this reason, we have the following example.
Example 1. 
Suppose that v is defined on N ξ 0 , h . Then, in the conclusion of (21), Lemmas 4 and 7, we have
L ξ 0 , h { ( ξ 0 A B L C Δ h α v ) ( τ ) } ( s ) = D ( α ) 1 α ( s V ( s ) v ( ξ 0 ) ) s α 1 ( h s + 1 ) 1 α s α ( h s + 1 ) 1 α λ = D ( α ) 1 α s α ( h s + 1 ) 1 α s α ( h s + 1 ) 1 α λ V ( s ) D ( α ) 1 α s α 1 ( h s + 1 ) 1 α s α ( h s + 1 ) 1 α λ v ( ξ 0 ) .
By applying L ξ 0 , h 1 on both sides and considering (26), we can conclude that we have the following relationships:
( ξ 0 A B L C Δ h α v ) ( τ ) = ( ξ 0 A B R L Δ h α v ) ( τ ) D ( α ) 1 α v ( ξ 0 ) E ( α ) h ( λ , τ ξ 0 ) .
By applying the operator of Q on (30), we can similarly deduce that
Δ ξ α h A B L C v ( τ ) = Δ ξ α h A B R L v ( τ ) D ( α ) 1 α v ( ξ ) E ( α ) h ( λ , ξ τ ) .
Remark 5. 
The application of our theoretical results could not be prepared easily in a simple example, which is why we have not prepared any application in this paper. For this reason, we refer the reader to see the recently accepted paper [31], in which the application of our new findings has been examined in the context of monotonic analysis.

4. Conclusions

This work emphasizes the development of novel special functions named as delta h -exponential and h -ML functions to address the Δ h -CF and Δ h -AB fractional difference operators including these functions in their kernels. We aim to reconstruct discrete fractional operators in the delta background sense. Afterwards, Laplace transformation is performed on the Δ h -CF and Δ h -AB fractional difference operators to obtain their corresponding Δ h -CF and Δ h -AB fractional sum operators. Relationships between the RL and LC cases are considered for the Δ h -CF and Δ h -AB fractional difference operators.
It can be noted that there may be some limitations in our study. This study only introduces one variate delta Mittag–Leffler function. Bivariate discrete Mittag–Leffler functions [32], higher-order differential equations [33] and monotonicity analysis are also needed to generalize the above new operators and address future work.

Author Contributions

Conceptualization, J.L.G.G.; Data curation, I.S.I.; Funding acquisition, M.V.-C.; Investigation, P.O.M.; Methodology, M.A.Y.; Project administration, J.L.G.G.; Resources, M.V.-C.; Software, I.S.I.; Supervision, N.C.; Visualization, N.C.; Writing—review and editing, P.O.M. and M.A.Y. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Acknowledgments

Researchers Supporting Project number (RSP2024R153), King Saud University, Riyadh, Saudi Arabia.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Goodrich, C.; Peterson, A.C. Discrete Fractional Calculus; Springer: New York, NY, USA, 2015. [Google Scholar]
  2. Cabada, A.; Dimitrov, N. Nontrivial solutions of non-autonomous Dirichlet fractional discrete problems. Fract. Calc. Appl. Anal. 2020, 23, 980–995. [Google Scholar] [CrossRef]
  3. Zhang, X.-L.; Yu, Y.; Wang, H.; Feng, J. Stability analysis of discrete-time tempered fractional-order neural networks with time delays. Fract. Calc. Appl. Anal. 2024, 27, 1972–1993. [Google Scholar] [CrossRef]
  4. Atici, F.M.; Atici, M.; Belcher, M.; Marshall, D. A new approach for modeling with discrete fractional equations. Fund. Inform. 2017, 151, 313–324. [Google Scholar]
  5. You, X.; Song, Q.; Zhao, Z. Existence and finite-time stability of discrete fractional-order complex-valued neural networks with time delays. Neural Netw. 2020, 123, 248–260. [Google Scholar] [CrossRef]
  6. Wei, Y.; Chen, Y.; Wei, Y.; Zhao, X. Lyapunov stability analysis for nonlinear nabla tempered fractional order systems. Asian J. Control 2023, 25, 3057–3066. [Google Scholar] [CrossRef]
  7. Ortigueira, M. Discrete-time fractional difference calculus: Origins, evolutions, and new formalisms. Fractal Fract. 2023, 7, 502. [Google Scholar] [CrossRef]
  8. Jia, B.; Du, F.; Erbe, L.; Peterson, A. Asymptotic behavior of nabla half order h-difference equations. J. Appl. Anal. Comput. 2018, 8, 1707–1726. [Google Scholar]
  9. Liu, X.; Yu, Y. Synchronization analysis for discrete fractional-order complex-valued neural networks with time delays. Neural Comput. Appl. 2021, 33, 10503–10514. [Google Scholar] [CrossRef]
  10. Mohammed, P.O.; Srivastava, H.M.; Baleanu, D.; Al-Sarairah, E.; Yousif, M.A.; Chorfi, N. Analytical and approximate monotone solutions of the mixed order fractional nabla operators subject to bounded conditions. Math. Comp. Model. Dyn. 2024, 30, 626–639. [Google Scholar] [CrossRef]
  11. Mohammed, P.O. An analysis of exponential kernel fractional difference operator for delta positivity. Nonlinear Eng. 2024, 13, 20220377. [Google Scholar] [CrossRef]
  12. Chen, C.R.; Bohner, M.; Jia, B.G. Ulam-Hyers stability of Caputo fractional difference equations. Math. Meth. Appl. Sci. 2019, 42, 7461–7470. [Google Scholar] [CrossRef]
  13. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; Wiley-Interscience Publisher: New York, NY, USA, 1993. [Google Scholar]
  14. Cesarano, C. Generalized special functions in the description of fractional diffusive equations. Commun. Appl. Ind. Math. 2019, 10, 31–40. [Google Scholar] [CrossRef]
  15. Jajarmi, A.; Baleanu, D.; Sajjadi, S.S.; Nieto, J.J. Analysis and some applications of a regularized Ψ-Hilfer fractional derivative. J. Comput. Appl. Math. 2022, 415, 114476. [Google Scholar] [CrossRef]
  16. Baleanu, D.; Jajarmi, A.; Hajipour, M. On the nonlinear dynamical systems within the generalized fractional derivatives with Mittag–Leffler Kernel. Nonlinear Dyn. 2018, 94, 397–414. [Google Scholar] [CrossRef]
  17. Abdeljawad, T. On Riemann and Caputo fractional differences. Comput. Math. Appl. 2011, 62, 1602–1611. [Google Scholar] [CrossRef]
  18. Abdeljawad, T.; Baleanu, D. Discrete fractional differences with nonsingular discrete Mittag–Leffler kernels. Adv. Differ. Equ. 2016, 2016, 232. [Google Scholar] [CrossRef]
  19. Abdeljawad, T.; Al-Mdallal, Q.M.; Hajji, M.A. Arbitrary order fractional difference operators with discrete exponential kernels and applications. Discret. Dyn. Nat. Soc. 2017, 2017, 4149320. [Google Scholar] [CrossRef]
  20. Mohammed, P.O.; Abdeljawad, T.; Hamasalh, F.K. On Riemann–Liouville and Caputo fractional forward difference monotonicity analysis. Mathematics 2021, 9, 1303. [Google Scholar] [CrossRef]
  21. Nagai, A. Discrete Mittag–Leffler function and its applications. Publ. Res. Inst. Math. Sci. Kyoto Univ. 2003, 1302, 1–20. [Google Scholar]
  22. Wu, G.C.; Baleanu, D.; Zeng, S.D.; Luo, W.H. Mittag–Leffler function for discrete fractional modelling. J. King Saud Univ. Sci. 2016, 28, 99–102. [Google Scholar] [CrossRef]
  23. Atici, F.M.; Chang, S.; Jonnalagadda, J.M. Mittag–Leffler functions in discrete time. Fractal Fract. 2023, 7, 254. [Google Scholar] [CrossRef]
  24. Mohammed, P.O.; Abdeljawad, T.; Hamasalh, F.K. Discrete Prabhakar fractional difference and sum operators. Chaos Soliton. Fract. 2021, 150, 111182. [Google Scholar] [CrossRef]
  25. Saenko, V.V. The calculation of the Mittag–Leffler function. Int. J. Comput. Math. 2022, 99, 1367–1394. [Google Scholar] [CrossRef]
  26. Abdeljawad, T. Different type kernel h–fractional differences and their fractional h–sums. Chaos Solit. Fract. 2018, 116, 146–156. [Google Scholar] [CrossRef]
  27. Mohammed, P.O.; Abdeljawad, T. Discrete generalized fractional operators defined using h-discrete Mittag–Leffler kernels and applications to AB fractional difference systems. Math. Meth. Appl. Sci. 2020, 46, 7688–7713. [Google Scholar] [CrossRef]
  28. Fereira, R.; Torres, D. Fractional h–difference equations arising from the calculus of variations. Appl. Anal. Discret. Math. 2011, 5, 110–121. [Google Scholar] [CrossRef]
  29. Abdeljawad, T. Dual identities in fractional difference calculus within Riemann. Adv. Differ. Equ. 2013, 2013, 36. [Google Scholar] [CrossRef]
  30. Abdeljawad, T. On delta and nabla caputo fractional differences and dual identities. Discr. Dynam. Nat. Soc. 2013, 2013, 12. [Google Scholar] [CrossRef]
  31. Mohammed, P.O. On the delta Mittag–Leffler functions and its application in monotonic analysis. J. Comput. Appl. Math. 2025. accepted. [Google Scholar]
  32. Mohammed, P.O.; Kürt, C.; Abdeljawad, T. Bivariate discrete Mittag–Leffler functions with associated discrete fractional operators. Chaos Solit. Fractals 2022, 165, 112848. [Google Scholar] [CrossRef]
  33. Zhang, X.; Feng, Y.; Luo, Z.; Liu, J. A spatial sixth-order numerical scheme for solving fractional partial differential equation. Appl. Math. Lett. 2025, 159, 109265. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Vivas-Cortez, M.; Mohammed, P.O.; Guirao, J.L.G.; Yousif, M.A.; Ibrahim, I.S.; Chorfi, N. Improved Fractional Differences with Kernels of Delta Mittag–Leffler and Exponential Functions. Symmetry 2024, 16, 1562. https://doi.org/10.3390/sym16121562

AMA Style

Vivas-Cortez M, Mohammed PO, Guirao JLG, Yousif MA, Ibrahim IS, Chorfi N. Improved Fractional Differences with Kernels of Delta Mittag–Leffler and Exponential Functions. Symmetry. 2024; 16(12):1562. https://doi.org/10.3390/sym16121562

Chicago/Turabian Style

Vivas-Cortez, Miguel, Pshtiwan Othman Mohammed, Juan L. G. Guirao, Majeed A. Yousif, Ibrahim S. Ibrahim, and Nejmeddine Chorfi. 2024. "Improved Fractional Differences with Kernels of Delta Mittag–Leffler and Exponential Functions" Symmetry 16, no. 12: 1562. https://doi.org/10.3390/sym16121562

APA Style

Vivas-Cortez, M., Mohammed, P. O., Guirao, J. L. G., Yousif, M. A., Ibrahim, I. S., & Chorfi, N. (2024). Improved Fractional Differences with Kernels of Delta Mittag–Leffler and Exponential Functions. Symmetry, 16(12), 1562. https://doi.org/10.3390/sym16121562

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