Lecture Notes for Laplace Transform

Wen Shen
April 2009
NB! These notes are used by myself. They are provided to students as a supplement to the
textbook. They can not substitute the textbook.
Laplace Transform is used to handle piecewise continuous or impulsive force.

6.1: Definition of the Laplace transform (1)

Topics:
Definition of Laplace transform,
Compute Laplace transform by definition, including piecewise continuous functions.
Definition: Given a function f (t), t 0, its Laplace transform F (s) = L{f (t)} is defined as
Z
Z A
.
.
st
F (s) = L{f (t)} =
e f (t) dt = lim
est f (t) dt
A

We say the transform converges if the limit exists, and diverges if not.
Next we will give examples on computing the Laplace transform of given functions by definition.
Example 1. f (t) = 1 for t 0.
A
Z A
1
1
1 st
st
F (s) = L{f (t)} = lim
e 1 dt = lim e = lim esA 1 = ,
A
A 0
A
s
s
s
0

(s > 0)

Example 2. f (t) = et .
Z
F (s) = L{f (t)} = lim

lim

e
0

st at

e dt = lim

1 (sa)A
1
e
1 =
,
sa
sa
1

(sa)t

(s > a)

A
1 (sa)t
dt = lim
e

A
sa
0

Example 3. f (t) = tn , for n 1 integer.

Z
F (s) =

lim

(
est tn dt =

n
= 0 + lim
s A

lim

est tn1 dt =

Z A n1 st )
st A
e
nt e

tn
dt

s 0
s
0
n
L{tn1 }.
s

So we get a recursive relation

L{tn } =
which means
L{tn1 } =

n
L{tn1 },
s

n1
L{tn2 },
s

n,

L{tn2 } =

n2
L{tn3 },
s

By induction, we get
n
n (n 1)
n (n 1) (n 2)
L{tn1 } =
L{tn2 } =
L{tn3 }
s
s
s
s
s
s
n (n 1) (n 2)
1
n! 1
n!
= =
L{1} = n = n+1 , (s > 0)
s
s
s
s
s s
s

L{tn } =

Example 4. Find the Laplace transform of sin at and cos at.

Method 1. Compute by definition, with integration-by-parts, twice. (lots of work...)
Method 2. Use the Eulers formula
eiat = cos at + i sin at,

L{eiat } = L{cos at} + iL{sin at}.

By Example 2 we have
L{eiat } =

1
1(s + ia)
s + ia
s
a
=
= 2
= 2
+i 2
.
2
2
s ia
(s ia)(s + ia)
s +a
s +a
s + a2

Comparing the real and imaginary parts, we get

s
a
L{cos at} = 2
,
L{sin at} = 2
, (s > 0).
2
s +a
s + a2
RA
R
Remark: Now we will use 0 instead of limA 0 , without causing confusion.
For piecewise continuous functions, Laplace transform can be computed by integrating each
integral and add up at the end.
Example 5. Find the Laplace transform of

1,
0 t < 2,
f (t) =
t 2, 2 t.
2

We do this by definition:
Z
Z 2
Z
st
st
F (s) =
e f (t) dt =
e dt +
(t 2)est dt
0
0
2
2 Z A

1
1 st
1 st
=
e + (t 2) est
e dt
s
s
2 s
t=0
t=2

1 2s
1 1 st
1 2s
1
=
(e 1) + (0 0) +
e
=
(e 1) + 2 e2s
s
s s
s
s
t=2

6.2: Solution of initial value problems (4)

