Next Article in Journal
Mathematical Calculation of Synchronous Electric Motors Dynamic Stability
Next Article in Special Issue
Hyers–Ulam Stability of Caputo Fractional Stochastic Delay Differential Systems with Poisson Jumps
Previous Article in Journal
A No-Chatter Single-Input Finite-Time PID Sliding Mode Control Technique for Stabilization of a Class of 4D Chaotic Fractional-Order Laser Systems
Previous Article in Special Issue
Pontryagin Maximum Principle for Incommensurate Fractional-Orders Optimal Control Problems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

On the Fractional Derivative Duality in Some Transforms

by
Manuel Duarte Ortigueira
1,* and
Gabriel Bengochea
2
1
NOVA School of Science and Technology, UNINOVA-CTS and LASI, NOVA University of Lisbon, Quinta da Torre, 2829-516 Caparica, Portugal
2
Academia de Matemáticas, Universidad Autónoma de la Ciudad de México, Prolongación San Isidro 151, Col. San Lorenzo Tezonco, Del. Iztapalapa, Ciudad de México C.P. 09790, Mexico
*
Author to whom correspondence should be addressed.
Mathematics 2023, 11(21), 4464; https://doi.org/10.3390/math11214464
Submission received: 3 October 2023 / Revised: 23 October 2023 / Accepted: 24 October 2023 / Published: 27 October 2023
(This article belongs to the Special Issue Recent Research on Fractional Calculus: Theory and Applications)

Abstract

:
Duality is one of the most interesting properties of the Laplace and Fourier transforms associated with the integer-order derivative. Here, we will generalize it for fractional derivatives and extend the results to the Mellin, Z and discrete-time Fourier transforms. The scale and nabla derivatives are used. Some consequences are described.

1. Introduction

Transform is a term that originated a long time ago, probably in the work of Euler. In fact, the Laplace transform (LT) was introduced by him in 1737 [1]. However, it is used in many applications, not only in mathematics, but in all sciences and engineering. We therefore need to restrict our definition so that it is in line with our objectives. We will use the designation transform for operators defined in R or Z that generate functions defined in C and that take on integral or series format. The classic integral transforms, namely the Fourier, Laplace, Mellin, Hankel and Stieltjes transforms [2,3,4,5], and some variants are the most useful in practical applications. In engineering, other transforms have been introduced, such as the fractional Fourier transform (FT) [6,7], the Wavelet transform [8,9] and the Radon transform [4]. In recent years, many transforms have been introduced, without their usefulness being clear [10]. Most are modified one-sided Laplace transforms.
Although many transforms are defined over continuous domains, there are some discrete transforms such as the Z transform [11], also called characteristic function, the discrete-time Fourier transform [12,13], the discrete Mellin transform [14] and the Fermat number transform [15].
Here, we will consider those transforms that have direct relations with linear systems:
  • the continuous-time (two-sided) Laplace and the Fourier transforms that are associated with the linear systems described by shift-invariant convolution [11,12,16,17];
  • the continuous-scale Mellin transform that is tied with the linear systems described by the Mellin convolution [17,18];
  • the discrete-time Z and Fourier transforms connected with the discrete shift-invariant convolution [11,12,16].
We have not considered the discrete Mellin transform here, as it poses some difficulties that have not been sufficiently studied [19].
These transforms have a duality property that can be expressed as a reversibility characteristic that can be stated as follows: if a given transform has a particular property, the inverse transform has a similar property. The extreme expression of this property is found in the Fourier transform (FT). To understand this idea, we consider the FT of an absolutely integrable function, f ( t ) ,
F ( i ω ) = F f ( t ) = f ( t ) e i ω t d t ,
and its inverse
f ( t ) = F 1 F ( i ω ) = 1 2 π F ( i ω ) e i ω t d ω .
Duality says that
F F ( i t ) = 2 π f ( ω ) ,
which is very useful in the computation of some transforms.
In this paper, we are interested in the derivative duality. Let D represent the usual derivative and n N . As D ω n e i ω t = ( i t ) n e i ω t , we have
F ( n ) ( i ω ) = ( i t ) n f ( t ) e i ω t d t ,
and, similarly, D t n e i ω t = ( i ω ) n e i ω t ,
f ( n ) ( t ) = 1 2 π ( i ω ) n F ( i ω ) e i ω t d ω .
These relations clearly express the duality property of FT and the associated usual derivative: deriving t (ω) in one domain corresponds to multiplying ω (t) in the other domain. This property is useful in transforming some variable coefficient ordinary differential equations.
In the current literature, only Fourier transform cases are dealt with. The fractional case is not considered anywhere. Here, we want to generalize the derivative duality to fractional orders and to other transforms. We will consider Laplace, Mellin (MT), Z (ZT) and discrete-time Fourier transforms (DTFT). To achieve this, we need to introduce the necessary fractional derivatives: Liouville, Hadamard and discrete, nabla and bilinear. The best-known fractional derivatives are not useful for our purposes, as we shall see.
The paper is outlined as follows. In Section 2, we introduce the required properties for the fractional derivatives to use. The integer order case of duality is dealt with in Section 3, and the fractional cases are studied in Section 4. In Section 5, we describe some consequences of the results we obtained in the previous sections. Finally, we present some conclusions.

