Laplace
Laplace
Start with an initial value Using Laplace method, the The solution is
ODE problem (IVP) IVP is transformed to an transferred back to the
(Time domain) algebraic equation that can physical domain
be solved by simple (Time domain)
algebraic manipulation
(s domain)
Laplace transform: basics and definition
• If 𝑓(𝑡) is a function defined for all 𝑡 ≥ 0, its Laplace transform is
Proof
∞
ℒ 𝑎(𝑓 𝑡 + 𝑏(𝑔 𝑡 = 𝑒 −𝑠𝑡 𝑎(𝑓 𝑡 + 𝑏(𝑔 𝑡 𝑑𝑡
0
∞ ∞
=𝑎 𝑒 −𝑠𝑡 𝑓 𝑡 𝑑𝑡 + 𝑏 𝑒 −𝑠𝑡 𝑔 𝑡 𝑑𝑡
0 0
= 𝑎ℒ(𝑓) + 𝑏ℒ(𝑔)
Application on the Linearity of the Laplace transforms
• Find the Laplace of the function 𝑓(𝑡) = sinh 𝑏𝑡
𝑒 𝑏𝑡 + 𝑒 −𝑏𝑡
sinh 𝑏𝑡 =
2
Solution
∞
ℒ sinh 𝑏𝑡 = 𝑒 −𝑠𝑡 sinh 𝑏𝑡 𝑑𝑡
0
∞ ∞ ∞
𝑒 𝑏𝑡 + 𝑒 −𝑏𝑡 1 1
= 𝑒 −𝑠𝑡 𝑑𝑡 = 𝑒 −𝑠𝑡 𝑒 𝑏𝑡 𝑑𝑡 + 𝑒 −𝑠𝑡 𝑒 −𝑏𝑡 𝑑𝑡
0 2 2 0 2 0
∞ ∞
1 1
= 𝑒 −𝑠𝑡 𝑒 𝑏𝑡 𝑑𝑡 + 𝑒 −𝑠𝑡 𝑒 −𝑏𝑡 𝑑𝑡
2 0 2 0
1 1 1 1
= +
2 𝑠−𝑏 2 𝑠+𝑏
𝑠
=
𝑠 2 − 𝑏2
Common Laplace transforms
First shifting theorem (s-Shifting)
Proof
∞
ℒ 𝑒 𝑎𝑡 𝑓(𝑡) = 𝑒 −𝑠𝑡 𝑒 𝑎𝑡 𝑓(𝑡)𝑑𝑡
0
∞ ∞
= 𝑒 −(𝑠−𝑎)𝑡 𝑓(𝑡)𝑑𝑡 = 𝑒 −𝑘𝑡 𝑓(𝑡)𝑑𝑡 , 𝑤ℎ𝑒𝑟𝑒 𝑘 = (𝑠 − 𝑎)
0 0
= 𝐹 𝑘 = 𝐹(𝑠 − 𝑎)
Shifting example
Knowing the following
Solution:
With problems like this, usually try to simplify using completing the squares and partial fractions
𝑠 2 + 2𝑠 + 401 = 𝑠 2 + 2𝑠 + 1 − 1 + 401
= 𝑠 + 1 𝑠 + 1 + 400 = 𝑠 + 1 2 + 202
3𝑠 − 137 s+1 20
ℒ −1 = 3ℒ −1 − 7ℒ −1
𝑠 2 + 2𝑠 + 401 𝑠 + 1 2 + 202 𝑠 + 1 2 + 202
= 3 𝑒 −𝑡 cos 20𝑡 − 7 𝑒 −𝑡 sin 20𝑡
Transforms of Derivatives and Integrals
Laplace Transform of Derivatives
First order
Second order
Proofs
1st order
2nd order
𝐵𝑢𝑡
ℒ 𝑓 ′′ 𝑡 = 𝑠 2 ℒ 𝑓 − 𝑠𝑓 0 − 𝑓′(0) we used the derivative rule
= 𝑠2ℒ 𝑓
𝑠
2𝜔 2 2 − 𝜔2 ℒ 𝑓 = 𝑠 2 ℒ 𝑓 solve using algebraic manipulation
𝑠 +𝜔
2𝜔𝑠
ℒ 𝑓 = 2
𝑠 + 𝜔2 2
Transforms of Derivatives and Integrals
Laplace Transform of Integrals
Proof
𝑡
𝐿𝑒𝑡 𝑔(𝑡) = 𝑓 𝜏 𝑑𝜏 𝑓 𝑡 = 𝑔′(𝑡)
0
ℒ 𝑓(𝑡)
ℒ 𝑔(𝑡) =
𝑆
𝑡
1
ℒ 𝑓 𝜏 𝑑𝜏 = 𝐹(𝑠)
0 𝑆
Solving Initial Value Problems Using Laplace
Solving IVP using Laplace Transforms
General Framework
• To solve an initial value problem
Solution
• Step 1: Transform it to an algebraic equation
𝑎𝑡 𝑥 = 0 𝐵 = −1
𝑎𝑡 𝑥 = 1 𝐶−𝐷 =1 𝐷 = −1
𝐶=0
𝑎𝑡 𝑥 = −1 −𝐶 − 𝐷 = 1
𝑎𝑡 𝑥 = 2 1 = 6𝐴 − 3 + 4 𝐴=0
1 𝐴 𝐵 𝐶𝑥 − 𝐷
= + +
𝑥 2 (𝑥 2 − 1) 𝑥 𝑥 2 (𝑥 2 − 1)
−1 1
= +
𝑥2 (𝑥 2 − 1)
Partial fraction expansion
Example: Find the partial fraction expansion of
1
𝑥 2 (3𝑥 + 1)
1 𝐴 𝐵 𝐶
General expansion = + +
𝑥 2 (3𝑥 + 1) 𝑥 𝑥 2 (3𝑥 + 1)
𝑎𝑡 𝑥 = 0 𝐵=1
𝑎𝑡 𝑥 = −1/3 𝐶=9
𝑎𝑡 𝑥 = 1 1 = 4𝐴 + 4 + 9 𝐴 = −3
1 −3 1 9
= + +
𝑥 2 (3𝑥 + 1) 𝑥 𝑥 2 (3𝑥 + 1)
Partial fraction expansion
Example: Find the partial fraction expansion of
1
𝑥 2 (3𝑥 + 1)2
1 𝐴 𝐵 𝐶 𝐷
General expansion = + + +
𝑥 2 (3𝑥 + 1)2 𝑥 𝑥 2 (3𝑥 + 1) (3𝑥 + 1)2
𝑎𝑡 𝑥 = 0 𝐵=1
𝑎𝑡 𝑥 = −1/3 𝐷=9
𝐶 = 18, 𝐴 = −6
1 6 1 18 9
=− + 2− +
𝑥 2 (3𝑥 + 1)2 𝑥 𝑥 (3𝑥 + 1) (3𝑥 + 1)2
Shifted Initial Value Problems
Solving Shifted IVP using Laplace
Initial value problems whose initial conditions are defined at a point
other than the origin (t=0) are referred to as shifted IVP
Example: 𝑦 ′′ + 𝑎𝑦 ′ + 𝑏𝑦 = 𝑅(𝑡) 𝑦 𝑐 = 𝐷, 𝑦 ′ 𝑐 = 𝐹
Solution:
• 𝐴𝑝𝑝𝑙𝑦 𝑡ℎ𝑒 𝐿𝑎𝑝𝑙𝑎𝑐𝑒 𝑡𝑟𝑎𝑛𝑠𝑓𝑜𝑟𝑚 𝑡𝑜 𝑏𝑜𝑡ℎ 𝑒𝑞𝑢𝑎𝑡𝑜𝑖𝑛𝑠
𝑠𝑌1 − 𝑦1 (0) = 2𝑌1 + 𝑌2
𝑠𝑌2 − 𝑦2 0 = 𝑌1 + 2𝑌2
• 𝑆𝑜𝑙𝑣𝑒 𝑓𝑜𝑟𝑌1 and 𝑌2
𝑠 − 2 𝑌1 − 𝑌2 = 𝑦1 (0)
𝑠 − 2 −1 𝑌1 1
=
−𝑌1 + 𝑠 − 2 𝑌2 = 𝑦2 (0) −1 𝑠 − 2 𝑌2 2
Systems of 1st order ODEs
𝑠 − 2 −1 𝑌1 1
=
−1 𝑠 − 2 𝑌2 2
𝑌1 −1
𝑠−2 −1 1
=
𝑌2 −1 𝑠−2 2
3 1 1 3
𝑌1 = − 𝑎𝑛𝑑 𝑌2 = +
2(𝑠 − 3) 2 𝑠 − 1 2(𝑠 − 1) 2(𝑠 − 3)
𝑦1 (𝑡) 3 1 3𝑡 1 −1 𝑡
= 𝑒 + 𝑒
𝑦2 (𝑡) 2 1 2 1
∞
∞
𝑒 −𝑠𝑡 𝑒 −𝑎𝑠
Proof: 𝑒 −𝑠𝑡 . 1 𝑑𝑡 = − =
𝑎 𝑠 𝑠
𝑎
Example: Unit step function
Write the following 𝑓(𝑡) using unit step function and find its Laplace
transform
𝑓(𝑡) 6
2, 𝑖𝑓 0 < 𝑡 < 1
𝑓 𝑡 = 6, 𝑖𝑓 1 < 𝑡 < 5 2
1, 𝑖𝑓 5 < 𝑡 < 8 1
𝑡
1 5 8
Solution:
𝑓 𝑡 = 2𝑢 𝑡 − 2𝑢 𝑡 − 1 + 6𝑢 𝑡 − 1 − 6𝑢 𝑡 − 5 + 1𝑢 𝑡 − 5 − 1𝑢(𝑡 − 8)
𝑓 𝑡 = 2 𝑢 𝑡 − 𝑢 𝑡 − 1 + 6 𝑢 𝑡 − 1 − 𝑢 𝑡 − 5 + 𝑢 𝑡 − 5 − 𝑢(𝑡 − 8)
Proof:
∞ ∞
= 𝑒 −𝑠𝑡 𝑓 𝑡 − 𝑎 𝑑𝑡 = 𝑒 −𝑠𝑡 𝑓 𝑡 − 𝑎 𝑢(𝑡 − 𝑎)𝑑𝑡
𝑎 0
Example: Unit step function and second shifting theorem
Write the following function using unit step functions and find its Laplace transform
Solution:
1 2 1 𝜋 𝜋
𝑓 𝑡 = 2𝑢 𝑡 − 2𝑢 𝑡 − 1 + 𝑡 𝑢 𝑡 − 1 − 𝑡 2𝑢 𝑡 − + cos 𝑡 𝑢 𝑡 −
2 2 2 2
To find the Laplace transform, cast all terms into the form 𝑓(𝑥 − 𝑎)𝑢(𝑥 − 𝑎)
1 2 1 1
𝑡 𝑢 𝑡−1 = 𝑡−1 2+ 𝑡−1 + 𝑢 𝑡−1
2 2 2
1 2 𝜋 1 𝜋 2 𝜋 𝜋 𝜋2 𝜋
𝑡 𝑢 𝑡− = 𝑡− + 𝑡− + 𝑢 𝑡−
2 2 2 2 2 2 8 2
𝜋 𝜋 𝜋
cos 𝑡 𝑢 𝑡 − = −sin 𝑡 − 𝑢 𝑡 −
2 2 2
Apply the Laplace transform to each term inside the three bracketed components
2 𝑒 −𝑠
• ℒ 2𝑢 𝑡 − 2𝑢 𝑡 − 1 = −𝑠
2 𝑠
Unit step function Example
1 2 1 2 1 𝑒 −𝑠 𝑒 −𝑠 𝑒 −𝑠
• ℒ 𝑡 𝑢 𝑡−1 =ℒ 𝑡−1 + 𝑡−1 + 𝑢(𝑡 − 1) = + +
2 2 2 𝑠3 𝑠2 2𝑠
1 2 𝜋 1 𝜋 2 𝜋 𝜋 𝜋2 𝜋
• ℒ 𝑡 𝑢(𝑡 − ) =ℒ 𝑡− + 𝑡− + 𝑢(𝑡 − )
2 2 2 2 2 2 8 2
𝜋 𝜋 𝜋
− 𝑠 − 𝑠 2 − 𝑠
𝑒 2 𝜋𝑒 2 𝜋 𝑒 2
= + +
𝑠3 2 𝑠2 8 𝑠
𝜋
− 𝑠
𝜋 𝜋 𝑒 2
• ℒ −sin 𝑡 − 𝑢(𝑡 − ) = −
2 2 𝑠 2 +1
1 1 𝜋 𝜋
ℒ 𝑓 𝑡 = ℒ 2𝑢 𝑡 − 2𝑢 𝑡 − 1 + 𝑡 2 𝑢 𝑡 − 1 − 𝑡 2 𝑢 𝑡 − + cos 𝑡 𝑢 𝑡 −
2 2 2 2
𝜋 𝜋 𝜋 𝜋
−2𝑠 −2𝑠 −2𝑠 −2𝑠
2 2𝑒 −𝑠 𝑒 −𝑠 𝑒 −𝑠 𝑒 −𝑠 𝑒 𝜋𝑒 2
𝜋 𝑒 𝑒
= − + + + − − − − 2
𝑠 𝑠 𝑠3 𝑠2 2𝑠 𝑠3 2 𝑠2 8 𝑠 𝑠 +1
Example: (second shifting theorem)
Find the Laplace inverse of the following
Solution:
1 1
So ℒ −1 𝐹 𝑠 = sin 𝜋 𝑡 − 1 𝑢 𝑡 − 1 + sin 𝜋 𝑡 − 2 𝑢 𝑡 − 2 + 𝑒 −2 𝑡−3
𝑡 − 3 𝑢(𝑡 − 3)
𝜋 𝜋
Delta Function
Delta function is used to represent short impulses. For example: to represent
mechanical impacts and sudden loadings (mechanical or electrical)
Impulse function definition • Dirac delta function is the impulse function
when 𝑘 → 0
Solution:
Apply the Laplace transform to the equation
𝑠 2 𝑌 + 3𝑠𝑌 + 2𝑌 = 𝑒 −𝑠
ℒ sin 𝜔𝑡 1
We know that = 2
𝜔 𝑠 + 𝜔2
sin 𝜔𝑡
So 𝑓 𝑡 = 𝑔 𝑡 =
𝜔
𝑡 𝑡
1 sin 𝜔𝜏 sin 𝜔(𝑡 − 𝜏)
ℒ −1 = 𝑓 𝜏 𝑔 𝑡 − 𝜏 𝑑𝜏 = 𝑑𝜏
𝑠2 + 𝜔2 2
0 0 𝜔 𝜔
𝑡
1 1 sin 𝜔𝑡
= 2 sin 𝜔𝜏 sin 𝜔(𝑡 − 𝜏) 𝑑𝜏 = −𝑡 cos 𝜔𝑡 +
𝜔 0 2 𝜔
Differentiation and Integration of Laplace transforms
• These relations are useful to solve IVP with variable coefficients
• Also they are useful in determining the Laplace or its inverse for complex cases