Next Article in Journal
Neural Network-Based Hybrid Forecasting Models for Time-Varying Passenger Flow of Intercity High-Speed Railways
Next Article in Special Issue
A Generalized Lyapunov Inequality for a Pantograph Boundary Value Problem Involving a Variable Order Hadamard Fractional Derivative
Previous Article in Journal
Plastic-Pollution Reduction and Bio-Resources Preservation Using Green-Packaging Game Coopetition
Previous Article in Special Issue
Fractional Systems’ Identification Based on Implicit Modulating Functions
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Unilateral Laplace Transforms on Time Scales

by
Müfit Şan
1 and
Manuel D. Ortigueira
2,*
1
Department of Mathematics, Faculty of Science, Çankırı Karatekin University, Çankırı 18100, Türkiye
2
Centre of Technology and Systems—UNINOVA and Department of Electrical Engineering, NOVA Faculty of Science and Technology of NOVA University of Lisbon, Quinta da Torre, 2829-516 Caparica, Portugal
*
Author to whom correspondence should be addressed.
Mathematics 2022, 10(23), 4552; https://doi.org/10.3390/math10234552
Submission received: 4 November 2022 / Revised: 23 November 2022 / Accepted: 25 November 2022 / Published: 1 December 2022
(This article belongs to the Special Issue Fractional Differential Equations and Control Problems)

Abstract

:
We review the direct and inverse Laplace transforms on non-uniform time scales. We introduce full backward-compatible unilateral Laplace transforms and studied their properties. We also present the corresponding inverse integrals and some examples.

1. Introduction

Mathematical developments do not solely depend on internal dynamics. The motivational effects of problems in science and engineering on these development processes may be fundamental. The attempts to solve such problems have triggered many new concepts in mathematics. Whenever notions and tools are not sufficient/appropriate to solve a given problem or do not yield the desired/fruitful results, we need to search for new ones to proceed to the generalization of the existing one. However, when performing such a generalization, it must satisfy the backward compatibility with the previously existing methodology. Here, we present an important example.
The classical theory of a differential/difference equation involves maintaining a parallelism between the formulation, solution, and properties; this implies great similarity and allows the establishment of equivalent rules [1]. These types of systems have uniform domains. However, some situations arise where we may have discrete or mixed discrete/continuous irregular domains in which classic ARMA-type equations are not suitable. To overcome such a discrepancy and allow obtaining similar results, Aulbach and Hilger [2] introduced a framework that unifies and extends the continuous and discrete existing methodologies: the “Calculus on Time Scales”.
A time scale, denoted by T , is defined as an arbitrary non-empty closed subset of the real line R and can be classified as uniform or non-uniform (or, irregular). Usually, uniform time scales are used [3,4], but non-uniform ones stand out in many applications, as pointed out in [5,6]. To solve the differential–difference equations corresponding to applications with underlying non-uniform time scales; it is obvious that the classic Laplace and Z transforms are not suitable and must be redefined for such time scales. A first generalization of the one-sided Laplace transform (LT) was proposed by Bohner and Peterson [7,8,9,10,11,12,13,14]. A second approach to such a goal was pursued by Davis et al. [8,15], where the domain of the proposed generalized LT was taken as the global time scale (bilateral transforms). However, such papers deal with transforms that are not exactly backward compatible, because they do not recover the Z transform when the time scale is discrete uniform. Moreover, they do not introduce an inverse transform. To overcome these difficulties, Ortigueira et al [4,6] proposed a different approach. They began by computing the exponentials, eigenfunctions of the nabla (causal) and delta (anti-causal) derivatives, and, from them, defined an inverse Laplace transform that was backward-compatible with the Bromwich integral for the classic Laplace transform (LT), as well as the Cauchy integral for the Z transform (ZT). Using the properties of the referred exponentials, they obtained nabla and delta bilateral LTs that were fully compatible with the classic LT and ZT.
In this work, following [6], we present unilateral nabla and delta Laplace transforms on non-uniform time scales. Unlike the unilateral delta LT introduced by Bohner and Peterson (see [16], page 118), both transforms will coincide with the Z transform if the time scale is Z 0 + . Moreover, we deduce the properties enjoyed by classical one-sided LT and ZT and present the fully compatible characteristics.
This new transform allows us to solve problems defined on irregular domains and implement control systems when the sampling is not uniform. Moreover, it will permit the introduction of fractional derivatives on irregular time scales by inverting the transform s α and using the convolution theorem.
The study is organized as follows. Section 2 introduces basic definitions of time scales, including graininess, jump operators, nabla and delta derivatives/anti-derivatives, and corresponding nabla and delta exponentials. Section 3 presents unilateral nabla and delta Laplace transforms on non-uniform time scales and their properties, and show the referred compatibility. We present some examples. Finally, in the last section we present the conclusions.

2. On-Time Scale Calculus

2.1. Basic Definitions

Let t T be the current time instance. Let ρ ( t ) and σ ( t ) be positive real-valued functions representing the previous and next instances, respectively. With n t h , the power of the ρ ( t ) and σ ( t ) , i.e., ρ n ( t ) and σ n ( t ) , we refer to the n t h previous and n t h next instances from the current one t , respectively, defined by the recursions (see Figure 1).
ρ 0 ( t ) : = t , ρ 1 ( t ) : = ρ ( t ) , ρ n ( t ) : = ρ ( ρ n 1 ( t ) ) and σ 0 ( t ) : = t , σ 1 ( t ) : = σ ( t ) , σ n ( t ) : = σ ( σ n 1 ( t ) )
where n 2 . These functions are useful for specifying the steps needed to go from one instance to another, allowing us to introduce the conceptual differences between two instances and later the translation used in the convolution.
These functions allow us to introduce the graininess functions ν ( t ) : = t ρ ( t ) and μ ( t ) : = σ ( t ) t . The n t h power, i.e., ν n ( t ) or μ n ( t ) , expresses the value of the difference between the current instance t and the n t h previous one ρ n ( t ) , as well as the difference between the current instance t and the n t h following one σ n ( t ) , respectively. They can be recursively generated according to:
μ n ( t ) = μ n 1 ( t ) + μ ( t + μ n 1 ( t ) ) , n = 1 , 2 ,
with μ 0 ( t ) : = 0 , and
ν n ( t ) = ν n 1 ( t ) + ν ( t ν n 1 ( t ) ) , n = 1 , 2 ,
where ν 0 ( t ) : = 0 . In the following, we can see an illustration of the recursions.
σ 0 ( t ) = t = t + μ 0 ( t ) , σ 1 ( t ) = σ ( t ) = t + μ ( t ) = t + μ 1 ( t ) , σ 2 ( t ) = σ ( σ ( t ) ) = σ ( t ) + μ ( σ ( t ) ) = t + μ 1 ( t ) + μ 1 ( t ) = t + μ 2 ( t ) , . . . σ n ( t ) = t + μ n ( t ) ,
moreover, analogously, the relation ρ n ( t ) = t ν n ( t ) can be obtained. At this stage, one may wonder about the meaning of a difference, such as t τ for t , τ T , and how to define it, so that it makes sense. Let t > τ . This means that one can reach t by moving from the instance τ into the future; otherwise, one can arrive at τ by moving from the instance t into the past. In a more formal way, there exists N N , such that t = σ N ( τ ) or τ = ρ N ( t ) . Set
ν n ( t ) = ρ n 1 ( t ) ρ n ( t ) and μ n ( t ) = σ n ( t ) σ n 1 ( t )
for n 1 . Then, one can write the difference t τ in the following two ways:
t τ = σ N ( τ ) τ = μ N ( τ ) = μ ( τ ) + μ ( τ + μ ( τ ) ) + = n = 1 N μ n ( τ )
or
t τ = t ρ N ( t ) = ν N ( t ) = ν ( t ) + ν ( t ν ( t ) ) + = n = 1 N ν n ( t ) .
In this study, the time scales that we deal with are assumed as a union of two sets, i.e., isolated points and closed intervals. As we will see in Section 2.4, with the isolated points and boundary points of each interval, we construct a discrete time scale T = t n : n Z . Each point T = t n : n Z is right dense, isolated, or left dense. For example, if t k is a left (or right) dense point, then it represents all points t k 1 + h (or t k + 1 h ), although we have to insert a limited computation when computing any derivative or integral. With this formulation, we can redefine direct and reverse graininess in a way that plays a vital role in representing the integral function in the summation forms. The direct graininess is given by ν n = ν ( t n ) = t n t n 1 , n Z , and the reverse graininess is defined by μ n = μ ( t n ) = t n + 1 t n , n Z , by avoiding representing the reference instance t 0 . In the following section, the nabla and delta exponentials will be formulated in terms of these two types of graininess.