2. Suitable Derivatives

2.1. Derivative Requirements

The most important transforms have exponentials or powers as kernels due to their relations with shift-invariant or scale-invariant systems [17,18]. Therefore, if we want to compute derivatives of the corresponding transforms and continue having a transform of the same type, we must have derivatives verifying (Liouville requirement)
D z α e a z = a α e a z ,
or (Hadamard requirement)
D z β z b = b β z b ,
for suitable orders α , β R and complex variable z C . For simplicity matters and later utility, we will assume that a , b R also.
Remark 1. 
It is important to note that the usual Riemann–Liouville and Caputo derivatives are not useful in this case since they do not check relationships (3) and (4).

2.2. Liouville-Type Derivatives

The Liouville-type derivatives can be expressed in a general unified way [20,21] assuming both summation and integral formulations.
Definition 1. 
We will consider the Grünwald–Letnikov (GL) derivatives that we define by:
D ± α f ( z ) : = lim h 0 + ( ± h ) α n = 0 + α n n ! f ( z n h ) ,
where we denoted by ( a ) n , n = 1 , 2 , the Pochhammer symbol for the rising factorial
( a ) 0 = 1 , ( a ) n = k = 0 n 1 ( a + k ) .
To avoid confusion of symbols for different derivatives, we will write frequently D ± α f ( t ) = f ± ( α ) ( t ) . If necessary, we will put the independent variable in the subscript.
Theorem 1. 
Let f ( z ) = e a z . Then, [21,22]
D ± α e a z = a α e a z ,
if ± R e ( a ) > 0 , while both diverge with ± R e ( a ) < 0 , unless α Z .
For applications to the Laplace transform, we need to compute derivatives of e s t , t R , s C . We obtain:
  • Derivative in t
    D t ± α e s t = s α e s t , ± R e ( s ) > 0 ;
    it is useful for the derivative computation of the inverse LT.
  • Derivative in s
    D s ± α e s t = ( t ) α e s t , t > 0 ;
    this derivative allows us to express the derivative of the direct LT.

2.3. Hadamard-Type Derivatives

Definition 2. 
The Hadamard-type derivatives are scale-invariant and verify the above requirement (4). Similarly to the Liouville’s type, we define the stretching (+) and shrinking (−) GL-type derivatives by [18]
D v ± α x ( v ) = lim q 1 + ln α ( q ± 1 ) n = 0 ( α ) n n ! x ( v q n ) ,
where q > 1 .
Theorem 2. 
D v ± α v b = b α v b ,
provided that ± R e ( v ) > 0 . If ± R e ( v ) < 0 , they diverge.
Similarly to the LT case, we need to compute derivatives of τ v , τ R + , v C . We obtain:
  • Derivative in τ
    D τ ± α τ v = v α τ v , ± R e ( v ) > 0 ;
    it is useful for the derivative computation of the inverse MT.
  • Derivative in v
    D v ± α τ v = ( ln τ ) α τ v , τ 1 > 1 ;
    this derivative allows us to express the derivative of the direct MT.

2.4. Discrete-Time Derivatives

2.4.1. Fractional Nabla and Delta Derivatives

In the following, we consider that our domain is the time scale or time sequence
T h = ( h Z ) = , n h , 2 h , h , 0 , h , 2 h , , n h , ,
with h R + , which is called the graininess or sampling interval [23,24].
Definition 3. 
Let f ( t ) be a function defined on T . Set t = n h . We define the nabla derivative by:
f ( t ) = f ( t ) f ( t h ) h
and the delta derivative by
Δ f ( t ) = f ( t + h ) f ( t ) h .
The corresponding fractional derivatives read [24]:
nabla
α f ( t ) = n = 0 ( α ) n n ! f ( t n h ) h α
and delta
Δ α f ( t ) = n = 0 ( α ) n n ! f ( t + n h ) h α .
Theorem 3. 
The eigenvalue of these derivatives is s α , and the corresponding eigenfunctions are the nabla and delta exponentials given by [24]
e ( k h , s ) = 1 ( 1 s h ) k
and
e Δ ( k h , s ) = ( 1 + s h ) k .
These results were generalized for irregular time sequences in [24]. In the following, we will continue with the nabla derivative.
Definition 4. 
With the nabla exponential, we can define the nabla Laplace transform [24] through
F ( s ) = h k = + f ( k h ) e k h , s ,
with its inverse transform being given by
f ( k h ) = 1 2 π j γ F ( s ) e ( k + 1 ) h , s d s ,
where the integration path, γ , is any simple closed contour in a region of analyticity of the integrand that includes the point s = 1 h . The simplest path is a circle with center at s = 1 h .
Corollary 1. 
Let z 1 = 1 s h . Then
α z n = 1 z 1 h α z n .
The proof is immediate.
With this change, we entered in the framework of the Z transform.

2.4.2. Forward and Backward Derivatives Based on the Bilinear Transformation

