ODE: Complex Roots and Oscillatory Solutions

When you design a suspension system for a car or filter noise from an electrical signal, you're not just solving equations—you're predicting and controlling oscillations. The transition from solving a quadratic equation to designing a comfortable ride hinges on one powerful mathematical idea: when the characteristic roots of a second-order linear ordinary differential equation (ODE) become complex, the solutions become oscillatory. This concept is the bridge between abstract math and the real-world vibratory behavior of countless engineering systems.

From Real to Complex: The Birth of Oscillations

Consider the canonical second-order, linear, homogeneous ODE with constant coefficients: $$a y'' + b y' + c y = 0.$$ The standard approach is to propose an exponential solution of the form $y = e^{rt}$, leading to the characteristic equation $a r^2 + b r + c = 0$. Using the quadratic formula, the roots are: $$r = \frac{-b \pm \sqrt{b^2 - 4ac}}{2a}.$$

Oscillations emerge when the discriminant is negative: $b^2 - 4ac < 0$. This makes the roots complex conjugates, which we write as $r = \alpha \pm i \beta$, where:

• $\alpha = -b/(2a)$ is the real part.
• $\beta = \sqrt{4ac - b^2}/(2a)$ is the imaginary part.

The parameter $\alpha$ dictates growth or decay, while $\beta$ is the key to oscillation. This condition, common in systems with low damping, is precisely what defines an underdamped response in engineering. The two linearly independent solutions are now complex exponentials: $e^{(\alpha + i\beta)t}$ and $e^{(\alpha - i\beta)t}$.

Euler's Formula: The Key to Real-Valued Solutions

Our fundamental solutions are complex, but we seek real-valued functions to describe physical quantities like displacement or voltage. The indispensable tool is Euler's formula: $$e^{i\theta} = \cos \theta + i \sin \theta.$$ This identity connects complex exponentials to the sinusoidal functions that model pure oscillation. Applying it to our complex roots, we expand: $$e^{(\alpha + i\beta)t} = e^{\alpha t} e^{i\beta t} = e^{\alpha t} \left( \cos(\beta t) + i \sin(\beta t) \right).$$

Since our ODE is linear and has real coefficients, both the real and imaginary parts of this complex solution are themselves real-valued solutions. Therefore, we can construct our fundamental real solution set by taking linear combinations of the two complex exponential solutions. The standard, most useful general real solution is: $$y(t) = e^{\alpha t} \left( C_1 \cos(\beta t) + C_2 \sin(\beta t) \right),$$ where $C_1$ and $C_2$ are real constants determined by initial conditions. This form explicitly reveals the solution's structure: an exponential envelope $e^{\alpha t}$ modulating an oscillation at angular frequency $\beta$.

Interpreting Amplitude, Phase, Frequency, and Period

The general solution with sine and cosine can be consolidated into a single sine (or cosine) function with a phase shift, offering a clearer physical interpretation. Using a trigonometric identity, we can rewrite: $$C_1 \cos(\beta t) + C_2 \sin(\beta t) = A \sin(\beta t + \phi),$$ where:

• Amplitude $A = \sqrt{C_1^2 + C_2^2}$ is the peak magnitude of the oscillatory component.
• Phase angle $\phi = \arctan(C_1 / C_2)$ determines the initial offset of the oscillation.

The complete solution becomes: $$y(t) = A e^{\alpha t} \sin(\beta t + \phi).$$

This is the most insightful form. The term $\beta$, the imaginary part of the characteristic root, is the angular frequency (or circular frequency) of oscillation, measured in radians per unit time. From $\beta$, we derive two critical temporal measures:

• The period $T = \frac{2\pi}{\beta}$, which is the time for one complete cycle.
• The frequency $f = \frac{\beta}{2\pi} = \frac{1}{T}$, which is the number of cycles per unit time.

The real part $\alpha$ governs the exponential envelope. If $\alpha < 0$ (the common case in stable systems), the oscillation decays; if $\alpha > 0$, it grows.

Physical Applications to Underdamped Systems

