Solution of Differential Equations Using Laplace Transforms

The Laplace transform is a powerful mathematical tool that can simplify the process of solving linear ordinary differential equations (ODEs), especially those with constant coefficients and specific initial conditions. The core idea is to transform the differential equation in the time domain (usually represented by 't') into an algebraic equation in the frequency domain (usually represented by 's'). Solving this algebraic equation is often much easier. Once the solution in the 's' domain is found, the inverse Laplace transform is used to convert it back to the time domain, yielding the solution to the original differential equation.

The Laplace Transform of Derivatives

The effectiveness of the Laplace transform in solving ODEs stems from its ability to transform derivatives of a function into algebraic expressions involving the transform of the function itself. This property allows us to convert a differential equation into an algebraic one.

Laplace Transform of the First Derivative

Let $f(t)$ be a function whose Laplace transform is $F(s) = \mathcal{L}\{f(t)\} = \int_0^\infty e^{-st} f(t) dt$. The Laplace transform of the first derivative, $f'(t)$, is given by:

$\mathcal{L}\{f'(t)\} = \int_0^\infty e^{-st} f'(t) dt$

We can solve this integral using integration by parts, with $u = e^{-st}$ and $dv = f'(t) dt$. Then $du = -se^{-st} dt$ and $v = f(t)$.

$\mathcal{L}\{f'(t)\} = [e^{-st} f(t)]_0^\infty - \int_0^\infty f(t) (-se^{-st}) dt$

Assuming $f(t)$ is of exponential order (meaning $|f(t)| \le Me^{at}$ for some constants M and a, and for all $t \ge 0$), the term $[e^{-st} f(t)]_0^\infty$ evaluates to $0 - f(0)$ as $s \to \infty$.

So, $\mathcal{L}\{f'(t)\} = 0 - f(0) + s \int_0^\infty e^{-st} f(t) dt$

$\mathcal{L}\{f'(t)\} = s F(s) - f(0)$

Laplace Transform of the Second Derivative

Similarly, we can find the Laplace transform of the second derivative, $f''(t)$. We know that $f''(t) = \frac{d}{dt}(f'(t))$. Let $g(t) = f'(t)$. Then $g'(t) = f''(t)$. Using the formula for the first derivative:

$\mathcal{L}\{f''(t)\} = \mathcal{L}\{g'(t)\} = s G(s) - g(0)$

Here, $G(s) = \mathcal{L}\{f'(t)\} = s F(s) - f(0)$. And $g(0) = f'(0)$.

Substituting these back:

$\mathcal{L}\{f''(t)\} = s [s F(s) - f(0)] - f'(0)$

$\mathcal{L}\{f''(t)\} = s^2 F(s) - s f(0) - f'(0)$

Laplace Transform of the nth Derivative

This pattern can be generalized for the nth derivative of a function $f(t)$, denoted as $f^{(n)}(t)$. The formula is:

$\mathcal{L}\{f^{(n)}(t)\} = s^n F(s) - s^{n-1} f(0) - s^{n-2} f'(0) - \dots - s f^{(n-2)}(0) - f^{(n-1)}(0)$

This formula is crucial because it incorporates the initial conditions $f(0), f'(0), \dots, f^{(n-1)}(0)$ directly into the transformed equation.

Solving Linear ODEs with Constant Coefficients

Consider a general nth-order linear ODE with constant coefficients:

$a_n y^{(n)}(t) + a_{n-1} y^{(n-1)}(t) + \dots + a_1 y'(t) + a_0 y(t) = R(t)$

with initial conditions $y(0), y'(0), \dots, y^{(n-1)}(0)$.

To solve this using Laplace transforms, we follow these steps:

Step 1: Take the Laplace Transform of the Entire Equation

Apply the Laplace transform operator $\mathcal{L}\{\cdot\}$ to both sides of the equation. Use the linearity property of the Laplace transform: $\mathcal{L}\{ay(t) + bz(t)\} = a\mathcal{L}\{y(t)\} + b\mathcal{L}\{z(t)\}$.

$a_n \mathcal{L}\{y^{(n)}(t)\} + a_{n-1} \mathcal{L}\{y^{(n-1)}(t)\} + \dots + a_1 \mathcal{L}\{y'(t)\} + a_0 \mathcal{L}\{y(t)\} = \mathcal{L}\{R(t)\}$

Step 2: Substitute the Laplace Transforms of Derivatives

Replace each term $\mathcal{L}\{y^{(k)}(t)\}$ with its corresponding formula involving $Y(s) = \mathcal{L}\{y(t)\}$ and the initial conditions $y(0), y'(0), \dots, y^{(k-1)}(0)$.

For example, for a second-order ODE $ay''(t) + by'(t) + cy(t) = R(t)$:

$a[s^2 Y(s) - s y(0) - y'(0)] + b[s Y(s) - y(0)] + c Y(s) = R(s)$

where $R(s) = \mathcal{L}\{R(t)\}$.

Step 3: Solve for Y(s)

Rearrange the equation obtained in Step 2 to isolate $Y(s)$. This will result in an algebraic equation for $Y(s)$.

In the second-order example:

$Y(s) [as^2 + bs + c] - a s y(0) - a y'(0) - b y(0) = R(s)$

$Y(s) [as^2 + bs + c] = R(s) + a s y(0) + a y'(0) + b y(0)$

$Y(s) = \frac{R(s) + a s y(0) + a y'(0) + b y(0)}{as^2 + bs + c}$

Step 4: Find the Inverse Laplace Transform

The final step is to find the inverse Laplace transform of $Y(s)$ to obtain the solution $y(t)$. This often involves techniques like partial fraction decomposition, especially when the denominator is a polynomial that can be factored.

$y(t) = \mathcal{L}^{-1}\{Y(s)\}$

Step 1: Take Laplace Transform $\mathcal{L}\{y''(t) - 3y'(t) + 2y(t)\} = \mathcal{L}\{0\}$ $\mathcal{L}\{y''(t)\} - 3\mathcal{L}\{y'(t)\} + 2\mathcal{L}\{y(t)\} = 0$

Step 2: Substitute Transforms of Derivatives Let $Y(s) = \mathcal{L}\{y(t)\}$. We know: $\mathcal{L}\{y'(t)\} = sY(s) - y(0) = sY(s) - 1$ $\mathcal{L}\{y''(t)\} = s^2Y(s) - s y(0) - y'(0) = s^2Y(s) - s(1) - 0 = s^2Y(s) - s$

Substituting these into the transformed equation: $(s^2Y(s) - s) - 3(sY(s) - 1) + 2Y(s) = 0$

Step 3: Solve for Y(s) $s^2Y(s) - s - 3sY(s) + 3 + 2Y(s) = 0$ Group terms with $Y(s)$: $Y(s)(s^2 - 3s + 2) - s + 3 = 0$ $Y(s)(s^2 - 3s + 2) = s - 3$ $Y(s) = \frac{s - 3}{s^2 - 3s + 2}$

Factor the denominator: $s^2 - 3s + 2 = (s-1)(s-2)$. $Y(s) = \frac{s - 3}{(s-1)(s-2)}$

Step 4: Inverse Laplace Transform We need to decompose $Y(s)$ using partial fractions. $\frac{s - 3}{(s-1)(s-2)} = \frac{A}{s-1} + \frac{B}{s-2}$ Multiply both sides by $(s-1)(s-2)$: $s - 3 = A(s-2) + B(s-1)$ To find A, set $s=1$: $1 - 3 = A(1-2) + B(1-1)$ $-2 = A(-1) + 0 \implies A = 2$ To find B, set $s=2$: $2 - 3 = A(2-2) + B(2-1)$ $-1 = 0 + B(1) \implies B = -1$ So, $Y(s) = \frac{2}{s-1} - \frac{1}{s-2}$.

Now, take the inverse Laplace transform using the standard formula $\mathcal{L}^{-1}\{\frac{1}{s-a}\} = e^{at}$. $y(t) = \mathcal{L}^{-1}\{Y(s)\} = \mathcal{L}^{-1}\{\frac{2}{s-1}\} - \mathcal{L}^{-1}\{\frac{1}{s-2}\}$ $y(t) = 2 \mathcal{L}^{-1}\{\frac{1}{s-1}\} - \mathcal{L}^{-1}\{\frac{1}{s-2}\}$ $y(t) = 2e^{1t} - e^{2t}$ $y(t) = 2e^t - e^{2t}$

This is the solution to the given ODE with the specified initial conditions.

Handling Non-Homogeneous Equations and Different RHS functions

The method extends seamlessly to non-homogeneous linear ODEs, where $R(t)$ is not identically zero. The Laplace transform of $R(t)$, denoted $R(s)$, is incorporated directly into the algebraic equation for $Y(s)$.

For example, if $R(t) = e^{kt}$, then $R(s) = \frac{1}{s-k}$. This $R(s)$ term will appear on the right-hand side when solving for $Y(s)$.

$Y(s) = \frac{R(s) + \text{terms from initial conditions}}{\text{characteristic polynomial}}$

The subsequent partial fraction decomposition might become more complex, potentially involving terms like $\frac{C}{s-k}$ if $k$ is not a root of the characteristic polynomial, or higher powers of $(s-k)$ if $k$ is a repeated root.

Example with a non-homogeneous term:

Solve $y'(t) + y(t) = \sin(t)$ with $y(0) = 0$.

Step 1 & 2: Transform and Substitute $\mathcal{L}\{y'(t)\} + \mathcal{L}\{y(t)\} = \mathcal{L}\{\sin(t)\}$ $[sY(s) - y(0)] + Y(s) = \frac{1}{s^2+1}$ Since $y(0)=0$: $sY(s) + Y(s) = \frac{1}{s^2+1}$

Step 3: Solve for Y(s) $Y(s)(s+1) = \frac{1}{s^2+1}$ $Y(s) = \frac{1}{(s+1)(s^2+1)}$

Step 4: Inverse Laplace Transform Use partial fractions: $\frac{1}{(s+1)(s^2+1)} = \frac{A}{s+1} + \frac{Bs+C}{s^2+1}$ $1 = A(s^2+1) + (Bs+C)(s+1)$ Set $s=-1$: $1 = A((-1)^2+1) + (B(-1)+C)(-1+1)$ $1 = A(2) + 0 \implies A = 1/2$ Expand the equation: $1 = A s^2 + A + B s^2 + Bs + Cs + C$ $1 = (A+B)s^2 + (B+C)s + (A+C)$ Equate coefficients: For $s^2$: $A+B = 0 \implies B = -A = -1/2$ For $s$: $B+C = 0 \implies C = -B = 1/2$ For constant: $A+C = 1/2 + 1/2 = 1$ (This checks out) So, $Y(s) = \frac{1/2}{s+1} + \frac{-1/2 s + 1/2}{s^2+1} = \frac{1}{2} \frac{1}{s+1} - \frac{1}{2} \frac{s}{s^2+1} + \frac{1}{2} \frac{1}{s^2+1}$

Now, take the inverse Laplace transform: $y(t) = \mathcal{L}^{-1}\{\frac{1}{2} \frac{1}{s+1}\} - \mathcal{L}^{-1}\{\frac{1}{2} \frac{s}{s^2+1}\} + \mathcal{L}^{-1}\{\frac{1}{2} \frac{1}{s^2+1}\}$ $y(t) = \frac{1}{2} e^{-t} - \frac{1}{2} \cos(t) + \frac{1}{2} \sin(t)$

This is the solution $y(t)$ for the given non-homogeneous ODE with the initial condition.

Advantages of Using Laplace Transforms

The Laplace transform method offers several significant advantages for solving ODEs:

• Handles Initial Conditions Naturally: Initial conditions are incorporated directly into the transformed equation, eliminating the need for separate steps to determine arbitrary constants.
• Converts Differential Equations to Algebraic Equations: The core transformation simplifies the problem from calculus (differentiation and integration) to algebra (solving for $Y(s)$).
• Effective for Discontinuous Forcing Functions: It can handle forcing functions $R(t)$ that are piecewise continuous, such as those involving step functions or impulses, which are more cumbersome with other methods.
• Systematic Approach: The procedure is systematic and generally applicable to any linear ODE with constant coefficients.

Limitations and Considerations

Despite its power, there are some limitations:

• Linear ODEs Only: The standard Laplace transform method is primarily for linear differential equations.
• Constant Coefficients Preferred: While extensions exist, the method is most straightforward for equations with constant coefficients. Variable coefficients lead to more complex transforms.
• Inverse Transform Complexity: The most challenging part can be finding the inverse Laplace transform, especially for complicated expressions of $Y(s)$. This often requires significant skill in partial fraction decomposition and knowledge of transform pairs.
• Requires Functions of Exponential Order: The derivation assumes that the function and its derivatives are of exponential order, which is true for most functions encountered in practical ODE problems.