2.2. Nabla and Delta Derivatives

Let T be a given time scale in which will introduce two derivatives.
Definition 1.
The nabla derivative is a causal operator defined by
f ( t ) = f ( t ) f ( ρ ( t ) ) ν ( t ) if ν ( t ) 0 lim h 0 + f ( t ) f ( t h ) h if ν ( t ) = 0 ,
while the delta derivative is an anti-causal operator given by
f Δ ( t ) = f ( σ ( t ) ) f ( t ) μ ( t ) if μ ( t ) 0 lim h 0 + f ( t + h ) f ( t ) h if μ ( t ) = 0 .
If we consider the time scales T = t n : n Z specialized in the previous section, consisting of the union of a numerable discrete set t n : n Z and a sequence of closed intervals, then these above definitions assume slightly different forms
f ( t n ) = f ( t n ) f ( t n 1 ) ν n for any t n that is not left dense lim h 0 + f ( t n ) f ( t n h ) h for any t n that is left dense ,
where ν n = t n t n 1 and the delta derivative is given by
f Δ ( t n ) = f ( t n + 1 ) f ( t n ) μ n for any t n that is not right dense lim h 0 + f ( t n + h ) f ( t n ) h for any t n that is right dense
where μ n = t n + 1 t n .
Example 1.
Consider a time scale defined by T , such that t n = n + r ( n ) , n = 1 , 2 , , L . The sequence r ( n ) , n = 1 , 2 , , L , is randomly uniformly distributed in the interval ( 0.5 , 0.5 ) . With this sequence, we sampled a sinusoid to obtain a signal f ( t n ) = c o s ( 2 π T t n + ϕ ) , where ϕ is a random initial phase, L = 250 , and T = 50 . In Figure 2, we depict the nabla derivative of f ( t n ) (to make the visualization easier, we used an interpolation on the plot).
The delta derivative gives a similar function.

2.3. The Nabla and Delta Anti-Derivatives

The relations (3)–(4) and (5)–(6) can be inverted to give the corresponding anti-derivatives [6]. The inverses of nabla and delta derivatives of (3) and (4) are given by
f 1 ( t ) = n = 0 ν n + 1 ( t ) f ( t ν n ( t ) ) ,
where ν n ( t ) was introduced in (2), and
f Δ 1 ( t ) = n = 0 μ n + 1 ( t ) f ( t + μ n ( t ) ) ,
where μ n ( t ) was given in (1). Here, f 1 ( ) = f Δ 1 ( + ) = 0 is assumed.
These anti-derivatives allow us to introduce two definite nabla and delta integrals over [ a , b ] , respectively, by
a b f ( t ) t and a b f ( t ) Δ t .
If T consists of the isolated points, then there exists a j N , such that b = σ j ( a ) = a + μ j ( a ) or a = ρ j ( b ) = b ν j ( b ) . By considering (7) we will write
a b f ( t ) t = n = 0 j 1 ν n + 1 ( b ) f ( b ν n ( b ) ) ,
The lower terminal point of the integral will be taken as a , even if the summation on the right-hand side does not involve a multiplier with f ( a ) and involves f ( σ ( a ) ) in the last multiplier. When we pass from the integral to the summation, we always keep this remark in mind. If f 1 is the nabla inverse of the function f , then the Equation (9) can be written as follows
a b f ( t ) t = n = 0 j 1 ν n + 1 ( b ) f ( b ν n ( b ) ) = f 1 ( b ) f 1 ( a ) ,
which expresses the nabla definite integral in a Barrow-like formula.
In the next section, we will define unilateral (right and left) -LTs by splitting the integral representing the bilateral LT at the reference point t 0 . Moreover, we will explain the relationship between -LT and -ZT. Therefore, it will be more useful to formulate the nabla integral on a discrete set a = t 0 , t 1 , , b = t j as follows:
a b f ( t ) t = t 0 t j f ( t ) t = n = 0 j 1 ν n + 1 f ( t n + 1 ) = n = 1 j ν n f ( t n )
where ν n = t n t n 1 for n = 1 , , j 1 . and the summation do not involve ν 0 f ( t 0 ) .
Analogously, from (8), we will write
a b f ( t ) Δ t = n = 0 j 1 μ n + 1 ( a ) f ( a + μ n ( a ) ) .
The upper terminal point of the integral will be assumed as b even if the summation on the right-hand side does not involve a multiplier with f ( b ) . Nevertheless, the last sum in this summation contains μ j ( a ) f ( ρ ( b ) ) as a last multiplier. When we pass from the integral to the summation, we always recall this note. If f Δ 1 is the delta inverse of the function f , then Equation (12) can be written as follows
a b f ( t ) Δ t = n = 0 j 1 μ n + 1 ( a ) f ( a + μ n ( a ) ) = f Δ 1 ( b ) f Δ 1 ( a ) ,
that resembles again the classic Barrow formula.
For the reasons we explained above for the nabla case, the formulation (13) for the discrete set a = t 0 , t 1 , , b = t j can be given
a b f ( t ) Δ t = n = 0 j 1 μ ( t n ) f ( t n ) = n = 0 j 1 μ n f ( t n ) ,
where μ n = μ ( t n ) = t n + 1 t n for all n = 0 , 1 , , j 1 . This sum already involves μ ( t 0 ) f ( t 0 ) as a parcel. Moreover, it is known that [16]
f ( t ) g ( t ) = f ( t ) g ( t ) + f ( ρ ( t ) ) g ( t ) 1
so that we can obtain the integration of the parts formula for the nabla derivative:
Theorem 1.
Let f and g be -differentiable functions, then
a b f ( t ) g ( t ) t = f ( b ) g ( b ) f ( a ) g ( a ) a b f ( ρ ( t ) ) g ( t ) t .
is fulfilled.
Similarly, by using (13), one can obtain the integration of the parts rule for the delta derivative
f ( t ) g ( t ) = f Δ ( t ) g ( t ) + f ( σ ( t ) ) g Δ ( t ) Δ 1 ,
which leads us to:
Theorem 2.
Let f and g be Δ -differentiable functions, then
a b f Δ ( t ) g ( t ) Δ t = f ( b ) g ( b ) f ( a ) g ( a ) a b f ( σ ( t ) ) g Δ ( t ) Δ t
is satisfied.
Relations (15) and (16) originate from the initial conditions in the nabla and delta Laplace transforms, respectively.

2.4. The Nabla and Delta Unit Step Functions

We define the nabla and delta unit step functions, respectively, as follows [4]:
u t , t 0 = 1 if t = t n t 0 0 if t = t n < t 0 ,
and
u Δ t , t 0 = 0 if t = t n > t 0 1 if t = t n t 0 .
It is obvious that u t , t 0 = u Δ t , t 0 . Moreover, the derivatives of these functions give the nabla and delta impulse functions:
D [ u t , t 0 ] = δ t , t 0 = 1 v 1 t 0 ( = 1 v 0 ) if t = t n = t 0 0 if t = t n t 0 .
and
D Δ [ u Δ t , t 0 ] = δ Δ t , t 0 = 1 μ 1 t 0 ( = 1 μ 0 ) if t = t n = t 0 0 if t = t n t 0 .
If t h Z , then these two impulses coincide with the Kronecker impulses given by [4]. Moreover, if t R , is a non-isolated point, we have
δ t , t 0 = δ Δ t , t 0 = δ t t 0
where δ represents the classical Dirac distribution.

