Next Article in Journal
Some Estimates of k-Fractional Integrals for Various Kinds of Exponentially Convex Functions
Next Article in Special Issue
Boundary Value Problem for Impulsive Delay Fractional Differential Equations with Several Generalized Proportional Caputo Fractional Derivatives
Previous Article in Journal
Fekete–Szegö Problem and Second Hankel Determinant for a Class of Bi-Univalent Functions Involving Euler Polynomials
Previous Article in Special Issue
Extended Stability and Control Strategies for Impulsive and Fractional Neural Networks: A Review of the Recent Results
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Fractional Scale Calculus: Hadamard vs. Liouville

by
Manuel D. Ortigueira
1,* and
Gary W. Bohannan
2
1
Centre of Technology and Systems–UNINOVA and Department of Electrical Engineering, NOVA School of Science and Technology, NOVA University of Lisbon, Quinta da Torre, 2829-516 Caparica, Portugal
2
9248 Liggon Green Lane, Germantown, TN 38139, USA
*
Author to whom correspondence should be addressed.
Fractal Fract. 2023, 7(4), 296; https://doi.org/10.3390/fractalfract7040296
Submission received: 7 February 2023 / Revised: 23 March 2023 / Accepted: 27 March 2023 / Published: 29 March 2023
(This article belongs to the Special Issue Feature Papers for the 'General Mathematics, Analysis' Section)

Abstract

:
A general fractional scale derivative is introduced and studied. Its relation with the Hadamard derivatives is established and reformulated. A new derivative similar to the Grünwald–Letnikov’s is deduced. Tempered versions are also introduced. Scale-invariant systems are described and exemplified. For solving the corresponding differential equations, a new logarithmic Mittag-Leffler series is proposed.

1. Introduction

Scale is a term we find often but that we do not know how to define. L. Cohen [1] considered it as a physical attribute similar to frequency; he proposed a scale operator for representing it and introduced a scale transform. However, L. Nottale [2] defined scale as the resolution with which measurements are performed and considered it as a relative state of the reference system. He paid attention mainly to the scale transformations he obtained with the operators based on the scale derivative obtained from fractal considerations [3,4]. A similar procedure was developed by J. Cresson [5,6]. For these approaches, the important factor is the scale invariance [7,8], which led to establishing bridges to other formulations such as the Lamperti transformations for stochastic processes [9,10] and the operator theory [11]. However, the importance of scale was fully recognized in the multiscale/multiresolution analysis introduced during the eighties of the last century in signal processing, mainly in image processing, obtained from wavelet transform [12,13,14,15,16]. In applications, the discrete version is normally used [17,18,19]. However, this approach goes further in the sense that it simultaneously explores the shift and scale characteristics of the wavelet used to define the transform such that it keeps the form but changes the position and scale. Contrary to this route, we are interested in the study of tools and systems that are scale invariant. It is not surprising that the Mellin transform has appeared in this context (see J. Bertrand et al. [20]) and, as consequence, in fractional calculus through the Hadamard derivative [21,22,23,24].
In this paper, we join Mellin transform and Hadamard’s idea to introduce a general framework for the scale fractional derivative. We follow a path identical to the one proposed by Liouville for the shift-invariant fractional derivatives [25,26].
The first reference to the possibility of non-integer order derivatives was found in a letter from G. Leibniz to J. Bernoulli [27]. Although Euler (1730), Fourier (1822), and Abel (1823) touched on the problem, we must consider Liouville as the true father of fractional calculus (FC) through two papers he published in 1832 [25,26]. His approach was based on generalizing the formula for the derivative of the exponential. This brought many difficulties in imposing his vision, since, at that time, the inverse Laplace integral was unknown; so, Liouville could not find a simple way of expressing a function in terms of exponentials. The main definitions of fractional derivatives (FD) we find nowadays are based on the formulæ presented by Liouville, mainly the Riemann–Liouville [22,28,29], (Dzherbashian)–Caputo [28,29,30], and Grünwald–Letnikov [22,26,31] definitions. However, and based on these derivatives, other ones were proposed. Among the most important, we refer to Hadamard’s [21,28] and Marchaud’s [22].
It is important to remark that Liouville’s approach was based on the fact that the order α fractional derivative, D α , of
f ( t ) = n 0 a n e γ n t ,
with a n , γ n R , would be given by
D α f ( t ) = n 0 a n γ n α e γ n t ,
Hadamard [21] arrived at FC through Riemann’s formulation [32] and had in mind the search for a derivative, D , such that if the function at hand has a MacLaurin expansion,
g ( t ) = n 0 a n t n ,
then
D g ( t ) = n 0 a n n α t n .
With a straightforward modification to the Riemann formulation followed by a logarithmic transformation, he could devise the first formula for what is called the Hadamard integral. If α > 0 , it reads
D α f ( t ) = 1 Γ ( α ) 0 1 ln ( t ) α 1 f ( t η ) d η η
It is interesting to remark that, in passing, he introduced the operator t α D α f ( t ) that we will recover later; however, he did not formally present the derivative expression. This was achieved, for the first time, by Samko et al. [22].
The framework we will propose starts from two different types of linear systems, shift and dilation (scale)-invariant, such that we establish a parallel way into the Liouville and Hadamard derivatives, having their eigenfunctions as “seeds”. The derivatives are defined from the corresponding eigenvalues by Laplace and Mellin inverse transforms. The shift-invariant systems led to the Grünwald–Letnikov (GL) and Liouville (L) derivatives [33]. Similar procedures here will lead to the corresponding Hadamard and new GL-like scale derivatives. The formulation clarifies the difference between left and right derivatives through their regions of convergence: right or left complex half-planes.
The paper is outlined as follows. In Section 2, we make a brief introduction to the linear systems. In Section 3, the derivatives based on the Liouville formulation and referred in [33,34,35] are then obtained sequentially; these result from the shift-invariant systems. Similar development is achieved from the scale-invariant systems in Section 4, where we obtain their eigenfunctions and corresponding eigenvalue that is the base for the definition of scale derivative. In passing, we established its relation with the classic and quantum derivatives and deduced a logarithmic series analog to the Taylor series. The properties of the scale derivative allowed us to reformulate the Hadamard derivatives that we generalized to obtain the tempered scale derivatives. Finally, we conclude the paper with a reflection on the results we obtained.
Remark 1.
The questions associated with the existence and uniqueness of the fractional derivatives are very important. They have been treated in many good texts: [22,28,30] for Liouville type derivatives, [22,31] for the GL, and [22,28,36,37,38,39,40,41,42,43] for the Hadamard derivatives. Therefore, we will not pay attention to the subject. We will assume always to work with functions having Laplace or Mellin transform.