The Tustin transformation is usually expressed by [11,13]
s = 2 h 1 z 1 1 + z 1 ,
where s is the derivative operator associated with the (continuous-time) LT and z 1 the delay operator tied to the Z transform.
Definition 5. 
Let x ( n h ) be a discrete-time function. We define the order 1 forward or nabla bilinear derivative b x ( n h ) of x ( n h ) as the solution of the difference equation
b x ( n h ) + b x ( n h h ) = 2 h x ( n h ) x ( n h h ) .
Similarly, we define the order 1 backward or delta bilinear derivative Δ b x ( n h ) of x ( n h ) as the solution of
Δ b x ( n h + h ) + Δ b x ( n h ) = 2 h x ( n h + h ) x ( n h ) .
Definition 6. 
We define the nabla bilinear derivative ( b ) as an elementary DT system such that
b z n = 2 h 1 z 1 1 + z 1 z n .
The transfer function of such derivative, H b ( z ) , is defined by
H ( z ) = 2 h 1 z 1 1 + z 1 , z > 1 .
For the backward bilinear derivative, a transfer function is defined similarly.
Let α R . The α -order nabla bilinear fractional derivative is a discrete-time linear system with transfer function
H ( z ) = 2 h 1 z 1 1 + z 1 α , | z | > 1 ,
such that
b α z n = 2 h 1 z 1 1 + z 1 α z n , | z | > 1 .
With this formulation, we entered again in the context of the Z transform.

3. Main Transforms and Integer-Order Derivatives

3.1. Continuous-Time Laplace and Fourier Transforms

Definition 7. 
The direct bilateral LT (BLT) is given by [25]
L f ( t ) = F ( s ) = f ( t ) e s t d t , s C R c ,
while the inverse LT (synthesis equation) reads
f ( t ) = L 1 F ( s ) = 1 2 π i a i a + i F ( s ) e s t d s , t R ,
where a R must be located inside the region of convergence, R c , (ROC) of F ( s ) . The right-hand side represents the Bromwich integral. In the following, we will denote by γ the integration path.
We can obtain existence conditions for the BLT from those of the FT [12,26,27]. Let f ( t ) be a function
  • piecewise continuous,
  • with bounded variation,
  • locally integrable (in the sense that the function is absolutely integrable in any real interval [ a , b ] , so that a b | f ( t ) | d t < ),
  • of exponential order,
then there exists the BLT of f ( t ) .
Remark 2. 
Loosely speaking, a function of exponential order is the one that does not “grow faster” than given exponentials as t ± . This means two things. First, that there are real constants A, a > 0 such that | f ( t ) | < A · e a t , when t is large and negative (say, for t < t 1 R ). Second, that there are real constants B, b > 0 such that | f ( t ) | < B · e b t , when t is large (say, for t > t 2 R ). It also has to be true that b < a . We are interested in dealing with functions for which b < 0 and a > 0 so that the function has Fourier transform.
Under the conditions indicated, the integral in (29) converges absolutely and uniformly in a vertical band in the complex plane defined by b < R e ( s ) < a , where F ( s ) is analytic. This band is called the region of convergence (ROC), and the values of the constants a and b are the abscissae of convergence. It can be shown that:
  • If f ( t ) is absolutely integrable and of finite duration, then the ROC is the entire s-plane since the Laplace transform is finite and F ( s ) exists for any s.
  • If f ( t ) is right-handed (i.e., it exists with t t 0 R ) and R e ( s ) = a R c , then any s to the right of a is also in R c .
  • If f ( t ) is left-handed (i.e., exists with t t 0 R ) and R e ( s ) = a R c , then any s to the left of a is also in R c .
  • A function f ( t ) is absolutely integrable (satisfying the Dirichlet conditions and having the Fourier transform) if and only if the ROC of the corresponding Laplace transform F ( s ) includes the imaginary axis since R e ( s ) = 0 and s = i ω .
  • A given complex variable function only can define univocally an LT if it has attached a suitable ROC.
  • If F ( s ) = L f ( t ) , then L f ( t ) = F ( s ) .
  • If the region of convergence of F ( s ) includes the frontiers b R e ( s ) a , then F ( s ) is completely defined in that region by the values at the lines F ( a + i τ ) and F ( b + i τ ) , τ R .
  • F ( s ) is bounded in the strip a + ϵ R e ( s ) b ϵ , with ϵ > 0 .
As the integer-order GL derivative of an exponential exists for any value, provided that s 0 , we obtain easily
Theorem 4. 
F ( n ) ( s ) = ( t ) n f ( t ) e s t d t , s C
and
f ( n ) ( t ) = 1 2 π i a i a + i s n F ( s ) e s t d s , t R .
From these relations, we obtain the corresponding properties that we introduced above for the FT.

3.2. The Mellin Transform

