Analogues of the Laplace Transform and Z -Transform with Piecewise Linear Kernels

: Two new transforms with piecewise linear kernels are introduced. These transforms are analogues of the classical Laplace transform and Z -transform. Properties of these transforms are investigated and applications to ordinary differential equations and integral equations are provided. This article is ideal for study as a foundational project in an undergraduate course in differential and/or integral equations.


Introduction
An integral transform maps a function from its original function space into another function space via integration. In many cases, properties of the original function might be more easily characterised in the transformed space rather than in the original space. Integral transforms arise in many areas of mathematics, e.g., differential and integral equations, probability, number theory and computer science (see for instance [1][2][3][4][5][6][7] and the comprehensive references therein). The reader is also referred to the recent article [8] that introduced and studied a broad class of integral transforms and which includes many well-known integral transforms as special cases. A discrete transform is one where the input function is a sequence and the integral is typically replaced by a summation [6,9].
Two of the most well-known and also widely used transforms are the Laplace transform [3,6,10] L{ f (t); s} = The Laplace transform and Z-transform have many nice properties, hence they are ubiquitous in many areas of applied mathematics and engineering [6,12,13]. Note that here for simplicity we assume that s and z are positive real numbers although they can also be considered in the complex plane.
Suppose that we replace the above two kernels by similarly behaved but piecewise linear functions. More precisely, let us define two new transforms and provided the improper integral and series converge. More specific assumptions on the appropriate function spaces will be given later. Then it is easy to see that K s and K z are nonnegative, nonincreasing and satisfy (1) and (2), respectively. Moreover, R c and R d are linear operators. The goal of this article is to investigate some properties and applications of these analogues of the Laplace transform and Z-transform with piecewise linear kernels. We consider the continuous transform R c in Section 2, while Section 3 is devoted to the discrete transform R d . A brief discussion is given in Section 4. The results in this article are useful as a springboard for further investigation of other continuous and discrete transforms, and are particularly useful as a foundational project for a first course in differential and/or integral equations.

Properties of the Continuous Transform
By breaking up the interval of integration, we see that (3) can be simplified to ; s} exists for all s ∈ R + . In fact, the mapping s → R c { f (t); s} belongs to C(R + ). Hence we may write R c : U → V, where U = C(R + ) and V = C(R + ).

Remark 1.
Note, however, that the assumption U = C(R + ) is a sufficient, but not a necessary, condition for the R c -transform to exist. For example, if f ∈ C(R + ) but f (t) becomes unbounded as t tends to zero from the right (e.g., f (t) = 1/ √ t), then the R c -transform may still exist; see Example 6, where it is also pointed out that a similar situation arises for the Laplace transform. As will be seen in following propositions, depending on the transform property that one wishes to prove, stricter conditions on U may need to be imposed (e.g., U = C n (R + ) ⊂ C(R + ), where n ∈ N). The next result shows how R c transforms derivatives. We omit the proof since the result follows from straightforward integration by parts. Let n ∈ N be the order of the derivative. For comparison recall that L{ f (n) (t); s} = s n L{ f (t); s} − n ∑ j=1 s n−j f (j−1) (0), n ≥ 1.
Now, we look at how R c transforms integrals. Comparing with the Laplace transform, it is known that Proof. For notational convenience define hence g is uniformly continuous on [0, t] and differentiable on (0, t). From integration by parts we obtain Observing (5), we can rewrite the above equation as which proves the assertion.
The next proposition shows that R c has a scaling property. An analogous scaling property for the Laplace transform is Proposition 3. Let a > 0 and f ∈ C(R + ). Then Proof. The result follows from making the substitution u = at and using (5).
Let F(s) = L{ f (t); s}. Then the inverse Laplace transform is given by where c is a real number such that the contour path of integration is in the region of convergence of F(s). We now derive the inverse R c -transform. Unlike the inverse Laplace transform, the inverse R c -transform does not involve contour integration but differentiation of the transformed function.
; s} is twice continuously differentiable with respect to s, then We have from (5) and the Leibniz Integral Rule that The Fundamental Theorem of Calculus implies that Hence we see that The Laplace convolution of two functions f : R + → R and g : R + → R is defined as It can be shown that the convolution operator is commutative and the convolution property gives where F(s) = L{ f (t); s} and G(s) = L{g(t); s}. We wish to define a convolution operator associated with the R c -transform such that the convolution property holds, namely Thus we take (7) as the definition of the R c -convolution so that the convolution property necessarily holds. The R c -convolution f * g defined in (7) is expressed in terms of the R c -transforms of f and g. However, if we look at (6), we note that the Laplace convolution is not given in terms of the Laplace transforms of f and g. Hence in the next result we derive an alternative formula for the R c -convolution.
Proposition 5. The R c -convolution (7) of f and g is formally equivalent to Then we can compute Thus (7) gives Note that the R c -convolution operator is also commutative.

Examples
As a special case, when p = 0, we get R c {1; s} = s/2. Note that if −1 < p < 0, then the function t → t p does not belong to C(R + ) but R c {t p ; s} may still exist as an improper Riemann integral. This can be seen from and therefore , p > −1.
Of course, to be able to apply the results of the previous section for this case, we have to assume that p ≥ 0. It is known that Observe that L{t p ; s} will still exist as an improper Riemann integral when −1 < p < 0, just as for R c {t p ; s}.
By comparison, Replacing a by ia in (9) gives Thus, equating real and imaginary parts, we obtain Recalling that cosh(x) = cos(ix) and sinh(x) = −i sin(ix) for x ∈ R, we deduce that The above formulas can also be compared with their Laplace transform counterparts.   (8) and (9), we see that

Example 5.
Integral transforms are usually used to solve linear ordinary differential equations (ODEs). This example shows that the R c -transform can be used to solve some nonlinear ODEs. Let f : R + → R be a given continuous function of t, and consider where y = y(t) is to be determined. With a slight abuse of notation, if f : R + → R and g : R + → R are any continuously differentiable functions of t, then from integration by parts. Choosing f (t) = y(t) and g(t) = y (t), we see from (12) that Taking the R c -transform of (11), we obtain −y(0)y (0) where F(s) = R c { f (t); s}. Solving for y(s) 2 yields y(s) 2 = 2sF(s) + 2sy(0)y (0) + y(0) 2 .
Replacing s by t, Thus the solution of the nonlinear ODE (11) using the R c -transform is where the initial conditions y(0) and y (0) are assumed to be given. Alternatively, (11) can be solved by rewriting it as (yy ) = f (t) and possibly making a change of variables to reduce the order.

Example 6.
Let us solve the integral equation for y = y(t). Suppose that we take f (t) = −1 and g(t) = y(t) in Proposition 5. We obtain Therefore the integral Equation (13) can be expressed as Taking the R c -transform, applying the R c -convolution property and using (8), we have Hence (7) gives as the solution of the integral Equation (13). For comparison, we observe that (13) can be converted to an ODE by differentiating with respect to t. Thus This is a nonautonomous linear first-order ODE, which can be solved using the Laplace transform only if we rewrite it as 2y (t) + ty (t) = −3y(t). Even then, the term L{ty (t); s} = −Y (s) produces another ODE in Laplace transform space, which is of the same degree of difficulty as the original problem. However, using the method of integrating factors, the general solution can be expressed as where c is an arbitrary constant. Note also from (13) that y(0) = 1, thus c = 8 and we recover the solution (14) obtained using the R c -transform.

Example 7.
More generally, let us solve the linear integral equation for y = y(t), where f :R + → R and g :R + → R are arbitrary but given continuous functions.
Assume further that f (t) ≤ 0 for t ≥ 0. It follows that Note that in general the integral Equation (15) cannot be converted to a linear ODE by differentiation with respect to t, unlike in the previous example. Expanding the R c -convolution operator in (15), we have An application of the R c -convolution property to (15) gives where Y(s) = R c {y(t); s}, F(s) = R c { f (t); s} and G(s) = R c {g(t); s}. If we define Therefore from (7) we deduce that Finally, y can be expressed as a R c -convolution, i.e., The function g is given, while F (and therefore h) is obtained from the given function f by taking its R c -transform and replacing s by t.

Example 8.
Consider the nonlinear integral equation from Proposition 5 and g :R + → R is a given continuous function such that g(t) ≤ 0 for t ≥ 0. This implies that As before, we want to determine y = y(t).
Taking the R c -transform of (16) and invoking the convolution property, we get where Y(s) = R c {y(t); s} and G(s) = R c {g(t); s}. Proposition 4 implies that An auxiliary condition is needed to determine the correct sign above.

Properties of the Discrete Transform
Throughout this section we assume that z > 0. If z denotes the greatest integer less than or equal to z, then z ≤ z < z + 1. For all j ≥ z + 1 we see that j ≥ z + 1 > z and 1 − j/z < 0; hence (4) reduces to Let U be the collection of all real-valued sequences defined on N 0 . As (17) is a finite sum, we see that R d { f (j); z} ∈ R is always defined for any z > 0. Because of the presence of z , we see that the R d -transform is not necessarily continuous in z. Let V be the collection of all real-valued functions defined on R + . Then R d : U → V.

Remark 2.
We shall see later that the image of R d is in fact a proper subset of V.
The Z-transform satisfies the backward shift and forward shift properties respectively, where k ∈ N. Similar to the Z-transform, the R d -transform also has backward shift and forward shift properties although they are not as simple since the kernel of the Z-transform satisfies a semigroup property while the kernel of the R d -transform does not. In the following computations we will assume that f (j) = 0 if j < 0. Now suppose that 1 ≤ z < 2. Then (20) gives We deduce that a(2) = f (0) The pattern is now apparent, i.e., provided b(0) = 0. Furthermore, if F(z) is a R d -transform, not only must (20) hold but a : N 0 → R and b : N 0 → R are such that (21) is true.
We are now ready to state and prove the following result.

Proposition 8.
Let be the R d -transform of some sequence f : N 0 → R, where a : N 0 → R is arbitrary and b : Then f (j) = R −1 d {F(z); j} is given by Proof. For all j ≥ 1 we have We need to show that Recalling (17) and (23), we get

Discussion
In this article, we introduced two new transforms R c and R d whose kernels are piecewise linear analogues of the Laplace transform and Z-transform, respectively. We gave several examples and derived some properties of these two transforms.
For the continuous case, we showed that the inverse R c -transform has a form that is relatively straightforward to calculate. We defined a R c -convolution operator with the aid of the R c -inverse transform, and then derived an alternative formula. We showed through examples that the R c -transform can be used to solve certain (linear and nonlinear) ODEs and integral equations. As to be expected, since the transform properties for the derivatives are different, the R c -transform is not to be used to solve linear ODEs with constant coefficients since the Laplace transform is more efficient for this.
For the discrete case, we derived formulas for the backward shift and forward shift properties for the R d -transform. Since the kernel of the Z-transform has the semigroup property that K z (j + k) = K z (j)K z (k) for j, k ∈ N 0 , which does not hold for the R d -transform kernel, the formulas for the R d -transform tend to be more complicated than those for the Z-transform. Moreover, a R d -transform is necessarily of the form a( z ) − b( z )/z for some sequences a : N 0 → R and b : N 0 → R related through (23). This imposes a restriction on functions whose inverse R d -transforms exist.
Finding further interesting applications for the new transforms introduced in this article is still an open problem. One possible direction is in model fitting. For example, suppose that f (j) represents the population of some species at time j and is assumed to be modelled by a discrete logistic equation The algebraic system (29) is linear in r and r/K. Hence explicit analytical formulas for r and K (e.g., using Cramer's Rule) can be derived in terms of the measured population sizes given by f (0), f (1), . . . , f (n). As (27) is nonlinear and its analytical solution is unknown, parameter estimation techniques such as those based on least squares, for example, to estimate r and K are not straightforward to implement. The parameter estimation technique outlined above can be viewed as a discrete version of the integration-based techniques introduced in [14,15].