374 Lecture 26
Chapter 9: Laplace Transforms
374 Lecture 26
Chapter 9: Laplace Transforms
Chapter 9: Laplace Transforms
9.1: An Introduction to the Laplace Transform
9.1: An Introduction to the Laplace Transform
Definition of Laplace Transform
Definition of Laplace Transform
Def: E is the class of continuous or piecewise continuous functions f(t), defined for at least t > 0, such that |f(t)| < for some A, B. If f ϵ E , then we say f is of exponential type.Note that not all functions are of exponential type. For example, f(t) = is not of exponential type. Def: If f ϵE, then the Laplace transform of f(t) is defined to be the function F(s) f(t)t, defined at least for all sufficiently large t. We will write F(s) = L{f} = L{f(t)}.
Bt
A
()
2
t
:=
∞
∫
0
-st
e
Laplace Transforms of Common Functions
Laplace Transforms of Common Functions
FrameBox[GridBox[Table[{{"f(t)","F(s)"},{,LaplaceTransform[,t,s]},{t*,LaplaceTransform[t*,t,s]},{*,LaplaceTransform[*,t,s]},{*,LaplaceTransform[*,t,s]},{Sin[at],LaplaceTransform[Sin[at],t,s]},{Cos[at],LaplaceTransform[Cos[at],t,s]},{t*Sin[at],LaplaceTransform[tSin[at],t,s]},{t*Cos[at],LaplaceTransform[tCos[at],t,s]},{Sinh[at],LaplaceTransform[Sinh[at],t,s]},{Cosh[at],LaplaceTransform[Cosh[at],t,s]},{1,LaplaceTransform[1,t,s]},{t,LaplaceTransform[t,t,s]},{,LaplaceTransform[,t,s]},{,LaplaceTransform[,t,s]}}],RowLinesTrue,ColumnLinesTrue]]//DisplayForm
at
at
at
at
2
t
at
2
t
at
n
t
at
n
t
at
2
t
2
t
n
t
n
t
Out[]//DisplayForm=
Properties of Laplace Transforms
Properties of Laplace Transforms
The Laplace transform is linear.
The Laplace transform is linear.
Let F(s) and G(s) be the Laplace transforms of f(t) and g(t), respectively. Then L{α f(t)+β g(t)} = α F(s) + β G(s).
Try applying this idea, using the table above, with h(t) =. Check with Mathematica!
Try applying this idea, using the table above, with h(t) =
4t-3Cos(t)
t
LaplaceTransform[4t-3Cos[t],t,s]
t
Out[]=
4
2
(-1+s)
3s
1+
2
s
Laplace Transforms of Derivatives
Laplace Transforms of Derivatives
One of the main reasons that Laplace transforms are used is the fact that they turn differentiation into multiplication!
Theorem 9.1: Let f(t) and f '(t) be inE. Then L{f ' (t)} =.
Proof: Use integration by parts.
Corollary 1: If f and its first n derivatives belong to E, then
Check these results hold in the following examples!
Theorem 9.1: Let f(t) and f '(t) be inE. Then L{f ' (t)} =
-f(0)+sL{f(t)}=-f(0)+sF(s)
Proof: Use integration by parts.
Corollary 1: If f and its first n derivatives belong to E, then
L{(t)}=-[f(0)+...+s(0)+(0)]+L{f(t)}=-[f(0)+...+s(0)+(0)]+F(s)
(n)
f
n-1
s
(n-2)
f
(n-1)
f
n
s
n-1
s
(n-2)
f
(n-1)
f
n
s
Check these results hold in the following examples!
LaplaceTransform[f'[t],t,s]/.LaplaceTransform[f[t],t,s]F[s]
Out[]=
-f[0]+sF[s]
LaplaceTransform[f''[t],t,s]/.LaplaceTransform[f[t],t,s]F[s]
Out[]=
-sf[0]+F[s]-[0]
2
s
′
f
LaplaceTransform[f'''[t],t,s]/.LaplaceTransform[f[t],t,s]F[s]
Out[]=
-f[0]+F[s]-s[0]-[0]
2
s
3
s
′
f
′′
f
Note that if f(t) is such that f and all of its derivatives up to order n-1 are zero at t = 0, then Corollary 1 implies .
Corollary 2: Let T(D) be a differential operator with constant coefficients: T(D) =++...+. If f and all of its derivatives up to order n belong to E and f(0) = f '(0) = ... = (0)=0, then L{T(D)f} = T(s) L{f} where p(s) = ++...+.
Try this out with the following command!
L{(t)}=L{f(t)}=F(s)
(n)
f
n
s
n
s
Corollary 2: Let T(D) be a differential operator with constant coefficients: T(D) =
a
0
n
D
a
1
n-1
D
a
n
(n-1)
f
a
0
n
s
a
1
n-1
s
a
n
Try this out with the following command!
Simplify[LaplaceTransform[5f''''[t]+3f''[t]-f'[t]+4f[t],t,s]/.{LaplaceTransform[f[t],t,s]F[s],f[0]0,f'[0]0,f''[0]0,f'''[0]0}]
Out[]=
(4-s+3+5)F[s]
2
s
4
s
The Convolution Theorem
The Convolution Theorem
Def: Let f and g be functions. The convolution of f and g, denoted f*g, is the function defined by f[t-u]g[u]u.Theorem 9.2 (Convolution Theorem): Let f(t) and g (t) be inE. Then L{f*g} = L{f} L{g}=F(s) G(s).Here is a way to define the convolution of two functions f and g via Mathematica. (Note that there is also a built-in function, Convolve.)
t
∫
0
In[]:=
Convolution[f_,g_]:=Integrate[f[t-u]g[u],{u,0,t}]
Convolution[f,g]
Out[]=
t
∫
0
Convolution[Sin,Cos]
Out[]=
1
2
Simplify[%]
Out[]=
1
2
Plot[Convolution[Sin,Cos],{t,0,2π}]
Out[]=
Test out the Convolution Theorem with Mathematica!
LaplaceTransform[Convolution[Sin,Cos],t,s]
Out[]=
s
2
(1+)
2
s
Simplify[%]
Out[]=
s
2
(1+)
2
s
LaplaceTransform[Sin[t],t,s]
Out[]=
1
1+
2
s
LaplaceTransform[Cos[t],t,s]
Out[]=
s
1+
2
s
Inverse Laplace Transforms
Inverse Laplace Transforms
Definition
Definition
Def: Given function F(s), the inverse Laplace transform of F(s) is the function f(t) such that L{f(t)} = F(s). Notation: {F(s)} = f(t).Since the Laplace transform is linear, it follows that the inverse Laplace transform is linear, i.e. {F(s)+G(s)}={F(s)}+{G(s)}.To find inverse Laplace transforms, one can use tables such as those above. For example, =. Mathematica has a built in command to find inverse Laplace transforms. Here is an example:
-1
L
-1
L
c
1
c
2
c
1
-1
L
-1
L
-1
L
1
s-4
4t
Partial Fractions and Inverse Laplace Transforms
Partial Fractions and Inverse Laplace Transforms
The Convolution Theorem and Inverse Laplace Transforms
The Convolution Theorem and Inverse Laplace Transforms
Heaviside's Formula
Heaviside's Formula
A useful way to find inverse Laplace transforms is via the following theorem, which is really just a special case of finding inverse Laplace Transforms via partial fractions.
9.2: Solving ODE with Laplace Transforms
9.2: Solving ODE with Laplace Transforms
Example 1: Finding a particular solution to a differential equation.
Example 1: Finding a particular solution to a differential equation.
To simplify things, use the ReplaceAll command (/.). Assume u[0]=0, u'[0]=0.
Check that the solution solves the ODE!
Try solving this problem using the Convolution Theorem!
Example 2: Solving an IVP with initial data specified at t = 0
Example 2: Solving an IVP with initial data specified at t = 0
Example 2 Solution
Example 2 Solution
Check our solution!