Next Article in Journal
Further Generalization and Approximation of Fractional-Order Filters and Their Inverse Functions of the Second-Order Limiting Form
Next Article in Special Issue
Mixed Convection Flow over an Elastic, Porous Surface with Viscous Dissipation: A Robust Spectral Computational Approach
Previous Article in Journal
Numerical Solutions of Variable-Coefficient Fractional-in-Space KdV Equation with the Caputo Fractional Derivative
Previous Article in Special Issue
Why Controlling the Asymptomatic Infection Is Important: A Modelling Study with Stability and Sensitivity Analysis
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Weighted Fractional Calculus: A General Class of Operators

by
Arran Fernandez
* and
Hafiz Muhammad Fahad
Department of Mathematics, Eastern Mediterranean University, 99628 Famagusta, Northern Cyprus, via Mersin 10, Turkey
*
Author to whom correspondence should be addressed.
Fractal Fract. 2022, 6(4), 208; https://doi.org/10.3390/fractalfract6040208
Submission received: 17 February 2022 / Revised: 16 March 2022 / Accepted: 4 April 2022 / Published: 7 April 2022

Abstract

:
We conduct a formal study of a particular class of fractional operators, namely weighted fractional calculus, and its extension to the more general class known as weighted fractional calculus with respect to functions. We emphasise the importance of the conjugation relationships with the classical Riemann–Liouville fractional calculus, and use them to prove many fundamental properties of these operators. As examples, we consider special cases such as tempered, Hadamard-type, and Erdélyi–Kober operators. We also define appropriate modifications of the Laplace transform and convolution operations, and solve some ordinary differential equations in the setting of these general classes of operators.

1. Introduction

Fractional calculus studies generalisations of the basic operators of calculus, namely differentiation and integration, defined by allowing their order to roam outside of Z to more general domains such as R and C [1,2,3]. Such generalisations are not only a pure mathematical curiosity, as differential equations involving the generalised operators have been applied in many scientific fields [4,5] ranging from viscoelasticity [6] to epidemiology [7] to economics [8,9] to electrical circuits [10].
The available methods for solving fractional partial differential equations are almost as diverse as their range of applications. To mention but a few examples, there are analytical approaches using transform methods [11,12], distribution theory [13,14,15], or weak solutions [16,17], as well as various numerical methods [18,19,20].
The classical fractional derivatives and integrals, named after Riemann and Liouville, are no longer the only definitions in the literature. As well as other classical definitions such as those of Riesz and Caputo, so many new definitions have been proposed in the 21st century [21,22] that it has become necessary to categorise them into general classes of operators for mathematical study [23]. The reasons for so many new operators being introduced are twofold, much like the reasons for fractional calculus itself being introduced: on the one hand, from the theoretical viewpoint, mathematical generalisation is its own motivation; on the other hand, from the applications viewpoint, different operators with different structures are needed in order to successfully model all the diverse behaviours and processes that exist in the real world.
One example of a general class of operators is given by fractional calculus with analytic kernel functions [24], defined in 2019, which can be related to the classical Riemann–Liouville operators via infinite series. Another example is given by fractional calculus with respect to functions, first theorised by Osler in 1970 [25], and later studied in more detail in the textbooks [3] (§18.2) and [26] (§2.5). The operators in this class can be related to the classical Riemann–Liouville operators via an algebraic conjugation relation, which we describe in more detail below. The two general classes just described can also be merged into a larger class, fractional calculus with analytic kernels with respect to functions [27], which is large enough to contain as special cases all of the operators in both of the above classes, as well as others such as the Hadamard-type fractional calculus which falls into neither of the above classes.
Another general class of operators, studied by Agrawal [28,29] in 2012, is known as weighted (or scaled) fractional calculus with respect to functions, and its applications in variational calculus [29] and probabilistic modelling [30] have been thoroughly explored. This can also be seen as a merger of two classes, namely fractional calculus with respect to functions (as mentioned above) and weighted fractional calculus (to be defined below). A few recent papers have also studied mathematical properties of these weighted operators [31,32,33] and of the associated differential equations [34,35,36].
Noticing this new pattern emerging in the literature, the authors were inspired to conduct a detailed study of the so-called weighted fractional calculus, in the spirit of Agrawal [28,29] and Kolokoltsov [30] who appreciated its operational properties, in order to understand its mathematical structure and provide results to make future studies much easier for anyone using these operators. In particular, by observing that a conjugation formula (similar to the one for fractional calculus with respect to functions) relates the operators of weighted fractional calculus with those of Riemann–Liouville fractional calculus, we can immediately deduce many results in the weighted theory from classical results on Riemann–Liouville. Our work, therefore, represents an undervalued perspective on weighted fractional calculus, which has not been noted in much of the recent literature on this topic.
The current work is structured as follows. Section 2 consists of preliminary definitions and results concerning the classical fractional calculus and fractional calculus with respect to functions. Section 3 consists of a detailed study of weighted fractional calculus, divided into subsections devoted to conjugation relations and their consequences, examples including the tempered and Kober–Erdélyi fractional calculi, and weighted versions of Laplace transforms and convolutions. Section 4 extends the results of the previous section to the larger class of operators given by weighted fractional calculus with respect to functions; this larger class can again be characterised by conjugation relations, and it includes examples such as the Hadamard-type and Erdélyi–Kober fractional calculi. In Section 5, we conclude with a summary and some pointers for future work on weighted fractional calculus.

2. Preliminaries

We firstly provide the definitions of the classical fractional integrals and derivatives of Riemann–Liouville and Caputo, then state some basic composition results and examples.
Definition 1
([1,2,3,37]). The Riemann–Liouville (RL) fractional integral of a function f L 1 ( a , b ) , to order α in R or C , is defined as
I a R L x α f ( x ) = 1 Γ ( α ) a x ( x t ) α 1 f ( t ) d t , x ( a , b ) ,
where we require Re ( α ) > 0 , or simply α > 0 if we assume real order.
Associated with this integral operator are two possible fractional derivative operators, named, respectively, after Riemann–Liouville (RL) and Caputo:
D a R L x α f ( x ) = d n d x n I a R L x n α f ( x ) , x ( a , b ) ,
D a C x α f ( x ) = I a R L x n α d n d x n f ( x ) , x ( a , b ) ,
where this time Re ( α ) 0 , or simply α > 0 if we assume real order, and f A C n [ a , b ] where n : = Re ( α ) + 1 so that n 1 Re ( α ) < n .
It is important to be aware that the Riemann–Liouville derivative (2) is the analytic continuation of the Riemann–Liouville integral (1) in the complex variable α, under the convention that integrals of negative order are derivatives of positive order:
D a R L x α f ( x ) = I a R L x α f ( x ) , Re ( α ) 0 .
This fact allows us to define both I a R L x α f ( x ) and D a R L x α f ( x ) for all values of α C , and apply principles of analytic continuation to extend many results from fractional integrals to fractional derivatives. In general, we may use the term “differintegral” to cover fractional integral and derivative operators together.
Lemma 1
([2,3]). Riemann–Liouville differintegrals have semigroup properties as follows:
I a R L x α I a R L x β f ( x ) = I a R L x α + β f ( x ) , α C , Re ( β ) > 0 ; d n d x n D a R L x α f ( x ) = D a R L x α + n f ( x ) , α C , n N ,
where in both cases f is any function such that the relevant expressions are well-defined.
Note that the operators labelled by α in both of these relations may be either fractional integrals or fractional derivatives, while the one labelled by β must be a fractional integral.
Lemma 2
([3,37]). The following composition properties are valid for Riemann–Liouville and Caputo differintegrals in cases where semigroup properties are not:
I a R L x α D a R L x α f ( x ) = f ( x ) k = 1 n ( x a ) α k Γ ( α k + 1 ) · lim x a + D a R L x α k f ( x ) ; I a R L x α D a C x α f ( x ) = f ( x ) k = 0 n 1 ( x a ) k k ! · lim x a + d k d x k f ( x ) ,
where in both cases α C with Re ( α ) > 0 and n = Re ( α ) + 1 while f is any function such that the relevant expressions are well-defined. We also have the following relationship between the Riemann–Liouville and Caputo derivatives:
D a C x α f ( x ) = D a R L x α f ( x ) k = 0 n 1 ( x a ) k α Γ ( k α + 1 ) · lim x a + d k d x k f ( x ) = D a R L x α f ( x ) k = 0 n 1 ( x a ) k k ! · lim x a + d k d x k f ( x ) ,
where f A C n [ a , b ] and α , n are as before.
Lemma 3
([1,2,37]). The Riemann–Liouville and Caputo differintegrals of certain functions are given as follows:
D a R L x α ( x a ) β = Γ ( β + 1 ) Γ ( β α + 1 ) ( x a ) β α , α C , Re ( β ) > 1 ; D a C x α E α ω ( x a ) α = ω E α ω ( x a ) α , ω C , Re ( α ) > 0 ,
where E α is the Mittag-Leffler function. Note that the operator in the first identity can be either a fractional integral or a fractional derivative, according to the sign of Re ( α ) .
Following the above basic introduction to RL and Caputo fractional calculus, we continue with a similar introduction to fractional calculus with respect to a function. Let ϕ : [ a , b ] R be a strictly increasing C 1 function, so that ϕ > 0 everywhere.
Definition 2
([3,25,26]). The Riemann–Liouville fractional integral of a function f L ϕ 1 ( a , b ) with respect to ϕ, to order α in R or C , is defined as
I a R L ϕ ( x ) α f ( x ) = 1 Γ ( α ) a x ϕ ( x ) ϕ ( t ) α 1 f ( t ) ϕ ( t ) d t , x ( a , b ) ,
where Re ( α ) > 0 , or simply α > 0 if we assume real order.
As before, associated with this integral operator are two possible fractional derivative operators, respectively, the Riemann–Liouville derivative with respect to ϕ and the Caputo derivative with respect to ϕ, and this time assuming that ϕ is a C function:
D a R L ϕ ( x ) α f ( x ) = 1 ϕ ( x ) · d d x n I a R L ϕ ( x ) n α f ( x ) , x ( a , b ) ,
D a C ϕ ( x ) α f ( x ) = I a R L ϕ ( x ) n α 1 ϕ ( x ) · d d x n f ( x ) , x ( a , b ) ,
where Re ( α ) 0 , or simply α > 0 if we assume real order, and f A C ϕ n [ a , b ] where n : = Re ( α ) + 1 so that n 1 Re ( α ) < n .
Here, we use the weighted L 1 space L ϕ 1 ( a , b ) = L 1 ( ( a , b ) , d ϕ ) and the spaces C ϕ n [ a , b ] and A C ϕ n [ a , b ] , which are, respectively, n times continuously differentiable or absolutely continuous with respect to the function ϕ.
Again, the Riemann–Liouville derivative (6) is the analytic continuation of the Riemann–Liouville integral (5) in the complex variable α, under the convention that integrals of negative order are derivatives of positive order. So, we can define both I a R L ϕ ( x ) α f ( x ) and D a R L ϕ ( x ) α f ( x ) for all values of α C , as before.
Lemma 4
([3,26]). Let ϕ C [ a , b ] be a function as above, and let Q ϕ be the functional operator of right-composition with ϕ: Q ϕ g = g ϕ for any function g defined on the interval [ ϕ ( a ) , ϕ ( b ) ] . The operators of fractional calculus with respect to functions can be expressed as algebraic conjugations, via the operator Q ϕ , of the classical fractional calculus operators:
I a R L ϕ ( x ) α = Q ϕ I ϕ ( a ) R L x α Q ϕ 1 , D a R L ϕ ( x ) α = Q ϕ D ϕ ( a ) R L x α Q ϕ 1 , D a C ϕ ( x ) α = Q ϕ D ϕ ( a ) C x α Q ϕ 1 .
Many results about fractional calculus with respect to functions can be obtained as immediate consequences of Lemma 4 together with corresponding results in classical fractional calculus. For example, the semigroup properties of Lemma 1 immediately give rise to corresponding semigroup properties for Riemann–Liouville differintegrals with respect to functions, simply by conjugation of the algebraic composition relations. Fractional differential equations with respect to functions can be transformed into classical fractional differential equations by simple substitutions equivalent to the conjugation relation [38]. Additionally, the examples of Lemma 3 immediately give rise to the following examples of Riemann–Liouville and Caputo differintegrals with respect to functions:
D a R L ϕ ( x ) α ϕ ( x ) ϕ ( a ) β = Γ ( β + 1 ) Γ ( β α + 1 ) ϕ ( x ) ϕ ( a ) β α , α C , Re ( β ) > 1 ;
D a C ϕ ( x ) α E α ω ϕ ( x ) ϕ ( a ) α = ω E α ω ϕ ( x ) ϕ ( a ) α , ω C , Re ( α ) > 0 .
We conclude this preliminaries section by showing the definition of fractional integrals with general analytic kernels.
Definition 3
([24]). For α , β C with positive real parts, let A be a kernel function given by a power series about 0 with radius of convergence at least ( b a ) Re ( β ) . The fractional integral with kernel A and parameters α , β of a function f L 1 ( a , b ) is defined as
I a A x α , β f ( x ) = a x ( x t ) α 1 A ( x t ) β f ( t ) d t , x ( a , b ) .
Fractional derivatives can also be defined in this general class, both of Riemann–Liouville type and of Caputo type, but we refer to [24] for the details of these.

