Laplace transform IV - Convolution redux

Recall that for the Fourier transform we had

F(fβˆ—g)=Fgβ‹…Fg,

where fβˆ—g was the convolution of f and g, defined by

(fβˆ—g)(t)=βˆ«βˆ’βˆžβˆžf(tβˆ’x)g(x)dx.

Something very similar holds for the Laplace transform, but we have to make a few small adjustments.

Definition of convolution


First note that, when dealing with the Laplace transform, we are now dealing with functions defined on [0,∞). In particular, in the integral above the function g(x) would not be defined for x<0. Similar, the function f(tβˆ’x) would not be defined for tβˆ’x<0, i.e., for x>t. So we need to restrict the bounds on the definite integral in order for both terms in the integrand to even be defined.

Fortunately, it turns out that is the only thing we need to do:

Definition of convolution for functions defined on [0,∞)

Fr functions f,g defined on [0,∞), their convolution is the function denoted fβˆ—g defined by the integral formula

(fβˆ—g)(t)=∫0tf(tβˆ’x)g(x)dx.

Example

Define functions f,g on [0,∞) by f(t)=t and g(t)=et. We then compute

(fβˆ—g)(t)=∫0tf(tβˆ’x)g(x)dx=∫0t(tβˆ’x)exdx=t∫0texdxβˆ’βˆ«0txexdx=t(etβˆ’e0)βˆ’[xexβˆ’ex]x=0x=t=t(etβˆ’1)βˆ’((tetβˆ’et)βˆ’(0β‹…e0βˆ’e0))=βˆ’1βˆ’t+et.

Convolution and the Laplace transform


Just as with the Fourier transform, it turns out that the Laplace transform converts convolution product into usual product:

L(fβˆ—g)=Lfβ‹…Lg.

Read "in reverse", this same property tells us

Lβˆ’1(Fβ‹…G)=Lβˆ’1Fβˆ—Lβˆ’1G,

and this is the situation in which we'll most often find ourselves running into the convolution.

Example

Suppose we wished/needed to compute the inverse Laplace transform of the function F(s)=1s2(sβˆ’1).

One option is to first use partial fractions to rewrite F(s) as

F(s)=βˆ’1sβˆ’1s2+1sβˆ’1.

The above calculation took a few minutes behind the scenes, but the payoff is that we can now instantly conclude

Lβˆ’1F=βˆ’Lβˆ’1(1s)βˆ’Lβˆ’1(1s2)+Lβˆ’1(1sβˆ’1)=βˆ’1βˆ’t+et.

A second (new) option is to view F(s) as a product of two functions with known inverse transforms, namely

F(s)=1s2β‹…1sβˆ’1

Then

Lβˆ’1F=Lβˆ’1(1s2β‹…1sβˆ’1)=Lβˆ’1(1s2)βˆ—Lβˆ’1(1sβˆ’1)=tβˆ—et=βˆ’1βˆ’t+et.

Although very different, these two methods (at least in this example) take roughly the same amount of time. In the first method, we needed to compute a partial fraction decomposition, while in the second we needed to compute a convolution integral.

Suggested next notes


Laplace transform V - Solving differential equations