Examples of the Solution Methods for Constant Coefficient Linear Differential Equations

Introduction

Linear differential equations with constant coefficients are fundamental in mathematical analysis, offering structured frameworks for modeling physical phenomena such as oscillatory motion, heat conduction, and electrical circuits. These equations are characterized by their ability to be solved using systematic methods rooted in algebraic and analytical techniques. The core challenge lies in determining the general solution, which often involves identifying the homogeneous solution and a particular solution. This article explores the most common solution methods for such equations, emphasizing their theoretical foundations and practical applications.

Core Content

Characteristic Equation Method

The characteristic equation method is a cornerstone of solving linear differential equations with constant coefficients. For a second-order equation of the form $ ay'' + by' + cy = 0 $, the characteristic equation $ ar^2 + br + c = 0 $ is derived by assuming a solution of the form $ y = e^{rx} $. The roots of this equation determine the nature of the solution:

Real and distinct roots ($ r_1 $ and $ r_2 $): The general solution is $ y(x) = C_1 e^{r_1 x} + C_2 e^{r_2 x} $.
Complex conjugate roots ($ r = \alpha \pm i\beta $): The solution involves sines and cosines, $ y(x) = C_1 \cos(\beta x) + C_2 \sin(\beta x) $.
Real and repeated roots ($ r_1 = r_2 $): The solution is $ y(x) = C_1 e^{r x} + C_2 x e^{r x} $.

This method relies on the exponential behavior of solutions and is particularly effective for equations with smooth, well-defined behavior. For example, the equation $ y'' + y = 0 $ has roots $ \pm i $, leading to the oscillatory solution $ y(x) = C_1 \cos(x) + C_2 \sin(x) $.

Method of Undetermined Coefficients

The method of undetermined coefficients is applicable to linear differential equations with constant coefficients when the nonhomogeneous term is a polynomial, exponential, sine, or cosine function. This method involves assuming a form for the particular solution based on the nonhomogeneous term and then solving for the coefficients.

For instance, consider the equation $ y'' + 3y' + 2y = 5 $. The homogeneous solution is $ y_h = C_1 e^{x} + C_2 e^{2x} $, and the nonhomogeneous term $ 5 $ is a constant. A particular solution is assumed to be a constant, $ y_p = A $. Substituting into the equation:
$$ A'' + 3A' + 2A = 5 \implies 0 + 0 + 2A = 5 \implies A = \frac{5}{2}. $$
Thus, the general solution is $ y(x) = C_1 e^{x} + C_2 e^{2x} + \frac{5}{2} $. This method is efficient for equations with simple nonhomogeneous terms but requires careful selection of the assumed form to avoid duplication with the homogeneous solution.

Variation of Parameters

When the nonhomogeneous term is not easily expressible via the method of undetermined coefficients, variation of parameters provides an alternative approach. This method involves expressing the particular solution as a linear combination of the homogeneous solutions, with coefficients that are functions of $ x $.

For example, consider the equation $ y'' + 4y = \sin(x) $. The homogeneous solution is $ y_h = C_1 \cos(2x) + C_2 \sin(2x) $. To find a particular solution, assume $ y_p = u_1 \cos(2x) + u_2 \sin(2x) $. Using the method of variation of parameters, the coefficients $ u_1 $ and $ u_2 $ are determined by solving the system:
$$ \begin{cases} -2u_1 \sin(2x) + 2u_2 \cos(2x) = \sin(x) \ 2u_1 \cos(2x) - 2u_2 \sin(2x) = 0 \end{cases} $$
Solving this system yields the particular solution $ y_p = \frac{1}{4} \sin(2x) + \frac{1}{4} \cos(2x) $. The general solution combines $ y_h $ and $ y_p $, demonstrating the method's versatility in handling nonhomogeneous terms with trigonometric functions.

Laplace Transforms

Laplace transforms offer a powerful tool for solving linear differential equations with constant coefficients, particularly when initial conditions are provided. The Laplace transform converts differential equations into algebraic equations, simplifying their solution.

For example, consider the equation $ y'' + 4y = \sin(x) $ with initial conditions $ y(0) = 0 $ and $ y'(0) = 0 $. Taking the Laplace transform of both sides:
$$ s^2 Y(s) - sY(0) - Y'(0) + 4Y(s) = \frac{1}{s^2 + 1} $$
Substituting the initial conditions:
$$ s^2 Y(s) + 4Y(s) = \frac{1}{s^2