3. Weighted Fractional Calculus

Definition 4.
The weighted Riemann–Liouville fractional integral of a given function f L 1 ( a , b ) , with a weight function w L ( a , b ) and order α in R or C , is defined by
I a R L x ; w ( x ) α f ( x ) = 1 Γ ( α ) w ( x ) a x ( x t ) α 1 w ( t ) f ( t ) d t , x ( a , b ) ,
where we require Re ( α ) > 0 , or simply α > 0 if we assume real order.
The weighted Riemann–Liouville fractional derivative of a given function f A C n [ a , b ] , with a weight function w A C n [ a , b ] and order α in R or C , is defined by
D a R L x ; w ( x ) α f ( x ) = d d x + w ( x ) w ( x ) n I a R L x ; w ( x ) n α f ( x ) , x ( a , b ) ,
where Re ( α ) 0 , or simply α > 0 if we assume real order, and n : = Re ( α ) + 1 so that n 1 Re ( α ) < n .
The weighted Caputo fractional derivative of a given function f A C n [ a , b ] , with a weight function w A C n [ a , b ] and order α in R or C , is defined by
D a C x ; w ( x ) α f ( x ) = I a R L x ; w ( x ) n α d d x + w ( x ) w ( x ) n f ( x ) ,
where Re ( α ) 0 , or simply α > 0 if we assume real order, and n : = Re ( α ) + 1 so that n 1 Re ( α ) < n .
Remark 1.
There are many possible function spaces that could be used as a setting or domain for the operators defined above. By Proposition 1 below, any suitable function space for the classical (Riemann–Liouville or Caputo) operators gives rise to a suitable function space for the corresponding weighted operators, simply by multiplying by the function w.
This explains the choice of function spaces in the above definition: if f L 1 ( a , b ) and w L ( a , b ) , then w · f L 1 ( a , b ) has a well-defined Riemann–Liouville integral; if f A C n [ a , b ] and w A C n [ a , b ] , then w · f A C n [ a , b ] has well-defined Riemann–Liouville and Caputo derivatives to order α.
Other function spaces can be used, but the condition for the w-weighted fractional operators to exist when applied to f is always precisely the condition for the corresponding classical operators to exist when applied to w · f .

3.1. Conjugation Relations

It is clear that the weighted Riemann–Liouville fractional integral with weight function w is given by multiplying by w, applying the original Riemann–Liouville fractional integral to the same order, and then dividing by w again. This gives rise to a conjugation relation for the weighted fractional integral, which we can also extend to weighted fractional derivatives in order to achieve the following result.
Proposition 1.
The weighted fractional differintegrals are conjugations of the original fractional differintegrals, as follows:
I a R L x ; w ( x ) α = M w ( x ) 1 I a R L x α M w ( x ) , D a R L x ; w ( x ) α = M w ( x ) 1 D a R L x α M w ( x ) , D a C x ; w ( x ) α = M w ( x ) 1 D a C x α M w ( x ) ,
where the operator M w ( x ) acting on functions is defined by multiplication:
M w ( x ) f ( x ) = w ( x ) f ( x ) .
Proof. 
The first result I a R L x ; w ( x ) α = M w ( x ) 1 I a R L x α M w ( x ) is clear. Since both types of weighted fractional derivatives are compositions of the weighted fractional integral with the operator d d x + w ( x ) w ( x ) repeated n times, it will suffice to show that this non-fractional operator also satisfies a conjugation relation, d d x + w ( x ) w ( x ) = M w ( x ) 1 d d x M w ( x ) . This is easily proved using the product rule:
M w ( x ) 1 d d x M w ( x ) f ( x ) = M w ( x ) 1 d d x w ( x ) f ( x ) = M w ( x ) 1 w ( x ) f ( x ) + w ( x ) f ( x ) = f ( x ) + w ( x ) w ( x ) f ( x ) = d d x + w ( x ) w ( x ) f ( x ) .
The results follow from composition of conjugation relations. □
The results of Proposition 1 are extremely useful in the study of weighted fractional calculus, because now many fundamental results from the original fractional calculus can be extended immediately, simply by composition, to the corresponding results on the weighted fractional operators. It is often true that the best way to prove facts about a new mathematical object is to express it in terms of an older one: this is the same principle that gave rise to the series formula for fractional calculus with general analytic kernels and to the conjugation relation for fractional calculus with respect to functions.
Proposition 2.
The weighted Riemann–Liouville derivative is the analytic continuation of the weighted Riemann–Liouville integral in the complex variable α, under the convention that integrals of negative order are derivatives of positive order:
D a R L x ; w ( x ) α = I a R L x ; w ( x ) α , Re ( α ) 0 .
This fact allows both I a R L x ; w ( x ) α and D a R L x ; w ( x ) α to be defined for all values of α C , in the same way as for the original Riemann–Liouville differintegrals.
Proof. 
This follows immediately from the conjugation relations of Proposition 1 together with the corresponding analytic continuation result for Riemann–Liouville differintegrals given at Equation (4). □
Proposition 3.
The weighted fractional differintegrals have semigroup properties as follows:
I a R L x ; w ( x ) α I a R L x ; w ( x ) β f ( x ) = I a R L x ; w ( x ) α + β f ( x ) , α C , Re ( β ) > 0 ; d d x + w ( x ) w ( x ) n D a R L x ; w ( x ) α f ( x ) = D a R L x ; w ( x ) n + α f ( x ) , α C , n N ,
where in both cases f is any function such that the relevant expressions are well-defined: for example, f L 1 ( a , b ) for the first identity, f A C n + Re α + 1 [ a , b ] for the second identity.
Note that the operators labelled by α in both of these relations may be either fractional integrals or fractional derivatives, while the one labelled by β must be a fractional integral.
Proof. 
This is an immediate consequence of Proposition 1 (conjugation relations) with Lemma 1 (Riemann–Liouville semigroup properties). □
Proposition 4.
The following composition properties are valid for weighted Riemann–Liouville and Caputo differintegrals in cases where semigroup properties are not:
I a R L x ; w ( x ) α D a R L x ; w ( x ) α f ( x ) = f ( x ) k = 1 n ( x a ) α k Γ ( α k + 1 ) · w ( a + ) w ( x ) · lim x a + D a R L x ; w ( x ) α k f ( x ) ; I a R L x ; w ( x ) α D a C x ; w ( x ) α f ( x ) = f ( x ) k = 0 n 1 ( x a ) k k ! · w ( a + ) w ( x ) · lim x a + d d x + w ( x ) w ( x ) k f ( x ) ,
where in both cases α C with Re ( α ) > 0 and n = Re ( α ) + 1 , while f is any function such that the relevant expressions are well-defined: specifically, f C n [ a , b ] is sufficient. We also have the following relationship between the weighted Riemann–Liouville and Caputo derivatives:
D a C x ; w ( x ) α f ( x ) = D a R L x ; w ( x ) α f ( x ) k = 0 n 1 ( x a ) k α Γ ( k α + 1 ) · w ( a + ) w ( x ) · lim x a + d d x + w ( x ) w ( x ) k f ( x ) = D a R L x ; w ( x ) α f ( x ) k = 0 n 1 ( x a ) k k ! · w ( a + ) w ( x ) · lim x a + d d x + w ( x ) w ( x ) k f ( x ) ,
where f A C n [ a , b ] and α , n are as before.
Proof. 
All of these relations follow directly from combining the conjugation relations of Proposition 1 with the previous relations of Lemma 2. However, the conjugation process is not quite so easy this time, as the final term (finite sum) is more tricky to deal with than single weighted fractional operators or compositions thereof. We verify the first of the stated relations carefully as follows:
I a R L x ; w ( x ) α D a R L x ; w ( x ) α f ( x ) = 1 w ( x ) I a R L x α D a R L x α w ( x ) f ( x ) = 1 w ( x ) w ( x ) f ( x ) k = 1 n ( x a ) α k Γ ( α k + 1 ) · lim x a + D a R L x α k w ( x ) f ( x ) = f ( x ) k = 1 n ( x a ) α k Γ ( α k + 1 ) w ( x ) · lim x a + w ( x ) · D a R L x ; w ( x ) α k f ( x ) = f ( x ) k = 1 n ( x a ) α k w ( a + ) Γ ( α k + 1 ) w ( x ) · lim x a + D a R L x ; w ( x ) α k f ( x ) ,
which is the required result. Similarly for the other relations stated in the proposition. □
Proposition 5.
The weighted Riemann–Liouville and Caputo differintegrals of certain functions are given as follows:
D a R L x ; w ( x ) α ( x a ) β w ( x ) = Γ ( β + 1 ) Γ ( β α + 1 ) ( x a ) β α w ( x ) , α C , Re ( β ) > 1 ;
D a C x ; w ( x ) α E α ω ( x a ) α w ( x ) = ω · E α ω ( x a ) α w ( x ) , ω C , Re ( α ) > 0 ,
where E α is the Mittag-Leffler function. Note that the operator in the first identity can be either a fractional integral or a fractional derivative, according to the sign of Re ( α ) .
Proof. 
The result follows directly from the conjugation relations given by Proposition 1, combined with the results given in Lemma 3. □
Remark 2.
Note that the functions used in Proposition 5 are just two possible examples that could have been chosen. Any known result for Riemann–Liouville or Caputo differintegrals of any particular functions can now easily be extended to an analogous result on weighted differintegrals, with the functions divided by w ( x ) on left and right sides of the identity.
Note also that, for any given smooth non-zero function h ( x ) on ( a , b ) , we can now find a weighted Caputo-type fractional derivative, which has this function as an eigenfunction, by choosing the weight function to be w ( x ) = E α ( ω ( x a ) α ) h ( x ) and using the second result of Proposition 5.