2. On the Linear Systems

In a simple way, we can define a system as a combination of components (mathematical, physical, biological, social, etc.) acting together to fulfill certain predefined goals. Interaction with systems is performed by means of functions called signals [44]. Mathematically, a system is defined as an application on the set of signals—that is, a transformation of one signal, x ( t ) , into another, y ( t ) . Let T [ . ] be an operator that symbolically represents such a transformation,
y ( t ) = T x ( t ) ,
where x ( t ) is the input or excitation and y ( t ) is the output or response. A system is linear if it is formally represented by a linear operator, which means that it satisfies the properties of additivity and homogeneity (superposition principle) [44,45]
y ( t ) = T a 1 x 1 ( t ) + a 2 x 2 ( t ) = a 1 T x 1 ( t ) + a 2 T x 2 ( t ) = a 1 y 1 ( t ) + a 2 y 2 ( t )
There are many classes of linear systems (LS). Among them, we consider the shift-invariant and dilation (scale)-invariant.
Definition 1.
We call a linear system shift-invariant (SI) if its input–output relation is given by the usual convolution:
y ( t ) = x ( t ) g ( t ) = + x ( t η ) g ( η ) d η ,
where t R and g ( t ) is its impulse response (or Green function): the response to x ( t ) = δ ( t ) .
Definition 2.
We call a linear system scale-invariant or dilation-invariant (DI) if its input–output relation is given by the Mellin convolution
y ( τ ) = x ( τ ) g ( τ ) = 0 x τ η g ( η ) d η η ,
where τ R + and g ( τ ) is the impulse response: the response to x ( τ ) = δ ( τ 1 ) .
We demand that the impulse response, g ( t ) , be at least
  • piecewise continuous,
  • with bounded variation.
The specific applications impose constraints on the signals. However, we will assume they are square integrable: “signals type energy” [46]. For “enough good” functions, the derivative (suitably defined) of a convolution is equal to the convolution of one function with the derivative of the other,
D α x ( t ) y ( t ) = x ( t ) D α y ( t ) .
Remark 2.
It is important to emphasize that, if D represents a SI derivative, from a numerical point of view, D α ( x y ) is preferable relative to D α x y or x D α y . For the Mellin convolution, the situation is similar.

3. Shift-Invariant Systems: The Liouville Derivatives

Consider a linear shift-invariant system (LTIS) defined by (4). Then, the following apply [44]:
1.
The exponentials are the eigenfunctions of LTIS
y ( t ) = G ( s ) e s t , t R ,
with eigenvalues given by G ( s ) given by
G ( s ) = + g ( t ) e s t d t
which is the transfer function and the (bilateral) Laplace transform (LT) of the impulse response of the system.
2.
If the region of convergence (ROC) of G ( s ) contains the imaginary axis, we can set s = i ω , ω R , in such a way that the response of an LS to a sinusoid is also a sinusoid with the same frequency. In such a situation, the LT degenerates into the Fourier transform and we say that the system is stable.
Definition 3.
Let α R . The α-order SI derivative, D α , is an operator such that
D α e s t = s α e s t , t R , s C ,
for R e ( s ) > 0 (causal case) or R e ( s ) < 0 (anti-causal case). We can extend the validity of this definition to R e ( s ) = 0 if s 0 .
This definition dates back to Liouville who introduced it in the first paper on fractional calculus [25].
The inverse LT theorem (Bromwich integral) [47] allows us to state that if x ( t ) is a function with LT, X ( s ) , then x ( t ) has a derivative in the Liouville sense and
D α x ( t ) = 1 2 π i γ s α X ( s ) e s t d s ,
where γ is a vertical straight line located in the ROC of X ( s ) .
Remark 3.
In the following, we shall be working with derivative orders that are always real numbers. The results remain valid if they are complex numbers, but the resulting derivatives are not hermitian operators, implying that they are not useful in applications [48].
Liouville’s definition highlights the importance of the transfer function, G ( s ) = s α . The way how we express it in the inverse time domain leads to several expressions for the derivative. To start, we can obtain direct generalizations of classic derivatives by noting that
s α = lim h 0 + ( 1 e s h ) h α R e ( s ) > 0 , ( e s h 1 ) h α R e ( s ) < 0 ,
that lead to the forward and backward Grünwald–Letnikov (GL) derivatives, respectively [22,31]. However, the first proposal of such a derivative was also presented by Liouville [26]. Let ( a ) n , n = 1 , 2 , define the Pochhammer symbol for the rising factorial
( a ) 0 = 1 , ( a ) n = k = 0 n 1 ( a + k ) .
As
( 1 e s h ) α = n = 0 ( α ) n n ! e n s h
it is a simple matter to obtain the GL derivatives using the translation property of the LT. For the causal (forward) case, we have
D + α f ( t ) : = lim h 0 + h α n = 0 + α n n ! f ( t n h ) .
The backward is obtained through the substitution h h .
However, we may be more interested here in obtaining the convolutional versions of the SI derivatives. As is well known, the impulse response, g ( t ) , of the SI system with transfer function s α is given by [47]
g ( t ) = ± t α 1 Γ ( α ) ε ( ± t )
where ε ( t ) is the Heaviside unit step; the + sign corresponds to the causal case ( R e ( s ) > 0 ) and − to the anti-causal ( R e ( s ) < 0 ) [31]. Consequently, attending to (7) and to the convolution property of the LT, we conclude that
D ± α x ( t ) = x ( t ) g ( t ) = ± 1 Γ ( α ) 0 η α 1 f ( t η ) d η ,
The formula corresponding to the − sign was introduced first by Liouville [25] who noted that the derivative case ( α > 0 ) leads to a singular convolution. However, instead of a regularization, he devised a trick for solving the singularity problem that we can describe as
s α = s N s α N = s α N s N ,
where N > α . The most usual choice is the first integer greater than α , N 1 < α N , that we denote by N = α . Basically, it consists of transferring the singular behavior to an integer order derivative. The first approach leads to the so-called (Riemann–)Liouville derivatives [22],
D + α f ( t ) = 1 Γ ( N α ) d N d t N 0 f ( t η ) η N α 1 d η = 1 Γ ( N α ) d N d t N t f ( η ) ( t η ) N α 1 d η
and
D α f ( t ) = ( 1 ) N Γ ( N α ) d N d t N 0 f ( t + η ) η N α 1 d η = ( 1 ) N Γ ( N α ) d N d t N t f ( η ) ( t η ) N α 1 d η .
Liouville’s second procedure leads to what can be called Liouville–Caputo derivatives [44,49] that read
D + α f ( t ) = 1 Γ ( N α ) t f ( N ) ( η ) ( t η ) N α 1 d η
and
D α f ( t ) = ( 1 ) N Γ ( N α ) t f ( N ) ( η ) ( t η ) N α 1 d η .
Remark 4.
The usually called Riemann–Liouville and (Dzherbashian–)Caputo derivatives are particular cases of these, obtained for functions defined in intervals of R .
The problem raised by the singularity of the integral (12) can be solved in an alternative way through a regularization [31,44]. For the forward case, it reads
D + α f ( t ) = 1 Γ ( α ) 0 η α 1 f ( t η ) 0 N 1 ( ) m f ( m ) ( t ) m ! η m d η ,
which we will call the regularized forward Liouville derivative. The backward case is similar
D + α f ( t ) = 1 Γ ( α ) 0 η α 1 f ( t + η ) 0 N 1 f ( m ) ( t ) m ! η m d η ,
Remark 5.
The method introduced in (13) can be generalized. For example,
s α = s N k s α N s k , α R + , N Z 0 + ,
with k < N and N > α . This procedure was used by Davidson, Essex, and Cavanati [35]. Similarly, Hilfer [50] proposed the decomposition
s α = s μ ( 1 α ) s s ( 1 μ ) ( 1 α ) ,
where 0 < α , μ < 1 . This approach can be generalized and can lead to the invention of many alternative derivatives [33].
The derivatives that we described in this section are suitable for modeling SI linear systems. However, they are useful in dealing with nonlinear systems too [51].