Definition 8. 
Let us define the Mellin transform by
G ( v ) = 0 g ( u ) u v 1 d u ,
which is a modified version of the usual Mellin transform. This has a parameter sign change v v relative to the current [14]. The idea is to keep the parallelism with the LT so that the properties related to the transform domain are essentially the same. In fact, it results from the change e t u . The inverse Mellin transform related to (10) is
g ( τ ) = M 1 X ( v ) = 1 2 π i γ G ( v ) τ v d v , τ R + ,
where γ is a vertical straight line in the ROC of the transform.
As the MT results from the LT through an exponential variable change, the convergence properties are easily deduced. In particular, the integral in (33) converges absolutely and uniformly in a vertical strip in the complex plane defined by b < R e ( v ) < a , where G ( v ) is analytic.
Theorem 5. 
The duality in the MT is expressed by
D ± n g ( τ ) = 1 2 π i γ ( ± v ) n G ( v ) τ v d v , τ R + ,
and
G ( n ) ( v ) = 0 ( ln ( u ) ) n g ( u ) u v 1 d u , v C .
The proof is immediate using the scale-derivative in the inverse and the classical derivative in the direct MT. Although this result may seem strange, it expresses the effect of the exponential transformation mentioned above.

3.3. On the Z and Discrete-Time Fourier Transforms

Definition 9. 
Let x ( n ) denote any function defined on T , leaving implicit the graininess, unless it is convenient to display it. The Z transform (ZT) is defined by
X ( z ) = Z x ( n ) = n = x ( n ) z n , z C .
The inverse ZT can be obtained by the integral defined by
x ( n ) = 1 2 π i γ X ( z ) z n 1 d z ,
where γ is a circle centered at the origin, located in the ROC of the transform and taken in a counterclock wise direction.
In some scientific fields, such as geophysics, z is used instead of z 1 , and sometimes the ZT is called the “generating function” or “characteristic function”. The existence conditions of the ZT are similar to those of the bilateral LT [11,13,16]. They can be stated as follows.
Definition 10. 
A discrete-time signal x ( n ) is called an exponential order signal if there exist integers n 1 and n 2 and positive real numbers a , b , A , and B such that A a n 1 < | x ( n ) | < B b n 2 for n 1 < n < n 2 .
For these signals, the ZT exists, and the ROC is an annulus centered at the origin, generally delimited by two circles of radius r and r + such that r < | z | < r + . However, there are some cases where the annulus can become infinite:
  • If the signal is right (i.e., x ( n ) = 0 , n < n 0 Z ), then the ROC is the exterior of a circle centered at the origin ( r + = ): | z | > r .
  • If the signal is left (i.e., x ( n ) = 0 , n > n 0 Z ), then the ROC is the interior of a circle centered at the origin ( r = 0 ): | z | < r + .
  • If the signal is a pulse (i.e., non null only on a finite set), then the ROC is the whole complex plane, possibly with the exception of the origin. In the ROC, the ZT defines an analytical function.
If the ROC contains the unit circle, then by making z = e i ω , ω < π , i = 1 , we obtain the discrete-time Fourier transform, which we will shortly call Fourier transform. This means that not all signals with ZT have FT. The signals with ZT and FT are those for which the ROC is non-degenerate and contains the unit circle ( r < 1 , r + > 1 ).
In this situation, the integral in (37) converges uniformly. The calculation uses Cauchy’s theorem for functions of complex variable [16].
To treat the duality, we must note the importance of the unit circle that suggests the use of the scale-derivative D ± according to the ROC: (+) for | z | > 1 and (−) for | z | < 1 .
Theorem 6. 
Using the stretching and shrinking derivatives, we obtain
D z ± n X ( z ) = k = ( k ) n x ( k h ) z k ,
for | z | ± 1 > 1 . Concerning the inverse ZT, we use the nabla derivative to obtain
n x ( k h ) = 1 2 π i γ X ( z ) 1 z 1 h n z k 1 d z
For both, the proofs are immediate.
Definition 11. 
For functions that have an ROC including the unit circle or for functions having a degenerate ROC, as it is the case of the periodic signals, it is preferable to work with the discrete-time Fourier transform that can be obtained from the ZT through the transformation z = e i ω h , ω h < π
X ( e i ω ) = k = x ( k h ) e i ω h k ,
with the inversion integral
x ( k ) = 1 2 π h π / h π / h X ( e i ω ) e i ω h k d ω .
To obtain the duality, we must note that we made an exponential transformation to pass from (36) to (40). Therefore, the Liouville derivative must be used. We can state:
Theorem 7. 
X ( n ) ( e i ω ) = k = ( i k h ) n x ( k h ) e i ω h k
and
n x ( k h ) = 1 2 π h π / h π / h 1 e i ω h h n X ( e i ω ) e i ω h k d ω .
Remark 3. 
The change from the nabla derivative stated in (13) to the corresponding bilinear implies a change of the factor 1 z 1 h by 2 h 1 z 1 1 + z 1 so that we obtain alternative derivative properties in (39) and (43).

4. Main Transforms and Non Integer-Order Derivatives

4.1. Laplace Transform