3.2. Examples

In this subsection, we discuss some particular choices of weight function w ( x ) , which lead to interesting special cases in the theory of weighted fractional calculus. We start with a trivial example and then proceed to two more examples, which give rise to models of fractional calculus that are already well-known and studied in the literature.
Example 1.
If w ( x ) = k is a constant, then the operators of weighted fractional calculus, as given in Definition 4, are exactly the same as the original Riemann–Liouville and Caputo operators given in Definition 1. This is because the operator M w ( x ) in this case commutes with the operators of fractional differentiation and integration, and indeed with any linear functional operator.
Note that this is the only possible type of function w ( x ) which returns the original Riemann–Liouville and Caputo operators, because it is the only case when the weighted 1st-order derivative d d x + w ( x ) w ( x ) is the same as the original 1st-order derivative d d x .
Remark 3.
More generally, we can conclude that, for a given weight function w ( x ) = w 0 ( x ) , any constant multiple w ( x ) = k w 0 ( x ) will give rise to exactly the same operators of weighted fractional calculus as w 0 ( x ) . In fact, this implication goes both ways: two weight functions w 1 ( x ) and w 2 ( x ) give rise to the same weighted operators if and only if they are constant multiples of each other. Therefore, we can say that the class of weighted fractional calculi bijects with the space of possible weight functions quotiented by the equivalence relation of constant multiplication.
Example 2.
If w ( x ) = e β x is an exponential function, then the operators of weighted fractional calculus, as given in Definition 4, are precisely those of tempered fractional calculus, defined [39,40] as follows:
I a T x α , β f ( x ) = 1 Γ ( α ) a x ( x t ) α 1 e β ( x t ) f ( t ) d t , β C , Re ( α ) > 0 ; D a T R x α , β f ( x ) = d d x + β n I a T x n α , β f ( x ) , β C , Re ( α ) 0 ; D a T C x α , β f ( x ) = I a T x n α , β d d x + β n f ( x ) , β C , Re ( α ) 0 .
Therefore, tempered fractional calculus is within the class of weighted fractional calculus. It is also known [41] that tempered fractional calculus is within the class of fractional calculus with analytic kernels. In fact, the intersection between these two general classes of fractional operators consists only of tempered fractional calculus, as we show in the following theorem.
Theorem 1.
The only operator which is a member of both the class of fractional integrals with general analytic kernels and the class of weighted fractional integrals is the tempered fractional integral. In other words, if
O = I a A x α , β = I a R L x ; w ( x ) α ,
with an analytic function A and a real continuous function w, then
O = I a T x α , κ
for some constant κ R .
Proof. 
We start with an operator O which is in both classes, and write the integrals explicitly as follows:
O f ( x ) = a x ( x t ) α 1 A ( x t ) β f ( t ) d t = 1 Γ ( α ) w ( x ) a x ( x t ) α 1 w ( t ) f ( t ) d t ,
for all f L 1 ( a , b ) and x ( a , b ) . Therefore, we must have
Γ ( α ) A ( x t ) β = w ( t ) w ( x ) , a < t < x < b .
The left-hand side is a function of x t ; let us write q ( z ) = Γ ( α ) A ( z β ) , noting that this is a smooth function on the real interval z [ 0 , b a ] . Without loss of generality, let us assume a = 0 as the lower end of the interval. Now we have q ( x t ) = w ( t ) w ( x ) for 0 < t < x < b . Letting t 0 , we find that q ( x ) = w ( 0 ) w ( x ) , so the equation becomes
w ( 0 ) w ( x t ) = w ( t ) w ( x ) , 0 < t < x < b .
Writing W ( z ) = w ( z ) w ( 0 ) (a constant multiple of the original function w), we obtain Cauchy’s functional equation:
W ( x t ) W ( t ) = W ( x ) , 0 < t < x < b .
Since we assumed w real and continuous, this means it must be a constant times an exponential function, w ( x ) = w ( 0 ) e κ x for some constant κ R , which gives the result. □
Example 3.
If w ( x ) = x η is a power function, then the operators of weighted fractional calculus, as given in Definition 4, are almost exactly those of the so-called Kober–Erdélyi (or simply Kober) fractional calculus, defined in [3] (§18.1) and [26] (§2.6) as follows:
I a K x α ; η f ( x ) = x α η Γ ( α ) a x ( x t ) α 1 t η f ( t ) d t , Re ( α ) > 0 , Re ( η ) > 0 .
(Usually, the Kober–Erdélyi operator is written with lower bound a = 0 , but we see no reason here to restrict to this particular case, as it can easily be written with the general constant of integration a.) It is clear that we have
I a K x α ; η f ( x ) = x α · I a R L x ; x η α f ( x ) .
This fact was already noted in [26] (Equation (2.6.15)), using a different notation of M η operators to give a relationship between the Kober–Erdélyi integral and the Riemann–Liouville integral.
A closely related, but more general, operator is the so-called Erdélyi–Kober fractional integral, which is also discussed in [3] (§18.1) and [26] (§2.6). However, this involves a fractional power of ( x σ t σ ) , not only ( x t ) , so it is not easily related to a special case of weighted fractional calculus as given in Definition 4. In fact, the appearance of this operator has elements of both weighted fractional calculus and fractional calculus with respect to functions, as given in Definition 2, and later we shall see that it is a special case of the larger class of operators described in Section 4 as weighted fractional calculus with respect to functions.

3.3. Laplace Transform and Convolution