Topics:
Properties of Laplace transform, with proofs and examples
Inverse Laplace transform, with examples, review of partial fraction,
Solution of initial value problems, with examples covering various cases.
Properties of Laplace transform:
1. Linearity: L{c1 f (t) + c2 g(t)} = c1 L{f (t)} + c2 L{g(t)}.
2. First derivative: L{f 0 (t)} = sL{f (t)} f (0).
3. Second derivative: L{f 00 (t)} = s2 L{f (t)} sf (0) f 0 (0).
4. Higher order derivative:
L{f (n) (t)} = sn L{f (t)} sn1 f (0) sn2 f 0 (0) sf (n2) (0) f (n1) (0).
5. L{tf (t)} = F 0 (s) where F (s) = L{f (t)}. This also implies L{tf (t)} = F 0 (s).
6. L{eat f (t)} = F (s a) where F (s) = L{f (t)}. This implies eat f (t) = L1 {F (s a)}.
Remarks:
Note property 2 and 3 are useful in differential equations. It shows that each derivative
in t caused a multiplication of s in the Laplace transform.
Property 5 is the counter part for Property 2. It shows that each derivative in s causes
a multiplication of t in the inverse Laplace transform.
Property 6 is also known as the Shift Theorem. A counter part of it will come later in
chapter 6.3.
Proof:
1. This follows by definition.
2. By definition
Z

L{f (t)} =

e
0

st 0

st

f (t)dt = e

f (t)
0

(s)est f (t)dt = f (0) + sL{f (t)}.

3. This one follows from Property 2. Set f to be f 0 we get

L{f 00 (t)} = sL{f 0 (t)} f 0 (0) = s(sL{f (t)} f (0)) f 0 (0) = s2 L{f (t)} sf (0) f 0 (0).
4. This follows by induction, using property 2.
5. The proof follows from the definition:
Z
Z
Z
d st
st
0
F (s) =
e f (t)dt =
(e )f (t)dt =
(t)est f (t)dt = L{tf (t)}.
ds 0
s
0
0
6. This proof also follows from definition:
Z
Z
at
st at
L{e f (t)}
e e f (t)dt =
0

e(sa)t f (t)dt = F (s a).

By using these properties, we could find more easily Laplace transforms of many other functions.
Example 1.
From

L{tn } =

n!

n!
.
(s a)n+1

we get

L{eat tn } =

s2

b
,
+ b2

we get

L{eat sin bt} =

b
.
(s a)2 + b2

s2

s
,
+ b2

we get

L{eat cos bt} =

sa
.
(s a)2 + b2

sn+1

Example 2.
From

L{sin bt} =

Example 3.
From

L{cos bt} =

Example 4.
L{t3 + 5t 2} = L{t3 } + 5L{t} 2L{1} =

Example 5.
L{e2t (t3 + 5t 2)} =

3!
1
1
+5 2 2 .
4
s
s
s

1
1
3!
+5
2
.
4
2
(s 2)
(s 2)
s2
5

Example 6.
L{(t2 + 4)e2t et cos t} =
because
L{t2 + 4} =

2
4
+ ,
3
s
s

2
4
s+1
+

,
3
(s 2)
s 2 (s + 1)2 + 1

L{(t2 + 4)e2t } =

2
4
+
.
3
(s 2)
s2

Next are a few examples for Property 5.

Example 7.
1
Given L{e } =
,
sa
at

at

we get

Example 8.

L{te } =

L{t sin bt} =

b
2
s + b2

Example 9.
L{t cos bt} =

s
2
s + b2

0
=

1
sa

2bs
.
+ b2 )2

