Modified Fractional Difference Operators Defined Using Mittag-Leffler Kernels

Abstract

The discrete fractional operators of Riemann–Liouville and Liouville–Caputo are omnipresent due to the singularity of the kernels. Therefore, convexity analysis of discrete fractional differences of these types plays a vital role in maintaining the safe operation of kernels and symmetry of discrete delta and nabla distribution. In their discrete version, the generalized or modified forms of various operators of fractional calculus are becoming increasingly important from the viewpoints of both pure and applied mathematical sciences. In this paper, we present the discrete version of the recently modified fractional calculus operator with the Mittag-Leffler-type kernel. Here, in this article, the expressions of both the discrete nabla derivative and its counterpart nabla integral are obtained. Some applications and illustrative examples are given to support the theoretical results.

1. Introduction

Fractional calculus is a 327-year-old interdisciplinary field [1,2] in which the integral and derivatives of a fractional real or complex order have been investigated (see [3,4,5,6]) as well as their applications (see [7,8,9,10]). It is well-known that there are many definitions of fractional calculus operators, mainly due to the fact that there is no single extension of meaning in this area. We can also mention that the Liouville–Caputo fractional derivative has appeared in several investigations in the history of fractional calculus (see, for example, the works of Abel [11], Liouville [12], Caputo [13], and Dzherbashian and Nersesian [14]). Based on the works of Boltzmann [15] and continuing with the contemporary works, we conclude that finding some generalized (see, for example, [16,17,18]) or modified versions [19] of some existing operators [20] is an interesting issue in the field of fractional calculus. In passing, we recall Wright’s general function ${\mathfrak{E}}_{\alpha ,\beta }(\varphi ;\mathrm{z})$, which occurred in his study of the asymptotic behavior of a certain Taylor-Maclaurin series, which obviously provides a remarkably deep generalization of the Mittag-Leffler function $E_{\alpha ,\beta }\left(\mathrm{z}\right)$ of two parameters $\alpha$ and $\beta$ (see [21], p. 424):

$${\mathfrak{E}}_{\alpha ,\beta }(\varphi ;\mathrm{z}):=\sum _{\ell =0}^{\infty }\frac{\varphi (\ell )}{\Gamma (\alpha \ell +\beta )}{\mathrm{z}}^{\ell }\left(\alpha ,\beta \in \mathbb{C};\mathrm{Re}\left(\alpha \right)>0\right),$$

(1)

where $\varphi \left(\tau \right)$ is a suitably restricted function of $\tau$. Recent works [22,23] (see also [24,25]) provide other historical and important backgrounds in detail regarding an interesting unification of Equation (1) and several multiparameter extensions of many functions happening in analytic number theory. It is defined by

$${\epsilon }_{\alpha ,\beta }(\phi ;\mathrm{z},s,a):=\sum _{\ell =0}^{\infty }\frac{\phi (\ell )}{{(\ell +a)}^s\Gamma (\alpha \ell +\beta )}{\mathrm{z}}^{\ell },$$

(2)

for $\beta ,\alpha \in \mathbb{C}$ with $\mathrm{Re}\left(\alpha \right)>0$, where for a suitably restricted function $\phi$, the parameters $\alpha ,$ $\beta ,$ s, and a are appropriately constrained. In its very specialized case where

$$\phi (\ell )=\Gamma (\alpha \ell +\beta )(\ell \in {\mathbb{N}}_0),$$

the Srivastava function ${\epsilon }_{\alpha ,\beta }(\phi ;\mathrm{z},s,a)$ reduces immediately to the Hurwitz–Lerch zeta function $\Phi (\mathrm{z},s,a)$. More importantly, the Srivastava function ${\epsilon }_{\alpha ,\beta }(\phi ;\mathrm{z},s,a)$, defined by Equation (2), has already been successfully used as the kernel of some general families of Riemann–Liouville-type fractional integrals and fractional derivatives (for details, see [22,23]). Appropriately defined operators of discrete fractional calculus, which are based upon the general operators of fractional calculus with the Srivastava function ${\epsilon }_{\alpha ,\beta }(\phi ;\mathrm{z},s,a)$ in the kernel, can possibly lead to some interesting future investigations.

The discrete fractional calculus [26,27,28,29] started to be an interesting part of fractional calculus [30,31,32,33] due to its multiple important applications in solving the complex dynamics of several complicated systems [34,35,36] arising from several fields of science and engineering [37,38,39]. Aside from that, the possibility of applying the discrete fractional calculus to improve some methods and techniques from artificial intelligence (AI) makes this branch of fractional calculus of great interest to readers. In addition, finding the discrete counterpart of a modified or generalised operator fractional operator is always an interesting topic, mainly because the discrete version of the continuous non-local operators possesses qualitatively different properties. In particular, the fractional continuous and discrete operators involving the Mittag-Leffler kernels (see [17] for the continuous version and [30,40,41,42] for the discrete version) present their own interest, due to the fact that this special function is considered to be the queen of fractional calculus.

Motivated by the aforementioned works, in the current study, we establish the discrete version of the recently modified fractional calculus operator defined using the Mittag-Leffler kernel in [19]. Furthermore, we obtain its series representation formula and some related examples.