In this subsection, we discuss a weighted integral transform, a very simple variation on the Laplace transform, which can be used to solve linear fractional differential equations involving weighted Riemann–Liouville and Caputo fractional derivatives. Related to this transform, we shall introduce a weighted convolution operation, which has a conjugation relation with the classical Laplace-type convolution. Using the conjugation approach and well-known facts about the classical Laplace transform and convolution, we can easily prove the fundamental properties of the weighted Laplace transform and weighted convolution.
Definition 5.
Let f : [ 0 , ) C be a real-valued or complex-valued function. Then, the weighted Laplace transform of f with a weight function w is defined by
L w ( x ) f ( x ) = F ( s ) = 0 e s x w ( x ) f ( x ) d x ,
for any s C and any function f such that this is a convergent integral (specific appropriate conditions on s and f are discussed below).
Remark 4.
It is clear that the w-weighted Laplace transform of f is precisely the usual Laplace transform of the function w ( x ) f ( x ) : in operational calculus terms, this means that
L w ( x ) = L M w ( x ) ,
where M w ( x ) is the multiplication operator defined by (9) above. This relation means that the w-weighted Laplace transform is essentially a trivial variation of the classical Laplace transform; there is no separate theory of weighted Laplace transforms, and no need to spend time proving anew any results on weighted Laplace transforms rather than deducing them directly from the theory of classical Laplace transforms.
We note that the function F ( s ) exists for Re ( s ) > c if w and f are such that w ( x ) f ( x ) is of exponential order c: for example, if w L [ 0 , ) and f is of exponential order c, or if w is of exponential order c 1 and f is of exponential order c 2 with c 1 + c 2 = c . In such circumstances, we say that f is w-exponentially bounded, which is the same as saying that w · f is exponentially bounded. This terminology is a convenience that will be useful later when we combine these concepts with those of Laplace transforms with respect to functions [42,43] and the corresponding function spaces.
As corollaries of the relation (11), we mention the following results.
Corollary 1.
The inverse weighted Laplace transform exists for any function which has a classical inverse Laplace transform, and it may be written as follows:
L w ( x ) 1 = M w ( x ) 1 L 1 ,
or in other words
L w ( x ) 1 F ( s ) = 1 2 π i w ( x ) c i c + i e s x F ( s ) d s .
Corollary 2.
If f is a function which has a classical Laplace transform F ( s ) , then the weighted Laplace transform of M w ( x ) 1 f ( x ) = f ( x ) w ( x ) is also F ( s ) :
L { f ( x ) } = F ( s ) L w ( x ) f ( x ) w ( x ) = F ( s ) .
The reason for introducing a weighted Laplace transform, even when it is so similar to the classical Laplace transform that all results concerning it are trivially proved, is that it has a natural relationship with the operators of weighted fractional calculus. In the following theorem, we find the weighted Laplace transforms of the weighted Riemann–Liouville and Caputo fractional differintegrals.
Theorem 2.
Let α > 0 and let f be a continuous function on [ 0 , ) which is of w-weighted exponential order, where w is a continuous weight function. Then, we have the following results.
1. 
L w ( x ) I 0 R L x ; w ( x ) α f ( x ) = s α L w ( x ) f ( x ) .
2. 
Let n 1 Re ( α ) < n Z + , and assume that D 0 R L x ; w ( x ) α f is continuous on [ 0 , ) and of w-weighted exponential order. Then,
L w ( x ) D 0 R L x ; w ( x ) α f ( x ) = s α L w ( x ) { f ( x ) } w ( 0 + ) i = 0 n 1 s n i 1 I 0 R L x ; w ( x ) n i α f ( 0 + ) .
3. 
Let n 1 Re ( α ) < n Z + , and assume that D x ; w ( x ) n f is continuous on [ 0 , ) and of w-weighted exponential order. Then,
L w ( x ) D 0 C x ; w ( x ) α f ( x ) = s α L w ( x ) { f ( x ) } w ( 0 + ) i = 0 n 1 s α i 1 d d x + w ( x ) w ( x ) i f ( 0 + ) .
Proof. 
All of these results follow from combining the results of Proposition 1 and Equation (11) with the following well-known facts on classical Laplace transforms of fractional integrals and derivatives:
L I 0 R L x α f ( x ) = s α L f ( x ) , L D 0 R L x α f ( x ) = s α L { f ( x ) } i = 0 n 1 s n i 1 I 0 x n i α f ( 0 + ) , L D 0 C x α f ( x ) = s α L { f ( x ) } i = 0 n 1 s α i 1 f ( i ) ( 0 + ) .
Note that the given conditions on f imply that (in the second case, Riemann–Liouville derivatives) the function I 0 R L x ; w ( x ) n α f and all its weighted derivatives from 1st to nth order, namely D 0 x ; w ( x ) α n + 1 f , , D 0 R L x ; w ( x ) α 1 f , D 0 R L x ; w ( x ) α f , are continuous on [ 0 , ) and of w-weighted exponential order, while (in the third case, Caputo derivatives) the function f and all its weighted derivatives from 1st to nth order, namely D 0 R L x ; w ( x ) 1 f , D 0 R L x ; w ( x ) 2 f , , D 0 R L x ; w ( x ) n f , are continuous on [ 0 , ) and of w-weighted exponential order, and so is D 0 C x ; w ( x ) α f . This is because continuity and exponential boundedness are preserved by integration; we have fixed our assumptions to be only at the highest order of differentiation. □
Now, we consider a w-weighted convolution operation between two functions, which relates naturally to the weighted Laplace transform.
Definition 6.
The w-weighted convolution of two real-valued or complex-valued functions f , g : [ 0 , ) C is the function f w g defined by
f w g ( x ) = 1 w ( x ) 0 x w ( x t ) f x t w ( t ) g ( t ) d t .
It is clear that this weighted convolution is closely related to classical convolution via the following formula:
f w g ( x ) = 1 w ( x ) w f w g ( x ) ,
or in other words, using the multiplication operator defined in (9),
f w g = M w ( x ) 1 M w ( x ) f M w ( x ) g .
Since both classical and w-weighted convolutions are binary operations acting on pairs of functions, we can use the alternative notations ( f , g ) and w ( f , g ) , instead of f g and f w g , in order to rewrite (15) as a conjugation relation of operators:
w = M w ( x ) 1 M w ( x ) , M w ( x ) ,
i.e.,
w f , g = M w ( x ) 1 M w ( x ) f , M w ( x ) g .
These relations, combined with the operational relation (11) for the weighted Laplace transforms, give rise to the following important corollary, an analogue of the classical Laplace convolution theorem.
Corollary 3.
If f , g : [ 0 , X ] C are piecewise continuous and their products with w are of exponential order c > 0 , then
L w ( x ) f w g = L w ( x ) { f } L w ( x ) { g } .
In the following theorem, we use the weighted Laplace transform method to establish a regularity condition for the solutions to weighted fractional differential equations of the following type:
D 0 C x ; w ( x ) α y ( x ) = A y ( x ) + g ( x ) , x 0 ,
y ( 0 ) = η ,
where 0 < α < 1 is fixed and A = a i j is an n × n constant matrix and g is a continuous n-dimensional vector-valued function and η is a constant n-dimensional vector.
Theorem 3.
Assume that the system (16)–(17) has a unique continuous solution y . If g is continuous on [ 0 , ) and w-weighted exponentially bounded, then y and D 0 C x ; w ( x ) α y are both w-weighted exponentially bounded too.
Proof. 
Using the result proved in [44] (Theorem 3.1) together with the conjugation relations given by Proposition 1, we obtain the required result. □

4. Weighted Fractional Calculus with Respect to Functions

Definition 7
([28,33]). Let ϕ : [ a , b ] R be a strictly increasing C 1 function, so that ϕ > 0 everywhere, and let w L ( a , b ) be a weight function. The w-weighted Riemann–Liouville fractional integral with respect to ϕ of a given function f L ϕ 1 ( a , b ) , to order α in R or C , is defined by
I a R L ϕ ( x ) ; w ( x ) α f ( x ) = 1 Γ ( α ) w ( x ) a x ϕ ( x ) ϕ ( t ) α 1 w ( t ) f ( t ) ϕ ( t ) d t , x ( a , b ) ,
where we require Re ( α ) > 0 , or simply α > 0 if we assume real order.
Assuming sufficient differentiability conditions on ϕ and w, we can also define the w-weighted Riemann–Liouville fractional derivative with respect to ϕ of a given function f A C ϕ n [ a , b ] and the w-weighted Caputo fractional derivative with respect to ϕ of a given function f C ϕ n [ a , b ] , to order α in R or C , as follows:
D a R L ϕ ( x ) ; w ( x ) α f ( x ) = D ϕ ( x ) ; w ( x ) n I a R L ϕ ( x ) ; w ( x ) n α f ( x ) , D a C ϕ ( x ) ; w ( x ) α f ( x ) = I a R L ϕ ( x ) ; w ( x ) n α D ϕ ( x ) ; w ( x ) n f ( x ) ,
where Re ( α ) 0 , or simply α > 0 if we assume real order, and n : = Re ( α ) + 1 so that n 1 Re ( α ) < n , and where the first-order operator D ϕ ( x ) ; w ( x ) is defined by
D ϕ ( x ) ; w ( x ) f ( x ) = 1 w ( x ) ϕ ( x ) · d d x w ( x ) f ( x ) ,
or equivalently D ϕ ( x ) ; w ( x ) = 1 ϕ ( x ) · d d x + w ( x ) w ( x ) .
Sufficient differentiability conditions on ϕ and w could be, for example, that both ϕ and w are in C n [ a , b ] , but the important condition overall is that the functions D ϕ ( x ) ; w ( x ) k f ( x ) should exist for k = 1 , 2 , , n and that the k = n case should be, respectively, continuous or L 1 with respect to ϕ.
These operators were introduced by Agrawal [28,29] and also studied by Jarad et al. [33]. Here, we shall consider them from the operational viewpoint, using the notion of conjugation to connect them directly to the original Riemann–Liouville and Caputo operators and thence to derive many results in an easy way without the unnecessary calculations of some previous works [33].

4.1. Conjugation Relations

