Tema 7
Tema 7
1
Z ∞ Z ∞ Z ∞
−st −st u=t− a
L( H (t − a)) = H (t − a)e dt = e dt = e−s(u+a) du =
0 a 0
e−as
= e−as L(1) = ( s ≥ a ).
s
R∞
Example 1.5. Remember the Gamma function Γ( x ) = 0 e−t t x−1 dt for x > 0. Then for
a > −1
Γ ( a + 1)
L (t a ) = , s > 0.
s a +1
R∞ u=st R ∞ a
Indeed L (t a ) = 0 t a e−st dt = 0 usa e−u du s = s a +1 Γ ( a + 1 ) .
1
and the problem is to determine the classes of functions f in such a way that the convergence
is guaranteed.
The next theorem gives us a sufficient condition for existence for the Laplace transforms.
Theorem 1.6 (Existence of L( f )). Let f (t) be piecewise continuous on every finite interval in
t ≥ 0 and satisfy
Proof. It has to be shown that the Laplace integral of f is finite for s > α. Advanced calculus
implies that it is sufficient to show that the integrand is absolutely bounded above by an
integrable function g(t). Take g(t) = Me−(s−α)t . Then g(t) ≥ 0. Furthermore, g is
integrable, because Z ∞
M
g(t) dt = .
0 s − α
Inequality | f (t)| ≤ Meαt implies the absolute value of the Laplace transform integrand
f (t)e−st is estimated by
Obviously (1) gives only a sufficient condition for the existence of a Laplace transform,
and not a necessary one, since Laplace transforms can be found for functions that do not
satisfy condition (1). For example, f (t) = t−1/4 (see Example 1.5).
Furthermore, the property (2) in the previous theorem gives us a criterion to determine
when a function is the Laplace transform of another. For example, polynomial functions
are not Laplace transforms. Instead, the function F (s) = arctan(1/s) for s > 0 could be a
Laplace transform, as we confirm in the Example 2.6.
2
1.2 Properties of the Laplace Operator
Linearity
Proposition 1.7. Let f 1 (t) and f 2 (t) be functions which Laplace transform exists and let c1 and c2
be constant real numbers, then
Translations
Proposition 1.8 (t-shift rule). Let f (t) be a function and g(t) = H (t − a) f (t − a), i.e.
(
f (t − a) for t > a
g(t) =
0 for t < a
with a ≥ 0, then
doing u = t − a
Z ∞ Z ∞ Z ∞
L( g(t)) = f (t − a)e−st dt = f (u)e−s(u+a) du = e−as f (u)e−su du =
a 0 0
− as
=e L( f (t))
Remark. The hypothesis a ≥ 0 in the previous Proposition is necessary, other case, this
result is not true. For a counterexample, given f (t) = t, we obtain
Z ∞
1 1
L [ H (t + 1) f (t + 1)] = (t + 1)e−st dt = L(t) + L(1) = 2
+
0 s s
1+s
= 2 6= es L(t)
s
Example 1.9. For calculating the Laplace transform for rectangular pulse function
1 −
(
1 for a ≤ t < b
f (t) = (with a > 0) | |
0 othercase
a b
3
observe what f (t) = H (t − a) − H (t − b) where H (t) is the Heaviside unit step function.
Then
Rescaling
1 s
Proposition 1.11. If L( f (t)) = F (s) and a > 0, then L( f ( at)) = F .
a a
Proof.
Z ∞ Z ∞ Z ∞
−st at=u −s ua du 1 s
L( f ( at)) = f ( at)e dt = f (u)e = f (u)e− a u du =
0 0 a a 0
1 s
= F
a a
On any interval [0, N ], there are finitely many intervals [ a, b] on each of which f 0 is continu-
ous. Add the above equality across these finitely many intervals [ a, b]. The boundary values
on adjacent intervals match and the integrals add to give
Z N Z N
f 0 (t)e−st dt = f ( N )e− Ns − f (0)e0 + s f (t)e−st dt
0 0
Take the limit across this equality as N → ∞. Then the right side has limit − f (0) +
sL( f (t)), because of the existence of L( f (t)) and limt→∞ f (t)e−st = 0 for large s. There-
fore, the left side has a limit, and by definition L( f 0 (t)) exists and L( f 0 (t)) = − f (0) +
sL( f (t)).
4
Similarly we have:
L( f 00 (t)) = sL( f 0 (t)) − f 0 (0) = s (sL( f (t)) − f (0)) − f 0 (0) =
= s2 L( f (t)) − s f (0) − f 0 (0)
and furthermore L( f 000 (t)) = s2 L( f (t)) − s2 f (0) − s f 0 (0) − f 00 (0). In general,
L f (n) (t) = sn L( f (t)) − sn−1 f (0) − sn−2 f 0 (0) − · · · − f (n−1) (0)
dn
• Hypothesis: L( f (t)) = (−1)n L (tn f (t)). Then
dsn
d n +1
n
d d d
n + 1
L( f (t)) = n
L( f (t)) = [(−1)n L (tn f (t))] =
ds ds ds ds
Z ∞ Z ∞
d n −st
n n
−tn+1 f (t)e−st dt =
= (−1) t f (t)e dt = (−1)
0 ds 0
n +1 n +1
= (−1) L t f (t) .
.
What proves the thesis.
4−s
Example 1.14. Use the above propositions to prove that L (t − 1)e3t =
for s > 3.
( s − 3)2
f (t)
Proposition 1.16. If limt→∞ t exists and L( f (t)) = F (s), then
Z ∞
f (t)
L = F (t) dt.
t s
Proof. Omitted.
5
1.3 Laplace Transform Table
1
Proposition 1.17. L(e at ) = assumed s > a.
s−a
Proof.
∞
Z ∞
e( a−s)t 1 for s > a
L(e at ) = e(a−s)t dt = = s−a
0 a−s ∞ for s ≤ a
t =0
a
Proposition 1.18. L(sin at) = assumed s > 0.
s2 + a2
Proof. First we calculate L(sin t)
d cos t d sin t
L(sin t) = L − = −sL(cos t) + 1 = −sL + 1 = −s2 L(sin t) + 1.
dt dt
1
Hence L(sin t) = . Rescaling (Proposition 1.11)
s2 + 1
1 1 a
L(sin at) = s2
= .
a +1 s2 + a2
a2
s
Proposition 1.19. L(cos at) = assumed s > 0.
s2 + a2
Proof. Analogous.
s
Proposition 1.20. L(cosh at) = assumed s > | a|.
s2 − a2
e at + e−at
Proof. Exercise. Hint: use cosh at = .
2
a
Proposition 1.21. L(sinh at) = 2 assumed s > | a|.
s − a2
Proof. Analogous.
The Dirac Delta is not a function but a concept called distribution. Roughly speaking, it can
be thought of as a function that is very, very high and very, very thin. δ(t − a) is usually
used to move the Dirac delta for a > 0 (see figure 1)..
Often this distribution is defined as the “function” which do
Z ∞
f (t)δ(t) dt = f (0)
−∞
and it can also be see as the limit of families of functions with certain properties, for example
6
`
a
r
n −nx2
• Gaussian functions : δn (t) = e for n = 1, 2, 3, . . .
π
n 1
• Lorentz functions : δn (t) = for n = 1, 2, 3, . . .
π 1 + n2 x 2
• and others.
From
Z t Z t
(
0 if t < a
δ( x − a) dx = lim δn ( x − a) dx = = H (t − a)
−∞ n→∞ −∞ 1 if t ≥ a
d
H (t − a) = δ(t − a)
dt
and so, using t-derivative rule, theorem 1.12, we get the Laplace transform for the Dirac
Delta:
7
(
0 for t = 2
Example 2.2. Functions f 1 (t) = et and f 2 (t) = .
et for t 6= 2
Check
1
L( f 1 (t)) = L( f 2 (t)) = ,
s−1
1
so both functions are the inverse Laplace transform of the same function F (s) = s −1 .
However there are conditions for the uniqueness of the inverse transform as established
next theorem we give without proof.
Theorem 2.3 (Lerch). If f 1 (t) and f 2 (t) are continuous, of exponential order and L( f 1 (t)) =
L( f 2 (t)) for all s > s0 then f 1 (t) = f 2 (t) for all t ≥ 0.
Table 1b shows most important Inverse Laplace transforms, immediate consequence of
table 1a.
Linearity
Let F1 (s) and F2 (s) be functions and let c1 and c2 be constant real numbers, then
Translation 1
If L−1 ( F (s)) = f (t) then L−1 ( F (s − a)) = e at f (t).
Translation 2
If L−1 ( F (s)) = f (t) then L−1 (e−as F (s)) = f (t − a) H (t − a) assuming a ≥ 0.
Rescaling
1 t
If L−1 ( F (s)) = f (t) then L−1 ( F ( as)) = f assuming a ≥ 0.
a a
Derivative rule
If L−1 ( F (s)) = f (t) then L−1 ( F (n) (s)) = (−1)n tn f (t).
Integral rule
R∞ f (t)
If L−1 ( F (s)) = f (t) then L−1
s
F (u) du = .
t
s sin φ + ω cos φ
Example 2.4. The inverse Laplace transform of the function X (s) = is
s2 + ω 2
x (t) = sin(ωt + φ).
8
We can now take the inverse Laplace transform of table 1(b) on page 13:
−1 s −1 ω
x (t) = (sin φ)L + (cos φ)L
s2 + ω 2 s2 + ω 2
= (sin φ)(cos ωt) + (sin ωt)(cos φ)
= sin(ωt + φ).
s+b
Example 2.5. Calculate the inverse Laplace transform of F (s) = .
( s + a )2 + ω 2
s+b s+a 1
F (s) = 2 2
= 2 2
+ (b − a) ,
(s + a) + ω (s + a) + ω ( s + a )2 + ω 2
so
b−a
f (t) = F −1 ( F (s)) = e−at cos ωt + sin ωt .
ω
1 sin t
Example 2.6. The inverse Laplace transform of F (s) = arctan( ) is f (t) = .
s t
−1
The derivative is F 0 (s) = 2 and using derivative rule L−1 ( F 0 (s)) = −t f (t), we
s +1
obtain
1 −1 1 sin t
f (t) = L 2
=
t s +1 t
3 Convolution
Definition 3.1 (Convolution). Let f (t) and g(t) be functions piecewise continuous of ex-
ponential order with f (t) = 0 and g(t) = 0 for t < 0. We call convolution product (or
simply convolution) of f and g to
Z ∞ Z t
( f ∗ g)(t) = f (u) g(t − u) du = f (u) g(t − u) du.
−∞ 0
9
Theorem 3.4. If L−1 ( F (s)) = f (t) and L−1 ( G (s)) = g(t) then
u t
Hence
Z ∞ Z t
F (s) G (s) = f (y) g(t − y)e−st dydt =
t =0 y =0
Z ∞ Z t Z ∞
−st
= e f (y) g(t − y) dydt = e−st ( f ∗ g)(t) dt =
t =0 y =0 0
= L(( f ∗ g)(t))
1
F (s) =
(s + a)(s + b)
The impulse response is simply the inverse Laplace transform of this transfer function f (t) =
L−1 ( F (s)).
To evaluate this inverse transform, we use the convolution property. That is, the inverse
of
1 1 1
F (s) = = ·
(s + a)(s + b) s+a s+b
is
e−at − e−bt
Z t
−1 1 −1 1 − at −bt
f (t) = L ∗L =e ∗e = e−ax e−b(t− x) dx = .
s+a s+b 0 b−a
Exercise. Use method of partial fraction expansion to evaluate the inverse Laplace transform f (s) =
L−1 ( F (s)) being
1 A B
F (s) = = +
(s + a)(s + b) s+a s+b
used in the Example 3.5 above.
10
Example 3.6. Convolution could be used to prove the rule 2.1:
Indeed,
Example 4.1. We use Laplace method for solving the linear ODE
1
(s2 L(y) − sy(0) − y0 (0)) + (sL(y) − y(0)) − 2L(y) =
s2
1
(s2 L(y) − 2s + 1) + (sL(y) − 2) − 2L(y) =
s2
1
(s2 + s − 2)L(y) − 2s − 1 =
s2
1 2s3 + s2 + 1
(s2 + s − 2)L(y) = 2 + 2s + 1 =
s s2
Hence
2s3 + s2 + 1 2s3 + s2 + 1
L(y) = =
s2 ( s2 + s − 2) s2 (s − 1)(s + 2)
Using partial fraction method
11
Example 4.2 (Damped oscillator). Solve by Laplace’s method the initial value problem
From here
s+1 s+1
L( x ) = =
s2 + 2s + 2 ( s + 1)2 + 1
and
−1 s+1 −t −1 s
x (t) = L =e L = e−t cos t
( s + 1)2 + 1 2
s +1
12
4.1 Laplace Transform Table
13
Exercises
Exercise 1
Find the Laplace transform of each of the following functions:
(
3 for 0 < t < 5,
a) f (t) =
0 for t > 5.
b) f (t) = e−2t cos2 3t − 3t2 e3t .
Hint: You can use the equality 2 cos2 a = 1 + cos 2a.
(
cos(t − 2π
3 ) for t > 2π
3 ,
c) f (t) = cos(t − 2π3 ) H ( t − 2π
3 ) = 2π
0 for t < 3 .
Solutions:
2π
3 e5s −1
s +2 se− 3 s
a) F (s) = se5s
. b) F (s) = 2(s1+2) + 1
2 s2 +4s+40 − 6
( s −3)3
. c) F (s) = s2 +1
.
Exercise 2
t
Z
sin u 1 1
Prove that L du = arctan .
0 u s s
Exercise 3
s2 − 2s + 1
If L( f (t)) = compute L( f (2t)). Solution: 4(ss− 2
+1)2
.
(2s + 1)2 (s − 1)
Exercise 4
s2 − a2
Prove L (t cos at) = .
( s2 + a2 )2
Exercise 5
Find L( f (t)), knowing L( f 00 (t)) = arctan 1
and f (0) = 2, f 0 (0) = −1.
s
1
2s+arctan −1
Solution: L( f (t)) = s2
s
.
Exercise 6
Prove that convolution is commutative, i.e. ( f ∗ g)(t) = ( g ∗ f )(t).
Exercise 7
Calculate the inverse Laplace transform of
6s − 4 s+5
a) F (s) = b) F (s) =
s2 − 4s + 20 ( s − 2)3 ( s + 3)
1 e−2s
c) F (s) = 2 2
d) F (s) = .
s (s + 3s − 4) s2 − 4
s
e) F (s) = 2 2
(s − 1) (s + 2s + 5)
14
Solutions:
2e2t 2te2t 7t2 e2t 2e−3t
a) f (t) = e2t (2 sin (4t) + 6 cos (4t)). b) f (t) = 125 − 25 + 10 − 125 .
et e−4t t 3
c) f (t) = 5 − 80 − 4 − 16 . d) f (t) = 12 H (t − 2) sinh(2t − 4).
t et et sin (2t) cos (2t)
e) f (t) = 8 + 16 − e−t 8 + 16 .
Exercise 8
Use the convolution rule to solve the following inverse Laplace transform:
" #
− 1 s − 1 1
a) L b) L
( s2 + a2 )2 s2 ( s + 1)2
Solutions:
t sin ( at)
a) . b) t e−t + 2e−t + t − 2.
2a
Exercise 9
2
Use convolution to prove L sin2 (t) =
.
s ( s2 + 4)
Exercise 10
Solve the following ODEs using Laplace method:
15