We reproduce here the results stated in Theorem 4
F ( n ) ( s ) = ( t ) n f ( t ) e s t d t , s C
and
h ( n ) ( t ) = 1 2 π i a i a + i F ( s ) s n e s t d s , t R .
As it is clear, the substitution n α creates problems, since the complex variable expression s α is no longer a function, it involves setting a branchcut line. The simplest way is to choose the negative or positive real semi-axis. The results in Section 2 allow us to write:
  • Right function case ( R e ( s ) > 0 )
    F ( α ) ( s ) = 0 ( t ) α f ( t ) e s t d t ,
    and
    f + ( α ) ( t ) = 1 2 π i a i a + i F ( s ) s α e s t d s , a > 0 .
  • Left function case ( R e ( s ) < 0 )
    F + ( α ) ( s ) = 0 ( t ) α f ( t ) e s t d t ,
    and
    f ( α ) ( t ) = 1 2 π i a i a + i F ( s ) s α e s t d s , a < 0 .
  • Two-sided function case ( | R e ( s ) | < b )
    The above relations suggest we introduce the two-sided fractional derivative [28]:
    D θ α f ( t ) : = lim h 0 + h α n = + ( 1 ) n Γ α + 1 Γ α + θ 2 n + 1 Γ α θ 2 + n + 1 f ( t n h ) ,
    for which
    D θ α e ± i ω t = | ω | α e i θ π 2 sgn ( ω ) e ± i ω t .
    Assuming that | θ | | α | , then
    D θ α F ( i ω ) = | t | α e i θ π 2 sgn ( t ) f ( t ) e i ω t d t ,
    and
    D θ α f ( t ) = 1 2 π | ω | α e i θ π 2 sgn ( ω ) F ( i ω ) e i ω t d ω .
    In particular, we obtain
    D 0 α F ( i ω ) = | t | α f ( t ) e i ω t d t ,
    and
    D 0 α f ( t ) = 1 2 π | ω | α F ( i ω ) e i ω t d ω .

4.2. Mellin Transform

The MT case is very similar to the LT case. The relationships are easily obtained from those in the previous subsection, taking into account the situations in which the Liouville derivative is replaced by the Hadamard derivative. The results in Section 2 allow us to write:
  • Stretching case ( R e ( v ) > 0 )
    G ( α ) ( v ) = 1 ( ln τ ) α g ( τ ) τ v 1 d τ ,
    and
    D τ + α g ( τ ) = 1 2 π i a i a + i v α G ( v ) τ v d s , a > 0 .
  • Shrinking case ( R e ( v ) < 0 )
    G + ( α ) ( v ) = 0 1 ( ln τ ) α g ( τ ) τ v 1 d τ ,
    and
    D τ α g ( τ ) = 1 2 π i a i a + i v α G ( v ) τ v d v , a < 0 .
  • Bilateral scale case ( | R e ( v ) | < b )
    The above relations suggest we introduce the two-sided scale derivative by:
    D τ , θ α g ( τ ) = lim q 1 + ln ( q ) α n = + ( 1 ) n Γ α + 1 Γ α + θ 2 n + 1 Γ α θ 2 + n + 1 g ( τ q n ) ,
    for which
    D τ , θ α τ i ω = | ω | α e i θ π 2 sgn ( ω ) τ i ω .
    Assuming that | θ | | α | , then
    D θ α G ( i ω ) = 0 | ln τ | α e i θ π 2 sgn ( ln τ ) g ( τ ) τ i ω 1 d τ ,
    and
    D τ , θ α g ( τ ) = 1 2 π | ω | α e i θ π 2 sgn ( ω ) G ( i ω ) τ i ω t d ω .
    In particular, we obtain
    D 0 α G ( i ω ) = 0 | ln τ | α g ( τ ) τ i ω 1 d τ ,
    and
    D 0 α g ( t ) = 1 2 π | ω | α G ( i ω ) e i ω t d ω .

4.3. Z and Discrete-Time Fourier Transforms

Theorem 6 expresses the duality of the ZT for the integer-order case. The situation here is similar to the one we found in the LT and MT cases, having to consider separately the three cases, corresponding to the exterior of the unit circle, | z | > 1 , the unity disk, | z | < 1 , and the annulus containing the unit circle.
  • Right sequence case ( | z | > 1 )
    D z + α X ( z ) = k = 0 ( k ) α x ( k h ) z k
    and
    α x ( k h ) = 1 2 π i γ X ( z ) 1 z 1 h α z k 1 d z ,
    where γ is a closed path lying outside the unity circle.
  • Left sequence case ( | z | < 1 )
    D z α X ( z ) = k = 1 ( k ) α x ( k h ) z k
    and
    α x ( k h ) = 1 2 π i γ X ( z ) 1 z 1 h α z k 1 d z ,
    where γ is now inside the unity circle.
  • Two-sided function case ( r < | z | < r + )
    Definition 12. 
    As above, attending to the relation [29,30]:
    ( 1 z ) a ( 1 z 1 ) b = n = + ( 1 ) n Γ a + b + 1 Γ a n + 1 Γ b + n + 1 z n ,
    and considering again two real parameters α, the derivative order, and θ, the asymmetry parameter, such that α > 1 if θ ± α , or α R if θ = ± α , we define a discrete-time two-sided derivative by
    Θ θ α x ( k h ) : = h α n = + ( 1 ) n Γ α + 1 Γ α + θ 2 n + 1 Γ α θ 2 + n + 1 x ( k h n h ) .
    Theorem 8. 
    Let x ( k h ) = e i ω h k . Then, [29]
    Θ θ α e i ω h k = 2 sin ( ω h / 2 ) α e i θ π 2 sgn ( ω ) e i ω h k .
    The proof comes from the left side in (66) by letting z = e i ω h k .
    Assuming that | θ | | α | , then
    Theorem 9. 
    D θ α X ( e i ω ) = k = | k h | α e i θ π 2 sgn ( k ) x ( k h ) e i ω h k
    and
    Θ θ α x ( k h ) = 1 2 π h π / h π / h 2 sin ( ω h / 2 ) α e i θ π 2 sgn ( ω ) X ( e i ω ) e i ω h k d ω .
    If h is very small, 2 sin ( ω h / 2 ) ω h , making (70) more similar to (69).
    In particular, we obtain
    D 0 α X ( e i ω ) = k = | k h | α x ( k h ) e i ω h k
    and
    Θ 0 α x ( k h ) = 1 2 π h π / h π / h 2 sin ( ω h / 2 ) α e i ω h k d ω .
