ODE: Solving IVPs with Laplace Transforms

When an engineering system is subjected to a sudden shock, a pulsed input, or any non-smooth forcing function, classical methods for solving differential equations become cumbersome. The Laplace transform provides a powerful algebraic alternative, converting differential equations in the time domain into simpler algebraic equations in the complex frequency domain (the s-domain). This method is exceptionally efficient for solving initial value problems (IVPs), especially those with discontinuous or impulsive forcing terms common in control systems, circuit analysis, and vibrations. By transforming the entire IVP, you handle the initial conditions seamlessly from the start, streamlining the solution process for both single equations and coupled systems.

The Core Laplace Transform Method for IVPs

The Laplace transform is defined as an integral that converts a function of time, $f(t)$, into a function of a complex variable, $s$. The key property for solving ODEs is its action on derivatives. For a function $y(t)$ with Laplace transform $Y(s)=\mathcal{L}\{y(t)\}$, the transform of its derivative is: $$\mathcal{L}\{y^{\prime }(t)\}=sY(s)-y(0)$$ Similarly, for a second derivative: $$\mathcal{L}\{y^{\prime \prime }(t)\}=s^{2}Y(s)-sy(0)-y^{\prime }(0)$$ This property systematically incorporates the initial conditions $y(0)$, $y^{\prime }(0)$, etc., into the algebraic equation.

The step-by-step solution method for an IVP is a structured four-step algorithm:

1. Apply the Laplace Transform to both sides of the entire differential equation. Use the derivative formulas to transform terms like $y^{\prime \prime }$ and $y^{\prime }$, which introduces the initial conditions and transforms the ODE into an algebraic equation for $Y(s)$.
2. Solve the Algebraic Equation for the unknown $Y(s)$. This typically involves rearranging terms and using algebra to isolate $Y(s)$ on one side of the equation.
3. Apply the Inverse Laplace Transform to find $y(t)={\mathcal{L}}^{-1}\{Y(s)\}$. This is often the most challenging step, requiring you to manipulate $Y(s)$ into a form recognizable from a table of Laplace transform pairs. Techniques like partial fraction decomposition are essential here.
4. Simplify and Interpret the resulting time-domain function $y(t)$, which is the unique solution to the original IVP.

Consider the IVP: $y^{\prime \prime }+3y^{\prime }+2y=0$, with $y(0)=1$, $y^{\prime }(0)=0$.

• Step 1: Transform: $(s^{2}Y(s)-s\cdot 1-0)+3(sY(s)-1)+2Y(s)=0$.
• Step 2: Solve for $Y(s)$: $s^{2}Y(s)-s+3sY(s)-3+2Y(s)=0\Rightarrow Y(s)(s^2+3s+2)=s+3$. Thus, $Y(s)=\frac{s+3}{s^2+3s+2}$.
• Step 3: Factor and use partial fractions: $Y(s)=\frac{s+3}{(s+1)(s+2)}=\frac{2}{s+1}-\frac{1}{s+2}$. Looking up inverse transforms: $y(t)=2e^{-t}-e^{-2t}$.
• Step 4: This is the solution, easily verified against the initial conditions.

Handling Various Forcing Functions and Discontinuities

The true power of the Laplace method shines when the non-homogeneous term, or forcing function $f(t)$, is not a simple exponential or sinusoid. Classical methods like undetermined coefficients require guessing the solution's form, which fails for discontinuous inputs. The Laplace transform handles them with ease.

Common forcing functions and their transforms include:

• Constants and Polynomials: $\mathcal{L}\{t^n\}=\frac{n!}{s^{n+1}}$
• Exponentials: $\mathcal{L}\{e^{at}\}=\frac{1}{s-a}$
• Sinusoids: $\mathcal{L}\{\sin (\omega t)\}=\frac{\omega }{s^2+{\omega }^2}$
• The Heaviside step function $u_c(t)$, which is 0 for $t<c$ and 1 for $t\ge c$: $\mathcal{L}\{u_c(t)\}=\frac{e^{-cs}}{s}$.
• The Dirac delta function $\delta (t-c)$, modeling an impulse: $\mathcal{L}\{\delta (t-c)\}=e^{-cs}$.