Theorem 4
([30] (§3.3)). The weighted fractional differintegrals with respect to functions may be written as conjugations of the original fractional differintegrals as follows:
I a R L ϕ ( x ) ; w ( x ) α = M w ( x ) 1 Q ϕ I ϕ ( a ) R L x α Q ϕ 1 M w ( x ) , D a R L ϕ ( x ) ; w ( x ) α = M w ( x ) 1 Q ϕ D ϕ ( a ) R L x α Q ϕ 1 M w ( x ) , D a C ϕ ( x ) ; w ( x ) α = M w ( x ) 1 Q ϕ D ϕ ( a ) C x α Q ϕ 1 M w ( x ) ,
where the multiplication operator M w ( x ) is defined by (9) and the composition operator Q ϕ is defined in Lemma 4.
Proof. 
It is clear that the w-weighted Riemann–Liouville fractional integral with respect to ϕ is given by multiplying by w, applying the (unweighted) Riemann–Liouville fractional integral with respect to ϕ to the same order, and then dividing by w again. This gives a conjugation relation between the weighted fractional integral with respect to ϕ and the original fractional integral with respect to ϕ , which can be combined with Lemma 4 as follows:
I a R L ϕ ( x ) ; w ( x ) α = M w ( x ) 1 I a R L ϕ ( x ) α M w ( x ) = M w ( x ) 1 Q ϕ I ϕ ( a ) R L x α Q ϕ 1 M w ( x ) .
Both types of fractional derivatives are compositions of this fractional integral operator with the first-order operator D ϕ ( x ) ; w ( x ) = 1 ϕ ( x ) · d d x + w ( x ) w ( x ) repeated n times, so it will suffice to show that D ϕ ( x ) ; w ( x ) also satisfies a conjugation relation, which is easily proved using the product rule and chain rule:
M w ( x ) 1 Q ϕ d d x Q ϕ 1 M w ( x ) f ( x ) = M w ( x ) 1 1 ϕ ( x ) · d d x M w ( x ) f ( x ) = M w ( x ) 1 1 ϕ ( x ) · d d x w ( x ) f ( x ) = 1 w ( x ) ϕ ( x ) w ( x ) f ( x ) + w ( x ) f ( x ) = 1 ϕ ( x ) d d x + w ( x ) w ( x ) f ( x ) .
The results follow from composition of conjugation relations. □
Remark 5.
The above result shows that the operators of Definition 7 are obtained by starting with the original operators of fractional calculus, applying Q ϕ -conjugation to make the operators with respect to functions, and then applying M w -conjugation to make them weighted. There is a definite order to the modifications: we are using the weighted versions of operators with respect to functions. What about the with-respect-to-functions versions of weighted operators, as mentioned briefly in [30] (Equations (3.36)–(3.38))?
In fact, these two potential general classes of operators come to the same thing, because of the relation M w ( x ) Q ϕ 1 = Q ϕ 1 M w ( ϕ ( x ) ) . It can be checked by some elementary calculations that
Q ϕ M w ( x ) 1 d d x M w ( x ) Q ϕ 1 f ( x ) = 1 ϕ ( x ) d d x + ( w ϕ ) ( x ) w ϕ ( x ) f ( x ) ,
confirming that the class of weighted fractional operators with respect to functions is the same class regardless of which order we apply the weighted and with-respect-to-functions aspects of the operators. To write the alternative conjugation relations explicitly, we have
I a R L ϕ ( x ) ; w ( x ) α = Q ϕ M w ( ϕ 1 ( x ) ) 1 I ϕ ( a ) R L x α M w ( ϕ 1 ( x ) ) Q ϕ 1 , D a R L ϕ ( x ) ; w ( x ) α = Q ϕ M w ( ϕ 1 ( x ) ) 1 D ϕ ( a ) R L x α M w ( ϕ 1 ( x ) ) Q ϕ 1 , D a C ϕ ( x ) ; w ( x ) α = Q ϕ M w ( ϕ 1 ( x ) ) 1 D ϕ ( a ) C x α M w ( ϕ 1 ( x ) ) Q ϕ 1 ,
derived from the results of Theorem 4.
The above conjugation results are very useful in understanding weighted fractional calculus with respect to functions. They were mentioned in [30] (Equations (3.31)–(3.35)), but they were not used in [33], many of whose results can be proved much more quickly by directly using the corresponding results from classical fractional calculus together with the conjugation relations.
Proposition 6.
The w-weighted RL derivative with respect to ϕ is the analytic continuation, in the complex variable α, of the w-weighted RL integral with respect to ϕ, under the convention that integrals of negative order are derivatives of positive order:
D a R L ϕ ( x ) ; w ( x ) α = I a R L ϕ ( x ) ; w ( x ) α , Re ( α ) 0 .
This fact allows both I a R L ϕ ( x ) ; w ( x ) α and D a R L ϕ ( x ) ; w ( x ) α to be defined for all values of α C , in the same way as for the original Riemann–Liouville differintegrals.
Proof. 
This follows directly from the conjugation relations of Theorem 4 together with the corresponding analytic continuation result for Riemann–Liouville differintegrals given at Equation (4). □
Proposition 7.
The weighted fractional differintegrals with respect to functions have semigroup properties as follows:
I a R L ϕ ( x ) ; w ( x ) α I a R L ϕ ( x ) ; w ( x ) β f ( x ) = I a R L ϕ ( x ) ; w ( x ) α + β f ( x ) , α C , Re ( β ) > 0 ; D ϕ ( x ) ; w ( x ) n D a R L ϕ ( x ) ; w ( x ) α f ( x ) = D a R L ϕ ( x ) ; w ( x ) n + α f ( x ) , α C , n N ,
where in both cases f is any function such that the relevant expressions are well-defined: for example, f L ϕ 1 ( a , b ) for the first identity, f A C ϕ n + Re α + 1 [ a , b ] for the second identity.
Note that the operators labelled by α in both of these relations may be either fractional integrals or fractional derivatives, while the one labelled by β must be a fractional integral.
Proof. 
This is an immediate consequence of Theorem 4 (conjugation relations) with Lemma 1 (Riemann–Liouville semigroup properties). □
Proposition 8.
The following composition properties are valid for weighted Riemann–Liouville and Caputo differintegrals with respect to functions in cases where semigroup properties are not:
I a R L ϕ ( x ) ; w ( x ) α D a R L ϕ ( x ) ; w ( x ) α f ( x ) = f ( x ) k = 1 n ϕ ( x ) ϕ ( a ) α k Γ ( α k + 1 ) · w ( a + ) w ( x ) · lim x a + D a R L ϕ ( x ) ; w ( x ) α k f ( x ) ;
I a R L ϕ ( x ) ; w ( x ) α D a C ϕ ( x ) ; w ( x ) α f ( x ) = f ( x ) k = 0 n 1 ϕ ( x ) ϕ ( a ) k k ! · w ( a + ) w ( x ) · lim x a + D ϕ ( x ) ; w ( x ) k f ( x ) ,
where in both cases α C with Re ( α ) > 0 and n = Re ( α ) + 1 while f is any function such that the relevant expressions are well-defined: specifically, f C ϕ n [ a , b ] is sufficient. We also have the following relationship between the weighted Riemann–Liouville and Caputo derivatives with respect to functions:
D a C ϕ ( x ) ; w ( x ) α f ( x ) = D a R L ϕ ( x ) ; w ( x ) α f ( x ) k = 0 n 1 ϕ ( x ) ϕ ( a ) k α Γ ( k α + 1 ) · w ( a + ) w ( x ) · lim x a + D ϕ ( x ) ; w ( x ) k f ( x ) = D a R L ϕ ( x ) ; w ( x ) α f ( x ) k = 0 n 1 ϕ ( x ) ϕ ( a ) k k ! · w ( a + ) w ( x ) · lim x a + D ϕ ( x ) ; w ( x ) k f ( x ) ,
where f A C ϕ n [ a , b ] and α , n are as before.
Proof. 
We verify the first of the stated relations as follows:
I a R L ϕ ( x ) ; w ( x ) α D a R L ϕ ( x ) ; w ( x ) α f ( x ) = M w ( x ) 1 Q ϕ I ϕ ( a ) R L x α D ϕ ( a ) R L x α Q ϕ 1 M w ( x ) f ( x ) = M w ( x ) 1 Q ϕ I ϕ ( a ) R L x α D ϕ ( a ) R L x α w ϕ 1 ( x ) f ϕ 1 ( x ) = M w ( x ) 1 Q ϕ [ w ϕ 1 ( x ) f ϕ 1 ( x ) k = 1 n x ϕ ( a ) α k Γ ( α k + 1 ) · lim x ϕ ( a ) + D ϕ ( a ) R L x α k w ϕ 1 ( x ) f ϕ 1 ( x ) ] = 1 w ( x ) [ w ( x ) f ( x ) k = 1 n ϕ ( x ) ϕ ( a ) α k Γ ( α k + 1 ) × lim x ϕ ( a ) + D ϕ ( a ) R L x α k w ϕ 1 ( x ) f ϕ 1 ( x ) ] = f ( x ) k = 1 n ϕ ( x ) ϕ ( a ) α k Γ ( α k + 1 ) w ( x ) · lim x ϕ ( a ) + D ϕ ( a ) R L x α k Q ϕ 1 M w ( x ) f ( x ) = f ( x ) k = 1 n ϕ ( x ) ϕ ( a ) α k Γ ( α k + 1 ) w ( x ) · lim x a + Q ϕ D ϕ ( a ) R L x α k Q ϕ 1 M w ( x ) f ( x ) = f ( x ) k = 1 n ϕ ( x ) ϕ ( a ) α k Γ ( α k + 1 ) · w ( a + ) w ( x ) · lim x a + D a R L ϕ ( x ) ; w ( x ) α k f ( x ) ,
which is the required result. Similar manipulations can be used to prove all of the other stated relations, starting from the corresponding relations for Riemann–Liouville and Caputo differintegrals; we omit the straightforward details. □
Proposition 9.
The weighted Riemann–Liouville and Caputo differintegrals of certain functions with respect to another function are given as follows:
D a R L ϕ ( x ) ; w ( x ) α ϕ ( x ) ϕ ( a ) β w ( x ) = Γ ( β + 1 ) Γ ( β α + 1 ) ϕ ( x ) ϕ ( a ) β α w ( x ) , α C , Re ( β ) > 1 ;
D a C ϕ ( x ) ; w ( x ) α E α ω ϕ ( x ) ϕ ( a ) α w ( x ) = ω · E α ω ϕ ( x ) ϕ ( a ) α w ( x ) , ω C , Re ( α ) > 0 ,
where E α is the Mittag-Leffler function. Note that the operator in the first identity can be either a fractional integral or a fractional derivative, according to the sign of Re ( α ) .
Proof. 
The result follows directly from the conjugation relations of Theorem 4, combined with the results of Lemma 3. □
Remark 6.
As before, the functions used in Proposition 9 are just two possible examples that could have been chosen. Any known result for Riemann–Liouville or Caputo differintegrals of any particular functions can now easily be extended to an analogous result on weighted differintegrals with respect to functions, with the functions divided by w ( x ) on left and right sides of the identity.

4.2. Examples