Remark 4. 
As stated in the previous section, the change from the nabla derivative stated in (13) to the corresponding bilinear implies to substitute the factor 1 z 1 h by 2 h 1 z 1 1 + z 1 so that we obtain alternative derivative properties in (63) and (70). Relative to (70), the substitution consists of 2 sin ( ω h / 2 ) tan ( ω h / 2 ) .

5. Some Consequences

The above results allow us to obtain some consequences that express a full coherence with known theory. We start by considering the case of the fractional derivative in the complex plane [30]. Return back to (44) and (45), where we treat the right-handed functions; the corresponding left case is similar, so we will not deal with it here.
With simple substitutions, we obtain:
  • Substitute the inverse LT (30) in (44) and note that the LT converges uniformly in the ROC to commute the integrations. Using the definition of the gamma function [31], we obtain
    F ( α ) ( s ) = ( 1 ) α Γ ( α + 1 ) 2 π i a i a + i F ( u ) ( s u ) α 1 d u , R e ( s > 0 ) .
    This last formula is another way of expressing the left derivative in the complex plane, suitable for dealing with LT, but in agreement with previous formulations [30,32,33,34]. It is important to highlight an interesting fact: (73) is defined only in the right-handed complex plane. In general, we do not know what happens in the left half plane since it is out of the ROC. We profit on this fact to define there a branchcut line implicit in the definition. For this reason, we choose ( 1 ) α = e i α π .
  • Let f ( t ) = ε ( t ) be the Heaviside unit step. Insert it into (44) and use again the definition of gamma function. Assume that α > 0 . We arrive at
    F ( α ) ( s ) = ( 1 ) α Γ ( α + 1 ) s α + 1 , R e ( s ) > 0 .
    From this relation, we conclude that
    L t α Γ ( α + 1 ) ε ( t ) = 1 s α + 1 , R e ( s ) > 0 ,
    and, as L ε ( t ) = 1 / s ,
    D α 1 s = ( 1 ) α Γ ( α + 1 ) s α + 1 , R e ( s ) > 0 ,
    from where we deduce the interesting result [34]
    D α s β = ( 1 ) α Γ ( α + β ) Γ ( β ) s β α , R e ( s ) > 0 ,
    where β > 0 .
We established a correspondence between two sequences of powers, positive in time ( α > 0 ), negative in transform domain:
t α Γ ( α + 1 ) ε ( t ) s α 1 , R e ( s ) > 0 ,
and we ask ourselves what happens when α < 0 . The answer can be searched in (45). However, we must be careful since we are going to enter in the field of the generalized function [35,36]. In fact, letting F ( s ) = 1 , we obtain the formal relation
D + α δ ( t ) = 1 2 π i a i a + i s α e s t d s , a > 0 .
Let us try another way [37] by calculating the derivative of the Heaviside function using the Grünwald–Letnikov definition (5). Starting from ([20]):
k = 0 n α k ( 1 ) k = α 1 n ( 1 ) n = 1 Γ ( 1 α ) Γ ( α + n + 1 ) Γ ( n + 1 ) ,
and setting n = t / h , with h 0 + , we can show that
D α ε ( t ) = t α Γ ( 1 α ) ε ( t ) .
Remark 5. 
An important fact must be highlighted: The Caputo derivative of the unit step is null. This implies the same for the derivatives we are going to compute in the following.
As δ ( t ) = D ε ( t ) , we deduce using (79) that
D α δ ( t ) = t α 1 Γ ( α + 1 ) ε ( t ) .
Attending to (78), we conclude that
L t α 1 Γ ( α + 1 ) ε ( t ) = s α ,
which suggests the inversion of the exponent roles relative to (77)
t α 1 Γ ( α + 1 ) ε ( t ) s α , R e ( s ) > 0 .
To test the coherence of this result and the agreement with (78), we set α = n , n N . As it is easy to show,
L D n δ ( t ) = s n .
Therefore, we must have
t n 1 Γ ( n + 1 ) = D n δ ( t ) .
As strange as this relationship may seem, it was deduced by Gelf and Shilov [35].
With (79), it is possible to obtain the derivative of any order of the continuous function p ( t ) = t β u ( t ) , with β > 0 . The L T of p ( t ) yields P ( s ) = Γ ( β + 1 ) s β + 1 , for R e ( s ) > 0 and the FD of order α is given by s α Γ ( β + 1 ) s β + 1 . Therefore, the expression
D + α t β ε ( t ) = Γ ( β + 1 ) Γ ( β α + 1 ) t β α ε ( t )
is valid for any α R [38].
These results are readily reproduced for the Mellin transform. We only have to remember that we pass from the LT to the MT with the substitution e t = τ . Therefore, it is enough to perform such substitution and set v s .
The above presented application to the LT can be repeated with the Z transform. We will consider the | z | > 1 , case, since the other is similar.
Take relation (62) and substitute there x ( k h ) by its expression in terms of the inverse ZT and compute the summation and integral to obtain
D z + α X ( z ) = 1 2 π i γ X ( u ) k = 0 ( k ) α ( z / u ) k d u u .
If we introduce the kernel (a discrete Mellin transform),
K ( u ) = k = 0 ( k ) α u k ,
we can write
D z + α X ( z ) = 1 2 π i γ X ( u ) K ( z / u ) d u u .
This expression represents a new scale-derivative defined in the complex plane. However, it has one major difficulty: the lack of a closed form for K ( z ) . The same is true in most situations that we can obtain from (62) or (63), which prevents us from obtaining simple expressions similar to those obtained in the case of LT. In any case, we can fall back on the Fourier transform and implement the expressions using the fast Fourier transform.