(s2

s2 b2
= = 2
.
(s + b2 )2

0
=

1
(s a)2

Inverse Laplace transform. Definition:

L1 {F (s)} = f (t),

if

F (s) = L{f (t)}.

Technique: find the way back.

Some simple examples:
Example 10.

3
s2 + 4

=L

3
2
2
2 s + 22

3
= L1
2

2
s2 + 22

3
sin 2t.
2

Example 11.

2
1
6
1 1
3!
1 5t 1 3!
1
1
1
L
=L

= L
= e L
= e5t t3 .
4
4
4
4
(s + 5)
3 (s + 5)
3
(s + 5)
3
s
3
Example 12.

s+1
s2 + 4

=L

s
2
s +4

1
+ L1
2

2
2
s +4

1
= cos 2t sin 2t.
2

Example 13.

s+1
s+1
3/4
1/4
3
1
1
1
1
L
=L
=L
+
= e2t + e2t .
2
s 4
(s 2)(s + 2)
s2 s+2
4
4
Here we used partial fraction to find out:
s+1
A
B
=
+
,
(s 2)(s + 2)
s2 s+2

A = 3/4,

B = 1/4.

Solutions of initial value problems.

We will go through one example first.
Example 14. (Two distinct real roots.) Solve the initial value problem by Laplace transform,
y 00 3y 0 10y = 2,

y(0) = 1, y 0 (0) = 2.

Step 1. Take Laplace transform on both sides: Let L{y(t)} = Y (s), and then
L{y 0 (t)} = sY (s) y(0) = sY 1,

L{y 00 (t)} = s2 Y (s) sy(0) y 0 (0) = s2 Y s 2.

Note the initial conditions are the first thing to go in!

L{y 00 (t)} 3L{y 0 (t)} 10L{y(t)} = L{2},

2
s2 Y s 2 3(sY 1) 10Y = .
s

Now we get an algebraic equation for Y (s).

Step 2: Solve it for Y (s):
(s2 3s 10)Y (s) =

2
s2 s + 2
+s+23=
,
s
s

Y (s) =

s2 s + 2
.
s(s 5)(s + 2)

Step 3: Take inverse Laplace transform to get y(t) = L1 {Y (s)}. The main technique here is
partial fraction.
Y (s) =

s2 s + 2
A
B
C
A(s 5)(s + 2) + Bs(s + 2) + Cs(s 5)
= +
+
=
.
s(s 5)(s + 2)
s
s5 s+2
s(s 5)(s + 2)

Compare the numerators:

s2 s + 2 = A(s 5)(s + 2) + Bs(s + 2) + Cs(s 5).
The previous equation holds for all values of s.
Set s = 0: we get 10A = 2, so A = 15 .
Set s = 5: we get 35B = 22, so B = 22
.
35
4
Set s = 2: we get 14C = 8, so C = 7 .
Now, Y (s) is written into sum of terms which we can find the inverse transform:
1
1
1
1 22
4
y(t) = AL1 { } + BL1 {
} + CL1 {
} = + e5t + e2t .
s
s5
s+2
5 35
7

Structure of solutions:
Take Laplace transform on both sides. You will get an algebraic equation for Y .
Solve this equation to get Y (s).
Take inverse transform to get y(t) = L1 {y}.
Example 15. (Distinct real roots, but one matches the source term.) Solve the initial value
problem by Laplace transform,
y 00 y 0 2y = e2t ,

y(0) = 0, y 0 (0) = 1.

Take Laplace transform on both sides of the equation, we get

L{y 00 } L{y 0 } L{2y} = L{e2t },

s2 Y (s) 1 sY (s) 2Y (s) =

1
.
s2

Solve it for Y :
(s2 s2)Y (s) =

1
s1
+1 =
,
s2
s2

Y (s) =

s1
s1
=
.
(s 2)(s2 s 2)
(s 2)2 (s + 1)

Use partial fraction:

s1
A
B
C
=
+
+
.
2
(s 2) (s + 1)
s + 1 s 2 (s 2)2
Compare the numerators:
s 1 = A(s 2)2 + B(s + 1)(s 2) + C(s + 1)
Set s = 1, we get A = 29 . Set s = 2, we get C = 13 . Set s = 0 (any convenient values of s
can be used in this step), we get B = 92 . So
Y (s) =
and

2 1
2 1
1
1
+
+
9 s + 1 9 s 2 3 (s 2)2

2
2
1
y(t) = L1 {Y } = et + e2t + te2t .
9
9
3

Compare this to the method of undetermined coefficient: general solution of the equation
should be y = yH + Y , where yH is the general solution to the homogeneous equation and
Y is a particular solution. The characteristic equation is r2 r 2 = (r + 1)(r 2) = 0,
so r1 = 1, r2 = 2, and yH = c1 et + c2 e2t . Since 2 is a root, so the form of the particular
solution is Y = Ate2t . This discussion concludes that the solution should be of the form
y = c1 et + c2 e2t + Ate2t
for some constants c1 , c2 , A. This fits well with our result.
9

Example 16. (Complex roots.)

y 00 2y 0 + 2y = et ,

y 0 (0) = 1.

y(0) = 0,

Before we solve it, lets use the method of undetermined coefficients to find out which terms
will be in the solution.
r2 2r + 2 = 0,

(r 1)1 + 1 = 0,

yH = c1 et cos t + c2 et sin t,

r1,2 = 1 i,

Y = Aet ,

so the solution should have the form:

y = yH + Y = c1 et cos t + c2 et sin t + Aet .
The Laplace transform would be
Y (s) = c1

s1
1
c1 (s 1) + c2
A
1
+
c
+
A
=
+
.
2
(s 1)2 + 1
(s 1)2 + 1
s+1
(s 1)2 + 1
s+1

This gives us idea on which terms to look for in partial fraction.

Now lets use the Laplace transform:
Y (s) = L{y},

L{y 0 } = sY y(0) = sY,

1
,
s+1

L{y 00 } = s2 Y sy(0) y(0) = s2 Y 1.

1
s+2
+1=
s+1
s+1
s+2
s+2
A
B(s 1) + C
Y (s) =
=
=
+
2
2
(s + 1)(s 2s + 2)
(s + 1)((s 1) + 1)
s+1
(s 1)2 + 1

s2 Y 1 2sY + 2Y =

(s2 2s + 2)Y (s) =

Compare the numerators:

s + 2 = A((s 1)2 + 1) + (B(s 1) + C)(s + 1).
Set s = 1: 5A = 1, A = 15 .
Compare coefficients of s2 -term: A + B = 0, B = A = 51 .
Set any value of s, say s = 0: 2 = 2A B + C, C = 2 2A + B = 95 .
Y (s) =

1 1
1
9
s1
1

+
2
5 s + 1 5 (s 1) + 1 5 (s 1)2 + 1
1
1
9
y(t) = et et cos t + et sin t.
5
5
5

We see this fits our prediction.

10

Example 17. (Pure imaginary roots.)

y 00 + y = cos 2t,

y(0) = 2,

y 0 (0) = 1.

Again, lets first predict the terms in the solution:

r2 + 1 = 0,

r1,2 = i,

yH = c1 cos t + c2 sin t,

Y = A cos 2t

so
y = yH + Y = c1 cos t + c2 sin t + A cos 2t,
and the Laplace transform would be
Y (s) = c1

s1

s
s
1
+ c2 2
+A 2
.
+1
s +1
s +4

Now, lets take Laplace transform on both sides:

s2 Y 2s 1 + Y = L{cos 2t} =

s2

s
+4

s
2s3 + s2 + 9s + 4
+
2s
+
1
=
s2 + 4
s2 + 4
2s3 + s2 + 9s + 4
As + B Cs + D
Y (s) =
= 2
+ 2
.
2
2
(s + 4)(s + 1)
s +1
s +4

(s2 + 1)Y (s) =

Comparing numerators, we get

2s3 + s2 + 9s + 4 = (As + B)(s2 + 4) + (Cs + D)(s2 + 1).
One may expand the right-hand side and compare terms to find A, B, C, D, but that takes
more work. Lets try by setting s into complex numbers.
Set s = i, and remember the facts i2 = 1 and i3 = i, we have
2i 1 + 9i + 4 = (Ai + B)(1 + 4),

3 + 7i = 3B + 3Ai,

7
B = 1, A = .
3

Set now s = 2i:

16i 4 + 18i + 4 = (2Ci + D)(3),
So
Y (s) =

0 + 2i = 3D 6Ci,

7 s
1
1 s
+ 2
2
2
3s +1 s +1 3s +4

and
y(t) =

7
1
cos t + sin t cos 2t.
3
3

11

1
D = 0, C = .
3

A very brief review on partial fraction, targeted towards inverse Laplace transform.
Goal: rewrite a fractional form PPmn (s)
(where Pn is a polynomial of degree n) into sum of
(s)
simpler terms. We assume n < m.
The type of terms appeared in the partial fraction is solely determined by the denominator
Pm (s). First, fact out Pm (s), write it into product of terms of (i) s a, (ii) s2 + a2 , (iii)
(sa )2 + b2 . The following table gives the terms in the partial fraction and their corresponding
inverse Laplace transform.
term in PM (s)

from where?
real root, or

sa

g(t) = eat

term in partial fraction

inverse L.T.

A
sa

Aeat

A
B
+
s a (s a)2

Aeat + Bteat

double roots,
(s a)2

or r = a and g(t) = eat

double roots,

(s a)

and g(t) = eat

A
B
C
+
+
2
s a (s a)
(s a)3

imaginary roots or
s 2 + 2

g(t) = cos t or sin t

Aeat + Bteat +

C 2 at
te
2

As + B
s2 + 2

A cos t + B sin t

A(s ) + B
(s )2 + 2

et (A cos t + B sin t)

complex roots, or
(s )2 + 2

g(t) = et cos t(or sin t)

In summary, this table can be written

Pn (s)
(s a)(s
c)3 ((s )2 + 2 )
B1
B2
C1
C2
C3
D1 (s ) + D2
A
+
+
+
+
+
+
.
=
2
2
3
s a s b (s b)
s c (s c)
(s c)
(s )2 + 2
b)2 (s

12

6.3: Step functions (2)

Topics:
Definition and basic application of unit step (Heaviside) function,
Laplace transform of step functions and functions involving step functions (piecewise
continuous functions),
Inverse transform involving step functions.
We use steps functions to form piecewise continuous functions.
Unit step function(Heaviside function):

uc t =

0, 0 t < c,
1, c t.

for c 0. A plot of uc (t) is below:

u

For a given function f (t), if it is multiplied with uc (t), then

0,
0 < t < c,
uc tf (t) =
f (t), c t.
We say uc picks up the interval [c, ).
Example 1. Consider

1 uc (t) =

1, 0 t < c,
0, c t.

A plot of this is given below

1 u

13

We see that this function picks up the interval [0, c).

Example 2. Rectangular pulse. The plot of the function looks like
u u
a b
1

for 0 a < b < . We see it can be expressed as

ua (t) ub (t)
and it picks up the interval [a, b).
Example 3. For the function

g(t) =

f (t), a t < b
0,
otherwise

We can rewrite it in terms of the unit step function as

g(t) = f (t) ua (t) ub (t) .

Example 4. For the function

sin t, 0 t < 1,
et ,
1 t < 5,
ft =
2
t
5 t,
we can rewrite it in terms of the unit step function as we did in Example 3, treat each interval
separately

t
f (t) = sin t u0 (t) u1 (t) + e u1 (t) u5 (t) + t2 u5 (t).

14

Laplace transform of uc (t): by definition

Z
Z
est
est
esc
st
st
L{uc (t)} =
e uc (t) dt =
e 1 dt =
=
0

=
,
s t=c
s
s
0
c

(s > 0).

Shift of a function: Given f (t), t > 0, then

f (t c), c t,
g(t) =
0,
0 t < c,
is the shift of f by c units. See figure below.
f

Let F (s) = L{f (t)} be the Laplace transform of f (t). Then, the Laplace transform of g(t) is
Z
Z
st
L{g(t)} = L{uc (t) f (t c)} =
e uc (t)f (t c) dt =
est f (t c) dt.
0

Let y = t c, so t = y + c, and dt = dy, and we continue

Z
Z
s(y+c)
sc
L{g(t)} =
e
f (y) dy = e
esy f (y) dy = ecs F (s).
0

So we conclude
L{uc (t)f (t c)} = ecs L{f (t)} = ecs F (s) ,
which is equivalent to
L1 {ecs F (s)} = uc (t)f (t c) .
Note now we are only considering the domain t 0. So u0 (t) = 1 for all t 0.

15

In following examples we will compute Laplace transform of piecewise continuous functions

with the help of the unit step function.
Example 5. Given

f (t) =

sin t,
0 t < 4 ,
sin t + cos(t 4 ), 4 t.

It can be rewritten in terms of the unit step function as

f (t) = sin t + u 4 (t) cos(t

).
4

(Or, if we write out each intervals

f (t) = sin t(1 u 4 (t)) + (sin t + cos(t

))u 4 (t) = sin t + u 4 (t) cos(t ).

4
4

which gives the same answer.)

And the Laplace transform of f is
F (s) = L{sin t} + L{u 4 (t) cos(t

Example 6. Given

f (t) =

1
s

)} = 2
+ e 4 s 2
.
4
s +1
s +1

t, 0 t < 1,
1, 1 t.

It can be rewritten in terms of the unit step function as

f (t) = t(1 u1 (t)) + 1 u1 (t) = t (t 1)u1 (t) .
The Laplace transform is
L{f (t)} = L{t} L{(t 1)u1 (t)} =

Example 7. Given

f (t) =

1
s 1

e
.
s2
s2

0,
0 t < 2,
t + 3, 2 t.

We can rewrite it in terms of the unit step function as

f (t) = (t + 3)u2 (t) = (t 2 + 5)u2 (t) = (t 2)u2 (t) + 5u2 (t) .
The Laplace transform is
L{f (t)} = L{(t 2)u2 (t)} + 5L{u2 (t)} = e2s
16

1
1
+ 5e2s .
2
s
s

Example 8. Given

g(t) =

1, 0 t < 2,
t2 , 2 t.

We can rewrite it in terms of the unit step function as

g(t) = 1 (1 u2 (t)) + t2 u2 (t) = 1 + (t2 1)u2 (t) .
Observe that
t2 1 = (t 2 + 2)2 1 = (t 2)2 + 4(t 2) + 4 1 = (t 2)2 + 4(t 2) + 3 ,
we have

g(t) = 1 + (t 2)2 + 4(t 2) + 3 u2 (t) .

The Laplace transform is

1
L{g(t)} = + e2s
s
Example 9. Given

4
3
2
+ 2+
3
s
s
s

0, 0 t < 3,
et , 3 t < 4,
f (t) =

0, 4 t.

We can rewrite it in terms of the unit step function as

f (t) = et (u3 (t) u4 (t)) = u3 (t)et3 e3 u4 (t)et4 e4 .
The Laplace transform is
L{g(t)} = e3 e3s

1
1
1 3(s1)
e4 e4s
=
e
e4(s1) .
s1
s1
s1

17

Inverse transform: We use two properties:

1
L{uc (t)} = ecs ,
s

and

L{uc (t)f (t c)} = ecs L{f (t)} .

In the following examples we want to find f (t) = L1 {F (s)}.

Example 10.

1 e2s
1
1
= 3 e2s 3 .
3
s
s
s
1 2
1 1
We know that L { s3 } = 2 t , so we have
( 1 2
t,
0 t < 2,
1
1
2
f (t) = L1 {F (s)} = t2 u2 (t) (t 2)2 =
.
1 2
2
2
t 12 (t 2)2 , 2 t.
2
F (s) =

Example 11. Given

e3s
1
F (s) = 2
= e3s
= e3s
s + s 12
(s + 4)(s + 3)

A
B
+
s+4 s3

By partial fraction, we find A = 17 and B = 17 . So

f (t) = L1 {F (s)} = u3 (t) Ae4(t3) + Be3(t3) = u3 (t) e4(t3) + e3(t3)

7
which can be written as a p/w continuous function

0,
f (t) =
17 e4(t3) + 17 e3(t3) ,

0 t < 3,
3 t.

Example 12. Given

s+22
ses
s+22
s s + 2 2
s
F (s) = 2
=e
=s
+
.
s + 4s + 5
(s + 2)2 + 1
(s + 2)2 + 1 (s + 2)2 + 1
So

f (t) = L1 {F (s)} = u1 (t) e2(t1) cos(t 1) 2e2(t1) sin(t 1)

which can be written as a p/w continuous function

0,
f (t) =
e2(t1) [cos(t 1) 2 sin(t 1)] ,

18

0 t < 1,
1 t.

6.4: Differential equations with discontinuous forcing

functions (1)
Topics:
Solve initial value problems with discontinuous force, examples of various cases,
Describe behavior of solutions, and make physical sense of them.
Next we study initial value problems with discontinuous force. We will start with an example.
Example 1. (Damped system with force, complex roots) Solve the following initial value
problem

0, 0 t < 1,
00
0
y + y + y = g(t),
g(t) =
,
y(0) = 1, y 0 (0) = 0 .
1, 1 t,
Let L{y(t)} = Y (s), so L{y 0 } = sY 1 and L{y 00 } = s2 Y s. Also we have L{g(t)} =
L{u1 (t)} = es 1s . Then
1
s2 Y s + sY 1 + Y = es ,
s
which gives
s+1
es
+ 2
.
Y (s) =
2
s(s + s + 1) s + s + 1
Now we need to find the inverse Laplace transform for Y (s). We have to do partial fraction
first. We have
1
A
Bs + C
=
+
.
s(s2 + s + 1)
s
s2 + s + 1
Compare the numerators on both sides:
1 = A(s2 + s + 1) + (Bs + C) s
Set s = 0, we get A = 1.
Compare s2 -term: 0 = A + B, so B = A = 1.
Compare s-term: 0 = A + C, so C = A = 1.
So

Y (s) = e

s+1
1
2
s s +s+1

s2

s+1
.
+s+1

We work out some detail

(s + 21 ) + 13 23
s+1
s+1

=
=
,
s2 + s + 1
(s + 12 )2 + ( 23 )2
(s + 12 )2 + ( 23 )2
19

so

L1

We conclude
"
y(t) = u1 (t) 1 e

s+1
2
s +s+1

21 (t1)

12 t

=e

!
3
1
3
cos
t + sin
t .
2
2
3

!#
"

1
3
3
3
1
3
cos
(t 1) sin
(t 1) +e 2 t cos
t + sin
t .
2
2
2
2
3

Remark: There are other ways to work out the partial fractions.
Extra question: What happens when t ?
Answer: We see all the terms with the exponential function will go to zero, so y 1 in the
limit. We can view this system as the spring-mass system with damping. Since g(t) becomes
constant 1 for large t, and the particular solution (which is also the steady state) with 1 on
the right hand side is 1, which provides the limit for y.
Further observation:
We see that the solution to the homogeneous equation is
"

#
1
3
3
e 2 t c1 cos
t + c2 sin
t ,
2
2
and these terms do appear in the solution.
Actually the solution consists of two part: the forced response and the homogeneous
solution.
Furthermore, the g has a discontinuity at t = 1, and we see a jump in the solution also
for t = 1, as in the term u1 (t).

20

Example 2. (Undamped system with force, pure imaginary roots) Solve the following initial
value problem

0, 0 t < ,
1, t < 2,
y 00 + 4y = g(t) =
y(0) = 1, y 0 (0) = 0 .

0, 2 t,
Rewrite

1
1
L{g} = es e2 .
s
s

g(t) = u (t) u2 (t),

So
s2 Y s + 4Y =

1
e e2 .
s

Solve it for Y :
Y (s) =

e e2
s
e
e2
s
+
=

+
.
s(s2 + 4)
s2 + 4
s(s2 + 4) s(s2 + 4) s2 + 4

Work out partial fraction

1
A Bs + C
=
+ 2
,
s(s2 + 4)
s
s +4
So
L1 {

1
A= ,
4

1
B= ,
4

C = 0.

1
1 1
} = cos 2t .
+ 4)
4 4

s(s2

Now we take inverse Laplace transform of Y

1 1
1 1
y(t) = u (t)( cos 2(t )) u2 (t)( cos 2(t 2)) + cos 2t
4 4
4 4
1
= (u (t) u2 ) (1 cos 2t) + cos 2t
1 4
(1 cos 2t), t < 2,
= cos 2t + 4
0,
otherwise,
= homogeneous solution + forced response

21

Example 3. In Example 14, let

0, 0 t < 4,
et , 4 5 < 2,
g(t) =

0, 5 t.
Find Y (s).
Rewrite
g(t) = et (u4 (t) u5 (t)) = u4 (t)et4 e4 u5 (t)et5 e5 ,
so
G(s) = L{g(t)} = e4 e4s

1
1
e5 e5s
.
s1
s1

Take Laplace transform of the equation, we get

(s2 + 4)Y (s) = G(s) + s,

Y (s) = e4 e4s e5 e5s

1
s
+ 2
.
2
(s 1)(s + 4) s + 4

Remark: We see that the first term will give the forced response, and the second term is from
the homogeneous equation.
The students may work out the inverse transform as a practice.

22

Example 4. (Undamped system with force, example 2 from the book p. 334)

0 t < 5,
0,
00
0
(t 5)/5, 5 5 < 10,
y + 4y = g(t), y(0) = 0, y (0) = 0, g(t) =

1,
10 t.
Lets first work on g(t) and its Laplace transform
g(t) =

t5
1
1
(u5 (t) u10 (t)) + u10 (t) = u5 (t)(t 5) u10 (t)(t 10),
5
5
5
1
1
1
1
G(s) = L{g} = e5s 2 e10s 2
5
s
5
s

Let Y (s) = L{y}, then

(s2 + 4)Y (s) = G(s),

Y (s) =

G(s)
1 5s
1
1 10s
1
=
e

e
s2 + 4
5
s2 (s2 + 4) 5
s2 (s2 + 4)

Work out the partial fraction:

.
H(s) =

1
A B Cs + 2D
= + 2+ 2
+ 4)
s
s
s +4

s2 (s2

one gets A = 0, B = 41 , C = 0, D = 81 . So

1
1 1
1
2
1
1
. 1
1
h(t) = L
=L
2 2
= t sin 2t.
2
2
2
s (s + 4)
4 s
8 s +2
4
8
Go back to y(t)
1
1
y(t) = L1 {Y } = u5 (t)h(t 5) u10 (t)h(t 10)
5
5

1
1
1
1
1
1
=
u5 (t) (t 5) sin 2(t 5) u10 (t) (t 10) sin 2(t 10)
5
4
8
5
4
8

0 t < 5,
0,
1
1
(t 5) 40 sin 2(t 5),
5 5 < 10,
=
20
1
1

(sin
2(t

5)

sin
2(t

10)),
10 t.
4
40
Note that for t 10, we have y(t) =

1
4

+ R cos(2t + ) for some amplitude R and phase .

The plots of g and y are given in the book. Physical meaning and qualitative nature of the
solution:
The source g(t) is known as ramp loading. During the interval 0 < t < 5, g = 0 and initial
conditions are all 0. So solution remains 0. For large time t, g = 1. A particular solution is
Y = 41 . Adding the homogeneous solution, we should have y = 14 + c1 sin 2t + c2 cos 2t for t
large. We see this is actually the case, the solution is an oscillation around the constant 14 for
large t.
23