4. On the Scale-Invariant Systems: Hadamard Derivatives

4.1. From the System to the Derivative

In this section, we are going to reproduce the approach used in the SI case to the DI. We start from (5), which we repeat here
y ( τ ) = x ( τ ) g ( τ ) = 0 x τ η g ( η ) d η η .
We must remark that the procedure in the following only makes sense if the working domain is R + . Therefore, we will assume it implicitly in the following.
Similarly to the SI case, the powers x ( τ ) = τ v , τ R + , v C are the eigenfunctions of the DI systems. In fact, if the input is x ( τ ) = τ v , then the output y ( τ ) = x ( τ ) g ( τ ) is
y ( τ ) = G ( v ) τ v
where G ( v ) is the transfer function given by
G ( v ) = 0 g ( u ) u v 1 d u ,
which is the Mellin transform (MT) of the impulse response.
The MT in (23) has a parameter sign change v v relatively to the usual MT [20,52,53]. However, we will keep this definition, since it establishes a better parallelism with the LT, which concerns the region of convergence. The corresponding inverse Mellin transform is given by the integral
x ( τ ) = M 1 X ( v ) = 1 2 π i γ X ( v ) τ v d v , τ R +
where γ is a vertical straight line in the ROC of the transform.
Remark 6.
As in the SI case, if the ROC of the transfer function includes the imaginary axis, we will say that the system is stable. This corresponds to write
0 | g ( u ) | u 1 d u < .
Consider the relation (22). We are going to use it as a starting point to obtain the Hadamard derivatives and two new ones, analog to the GL derivatives.
Definition 4.
Let α R . We define the α-order scale derivative (SD) as the operator D s obeying the rule
D s α τ v = v α τ v , τ R + , v C ,
for R e ( v ) > 0 (expansion case) or R e ( v ) < 0 (shrinkage case).
If a function, x ( τ ) , has MT, X ( v ) , then it has a fractional scale-derivative that is given by
D s α x ( τ ) = 1 2 π i γ v α X ( v ) τ v d v , t R +
where γ is a vertical straight line located in the ROC of X ( v ) . Similarly to the shift-invariant case, (25) shows that the scale derivative is an elemental system with transfer function v α . Again, the way we express it in the inverse scale domain leads to various expressions for the derivative.
Theorem 1.
Let q > 1 . The following expressions,
D s + α x ( τ ) = lim q 1 + ln α ( q ) n = 0 ( α ) n n ! x ( τ q n ) ,
and
D s α x ( τ ) = lim q 1 + ( 1 ) α ln α ( q ) n = 0 ( α ) n n ! x ( τ q n ) ,
represent scale-derivatives that we can call stretching and shrinking GL-type derivatives, respectively.
Proof. 
As it is not hard to deduce,
v α = lim q 1 + ( 1 q v ) ln q α R e ( v ) > 0 , ( q v 1 ) ln q α R e ( v ) < 0 ,
Using the binomial theorem,
( 1 q ± v ) α = n = 0 ( α ) n n ! q ± n v ,
and the dilation property of the MT,
M x ( a τ ) = a v X ( v ) a R + ,
we obtain (27) and (28). □
We can extend these results to R e ( v ) = 0 , v 0 .
Particular cases are interesting. Let α = ± 1 .
1.
From (27), we obtain
D s + x ( τ ) = lim q 1 + x ( τ ) x ( τ q 1 ) ln q
and
D s + 1 x ( τ ) = lim q 1 + ln q n = 0 x ( τ q n ) ,
2.
Ref. (28) gives
D s x ( τ ) = lim q 1 + x ( τ q ) x ( t ) ln q
and
D s 1 x ( τ ) = lim q 1 + ln q n = 0 x ( τ q n ) ,
Example 1.
1.
Power functions: τ a
We have
D s + τ a = lim q 1 + τ a τ a q a ln q = τ a lim q 1 + 1 1 a ln ( q ) + a 2 a ln 2 ( q ) / 2 ln q = a τ a
and
D s τ a = lim q 1 + τ a q a τ a ln q = τ a lim q 1 + q a 1 ln q = a τ a
2.
Logarithm: ln a ( τ )
As above, we obtain
D s + ln a ( τ ) = lim q 1 + ln a ( τ ) ln ( τ ) ln ( q ) a ln ( q ) = ln a ( τ ) lim q 1 + a ln ( q ) ln ( τ ) a 2 2 ln ( q ) ln ( τ ) 2 + ln ( q ) = a ln a 1 ( τ )
and
D s ln a ( τ ) = lim q 1 + ln ( τ ) + ln ( q ) a ln a ( τ ) ln ( q ) = ln a ( τ ) lim q 1 + a ln ( q ) ln ( τ ) + a 2 2 ln ( q ) ln ( τ ) 2 + ln ( q ) = a ln a 1 ( τ )

4.2. Properties of the Scale Derivatives

