Definition and General Form
Homogeneous Second Order ODEs
Form: y'' + p(x)y' + q(x)y = 0. Homogeneous: zero right-hand side. Linear: dependent variable and derivatives appear to first power. Order: highest derivative is second.
Linearity and Homogeneity
Superposition principle applies: if y1, y2 solutions, then c1y1 + c2y2 solution. No external forcing term: purely natural system response.
Initial and Boundary Conditions
Two conditions required for unique solution. Examples: y(x0) = y0, y'(x0) = y0'. Essential for physical interpretation and numerical methods.
Characteristic Equation
Definition
Algebraic equation derived by substituting y = e^{rx} into linear constant coefficient ODE. Converts differential equation to polynomial equation: ar^2 + br + c = 0.
Derivation
Substitute y = e^{rx}, then y' = re^{rx}, y'' = r^2 e^{rx}, divide by e^{rx} (nonzero), get quadratic in r.
Role in Solutions
Roots determine solution type: distinct real roots, repeated roots, or complex conjugates. Root nature dictates exponential, polynomial, or trigonometric solutions.
Solution Structure
General Solution
General solution: linear combination of two linearly independent solutions y = c1 y1 + c2 y2. Constant coefficients yield closed-form expressions.
Complementary Function
Also called homogeneous solution. Forms basis for general solutions of nonhomogeneous equations when combined with particular integral.
Dependence on Roots
Root types (real distinct, repeated, complex) yield different forms of y1, y2. Ensures completeness and linear independence.
Constant Coefficients Case
Equation Form
Standard form: a y'' + b y' + c y = 0 with a, b, c constants. Simplifies analysis and solution via characteristic equation.
Characteristic Polynomial
Polynomial: a r^2 + b r + c = 0. Solutions r1, r2 found using quadratic formula: roots dictate solution form.
Solution Methods
Direct substitution, quadratic formula, and verification of linear independence. Closed form solutions available.
Distinct Real Roots
Condition
Discriminant D = b^2 - 4ac > 0. Two distinct real roots r1 ≠ r2.
General Solution
y = c1 e^{r1 x} + c2 e^{r2 x}. Solutions exponential with different growth/decay rates.
Behavior
System response: sum of two exponentials. Stability depends on sign of roots. Real negative roots indicate decay.
Characteristic equation: a r^2 + b r + c = 0Roots: r1, r2 = (-b ± √(b² - 4ac)) / 2aGeneral solution: y = c1 e^{r1 x} + c2 e^{r2 x} Repeated Roots
Condition
Discriminant D = b^2 - 4ac = 0. Double root r = -b/(2a).
General Solution
y = (c1 + c2 x) e^{r x}. Second solution multiplied by x to ensure linear independence.
Derivation
Reduction of order or limit of distinct roots solution as roots converge. Ensures complete basis.
If r1 = r2 = r:y1 = e^{r x}y2 = x e^{r x}General solution: y = (c1 + c2 x) e^{r x} Complex Conjugate Roots
Condition
Discriminant D = b^2 - 4ac < 0. Roots: r = α ± i β, α, β ∈ ℝ, β ≠ 0.
General Solution
y = e^{α x} (c1 cos β x + c2 sin β x). Oscillatory solutions modulated by exponential envelope.
Physical Interpretation
Damped oscillations: α < 0 damping, α = 0 undamped, α > 0 growth. Common in mechanical vibrations, circuits.
Reduction of Order Method
Purpose
Find second solution y2 given one known solution y1. Applicable when standard methods fail or coefficients variable.
Procedure
Assume y2 = v(x) y1(x), substitute into ODE, derive first order ODE for v'. Integrate to find v and thus y2.
Formula
v' = C exp(-∫P(x) dx) / (y1)^2 where ODE in standard form y'' + P(x) y' + Q(x) y = 0.
Given y1, set y2 = v y1Substitute into ODE, obtain:v'' y1 + 2 v' y1' + v y1'' + P(x)(v' y1 + v y1') + Q(x) v y1 = 0Simplify, solve for v', integrate to find y2 Euler-Cauchy Equations
Form
x^2 y'' + a x y' + b y = 0, a, b constants. Variable coefficients but solvable by substitution.
Solution Method
Try y = x^m, reduce to characteristic equation m^2 + (a - 1)m + b = 0. Roots m1, m2 determine solution.
Cases
Distinct roots: y = c1 x^{m1} + c2 x^{m2}. Repeated roots: y = c1 x^{m} + c2 x^{m} ln x. Complex roots: y = x^{α} (c1 cos β ln x + c2 sin β ln x).
Applications and Examples
Mechanical Vibrations
Mass-spring-damper modeled by m y'' + c y' + k y = 0. Solution type determines damping regime: underdamped, overdamped, critically damped.
Electrical Circuits
RLC circuit differential equation analogous to mechanical system. Current or voltage follow second order homogeneous ODEs.
Population Dynamics
Some simplified models reduce to linear homogeneous ODEs for growth or decay phases.
| Application | ODE Form | Physical Interpretation |
|---|---|---|
| Mass-Spring-Damper | m y'' + c y' + k y = 0 | Vibration and damping response |
| RLC Circuit | L y'' + R y' + (1/C) y = 0 | Current/voltage oscillations |
Summary Tables
Characteristic Roots and Solutions
| Root Type | Condition | General Solution |
|---|---|---|
| Distinct Real Roots | D > 0 | y = c1 e^{r1 x} + c2 e^{r2 x} |
| Repeated Roots | D = 0 | y = (c1 + c2 x) e^{r x} |
| Complex Roots | D < 0 | y = e^{α x} (c1 cos β x + c2 sin β x) |
Euler-Cauchy Equation Cases
| Case | Characteristic Roots | Solution Form |
|---|---|---|
| Distinct Real Roots | m1 ≠ m2 ∈ ℝ | y = c1 x^{m1} + c2 x^{m2} |
| Repeated Root | m1 = m2 = m | y = c1 x^{m} + c2 x^{m} ln x |
| Complex Roots | m = α ± i β | y = x^{α} (c1 cos β ln x + c2 sin β ln x) |
References
- Coddington, E. A., & Levinson, N. "Theory of Ordinary Differential Equations." McGraw-Hill, 1955, pp. 1-320.
- Boyce, W. E., & DiPrima, R. C. "Elementary Differential Equations and Boundary Value Problems," 10th ed., Wiley, 2012, pp. 100-150.
- Ince, E. L. "Ordinary Differential Equations." Dover Publications, 1956, pp. 200-280.
- Butcher, J. C. "Numerical Methods for Ordinary Differential Equations," 2nd ed., Wiley, 2008, pp. 50-90.
- Smith, H. L. "An Introduction to Ordinary Differential Equations." Cambridge University Press, 2013, pp. 75-130.