Loading [MathJax]/jax/output/HTML-CSS/jax.js
Skip to main content
Library homepage
 

Text Color

Text Size

 

Margin Size

 

Font Type

Enable Dyslexic Font
Mathematics LibreTexts

9.10: The Inverse Laplace Transform

( \newcommand{\kernel}{\mathrm{null}\,}\)

Until this point we have seen that the inverse Laplace transform can be found by making use of Laplace transform tables and properties of Laplace transforms. This is typically the way Laplace transforms are taught and used in a differential equations course. One can do the same for Fourier transforms. However, in the case of Fourier transforms we introduced an inverse transform in the form of an integral. Does such an inverse integral transform exist for the Laplace transform? Yes, it does! In this section we will derive the inverse Laplace transform integral and show how it is used.

We begin by considering a causal function f(t) which vanishes for t<0 and define the function g(t)=f(t)ect with c>0. For g(t) absolutely integrable,

|g(t)|dt=0|f(t)|ectdt<,

we can write the Fourier transform,

ˆg(ω)=g(t)eiωtdt=0f(t)eiωtctdt

and the inverse Fourier transform,

g(t)=f(t)ect=12πˆg(ω)eiωtdω.

Multiplying by ect and inserting ˆg(ω) into the integral for g(t), we find

f(t)=12π0f(τ)e(iωc)τdτe(iωc)tdω.

Letting s=ciω( so dω=ids), we have

f(t)=i2πcic+i0f(τ)esτdτestds

Note that the inside integral is simply F(s). So, we have

f(t)=12πic+iciF(s)estds.

The integral in the last equation is the inverse Laplace transform, called the Bromwich integral and is named after Thomas John I’Anson Bromwich (1875-1929). This inverse transform is not usually covered in differential equations courses because the integration takes place in the complex plane. This integral is evaluated along a path in the complex plane called the Bromwich contour. The typical way to compute this integral is to first chose c so that all poles are to the left of the contour. This guarantees that f(t) is of exponential type. The contour is closed a semicircle enclosing all of the poles. One then relies on a generalization of Jordan’s lemma to the second and third quadrants.1

Note

Closing the contour to the left of the contour can be reasoned in a manner similar to what we saw in Jordan’s Lemma. Write the exponential as est= e(sR+isI)t=esRteislt. The second factor is an oscillating factor and the growth in the exponential can only come from the first factor. In order for the exponential to decay as the radius of the semicircle grows, sRt<0. Since t>0, we need s<0 which is done by closing the contour to the left. If t<0, then the contour to the right would enclose no singularities and preserve the causality of f(t).

Example 9.10.1

Find the inverse Laplace transform of F(s)=1s(s+1).

Solution

The integral we have to compute is

f(t)=12πic+iciests(s+1)ds

This integral has poles at s=0 and s=1. The contour we will use is shown in Figure 9.10.1. We enclose the contour with a semicircle to the left of the path in the complex s-plane. One has to verify that the integral over the semicircle vanishes as the radius goes to infinity. Assuming that we have done this, then the result is simply obtained as 2πi times the sum of the residues. The residues in this case are:

Res[eztz(z+1);z=0]=limz0ezt(z+1)=1

and

Res[eztz(z+1);z=1]=limz1eztz=et.

Therefore, we have

f(t)=2πi[12πi(1)+12πi(et)]=1et.

clipboard_e5311c7e9377a0ed9390d9048bc278277.png
Figure 9.10.1: The contour used for applying the Bromwich integral to the Laplace transform F(s)=1s(s+1).

We can verify this result using the Convolution Theorem or using a partial fraction decomposition. The latter method is simplest. We note that

1s(s+1)=1s1s+1

The first term leads to an inverse transform of 1 and the second term gives et. So,

L1[1s1s+1]=1et.

Thus, we have verified the result from doing contour integration.

Example 9.10.2

Find the inverse Laplace transform of F(s)=1s(1+es).

Solution

In this case, we need to compute

f(t)=12πic+iciests(1+es)ds.

This integral has poles at complex values of s such that 1+es=0, or es=1. Letting s=x+iy, we see that

es=ex+iy=ex(cosy+isiny)=1.

We see x=0 and y satisfies cosy=1 and siny=0. Therefore, y=nπ for n an odd integer. Therefore, the integrand has an infinite number of simple poles at s=nπi,n=±1,±3,. It also has a simple pole at s=0.

In Figure 9.10.2 we indicate the poles. We need to compute the resides at each pole. At s=nπi we have

Res[ests(1+es);s=nπi]=limsnπi(snπi)ests(1+es)=limsnπiestses=enπitnπi,n odd.

At s=0, the residue is

Res[ests(1+es);s=0]=lims0est1+es=12.

clipboard_eca5764decfb9834d7a8cebfe55566dd5.png
Figure 9.10.2: The contour used for applying the Bromwich integral to the Laplace transform F(s)=11+es.

Summing the residues and noting the exponentials for ±n can be combined to form sine functions, we arrive at the inverse transform.

f(t)=12n odd enπitnπi=122k=1sin(2k1)πt(2k1)π.

The series in this example might look familiar. It is a Fourier sine series with odd harmonics whose amplitudes decay like 1/n. It is a vertically shifted square wave. In fact, we had computed the Laplace transform of a general square wave in Example 9.8.7.

In that example we found

L[n=0[H(t2na)H(t(2n+1)a)]]=1eass(1e2as)=1s(1+eas).

In this example, one can show that

f(t)=n=0[H(t2n+1)H(t2n)].

The reader should verify that this result is indeed the square wave shown in Figure 9.10.3.

clipboard_e69366f97a8a65dc3c0f93ea7f6e7e47d.png
Figure 9.10.3: Plot of the square wave result as the inverse Laplace transform of F(s)=1s(1+es with 50 terms.

This page titled 9.10: The Inverse Laplace Transform is shared under a CC BY-NC-SA 3.0 license and was authored, remixed, and/or curated by Russell Herman via source content that was edited to the style and standards of the LibreTexts platform.

Support Center

How can we help?