Understanding Laplace Transforms for Solving Differential Equations of Dynamic Systems
3Blue1Brown
Summary:
The video introduces the Laplace transform as a powerful tool for analyzing dynamic systems, using the example of a forced harmonic oscillator, where initial irregular motion transitions to a consistent rhythm.
It recaps key Laplace transform properties: how exponential functions transform into simple fractions with poles in the s-plane, and the transform's linearity, where poles in the s-plane directly reveal the exponential components of a function.
A crucial third property is explained: Laplace transforms convert differentiation in the time domain into multiplication by 's' in the s-domain, simplifying the solution of differential equations by incorporating initial conditions.
0].">
The video demonstrates applying this property to the forced oscillator equation, transforming it into an algebraic problem whose solution, represented by poles in the s-plane, reveals the system's dynamic behavior. These poles show both the decaying natural oscillation (transient response) and the sustained oscillation from the external force (steady-state response).
The process of inverting the transform back to the time domain using partial fraction decomposition is shown, providing an exact analytical solution and highlighting how the amplitude of steady-state oscillation relates to the system's resonant frequency and external forcing frequency, which is crucial for understanding phenomena like resonance in bridges.
Opening Puzzle [0:00]
The video begins by presenting a simulation of a mass on a spring influenced by an oscillating external force. Initially, the system exhibits irregular behavior before settling into a consistent rhythm. The central questions posed are:
- How can this irregular startup trajectory be mathematically analyzed?
- How long does it take for the system to reach its steady rhythm?
- What determines the amplitude of the oscillations once it stabilizes?
The Laplace transform is introduced as a powerful tool to study such dynamic systems and differential equations.
Key Properties of a Laplace Transform [1:06]
This section recaps fundamental concepts from previous chapters on Laplace transforms:
- S-plane and
e^st functions: The s-plane is a complex plane representing values of 's' in e^st.
- Larger imaginary values of 's' correspond to more oscillation.
- Negative real parts of 's' indicate decay, while positive real parts indicate growth.
- Many natural functions can be decomposed into exponential pieces.
- Laplace Transform Definition: The transform converts a time-domain function
f(t) into an s-domain function F(s).
- Poles in the
F(s) plot correspond to the exponential pieces hidden within f(t).
- Property 1: Exponential to Fraction:
L{e^at} = 1/(s-a).
- This means an exponential function in the time domain corresponds to a pole in the s-plane at value
a.
- Property 2: Linearity:
L{a*f(t) + b*g(t)} = a*F(s) + b*G(s).
- A scaled sum of functions transforms into the same scaled sum of their individual transforms.
- This implies that a function expressible as a sum of exponentials will transform into a sum of fractions, with each fraction having a pole corresponding to an exponential piece.
Qualitative Analysis with Laplace Transforms [3:29]
The properties of Laplace transforms enable qualitative analysis of a system's dynamics:
- Imaginary poles: Indicate oscillation in the system.
- Negative real parts of poles: Suggest a tendency for the system to decay towards zero (stability).
- Positive real parts of poles: Point to instability, where the system tends to explode away from zero (growth).
The Laplace Transforms of a Derivative [4:29]
A third crucial property allows for the transformation of differential equations into algebraic ones:
- Property 3:
L{f'(t)} = sF(s) - f[<a href="https://youtube.com/watch?v=FE-hM1kRK4Y&t=0">0</a>].
- This transforms differentiation in the time domain (
d/dt) into multiplication by s in the s-domain.
- The
f[<a href="https://youtube.com/watch?v=FE-hM1kRK4Y&t=0">0</a>] term represents the initial condition, which is automatically incorporated into the transformed equation.
The Forced Oscillator [6:06]
The video demonstrates applying Laplace transforms to solve the forced harmonic oscillator problem:
- System Description: A mass on a spring with damping is subjected to an external oscillating force
F0 cos(ωt).
- The differential equation is
mx''(t) + μx'(t) + kx(t) = F0 cos(ωt).
- General Strategy:
- Take the Laplace transform of all terms in the differential equation.
- Solve the resulting algebraic equation to find the transformed solution
X(s).
- Invert the Laplace transform to recover the solution
x(t) in the time domain.
- Transforming the Equation:
L{x(t)} = X(s)
L{x'(t)} = sX(s) - x0 (where x0 is initial position)
L{x''(t)} = s^2 X(s) - sx0 - v0 (where v0 is initial velocity)
- For simplicity, assume
x0 = 0 and v0 = 0.
- The left side
m x''(t) + μx'(t) + kx(t) transforms into X(s)(ms^2 + μs + k).
- The right side
F0 cos(ωt) transforms into F0s / (s^2 + ω^2).
- Combining these yields the transformed solution:
L{x(t)} = X(s) = (F0s) / ((s^2 + ω^2)(ms^2 + μs + k)).
Intuition from the Transformed Solution [11:59]
Analyzing the poles of the transformed solution X(s) provides insight into the system's behavior:
- Poles from
(ms^2 + μs + k) = 0: These roots typically have negative real parts and imaginary parts (assuming damping and oscillation).
- They represent the natural resonant frequency of the unforced oscillator, indicating a decaying oscillation.
- Poles from
(s^2 + ω^2) = 0 (i.e., s = ±iω): These are pure imaginary poles.
- They correspond to the frequency of the external cosine force, representing a pure, undamped oscillation.
- Combined Dynamics: The total solution
x(t) is a sum of components corresponding to these poles.
- The decaying oscillation from the system's natural frequency (transient response) combined with the pure oscillation from the external force (steady-state response).
- This explains the initial "wibbling" period in the simulation, where both frequencies interact, eventually settling into the steady-state rhythm driven by the external force as the natural oscillation decays.
Inverting to Find a Final Answer [15:15]
To obtain the exact time-domain solution x(t), the inverse Laplace transform is applied:
- Partial Fraction Decomposition: The transformed solution
X(s) is broken down into a sum of simpler fractions, each corresponding to one of its poles (roots of the denominator).
- Example:
X(s) = c1/(s-r1) + c2/(s-r2) + c3/(s-r3) + c4/(s-r4).
- Solving for the constants (c_n) is required.
- Inverse Transform of Each Term: Each simple fraction
c_n/(s-r_n) inverts back into an exponential term c_n * e^(r_n * t).
- Steady-State Amplitude: Focusing on the poles
±iω (from the external force), their corresponding terms combine to form the steady-state cosine rhythm: (F0 / (m(k/m - ω^2) + 2μiω)) * 2cos(ωt).
- The homework exercise is to analyze how the amplitude of this final expression depends on the difference between the resonant frequency of the spring and the frequency of the external force (
k/m - ω^2), which relates to phenomena like resonance in bridges.
Explaining the Derivative Property [17:40]
The video provides three ways to understand why L{f'(t)} = sF(s) - f[<a href="https://youtube.com/watch?v=FE-hM1kRK4Y&t=0">0</a>] is true:
- 1. Elementary but Limited (Testing an Example) [18:10]:
- Demonstrates that for
f(t) = e^at, applying the derivative and then the Laplace transform yields a/(s-a).
- Algebraically manipulating
sF(s) - f[<a href="https://youtube.com/watch?v=FE-hM1kRK4Y&t=0">0</a>] for this function also results in a/(s-a), proving consistency for this specific case.
- This approach provides intuition but is not a general proof; however, it can be extended to combinations of exponentials due to linearity.
- 2. General but Opaque (Integration by Parts) [20:25]:
- The formal derivation uses the definition of the Laplace transform as an integral (
L{f'(t)} = ∫[0,∞] f'(t)e^(-st) dt) and applies integration by parts.
- This method rigorously proves the property but can obscure the underlying intuition.
- 3. Favorite but Theoretical (Reinventing the Transform) [21:13]:
- This explanation, to be covered in the next chapter, involves understanding the inverse Laplace transform and contour integrals.
- The idea is to logically "reinvent" the Laplace transform and its inverse from the goal of converting derivatives into multiplication, showing the property is inherent to the transform's design. This connects to Fourier series and transforms.