Definition 4 together with (26), which serve as base for the introduction of the scale derivatives, allow us to obtain their corresponding properties that are similar to those of the Liouville derivatives [31]. Let α , β R . The scale derivative enjoys the following properties:
  • Linearity
    It is obvious from (26).
  • Additivity and Commutativity of the orders
    D s α D s β x ( τ ) = D s α + β x ( τ ) .
    This comes from (25).
  • Neutral and inverse elements
    Let α = β . Then,
    D s α D s α x ( τ ) = D s 0 x ( τ ) = x ( τ ) .
    From (36), we conclude that there is always an inverse element—that is, for every α there is always the α order that we call anti-derivative.
  • The generalized Leibniz rule
    This rule gives the FD of the product of two functions and assumes the format of other fractional derivatives [31]
    D s α x ( τ ) y ( τ ) = k = 0 α k D s k x ( τ ) D s α k y ( τ ) .
    To prove this relation, we note first that
    M x ( τ ) y ( τ ) = X ( v ) Y ( v ) .
    Using the Bromwich inverse Mellin transform, we can write
    D s α x ( τ ) y ( τ ) = 1 2 π i γ 1 v α γ 2 X ( u ) Y ( v u ) d u τ v d v ,
    where γ 1 and γ 2 are vertical straight lines in the intersection of the region of convergence of both transforms. With a trick, with which we must be careful,
    v α = ( v u + u ) α = ( v u ) α 1 + u v u α = k = 0 α k k k ( v u ) α k ,
    and evident manipulations, we obtain
    D s α x ( τ ) y ( τ ) = k = 0 α k 1 2 π i γ 1 γ 2 u k X ( u ) ( v u ) α k Y ( v u ) d u τ v d v ,
    from where the property (37) results.
These properties refer to derivatives defined on R + . In mathematical texts, it is usual to restrict the definitions by including the domain of the functions that are assumed to have bounded support. In such situations, the properties are not exactly as those we presented [37].

4.3. Relation with Classic and Quantum Derivatives

Consider derivative (31). Assume that we apply the l’Hôpital rule for the indeterminacy. We obtain
D s + x ( τ ) = lim q 1 + x ( τ ) x ( τ q 1 ) ln q = lim q 1 + τ x ( τ q 1 ) 1 / q = τ x ( τ ) ,
which expresses the usual way of representing D s + x ( τ ) [22]. Therefore, the scale derivative of order n, D s + n is the result of applying τ x ( τ ) n times. Such interpretation is not so obvious in the fractional case. Anyway, we can write D s + α = τ D α . As for D s x ( τ ) , the conclusion is the same.
These results suggest we go into the usual derivative definition
x ( t ) = lim h 0 + x ( t ) x ( t h ) h
and make the substitution t h t q 1 to obtain
x ( t ) = lim q 1 + x ( t ) x ( t q 1 ) 1 q 1 t .
Similarly, with
x ( t ) = lim h 0 + x ( t + h ) x ( t ) h
and the substitution t + h t q , we obtain
x ( t ) = lim q 1 + x ( q t ) x ( t ) q 1 t .
The Formulae (38) and (39) express the so-called “quantum derivatives”, frequently written without the limit operation [54,55]. These results lead to other formulations for the scale derivative alternative to (31) and (32). We may write
D q + x ( τ ) = lim q 1 + x ( τ ) x ( τ q 1 ) 1 q 1 .
It must be noted that D q + x ( τ ) = D s + x ( τ ) , since lim q 1 + ln q = lim q 1 + 1 q 1 = 0 . Similarly, we define
D q x ( τ ) = lim q 1 + x ( q τ ) x ( τ ) q 1 .
The fractionalization of derivatives (38) and (39) was performed in [56]. For the first case, we obtained a Grünwald–Letnikov like fractional quantum derivative
D q + α x ( τ ) = τ α lim q 1 + j = 0 α j q 1 ( 1 ) j q j ( j + 1 ) / 2 q j α x q j τ ( 1 q 1 ) α
where the q-binomial coefficients are given by
α i q = [ α ] q ! [ j ] q ! [ α i ] q !
and
[ α ] q = 1 q α 1 q
Identically, we obtain
D q α x ( τ ) = τ α lim q 1 + j = 0 α j q 1 ( 1 ) j q j ( j 1 ) / 2 x q j τ ( q 1 ) α
Formulas (42) and (43) represent derivatives that are essentially Liouville derivatives. From them, we obtain new scale derivatives through
D q ± α x ( τ ) = τ α D q ± α x ( τ ) .
However, these are different from those we presented above, (27) and (28). Its usefulness is in solving the Euler–Cauchy differential equations [57].

4.4. Scale Conversion: Logarithmic Series

The previous results show that the logarithm here plays the role performed by the power in the Liouville derivatives. Therefore, we expect to obtain a series analog to Taylor/MacLaurin’s.
On the other hand, we considered a DI system, but nothing was said about the scale modification. With the tools developed above, we can obtain a scale conversion.
Theorem 2.
Let f ( τ ) be an indefinitely scale differentiable bounded function and a > 0 . Then,
f ( a τ ) = n = 0 D s n f ( a ) n ! ln n ( τ )
Proof. 
Computing the MT, (23), of f ( a τ ) relative to a (the two variables play the same role), it is not hard to show that
M f ( a τ ) ( v ) = τ v F ( v ) .
Therefore, there is an operator with a transfer function τ v that performs a scale conversion. We can write
τ v = n = 0 v n n ! ln n ( τ ) ,
which leads immediately to the result. We can obtain another justification by performing successive scale derivations using (33) followed by τ 1 . We must note that D s n ln n ( τ ) = n ln n 1 ( τ ) .
The convergence of the series can be stated by the traditional methods used in the study of power series.
Expression (44) shows how ideally we can obtain a scale conversion.
Example 2.
Let N Z + and f ( τ ) = τ N . According to Example 1,
D s n f ( a ) = N n f ( a ) = N n a N
and
f ( a τ ) = a N n = 0 N n n ! ln n ( τ ) = a N e N ln ( τ ) = a N τ N
as expected.
Corollary 1.
From (44), we obtain
f ( a / τ ) = n = 0 ( 1 ) n D s + n f ( a ) n ! ln n ( τ )
The proof is obvious.
Relations (44) and (46) are valid for both D s + n and D s n The logarithmic series will be useful in the next sub-section to regularize a fractional derivative.

4.5. Hadamard Derivatives