The organization of this manuscript is as follows. Section 2 deals with the basic definitions and properties of the discrete fractional calculus used in this manuscript. Section 2.1 is dedicated to discrete Mittag-Leffler functions, and Section 2.2 is for the discrete Laplace transform, including some of its properties. Section 3 deals with recalling the discrete Atangana–Baleanu derivative of the Liouville–Caputo-type fractional operators and the new finding of discrete modifications of the discrete Atangana–Baleanu derivative of the Liouville–Caputo-type fractional operators. Section 4 contains two illustrative examples, and the last section is devoted to our conclusions.

2. Preliminary Tools

This section recalls some basic concepts of discrete fractional calculus, such as discrete Mittag-Leffler functions and discrete Laplace transformations on the time set ${\mathbb{N}}_a:=\{a,a+1,a+2,\dots \}$ for $a\in \mathbb{R}$.

2.1. Riemann–Liouville Fractional Sums and Mittag-Leffler Functions

Definition 1

(see [30,40,43]). For $\kappa \in {\mathbb{N}}_1$, the κ-rising factorial function can be expressed by

$${\mathrm{z}}^{\overline{\kappa }}=\prod _{\kappa =0}^{\kappa -1}(\mathrm{z}+\kappa ),{\mathrm{z}}^{\overline{0}}=1.$$

In general, it can be expressed as follows:

$${\mathrm{z}}^{\overline{\alpha }}=\frac{\Gamma \left(\mathrm{z}+\alpha \right)}{\Gamma \left(\mathrm{z}\right)},$$

(3)

for $\mathrm{z},\alpha \in \mathbb{R}$ such that neither $\mathrm{z}$ nor $\mathrm{z}+\alpha$ is a pole of the Gamma function. One of the major properties of this function is given by

$$\nabla \left({\mathrm{z}}^{\overline{\alpha }}\right)=\alpha {\mathrm{z}}^{\overline{\alpha -1}}.$$

Definition 2

(see [30,40,43]). For $f:{\mathbb{N}}_a⟶\mathbb{R}$ and $\alpha >0$, the nabla left-sided Riemann–Liouville fractional sum of the order α can be expressed as follows:

$$\left({}_a\nabla ^{-\alpha }f\right)\left(\mathrm{z}\right)=\frac{1}{\Gamma \left(\alpha \right)}\sum _{r=a+1}^{\mathrm{z}}{(\mathrm{z}+1-r)}^{\overline{\alpha -1}}f\left(r\right)(\forall \mathrm{z}\in {\mathbb{N}}_{a+1}),$$

(4)

In addition, for $f:{}_b\mathbb{N}:=\{b,b-1,b-2,\dots \}⟶\mathbb{R}$, the right-sided one can be expressed as follows:

$$\left(\nabla _b^{-\alpha }f\right)\left(\mathrm{z}\right)=\frac{1}{\Gamma \left(\alpha \right)}\sum _{r=\mathrm{z}}^{b-1}{(r-\rho \left(\mathrm{z}\right))}^{\overline{\alpha -1}}f\left(r\right)(\forall \mathrm{z}\in {}_{b-1}\mathbb{N}).$$

(5)

Lemma 1

(see [30,40,43]). Suppose that $a,b,\beta \in \mathbb{R}$, and $\alpha \in {\mathbb{R}}^{+}$. Then, one can have

$$\begin{array}{c}\hfill {}_a\nabla ^{-\alpha }\frac{{(\mathrm{z}-a)}^{\overline{\beta }}}{\Gamma \left(\beta +1\right)}=\frac{1}{\Gamma \left(\beta +\alpha +1\right)}{(\mathrm{z}-a)}^{\overline{\beta +\alpha }},\\ \hfill \nabla _b^{-\alpha }\frac{{(b-\mathrm{z})}^{\overline{\beta }}}{\Gamma \left(\beta +1\right)}=\frac{1}{\Gamma \left(\beta +\alpha +1\right)}{(b-\mathrm{z})}^{\overline{\beta +\alpha }}.\end{array}$$

Definition 3

(see [43]). For $\alpha ,\beta ,\gamma ,\mathrm{z}\in \mathbb{C}$ with $\mathrm{Re}\left(\alpha \right)>0$, the discrete Mittag-Leffler function of a generalized form can be expressed as follows:

$${\mathrm{E}}_{\overline{\alpha ,\beta }}^{\gamma }(\lambda ,\mathrm{z}):=\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa }\frac{{\mathrm{z}}^{\overline{\kappa \alpha +\beta -1}}{\left(\gamma \right)}_{\kappa }}{\Gamma (\kappa \alpha +\beta )}(\forall \lambda \in \mathbb{R}suchthat|\lambda |<1),$$

(6)

where ${\left(\gamma \right)}_{\kappa }=\frac{\Gamma (\gamma +\kappa )}{\Gamma \left(\kappa \right)}$. Particularly, the discrete Mittag-Leffler function of two parameters can be deduced when $\gamma =1$ as follows:

$${\mathrm{E}}_{\overline{\alpha ,\beta }}(\lambda ,\mathrm{z})={\mathrm{E}}_{\overline{\alpha ,\beta }}^1(\lambda ,\mathrm{z}):=\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa }\frac{{\mathrm{z}}^{\overline{\kappa \alpha +\beta -1}}}{\Gamma (\kappa \alpha +\beta )},$$