In this subsection, we discuss some particular choices of functions for both w ( x ) and ϕ ( x ) which will lead to interesting special cases of weighted fractional calculus with respect to functions, some of which are already well known and studied in the literature.
Example 4.
If w ( x ) = k is a constant, then the operators of w-weighted fractional calculus with respect to ϕ, as given in Definition 7, are exactly the same as the operators of the original fractional calculus with respect to ϕ, as given in Definition 2. Indeed, this is the only case where w-weighted fractional calculus with respect to ϕ reduces to fractional calculus with respect to a function, because it is the only case when the first-order operator D ϕ ( x ) ; w ( x ) = 1 ϕ ( x ) · d d x + w ( x ) w ( x ) becomes simply a function times d d x .
If ϕ ( x ) = x , then the operators of w-weighted fractional calculus with respect to ϕ, as given in Definition 7, are exactly the same as the operators of w-weighted fractional calculus, as given in Definition 4. Indeed, this is the only case where w-weighted fractional calculus with respect to ϕ reduces to a case of weighted fractional calculus, because it is the only case when the first-order operator D ϕ ( x ) ; w ( x ) = 1 ϕ ( x ) · d d x + w ( x ) w ( x ) contains simply the operator d d x without any function multiplier.
Example 5.
If w ( x ) = e β ϕ ( x ) is an exponential function of ϕ, then the operators of w-weighted fractional calculus with respect to ϕ, as given in Definition 7, are precisely those of tempered fractional calculus with respect to a function, as defined in [45]. Therefore, tempered fractional calculus with respect to a function forms an overlap between the general class of weighted fractional operators with respect to functions [33] and the general class of fractional operators with analytic kernels with respect to functions [27].
Example 6.
If w ( x ) = x β is a power function and ϕ ( x ) = log ( x ) is the natural logarithm function, then the operators of w-weighted fractional calculus with respect to ϕ become precisely those of Hadamard-type fractional calculus, defined [46,47] as follows:
I a H x α , β f ( x ) = 1 Γ ( α ) a x t x β log x t α 1 f ( t ) t d t , β C , Re ( α ) > 0 ; D a H R x α , β f ( x ) = x · d d x + β n I a H x n α , β f ( x ) , β C , Re ( α ) 0 ; D a H C x α , β f ( x ) = I a H x n α , β x · d d x + β n f ( x ) , β C , Re ( α ) 0 .
This model of fractional calculus was also studied in [45], where it was determined to be a special case of tempered fractional calculus with respect to a function, namely the case with respect to the logarithm function.
Remark 7.
Note that choosing w ( x ) = e β ϕ ( x ) in Example 5, rather than w ( x ) = e β x as we did to obtain tempered fractional calculus as a special case of weighted fractional calculus in Example 2, is necessary because we are defining weighted fractional calculus with respect to functions by applying the M w conjugation after the Q ϕ conjugation. As discussed in Remark 5, this means that our operators of Definition 7 are
weighted ( fractional calculus with respect to functions )
and not
( weighted fractional calculus ) with respect to functions .
Thus, choosing w ( x ) = e β x in Definition 4 gives tempered fractional calculus, but choosing w ( x ) = e β x in Definition 7 does not give the with-respect-to-functions version of tempered fractional calculus: instead, it gives the tempered version of fractional calculus with respect to functions. On the other hand, if we swapped the order of the M w and Q ϕ operators in the definition, then Hadamard-type fractional calculus would be given by w ( x ) = e β x and ϕ ( x ) = log ( x ) , as it is tempered fractional calculus taken with respect to the natural logarithm function.
Explicitly, the operators of Hadamard-type fractional calculus are given by the following conjugation relations:
I a H x α , β f ( x ) = M x β 1 Q log I log ( a ) R L x α Q log 1 M x β = Q log M e β x 1 I log ( a ) R L x α M e β x Q log 1 ; D a H R x α , β f ( x ) = M x β 1 Q log D log ( a ) R L x α Q log 1 M x β = Q log M e β x 1 D log ( a ) R L x α M e β x Q log 1 ; D a H C x α , β f ( x ) = M x β 1 Q log D log ( a ) C x α Q log 1 M x β = Q log M e β x 1 D log ( a ) C x α M e β x Q log 1 ,
which was already known from [45] (Theorem 3.2).
Example 7.
If w ( x ) = x σ η and ϕ ( x ) = x σ are power functions, with Re ( η ) > 0 and σ > 0 , then the operators of w-weighted fractional calculus with respect to ϕ, as given in Definition 7, are almost exactly those of the so-called Erdélyi–Kober fractional calculus, which we define following [3] (§18.1) and [26] (§2.6) as follows. The fractional integral is
I a E x α ; σ , η f ( x ) = σ x σ ( α + η ) Γ ( α ) a x x σ t σ α 1 t σ η + σ 1 f ( t ) d t , Re ( α ) > 0 ,
while the fractional derivative (of Riemann–Liouville type) is
D a E R x α ; σ , η f ( x ) = x σ η 1 σ x σ 1 · d d x n x σ ( η + n ) I a E x n α ; σ , η + α f ( x ) , Re ( α ) 0 ,
and the fractional derivative of Caputo type, defined more recently in [48], is:
D a E C x α ; σ , η f ( x ) = x σ n I a E x n α ; σ , η + α x σ ( η + α ) 1 σ x σ 1 · d d x n x σ ( η + α ) f ( x ) , Re ( α ) 0 ,
where in both of the last two cases the natural number n is defined by n 1 < Re ( α ) n , or in other words n = Re ( α ) + 1 .
It is clear that the Erdélyi–Kober integral is related to the w-weighted fractional integral with respect to ϕ as follows:
I a E x α ; σ , η f ( x ) = x σ α · I a R L x σ ; x σ η α f ( x ) .
The original Erdélyi–Kober derivative can also be related to the w-weighted fractional derivative with respect to ϕ of Riemann–Liouville type, as follows:
D a E R x α ; σ , η f ( x ) = x σ α · D a R L x σ ; x σ ( η + α ) α f ( x ) .
And the Caputo-type Erdélyi–Kober derivative can be related similarly to the w-weighted fractional derivative with respect to ϕ of Caputo type, as follows:
D a E C x α ; σ , η f ( x ) = x σ α · D a C x σ ; x σ ( η + α ) α f ( x ) .
These relationships were already noted in [26] (Equation (2.6.9)), using a different notation of M η and N σ operators to give a relationship between the Erdélyi–Kober integral and the Riemann–Liouville integral, and in [48] (Equations (16)–(24)), using a direct formulation involving substitutions of power functions.
It is also interesting to note that, by our Proposition 6, the Erdélyi–Kober fractional derivative is the unique analytic continuation of the Erdélyi–Kober fractional integral I a E x α ; σ , η f ( x ) from the original domain Re ( α ) > 0 to the whole complex plane for α, under the convention that
I a E x α ; σ , η f ( x ) = D a E R x α ; σ , η + α f ( x ) , Re ( α ) 0 .
As usual with analytic continuation relations, this fact will be generally useful in proving results about Erdélyi–Kober derivatives when the corresponding results for Erdélyi–Kober integrals are already known, simply by extension using the fact of analytic continuation.

4.3. Laplace Transform and Convolution

