ODE-Project The Laplace Transform

Section 6.1 The Laplace Transform

Objectives

To understand and be able to compute the Laplace transform of a function $f (t),$ which can be accomplished by computing

$L (f) (s) = F (s) = \int_{0}^{\infty} e^{- s t} f (t) d t,$

provided the integral converges.
To understand that if we know the Laplace transform $F (s)$ of a function, we can recover the original function using the inverse Laplace transform, $L^{- 1} (F (s)) (t) .$
To understand that both the Laplace transform and inverse Laplace transform are linear operators.

Consider the electrical circuit governed by the differential equation

L I^{″} + R I^{'} + \frac{1}{C} I = E^{'} (t) .

The voltage function,

E^{'} (t),

might have discontinuities. For example, the voltage in the circuit can be periodically turned on and off. The previous methods that we have used to solve second order linear differential equations may not apply here. However, the Laplace transform, an integral transform, gives a method of solving such equations.

🔗

As a second example, let us consider a population of fish that is governed by exponential growth,

🔗

\begin{aligned} \frac{d P}{d t} & = k P \\ P (0) & = P_{0}, \end{aligned}

and suppose that we wish to determine the effects of seasonal fishing. In other words, harvesting will not be continuous. For example, we might only allow fishing at a constant rate

r

during the first half of the year,

🔗

H (t) = {\begin{cases} r & 0 \leq t \leq 1 / 2, \\ 0 & 1 / 2 < t < 1, \\ r & 1 \leq t \leq 3 / 2, \\ 0 & 3 / 2 < t < 2, \\ ⋮ \end{cases}

Our initial value problem now becomes

🔗

\begin{aligned} \frac{d P}{d t} & = k P - H \\ P (0) & = P_{0}, \end{aligned}

It should be clear that we need some additional tools to analyze differential equations possessing discontinuous terms.

🔗

Given an initial value problem

🔗

\begin{aligned} a y^{″} + b y^{'} + c y & = g (t) \\ y (0) & = y_{0} \\ y^{'} (0) & = y_{0}^{'}, \end{aligned}

the idea is to use the Laplace transform to change the differential equation into an equation that can be solved algebraically and then transform the algebraic solution back into a solution of the differential equation. Surprisingly, this method will even work when

g

is a discontinuous function, provided the discontinuities are not too bad.

🔗

Subsection 6.1.1 Definition of the Laplace Transform

We shall define the Laplace transform of a function

f (t)

🔗

L (f) (s) = F (s) = \int_{0}^{\infty} e^{- s t} f (t) d t,

provided this integral converges. The Laplace transform of a function has many nice properties, especially with respect to the derivative of

f .

However, before we investigate these properties, let us compute several Laplace transforms.

🔗

Example 6.1.1.

Let $f (t) = k,$ where $k$ is a constant. Then

$L (f) (s) = \int_{0}^{\infty} e^{- s t} k d t = \frac{k}{s},$

for $s > 0 .$
For $f (t) = t,$ the Laplace transform is

$L (f) (s) = \int_{0}^{\infty} e^{- s t} t d t = \frac{1}{s^{2}},$

for $s > 0 .$
If $f (t) = e^{a t}$ with $t > 0,$ then

$L (f) (s) = \int_{0}^{\infty} e^{- s t} e^{a t} d t = \int_{0}^{\infty} e^{- (s - a) t} d t = {[- \frac{e^{- (s - a) t}}{s - a}]}_{0}^{\infty} .$

Noting that $e^{- (s - a) t} \to 0$ as $t \to \infty$ for $s > a,$ we find that

$L (e^{a t}) = \frac{1}{s - a}$

for $s > a .$

🔗

Example 6.1.2.

The Laplace transform of a function does not always exist, even for functions that are infinitely differentiable. For example, let

f (t) = e^{t^{2}} .

Then for any real number

s,

we know that

🔗

e^{t^{2}} e^{- s t} = e^{t (t - s)} > 1,

for

t > s

with

t \geq 0;

hence,

🔗

\int_{0}^{b} e^{t^{2}} e^{- s t} d t = \int_{0}^{b} e^{t (t - s)} d t \geq b .

Thus, the integral

🔗

\int_{0}^{\infty} e^{t^{2}} e^{- s t} d t

does not converge.

🔗

The Laplace transform, however, does exist in many cases. In Subsection 6.1.3, we will show that the Laplace transform of a function exists provided the function does not grow too quickly and does not possess bad discontinuities.

🔗

Subsection 6.1.2 Properties of the Laplace Transform

One of the most important properties of the Laplace transform is linearity. That is,

🔗

L (α f + β g) (s) = α L (f) (s) + β L (g) (s),

where

α, β \in R

provided the Laplace transforms of

f

and

g

exist. The proof of this statement follows directly from the definition of the Laplace transform and the properties of integration,

🔗

\begin{aligned} L (α f + β g) (s) & = \int_{0}^{\infty} e^{- s t} [α f (t) + β g (t)] d t \\ = lim_{b \to \infty} \int_{0}^{b} e^{- s t} [α f (t) + β g (t)] d t \\ = α lim_{b \to \infty} \int_{0}^{b} e^{- s t} f (t) d t + β lim_{b \to \infty} \int_{0}^{b} e^{- s t} β g (t) d t \\ = α L (f) (s) + β L (g) (s) \end{aligned}

We state the linearity of the Laplace transform as a theorem.

Theorem 6.1.3.

Let

f

and

g

be two functions whose Laplace transforms

L (f) (s)

and

L (g) (s)

exist for

s > a

and

s > b,

respectively. Then for any real constants

α

and

β

and

s > max {a, b}

🔗

L (α f + β g) (s) = α L (f) (s) + β L (g) (s) .

🔗

Transforms of functions having the linearity property are called linear operators.

🔗

Example 6.1.4.

Recall that

\sinh a t = (e^{a t} - e^{- a t}) / 2 .

Provided

a > 0,

linearity makes it very easy to compute the Laplace transform of

\sinh a t .

By Example 6.1.1,

🔗

L (\sinh a t) = \frac{1}{2} L (e^{a t}) - \frac{1}{2} L (e^{- a t}) = \frac{1}{2} (\frac{1}{s - k} - \frac{1}{s + k}) .

Thus, for

s > k > 0,

🔗

L (\sinh a t) = \frac{k}{s^{2} - k^{2}} .

🔗

The Laplace transform of discontinuous functions also exist, provided the discontinuities are not too bad. We say that a function

f

is piecewise continuous on an interval

[a, b]

f

satisfies the following conditions.

🔗

There are a finite number of discontinuities of $f$ on $[a, b] .$
If $c \in [a, b]$ is a discontinuity, then the one-sided limits

$lim_{t \to c^{-}} f (t) and lim_{t \to c^{+}} f (t)$

both exist. We also assume that $lim_{t \to a^{+}} f (t)$ and $lim_{t \to b^{-}} f (t)$ exist.

We say that a function is piecewise continuous on an interval

[0, \infty),

if it is piecewise continuous on the interval

[0, b]

for all

b > 0 .

The types of discontinuities that we are describing are sometimes called jump discontinuities. If a function is piecewise continuous, then the continuities are not too bad.

🔗

Example 6.1.5.

One of the simplest piecewise continuous functions is

🔗

u (t) = {\begin{cases} 0, & t < 0 \\ 1, & t \geq 0. \end{cases}

The function

u (t)

has a jump discontinuity at

t = 0 .

It follows very quickly that

L (u (t)) = 1 / s

for

s > 0 .

🔗

If we define the unit step function at

a

to be

🔗

u_{a} (t) = u (t - a) = {\begin{cases} 0, & t < a \\ 1, & t \geq a, \end{cases}

then we have a jump discontinuity at

t = a

(Figure 6.1.6). If

a > 0,

then

🔗

L (u_{a} (t)) = \int_{0}^{\infty} e^{- s t} u_{a} (t) d t = \int_{a}^{\infty} e^{- s t} d t = lim_{b \to \infty} {[- \frac{e^{- s t}}{s}]}_{a}^{b} = \frac{e^{- s a}}{s},

where

s > 0

and

a > 0 .

A function of the form

u_{a} (t) = u (t - a)

is called a step function or Heaviside function, named for the British engineer Oliver Heaviside.

🔗

🔗

Example 6.1.7.

Given a function

f (t),

consider the new function

🔗

g (t) = u_{a} (t) f (t - a) .

To obtain the function

g

from

f,

we shift the graph of the

f

to the right by

a

units and let

g (t) = 0

for

0 \leq t < a

(see Figure 6.1.8). Let us compute the Laplace transform of

g,

🔗

L (g) = \int_{0}^{\infty} g (t) e^{- s t} d t = \int_{a}^{\infty} f (t - a) e^{- s t} d t .

If we make the substitution

w = t - a,

then

🔗

L (g) = \int_{0}^{\infty} f (w) e^{- s (w + a)} d t = e^{- s a} \int_{0}^{\infty} f (w) e^{- s w} d t = e^{- s a} L (f) .

In other words, if

L (f) = F (s),

then

L (u_{a} (t) f (t - a)) = e^{- a s} F (s) .

🔗

🔗

Activity 6.1.1. Finding Laplace Transforms.

Find the Laplace transform for each of the functions below. Recall that

\cosh t = (e^{t} + e^{- t}) / 2

and

\sinh t = (e^{t} - e^{- t}) / 2

🔗

(a)

f (t) = t^{6}

🔗

(b)

f (t) = t^{7} e^{- t}

🔗

(c)

f (t) = 3 \sin^{2} (\sqrt{2} t)

🔗

(d)

f (t) = e^{3 t} \sin 2 t

🔗

(e)

f (t) = 2 t^{2} \cosh t

🔗

Subsection 6.1.3 Existence and Uniqueness of the Laplace Transform

If we are to use Laplace transforms to study differential equations, we would like to know which functions actually have Laplace transforms. Furthermore, if two functions have the same Laplace transform, we can ask if the functions must be the same. In other words, we wish to know if the Laplace transform of a function exists and is unique. We can answer both of these questions affirmatively if the function

f (t)

is piecewise continuous on

[0, \infty)

and does not grow too quickly as

t \to \infty .

We say a function

f (t)

is exponentially bounded on

[0, \infty)

if there exist constants

M \geq 0

and

a

such that

🔗

| f (t) | \leq M e^{a t},

for all

t

[0, \infty) .

In other words, the graph of

f

must lie between the curves

y = M e^{a t}

and

y = - M e^{a t} .

The following two theorems tell us that Laplace transforms exist and are unique for piecewise continuous, exponentially bounded functions on the interval

[0, \infty) .

We leave the proofs of these theorem as exercises.

🔗

Theorem 6.1.9.

f

is a piecewise continuous, exponentially bounded function defined

[0, \infty),

then the Laplace transform of

f,

🔗

L (f) (s) = F (s) = \int_{0}^{\infty} e^{- s t} f (t) d t

exists.

🔗

Theorem 6.1.10.

Let

f

and

g

be two piecewise continuous exponentially bounded functions with Laplace transforms,

F (s)

and

G (s),

respectively. Suppose that

F (s) = G (s)

for all

s > c,

where

c

is some positive number. Then

f (t) = g (t)

for all

t \geq 0 .

🔗

In light of Theorem 6.1.10, it now makes sense to define the inverse Laplace transform of a function

F (s),

which we will denote by

L^{- 1} (F (s)) (t),

as the function

f (t)

whose unique Laplace transform is

F (s) .

Furthermore, the inverse Laplace transform is linear,

🔗

L^{- 1} (α F + β G) = α L^{- 1} (F) + β L^{- 1} (G) .

🔗

Subsection 6.1.4 Finding Laplace Transforms and Inverse Transforms

To use the method of Laplace transforms effectively, we need either tables or computer software such as Sage so that we can easily find Laplace transforms and inverse Laplace transforms.

🔗

Table 6.1.11. Table of Laplace Transforms

$f (t) = L^{- 1} (F (s))$	$F (s) = L (f (t))$
1	$1 / s,$ $s > 0$
$e^{a t}$	$1 / (s - a),$ $s > a$
$t^{n},$ $n \in N$	$n! / s^{n + 1},$ $s > 0$
$t^{a},$ $a > - 1$	$Γ (a + 1) / s^{a + 1},$ $s > 0$
$\sin a t$	$a / (s^{2} + a^{2}),$ $s > 0$
$\cos a t$	$s / (s^{2} + a^{2}),$ $s > 0$
$\sinh a t$	$a / (s^{2} - a^{2}),$ $s > \| a \|$
$\cosh a t$	$s / (s^{2} - a^{2}),$ $s > \| a \|$
$e^{a t} \sin b t$	$b / [(s - a)^{2} + b^{2}],$ $s > a$
$e^{a t} \cos b t$	$(s - a) / [(s - a)^{2} + b^{2}],$ $s > a$
$t^{n} e^{a t},$ $n \in N$	$n! / (s - a)^{n + 1},$ $s > a$
$u_{c} (t)$	$e^{- c s} / s,$ $s > 0$
$u_{c} (t) f (t - c)$	$e^{- c s} F (s)$
$e^{c t} f (t)$	$F (s - c)$
$f (c t)$	$(1 / c) F (s / c),$ $c > 0$
$δ (t - c) = δ_{c} (t)$	$e^{- c s}$
$f^{'} (t)$	$s F (s) - f (0)$
$f^{″} (t)$	$s^{2} F (s) - s f (0) - f^{'} (0)$

🔗

We can also use Sage to find Laplace transforms and inverse transforms (see Subsection 6.1.8).

🔗

Activity 6.1.2. Finding Inverse Transforms.

Find the inverse Laplace transform for each of the following functions.

🔗

(a)

F (s) = \frac{2}{s^{2}}

🔗

(b)

F (s) = \frac{1}{s^{4}}

🔗

(c)

F (s) = \frac{2}{(s - 3)^{2} + 4}

🔗

(d)

F (s) = \frac{2 s^{2}}{(s^{2} + 1) (s - 1)}

🔗

(e)

F (s) = \frac{2 s - 13}{s (s^{2} - 4 s + 13)}

🔗

Subsection 6.1.5 Important Lessons

Many initial value problems have discontinuous forcing terms.
The Laplace transform of a function $f (t)$ by

$L (f) (s) = F (s) = \int_{0}^{\infty} e^{- s t} f (t) d t,$

provided the integral converges.
If $f$ is a piecewise continuous, exponentially bounded function defined $[0, \infty),$ then the Laplace transform of $f,$

$L (f) (s) = F (s) = \int_{0}^{\infty} e^{- s t} f (t) d t$

exists.
The Laplace transform is a linear operator; i.e.,

$L (α f + β g) (s) = α L (f) (s) + β L (g) (s),$

where $α, β \in R$ provided the Laplace transforms of $f$ and $g$ exist.
If we know the Laplace transform $F (s)$ of a function, we can recover the original function using the inverse Laplace transform of a function $L^{- 1} (F (s)) (t) .$
The inverse Laplace transform is linear,

$L^{- 1} (α F + β G) = α L^{- 1} (F) + β L^{- 1} (G) .$
If $L (f) = F (s),$ then $L (u_{c} (t) f (t - c)) = e^{- s c} F (s) .$
If $c \geq 0,$ we define the Heaviside function to be

$u_{c} (t) = {\begin{cases} 0 & t < c \\ 1 & t \geq c . \end{cases}$

The Laplace transform of $u_{c}$ is

$L (u_{c} (t)) (s) = \frac{e^{- c s}}{s} .$

🔗

Reading Questions 6.1.6 Reading Questions

1.

Does the Laplace transform of a function always exist? If so, why? If not, given an example?

🔗

2.

What is the Heaviside function? Explain.

🔗

Exercises 6.1.7 Exercises

Finding Laplace transforms.

Find the Laplace transform for each of the functions in Exercise Group 6.1.7.1–8. Recall that

\cosh t = (e^{t} + e^{- t}) / 2

and

\sinh t = (e^{t} - e^{- t}) / 2 .

🔗

1.

f (t) = \sin a t

🔗

2.

f (t) = \cos a t

🔗

3.

f (t) = e^{a t} \sin b t

🔗

4.

f (t) = \cosh a t

🔗

5.

f (t) = 2 \cos a t - 3 \sinh a t

🔗

6.

f (t) = t e^{a t}

🔗

7.

f (t) = t \sin a t

🔗

8.

f (t) = t^{2} \cos a t

🔗

Finding inverse Laplace transforms.

Find the inverse Laplace transform for each of the functions in Exercise Group 6.1.7.9–16. You will find partial fraction decomposition very useful.

🔗

9.

F (s) = \frac{1}{s (s + 1)}

🔗

10.

F (s) = s^{- 3 / 2}

🔗

11.

F (s) = \frac{1}{s + 7}

🔗

12.

F (s) = \frac{2}{s^{2} + 9}

🔗

13.

F (s) = \frac{5 s - 4}{2 s^{2} + s - 1}

🔗

14.

F (s) = \frac{2 s^{2} - s + 4}{s^{3} + 4 s}

🔗

15.

F (s) = \frac{1}{(s + 2)^{3}}

🔗

16.

F (s) = \frac{1}{(s + 1)^{2} (s^{2} - 9)}

🔗

17.

Prove Theorem 6.1.9: If

f

is a piecewise continuous, exponentially bounded function defined

[0, \infty),

then the Laplace transform of

f,

🔗

L (f) (s) = F (s) = \int_{0}^{\infty} e^{- s t} f (t) d t

exists.

🔗

18.

Let

f

be a function whose Laplace transform

L (f) (s) = F (s)

exists for

s > a .

Prove that

🔗

L (e^{c t} f (t)) (s) = F (s - c)

for

s > a + c .

🔗

19.

Prove Theorem 6.1.10: Suppose that

f

and

g

are piecewise continuous, exponentially bounded functions defined

[0, \infty),

with Laplace transforms of

F (s)

and

G (s),

respectively. If

F (s) = G (s),

then

f (t) = g (t) .

🔗

20.

f

is a piecewise continuous, exponentially bounded function defined

[0, \infty),

prove that the Laplace transform of

f,

🔗

L (f) (s) = F (s) = \int_{0}^{\infty} e^{- s t} f (t) d t

exists.

🔗

21.

Define the gamma function to be

🔗

Γ (x) = \int_{0}^{\infty} e^{- t} t^{x - 1} d t

for

x > 0 .

The gamma function is very useful for expressing the Laplace transform of the function

t^{α} .

🔗

Show that $Γ (1) = 1 .$
Prove that $Γ (x + 1) = x Γ (x)$ and deduce that $Γ (n + 1) = n!$ for $n \in N .$
Show that the Laplace transform of $f (t) = t^{α}$ is

$L (f) (s) = \frac{Γ (α + 1)}{s^{α + 1}} .$

If $n$ is a positive integer, then $L (t^{n}) (s) = n! / s^{n + 1} .$

🔗

Subsection 6.1.8 Laplace Transforms with Sage

Computer algebra systems have now replaced tables of Laplace transforms just as the calculator has replaced the slide rule. It is easy to calculate Laplace transforms with Sage. For example, suppose that we wish to compute the Laplace transform of

f (x) = t^{3} e^{t} - \cos t .

We can use the Sage command laplace.

🔗

xxxxxxxxxx
 
s = var("s")
t = var("t")
f = t^3*exp(t) - cos(t)
f.laplace(t,s)

To recover the original function, we can use the Sage command inverse_laplace. Suppose that

🔗

F (s) = - \frac{s}{s^{2} + 1} + \frac{6}{(s - 1)^{4}} .

🔗

xxxxxxxxxx
 
s = var("s")
t = var("t")
F = -s/(s^2 + 1) + 6/(s - 1)^4
F.inverse_laplace(s,t)

We can even use Sage to find the Laplace transform of piecewise-defined functions. Suppose that

🔗

f (t) = {\begin{cases} 1, & 0 \leq t \leq 3 \\ 0, & t > 3. \end{cases}

The Sage command to define a piecewise-defined function is piecewise.

🔗

xxxxxxxxxx
 
t = var("t")
f1(t) = 1
f2(t) = 0
f = piecewise([[[0, 3], f1], [(3, infinity), f2]])
f.plot((t, 0, 6), thickness=3, color="red", axes_labels=['$t$','$f(t)$'])

Sage can be used to find the Laplace transform of a piecewise-defined function.

🔗

xxxxxxxxxx
 
s = var("s")
t = var("t")
f1(t) = 1
f2(t) = 0
f = piecewise([[[0, 3], f1], [(3, infinity), f2]])
f.laplace(t,s)

🔗

Exercises 6.1.9 Exercises

1.

Use Sage to find the Laplace transform of

🔗

f (t) = {\begin{cases} - t, & 0 \leq t \leq 1 \\ 2, & t > 1. \end{cases}

🔗

xxxxxxxxxx
 
​

🔗

2.

Use Sage to find the Laplace transforms for each of the functions in Exercise Group 6.1.7.1–8.

🔗

xxxxxxxxxx
 
​

🔗

3.

Use Sage to find the inverse Laplace transform for each of the functions in Exercise Group 6.1.7.9–16.

🔗

xxxxxxxxxx
 
​

🔗

Prev Top Next