The second translation theorem (or time-shift property) is critical for discontinuous inputs: $\mathcal{L}\{f(t-c)u_c(t)\}=e^{-cs}F(s)$, where $F(s)=\mathcal{L}\{f(t)\}$. This allows you to transform functions that "turn on" at time $t=c$. For example, to solve $y^{\prime \prime }+y=u_3(t)$ with $y(0)=y^{\prime }(0)=0$, the forcing function $f(t)=u_3(t)$ transforms to $\frac{e^{-3s}}{s}$. After solving for $Y(s)$, the term $e^{-3s}\cdot \frac{1}{s(s^2+1)}$ requires an inverse transform using the time-shift property, yielding a solution $y(t)=[1-\cos (t-3)]u_3(t)$, which is zero until $t=3$.

Solving Systems of ODEs Simultaneously

The Laplace transform method extends elegantly to systems of coupled linear ODEs with constant coefficients. The procedure is a direct generalization of the single-equation case: you transform every equation in the system simultaneously. This yields a system of algebraic equations in the Laplace transforms of the unknown functions (e.g., $X(s)$ and $Y(s)$). You then solve this algebraic system using linear algebra techniques like substitution or Cramer's rule, find the inverse transforms of each result, and obtain the solutions $x(t)$ and $y(t)$.

Consider the coupled system describing two interconnected masses or circuits: $$\begin{array}{rlr}{x}^{\prime }& =-2x+y,& x(0)=1\\ y^{\prime }& =x-2y,& y(0)=0\end{array}$$ Applying the Laplace transform: $$\begin{array}{rl}sX(s)-1& =-2X(s)+Y(s)\\ sY(s)-0& =X(s)-2Y(s)\end{array}$$ Rearranging gives the linear system: $$\begin{array}{rl}(s+2)X(s)-Y(s)& =1\\ -X(s)+(s+2)Y(s)& =0\end{array}$$ Solving for $X(s)$ and $Y(s)$ (e.g., using Cramer's rule) yields: $$X(s)=\frac{s+2}{(s+1)(s+3)},Y(s)=\frac{1}{(s+1)(s+3)}$$ Partial fractions and inverse transforms then give the solution: $x(t)=\frac{1}{2}{e}^{-t}+\frac{1}{2}{e}^{-3t}$ and $y(t)=\frac{1}{2}{e}^{-t}-\frac{1}{2}{e}^{-3t}$. This simultaneous solution elegantly handles the coupling without first needing to decouple the system.

Common Pitfalls

1. Incorrectly Transforming Derivatives: The most frequent error is misapplying the derivative formula, especially forgetting to subtract the initial condition. Remember: $\mathcal{L}\{y^{\prime }\}=sY(s)-y(0)$, not simply $sY(s)$. Always write out the formula explicitly.
2. Algebraic Errors in Solving for $Y(s)$: The intermediate algebra between Step 1 and Step 2 is a common failure point. Move slowly, distribute terms carefully, and group all $Y(s)$ terms together before factoring. A single sign error here will propagate through the entire solution.
3. Misapplying the Inverse Transform: Attempting to find ${\mathcal{L}}^{-1}\{F(s)G(s)\}$ as ${\mathcal{L}}^{-1}\{F(s)\}\cdot {\mathcal{L}}^{-1}\{G(s)\}$ is incorrect—this ignores the convolution property. Always break $Y(s)$ into a sum of simpler terms via partial fractions before using the transform table. Also, when using the time-shift theorem, ensure the function is properly expressed as $e^{-cs}F(s)$ where $F(s)$ is the transform of the "unshifted" function.
4. Neglecting the Domain of the Solution for Discontinuous Inputs: When your solution involves a Heaviside function $u_c(t)$, your final answer must include it to explicitly show the solution is zero for $t<c$. Writing $y(t)=1-\cos (t-3)$ for the earlier example would be incorrect; it must be $y(t)=[1-\cos (t-3)]u_3(t)$.

Summary

• The Laplace transform converts an initial value problem from the time domain into an algebraic problem in the s-domain, where initial conditions are incorporated automatically from the outset.
• The method is superior to classical approaches for equations with discontinuous or impulsive forcing functions (modeled by Heaviside and Delta functions), thanks to the second translation theorem.
• The solution process is algorithmic: Transform the ODE, Solve the algebraic equation for $Y(s)$, Invert the transform to find $y(t)$, and Simplify. Partial fraction decomposition is a crucial tool for the inversion step.
• The technique generalizes directly to solving systems of ODEs by transforming all equations simultaneously, solving the resulting algebraic system for multiple transforms, and inverting each one.
• Common mistakes to avoid include errors in the derivative transform formula, algebraic slips when isolating $Y(s)$, and misapplication of inverse transform rules, particularly with convolutions and time-shifts.