The power of this mathematics is realized in modeling real underdamped oscillations. This occurs when damping is present but insufficient to prevent oscillation ($b^2 < 4ac$). Let's examine two quintessential engineering applications.

Mechanical Systems (Mass-Spring-Damper): For a system with mass $m$, damping coefficient $c$, and spring constant $k$, the equation of motion is $m x'' + c x' + k x = 0$. The roots are complex if $c^2 < 4mk$. Here:

• $\alpha = -c/(2m)$: The decay rate is proportional to the damping.
• $\beta = \sqrt{4mk - c^2}/(2m) = \sqrt{\frac{k}{m} - \left(\frac{c}{2m}\right)^2}$: The oscillation frequency, often called the damped natural frequency. Notice it is less than the undamped natural frequency $\omega_n = \sqrt{k/m}$ due to damping.

The solution $x(t) = A e^{-(c/(2m)) t} \sin(\beta t + \phi)$ models the displacement of a car's suspension after hitting a bump or a building swaying in the wind.

Electrical Systems (RLC Circuit): For a series RLC circuit with resistance $R$, inductance $L$, and capacitance $C$, the charge $q$ on the capacitor obeys $L q'' + R q' + (1/C) q = 0$. Complex roots arise when $R^2 < 4L/C$.

• $\alpha = -R/(2L)$: Decay depends on resistance.
• $\beta = \sqrt{4L/C - R^2}/(2L) = \sqrt{\frac{1}{LC} - \left(\frac{R}{2L}\right)^2}$: The damped oscillation frequency of the charge and current.

This models the ringing response of a tuned circuit or the transient behavior in power systems after a fault is cleared.

Common Pitfalls

1. Misidentifying the Oscillation Frequency: A frequent error is using the undamped natural frequency $\omega_n$ instead of the damped frequency $\beta$ in the sine/cosine terms. Remember, $\beta = \sqrt{\omega_n^2 - \alpha^2}$ for standard forms. The oscillation you observe in an underdamped system is at frequency $\beta$, which is always less than $\omega_n$.

1. Forgetting the Exponential Envelope: When sketching or interpreting solutions, students sometimes draw a pure sine wave. You must remember the $e^{\alpha t}$ multiplier. For stable systems ($\alpha < 0$), the oscillations are contained within the decaying envelope curves $\pm A e^{\alpha t}$.

1. Incorrect Phase Angle Calculation: When converting $C_1 \cos(\beta t) + C_2 \sin(\beta t)$ to $A \sin(\beta t + \phi)$, the formula $\phi = \arctan(C_1 / C_2)$ is correct, but you must consider the quadrant. $C_1$ is the coefficient of cosine and corresponds to the sine of the phase in the identity. Using the atan2(C1, C2) function in computation avoids sign errors.

1. Applying to the Wrong Condition: The oscillatory solution form is only valid for the underdamped case ($b^2 - 4ac < 0$). Applying it to critically damped or overdamped systems (with real roots) is incorrect. Always check the discriminant first.

Summary

• Complex conjugate roots $r = \alpha \pm i \beta$ to the characteristic equation of a constant-coefficient linear ODE are the direct mathematical source of oscillatory solutions, describing underdamped system behavior.
• Euler's formula $e^{i\theta} = \cos \theta + i \sin \theta$ is essential for converting complex-exponential solutions into the real-valued form $y(t) = e^{\alpha t} (C_1 \cos \beta t + C_2 \sin \beta t)$.
• The solution is best interpreted as $y(t) = A e^{\alpha t} \sin(\beta t + \phi)$, where $A$ is the initial amplitude, $\phi$ is the phase angle, $\beta$ is the damped angular frequency governing oscillation, and $\alpha$ (typically negative) controls the exponential decay.
• The period of oscillation is $T = 2\pi / \beta$, and the frequency is $f = \beta / (2\pi)$.
• These solutions are ubiquitous in engineering, modeling phenomena from the decaying vibration of a mass-spring-damper system to the ringing current in an RLC circuit, providing the predictive power needed for design and analysis.