In Section 4.1, we presented two GL-like scale derivatives. If we apply the MT to (27) and (28), and attend to (26), we are led to the generalization of a well-known property:
M D s ± α x ( τ ) = v α X ( v ) , ± R e ( v ) > 0 .
This result means that there should exist two, g ± ( t ) , MT inverses of v α , one for each ROC.
Theorem 3.
Let us assume that α < 0 . The MT inverse of v α is given by
M 1 v α ( τ ) = ln α 1 ( τ ) Γ ( α ) ε ( τ 1 ) R e ( v ) > 0 , ln α 1 ( 1 / τ ) Γ ( α ) ε ( 1 τ ) R e ( v ) < 0 .
Proof. 
We can write for any β < 0 ,
1 v β = 0 η β 1 Γ ( β ) e v η d η = 1 ln β 1 ( τ ) Γ ( β ) τ v 1 d τ ,
from where the first expression results. The second is obtained with the substitution τ 1 / τ . □
These MT inverses allow us to obtain two new scale anti-derivatives.
Corollary 2.
Let α < 0 and τ R + . The relations stated in (47) and (48) used in the Mellin convolution lead to the scale anti-derivatives given by
D s + α x ( τ ) = 1 Γ ( α ) 1 x ( τ / η ) ln α 1 ( η ) d η η = 1 Γ ( α ) 0 τ x ( u ) ln α 1 ( τ / u ) d u u
and
D s α x ( τ ) = 1 Γ ( α ) 0 1 x ( τ / η ) ln α 1 ( 1 / η ) d η η = 1 Γ ( α ) τ x ( u ) ln α 1 ( u / τ ) d u u
These anti-derivatives are the Hadamard integrals ([22,28]).
Theorem 4.
As in the shift-invariant cases, these integrals become singular when the derivative order is positive. We can proceed as in (13) to obtain the following:
1.
Hadamard right derivative [22,28]
D s + α x ( τ ) = 1 Γ ( N α ) D s + N 1 x ( τ / η ) ln N α 1 ( η ) d η η
2.
Hadamard left derivative [22,28]
D s α x ( τ ) = 1 Γ ( N α ) D s N 0 1 x ( τ / η ) ln N α 1 ( 1 / η ) d η η
3.
Hadamard–Liouville right derivative
D s + α x ( τ ) = 1 Γ ( N α ) 1 D s + N x ( τ ) ln N α 1 ( τ / u ) d u u
4.
Hadamard–Liouville left derivative
D s α x ( τ ) = 1 Γ ( N α ) 0 1 D s N x ( τ ) ln N α 1 ( u / τ ) d u u
To prove these results, it is enough to note that
v α = v N v α N = v α N v N , α N Z 0 + .
Remark 7.
The designation “Hadamard–Liouville” is more correct than “Hadamard–Caputo” or “Caputo–Hadamard” [58,59,60], since the idea of performing an integration after a derivation was proposed by Liouville [25]. On the other hand, M. Caputo did not publish anything about this derivative.
Similarly to the procedure used in (18), we can regularize the integrals in (49) and (50) (see [61] for an alternative). To do it, we use the logarithmic series (44) and (46).
Definition 5.
We define regularized Hadamard derivatives by
D s + α x ( τ ) = 1 Γ ( α ) 0 τ x ( u ) n = 0 N 1 ( 1 ) n D s + N x ( τ ) n ! ln n ( τ / u ) ln α 1 ( τ / u ) d u u
and
D s α x ( τ ) = 1 Γ ( α ) τ x ( u ) n = 0 N 1 D s + N x ( τ ) n ! ln n ( u / τ ) ln α 1 ( u / τ ) d u u
We can consider these two expressions valid for any real order, provided that we assume the summation to be null for N 0 .

4.6. Tempered Scale-Invariant Derivatives

The tempered scale-invariant derivatives were already studied in the context of a generalization of the Hadamard derivatives [36,37,62] but without referring directly to the concept. In fact, the SI tempered fractional derivatives are obtained from the SI Liouville derivatives, described previously (Section 3), through a translation in their transfer function. As shown there, the transfer function of the Liouville derivatives is given by G ( s ) = s α with ± R e ( s ) > 0 , according to the causality. Let λ R 0 + . The tempered Liouville derivatives have transfer functions defined by [63]
G ( s ) = ( s + λ ) α ± ( R e ( s ) + λ ) > 0 .
The corresponding impulse functions are given by
L 1 ( s + λ ) α = ± e λ t t α 1 Γ ( α ) ε ( ± t ) .
With (57) and (58), and proceeding as in Section 3, we obtain tempered versions of the Liouville derivatives. The causal derivatives (the anti-causal are readily obtained) read [63]
1.
Forward Grünwald–Letnikov
D λ , f α f ( t ) = lim h 0 + h α n = 0 ( α ) n n ! e n λ h f ( t n h ) ,
for α R .
2.
Forward regularized derivative
D λ , f α f ( t ) = 0 f ( t τ ) 0 N 1 ( 1 ) m f ( m ) ( t ) m ! τ m e λ τ τ α 1 Γ ( α ) d τ
where α N , again. If N 0 , the summation is null.
It can be shown that [63]
D λ , f β f ( t ) = e λ t D 0 , f β e λ t f ( t )
and
D λ , b β f ( t ) = e λ t D 0 , b β e λ t f ( t ) .
The approach to introducing tempered scale derivatives is similar.
Definition 6.
Let μ R 0 + . The tempered SD have transfer functions defined by
G ( v ) = ( v + μ ) α ± ( R e ( v ) + μ ) > 0 .
Theorem 5.
Let q > 1 . The following expressions
D μ , s + α x ( τ ) = lim q 1 + ln α ( q ) n = 0 ( α ) n n ! q n μ x ( τ q n ) ,
and
D μ , s α x ( τ ) = lim q 1 + ( 1 ) α ln α ( q ) n = 0 ( α ) n n ! q n μ x ( τ q n )
represent the tempered scale-derivatives corresponding to the stretching (27) and shrinking (28) GL-type derivatives, respectively.
The proof is readily found from
( v + μ ) α = lim q 1 + ( 1 q v μ ) ln q α R e ( v ) > μ , ( q v + μ 1 ) ln q α R e ( v ) < μ ,
To obtain the Hadamard tempered derivatives, we have to compute the corresponding impulse response functions, a simple task attending to the properties of the MT. We obtain
M 1 ( v + μ ) α ( τ ) = τ μ ln α 1 ( τ ) Γ ( α ) ε ( τ 1 ) ( R e ( v ) + μ ) > 0 τ μ ln α 1 ( 1 / τ ) Γ ( α ) ε ( 1 τ ) ( R e ( v ) + μ ) < 0
With this relation, we can obtain the tempered Hadamard derivatives. In particular, the tempered regularized Hadamard derivatives are given by
D s + α x ( τ ) = τ μ Γ ( α ) 0 τ u μ x ( u ) n = 0 N 1 ( 1 ) n D s + N x ( τ ) n ! ln n ( τ / u ) ln α 1 ( τ / u ) d u u
and
D s α x ( τ ) = τ μ Γ ( α ) τ u μ x ( u ) n = 0 N 1 D s + N x ( τ ) n ! ln n ( u / τ ) ln α 1 ( u / τ ) d u u .
From (67) and (68), we obtain
D μ , s + α x ( τ ) = τ μ D 0 , s + α x ( τ ) τ μ f ( t )
and
D μ , s α x ( τ ) = τ μ D 0 , s α x ( τ ) τ μ f ( t ) ,
similar to (69) and (70).

