Homogeneous Linear Equations with Constant Coefficients & Cauchy-Euler Equations
Complete algebraic theory of homogeneous equations with constant coefficients: characteristic polynomials, distinct real roots, repeated root multipliers, complex conjugate roots and Euler's harmonic representation, and Cauchy-Euler (equidimensional) equations via the logarithmic transformation x = e^t.
ยง6.1 Characteristic Polynomials & Roots Classification
1. The Second-Order Constant Coefficient Equation
Consider the homogeneous linear ODE with constant real coefficients:
Assuming a trial solution of exponential form $y(x) = e^{rx}$, we compute $y' = r e^{rx}$ and $y'' = r^2 e^{rx}$. Substituting into the ODE:
Since $e^{rx} \neq 0$ for all real $x$, $r$ must satisfy the characteristic (auxiliary) equation:
2. The Three Root Regimes
- Case 1: Real and Distinct Roots ($\Delta = b^2 - 4ac > 0$)
The roots $r_1 \neq r_2$ are real numbers. Two linearly independent solutions are $y_1 = e^{r_1 x}$ and $y_2 = e^{r_2 x}$. The general solution is: $$y(x) = c_1 e^{r_1 x} + c_2 e^{r_2 x}$$ - Case 2: Real and Repeated Roots ($\Delta = b^2 - 4ac = 0$)
There is a single root of multiplicity two: $r = -b / (2a)$. One solution is $y_1 = e^{rx}$. Using reduction of order, the second independent solution is $y_2 = x e^{rx}$. The general solution is: $$y(x) = (c_1 + c_2 x) e^{rx}$$ - Case 3: Complex Conjugate Roots ($\Delta = b^2 - 4ac < 0$)
The roots are $r = \alpha \pm i \beta$, where $\alpha = -b/(2a)$ and $\beta = \frac{\sqrt{4ac - b^2}}{2a} > 0$. Using Euler's formula $e^{(\alpha \pm i\beta)x} = e^{\alpha x}(\cos\beta x \pm i \sin\beta x)$, the real linearly independent fundamental solutions are $y_1 = e^{\alpha x} \cos(\beta x)$ and $y_2 = e^{\alpha x} \sin(\beta x)$. The general solution is: $$y(x) = e^{\alpha x} \left[ c_1 \cos(\beta x) + c_2 \sin(\beta x) \right] = R e^{\alpha x} \cos(\beta x - \delta)$$
ยง6.2 Higher-Order Constant Coefficient Equations ($n$-th Order)
1. Generalization to $n$-th Order Equations
For an $n$-th order linear ODE $a_n y^{(n)} + a_{n-1} y^{(n-1)} + \dots + a_1 y' + a_0 y = 0$, the characteristic polynomial is:
By the Fundamental Theorem of Algebra, $P(r)$ has exactly $n$ complex roots (counting multiplicities):
- Each real root $r$ of multiplicity $k$ contributes $k$ linearly independent solutions: $$e^{rx}, \quad x e^{rx}, \quad x^2 e^{rx}, \quad \dots, \quad x^{k-1} e^{rx}$$
- Each complex conjugate pair $\alpha \pm i\beta$ of multiplicity $k$ contributes $2k$ linearly independent solutions: $$\begin{aligned} e^{\alpha x} \cos\beta x, \quad x e^{\alpha x} \cos\beta x, \quad \dots, \quad x^{k-1} e^{\alpha x} \cos\beta x \\ e^{\alpha x} \sin\beta x, \quad x e^{\alpha x} \sin\beta x, \quad \dots, \quad x^{k-1} e^{\alpha x} \sin\beta x \end{aligned}$$
The sum of all these solutions multiplied by arbitrary constants $c_1, \dots, c_n$ forms the complete general solution.
ยง6.3 Cauchy-Euler (Equidimensional) Differential Equations
1. Standard Form of Cauchy-Euler Equations
A linear differential equation of the form:
where the power of $x$ matches the order of the derivative in each term, is called a Cauchy-Euler (or equidimensional) equation.
2. Second-Order Homogeneous Cauchy-Euler Equation
Consider $a x^2 y'' + b x y' + c y = 0$ for $x > 0$. We seek solutions of the form $y = x^m$. Differentiating:
Substituting into the ODE:
This yields the indicial (auxiliary) equation:
3. Three Cases of Solutions for Cauchy-Euler Equations
- Distinct Real Roots ($m_1 \neq m_2$): $y(x) = c_1 x^{m_1} + c_2 x^{m_2}$.
- Repeated Real Root ($m_1 = m_2 = m$): The second solution is obtained by logarithmic scaling: $$y(x) = x^m (c_1 + c_2 \ln x)$$
- Complex Conjugate Roots ($m = \alpha \pm i \beta$): Using $x^{\alpha \pm i\beta} = x^\alpha e^{\pm i \beta \ln x} = x^\alpha [\cos(\beta \ln x) \pm i \sin(\beta \ln x)]$: $$y(x) = x^\alpha \left[ c_1 \cos(\beta \ln x) + c_2 \sin(\beta \ln x) \right]$$
Step-by-Step Solved Examination Problems
Comprehensive analytical derivations, multi-tier solutions (Foundational, Intermediate Exam, and Honors/Proof Challenge) with complete line-by-line verification.
Find the general solution of the fourth-order differential equation:
Step 1: Write the characteristic equation
Step 2: Factor the polynomial
Step 3: Identify the roots
The four roots are:
Step 4: Formulate the fundamental solutions
For real roots: $y_1 = e^{2x}, y_2 = e^{-2x}$.
For imaginary roots $\alpha = 0, \beta = 2$: $y_3 = \cos(2x), y_4 = \sin(2x)$.
Step 5: Write the general solution
$y(x) = c_1 e^{2x} + c_2 e^{-2x} + c_3 \cos(2x) + c_4 \sin(2x)$ (or $A \cosh(2x) + B \sinh(2x) + c_3 \cos(2x) + c_4 \sin(2x)$).
Solve the Cauchy-Euler boundary value problem:
Step 1: Set up the indicial equation
With $y = x^m$, the equation becomes:
Step 2: Solve for $m$
Here $\alpha = 2, \beta = 3$.
Step 3: General solution
Step 4: Apply boundary conditions
At $x = 1$ ($\ln 1 = 0$):
At $x = e^{\pi/6}$ ($\ln(e^{\pi/6}) = \pi/6$):
Step 5: Final solution
$y(x) = 2x^2 \cos(3 \ln x)$.
An unforced RLC circuit is described by $L q'' + R q' + \frac{1}{C} q = 0$. For fixed $L = 1\text{ H}$ and $C = 0.25\text{ F}$, determine the critical damping resistance $R_{\text{crit}}$. If the circuit starts with initial charge $q(0) = Q_0$ and zero current $q'(0) = 0$, solve the critical IVP and prove that the charge never crosses zero for $t > 0$.
Step 1: Indicial equation and critical damping
Characteristic equation: $L r^2 + R r + \frac{1}{C} = 0 \implies r^2 + R r + 4 = 0$.
The discriminant is $\Delta = R^2 - 4(1)(4) = R^2 - 16$.
Critical damping occurs when $\Delta = 0 \implies R_{\\text{crit}} = \sqrt{16} = 4\\,\\Omega$.
Step 2: Repeated root solution
With $R = 4\\,\\Omega$, the repeated root is $r = -R/(2L) = -4/2 = -2$.
The general solution for the charge is:
Step 3: Apply initial conditions
$q(0) = Q_0 \implies c_1 = Q_0$.
Compute the current (derivative):
At $t = 0$: $q'(0) = c_2 - 2c_1 = 0 \implies c_2 = 2c_1 = 2Q_0$.
Thus:
Step 4: Zero-crossing proof
For $t > 0$, since $Q_0 > 0$, $1 + 2t > 1 > 0$ and $e^{-2t} > 0$. Thus $q(t) > 0$ for all $t \ge 0$. The charge monotonically decays toward zero as $t \to \infty$ without ever oscillating or crossing zero, proving the defining property of critical damping.
$R_{\text{crit}} = 4\,\Omega$. Solution: $q(t) = Q_0 (1 + 2t) e^{-2t}$. Since $1 + 2t > 0$ for all $t > 0$, $q(t) > 0$ strictly, verifying that the charge never crosses zero.