We shall now study the w-weighted Laplace transform with respect to ϕ , an integral transform which is ideally suited for studying fractional differential equations which are both weighted and with respect to functions. Most of the results of this subsection were already seen in [33], but we shall now see how to prove them much more quickly and easily by using operational calculus.
Definition 8
([33]). Let f : [ a , ) C be a real-valued or complex-valued function, and let w and ϕ be functions as above. Then, the w-weighted Laplace transform of f with respect to ϕ is defined by
L ϕ ; w ( x ) f ( x ) = F ( s ) = a e s [ ϕ ( x ) ϕ ( a ) ] w ( x ) f ( x ) ϕ ( x ) d x ,
for any s C and any function f such that this is a convergent integral (specific appropriate conditions on s and f are discussed below).
Theorem 5.
The w-weighted Laplace transform with respect to ϕ can be written as a composition of the usual Laplace transform with multiplication and composition operators, as follows:
L ϕ ; w ( x ) = L Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) ,
where M and Q are the operators defined above in (9) and in Lemma 4, respectively, and where we assume the increasing function ϕ satisfies ϕ ( x ) as x .
Proof. 
We check the effect of applying the three operators from the right-hand side, one by one, on an appropriate function f:
M w ( x ) f ( x ) = w ( x ) f ( x ) ; Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) f ( x ) = w ϕ 1 x + ϕ ( a ) f ϕ 1 x + ϕ ( a ) ; L Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) f ( x ) = 0 e s x w ϕ 1 x + ϕ ( a ) f ϕ 1 x + ϕ ( a ) d x = a e s [ ϕ ( t ) ϕ ( a ) ] w ( t ) f ( t ) ϕ ( t ) d t ,
where in the last step we substituted x = ϕ ( t ) ϕ ( a ) and used the assumption on the infinite limiting behaviour of ϕ . □
Corollary 4.
If ϕ : [ 0 , ) [ 0 , ) is an increasing bijection, then the w-weighted Laplace transform with respect to ϕ can be written as a composition of the usual Laplace transform with multiplication and composition operators, as follows:
L ϕ ; w ( x ) = L Q ϕ 1 M w ( x ) ,
where M w ( x ) and Q ϕ are the operators defined above in (9) and in Lemma 4, respectively.
Proof. 
This is the case a = 0 , ϕ ( 0 ) = 0 of the preceding theorem. □
Corollary 5.
If a R is any real number and ϕ : [ a , ) [ ϕ ( a ) , ) is a bijection, then the Laplace transform with respect to ϕ can be written in terms of the usual Laplace transform as follows:
L ϕ = L Q ϕ ( x ) ϕ ( a ) 1 .
Proof. 
This is the case w ( x ) = 1 of the preceding theorem. Note that it also provides a generalisation of the results of [45], where it was assumed that ϕ ( 0 ) = 0 . The operational-calculus approach in [45] was therefore used only with fractional differintegrals having lower limit 0, but now those results can be straightforwardly extended to fractional differintegral operators with general lower limit a R . □
As corollaries of the relation (19), we obtain the following results.
Corollary 6.
If w and f are such that w ( x ) f ( x ) is of ϕ-exponential order c (as defined in [45]), then the w-weighted Laplace transform F ( s ) of f with respect to ϕ exists for Re ( s ) > c . In this case, we say that f is of w-weighted ϕ-exponential order c.
Corollary 7.
The inverse w-weighted Laplace transform with respect to ϕ exists for any function which has a classical inverse Laplace transform, and it may be written as follows:
L ϕ ; w ( x ) 1 = M w ( x ) 1 Q ϕ ( x ) ϕ ( a ) L 1 ,
or in other words
L ϕ ; w ( x ) 1 F ( s ) = 1 2 π i w ( x ) c i c + i e s [ ϕ ( x ) ϕ ( a ) ] F ( s ) d s .
Corollary 8.
If f is a function which has a classical Laplace transform F ( s ) , then the w-weighted Laplace transform with respect to ϕ of the function M w ( x ) 1 Q ϕ ( x ) ϕ ( a ) f ( x ) = f ϕ ( x ) ϕ ( a ) w ( x ) is also F ( s ) :
L { f ( x ) } = F ( s ) L ϕ ; w ( x ) f ϕ ( x ) ϕ ( a ) w ( x ) = F ( s ) .
The following theorem gives the natural relationship between the weighted Laplace transform with respect to a function and the operators of weighted fractional calculus with respect to a function.
Theorem 6.
Let α > 0 and let f be a continuous function on [ 0 , ) , which is of w-weighted ϕ-exponential order, where w is a continuous weight function and ϕ : [ a , ) [ ϕ ( a ) , ) is an increasing bijection. Then, we have the following results.
1. 
L ϕ ; w ( x ) I a R L ϕ ( x ) ; w ( x ) α f ( x ) = s α L ϕ ; w ( x ) f ( x ) .
2. 
Let n 1 Re ( α ) < n Z + , and assume that D a R L ϕ ( x ) ; w ( x ) α f is continuous on [ a , ) and of w-weighted ϕ-exponential order. Then,
L ϕ ; w ( x ) D a R L ϕ ( x ) ; w ( x ) α f ( x ) = s α L ϕ ; w ( x ) { f ( x ) } w ( a + ) i = 0 n 1 s n i 1 I a R L ϕ ( x ) ; w ( x ) n i α f ( a + ) .
3. 
Let n 1 Re ( α ) < n Z + , and assume that D ϕ ( x ) ; w ( x ) n f is continuous on [ a , ) and of w-weighted ϕ-exponential order. Then,
L ϕ ; w ( x ) D a C ϕ ( x ) ; w ( x ) α f ( x ) = s α L ϕ ; w ( x ) { f ( x ) } w ( a + ) i = 0 n 1 s α i 1 D ϕ ( x ) ; w ( x ) i f ( a + ) .
Proof. 
All of these results follow from combining the composition results of Theorem 4 and Theorem 5 with the classical facts on Laplace transforms of fractional integrals and derivatives which were quoted previously in the proof of Theorem 2. However, in this case, the Q operators are a little more tricky to deal with, since we have both Q ϕ and Q ϕ ( x ) ϕ ( a ) , closely related but not the same operator, involved in the same manipulation.
For the fractional integral, we have:
L ϕ ; w ( x ) I a R L ϕ ( x ) ; w ( x ) α = L Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) M w ( x ) 1 Q ϕ I ϕ ( a ) R L x α Q ϕ 1 M w ( x ) = L Q ϕ ( x ) ϕ ( a ) 1 Q ϕ I ϕ ( a ) R L x α Q ϕ 1 M w ( x ) ,
and the composition Q ϕ ( x ) ϕ ( a ) 1 Q ϕ is equivalent to a linear substitution:
Q ϕ ( x ) ϕ ( a ) 1 Q ϕ g ( x ) = g x + ϕ ( a ) ,
so this operator has the following effect on the fractional integral operator:
Q ϕ ( x ) ϕ ( a ) 1 Q ϕ I ϕ ( a ) R L x α g ( x ) = 1 Γ ( α ) ϕ ( a ) x + ϕ ( a ) x + ϕ ( a ) t α 1 g ( t ) d t = 1 Γ ( α ) 0 x ( x u ) α 1 g u + ϕ ( a ) d u = I 0 R L x α Q ϕ ( x ) ϕ ( a ) 1 Q ϕ g ( x ) ,
and therefore
L ϕ ; w ( x ) I a R L ϕ ( x ) ; w ( x ) α = L I 0 R L x α Q ϕ ( x ) ϕ ( a ) 1 Q ϕ Q ϕ 1 M w ( x ) = M s α 1 L Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) = M s α 1 L ϕ ; w ( x ) ,
which is the stated relation for fractional integrals.
For fractional derivatives, since the linear substitution Q ϕ ( x ) ϕ ( a ) 1 Q ϕ commutes with the d d x operator, the same manipulations as above give rise to:
L ϕ ; w ( x ) D a R L ϕ ( x ) ; w ( x ) α = L D 0 R L x α Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) , L ϕ ; w ( x ) D a C ϕ ( x ) ; w ( x ) α = L D 0 C x α Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) .
For the Riemann–Liouville case, we then have
L ϕ ; w ( x ) D a C ϕ ( x ) ; w ( x ) α f ( x ) = L D a C ϕ ( x ) ; w ( x ) α Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) f ( x ) = s α L { f ( x ) } i = 0 n 1 s n i 1 lim x 0 + I 0 x n i α Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) f ( x ) = s α L { f ( x ) } i = 0 n 1 s n i 1 lim x a + Q ϕ ( x ) ϕ ( a ) I 0 x n i α Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) f ( x ) = s α L { f ( x ) } w ( a + ) i = 0 n 1 s n i 1 I a R L ϕ ( x ) ; w ( x ) n i α f ( a + ) ,
which is the stated relation, and similarly for the Caputo case. □
Definition 9
([33]). The w-weighted ϕ-convolution of two real-valued or complex-valued functions f , g : [ a , ) C is the function f ϕ w g defined by
f ϕ w g ( x ) = 1 w ( x ) a x w ϕ 1 ϕ ( x ) + ϕ ( a ) ϕ ( t ) × f ϕ 1 ϕ ( x ) + ϕ ( a ) ϕ ( t ) w ( t ) g ( t ) ϕ ( t ) d t .
Theorem 7.
The w-weighted ϕ-convolution is related to classical convolution via the following formula:
f ϕ w g = M w ( x ) 1 Q ϕ ( x ) ϕ ( a ) Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) f Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) g ,
or equivalently, using the notation of binary operations and conjugation,
ϕ w = M w ( x ) 1 Q ϕ ( x ) ϕ ( a ) Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) , Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) ,
where the M and Q operators are as in (9) and Lemma 4, respectively.
Proof. 
Clearly, the convolution of Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) f and Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) g can be written as
Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) f Q ϕ ( x ) ϕ ( a ) 1 M w ( x ) g ( x ) = 0 x w ϕ 1 x u + ϕ ( a ) f ϕ 1 x u + ϕ ( a ) × w ϕ 1 u + ϕ ( a ) g ϕ 1 u + ϕ ( a ) d u = a ϕ 1 ( x + ϕ ( a ) ) w ϕ 1 x + 2 ϕ ( a ) ϕ ( t ) × f ϕ 1 x + 2 ϕ ( a ) ϕ ( t ) w ( t ) g ( t ) ϕ ( t ) d t ,
where we substituted u = ϕ ( t ) ϕ ( a ) in the integral. Applying Q ϕ ( x ) ϕ ( a ) to this expression turns it into
a x w ϕ 1 ϕ ( x ) + ϕ ( a ) ϕ ( t ) f ϕ 1 ϕ ( x ) + ϕ ( a ) ϕ ( t ) w ( t ) g ( t ) ϕ ( t ) d t ,
and then dividing by w ( x ) gives exactly the w-weighted ϕ -convolution f ϕ w g ( x ) , as required. □
The result of Theorem 7 for the w-weighted ϕ -convolution, combined with the previous result of Theorem 5 for the w-weighted ϕ -Laplace transform, enable the following result, already seen in [33] (Theorem 5.9), to be proved immediately by composition of operators.
Corollary 9.
If f , g : [ 0 , X ] C are piecewise continuous and of w-weighted ϕ-exponential order c > 0 , then
L ϕ ; w ( x ) f ϕ w g = L ϕ ; w ( x ) { f } L ϕ ; w ( x ) { g } .
Finally, the w-weighted ϕ -Laplace transform can be used to establish a regularity condition for the solutions to weighted fractional differential equations with respect to functions, for example the following initial value problem:
D 0 C ϕ ( x ) ; w ( x ) α y ( x ) = A y ( x ) + g ( x ) , x 0 ,
y ( 0 ) = η ,
where 0 < α < 1 is fixed and A = a i j is an n × n constant matrix and g is a continuous n-dimensional vector-valued function and η is a constant n-dimensional vector.
Theorem 8.
Assume that the system (22) and (23) has a unique continuous solution y . If g is continuous on [ 0 , ) and w-weighted ϕ-exponentially bounded, then y and D 0 C x ; w ( x ) α y are both w-weighted ϕ-exponentially bounded too.
Proof. 
Using the result proved in [44] (Theorem 3.1) together with the conjugation relations given by Theorem 4, we obtain the required result. □

5. Conclusions

This paper has constructed a formal mathematical analysis of the structure of the so-called weighted fractional calculus, and also of the same operators taken with respect to an increasing function. Both of these frameworks can be seen as general classes of operators, and within each general class there is an operational conjugation relation which connects every operator in the class back to the basic Riemann–Liouville and Caputo fractional calculi via multiplication and composition. Although these conjugation relations have been mentioned previously in the literature, some previous research on weighted fractional calculus with respect to functions had failed to take them into account. As well as the new results proved in this paper, one of our main contributions has been to promote the value of the conjugation relations, which make the whole theory more efficient by enabling shorter simpler proofs based on knowledge of the RL and Caputo fractional calculi.
The general class of weighted fractional calculus operators includes such well-known types of operators as those of tempered fractional calculus and Kober–Erdélyi fractional calculus. The even more general class of weighted fractional calculus operators with respect to functions includes operators such as those of Hadamard-type fractional calculus and Erdélyi–Kober fractional calculus, as well as their extensions to be taken with respect to an arbitrary increasing function. Thus, the theory developed in this paper can be useful in the understanding of various types of fractional calculus that have already appeared in the literature and been found with applications worthy of discussion. The whole general class itself, weighted fractional calculus with respect to functions, has been found useful in probability theory and variational calculus.
Corresponding to each fractional differintegral operator in these general classes, there is both a Laplace-type integral transform and a convolution operation. Some of these are defined for the first time here, and we emphasise (which has not been noted before) that these too can be expressed as simple operational modifications of the classical Laplace transform and convolution. It is important to be aware that the modified Laplace transform and convolution do not give rise to any new theory, as their connection with the classical versions allows all results concerning them to be deduced directly. Just like the classical Laplace transform, the weighted Laplace transform (with respect to a function) can be used to solve differential equations with appropriate differential operators, as we have demonstrated briefly in this paper.
The current work is purely theoretical, but much remains to be done in the direction of fractional differential equations and applications. It is presumed that many methods used for solving fractional differential equations, either analytically or numerically, can be extended to corresponding methods which apply to solve weighted fractional differential equations or weighted fractional differential equations with respect to functions. Our current study, particularly the conjugation relations proved herein, will be vital in establishing such extensions.
In the pure mathematical direction, it is possible to extend these general classes of operators still further. The class of fractional differintegrals with analytic kernels and the class of fractional differintegrals with respect to functions have already been combined into a single superclass; the same is done here with the class of weighted fractional differintegrals and the class of fractional differintegrals with respect to functions. Future work may focus on combining the class of weighted differintegrals with the class of differintegrals with analytic kernels, in order to obtain yet another superclass of operators in fractional calculus.

Author Contributions

Conceptualization, A.F.; formal analysis, A.F. and H.M.F.; investigation, A.F. and H.M.F.; writing—original draft preparation, H.M.F.; writing—review and editing, A.F.; supervision, A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Miller, K.S.; Ross, B. An Introduction to the Fractional Calculus and Fractional Differential Equations; John Wiley: New York, NY, USA, 1993. [Google Scholar]
  2. Oldham, K.B.; Spanier, J. The Fractional Calculus; Academic Press: New York, NY, USA; London, UK, 1974. [Google Scholar]
  3. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publishers: New York, NY, USA; London, UK, 1993. [Google Scholar]
  4. Sun, H.; Zhang, Y.; Baleanu, D.; Chen, W.; Chen, Y. A new collection of real world applications of fractional calculus in science and engineering. Commun. Nonlin. Sci. Numer. Simul. 2018, 64, 213–231. [Google Scholar] [CrossRef]
  5. Tarasov, V.E. Review of some promising fractional physical models. Int. J. Mod. Phys. B 2013, 27, 1330005. [Google Scholar] [CrossRef] [Green Version]
  6. Bonfanti, A.; Kaplan, J.L.; Charras, G.; Kabla, A. Fractional viscoelastic models for power-law materials. Soft Matt. 2020, 16, 6002–6020. [Google Scholar] [CrossRef] [PubMed]
  7. Skwara, U.; Martins, J.; Ghaffari, P.; Aguiar, M.; Boto, J.; Stollenwerk, N. Applications of fractional calculus to epidemiological models. In AIP Conference Proceedings; American Institute of Physics: College Park, MD, USA, 2012; Volume 1479, pp. 1339–1342. [Google Scholar]
  8. Acay, B.; Bas, E.; Abdeljawad, T. Fractional economic models based on market equilibrium in the frame of different type kernels. Chaos Solitons Fractals 2020, 130, 109438. [Google Scholar] [CrossRef]
  9. Tarasov, V.E. Mathematical economics: Application of fractional calculus. Mathematics 2020, 8, 660. [Google Scholar] [CrossRef]
  10. Acay, B.; Inc, M. Electrical circuits RC, LC, and RLC under generalized type non-local singular fractional operator. Fractal Fract. 2021, 5, 9. [Google Scholar] [CrossRef]
  11. El-Ajou, A. Adapting the Laplace transform to create solitary solutions for the nonlinear time-fractional dispersive PDEs via a new approach. Eur. Phys. J. Plus 2021, 136, 1–22. [Google Scholar] [CrossRef]
  12. Fernandez, A.; Baleanu, D.; Fokas, A.S. Solving PDEs of fractional order using the unified transform method. Appl. Math. Comput. 2018, 339, 738–749. [Google Scholar] [CrossRef] [Green Version]
  13. Baleanu, D.; Fernandez, A. A generalisation of the Malgrange–Ehrenpreis theorem to find fundamental solutions to fractional PDEs. Elect. J. Qualit. Theory Differ. Equ. 2017, 15, 1–12. [Google Scholar] [CrossRef]
  14. Fernandez, A. An elliptic regularity theorem for fractional partial differential operators. Comput. Appl. Math. 2018, 37, 5542–5553. [Google Scholar] [CrossRef] [Green Version]
  15. Williams, P. Fractional Calculus of Schwartz Distributions. Bachelor of Science Thesis, University of Melbourne, Melbourne, Australia, 2007. [Google Scholar]
  16. Djida, J.D.; Nieto, J.J.; Area, I. Nonlocal time-porous medium equation: Weak solutions and finite speed of propagation. Disc. Contin. Dynam. Syst. B 2019, 24, 4031–4053. [Google Scholar] [CrossRef] [Green Version]
  17. Djida, J.D.; Nieto, J.J.; Area, I. Nonlocal time porous medium equation with fractional time derivative. Rev. MatemáTica Complut. 2019, 32, 273–304. [Google Scholar] [CrossRef] [Green Version]
  18. Guo, B.; Pu, X.; Huang, F. Fractional Partial Differential Equations and Their Numerical Solutions; World Scientific: Singapore, 2011. [Google Scholar]
  19. Li, C.; Chen, A. Numerical methods for fractional partial differential equations. Int. J. Comput. Math. 2018, 95, 1048–1099. [Google Scholar] [CrossRef]
  20. Momani, S.; Odibat, Z.; Hashim, I. Algorithms for nonlinear fractional partial differential equations: A selection of numerical methods. Topol. Meth. Nonlinear Anal. 2008, 31, 211–226. [Google Scholar]
  21. Hilfer, R.; Luchko, Y. Desiderata for Fractional Derivatives and Integrals. Mathematics 2019, 7, 149. [Google Scholar] [CrossRef] [Green Version]
  22. Teodoro, G.S.; Machado, J.T.; de Oliveira, E.C. A review of definitions of fractional derivatives and other operators. J. Comput. Phys. 2019, 388, 195–208. [Google Scholar] [CrossRef]
  23. Baleanu, D.; Fernandez, A. On fractional operators and their classifications. Mathematics 2019, 7, 830. [Google Scholar] [CrossRef] [Green Version]
  24. Fernandez, A.; Özarslan, M.A.; Baleanu, D. On fractional calculus with general analytic kernels. Appl. Math. Comput. 2019, 354, 248–265. [Google Scholar] [CrossRef] [Green Version]
  25. Osler, T.J. Leibniz rule for fractional derivatives generalized and an application to infinite series. SIAM J. Appl. Math. 1970, 18, 658–674. [Google Scholar] [CrossRef]
  26. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  27. Oumarou, C.M.S.; Fahad, H.M.; Djida, J.D.; Fernandez, A. On fractional calculus with analytic kernels with respect to functions. Comput. Appl. Math. 2021, 40, 244. [Google Scholar] [CrossRef]
  28. Agrawal, O.P. Some generalized fractional calculus operators and their applications in integral equations. Fract. Calc. Appl. Anal. 2012, 15, 700–711. [Google Scholar] [CrossRef]
  29. Agrawal, O.P. Generalized multiparameters fractional variational calculus. Int. J. Differ. Equ. 2012, 2012, 521750. [Google Scholar] [CrossRef] [Green Version]
  30. Kolokoltsov, V.N. The probabilistic point of view on the generalized fractional partial differential equations. Fract. Calc. Appl. Anal. 2019, 22, 543–600. [Google Scholar] [CrossRef] [Green Version]
  31. Al-Refai, M. On weighted Atangana–Baleanu fractional operators. Adv. Differ. Equ. 2020, 2020, 3. [Google Scholar] [CrossRef] [Green Version]
  32. Al-Refai, M.; Jarrah, M. Fundamental results on weighted Caputo–Fabrizio fractional derivative. Chaos Solitons Fractals 2019, 126, 7–11. [Google Scholar] [CrossRef]
  33. Jarad, F.; Abdeljawad, T.; Shah, K. On the weighted fractional operators of a function with respect to another function. Fractals 2020, 28, 2040011. [Google Scholar] [CrossRef]
  34. Abdo, M.S.; Abdeljawad, T.; Ali, M.; Shah, K.; Jarad, F. Existence of positive solutions for weighted fractional order differential equations. Chaos Solitons Fractals 2020, 141, 110341. [Google Scholar] [CrossRef]
  35. Bayrak, M.A.; Demir, A.; Özbilge, E. On solution of fractional partial differential equation by the weighted fractional operator. Alex. Eng. J. 2020, 59, 4805–4819. [Google Scholar] [CrossRef]
  36. Liu, J.G.; Yang, X.J.; Feng, Y.Y.; Geng, L.L. Fundamental results to the weighted Caputo-type differential operator. Appl. Math. Lett. 2021, 121, 107421. [Google Scholar] [CrossRef]
  37. Diethelm, K. The Analysis of Fractional Differential Equations: An Application-Oriented Exposition Using Differential Operators of Caputo Type; Springer: Berlin/Heidelberg, Germany, 2010. [Google Scholar]
  38. Zaky, M.A.; Hendy, A.S.; Suragan, D. A note on a class of Caputo fractional differential equations with respect to another function. Math. Comput. Simul. 2022, 196, 289–295. [Google Scholar] [CrossRef]
  39. Li, C.; Deng, W.; Zhao, L. Well-posedness and numerical algorithm for the tempered fractional ordinary differential equations. Discrete Contin. Dyn. Syst. Ser. B 2019, 24, 1989–2015. [Google Scholar]
  40. Meerschaert, M.M.; Sabzikar, F.; Chen, J. Tempered fractional calculus. J. Comput. Phys. 2015, 293, 14–28. [Google Scholar] [PubMed] [Green Version]
  41. Fernandez, A.; Ustaoğlu, C. On some analytic properties of tempered fractional calculus. J. Comput. Appl. Math. 2020, 366, 112400. [Google Scholar] [CrossRef] [Green Version]
  42. Fahad, H.M.; Ur Rehman, M.; Fernandez, A. On Laplace transforms with respect to functions and their applications to fractional differential equations. Math. Methods Appl. Sci. 2021, 1–20. [Google Scholar] [CrossRef]
  43. Jarad, F.; Abdeljawad, T. Generalized fractional derivatives and Laplace transform. Discret. Contin. Dyn. Syst.-S 2019, 13, 709–722. [Google Scholar] [CrossRef] [Green Version]
  44. Li, K.; Peng, J. Laplace transform and fractional differential equations. Appl. Math. Lett. 2011, 24, 2019–2023. [Google Scholar]
  45. Fahad, H.M.; Fernandez, A.; ur Rehman, M.; Siddiqi, M. Tempered and Hadamard-type fractional calculus with respect to functions. Mediterr. J. Math. 2021, 18, 143. [Google Scholar] [CrossRef]
  46. Butzer, P.L.; Kilbas, A.A.; Trujillo, J.J. Fractional calculus in the Mellin setting and Hadamard-type fractional integrals. J. Math. Anal. Appl. 2002, 269, 1–27. [Google Scholar] [CrossRef] [Green Version]
  47. Kilbas, A.A. Hadamard-Type Fractional Calculus. J. Korean Math. Soc. 2001, 38, 1191–1204. [Google Scholar]
  48. Odibat, Z.; Baleanu, D. On a New Modification of the Erdélyi–Kober Fractional Derivative. Fractal Fract. 2021, 5, 121. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Fernandez, A.; Fahad, H.M. Weighted Fractional Calculus: A General Class of Operators. Fractal Fract. 2022, 6, 208. https://doi.org/10.3390/fractalfract6040208

AMA Style

Fernandez A, Fahad HM. Weighted Fractional Calculus: A General Class of Operators. Fractal and Fractional. 2022; 6(4):208. https://doi.org/10.3390/fractalfract6040208

Chicago/Turabian Style

Fernandez, Arran, and Hafiz Muhammad Fahad. 2022. "Weighted Fractional Calculus: A General Class of Operators" Fractal and Fractional 6, no. 4: 208. https://doi.org/10.3390/fractalfract6040208

APA Style

Fernandez, A., & Fahad, H. M. (2022). Weighted Fractional Calculus: A General Class of Operators. Fractal and Fractional, 6(4), 208. https://doi.org/10.3390/fractalfract6040208

Article Metrics

Back to TopTop