6. Conclusions

The duality property of the Laplace and Fourier transforms associated with the integer-order derivative was reviewed. We generalized it for fractional derivatives and extended the results to the Mellin, Z and discrete-time Fourier transforms. To do it, we used the scale and discrete derivatives. Some consequences of the proposed theory are described.

Author Contributions

Conceptualization, M.D.O.; methodology, M.D.O. and G.B.; formal analysis, M.D.O.; investigation, M.D.O. and G.B.; writing—original draft preparation, M.D.O.; writing—review and editing, M.D.O. and G.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. Deakin, M.A. The development of the Laplace transform, 1737–1937: I. Euler to Spitzer, 1737–1880. Arch. Hist. Exact Sci. 1981, 25, 343–390. [Google Scholar] [CrossRef]
  2. Bateman, H. Tables of Integral Transforms; McGraw-Hill Book Company: New York, NY, USA, 1954; Volume 1. [Google Scholar]
  3. Erdélyi, A.; Magnus, W.; Oberhettinger, F.; Tricomi, F.G. Tables of Integral Transforms; McGraw-Hill Book Company: New York, NY, USA, 1954; Volume 2. [Google Scholar]
  4. Poularikas, A.D. Handbook of Formulas and Tables for Signal Processing; CRC Press: Boca Raton, FL, USA, 1999. [Google Scholar]
  5. Bracewell, R.N. The Fourier Transform and Its Applications; McGraw-Hil Higher Education: New York, NY, USA, 2000. [Google Scholar]
  6. Bailey, D.H.; Swarztrauber, P.N. The fractional Fourier transform and applications. SIAM Rev. 1991, 33, 389–404. [Google Scholar] [CrossRef]
  7. Almeida, L.B. The fractional Fourier transform and time-frequency representations. IEEE Trans. Signal Process. 1994, 42, 3084–3091. [Google Scholar] [CrossRef]
  8. 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]
  9. Daubechies, I. The wavelet transform, time-frequency localization and signal analysis. IEEE Trans. Inf. Theory 1990, 36, 961–1005. [Google Scholar] [CrossRef]
  10. Aruldoss, R.; Devi, R.A. A generalized fractional integral transform with exponential type kernel. Malaya J. Mat. (MJM) 2020, 8, 544–550. [Google Scholar] [CrossRef] [PubMed]
  11. Proakis, J.G.; Manolakis, D.G. Digital Signal Processing: Principles, Algorithms, and Applications; Prentice Hall: Upper Saddle River, NJ, USA, 2007. [Google Scholar]
  12. Oppenheim, A.V.; Willsky, A.S.; Hamid, S. Signals and Systems, 2nd ed.; Prentice-Hall: Upper Saddle River, NJ, USA, 1997. [Google Scholar]
  13. Oppenheim, A.V.; Schafer, R.W. Discrete-Time Signal Processing, 3rd ed.; Prentice Hall Press: Upper Saddle River, NJ, USA, 2009. [Google Scholar]
  14. Bertrand, J.; Bertrand, P.; Ovarlez, J. The Mellin Transform. In The Transforms and Applications Handbook, 2nd ed.; Poularikas, A.D., Ed.; Taylor Francis Group: Abingdon, UK, 2000. [Google Scholar]
  15. Leibowitz, L. A simplified binary arithmetic for the Fermat number transform. IEEE Trans. Acoust. Speech Signal Process. 1976, 24, 356–359. [Google Scholar] [CrossRef]
  16. Roberts, M. Signals and Systems: Analysis Using Transform Methods and Matlab, 2nd ed.; McGraw-Hill: New York, NY, USA, 2003. [Google Scholar]
  17. Ortigueira, M.D. Discrete-Time Fractional Difference Calculus: Origins, Evolutions, and New Formalisms. Fractal Fract. 2023, 7, 502. [Google Scholar] [CrossRef]
  18. Ortigueira, M.D.; Bohannan, G.W. Fractional scale calculus: Hadamard vs. Liouville. Fractal Fract. 2023, 7, 296. [Google Scholar] [CrossRef]
  19. Bertrand, J.; Bertrand, P.; Ovarlez, J.P. Discrete Mellin transform for signal analysis. In Proceedings of the International Conference on Acoustics, Speech, and Signal Processing, Albuquerque, NM, USA, 3–6 April 1990; pp. 1603–1606. [Google Scholar]
  20. Samko, S.; Kilbas, A.; Marichev, O. Fractional Integrals and Derivatives: Theory and Applications; Gordon and Breach Science Publishers: Amsterdam, The Netherlands, 1993. [Google Scholar]
  21. Valério, D.; Ortigueira, M.D.; Lopes, A.M. How many fractional derivatives are there? Mathematics 2022, 10, 737. [Google Scholar] [CrossRef]
  22. Ortigueira, M.D.; Valério, D. Fractional Signals and Systems; De Gruyter: Berlin, Germany; Boston, MA, USA, 2020. [Google Scholar]
  23. Bohner, M.; Peterson, A. Dynamic Equations on Time Scales: An Introduction with Applications; Springer Science & Business Media: New York, NY, USA, 2001. [Google Scholar]
  24. Ortigueira, M.D.; Coito, F.J.; Trujillo, J.J. Discrete-time differential systems. Signal Process. 2015, 107, 198–217. [Google Scholar] [CrossRef]
  25. Ortigueira, M.D.; Machado, J.T. Revisiting the 1D and 2D Laplace transforms. Mathematics 2020, 8, 1330. [Google Scholar] [CrossRef]
  26. van der Pol, B.; Bremmer, H. Operational Calculus: Based on the Two-Sided Laplace Integral; Cambridge University Press: Cambridge, UK, 1950. [Google Scholar]
  27. Dayal, S.; Singh, M.K. An analysis of convergence of Bi-lateral Laplace Transform. Int. J. Math. Its Appl. 2017, 5, 223–229. [Google Scholar]
  28. Ortigueira, M.D. Two-sided and regularised Riesz-Feller derivatives. In Mathematical Methods in the Applied Sciences; Wiley: Hoboken, NJ, USA, 2019. [Google Scholar] [CrossRef]
  29. Ortigueira, M.D. Fractional central differences and derivatives. IFAC Proc. Vol. 2006, 39, 58–63. [Google Scholar] [CrossRef]
  30. Ortigueira, M.D. Fractional Calculus for Scientists and Engineers; Lecture Notes in Electrical Engineering; Springer: Berlin/Heidelberg, Germany, 2011. [Google Scholar]
  31. Henrici, P. Applied and Computational Complex Analysis; Wiley-Interscience: New York, NY, USA, 1991; Volume 2. [Google Scholar]
  32. Campos, L. On the branchpoint operator and the annihilation of differintegrations. SIAM J. Math. Anal. 1989, 20, 439–453. [Google Scholar] [CrossRef]
  33. Campos, L. On a generalized Mittag–Leffler theorem and implicit differintegration. SIAM J. Math. Anal. 1989, 20, 454–467. [Google Scholar] [CrossRef]
  34. Ortigueira, M.D. A coherent approach to non-integer order derivatives. Signal Process. 2006, 86, 2505–2515. [Google Scholar] [CrossRef]
  35. Gelfand, I.M.; Shilov, G.P. Generalized Functions; Academic Press: New York, NY, USA, 1964; Volume 1. [Google Scholar]
  36. Zemanian, A.H. Distribution Theory and Transform Analysis: An Introduction to Generalized Functions, with Applications; Lecture Notes in Electrical Engineering, 84; Dover Publications: New York, NY, USA, 1987. [Google Scholar]
  37. Ortigueira, M.; Machado, J. Which derivative? Fractal Fract. 2017, 1, 3. [Google Scholar] [CrossRef]
  38. Kilbas, A.; Srivastava, H.; Trujillo, J. Theory and Applications of Fractional Differential Equations; North-Holland Mathematics Studies; Elsevier: Amsterdam, The Netherlands, 2006; Volume 204. [Google Scholar]
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.; Bengochea, G. On the Fractional Derivative Duality in Some Transforms. Mathematics 2023, 11, 4464. https://doi.org/10.3390/math11214464

AMA Style

Ortigueira MD, Bengochea G. On the Fractional Derivative Duality in Some Transforms. Mathematics. 2023; 11(21):4464. https://doi.org/10.3390/math11214464

Chicago/Turabian Style

Ortigueira, Manuel Duarte, and Gabriel Bengochea. 2023. "On the Fractional Derivative Duality in Some Transforms" Mathematics 11, no. 21: 4464. https://doi.org/10.3390/math11214464

APA Style

Ortigueira, M. D., & Bengochea, G. (2023). On the Fractional Derivative Duality in Some Transforms. Mathematics, 11(21), 4464. https://doi.org/10.3390/math11214464

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop