Discrete Fractional Calculus (2015)
3. Nabla Fractional Calculus
3.12. Convolution
We are now ready to investigate one of the most important properties in solving initial-value fractional nabla difference equations: convolution. This definition is motivated by the desire to express the fractional nabla sums and fractional nabla differences as convolutions of arbitrary functions and Taylor monomials. As a consequence, the resulting properties that stem from this definition are, in fact, consistent with the standard convolution. Many of the results in this section appear in Hein et al. [119] and Ahrendt et al. [3].
Definition 3.77.
For , we define the nabla convolution product of f and g by
Example 3.78.
Use the definition of the nabla convolution product to find 1 ∗Sin p (⋅ , a), p ≠ 0, ±i. By Definition 3.77,
for .
Example 3.79.
Use the definition of the (nabla) convolution product to find
Assume q ≠ 0. By Definition 3.77 and using the second integration by parts formula, we have that
Solving for , we obtain
for . We leave it to the reader to show that this last formula is also valid if q = 0.
Theorem 3.80.
Assume and . Then
for .
Proof.
The result follows from the following:
for . □
Theorem 3.81 (Nabla Convolution Theorem).
Assume and their nabla Laplace transforms converge for |s − 1| < r for some r > 0. Then
for |s − 1| < r.
Proof.
The following proves our result:
for | s − 1 | < r. □
With the above result and the uniqueness of the Laplace transform, it follows that the convolution product is commutative and associative (see Exercise 3.32).
We next establish properties of the Laplace transform that will be useful in solving initial value problems for integer nabla difference equations.
Theorem 3.82 (Transformation of Fractional Sums).
Assume ν > 0 and the nabla Laplace transform of converges for |s − 1| < r for some r > 0. Then
for .
Proof.
The result follows since
for . □
Assuming that ν is a positive integer, this result is consistent with the formula in the continuous case for the Laplace transform of the n-th iterated integral of a function. We want to establish similar properties for fractional differences; however, we will first establish integer-order difference properties.
Theorem 3.83 (Transform of Nabla Difference).
Assume is of exponential order r > 0. Then
for |s − 1| < r.
Proof.
Note that since we are assuming that , we have that and so we can consider . Since is of exponential order r > 0, it follows that is of exponential order r > 0. It follows that for | s − 1 | < r,
This completes the proof. □
The following is a simple example where we will use the Nabla Convolution Theorem and Theorem 3.83 to solve an initial value problem.
Example 3.84.
Use the Nabla Convolution Theorem to help you solve the IVP
If y(t) is the solution of this IVP and its Laplace transform, Y a (s), exists, then we have that
Using the initial condition and solving for Y a (s) we obtain
Using the Nabla Convolution Theorem we see that
Using Example 3.79 we find that
is the solution of our given IVP on . Of course, in this simple example one could also use partial fractions to find y(t).
We can then generalize this result for an arbitrary number of nabla differences.
Theorem 3.85 (Transform of n-th-Order Nabla Difference).
Assume is of exponential order r > 0. Then
(3.38)
for |s − 1| < r, for each .
Proof.
Note that since f is of exponential order r > 0, ∇ n f is of exponential order r > 0 for each . Hence converges for | s − 1 | < r for each . The proof of (3.38) is by induction for . The base case n = 1 follows from Theorem 3.83. Now assume n ≥ 1 and (3.38) holds for | s − 1 | < r. Then, using Theorem 3.83, we have that
Hence, (3.38) holds when n is replaced by n + 1 and the proof is complete. □
Example 3.86.
Solve the IVP
If y(t) is the solution of this IVP and we let we have that
Using the initial conditions we have
Solving for Y a (s) we have that
It follows that
for .