5. Scale-Invariant Systems

Definition 7.
In agreement with the concepts introduced above, we define the dilation (scale)-invariant fractional autoregressive-moving average (DI-FARMA) system through
k = 0 N 0 a k D s ± α k y ( τ ) = k = 0 M 0 b k D s ± β k x ( τ ) τ R +
where D s ± α k , ( β k ) , k = 0 , 1 , 2 , mean the fractional α k ( β k ) -order scale-derivatives and N 0 , M 0 are the system orders. The parameters a k , b k , k = 0 , 1 , are considered real numbers. Without losing generality, we set a N 0 = 1 .
We could use tempered SD but we will not do so here. The results in [63] can be easily adapted. The name “autoregressive-moving average” was borrowed from a current nomenclature used for shift-invariant systems [64]. As the power t v is the eigenfunction of (71), we obtain easily the transfer function
G ( v ) = k = 0 M 0 b k v β k k = 0 N 0 a k v α k .
which is a bit difficult to manipulate [65]. In the following, we shall be considering the so-called “commensurate” systems described by differential equations with the format
k = 0 N 0 a k D s ± k α y ( τ ) = k = 0 M 0 b k D s ± k α x ( τ ) τ R +
The corresponding transfer function is
G ( v ) = k = 0 M 0 b k v k α k = 0 N 0 a k v k α ,
where we assume that M 0 < N 0 for simplicity, and all the roots, p k , k = 1 , 2 , , of k = 0 N 0 a k w k are simple, which allows us to write
G ( v ) = k = 1 N A k v α p k ,
where the A k are the residues obtained by substituting w for v α in (72). The impulse response results from the inversion of a combination of partial fractions such as
F ( v ) = 1 v α p .
We will treat the case ( R e ( v ) > 0 ). We write
1 v α p = n = 0 p n v n α 1 , R e ( v ) > | p | .
The corresponding inverse MT is
f ( τ ) = n = 0 p n ln n α ( τ ) Γ ( n α + 1 ) , τ 1 ,
which is analog to the α exponential function [28] and can be expressed as
f ( τ ) = ln β 1 ( τ ) E α , β ( p ln α ( τ ) ) , τ 1 ,
where we introduced the logarithmic Mittag-Leffler function
E α , β ( ln α ( τ ) ) = n = 0 ln n α ( τ ) Γ ( n α + 1 ) , τ 1 .
The solution corresponding to R e ( v ) < 0 is obtained using a similar procedure, accounting for (48).
The α exponential function is more useful, since it allows us to obtain the inversion in the multiple-order pole case:
F ( v ) = 1 ( v α p ) n , n Z + .
This case can be treated easily from (76) through a n t h -order derivative assuming p as the variable.
Example 3.
Consider the simple AR(1) system
D y ( τ ) + y ( τ ) = x ( τ )
Its transfer function is
G ( v ) = 1 v + 1 = n = 0 ( 1 ) n v n 1 , R e ( v ) > 1 .
Using (48), the impulse response is
g ( τ ) = n = 0 ( 1 ) n ln n ( τ ) n ! ε ( τ 1 ) = e ln τ ε ( τ 1 ) = 1 τ ε ( τ 1 )
where we used (48) if R e ( v ) > 0 . In the R e ( v ) < 1 case, we obtain
g ( τ ) = n = 0 ln n ( τ ) n ! ε ( 1 τ ) = e ln ( τ ) ε ( 1 τ ) = τ ε ( 1 τ )
Remark 8.
It is important to note that, according to the theoretical framework that we developed in the previous section, the initial conditions must be taken at τ = 1 ± , not at τ = 0 ± . A similar scheme for introducing those initial conditions as the one in [48] can be found.
We could define two “unilateral” MTs with corresponding initial conditions, but that are not necessarily those in any system, as it happens with the one-sided Laplace transform [66]. For instance,
F u ( v ) = 1 f ( τ ) τ v 1 d τ .
Example 4.
In Section 4.4, we introduced the scale conversion operator that reads
S ( v ) = a v = n = 0 v n n ! ln n ( a )
Using signal processing language, we can call it an infinite impulse response system. We can obtain an interesting approximation. Assume that a = 1 + η with η being a small positive number. In such a case, ln ( a ) = ln ( 1 + η ) η . Then, we can write
S ( v ) = a v 2 a v 2 = e ln ( a ) 2 v e ln ( a ) 2 v 1 + η 2 v 1 η 2 v
Therefore, S ( v ) is the transfer function of a linear system described by the differential equation
η 2 v D y ( τ ) + y ( τ ) = η 2 D x ( τ ) + x ( τ )
If η is not small, we can decompose a into a product; for example, a = k = 1 n a 1 n and use n cascaded systems. We do not go further here.

6. Conclusions

Using the concept of scale-invariant linear systems, a general fractional scale derivative was introduced and studied. A Liouville-like framework for such systems was described. The relationship between scale derivatives and Hadamard’s was re-established and reformulated. A new derivative similar to Grünwald–Letnikov’s was deduced. Tempered versions were also introduced. Scale-invariant systems with an autoregressive-moving average form differential equation were introduced. For solving it, a new logarithmic Mittag-Leffler series was proposed.

Author Contributions

Conceptualization, M.D.O. and G.W.B.; methodology, M.D.O.; formal analysis, M.D.O.; investigation, M.D.O. and G.W.B.; writing—original draft preparation, M.D.O.; writing—review and editing, M.D.O. and G.W.B. All authors have read and agreed to the published version of the manuscript.

Funding