2.5. The Nabla and Delta General Exponentials

Consider a time scale as above. A recursive procedure allows us to obtain the eigenfunctions (exponentials) of nabla and delta derivatives [6].

2.5.1. The General Nabla Exponential

Let t > t 0 . Supposing f ( t ) = s f ( t ) ,
f ( t ) [ 1 s ν ( t ) ] = f ( t ν ( t ) ) ,
is obtained. If we write t = t ν ( t ) in the last equation, then we have
f ( t ν ( t ) ) [ 1 s ν ( t ν ( t ) ) ] = f ( t ν ( t ) ν ( t ν ( t ) ) .
Proceeding in the same way for n times and using the relations ν n t = ν n 1 t ν t ν n 1 t and ν ( t ) = ν 1 ( t ) , ν ( t ν ( t ) ) = ν 2 ( t ) , , ν ( t ν n 1 ( t ) ) = ν n ( t ) , we have
f ( t ) k = 1 n 1 s v k ( t ) = f t v n ( t ) = f ( t 0 ) .
By using f ( t 0 ) = 1 , ν n k + 1 = μ k ( t 0 ) and t k = t k t k 1 for k = 1 , 2 , , n in the last equation, we have the discrete and integral form of the nabla exponential
f ( t ) = k = 1 n 1 s μ k t 0 1 = e k = 1 n ln 1 s μ k ( t 0 ) μ k ( t 0 ) t k = e t 0 t l n ( 1 s μ ( τ ) ) μ ( τ ) τ
for t > t 0 .
Now let t < t 0 . By following the same way used for t > t 0 but by starting from the point t 0 , we have
f ( t 0 ) k = 1 n 1 s v k ( t 0 ) = f t 0 v n ( t 0 ) = f ( t n ) ,
where f ( t 0 ) = 1 . Applying the relation ν n k + 1 t 0 = Δ t n k + 1 for k = 1 , 2 , , n , to the last equation, for t < t 0 we have
f t n = k = 1 n 1 s ν n k + 1 t 0 = e k = 1 n ln 1 s ν n k + 1 t 0 ν n k + 1 t 0 Δ t n k + 1 = e t t 0 ln ( 1 s v ( τ ) ) v ( τ ) Δ τ .
Consequently, from (23) and (24), we have
e ( t , t 0 ; s ) = k = 1 n 1 s μ k ( t 0 ) 1 if t > t 0 1 if t = t 0 , k = 1 n 1 s ν k ( t 0 ) if t < t 0
where it is assumed that μ k 0 and ν k 0 for some k .
Moreover, taking into consideration ν n ( t 0 ) = μ n = ν n + 1 and μ n ( t 0 ) = ν n = μ n 1 for n 1 , we have the equivalent representations of the above formulation
e ( t n , t 0 ; s ) = k = 1 n 1 s ν k 1 if t n > t 0 1 if t n = t 0 , k = 0 n 1 1 s ν k if t n < t 0 , = k = 0 n 1 1 s μ k 1 if t n > t 0 1 if t n = t 0 , k = 1 n 1 s μ k if t n < t 0 .
Furthermore, from (23) and (25), the generalized nabla exponential for any time scale is as follows
e ( t , t 0 ; s ) = e x p ( t 0 t l n ( 1 s μ ( τ ) ) μ ( τ ) τ ) if t > t 0 1 if t = t 0 , e x p ( t t 0 l n ( 1 s ν ( τ ) ) ν ( τ ) Δ τ ) if t < t 0 .

2.5.2. The General Delta Exponential

For t > t 0 , assuming f Δ ( t ) = s f ( t ) , we have
f ( t + μ ( t ) ) = f ( t ) [ 1 + s μ ( t ) ] .
Set t = t 0 , then from the above equation we have
f ( t 0 + μ ( t 0 ) ) = f ( t 0 ) [ 1 + s μ ( t 0 ) ] .
Let t 0 = t 0 + μ ( t 0 ) . We have
f ( t 0 + μ ( t 0 ) + μ ( t 0 + μ ( t 0 ) ) ) = f ( t 0 + μ ( t 0 ) ) [ 1 + s μ ( t 0 + μ ( t 0 ) ) ] .
Repeating the above relation n times and using that f t 0 = 1 , we obtain
f ( t n ) = f t 0 + μ n t 0 = f t 0 k = 1 n 1 + s μ k t 0 = k = 1 n 1 + s μ k t 0 ,
where we used the relation μ n ( t 0 ) = μ n 1 ( t 0 ) + μ ( t 0 + μ n 1 ( t 0 ) ) , and μ ( t 0 + μ ( t 0 ) = μ 2 ( t 0 ) , …, μ ( t 0 + μ n 1 ( t 0 ) = μ n ( t 0 ) . By using μ k t 0 = t k = t k t k 1 , the last equation yields
f ( t n ) = k = 1 n 1 + s μ k t 0 = e k = 1 n l n ( 1 + s μ k t 0 ) μ k t 0 μ k t 0 = e t 0 t ln ( 1 + s μ ( τ ) ) μ ( τ ) τ
for t > t 0 .
Now we assume that t < t 0 , i.e., t 0 = t + μ m ( t ) . By following the same procedure applied for the case t > t 0 , we have
f ( t ) = k = 1 m 1 + s μ k t 1 .
Since μ 1 ( t ) = ν m ( t 0 ) , μ 2 ( t ) = ν m 1 ( t 0 ) , , μ m ( t ) = ν 1 ( t 0 ) , and ν m k + 1 t 0 = Δ t m k + 1 , from the last equation
f ( t ) = k = 1 m 1 + s ν m k + 1 t 0 1 = e k = 1 m l n ( 1 + s ν m k + 1 t 0 ) ν m k + 1 t 0 Δ t m k + 1 = e t t 0 ln ( 1 + s v ( τ ) ) v ( τ ) Δ τ
is obtained.
From (30) and (31), the delta exponential in the discrete form can be given as in [6]
e Δ ( t , t 0 ; s ) = k = 1 n 1 + s μ k ( t 0 ) if t > t 0 1 if t = t 0 k = 1 n 1 + s ν k ( t 0 ) 1 if t < t 0 .
where t = t 0 + μ n ( t 0 ) , when t > t 0 and t = t 0 ν n ( t 0 ) , when t < t 0 , and it is assumed that μ k 0 and ν k 0 for some k .
By using ν n ( t 0 ) = μ n = ν n + 1 and μ n ( t 0 ) = ν n = μ n 1 for n 1 , we have the equivalent representations of the above formulation
e Δ ( t n , t 0 ; s ) = k = 0 n 1 1 + s μ k if t n > t 0 1 if t n = t 0 k = 1 n 1 + s μ k 1 if t n < t 0 , = k = 1 n 1 + s ν k if t n > t 0 1 if t n = t 0 k = 0 n 1 1 + s ν k 1 if t n < t 0 .
From (30) and (32), we have the integral representations for the general delta exponential as follows:
e Δ ( t , t 0 ; s ) = e x p ( t 0 t l n ( 1 + s μ ( τ ) ) μ ( τ ) τ ) if t t 0 1 if t = t 0 , e x p ( t t 0 l n ( 1 + s ν ( τ ) ) ν ( τ ) Δ τ ) if t < t 0 .

2.5.3. Properties of the Exponentials

The above-introduced exponentials enjoy some properties that will be useful in the following section [6]. Here, we will consider those most interesting for our objectives. As it is easy to verify, there is a relation between both exponentials
e Δ ( t , t 0 ; s ) = 1 / e ( t , t 0 ; s ) .
Therefore, we do not need two exponentials. So, we will use the nabla exponential only and remove the subscript whenever it is not needed ( e ( t , t 0 ; s ) ).
In many situations, it is important to know how exponentials increase/decrease for s C . As said, consider the nabla exponential case. We start by noting that each graininess ν k defines a Hilger circle centered at its reciprocal, passing through s = 0 and being located in the right-hand half-complex plane. With μ k , the situation is similar, but the Hilger circles will be located in the left half complex plane. If we assume that h max = max ( ν k ) , k Z , h min = min ( ν k ) , and k Z , then the circle centered at 1 / h min will be the outermost one amongst all right Hilger circles while the circle centered at 1 / h max will be the innermost. Since we have infinite graininess values in the time scale, the number of corresponding circles is infinite. However, if the graininess is constant, then they reduce to one. Lastly, in the case of zero graininess, we will talk about the imaginary axis instead of circles. Considering all of these facts and the definition of the nabla exponentials e ( t , t 0 ; s ) , one can observe the following:
  • It is a real-valued function for any s R ;
  • It is a positive real-valued function for any s R , such that s < 1 / h m a x ;
  • oscillates for any s R , such that s > 1 / h m i n ;
  • It is a bounded function for s C in the innermost Hilger circle 1 s h m a x = 1 ;
  • It has an absolute value that increases as s increases outside the outermost Hilger circle 1 s h min = 1 , going to infinite as s .
  • Let T = h Z , h > 0 . We make t 0 = 0 and use ν ( t ) = μ ( t ) = h , t = n h , n Z , leading to e ( t , 0 ; s ) = 1 s h n . Using z 1 for 1 s h , we obtain the current discrete-time exponential, z n .
  • Let T = R . Return to the above case and use h = t n . As
    lim n 1 s t n n = e s t ,
    we obtain e ( t , 0 ) = e s t .
The main properties of the exponential read [6]
1.
Interchanging the role of instances.
e ( t 0 , t ; s ) = 1 / e ( t , t 0 ; s ) .
2.
Scale changing.
Let a R + . Then, the equality
e ( a t , t 0 ; s ) = e ( t , t 0 ; a s )
holds.
To prove this, we first assume that n > 0 , i.e., t > 0 . Since, ν k = ν ( t k ) = t k t k 1 for k Z , we have ν ( a t k ) = a t k a t k 1 = a ( t k t k 1 ) = a ν ( t k ) , where a t k a T . From here, we have
e a t , t 0 ; s = k = 1 n 1 s ν ( a t k ) 1 = k = 1 n 1 a s ν ( t k ) 1 = e t , t 0 ; a s ,
where e t , t 0 ; s = k = 1 n 1 s ν ( t k ) 1 as in (27) with ν k = ν ( t k ) .
Let n < 0 , i.e., t < 0 . Since ν k = ν ( t k ) = t k t k 1 for k N , we have ν ( a t k ) = a t k a t k 1 = a ( t k t k 1 ) = a ν ( t k ) . Then, we have
e a t , t 0 ; s = k = 1 n 1 1 s ν ( a t k ) 1 = k = 1 n 1 1 a s ν ( t k ) 1 = e t , t 0 ; a s ,
where e t , t 0 ; s = k = 0 n 1 1 s ν k for t < 0 as in (27) with ν k = ν ( t k ) .
3.
Product of exponentials.
The general products of exponentials may not be exponentials in the known sense, even if they are well-defined. However, the following relations are satisfied [6]:
(a)
By using the first and second properties given just above, we have
e ( t , t 0 ; s ) / e ( τ , t 0 ; s ) = e ( t , τ ; s ) , t τ .
(b)
Let t n > t m . Then, we obtain
e ( t n , t 0 ; s ) = e ( t m , t 0 ; s ) e ( t n , t m ; s ) .
4.
Shift Property.
Here, we deal with nonuniform time scales. Since t n t m may not be an element of the nonuniform time scale T , there is no guarantee of an existing mean for e ( t n t m , t j , ± s ) with t j T Instead of using t n t m to define the shift property, we use t n m , which is equal to t n t m = t τ when the time scales are uniform. Let n > m . From the definition of the (nabla) exponential in (27) we have
e ( t n m , t 0 ; s ) = k = 1 n m 1 s ν k 1 = k = 1 n 1 s ν k 1 k = n m + 1 n 1 s ν k 1 = e ( t n , t 0 ; s ) e ( t n , t n m ; s ) = e ( t n , t 0 ; s ) e ( t n m , t n ; s )
where we used property 1. For m > n , one can have the same equality in (37).

3. Laplace Transforms on Time Scales

The nabla and delta exponentials introduced in the last section allow us to define with generality two Laplace transforms defined over any time scale. The nabla and delta two-sided Laplace transforms were introduced and studied before in [6]. However, in agreement with the considerations we conducted in the previous sub-section, we will revise both and consider the corresponding one-sided transforms as well. However, we will pay attention to the -Laplace transform ( )-LT to distinguish it from the usual LT defined on R [17].

3.1. Inverse and Bilateral -Laplace Transform

In [6], Ortigueira et al. first presented the inverse nabla–Laplace transform before defining the corresponding direction.
Definition 2.
Let f ( t ) be defined on a given time scale T and assume it has a LT, F ( s ) = L f ( t ) . Then, f ( t ) can be synthesized through a continuous, infinite set of elemental exponentials with differential amplitudes:
f t = 1 2 π i C F ( s ) 1 e t ν ( t ) , t 0 ; s d s
where the integration path C , is a simple-closed contour in a region of analyticity of the integrand and encircles the poles of the delta exponential. Relation (38) is also called the synthesis equation of the LT. For this reason, the direct transform is also called the analysis formula [18].
The coherence of (38) can be verified through some tests
  • Let us take F ( s ) = 1 and calculate its inverse Laplace transform. By definition of e t , t 0 ; s in (27), one can see that the integrand is analytic when t > t 0 ; it is in the form of 1 P ( s ) (P is a polynomial of a degree greater than 1) when t < t 0 . So, its integral over C for both cases is null. When t = t 0 , we have only one pole, and
    f ( t ) = 1 2 π i C F ( s ) 1 e t ν ( t ) , t 0 ; s d s = 1 2 π i C 1 1 + s ν 1 ( t 0 ) d s = 1 ν 1 ( t 0 ) = 1 v 0 .
    So, the inverse Laplace transform of F ( s ) = 1 is equal to δ t n , t 0 where δ is a nabla impulse function given by (19). In view of (21), this result makes sense and is compatible with its continuous counterpart obtained for a non-isolated point t R .
  • Let F ( s ) = e t k , t 0 ; s , k Z . Then
    f ( t k ) = 1 2 π i C e t k , t 0 ; s 1 e t n 1 , t 0 ; s d s = 1 2 π i C e t k , t n 1 ; s d s = δ ( t k , t n )
    where the relationship between the nabla and delta exponentials [6]
    e t , t 0 ; s 1 e τ , t 0 ; s = e ( t , τ ; s )
    was used together with
    1 2 π i C e t k , t n 1 ; s d s = 0 if k n 1 ν n if k = n .
  • Consider a uniform time scale T = h Z , h R + . In this case, t = n h , n Z , ν ( t ) = h , and t 0 = 0 so that e t ν ( t ) , t 0 ; s = 1 s h n . With
    z 1 = 1 s h , we obtain the inverse Z transform;
    Setting h = t / n and letting h 0 we arrive to the Bromwich integral inverse of the usual continuous-time LT.
  • Moreover, by -differentiating on both sides of (38), we have L f ( t ) ( s ) = s L f ( ρ ( t ) ) ( s ) . This result will be confirmed later by direct transformation.
Definition 3.
Attending to the relation (40), the nabla–Laplace transform of a function f : T C is defined by
F ( s ) = L { f ( t ) } ( s ) = n = ν n f ( t n ) e t n , t 0 ; s
for those values s C for which the corresponding series converges.
Considering the relation between -integration and -summation given in (11), the above formulation can be formulated in terms of the integral, as follows:
L f ( t ) ( s ) = + f ( t ) e ( t , t 0 ; s ) t .
Similar to the above procedure, we can test the coherence of the definition by:
  • if f ( t ) = δ t n , t 0 , then from (42) and (43)
    L δ t n , t 0 ( s ) = + δ t , t 0 e ( t , t 0 ; s ) t = n = ν n δ t n , t 0 e t n , t 0 ; s = 1
    where we use δ t 0 , t 0 = 1 / ν 0 . This reproduces the above result L 1 1 = δ t n , t 0 .
Theorem 3
(Inverse -Laplace transform). Let F ( s ) be the Laplace transform of f : T C defined by (42) Then the formula in (38) represents the inverse -Laplace transform.
Proof. 
Inserting the nabla–Laplace transform formulation of f ( t ) given by (42) into (38), we have, attending to (40)
L 1 { F ( s ) } ( t n ) = 1 2 π i C k = ν k f ( t k ) e t k , t 0 ; s 1 e t n 1 , t 0 ; s d s = k = ν k f ( t k ) 1 2 π i C e t k , t 0 ; s 1 e t n 1 , t 0 ; s d s = k = ν k f ( t k ) 1 2 π i C e t k , t n 1 ; s d s = k = ν k f ( t k ) δ ( t k , t n ) = f ( t n )
for all n N , where we use the uniform convergence of the summation representing nabla–Laplace transform. □
  • Consider a uniform time scale T = h Z , h R + . In this case, t = n h , n Z , ν ( t ) = h , and t 0 = 0 so that e t ν ( t ) , 0 ; s = 1 s h n . With
    z 1 = 1 s h we obtain the usual Z transform:
    Setting h = t / n and letting h 0 we obtain the continuous-time bilateral LT.

3.2. The Unilateral -Laplace Transform

Definition 4.
We define unilateral (one-sided) -Laplace transform by:
L f ( t ) ( s ) = n = u ( ± ( t n t 0 ) ) ν n f ( t n ) e t n , t 0 ; s
where s C , t 0 is the reference point used for defining the nabla exponential and u is the nabla unit step function given by (17). The sign in the argument of u in the formulation (45) defines the left (−) and right (+) transforms.
To show the integral forms of the unilateral (one-sided) -Laplace transforms, we need to consider the relationship between the -integration and -summation in Section 2.3. According to the remark noted there, because the summation in (45) involves the multiplier v 0 f ( t 0 ) when n n 0 , the lower terminal point of the corresponding integral must be ρ ( t 0 ) . So, the integral form of the right unilateral -Laplace transform is given by
L + f ( t ) ( s ) = ρ ( t 0 ) + f ( t ) e ( t , t 0 ; s ) t .
If t R , then ρ ( t 0 ) t 0 and the above formula turns into the following
L + f ( t ) ( s ) = t 0 + f ( t ) e ( t , t 0 ; s ) t .
It means that t 0 is fully involved. In this respect, the nabla situation is different from the corresponding delta, as we will see later. It is very interesting to note that (47) agrees with a modified Laplace transform used in the control [19].
It is obvious that the left unilateral -Laplace transform is in the integral form of
L f ( t ) ( s ) = ρ ( t 0 ) f ( t ) e ( t , t 0 ; s ) t .
Theorem 4.
Let f ( t ) be a bounded function of the exponential type, i.e.,
(i) There exists A R + and a R , such that
f ( t n ) < A / e ( t n , t 0 ; a ) ( n < n 1 < 0 )
when 1 + a ν M * < 1 + s ν m , where ν M * is one of ν k so that 1 + a ν M * = max ν k 1 + a ν k , ν m = min k 1 , 2 , , n ν k with n > 0 , and
(ii) There exists B R + and b R such that
f ( t n ) < B / e ( t n , t 0 ; b ) ( n > n 2 > 0 )
when 1 + s ν M < 1 + b ν m * , where ν m * is one of ν k so that 1 + b ν m * = min ν k 1 + b ν k , ν M = max k { 0 , 1 , 2 , , n 1 } ν k with n < 0 .
If s C is inside the intersection of the circles defined by 1 + s ν M < 1 + b ν m * and 1 + s ν M < 1 + b ν m * , then the integral in (45) is convergent.
The set of values s C for which (45) exists and is finite is called the region of convergence (ROC).
Proof. 
By considering n 1 < 0 and n 2 > 0 , we split the summation corresponding to the nabla–Laplace transform into three summations; using the assumptions on the function on f, we have
n = + ν n f ( t n ) e ( t , t 0 , s ) V n = n 1 + 1 n 2 f ( t n ) e ( t , t 0 , s ) + A V n = n 1 e ( t , t 0 , s ) / e ( t , t 0 , a ) + B V n = n 2 + 1 + e ( t , t 0 , s ) / e ( t , t 0 , b )
where V = sup n N ν n .
It is obvious that the first summation on the right-hand side of the inequality is finite. The finiteness of the second summation can be obtained as follows
n = n 1 e ( t , t 0 , s ) / e ( t , t 0 , b ) n = n 1 k = 0 n 1 1 + s ν k k = 0 n 1 1 + b ν k n = n 1 k = 0 n 1 1 + s ν M k = 0 n 1 1 + b ν m * n = n 1 1 + s ν M 1 + b ν m * n <
when 1 + s ν M < 1 + b ν m * , where ν m * and ν M are as supposed in assumption (i).
For the third summation, we have
n = n 2 + 1 + e ( t , t 0 , s ) / e ( t , t 0 , a ) n = n 2 + 1 + k = 1 n 1 + a ν k k = 1 n 1 + s ν k n = n 2 + 1 + k = 1 n 1 + a ν M * k = 1 n 1 + s ν m n = n 2 + 1 + 1 + a ν M * 1 + s ν m n <
when 1 + a ν M * < 1 + s ν m , where ν M and ν m are as supposed in assumption (ii).
Consequently, (50) and (51) give the convergence of the series defined in (45).
It is clear that the verifications of the two inequalities in this theorem are sufficient conditions for the existence of the nabla bilateral LT.
Example 2.
As an application of the LT just introduced, we computed the transform of the nabla derivative obtained in Example 1, using N = 500 exponentials corresponding to equal values for s. Such values were generated according to the following formula: s ( k ) = i ω k = i π N ( k 1 ) , k = 1 , 2 , N . In Figure 3, we depict the absolute value of the transform as a function of ω k .
It is interesting to note the ringing at low frequencies due to the truncation of the signal and the similarity to the classic continuous-time case.

3.3. Backward Compatibility

We are going to show that one-sided -Laplace transforms given by (46) are backward compatible. Assume first that t t 0 , to obtain the right transform that enjoys the following properties:
  • Let T = R 0 + and t 0 = 0 . Then, we have ρ ( t 0 ) = 0 and e ( t , 0 ; s ) = e s t . From (46), we obtain the modified classical one-sided LT [19]
    L + f ( t ) ( s ) = 0 + e s t f ( t ) d t .
  • Let T = Z 0 + , i.e t = n and t 0 = 0 . Then, we have e ( n , 0 ; s ) = ( 1 + s ) n for n Z 0 + . By changing them with corresponding terms in (45) and using the relation between -integration and -summation (11), we have
    L + f ( t ) ( s ) = ρ ( 0 ) + f ( t ) e ( t , 0 ; s ) t = n = 0 f ( n ) z n
    where the transformation 1 + s = z was used. The last summation in (53) is the well-known Z transform [3].
Assume now that t < t 0 . Then, we have u ( ± ( t t 0 ) ) = u ( t 0 t ) , which leads to the left transform that verifies:
  • Let T = R 0 and t 0 = 0 . Then, we have ρ ( t 0 ) = 0 and e ( t , 0 ; s ) = e s t . By replacing these with corresponding terms in (48), we obtain the classical one-sided left LT
    L f ( t ) ( s ) = 0 e s t f ( t ) d t .
  • Let T = Z , i.e., t = n and t 0 = 0 . Then, we have ρ ( 0 ) = 1 and e ( n , 0 ; s ) = ( 1 + s ) n . By applying these in (48), we have
    L f ( t ) ( s ) = ρ ( 0 ) f ( t ) e ( t , t 0 ; s ) t = n = 1 f ( n ) z n
    where the transformation 1 + s = z is used. The last summation is the left Z transform.
Remark 1.
The left transforms are not very useful in applications. For this reason, we will stop considering it, in the following.

3.4. Some Properties of the Unilateral -Laplace Transform

The unilateral -LT enjoys some interesting properties:
  • Linearity
    The linearity of -LT is an obvious result of its definition (45).
  • Transform of the -derivative.
    We suppose that lim t + f ( t ) e ( t , t 0 ; s ) = 0 holds for a given function f : T R . By integrating the parts formula for the nabla derivative in (15) to the integral, we have:
    L + f ( t ) ( s ) = ρ ( t 0 ) + f ( t ) e ( t , t 0 ; s ) t = f ( ρ ( t 0 ) ) e ( ρ ( t 0 ) , t 0 ; s ) + s ρ ( t 0 ) + f ( ρ ( t ) ) e ( t , t 0 ; s ) t = f ( ρ ( t 0 ) ) e ( ρ ( t 0 ) , t 0 ; s ) + s L + f ( ρ ( t ) ) ( s ) .
    If T = R 0 + , then ρ ( t ) = t , and we have from (54)
    L + f ( t ) ( s ) = f ( t 0 ) + s L + f ( t ) ( s ) .
    which is the well-known result of the classical LT. It is important to mention that this is a property of the LT, not of any other operator or system.
  • The transform of the -anti-derivative.
    Let F ( t ) = t 0 t f ( τ ) τ with lim t + F ( t ) e ( t , t 0 ; s ) = 0 . By the nabla derivative of -exponential e ( t , t 0 ; s ) = s e ( t , t 0 ; s ) and the nabla derivative of the multiplication of two functions
    F ( t ) G ( t ) = F ( t ) G ( t ) + F ( ρ ( t ) ) G ( t ) ,
    we have
    L + F ( ρ ( t ) ) ( s ) = ρ ( t 0 ) + F ( ρ ( t ) ) e ( t , t 0 ; s ) t = 1 s lim t + F ( t ) e ( t , t 0 ; s ) F ( ρ ( t 0 ) ) e ( ρ ( t 0 ) , t 0 ; s ) ) + 1 s ρ ( t 0 ) + f ( t ) e ( t , t 0 ; s ) t = 1 s F ( ρ ( t 0 ) ) e ( ρ ( t 0 ) , t 0 ; s ) ) + 1 s L f ( t ) ( s ) .
    If t R , then this equality turns into the well-known relation
    L + F ( t ) ( s ) = 1 s L + f ( t ) ( s ) .
  • Initial value.
    Suppose that f ( t ) has -LT. Using the property “Transform of the -derivative”, we have
    lim s s L + f ( ρ ( t ) ) ( s ) f ( ρ ( t 0 ) ) e ( ρ ( t 0 ) , t 0 ; s ) = lim s + L + f ( t ) ( s ) = 0
    since lim s L + f ( t ) ( s ) = 0 is satisfied for a function f whose nabla derivative exists and of exponential type II. Hence,
    lim s + s L + f ( ρ ( t ) ) ( s ) = f ( ρ ( t 0 ) ) e ( ρ ( t 0 ) , t 0 ; s )
    holds.
  • Final value.
    Let us suppose that f ( t ) has -LT. From the property “Transform of the -derivative”, one can write
    lim s 0 s L + f ( ρ ( t ) ) ( s ) = lim s 0 L + f ( t ) ( s ) + f ( ρ ( t 0 ) ) e ( ρ ( t 0 ) , t 0 ; s ) = lim s 0 ρ ( t 0 ) + f ( t ) e ( t , ρ ( t 0 ) ; s ) t + f ( ρ ( t 0 ) ) = ρ ( t 0 ) + f ( t ) t + f ( ρ ( t 0 ) ) = lim t + f ( t ) .
  • Time scaling.
    Let ρ ( t 0 ) = 0 and a > 0 . By changing the variable a t = τ and applying the property of the scale changing in [6] for the nabla exponential in the definition of -LT, we have
    L + f ( a t ) ( s ) = 0 + f ( a t ) e ( t , t 0 ; s ) t = 1 a 0 + f ( τ ) e ( τ a , t 0 ; s ) τ = 1 a 0 + f ( τ ) e ( τ , t 0 ; s a ) τ = 1 a L + f ( t ) ( s a ) .

3.5. Inverse and Bilateral Δ -Laplace Transform

As such, in -case, for revealing the definition of the Δ -Laplace transform in the correct way, the inverse Δ -Laplace transform will be first introduced. As the inverse Laplace transform of F Δ ( s ) = 1 should equal the delta impulse given by (20), the formula for the inverse Laplace transform is obtained from the following formula
f t n = 1 2 π i C F Δ ( s ) e t , t 0 ; s d s
where the integration path, C , is a simple closed contour in a region of analyticity of the integrand and the poles of the (nabla) exponential.
The integrand of (56) is analytic for t > t 0 and has n poles at s = + 1 μ k ( t 0 ) ( k = 1 , 2 , , n ) for t < t 0 when we assume that F Δ ( s ) = 1 . For both cases, the integral is null. To obtain the expected result, a pole for n = 0 (or t = t 0 ) has to be introduced. For this, it is obvious that the exponential existing in the integrand have to be translated to the future by μ ( t ) . This leads the formulation in (56) to convert into the following
f ( t ) = 1 2 π i C F Δ ( s ) e t + μ ( t ) , t 0 ; s d s ,
which is the definition of the inverse Δ -Laplace transform. In addition to this, for n = 0 ,
f ( t ) = 1 2 π i C F Δ ( s ) e t + μ ( t ) , t 0 ; s d s = 1 2 π i C 1 1 s μ 1 ( t 0 ) = 1 μ 1 ( t 0 ) = 1 μ 0 .
So, the inverse Δ -Laplace transform of F Δ ( s ) = 1 is equal to δ Δ t , t 0 given by (20) as expected. Furthermore, by the Δ -differentiating on both sides of (57), we have L Δ f Δ ( t ) ( s ) = s L Δ f ( σ ( t ) ) ( s ) . This result will be reformulated by the one-sided delta direct transform in Section 3.7.
Definition 5.
The delta Laplace transform of a function f : T C is defined by
F Δ ( s ) = L Δ { f ( t ) } ( s ) = n = μ n f ( t n ) [ e t n , t 0 ; s ] 1
for those values s C for which the corresponding series converges.
The relationship between Δ -integral and Δ -summation given in Section 2.3 enable us to give the formula in (58) in the integral form as follows
L Δ f ( t ) ( s ) = + f ( t ) [ e t , t 0 ; s ] 1 Δ t .
Theorem 5.
(Inverse Δ -Laplace transform). Let F Δ ( s ) be the Laplace transform of f : T C defined by (58). Then the formula in (57) represents the inverse Δ -Laplace transform.
Proof. 
Inserting the delta–Laplace transform expression of f ( t ) given by (58) into (57), we have
L Δ 1 { F Δ ( s ) } ( t n ) = 1 2 π i C k = μ k f ( t k ) 1 e t k , t 0 ; s e t n + 1 , t 0 ; s d s = k = μ k f ( t k ) 1 2 π i C 1 e t k , t 0 ; s e t n + 1 , t 0 ; s d s = k = μ k f ( t k ) 1 2 π i C e t n + 1 , t k ; + s d s = k = μ k f ( t k ) δ Δ ( t k , t n ) = f ( t n )
for all n N , where we use the uniform convergence of the summation representing the delta Laplace transform, the relationship between nabla and delta exponentials given by [6] as follows:
1 e t , t 0 ; s e τ , t 0 ; s = e ( τ , t ; s )
and
1 2 π i C e t n + 1 , t k ; + s d s = 0 if k n 1 μ n if k = n .

3.6. The Unilateral Δ -Laplace Transform

Definition 6.
We define the unilateral Δ -LT
L Δ f ( t ) ( s ) = n = u Δ ( ( t n t 0 ) ) μ n f ( t n ) [ e t n , t 0 ; s ] 1
where s C , t 0 is the reference point used for defining the delta exponential and u Δ is the delta unit step function given by (18). The sign of the argument of u in formulation (45) defines the left (−) and right (+) transforms.
From the definition of the Δ -unit step function, the right unilateral Laplace transform has the summation form of
L Δ f ( t ) ( s ) = n = 1 μ n f ( t n ) [ e t n , t 0 ; s ] 1
which does not involve μ 0 f ( t 0 ) . By considering the relationship between the Δ -integration and Δ -summation in Section 2.3, the integral form of the unilateral right Δ -Laplace transform is given by
L Δ + f ( t ) ( s ) = σ ( t 0 ) + f ( t ) [ e t , t 0 ; s ] 1 Δ t .
If t R , then σ ( t 0 ) t 0 + and the above formula turns into the following
L + f ( t ) ( s ) = t 0 + + f ( t ) [ e t , t 0 ; s ] 1 Δ t .
It means that t 0 is not fully involved.
It is obvious that the unilateral left Δ -Laplace transform has the following integral form
L Δ f ( t ) ( s ) = n = 0 μ n f ( t n ) [ e t n , t 0 ; s ] 1 = σ ( t 0 ) f ( t ) [ e t , t 0 ; s ] 1 Δ t .
Now, let us discuss the existence and the region of convergence D Δ ( f ) . The proof of the following theorem will be omitted here since it can be proved by following Theorem 4.
Theorem 6.
Let f ( t ) be a bounded function of the exponential type, i.e.,
(i) There exists A R + and a R , such that
f ( t n ) < A e ( t n , t 0 ; a ) ( n < n 1 < 0 )
when 1 a μ M * < 1 s μ m , where μ M * is one of μ k , so that 1 a μ M * = max μ k 1 a μ k , μ m = min k 1 , 2 , , n μ k with n < 0 , and
(ii) There exists B R + and b R , such that
f ( t n ) < B e ( t n , t 0 ; b ) ( n > n 2 > 0 )
when 1 s μ M < 1 b μ m * , where μ m * is one of μ k , so that 1 b μ m * = min μ k 1 b μ k , μ M = max k { 0 , 1 , 2 , , n 1 } μ k with n > 0 .
If s C is inside the intersection of the circles defined by 1 a μ M * < 1 s μ m and 1 s μ M < 1 b μ m * , then the integral in (62) is convergent.
The set of values s C for which (62) exists and is finite is called the region of convergence (ROC).
Let us now prove that the one-sided Δ -LT defined by (62) ensures backward compatibility. At first, we suppose that t t 0 . Then, from (63) and (64), we determine the following relations:
  • Let T = R 0 + and t 0 = 0 . Then, we have σ ( t 0 ) = 0 and e ( t , 0 ; s ) = e s t . By interchanging these with the corresponding terms in (64), we obtain the classical one-sided Laplace transformation
    L Δ + f ( t ) ( s ) = 0 + + e s t f ( t ) d t .
  • Let T = Z 0 + and t 0 = 0 . Then, we have e ( t , t 0 ; s ) = ( 1 s ) n . By considering these in (63) and (64), we have
    L Δ + f ( t ) ( s ) = σ ( 0 ) + f ( t ) [ e t , t 0 ; s ] 1 Δ t = n = 1 f ( n ) ( 1 s ) n = n = 1 f ( n ) z n ,
    where the transformation 1 s = z 1 is used. The summation in (67) is the well-known right Z transform without the first term f ( 0 ) . From this aspect, this result differs from the nabla case.
At the moment, we assume t < t 0 . Then, we obtain the following relations:
  • Let T = R and t 0 = 0 . Then, we have σ ( t 0 ) = 0 and e ( t , 0 ; s ) = e s t . By applying these to (66), it yields the classical one-sided left LT
    L Δ f ( t ) ( s ) = 0 e s t f ( t ) d t .
  • Let T = Z and t 0 = 0 . Then, we have e ( t , t 0 ; s ) = ( 1 s ) n . By changing these with corresponding terms in (62) and considering the relation between Δ -integration and Δ -summation in (14) we have
    L Δ f ( t ) ( s ) = σ ( 0 ) f ( t ) [ e t , t 0 ; s ] 1 Δ t = n = ρ ( σ ( 0 ) ) f ( n ) ( 1 s ) n = n = 0 f ( n ) z n ,
    where the transformation 1 s = z 1 was used. The last summation consists of the sum of the well-known left Z transform and the term f ( 0 ) .

3.7. Some Properties of the Unilateral Δ -Laplace Transformation

As we have done before, here, we only take into account the properties of L Δ + . We will give some properties of the operator L Δ + in the following:
  • Linearity.
    The linearity of the Δ -Laplace transformation can be easily obtained from its integral representation.
  • Transform of the Δ -derivative.
    Let lim t + f ( t ) [ e t , t 0 ; s ] 1 = 0 for a given function f : T R . By integrating the parts formula for the delta derivative in (16) to (62), we have
    L Δ + f Δ ( t ) ( s ) = σ ( t 0 ) + f Δ ( t ) [ e t , t 0 ; s ] 1 Δ t = lim t + f ( t ) [ e t , t 0 ; s ] 1 f ( σ ( t 0 ) ) [ e σ ( t 0 ) , t 0 ; s ] 1 σ ( t 0 ) + f ( σ ( t ) ) [ e 1 t , t 0 ; s ] Δ Δ t = f ( σ ( t 0 ) ) [ e σ ( t 0 ) , t 0 ; s ] 1 + s σ ( t 0 ) + f ( σ ( t ) ) [ e t , t 0 ; s ] 1 Δ t
    which yields
    L Δ + f Δ ( t ) ( s ) = f ( σ ( t 0 ) ) [ e σ ( t 0 ) , t 0 ; s ] 1 + s L Δ + f ( σ ( t ) ) ( s ) .
    It is remarkable that if T = R 0 + , then we have the well-known result of the classical LT
    L + f ( t ) ( s ) = f ( t 0 + ) + s L + f ( t ) ( s ) .
  • The transform of the Δ -anti-derivative.
    Let F ( t ) = t 0 t f ( τ ) Δ τ with lim t + F ( t ) [ e t , t 0 ; s ] 1 = 0 . By the delta derivative of Δ -exponential [ e 1 t , t 0 ; s ] Δ = s [ e t , t 0 ; s ] 1 and the delta derivative of the multiplication of two functions F ( t ) G ( t ) Δ = F Δ ( t ) G ( t ) + F ( σ ) ) G Δ ( t ) , we have
    L Δ F ( σ ( t ) ) ( s ) = σ ( t 0 ) + F ( σ ( t ) ) [ e t , t 0 ; s ] 1 Δ t = 1 s lim t + F ( t ) [ e t , t 0 ; s ] 1 F ( σ ( t 0 ) ) [ e σ ( t 0 ) , t 0 ; s ] 1 + 1 s σ ( t 0 ) + f ( t ) [ e t , t 0 ; s ] 1 Δ t = 1 s F ( σ ( t 0 ) ) [ e σ ( t 0 ) , t 0 ; s ] 1 + 1 s L Δ f ( t ) ( s ) .
  • Initial value.
    Assume that f ( t ) has Δ -LT. Taking into account the property “Transform of the Δ -derivative” in (69), we have
    lim s + s L Δ + f ( σ ( t ) ) ( s ) + f ( σ ( t 0 ) ) [ e σ ( t 0 ) , t 0 ; s ] 1 = lim s + L Δ + f Δ ( t ) ( s ) = 0
    since lim s + L Δ + f Δ ( t ) ( s ) = 0 is satisfied for a function f whose delta derivative exists and is of exponential type II. Hence, we have
    lim s + s L Δ + f ( ρ ( t ) ) ( s ) = f ( σ ( t 0 ) ) [ e σ ( t 0 ) , t 0 ; s ] 1 .
  • Final value.
    Assume that f ( t ) has Δ -LT. From the property “Transform of the Δ -derivative” in (69) we have
    lim s 0 s L Δ + f ( σ ( t ) ) ( s ) = lim s 0 L + f Δ ( t ) ( s ) + f ( σ ( t 0 ) ) [ e σ ( t 0 ) , t 0 ; s ] 1 = lim s 0 σ ( t 0 ) + f Δ ( τ ) [ e τ , t 0 ; s ] 1 Δ τ + f ( σ ( t 0 ) ) = σ ( t 0 ) + f Δ ( τ ) Δ τ + f ( σ ( t 0 ) ) = lim t + f ( t ) .
  • Time scaling.
    Let t 0 = 0 and a > 0 . By changing the variable a t = τ and using the property of the scale changing for the nabla exponential in the definition of Δ -LT, we have
    L Δ + f ( a t ) ( s ) = 0 + f ( a t ) [ e t , t 0 ; s ] 1 Δ t = 1 a 0 + f ( τ ) e τ a , t 0 ; s 1 Δ τ = 1 a 0 + f ( τ ) e ( τ , t 0 ; s a ) 1 Δ τ = 1 a L Δ + f ( t ) ( s a ) .