(7)

for $\lambda \in \mathbb{R}$, satisfying $\left|\lambda \right|<1$.

Furthermore, the discrete Mittag-Leffler function of one parameter can be deduced when $\beta =\gamma =1$ as follows:

$${\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z})={\mathrm{E}}_{\overline{\alpha ,\beta }}^1(\lambda ,\mathrm{z}):=\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa }\frac{{\mathrm{z}}^{\overline{\kappa \alpha }}}{\Gamma (\kappa \alpha +1)},$$

(8)

where λ is as explained above.

Lemma 2

(see [43]). Assume that $\alpha >0,\beta >-1,\gamma ,\mathrm{z}\in \mathbb{C}$. Then, it is asserted that

$${}_a{}^{RL}\nabla ^{-\nu }{\mathrm{E}}_{\overline{\alpha ,\beta }}^{\gamma }(\lambda ,\mathrm{z})={\mathrm{E}}_{\overline{\alpha ,\beta +\nu }}^{\gamma }(\lambda ,\mathrm{z}),$$

for $\lambda \in \mathbb{R}$, satisfying $\left|\lambda \right|<1$.

2.2. Discrete Laplace Transformation

Definition 4

(see [30]). Let f and g be defined on ${\mathbb{N}}_a$. Then, the discrete Laplace transformation can be expressed as follows:

$$\left({\mathfrak{L}}_{a}f\left(\mathrm{z}\right)\right)\left(\mathrm{S}\right)=\sum _{t=a+1}^{\infty }\frac{f\left(\mathrm{z}\right)}{{(1-\mathrm{S})}^{a+1-\mathrm{z}}}.$$

(9)

Furthermore, the discrete convolution of f and g can be expressed as follows:

$$\left(f\ast g\right)\left(\mathrm{z}\right)=\sum _{\kappa =a+1}^{\mathrm{z}}f\left(\kappa \right)g\left(\mathrm{z}-\rho \left(\kappa \right)+a\right).$$

(10)

Lemma 3

(see [30]). For any $\alpha \in \mathbb{R}\backslash \{\dots ,-2,-1,0\},$ it is asserted that

$$\left({\mathfrak{L}}_a{(\mathrm{z}-a)}^{\overline{\alpha -1}}\right)\left(\mathrm{S}\right)=\frac{\Gamma \left(\alpha \right)}{{\mathrm{S}}^{\alpha }},$$

provided that $|1-\mathrm{S}|<1$.

Lemma 4

(see [30]). Let f and g be defined on ${\mathbb{N}}_a$. Then, it is asserted that

$$\left({\mathfrak{L}}_a(f\ast g)\left(\mathrm{z}\right)\right)\left(\mathrm{S}\right)=\left({\mathfrak{L}}_{a}f\right)\left(\mathrm{S}\right)\left({\mathfrak{L}}_{a}g\right)\left(\mathrm{S}\right).$$

Lemma 5

(see [30,40,43]). Let f be defined on ${\mathbb{N}}_a$. Then, it is asserted that

$$\left({\mathfrak{L}}_a(\nabla f)\left(\mathrm{z}\right)\right)\left(\mathrm{S}\right)=\mathrm{S}\left({\mathfrak{L}}_{a}f\right)\left(\mathrm{S}\right)-f\left(a\right).$$

More generally, the following result holds true:

$$\left({\mathfrak{L}}_a\left({\nabla }^{n}f\right)\left(\mathrm{z}\right)\right)\left(\mathrm{S}\right)={\mathrm{S}}^n\left({\mathfrak{L}}_{a}f\right)\left(\mathrm{S}\right)-\sum _{ı=0}^{n-1}{\mathrm{S}}^{n-1-ı}{\nabla }^{ı}f(a+1).$$

Lemma 6

(see [30]). Let α be any real numbers. Then, it is asserted that

$$\left({\mathfrak{L}}_a{}_a{}^{RL}\nabla ^{-\alpha }f\left(\mathrm{z}\right)\right)\left(\mathrm{S}\right)=\frac{1}{{\mathrm{S}}^{\alpha }}\left({\mathfrak{L}}_{a}f\right)\left(\mathrm{S}\right).$$

Lemma 7

(see [43]). For $\alpha ,\beta ,\lambda ,\mathrm{S}$ in $\mathbb{C}$ with $\mathrm{Re}\left(\beta \right)>0,$ and if $\left|\lambda {\mathrm{S}}^{-\alpha }\right|<1$ with $\mathrm{Re}\left(\mathrm{S}\right)>0,$ it is asserted that

$$\left({\mathfrak{L}}_a{\mathrm{E}}_{\overline{\alpha ,\beta }}(\lambda ,\mathrm{z}-a)\right)\left(\mathrm{S}\right)={\mathrm{S}}^{-\beta }{\left[1-\lambda {\mathrm{S}}^{-\alpha }\right]}^{-1}$$

(11)

and, in particular, that

$$\left({\mathfrak{L}}_a{\mathrm{E}}_{\overline{\alpha ,\alpha }}(\lambda ,\mathrm{z}-a)\right)\left(\mathrm{S}\right)=\frac{1}{{\mathrm{S}}^{\alpha }-\lambda }.$$

(12)

3. Discrete Atangana–Baleanu and the Modified Atangana–Baleanu of the Liouville–Caputo Fractional Differences

We start this section by briefly recalling the Atangana–Baleanu derivative of the Liouville–Caputo-type fractional difference operators, and then we introduce their modified versions, which are discrete fractional analogues of the continuous case in [19]. For the more salient details on these subjects, see [30,40,43].

Definition 5

(see [30]). For $\lambda =-\frac{\alpha }{1-\alpha }$ and $0<\alpha <1/2$, the left discrete generalized Atangana–Baleanu of the Liouville–Caputo-type fractional difference is given by

$$\left({}_a{}^{ABC}\nabla ^{\alpha }f\right)\left(\mathrm{z}\right)=\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\sum _{\mathrm{s}=a+1}^{\mathrm{z}}{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)\left(\nabla f\right)\left(\mathrm{s}\right)(\forall \mathrm{z}\in {\mathbb{N}}_a)$$

(13)

and the right one is given by

$$\left({}^{ABC}\nabla _b^{\alpha }f\right)\left(\mathrm{z}\right)=\frac{-\mathsf{A}\left(\alpha \right)}{1-\alpha }\sum _{\mathrm{s}=\mathrm{z}}^{b-1}{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{s}-\mathrm{z}+1)\left(\Delta f\right)\left(\mathrm{s}\right)(\forall \mathrm{z}\in {}_b\mathbb{N}).$$

(14)

The corresponding discrete Atangana–Baleanu fractional sum is given by

$$\left({}_a{}^{AB}\nabla ^{-\alpha }f\right)\left(\mathrm{z}\right)=\frac{1-\alpha }{\mathsf{A}\left(\alpha \right)}f\left(\mathrm{z}\right)+\frac{\alpha }{\mathsf{A}\left(\alpha \right)}\left({}_a\nabla ^{-\alpha }f\right)\left(\mathrm{z}\right)(\forall \mathrm{z}\in {\mathbb{N}}_a),$$

(15)

where $\mathsf{A}\left(\alpha \right)>0$ such that $\mathsf{A}\left(0\right)=\mathsf{A}\left(1\right)=1$.

It is worth mentioning that the above definition is the discrete analogue of the Atangana–Baleanu fractional operators in [20]. Now, we will proceed to obtain the modified version of Definition 5.

Considering the definition in Equation (13), we see that

$$\begin{array}{c}\left({}_a{}^{ABC}\nabla ^{\alpha }f\right)\left(\mathrm{z}\right)\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left\{\sum _{\mathrm{s}=a+1}^{\mathrm{z}}{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)f\left(\mathrm{s}\right)-\sum _{\mathrm{s}=a+1}^{\mathrm{z}}{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)f(\mathrm{s}-1)\right\}\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left\{\sum _{\mathrm{s}=a+1}^{\mathrm{z}}{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)f\left(\mathrm{s}\right)-\sum _{\mathrm{s}=a}^{\mathrm{z}-1}{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s})f\left(\mathrm{s}\right)\right\}\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left\{f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)f\left(a\right)+\sum _{\mathrm{s}=a+1}^{\mathrm{z}}\left\{{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s})\right\}f\left(\mathrm{s}\right)\right\}\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left\{f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)f\left(a\right)+\sum _{\mathrm{s}=a+1}^{\mathrm{z}}\nabla {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)f\left(\mathrm{s}\right)\right\}\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left\{f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)f\left(a\right)+\lambda \sum _{\mathrm{s}=a+1}^{\mathrm{z}}{\mathrm{E}}_{\overline{\alpha ,\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)f\left(\mathrm{s}\right)\right\},\hfill \end{array}$$

where we have used (see [31])

$${\mathrm{E}}_{\overline{\alpha }}(\lambda ,1)=1\mathrm{and}\nabla {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z})=\lambda {\mathrm{E}}_{\overline{\alpha ,\alpha }}(\lambda ,\mathrm{z}).$$

Therefore, we can define the left discrete modified Atangana–Baleanu of the Liouville–Caputo-type fractional difference for $\alpha \in \left(0,\frac{1}{2}\right)$ as follows:

Definition 6.

For $\lambda =-\frac{\alpha }{1-\alpha }$ and $0<\alpha <1/2$, the left discrete modified Atangana–Baleanu of the Liouville–Caputo-type fractional difference is given by

$$\begin{array}{c}\left({}_a{}^{MABC}\nabla ^{\alpha }f\right)\left(\mathrm{z}\right)\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)f\left(a\right)+\lambda \sum _{\mathrm{s}=a+1}^{\mathrm{z}}{\mathrm{E}}_{\overline{\alpha ,\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)f\left(\mathrm{s}\right)\right]\hfill \end{array}$$

(16)

for each $\mathrm{z}$ in ${\mathbb{N}}_a$. Furthermore, by applying the action of the $\mathcal{Q}$ operator to Equation (16) (for further information on this action, see [32,37]), we can deduce the right one as follows:

