
Nabla discrete fractional Mittag-Leffler (ML) functions are the key of discrete fractional calculus within nabla analysis since they extend nabla discrete exponential functions. In this article, we define two new nabla discrete ML functions depending on the Cayley-exponential function on time scales. While, the nabla discrete ML function E¯γ(λ,t) converges for |λ|<1, both of the defined discrete functions converge for more relaxed λ. The nabla discrete Laplace transforms of the newly defined functions are calculated and confirmed as well. Some illustrative graphs for the two extensions are provided.
Citation: Thabet Abdeljawad. Two discrete Mittag-Leffler extensions of the Cayley-exponential function[J]. AIMS Mathematics, 2023, 8(6): 13543-13555. doi: 10.3934/math.2023687
[1] | Zunyuan Hu, Can Li, Shimin Guo . Fast finite difference/Legendre spectral collocation approximations for a tempered time-fractional diffusion equation. AIMS Mathematics, 2024, 9(12): 34647-34673. doi: 10.3934/math.20241650 |
[2] | Sarkhel Akbar Mahmood, Pshtiwan Othman Mohammed, Dumitru Baleanu, Hassen Aydi, Yasser S. Hamed . Analysing discrete fractional operators with exponential kernel for positivity in lower boundedness. AIMS Mathematics, 2022, 7(6): 10387-10399. doi: 10.3934/math.2022579 |
[3] | Maryam Saddiqa, Ghulam Farid, Saleem Ullah, Chahn Yong Jung, Soo Hak Shim . On Bounds of fractional integral operators containing Mittag-Leffler functions for generalized exponentially convex functions. AIMS Mathematics, 2021, 6(6): 6454-6468. doi: 10.3934/math.2021379 |
[4] | Juan L. G. Guirao, Pshtiwan Othman Mohammed, Hari Mohan Srivastava, Dumitru Baleanu, Marwan S. Abualrub . Relationships between the discrete Riemann-Liouville and Liouville-Caputo fractional differences and their associated convexity results. AIMS Mathematics, 2022, 7(10): 18127-18141. doi: 10.3934/math.2022997 |
[5] | Erdal Bas, Ramazan Ozarslan . Theory of discrete fractional Sturm–Liouville equations and visual results. AIMS Mathematics, 2019, 4(3): 593-612. doi: 10.3934/math.2019.3.593 |
[6] | Pshtiwan Othman Mohammed, Donal O'Regan, Aram Bahroz Brzo, Khadijah M. Abualnaja, Dumitru Baleanu . Analysis of positivity results for discrete fractional operators by means of exponential kernels. AIMS Mathematics, 2022, 7(9): 15812-15823. doi: 10.3934/math.2022865 |
[7] | Hengxiao Qi, Muhammad Yussouf, Sajid Mehmood, Yu-Ming Chu, Ghulam Farid . Fractional integral versions of Hermite-Hadamard type inequality for generalized exponentially convexity. AIMS Mathematics, 2020, 5(6): 6030-6042. doi: 10.3934/math.2020386 |
[8] | Dong Pan, Huizhen Qu . Finite-time boundary synchronization of space-time discretized stochastic fuzzy genetic regulatory networks with time delays. AIMS Mathematics, 2025, 10(2): 2163-2190. doi: 10.3934/math.2025101 |
[9] | Erdal Bas, Ramazan Ozarslan, Resat Yilmazer . Spectral structure and solution of fractional hydrogen atom difference equations. AIMS Mathematics, 2020, 5(2): 1359-1371. doi: 10.3934/math.2020093 |
[10] | Mohammad Tamsir, Neeraj Dhiman, Deependra Nigam, Anand Chauhan . Approximation of Caputo time-fractional diffusion equation using redefined cubic exponential B-spline collocation technique. AIMS Mathematics, 2021, 6(4): 3805-3820. doi: 10.3934/math.2021226 |
Nabla discrete fractional Mittag-Leffler (ML) functions are the key of discrete fractional calculus within nabla analysis since they extend nabla discrete exponential functions. In this article, we define two new nabla discrete ML functions depending on the Cayley-exponential function on time scales. While, the nabla discrete ML function E¯γ(λ,t) converges for |λ|<1, both of the defined discrete functions converge for more relaxed λ. The nabla discrete Laplace transforms of the newly defined functions are calculated and confirmed as well. Some illustrative graphs for the two extensions are provided.
In the last decade of the last century, the time scale calculus appeared as a unification tool for the discrete, quantum and ordinary calculi [10]. Since then, it has affected the development of both the ordinary and fractional calculus. Many researchers have started to study fractional calculus benefiting from the time scales tools [3,9,20]. Exponential functions are considered to be the main set of functions to study linear dynamic equations. Due to this vitality, they occupied a considerable place in the time scale theory. Different version of exponential functions have appeared in this sequel. The author in [11] formulated a new type of exponential functions on time scales called Cayley-exponential function. Indeed, instead of assuming that the exponential function should satisfy the dynamic equation Δχ(t)=λχ(t),χ(0)=1, for example, he requested it to satisfy the dynamic equation
Δχ(t)=λ2[χ(t)+χ(σ(t))],χ(0)=1, | (1.1) |
where Δ is the forward type derivative on time scales and σ(t) is the forward jumping operator. Notice that when the time scale is R, then σ(l)=l and, hence, Eq (1.1) will be reduced to the ordinary linear equation ddtχ(t)=λχ(t),χ(0)=1 whose solution, then becomes the exponential function eλt and hence this new exponential function has a sense. The author in [11] used a bijection that relates the delta exponential function to Cayley-exponential function by replacing χ(σ(t)) by Δχ(t)+χ(t). In fractional calculus the ML functions play the role of or extended exponential functions [14,15,17]. On the other hand, researchers extended the delta and nabla discrete (on the time scale Z) exponential functions to nabla and delta discrete ML functions either in delta ot nabla senses [2,4,7,19]. For some recent developments in discrete fractional calculus and its importance in short memory, we refer to [18,20]. In this work, we define nabla discrete ML extensions to discrete Cayley-exponential functions. We obtain two types of extensions: one by using the bijection through replacing χ(t−1) by χ(t)−∇χ(t) and then allowing the nabla Caputo fractional difference to act in place of ∇χ(t). The first type, the one through the bijection, is easy to obtain by making use of the already defined nabla discrete ML function in literature mentioned above. However, the second looks harder and its compact formula consists of double summations instead of one summation. For that reason, we may call nabla discrete ML functions as first level extensions and the one with double summation a second level one. The two newly defined nabla discrete ML functions in this article have the advantage that they relax the interval of convergence of the parameter λ.
This article is organized as follows: Section 2 includes preliminaries about discrete fractional sums and differences in the nabla sense, where the nabla discrete exponential and ML functions are recalled, and provides some nabla discrete Laplace transforms tools. In Section 3, we give a discrete ML function extending the nabla exponential function of Cayley type via a bijection. Section 4 is about the Cayley ML extension without the use of the bijection. Section 5 is devoted to the discrete Laplace transforms of the obtained discrete ML extensions. Finally, in Section 6 we summarize our conclusions.
In this section, we review some basics about nabla discrete fractional calculus on a time scale N0={0,1,2,...}. For more details and more generally on a time scale N0,h={0,h,2h,3h,..},h>0, we refer the reader to [6]. The book [12] is a recommended reference as well.
In this section we set some essential terminology to proceed on.
Definition 2.1. [9] (The rising factorial)
(i) Let ς be a natural number; then the ς rising factorial of t has the form
t¯ς=ς−1∏k=0(t+k),t¯0=1. | (2.1) |
(ii) For any real number, the γ rising function becomes
t¯γ=Γ(t+γ)Γ(t),suchthatt∈C∖{...,−2,−1,0},0¯γ=0. | (2.2) |
In addition, we have
∇(t¯γ)=γt¯γ−1, | (2.3) |
hence t¯γ is increasing on N0.
The backward difference operator on Z is given by ∇χ(t)=χ(t)−χ(t−1), while the jumping operator is ρ(n)=n−1.
Consider the dynamic equation
∇χ(t)=ω(t)χ(t),t∈Na,χ(a)=1, | (2.4) |
where 1−ω(t)≠0 on Na={a,a+1,a+2,...}.
Then, the nabla exponential function ˆep(t,a) satisfies (2.4) and if ω is a positive constant with 1−ω≠0 then it has the form:
ˆeω(t,a)=(11−ω)t−a,t∈Na+1. | (2.5) |
Definition 2.2. (Discrete ML function in nabla sense) ([4,5]) For λ∈R,|λ|<1 and γ,z∈C with Re(γ)>0, the nabla discrete ML function is
E¯γ(λ,z)≜∞∑k=0λkz¯kγΓ(γk+1),|λ|<1. | (2.6) |
To know more details about the properties and the use of ML functions in fractional calculus, we refer to [14,15,17].
Definition 2.3. [6] For a function χ:Na={a,a+1,a+2…}→R, the nabla left fractional sum of order γ>0 is given by
(a∇−γχ)(t)=1Γ(γ)t∑k=a+1(t−ρ(k))¯γ−1χ(k),t∈Na+1. |
Definition 2.4. [6] (The Caputo fractional differences) Let γ∈(0,1], then, the left Caputo fractional difference of order γ starting at a is defined by
(Ca∇γχ)(t)=(a∇−(1−γ)∇χ)(t),t∈Na+1. | (2.7) |
Lemma 2.1. [6] Let γ>0,μ>0. Then,
a∇−γ(t−a)¯μ=Γ(μ+1)Γ(μ+1+γ)(t−a)¯γ+μ. | (2.8) |
The Laplace transform on time scales was studied in [10,13,21]. The Z−transform can be found in [10] as well. As a special isolated time scale, we can consider the Laplace transform for the discrete fractional calculus. Following the nabla time scale book [10], the discrete Laplace transform on Na={a,a+1,a+2,...} is given as follows.
Definition 2.5. The nabla discrete Laplace transform of a discrete function χ(t) defined on Na={a,a+1,a+2,...} has the form
Na{χ(t)}(s)=∞∑t=a+1(1−s)t−a−1χ(t). | (2.9) |
In case a=0 and we write
N{χ(t)}(s)=∞∑t=1(1−s)t−1χ(t). |
Definition 2.6. [5] Let s∈R, and χ1,χ2:Na→R be a functions. The discrete convolution of χ1 with χ2 within nabla is defined by
(χ1∗χ2)(t)=t∑k=a+1χ2(t−k+1+a)χ1(k). | (2.10) |
If a=0, then
(χ1∗χ2)(t)=t∑k=1χ2(t−k+1)χ1(k). | (2.11) |
Theorem 2.2. [6] (The discrete convolution theorem in nabla sense) For any γ∈R∖{…,−2,−1,0}, s∈R and χ1,χ2 defined on Na, we have
Na{(χ1∗χ2)(t)}(s)=Na{χ1(t)}(s)Na{χ2(t)}(s). | (2.12) |
Lemma 2.3. [5,6] Let χ be defined on Na. Then,
Na{∇χ(t)}(s)=sNa{χ(t)}(s)−χ(a). | (2.13) |
Lemma 2.4. (For a = 0 see [6]) For any γ∈R {…,−2,−1,0} and |1−s|<1, we have
Na{(t−a)¯γ−1}(s)=Γ(γ)sγ. |
Lemma 2.5. [6] For γ,λ∈C (Re(β)>0), and s∈C with Re(s)>0,|λs−γ|<1, we have
N{E¯γ(λ,t)}(s)=s−1[1−λs−γ]−1. |
The following discrete Laplace lemma is useful in what follows.
Lemma 2.6. Assume χ(t) is a real valued function defined on the set of integers Z. Then, we have
N{χ(t−k)}(s)=(1−s)k[0∑t=1−k(1−s)t−1χ(t)+X(s)],k=1,2,..., | (2.14) |
where N{χ(t)}(s)=X(s), and the empty sum leads to 0.
Proof. By Definition 2.5, we have
N{χ(t−k)}(s)=∞∑t=1(1−s)t−1χ(t−k)=(1−s)k[0∑t=1−k(1−s)t−1χ(t)+X(s)]. | (2.15) |
Of particular interest of Lemma 2.6 above are the following two cases:
N{χ(t−1)}(s)=χ(0)+(1−s)X(s), | (2.16) |
and
N{P(t−k)}(s)=(1−s)kX(s), | (2.17) |
where P(t)=t¯γ,γ>0. The latter case is due to the fact that dividing by poles of Gamma function leads to 0.
In this section, we study the bijection between Anderson etal-exponential function and the nabla Cayley-exponential function, and define a nabla discrete ML extension. Consider the equation
∇χ(t)=λχ(t)+χ(t−1)2,χ(0)=1,λ≠∓2. | (3.1) |
Using that χ(t−1)=χ(t)−∇χ(t), then (3.1) becomes
∇χ(t)=2λ2+λχ(t),χ(0)=1,λ≠∓2. | (3.2) |
Notice that 2λ2+λ≠1 implies that λ≠2.
Definition 3.1. (The nabla Cayley exponential function) The solution of the Eq (3.1), is called the nabla discrete Cayley function. It is given by
E¯1(λ,a)=ˆeω(t,a)=(11−ω)t−a,ω=2λ2+λt∈Na+1, | (3.3) |
or
E¯1(λ,a)=(2+λ2−λ)t−a,λ≠∓2. | (3.4) |
If we replace ∇χ(t) in (3.2) with the Caputo fractional difference (C∇γχ)(t) of order 0<γ≤1, and make use of Example 48 in [4], then the solution of (3.2) will be of the form:
χ(t)=E¯γ(ω,t)=∞∑l=0ωlt¯γlΓ(γl+1),|ω|<1,t∈N0,ω=2λ2+λ. | (3.5) |
Notice that the convergence condition |ω|<1 in (3.5) implies that λ∈(−23,2).
On the light of (3.5), we can define the following discrete ML extension of the Cayley-exponential (C-exponential) function obtained by making use of the bijection ω=2λ2+λ.
Definition 3.2. The nabla discrete ML function in the sense of Cayley via the bijection ω is defined by
CBE¯γ(λ,t)=E¯γ(ω,t),ω=2λ2+λ,λ∈(−23,2),and0<γ≤1, |
at least.
Remark 3.1. Notice that while the nabla discrete ML function E¯γ(λ,t) extending the discrete Anderson etal-exponential function converges for λ∈(−1,1), the ML extension CBE¯γ(λ,t) converges for λ∈(−23,2).
Below Figure 1 explains the graphs of the CBE¯γ(λ,t) function and its convergence behavior for γ=12 and different values of λ (0.5,0.7,0.9,1.2).
In general, it is not true that
0∇−γ(f(r−k))(t)=(0∇−γf(r))(t−k), | (4.1) |
where f:Z→R is a discrete function defined on the integers. However, for the power functions we have the following assertion.
Lemma 4.1. Assume γ>0 and η>−1. Then, for any k=0,1,2,... and t∈N0, we have
(0∇−γ(r−k)¯η)(t)=(0∇−γr¯η)(t−k). | (4.2) |
Proof. We follow the proof by induction on k. For k=1, we have
(0∇−γ(r−1)¯η)(t)=1Γ(γ)t∑s=1(t−s+1)¯γ−1(s−1)¯η | (4.3) |
=1Γ(γ)t−1∑s=0(t−s)¯γ−1s¯η. | (4.4) |
Now, by (4.3) we have
(0∇−γ(r−1)¯η)(t+1)=1Γ(γ)t∑s=0(t−s+1)¯γ−1s¯η | (4.5) |
=1Γ(γ)t∑s=1(t−s+1)¯γ−1s¯η | (4.6) |
=(0∇−γr¯η)(t). | (4.7) |
Hence, by replacing t by t−1 we obtain the case k=1. In fact, the assertion for any k=2,3,... can be proved by following similar to above and by using the fact that (−j)¯η=0 for any j=0,1,2,3,..., which is due to that division by a pole of the Gamma function results in 0.
Remark 4.1. ● If in Lemma 4.1, the step in the discrete time scale is any positive h, then we can have
(0∇−γh(r−kh)¯ηh)(t)=(0∇−γr¯ηh)(t−kh),k=0,1,2,.... | (4.8) |
● From above it can be seen that (4.1) is in general true for discrete functions f:Z→R satisfying f(−j)=0 for all j,0,−1,−2,−3,....
Consider the equation
C0∇γχ(t)=λχ(t)+χ(t−1)2,χ(0)=1,t=1,2,3,...,χ:Z→R,γ∈(0,1]. | (4.9) |
By applying the nabla fractional sum 0∇−γ of order γ on the left hand side of (4.9) and making use of Proposition 28 in [4], we get the solution representation:
χ(t)=χ(0)+λ20∇−γ[χ(r)+χ(r−1)](t). | (4.10) |
To proceed by the successive approximation (the existence and uniqueness theorem is straight forward), we set
χs(t)=1+λ20∇−γ[χs−1(r)+χs−1(r−1)](t),s=1,2,3,...,χ0(t)≡1. | (4.11) |
For s=1, and by the help of the identity (see for example Proposition 22 in [4])
0∇−γt¯η=Γ(η+1)Γ(η+1+γ)t¯η+γ,γ>−1, | (4.12) |
or by Lemma 2.1, with a=0, we have
χ1(t)=1+λ20∇−γ[1+1](t)=1+λt¯γΓ(γ+1). | (4.13) |
For s=2, we have
χ2(t)=1+λ20∇−γ[1+λr¯γΓ(γ+1)+1+λ(r−1)¯γΓ(γ+1)](t). |
Lemma 4.1 with k=1, then implies that
χ2(t)=1+λ2[2t¯γΓ(γ+1)+λt¯2γΓ(2γ+1)+λ(t−1)¯2γΓ(2γ+1)]. |
Similarly for s=3, Lemma 4.1 with k=1,2 implies
χ3(t)=1+λ2[2t¯γΓ(γ+1)+λt¯2γΓ(2γ+1)+λ22t¯3γΓ(3γ+1)+λ(t−1)¯2γΓ(2γ+1)+λ2(t−1)¯3γΓ(3γ+1)+λ22(t−2)¯3γΓ(3γ+1)]. |
If we proceed in this way, then we claim the following form:
χs(t)=1+λ2s∑i=1i−1∑j=012i−2Cji−1λi−1(t−j)¯iγΓ(iγ+1) if s=0,1,2,3,... | (4.14) |
where an empty sum is understood to be nil and Cji=i!j!(i−j)!. To prove that (4.14) is true we proceed inductively. Assume (4.14) is true for s=n, and we need to show that it is true for s=n+1. That is
xn+1(t)=1+λ2n+1∑i=1i−1∑j=012i−2Cji−1λi−1(t−j)¯iγΓ(iγ+1). | (4.15) |
To this end, we have
xn+1(t)=1+λ2.0∇−γ(xn(t)+xn(t−1))=1+λ2.0∇−γ(1+λ2n∑i=1i−1∑j=012i−2Cji−1λi−1(t−j)¯iγΓ(iγ+1)+1+λ2n∑i=1i−1∑j=012i−2Cji−1λi−1(t−1−j)¯iγΓ(iγ+1))=1+λ2.0∇−γ(2+λ2n∑i=1i−1∑j=012i−2Cji−1λi−1(t−j)¯iγΓ(iγ+1)+λ2n∑i=1i−1∑j=012i−2Cji−1λi−1(t−(1+j))¯iγΓ(iγ+1))=1+λ2.0∇−γ(2+λ2n∑i=1i−1∑j=012i−2Cji−1λi−1(t−j)¯iγΓ(iγ+1)+λ2n∑i=1i∑j=112i−2Cj−1i−1λi−1(t−j)¯iγΓ(iγ+1))=1+λ2.0∇−γ(2+λ2n∑i=112i−2C0i−1λi−1t¯iγΓ(iγ+1)+λ2n∑i=1i−1∑j=112i−2Cji−1λi−1(t−j)¯iγΓ(iγ+1)+λ2n∑i=1i−1∑j=112i−2Cj−1i−1λi−1(t−j)¯iγΓ(iγ+1)+λ2n∑i=112i−2Ci−1i−1λi−1(t−i)¯iγΓ(iγ+1))=1+λ2.0∇−γ(2+λ2n∑i=112i−2C0i−1λi−1t¯iγΓ(iγ+1)+λ2n∑i=1i−1∑j=112i−2(Cji−1+Cj−1i−1)λi−1(t−j)¯iγΓ(iγ+1)+λ2n∑i=112i−2Ci−1i−1λi−1(t−i)¯iγΓ(iγ+1)). |
That is,
χn+1(t)=1+λ2.0∇−γ(2+λ2n∑i=112i−2C0i−1λi−1t¯iγΓ(iγ+1)+λ2n∑i=1i−1∑j=112i−2Cjiλi−1(t−j)¯iγΓ(iγ+1)+λ2n∑i=112i−2Ci−1i−1λi−1(t−i)¯iγΓ(iγ+1))=1+λ2(2t¯γΓ(γ+1)+λ2n∑i=112i−2C0i−1λi−1t¯(i+1)γΓ((i+1)γ+1)+λ2n∑i=1i−1∑j=112i−2Cjiλi−1(t−j)¯(i+1)γΓ((i+1)γ+1)+λ2n∑i=112i−2Ci−1i−1λi−1(t−i)¯(i+1)γΓ((i+1)γ+1))=1+λ2(4λ1−1t¯γ2Γ(γ+1)+λ2n∑i=112i−2λi−1t¯(i+1)γΓ((i+1)γ+1)+λ2n∑i=1i−1∑j=112i−2Cjiλi−1(t−j)¯(i+1)γΓ((i+1)γ+1)+λ2n∑i=112i−2λi−1(t−i)¯(i+1)γΓ((i+1)γ+1))=1+λ2(λ2λ−1t¯γ2−2Γ(γ+1)+λ2n∑i=112i−2C0iλi−1t¯(i+1)γΓ((i+1)γ+1)+λ2n∑i=1i−1∑j=112i−2Cjiλi−1(t−j)¯(i+1)γΓ((i+1)γ+1)+λ2n∑i=112i−2λi−1(t−i)¯(i+1)γΓ((i+1)γ+1))=1+λ2(λ2λ0−1t¯γ20−2Γ(γ+1)+λ2n∑i=1i−1∑j=012i−2Cjiλi−1(t−j)¯(i+1)γΓ((i+1)γ+1)+λ2n∑i=112i−2λi−1(t−i)¯(i+1)γΓ((i+1)γ+1))=1+λ2(λ2n∑i=0i−1∑j=012i−2Cjiλi−1(t−j)¯(i+1)γΓ((i+1)γ+1)+λ2n∑i=012i−2Ciiλi−1(t−i)¯(i+1)γΓ((i+1)γ+1))=1+λ2(λ2n∑i=0i∑j=012i−2Cjiλi−1(t−j)¯(i+1)γΓ((i+1)γ+1))=1+λ2(λ2n∑i=0i∑j=012(i+1)−3Cjiλ(i+1)−2(t−j)¯(i+1)γΓ((i+1)γ+1))=1+λ2(λ2n+1∑i=1i−1∑j=012i−3Cji−1λi−2(t−j)¯iγΓ(iγ+1))=1+λ2(λ2n+1∑i=1i−1∑j=012i−2λ−12−1Cji−1λi−1(t−j)¯iγΓ(iγ+1))=1+λ2(λ−12−1λ2n+1∑i=1i−1∑j=012i−2Cji−1λi−1(t−j)¯iγΓ(iγ+1))=1+λ2(n+1∑i=1i−1∑j=012i−2Cji−1λi−1(t−j)¯iγΓ(iγ+1)). |
Which completes the proof.
Upon (4.14), we set the following definition.
Definition 4.1. The nabla discrete ML function in the sense of Cayley without bijection use, of order γ, is given by:
CAE¯γ(λ,t)=1+λ2∞∑i=1i−1∑j=012i−2Cji−1λi−1(t−j)¯iγΓ(iγ+1),t=0,1,2,3,.... | (4.16) |
Remark 4.2. By making use of Definition 9, the solution of (4.9) is given by
χ(t)=CAE¯γ(λ,t). |
Notice that it converges for |λ|<2.
Below Figure 2 explains the graph of the discrete function CAE¯γ(λ,t) and its convergence behavior for γ=12 and λ=0.7.
From Lemma 2.5, we recall that for |λs−γ|<1, we have
N{E¯γ(λ,t)}(s)=s−1[1−λs−γ]−1=sγ−1(sγ−λ). |
Hence, we have
N{CBE¯γ(λ,t)}(s)=N{E¯γ(ω,t)}(s)=sγ−1(sγ−ω). | (5.1) |
If we substitute ω=2λ2+λ, then for |ωs−γ|<1, we have
N{CBE¯γ(λ,t)}(s)=sγ−1(2+λ))(2+λ)sγ−2λ. | (5.2) |
To calculate the discrete Laplace transform of CAE¯γ(λ,t), we shall follow an indirect way. Namely, we shall solve the fractional difference initial value problem (4.9) by using the nabla discrete Laplace transform. Firstly, by writing (C0∇γχ)(t)=t¯−γΓ(1−γ)∗(∇χ)(t), and by the help Theorem 2.2, Lemma 2.3 and Lemma 2.4, we have
N{(C0∇γχ)(t)}(s)=sγ−1[sX(s)−1], | (5.3) |
where N{χ(t)}(s)=X(s). On the other hand, by (2.16) with χ(0)=1, we conclude that
sγ−1[sX(s)−1]=λ2[X(s)+1+(1−s)X(s)]. |
Hence,
X(s)=λ+2sγ−12sγ+λ(s−2). | (5.4) |
Finally, upon Remark 4.2, we conclude that
N{CAE¯γ(λ,t)}(s)=λ+2sγ−12sγ+λ(s−2). | (5.5) |
On the other hand, it seems complicated to calculate the nabla discrete Laplace of N{CAE¯γ(λ,t)}(s) by directly applying the N to the formula stated in Definition 4.1. Indeed, by making use of (2.17) or Lemma 2.6, we have
N{CAE¯γ(λ,t)}(s)=1s+2s∞∑i=1(λ2sγ)ii−1∑j=0Cji−1(1−s)j. | (5.6) |
Therefore, we conclude that
1s+2s∞∑i=1(λ2sγ)ii−1∑j=0Cji−1(1−s)j=λ+2sγ−12sγ+λ(s−2), | (5.7) |
where the left hand side of (5.7) is at least complicated to be evaluated. But if we set Λ(s)=[λs+2sγ4sγ+2λ(s−2)−12] and Di(s)=∑i−1j=0Cji−1(1−s)j, then under the condition that |λDi(s)1/i2sγ|<1, we can come to the conclusion that
Di(s)=[2Λ(s)sγ2+Λ(s)]i. | (5.8) |
We have defined two discrete ML functions CBE¯γ(λ,t) and CAE¯γ(λ,t) that extended the nabla discrete exponential functions in the sense Cayley. In fact, using the fractional order gives the chance to get the double series discrete ML extension CAE¯γ(λ,t), whereas the Cayley exponential function is only related to the classical one by a bijection. Both of the obtained discrete ML functions obtained will result in a discrete Cayley exponential function when γ=1.
We have noticed that, while the discrete ML function E¯γ(λ,t) converges for |λ|<1, CBE¯γ(λ,t) and CAE¯γ(λ,t) converge for λ∈(−23,2) and |λ|<2, respectively. Further, we have calculated the nabla discrete Laplace transforms of both of the defined discrete ML functions. The newly defined discrete Mittag-Leffler functions are allowed to express the solutions of linear fractional difference equations when we average the solution and its one step backward in the right hand side of the linear equation. Due to this averaging, the formula is expressed by means of double summation seeking for better representations. Hence, we believe that the second level extension CAE¯γ(λ,t) will play a central role in developing a new branch of discrete fractional calculus with better convergence criteria.
In this article, we have considered the nabla approach. Although, we can always construct dual identities between nabla and delta fractional difference operators, it is unclear to find direct dual identities that relate nabla and delta discrete exponential functions and more generally Mittag-Leffler functions. Hence, it will be of interest to consider the delta case in future works. For the classical delta discrete Mittag-Leffler case, we may refer to [22]. Also, it is of interest to consider the q−fractional case to study several theoretical aspects, such as in [1] and proceed in applications, such as in [8,16].
The author would like to thank Prince sultan University for paying the APC and for the support through the TAS research lab. The author would like to thank Prof. Hamid Boulares for the help in obtaining the compact formula (4.14).
The authors declare no conflicts of interest.
[1] |
B. Abdalla, On the oscillation of q-fractional difference equations, Adv. Differ. Equ., 2017 (2017), 254. https://doi.org/10.1186/s13662-017-1316-x doi: 10.1186/s13662-017-1316-x
![]() |
[2] |
T. Abdeljawad, On Riemann and Caputo fractional differences, Comput. Math. Appl., 62 (2011), 1602–1611. https://doi.org/10.1016/j.camwa.2011.03.036 doi: 10.1016/j.camwa.2011.03.036
![]() |
[3] |
T. Abdeljawad, F. M. Atici, On the definitions of nabla fractional differences, Abstr. Appl. Anal., 2012 (2012), 1–13. https://doi.org/10.1155/2012/406757 doi: 10.1155/2012/406757
![]() |
[4] |
T. Abdeljawad, On delta and nabla Caputo fractional differences and dual identities, Discrete Dyn. Nat. Soc., 2013 (2013), 1–12. https://doi.org/10.1155/2013/406910 doi: 10.1155/2013/406910
![]() |
[5] |
T. Abdeljawad, F. Jarad, D. Baleanu, A semigroup-like property for discrete Mittag-Leffler functions, Adv. Differ. Equ., 2012 (2012), 72. https://doi.org/10.1186/1687-1847-2012-72 doi: 10.1186/1687-1847-2012-72
![]() |
[6] |
T. Abdeljawad, Different type kernel h-fractional differences and their fractional h-sums, Chaos Solitons Fract., 116 (2018), 146–156. https://doi.org/10.1016/j.chaos.2018.09.022 doi: 10.1016/j.chaos.2018.09.022
![]() |
[7] |
N. Acar, F. M. Atıcı, Exponential functions of discrete fractional calculus, Appl. Anal. Discrete Math., 7 (2013), 343–353. https://doi.org/10.2298/AADM130828020A doi: 10.2298/AADM130828020A
![]() |
[8] |
J. Alzabut, S. R. Grace, J. M. Jonnalagadda, E. Thandapani, Bounded non-oscillatory solutions of nabla forced fractional difference equations with positive and negative terms, Qual. Theory Dyn. Syst., 22 (2023), 28. https://doi.org/10.1007/s12346-022-00729-0 doi: 10.1007/s12346-022-00729-0
![]() |
[9] | F. M. Atıcı, P. W. Eloe, Discrete fractional calculus with the nabla operator, Electr. J. Qual. Theory Differ. Equ., 2009, 1–12. |
[10] | M. Bohner, A. Peterson, Advances in dynamic equations on time scales, Birkhäuser, Boston, 2003. https://doi.org/10.1007/978-0-8176-8230-9 |
[11] |
J. L. Cieslinski, Some implications of a new approach to exponential functions on time scales, Discrete Cont. Dyn. Syst., 2011 (2011), 302–311. https://doi.org/10.3934/proc.2011.2011.302 doi: 10.3934/proc.2011.2011.302
![]() |
[12] | C. Goodrich, A. Peterson, Discrete fractional calculus, Springer, 2015. https://doi.org/10.1007/978-3-319-25562-0 |
[13] |
M. T. Holm, The Laplace transform in discrete fractional calculus, Comput. Math. Appl., 62 (2011), 1591–1601. https://doi.org/10.1016/j.camwa.2011.04.019 doi: 10.1016/j.camwa.2011.04.019
![]() |
[14] | A. Kilbas, H. M. Srivastava, J. J. Trujillo, Theory and application of fractional differential equations, North Holland Mathematics Studies, 2006. |
[15] |
A. A. Kilbas, M. Saigo, R. K. Saxena, Generalized Mittag-Leffler function and generalized fractional calculus operators, Integr. Transf. Spec. Funct., 15 (2004), 31–49. https://doi.org/10.1080/10652460310001600717 doi: 10.1080/10652460310001600717
![]() |
[16] |
N. Mlaiki, Integral-type fractional equations with a proportional Riemann-Liouville derivative, J. Math., 2021 (2021), 1–7. https://doi.org/10.1155/2021/9990439 doi: 10.1155/2021/9990439
![]() |
[17] | I. Podlubny, Fractional differential equations, Academic Press, San Diego CA, 1999. |
[18] | T. T. Song, G. C. Wu, J. L. Wei, Hadamard fractional calculus on time scales, Fractals, 30 (2022) 2250145. |
[19] |
G. C. Wu, D. Baleanu, S. Zeng, W. Luo, Mittag-Leffler function for discrete fractional modelling, J. King Saud Univ. Sci., 28 (2016), 99–102. https://doi.org/10.1016/j.jksus.2015.06.004 doi: 10.1016/j.jksus.2015.06.004
![]() |
[20] |
G. C. Wu, M. K. Luo, L. L. Huang, S. Banerjee, Short memory fractional differential equations for new neural network and memristor design, Nonlinear Dyn., 100 (2020), 3611–3623. https://doi.org/10.1007/s11071-020-05572-z doi: 10.1007/s11071-020-05572-z
![]() |
[21] |
D. Baleanu, G. C. Wu, Some further results of the Laplace transform for variable-order fractional difference equations, Fract. Calc. Appl. Anal., 22 (2019), 1641–1654. https://doi.org/10.1515/fca-2019-0084 doi: 10.1515/fca-2019-0084
![]() |
[22] | G. C. Wu, T. Abdeljawad, J. Liuc, D. Baleanud, K. T. Wu, Mittag-Leffler stability analysis of fractional discrete-time neural networks via fixed point technique, Nonlinear Anal.: Model. Contr., 24 (2019), 919–936. |