4. Conclusions

We propose unilateral nabla and delta Laplace transforms on non-uniform time scales by following the methodology introduced in [6]. We showed that both coincide with the Z transform if the time scale is Z 0 + and with unilateral LT if the time scale is R 0 + . Moreover, we obtained properties enjoyed by classical one-sided LT and ZT, showing the fully compatible character.

Author Contributions

M.Ş. and M.D.O. contributed equally to this article and approved the final manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partially funded by national funds through the Foundation for Science and Technology (FCT) within the scope of the CTS Research Unit—Center of Technology and Systems/UNINOVA/FCT/NOVA, under reference UIDB/00066/2020.

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. 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]
  2. Aulbach, B.; Hilger, S. A unified approach to continuous and discrete dynamics. In Qualitative Theory of Differential Equations; Colloquia Mathematica Sociefatis János Bolyai; North-Holland Publishing: Amsterdam, The Netherlands, 1990; Volume 53, pp. 37–56. [Google Scholar]
  3. Proakis, J.G.; Manolakis, D.G. Digital Signal Processing: Principles, Algorithms, and Applications; Prentice Hall: Hoboken, NJ, USA, 2007. [Google Scholar]
  4. Ortigueira, M.D.; Coito, F.J.; Trujillo, J.J. Discrete-time differential systems. Signal Process. 2015, 107, 198–217. [Google Scholar] [CrossRef]
  5. Ortigueira, M.D. The comb signal and its Fourier transform. Signal Process. 2001, 81, 581–592. [Google Scholar] [CrossRef] [Green Version]
  6. Ortigueira, M.D.; Torres, D.F.; Trujillo, J.J. Exponentials and Laplace transforms on nonuniform time scales. Commun. Nonlinear Sci. Numer. Simul. 2016, 39, 252–270. [Google Scholar] [CrossRef] [Green Version]
  7. Bohner, M.; Peterson, A. Laplace transform and Z-transform: Unification and extension. Methods Appl. Anal. 2002, 9, 151–158. [Google Scholar] [CrossRef] [Green Version]
  8. Davis, J.M.; Gravagne, I.A.; Jackson, B.J.; Marks II, R.J.; Ramos, A.A. The Laplace transform on time scales revisited. J. Math. Anal. Appl. 2007, 332, 1291–1307. [Google Scholar] [CrossRef] [Green Version]
  9. Ahrendt, C. The Laplace transform on time scales. Pan. Am. Math. J. 2009, 19, 1–36. [Google Scholar]
  10. Bohner, M.; Guseinov, G.S. The Laplace transform on isolated time scales. Comput. Math. Appl. 2010, 60, 1536–1547. [Google Scholar] [CrossRef] [Green Version]
  11. Davis, J.M.; Gravagne, I.A.; Marks, R.J. Convergence of unilateral Laplace transforms on time scales. Circuits Syst. Signal Process. 2010, 29, 971–997. [Google Scholar] [CrossRef]
  12. Bastos, N.R.; Mozyrska, D.; Torres, D.F. Fractional derivatives and integrals on time scales via the inverse generalized Laplace transform. arXiv 2010, arXiv:1012.1555. [Google Scholar]
  13. Karpuz, B. On uniqueness of the Laplace transform on time scales. arXiv 2011, arXiv:1102.5590. [Google Scholar]
  14. Rahmat, M.R.S. The (q, h)-Laplace transform on discrete time scales. Comput. Math. Appl. 2011, 62, 272–281. [Google Scholar] [CrossRef] [Green Version]
  15. Davis, J.M.; Gravagne, I.A.; Marks, R.J. Bilateral Laplace transforms on time scales: Convergence, convolution, and the characterization of stationary stochastic time series. Circuits Syst. Signal Process. 2010, 29, 1141–1165. [Google Scholar] [CrossRef]
  16. Bohner, M.; Peterson, A. Dynamic Equations on Time Scales: An Introduction with Applications; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2001. [Google Scholar]
  17. Ortigueira, M.D.; Machado, J.T. Revisiting the 1D and 2D Laplace transforms. Mathematics 2020, 8, 1330. [Google Scholar] [CrossRef]
  18. Ortigueira, M.D.; Valério, D. Fractional Signals and Systems; De Gruyter: Berlin, Germany; Boston, MA, USA, 2020. [Google Scholar]
  19. Roberts, M. Signals and Systems: Analysis Using Transform Methods and MATLAB, 2nd ed.; McGraw-Hill: New York, NY, USA, 2003. [Google Scholar]
Figure 1. Scheme for a time scale.
Figure 1. Scheme for a time scale.
Mathematics 10 04552 g001
Figure 2. Nabla derivative of a sinusoidal signal.
Figure 2. Nabla derivative of a sinusoidal signal.
Mathematics 10 04552 g002
Figure 3. Transforms of the signals used in Example 1.
Figure 3. Transforms of the signals used in Example 1.
Mathematics 10 04552 g003
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Şan, M.; Ortigueira, M.D. Unilateral Laplace Transforms on Time Scales. Mathematics 2022, 10, 4552. https://doi.org/10.3390/math10234552

AMA Style

Şan M, Ortigueira MD. Unilateral Laplace Transforms on Time Scales. Mathematics. 2022; 10(23):4552. https://doi.org/10.3390/math10234552

Chicago/Turabian Style

Şan, Müfit, and Manuel D. Ortigueira. 2022. "Unilateral Laplace Transforms on Time Scales" Mathematics 10, no. 23: 4552. https://doi.org/10.3390/math10234552

APA Style

Şan, M., & Ortigueira, M. D. (2022). Unilateral Laplace Transforms on Time Scales. Mathematics, 10(23), 4552. https://doi.org/10.3390/math10234552

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