$$\begin{array}{c}\left({}^{MABC}\nabla _b^{\alpha }f\right)\left(\mathrm{z}\right)\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,b-\mathrm{z})f\left(b\right)+\lambda \sum _{\mathrm{s}=\mathrm{z}}^{b-1}{\mathrm{E}}_{\overline{\alpha ,\alpha }}(\lambda ,\mathrm{s}-\mathrm{z}+1)f\left(\mathrm{s}\right)\right],\hfill \end{array}$$

(17)

for each $\mathrm{z}$ in ${}_b\mathbb{N}$.

The above can be generalized by the same technique used for $0<\alpha <1$, and we can deduce the following:

Definition 7.

For ${\lambda }_{\ell }=-\frac{\alpha -\ell }{\ell +1-\alpha }$ with $\ell <\alpha <\ell +\frac{1}{2}$ and $\ell \in {\mathbb{N}}_0$, the left discrete modified Atangana–Baleanu of the Liouville–Caputo-type fractional difference of a higher order can be expressed as follows:

$$\begin{array}{cc}\hfill \left({}_a{}^{MABC}\nabla ^{\alpha }f\right)\left(\mathrm{z}\right)& =\left({}_a{}^{ABC}\nabla ^{\alpha -\ell }{\nabla }^{\ell }f\right)\left(\mathrm{z}\right)\hfill \\ \hfill & =\frac{\mathsf{A}(\alpha -\ell )}{\ell +1-\alpha }[\left({\nabla }^{\ell }f\right)\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha -\ell }}({\lambda }_{\ell },\mathrm{z}-a)\left({\nabla }^{\ell }f\right)\left(a\right)\hfill \\ \hfill & +{\lambda }_{\ell }\sum _{\mathrm{s}=a+1}^{\mathrm{z}}{\mathrm{E}}_{\overline{\alpha -\ell ,\alpha -\ell }}({\lambda }_{\ell },\mathrm{z}-\mathrm{s}+1)\left({\nabla }^{\ell }f\right)\left(\mathrm{s}\right)]\hfill \end{array}$$

(18)

for each $\mathrm{z}$ in ${\mathbb{N}}_a$. Additionally, the right one can be expressed as follows:

$$\begin{array}{cc}\hfill \left({}^{MABC}\nabla _b^{\alpha }f\right)\left(\mathrm{z}\right)& =\frac{\mathsf{A}(\alpha -\ell )}{\ell +1-\alpha }[\left({\nabla }^{\ell }f\right)\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha -\ell }}({\lambda }_{\ell },b-\mathrm{z})\left({\nabla }^{\ell }f\right)\left(b\right)\hfill \\ \hfill & +{\lambda }_{\ell }\sum _{\mathrm{s}=\mathrm{z}}^{b-1}{\mathrm{E}}_{\overline{\alpha -\ell ,\alpha -\ell }}({\lambda }_{\ell },\mathrm{s}-\mathrm{z}+1)\left({\nabla }^{\ell }f\right)\left(\mathrm{s}\right)],\hfill \end{array}$$

(19)

for each $\mathrm{z}$ in ${}_b\mathbb{N}$.

Throughout the remainder of this article, we consider the left discrete modified Atangana–Baleanu of the Liouville–Caputo-type fractional differences. Interested readers can use the $\mathcal{Q}$ operator action on the left one’s results, and they will be able to obtain the corresponding results for the right one.

Theorem 1.

For $\alpha \in \left(0,\frac{1}{2}\right),$ each of the following results holds true:

$$\left({}_a{}^{AB}\nabla ^{-\alpha }{}_a{}^{MABC}\nabla ^{\alpha }f\right)\left(\mathrm{z}\right)=f\left(\mathrm{z}\right)-f\left(a\right)$$

(20)

and

$$\left({}_a{}^{MABC}\nabla ^{\alpha }{}_a{}^{AB}\nabla ^{-\alpha }f\right)\left(\mathrm{z}\right)=f\left(\mathrm{z}\right)-f\left(a\right){\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a),$$

(21)

for $\mathrm{z}$ in ${\mathbb{N}}_{a+1}$.

Proof. 

Denote the following:

$$\left({}_a{}^{MABC}\nabla ^{\alpha }f\right)\left(\mathrm{z}\right):=u\left(\mathrm{z}\right),\mathrm{for}\mathrm{all}\mathrm{z}\in {\mathbb{N}}_{a+1}.$$

Take the Laplace transform ${\mathfrak{L}}_a$ on both sides to obtain

$$\begin{array}{cc}\hfill \left({\mathfrak{L}}_a{}_a{}^{MABC}\nabla ^{\alpha }f\right)\left(\mathrm{S}\right)& \stackrel{by}{\underset{\mathrm{Lemma}7}{=}}\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[F\left(\mathrm{S}\right)-\frac{{\mathrm{S}}^{\alpha -1}{(1-\mathrm{S})}^a}{{\mathrm{S}}^{\alpha }-\lambda }f\left(a\right)+\frac{\lambda }{{\mathrm{S}}^{\alpha }-\lambda }F\left(\mathrm{S}\right)\right]\hfill \\ \hfill & =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[\frac{{\mathrm{S}}^{\alpha }}{{\mathrm{S}}^{\alpha }-\lambda }F\left(\mathrm{S}\right)-\frac{{\mathrm{S}}^{\alpha -1}{(1-\mathrm{S})}^a}{{\mathrm{S}}^{\alpha }-\lambda }f\left(a\right)\right]=U\left(\mathrm{S}\right),\hfill \end{array}$$

