Engineering Math Series 4: Second-Order Linear ODEs

While first-order differential equations dealt with "the rate of change of the current state," second-order differential equations handle **systems that include acceleration or curvature**. That is why they frequently appear when describing phenomena with inertia, such as vibration, circuits, and mechanical systems.

Standard Form

The most basic second-order linear differential equation is

$$ay'' + by' + cy = g(x)$$

where $a$, $b$, $c$ are constants and $g(x)$ is the external input. Looking first at the homogeneous equation without external input,

$$ay'' + by' + cy = 0$$

Why Think of Exponential Functions First

For constant-coefficient linear equations, substituting $y = e^{rx}$ maintains the form under differentiation. This leads to

$$ar^2 + br + c = 0$$

which is the **characteristic equation**, and the structure of the solution is determined by the roots of this quadratic equation.

Three Basic Cases

Two distinct real roots

If the roots are $r_1$ and $r_2$, then

$$y = C_1 e^{r_1x} + C_2 e^{r_2x}$$

Repeated root

If the root $r$ is repeated, then

$$y = (C_1 + C_2 x)e^{rx}$$

Complex conjugate roots

If the roots are $\alpha \pm i\beta$, then

$$y = e^{\alpha x}\left(C_1 \cos \beta x + C_2 \sin \beta x\right)$$

This last form is precisely what connects to damped oscillation.

Worked Example

Consider the following problem.

$$y'' - 3y' + 2y = 0$$

The characteristic equation is

$$r^2 - 3r + 2 = 0$$

which factors as

$$ (r-1)(r-2) = 0 $$

so the roots are $r=1$ and $r=2$. Therefore the general solution is

$$y = C_1 e^x + C_2 e^{2x}$$

If the initial conditions are

$$y(0)=1, \quad y'(0)=0$$

then

$$C_1 + C_2 = 1$$

$$C_1 + 2C_2 = 0$$

Solving gives

$$C_1 = 2, \quad C_2 = -1$$

so

$$y = 2e^x - e^{2x}$$

Engineering Applications

The natural response of a mass-spring-damper system is typically written as

$$m x'' + c x' + kx = 0$$

where

- $m$ is the mass

- $c$ is the damping coefficient

- $k$ is the spring constant

In circuits, the same structure appears as

$$Lq'' + Rq' + \frac{1}{C}q = 0$$

in an RLC circuit. In other words, mechanical systems and electrical systems share the same mathematical structure.

How to Interpret the Solution

In engineering mathematics, you must always read the behavior after finding the solution.

- Positive real roots in the exponential mean divergence

- Negative real roots mean decay

- Complex roots mean oscillation

- Complex roots with negative real part mean damped oscillation

In other words, the equations ultimately serve as tools that tell you about the stability and response speed of the system.

Common Mistakes

Writing the solution incorrectly after solving the characteristic equation

The solution forms for repeated roots and complex roots are different. If you mechanically write the general solution after only finding the roots, mistakes are easy.

Being careless with differentiation when applying initial conditions

If you do not compute $y'$ accurately and just try to match the constants, errors arise. It is safer to first organize the general solution and then differentiate carefully.

Not reading the meaning of the solution

In engineering problems, whether the solution diverges, decays, or oscillates is the key point. It is important to develop the habit of not skipping system analysis after completing the calculation.

One-Line Summary

Second-order linear ODEs are tools for reading the form of the natural response through the roots of the characteristic equation.

Next Post Preview

In the next post, we will connect second-order equations to real physical systems and examine how **damping, forced oscillation, and resonance** arise.

References

- Erwin Kreyszig, _Advanced Engineering Mathematics_, 10th Edition

- Steven H. Strogatz, _Nonlinear Dynamics and Chaos_

- MIT OpenCourseWare, Vibrations and Waves materials