The first author was partially funded by National Funds through the Foundation for Science and Technology of Portugal, under the projects UIDB/00066/2020.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Cohen, L. The scale representation. IEEE Trans. Signal Process. 1993, 41, 3275–3292. [Google Scholar] [CrossRef]
  2. Nottale, L. The theory of scale relativity. Int. J. Mod. Phys. A 1992, 7, 4899–4936. [Google Scholar] [CrossRef]
  3. Nottale, L. Non-differentiable space-time and scale relativity. In Proceedings of the International Colloquium Geometrie au XXe Siecle, Paris, France, 24–29 September 2001. [Google Scholar]
  4. Nottale, L. The Theory of Scale Relativity: Non-Differentiable Geometry and Fractal Space-Time. In Computing Anticipatory Systems. CASYS’03—Sixth International Conference. American Institute of Physics Conference Proceedings; Dubois, D.M., Ed.; American Institute of Physics: College Park, MA, USA, 2004; Volume 718, pp. 68–75. [Google Scholar]
  5. Cresson, J. Scale relativity theory for one-dimensional non-differentiable manifolds. Chaos Solitons Fractals 2002, 14, 553–562. [Google Scholar] [CrossRef]
  6. Cresson, J. Scale calculus and the Schrödinger equation. J. Math. Phys. 2003, 44, 4907–4938. [Google Scholar] [CrossRef]
  7. Proekt, A.; Banavar, J.R.; Maritan, A.; Pfaff, D.W. Scale invariance in the dynamics of spontaneous behavior. Proc. Natl. Acad. Sci. USA 2012, 109, 10564–10569. [Google Scholar] [CrossRef]
  8. Khaluf, Y.; Ferrante, E.; Simoens, P.; Huepe, C. Scale invariance in natural and artificial collective systems: A review. J. R. Soc. Interface 2017, 14, 20170662. [Google Scholar] [CrossRef]
  9. Lamperti, J. Semi-stable stochastic processes. Trans. Am. Math. Soc. 1962, 104, 62–78. [Google Scholar] [CrossRef]
  10. Borgnat, P.; Amblard, P.O.; Flandrin, P. Scale invariances and Lamperti transformations for stochastic processes. J. Phys. A Math. Gen. 2005, 38, 2081. [Google Scholar] [CrossRef]
  11. Belbahri, K. Scale invariant operators and combinatorial expansions. Adv. Appl. Math. 2010, 45, 548–563. [Google Scholar] [CrossRef]
  12. Grossmann, A.; Morlet, J. Decomposition of Hardy Functions into Square Integrable Wavelets of Constant Shape. SIAM J. Math. Anal. 1984, 15, 723–736. [Google Scholar] [CrossRef]
  13. Meyer, Y. Orthonormal wavelets. In Wavelets: Time-Frequency Methods and Phase Space Proceedings of the International Conference, Marseille, France, 14–18 December 1987; Springer: Berlin/Heidelberg, Germany, 1989; pp. 21–37. [Google Scholar]
  14. Mallat, S.G. Multiresolution Representations and Wavelets; University of Pennsylvania: Philadelphia, PA, USA, 1988. [Google Scholar]
  15. Mallat, S. A theory for multiresolution signal decomposition: The wavelet representation. IEEE Trans. Pattern Anal. Mach. Intell. 1989, 11, 674–693. [Google Scholar] [CrossRef]
  16. Daubechies, I. The wavelet transform, time-frequency localization and signal analysis. IEEE Trans. Inf. Theory 1990, 36, 961–1005. [Google Scholar] [CrossRef]
  17. Heil, C.E.; Walnut, D.F. Continuous and discrete wavelet transforms. SIAM Rev. 1989, 31, 628–666. [Google Scholar] [CrossRef]
  18. Edwards, T. Discrete Wavelet Transforms: Theory and Implementation; Stanford University: Stanford, CA, USA, 1991; pp. 28–35. [Google Scholar]
  19. Van Fleet, P.J. Discrete Wavelet Transformations: An Elementary Approach with Applications; John Wiley & Sons: Hoboken, NJ, USA, 2019. [Google Scholar]
  20. Poularikas, A.D. The Transforms and Applications Handbook; CRC Press LLC: Boca Raton, FL, USA, 2000. [Google Scholar]
  21. Hadamard, J. Essai sur L’étude des Fonctions, Données par leur Développement de Taylor; Gallica: Tokyo, Japan, 1892; pp. 101–186. [Google Scholar]
  22. Samko, S.G.; Kilbas, A.A.; Marichev, O.I. Fractional Integrals and Derivatives; Gordon and Breach: Yverdon, Switzerland, 1993. [Google Scholar]
  23. Garra, R.; Orsingher, E.; Polito, F. A Note on Hadamard Fractional Differential Equations with Varying Coefficients and Their Applications in Probability. Mathematics 2018, 6, 4. [Google Scholar] [CrossRef]
  24. Tarasov, V.E. Fractional dynamics with non-local scaling. Commun. Nonlinear Sci. Numer. Simul. 2021, 102, 105947. [Google Scholar] [CrossRef]
  25. Liouville, J. Memóire sur quelques questions de Géométrie et de Méchanique, et sur un nouveau genre de calcul pour résoudre ces questions. J. L’École Polytech. Paris 1832, 13, 1–69. [Google Scholar]
  26. Liouville, J. Memóire sur le calcul des différentielles à indices quelconques. J. L’École Polytech. Paris 1832, 13, 71–162. [Google Scholar]
  27. Dugowson, S. Les Différentielles Métaphysiques (Histoire et Philosophie de la Généralisation de L’ordre de Dérivation). Ph.D. Thesis, Université Paris Nord, Paris, France, 1994. [Google Scholar]
  28. Kilbas, A.A.; Srivastava, H.M.; Trujillo, J.J. Theory and Applications of Fractional Differential Equations; Elsevier: Amsterdam, The Netherlands, 2006. [Google Scholar]
  29. Mainardi, F. Fractional Calculus and Waves in Linear Viscoelasticity: An Introduction to Mathematical Models; Imperial College Press: London, UK, 2010. [Google Scholar]
  30. Podlubny, I. Fractional Differential Equations; Academic Press: San Diego, CA, USA, 1999; Volume 198, p. 340. [Google Scholar]
  31. Ortigueira, M.D. Fractional Calculus for Scientists and Engineers; Lecture Notes in Electrical Engineering; Springer: Dordrecht, The Netherlands, 2011. [Google Scholar]
  32. Riemann, B. Versuch einer allgemeinen Auffassung der Integration und Differentiation (Jan. 14, 1847). In The Collected Works of Bernard Riemann Edited by Heinrich Weber with the Assistance of Richard Dedekind; Dover Publications: New York, NY, USA, 1953; pp. 353–366. [Google Scholar]
  33. Valério, D.; Ortigueira, M.D.; Lopes, A.M. How Many Fractional Derivatives Are There? Mathematics 2022, 10, 737. [Google Scholar] [CrossRef]
  34. De Oliveira, E.C.; Tenreiro Machado, J.A. A review of definitions for fractional derivatives and integrals. Math. Probl. Eng. 2014, 2014, 238459. [Google Scholar] [CrossRef]
  35. 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]
  36. Kilbas, A.A. Hadamard-type fractional calculus. J. Korean Math. Soc. 2001, 38, 1191–1204. [Google Scholar]
  37. Butzer, P.L.; Kilbas, A.A.; Trujillo, J.J. Compositions of Hadamard-type fractional integration operators and the semigroup property. J. Math. Anal. Appl. 2002, 269, 387–400. [Google Scholar] [CrossRef]
  38. Klimek, M. Sequential fractional differential equations with Hadamard derivative. Commun. Nonlinear Sci. Numer. Simul. 2011, 16, 4689–4697. [Google Scholar] [CrossRef]
  39. Garra, R.; Polito, F. On some operators involving Hadamard derivatives. Integral Transform. Spec. Funct. 2013, 24, 773–782. [Google Scholar] [CrossRef]
  40. Kamocki, R. Necessary and sufficient conditions for the existence of the Hadamard-type fractional derivative. Integral Transform. Spec. Funct. 2015, 26, 442–450. [Google Scholar] [CrossRef]
  41. Chen, H.; Katugampola, U.N. Hermite–Hadamard and Hermite–Hadamard–Fejér type inequalities for generalized fractional integrals. J. Math. Anal. Appl. 2017, 446, 1274–1291. [Google Scholar] [CrossRef]
  42. Zheng, X. Logarithmic transformation between (variable-order) Caputo and Caputo–Hadamard fractional problems and applications. Appl. Math. Lett. 2021, 121, 107366. [Google Scholar] [CrossRef]
  43. Liu, W.; Liu, L. Properties of Hadamard Fractional Integral and Its Application. Fractal Fract. 2022, 6, 670. [Google Scholar] [CrossRef]
  44. Ortigueira, M.D.; Valério, D. Fractional Signals and Systems; De Gruyter: Berlin, Germany; Boston, MA, USA, 2020. [Google Scholar]
  45. Kailath, T. Linear Systems; Information and System Sciences Series; Prentice-Hall: Hoboken, NJ, USA, 1980. [Google Scholar]
  46. Bengochea, G.; Ortigueira, M.D. Fractional derivative of power type functions. Comput. Appl. Math. 2022, 41, 1–18. [Google Scholar]
  47. Ortigueira, M.D.; Machado, J.A.T. Fractional Derivatives: The Perspective of System Theory. Mathematics 2019, 7, 150. [Google Scholar] [CrossRef]
  48. Ortigueira, M.D. The complex order fractional derivatives and systems are non hermitian. In Proceedings of the International Conference on Fractional Differentiation and Its Applications (ICFDA’21), Online, 6–8 September 2021; Springer: Berlin/Heidelberg, Germany, 2022; pp. 38–44. [Google Scholar]
  49. Herrmann, R. Fractional Calculus, 3rd ed.; World Scientific: Singapore, 2018. [Google Scholar]
  50. Rudolf, H. Applications of Fractional Calculus in Physics; World Scientific: Singapore, 2000. [Google Scholar]
  51. Ortigueira, M.; Bengochea, G. A new look at the fractionalization of the logistic equation. Phys. A Stat. Mech. Its Appl. 2017, 467, 554–561. [Google Scholar] [CrossRef]
  52. Butzer, P.L.; Jansche, S. A direct approach to the Mellin transform. J. Fourier Anal. Appl. 1997, 3, 325–376. [Google Scholar] [CrossRef]
  53. Luchko, Y.; Kiryakova, V. The Mellin integral transform in fractional calculus. Fract. Calc. Appl. Anal. 2013, 16, 405–430. [Google Scholar] [CrossRef]
  54. Kac, V.G.; Cheung, P. Quantum Calculus; Springer: New York, NY, USA, 2002; Volume 113. [Google Scholar]
  55. Ernst, T. A Comprehensive Treatment of q-Calculus; Birkhäuser: Basel, Switzerland, 2012. [Google Scholar]
  56. Ortigueira, M.D. The fractional quantum derivative and its integral representations. Commun. Nonlinear Sci. Numer. Simul. 2010, 15, 956–962. [Google Scholar] [CrossRef]
  57. Ortigueira, M.D. On the Fractional Linear Scale Invariant Systems. IEEE Trans. Signal Process. 2010, 58, 6406–6410. [Google Scholar] [CrossRef]
  58. Jarad, F.; Abdeljawad, T.; Baleanu, D. Caputo-type modification of the Hadamard fractional derivatives. Adv. Differ. Equ. 2012, 2012, 1–8. [Google Scholar] [CrossRef]
  59. Ma, L.; Li, C. On Hadamard fractional calculus. Fractals 2017, 25, 1750033. [Google Scholar] [CrossRef]
  60. Almeida, R. Caputo–Hadamard fractional derivatives of variable order. Numer. Funct. Anal. Optim. 2017, 38, 1–19. [Google Scholar] [CrossRef]
  61. Ma, L.; Li, C. On finite part integrals and Hadamard-type fractional derivatives. J. Comput. Nonlinear Dyn. 2018, 13, 090905. [Google Scholar] [CrossRef]
  62. 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]
  63. Ortigueira, M.D.; Bengochea, G.; Machado, J.A.T. Substantial, tempered, and shifted fractional derivatives: Three faces of a tetrahedron. Math. Methods Appl. Sci. 2021, 44, 9191–9209. [Google Scholar] [CrossRef]
  64. Ortigueira, M.D.; Magin, R.L. On the Equivalence between Integer-and Fractional Order-Models of Continuous-Time and Discrete-Time ARMA Systems. Fractal Fract. 2022, 6, 242. [Google Scholar] [CrossRef]
  65. Bengochea, G.; Ortigueira, M.; Verde-Star, L. Operational calculus for the solution of fractional differential equations with noncommensurate orders. Math. Methods Appl. Sci. 2021, 44, 8088–8096. [Google Scholar] [CrossRef]
  66. Ortigueira, M.D.; Machado, J.T. Revisiting the 1D and 2D Laplace transforms. Mathematics 2020, 8, 1330. [Google Scholar] [CrossRef]
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ortigueira, M.D.; Bohannan, G.W. Fractional Scale Calculus: Hadamard vs. Liouville. Fractal Fract. 2023, 7, 296. https://doi.org/10.3390/fractalfract7040296

AMA Style

Ortigueira MD, Bohannan GW. Fractional Scale Calculus: Hadamard vs. Liouville. Fractal and Fractional. 2023; 7(4):296. https://doi.org/10.3390/fractalfract7040296

Chicago/Turabian Style

Ortigueira, Manuel D., and Gary W. Bohannan. 2023. "Fractional Scale Calculus: Hadamard vs. Liouville" Fractal and Fractional 7, no. 4: 296. https://doi.org/10.3390/fractalfract7040296

APA Style

Ortigueira, M. D., & Bohannan, G. W. (2023). Fractional Scale Calculus: Hadamard vs. Liouville. Fractal and Fractional, 7(4), 296. https://doi.org/10.3390/fractalfract7040296

Article Metrics

Back to TopTop