(22)

where $F\left(\mathrm{S}\right):=\left({\mathfrak{L}}_{a}f\right)\left(\mathrm{S}\right)$ and $U\left(\mathrm{S}\right):=\left({\mathfrak{L}}_{a}u\right)\left(\mathrm{S}\right)$. By solving for $F\left(\mathrm{S}\right)$, we find that

$$\begin{array}{cc}\hfill F\left(\mathrm{S}\right)& =\frac{1-\alpha }{\mathsf{A}\left(\alpha \right)}U\left(\mathrm{S}\right)+\frac{\alpha }{\mathsf{A}\left(\alpha \right)}\frac{1}{{\mathrm{S}}^{\alpha }}U\left(\mathrm{S}\right)+\frac{{(1-\mathrm{S})}^a}{\mathrm{S}}f\left(a\right)\hfill \\ \hfill & =({\mathfrak{L}}_a\left({}_a{}^{AB}\nabla ^{-\alpha }f)\left(\mathrm{z}\right)\right)\left(\mathrm{S}\right)+\frac{{(1-\mathrm{S})}^a}{\mathrm{S}}f\left(a\right),\hfill \end{array}$$

(23)

where we used the following (see [30], Equation (32)):

$$({\mathfrak{L}}_a\left({}_a{}^{AB}\nabla ^{-\alpha }f)\left(\mathrm{z}\right)\right)\left(\mathrm{S}\right)=\frac{1-\alpha }{\mathsf{A}\left(\alpha \right)}\left({\mathfrak{L}}_{a}f\right)\left(\mathrm{S}\right)+\frac{\alpha }{\mathsf{A}\left(\alpha \right)}\frac{1}{{\mathrm{S}}^{\alpha }}\left({\mathfrak{L}}_{a}f\right)\left(\mathrm{S}\right).$$

(24)

By taking the inverse Laplace transform ${\mathfrak{L}}_a^{-1}$ on both sides of (23), we obtain

$$f\left(\mathrm{z}\right)=\left({}_a{}^{AB}\nabla ^{-\alpha }u\right)\left(\mathrm{z}\right)+f\left(a\right),$$

which rearranges to the desired result of Equation (20).

To prove the second part, we set

$$\left({}_a{}^{AB}\nabla ^{-\alpha }f\right)\left(\mathrm{z}\right):=g\left(\mathrm{z}\right),\mathrm{for}\mathrm{all}\mathrm{z}\in {\mathbb{N}}_{a+1}.$$

Then, by taking the Laplace transform ${\mathfrak{L}}_a$ on $\left({}_a{}^{MABC}\nabla ^{\alpha }g\right)\left(\mathrm{z}\right)$, we obtain

$$\begin{array}{cc}\hfill \left({\mathfrak{L}}_a\left({}_a{}^{MABC}\nabla ^{\alpha }g\right)\left(\mathrm{z}\right)\right)\left(\mathrm{S}\right)& \stackrel{by}{\underset{\left(22\right)}{=}}\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[\frac{{\mathrm{S}}^{\alpha }}{{\mathrm{S}}^{\alpha }-\lambda }G\left(\mathrm{S}\right)-\frac{{\mathrm{S}}^{\alpha -1}{(1-\mathrm{S})}^a}{{\mathrm{S}}^{\alpha }-\lambda }g\left(a\right)\right]\hfill \\ \hfill & =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }[\frac{{\mathrm{S}}^{\alpha }}{{\mathrm{S}}^{\alpha }-\lambda }\left(\frac{1-\alpha }{\mathsf{A}\left(\alpha \right)}F\left(\mathrm{S}\right)+\frac{\alpha }{\mathsf{A}\left(\alpha \right)}\frac{1}{{\mathrm{S}}^{\alpha }}F\left(\mathrm{S}\right)\right)\hfill \\ \hfill & -\frac{{\mathrm{S}}^{\alpha -1}{(1-\mathrm{S})}^a}{{\mathrm{S}}^{\alpha }-\lambda }\frac{1-\alpha }{\mathsf{A}\left(\alpha \right)}f\left(a\right)]\hfill \\ \hfill & =F\left(\mathrm{S}\right)-\frac{{\mathrm{S}}^{\alpha -1}{(1-\mathrm{S})}^a}{{\mathrm{S}}^{\alpha }-\lambda }f\left(a\right),\hfill \end{array}$$

(25)

where $G\left(\mathrm{S}\right):=\left({\mathfrak{L}}_{a}g\right)\left(\mathrm{S}\right)$, and we first used

$$G(\mathrm{S})\stackrel{by}{\underset{\left(24\right)}{=}}\frac{1-\alpha }{\mathsf{A}\left(\alpha \right)}F\left(\mathrm{S}\right)+\frac{\alpha }{\mathsf{A}\left(\alpha \right)}\frac{1}{{\mathrm{S}}^{\alpha }}F\left(\mathrm{S}\right),$$

Later, we used

$$g(a)=\left({}_a{}^{AB}\nabla ^{-\alpha }f\right)\left(a\right)\stackrel{by}{\underset{\left(15\right)}{=}}\frac{1-\alpha }{\mathsf{A}\left(\alpha \right)}f\left(a\right).$$

By taking the Laplace transform ${\mathfrak{L}}_a$ on both sides of (25), we obtain

$$\begin{array}{cc}\hfill \left({}_a{}^{MABC}\nabla ^{\alpha }g\right)\left(\mathrm{z}\right)& =f\left(\mathrm{z}\right)-f\left(a\right){\mathfrak{L}}_a^{-1}\left\{\frac{{\mathrm{S}}^{\alpha -1}{(1-\mathrm{S})}^a}{{\mathrm{S}}^{\alpha }-\lambda }\right\}\hfill \\ \hfill & \stackrel{by}{\underset{\left(12\right)}{=}}f\left(\mathrm{z}\right)-f\left(a\right){\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a),\hfill \end{array}$$

which ends the second part our proof. The proof of Theorem 1 is thus completed.  □

Theorem 2.

For $\alpha \in \left(0,\frac{1}{2}\right),$ the following result provides an alternative series representation of the discrete modified Atangana–Baleanu of the Liouville–Caputo-type fractional difference:

$$\left({}_a{}^{MABC}\nabla ^{\alpha }f\right)\left(\mathrm{z}\right)=\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)f\left(a\right)+\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa +1}\left({}_a\nabla ^{-(\alpha \kappa +\alpha )}f\right)\left(\mathrm{z}\right)\right],$$

(26)

for $\mathrm{z}$ in ${\mathbb{N}}_{a+1}$.

Proof. 

According to Definitions 3 and 6, we can have

$$\begin{array}{c}\left({}_a{}^{MABC}\nabla ^{\alpha }f\right)\left(\mathrm{z}\right)\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)f\left(a\right)+\lambda \sum _{\mathrm{s}=a+1}^{\mathrm{z}}{\mathrm{E}}_{\overline{\alpha ,\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)f\left(\mathrm{s}\right)\right]\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)f\left(a\right)+\lambda \sum _{\mathrm{s}=a+1}^{\mathrm{z}}\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa }\frac{{(\mathrm{z}-\mathrm{s}+1)}^{\overline{\alpha \kappa +\alpha -1}}}{\Gamma (\alpha \kappa +\alpha )}f\left(\mathrm{s}\right)\right]\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)f\left(a\right)+\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa +1}\frac{1}{\Gamma (\alpha \kappa +\alpha )}\sum _{\mathrm{s}=a+1}^{\mathrm{z}}{(\mathrm{z}-\mathrm{s}+1)}^{\overline{\alpha \kappa +\alpha -1}}f\left(\mathrm{s}\right)\right]\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[f\left(\mathrm{z}\right)-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)f\left(a\right)+\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa +1}\left({}_a\nabla ^{-(\alpha \kappa +\alpha )}f\right)\left(\mathrm{z}\right)\right],\hfill \end{array}$$

which is the result asserted by Theorem 2.  □

4. Applications

This section presents some specific illustrative examples to verify the results which we obtained in the preceding sections. In addition, each of these examples shows the applicability of the alternative discrete modified Atangana–Baleanu of the Liouville–Caputo-type fractional difference series representation.

Example 1.

Let $f\left(\mathrm{z}\right)=\ell$ be any constant function. We can see from Definition 6 that

$$\begin{array}{cc}\hfill \left({}_a{}^{MABC}\nabla ^{\alpha }1\right)\left(\mathrm{z}\right)& =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[\ell -\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)+\ell \lambda \sum _{\mathrm{s}=a+1}^{\mathrm{z}}{\mathrm{E}}_{\overline{\alpha ,\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)\right]\hfill \\ \hfill & =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[\ell -\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)+\ell \sum _{\mathrm{s}=a+1}^{\mathrm{z}}\nabla {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}+1)\right]\hfill \\ \hfill & =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[\ell -\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)+\ell {\int }_a^{\mathrm{z}}\nabla {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s})\nabla \mathrm{s}\right]\hfill \\ \hfill & =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[\ell -\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)-\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-\mathrm{s}){|}_{\mathrm{s}=a}^{\mathrm{s}=\mathrm{z}}\right]=0.\hfill \end{array}$$

On the other hand, by using the series representation theorem (Theorem 2), we have

$$\begin{array}{cc}\hfill \left({}_a{}^{MABC}\nabla ^{\alpha }1\right)\left(\mathrm{z}\right)& =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[\ell -\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)+\ell \sum _{\kappa =0}^{\infty }{\lambda }^{\kappa +1}\left({}_a\nabla ^{-(\alpha \kappa +\alpha )}1\right)\left(\mathrm{z}\right)\right]\hfill \\ \hfill & =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[\ell -\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)+\ell \sum _{\kappa =0}^{\infty }{\lambda }^{\kappa +1}\underset{byLemma1}{\underbrace{\frac{{(\mathrm{z}-a)}^{\overline{\alpha \kappa +\alpha }}}{\Gamma \left(\alpha \kappa +\alpha +1\right)}}}\right]\hfill \\ \hfill & =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[-\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)+\ell \sum _{\kappa =0}^{\infty }{\lambda }^{\kappa }\frac{{(\mathrm{z}-a)}^{\overline{\alpha \kappa }}}{\Gamma \left(\alpha \kappa +1\right)}\right],\hfill \\ \hfill & =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }[-\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)+\ell {\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a)]=0,\hfill \end{array}$$

where we used ${(\mathrm{z}-a)}^{\overline{0}}=1$.

Example 2.

Consider

$$f\left(\mathrm{z}\right)={(\mathrm{z}-a)}^{\overline{\beta }},$$

for $\beta \in \mathbb{R}$ and $\alpha \in {\mathbb{R}}^{+}$. Then, from, Equation (26), we have

$$\begin{array}{c}\left({}_a{}^{MABC}\nabla ^{\alpha }{(\mathrm{z}-a)}^{\overline{\beta }}\right)\left(\mathrm{z}\right)\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[{(\mathrm{z}-a)}^{\overline{\beta }}-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a){\left(0\right)}^{\overline{\beta }}+\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa +1}\left({}_a\nabla ^{-(\alpha \kappa +\alpha )}{(\mathrm{z}-a)}^{\overline{\beta }}\right)\left(\mathrm{z}\right)\right]\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[{(\mathrm{z}-a)}^{\overline{\beta }}-{\mathrm{E}}_{\overline{\alpha }}(\lambda ,\mathrm{z}-a){\left(0\right)}^{\overline{\beta }}+\Gamma (\beta +1)\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa +1}\frac{{(\mathrm{z}-a)}^{\overline{\alpha \kappa +\alpha +\beta }}}{\Gamma (\alpha \kappa +\alpha +\beta +1)}\right]\hfill \\ =\frac{\mathsf{A}\left(\alpha \right)}{1-\alpha }\left[{(\mathrm{z}-a)}^{\overline{\beta }}+\Gamma (\beta +1)\sum _{\kappa =1}^{\infty }{\lambda }^{\kappa }\frac{{(\mathrm{z}-a)}^{\overline{\alpha \kappa +\beta }}}{\Gamma (\alpha \kappa +\beta +1)}\right]\hfill \\ =\frac{\Gamma (\beta +1)\mathsf{A}\left(\alpha \right)}{1-\alpha }\sum _{\kappa =0}^{\infty }{\lambda }^{\kappa }\frac{{(\mathrm{z}-a)}^{\overline{\alpha \kappa +\beta }}}{\Gamma (\alpha \kappa +\beta +1)}\hfill \\ =\frac{\Gamma (\beta +1)\mathsf{A}\left(\alpha \right)}{1-\alpha }{\mathrm{E}}_{\overline{\alpha ,\beta +1}}(\lambda ,\mathrm{z}-a),\hfill \end{array}$$

where we used ${\left(0\right)}^{\overline{\beta }}=0$.

In particular, for

$$f\left(\mathrm{z}\right)=(1-\alpha )\frac{{(\mathrm{z}-a)}^{\overline{\alpha -1}}}{\mathsf{A}\left(\alpha \right)\Gamma \left(\alpha \right)},$$

it becomes

$$\left({}_a{}^{MABC}\nabla ^{\alpha }\frac{(1-\alpha ){(\mathrm{z}-a)}^{\overline{\alpha -1}}}{\Gamma \left(\alpha \right)\mathsf{A}\left(\alpha \right)}\right)\left(\mathrm{z}\right)={\mathrm{E}}_{\overline{\alpha ,\alpha }}(\lambda ,\mathrm{z}-a).$$

This has been shown in Figure 1 for $\alpha =0.1,0.2,0.3,0.4$, and $\mathrm{z}\in {\mathbb{N}}_a^{a+20}$.

5. Concluding Remarks

Continuous fractional calculus, as an extension of its meaning, has a profound impact on its discrete version counterpart. During recent years, several singular and nonsingular fractional calculus operators were deeply scrutinized from the perspective and viewpoints of both pure and applied mathematical sciences. In particular, discrete fractional calculus has huge potential for applications in treating the extraction and modeling of hidden aspects from complicated real-world problems. Typically, the results provided by discrete fractional calculus can be easily adapted to improve and extend the classical methods and techniques within the artificial intelligence field. In particular, the operator possessing the Mittag-Leffler-type kernel, in the Liouville–Caputo sense, was extended in several ways in order to bypass the standard initialization procedure. In this paper, we constructed the modified nabla version of the Atangana–Baleanu discrete operator of the Liouville–Caputo type in Definitions 6 and 7, which was derived from the original discrete Atangana–Baleanu operator as defined in Definition 5. Furthermore, we proved the corresponding commutation relations and its series representation formula in Theorem 2. Aside from that, some examples were provided to see the similarities and differences with the classical discrete counterpart. We believe that the new expressions presented in this paper can be applied successfully in the modeling of complicated systems from various branches of engineering and science. We also briefly indicated some potential directions for further research by using much more general families of the appropriately defined operators of discrete fractional calculus, which are based upon the continuous fractional calculus operators with the kernel involving the Srivastava function ${\epsilon }_{\alpha ,\beta }(\phi ;\mathrm{z},s,a)$ defined by Equation (2) (for details, see [22,23]).