Higher-Order Non-Homogeneous DEs - Theory & Concepts

Higher-Order Non-Homogeneous DEs

A linear differential equation is non-homogeneous if the right-hand side is a non-zero function of

x

General Solution

The general solution is the sum of the complementary solution (

y_c

, solution to the associated homogeneous equation

L[y]=0

) and a particular solution (

y_p

, any specific solution to the full non-homogeneous equation

L[y]=f(x)

y = y_c + y_p

Superposition Principle

If the non-homogeneous term

f(x)

consists of a sum of multiple functions, say

f(x) = f_1(x) + f_2(x) + \dots + f_k(x)

, then we can find individual particular solutions

y_{p1}, y_{p2}, \dots, y_{pk}

corresponding to each individual forcing term

f_i(x)

. The overall particular solution is simply their linear combination:

y_p = y_{p1} + y_{p2} + \dots + y_{pk}

This drastically simplifies solving for complex

f(x)

terms.

Method 1: Undetermined Coefficients

This method is used when

f(x)

is a polynomial, an exponential, a sine, or a cosine function, or finite sums and products of these forms. We make an "educated guess" for the form of

y_p

based on the form of

f(x)

and its derivatives.

Guessing yp

Form of $f(x)$	Assumed Form of $y_p$
$k$ (Constant)	$A$
$ax^n + \dots + cx + d$ (Polynomial)	$A_n x^n + \dots + A_1 x + A_0$
$ce^{kx}$ (Exponential)	$A e^{kx}$
$c \sin(kx)$ or $c \cos(kx)$	$A \sin(kx) + B \cos(kx)$
$(x^n + \dots)e^{kx}$ (Product)	$(A_n x^n + \dots + A_0)e^{kx}$
$e^{kx} \sin(\omega x)$ or $e^{kx} \cos(\omega x)$	$e^{kx} (A \sin(\omega x) + B \cos(\omega x))$

The Multiplication Rule: Compare your assumed form for

y_p

with the complementary solution

y_c

. If any term in your initial guess for

y_p

is a duplicate of a term in

y_c

(i.e., it is a solution to the homogeneous equation), multiply the entire assumed

y_p

guess by

x^s

, where

s

is the smallest positive integer that eliminates the duplication.

Annihilator Approach

An alternative way to formally find

y_p

is using the Annihilator Approach with differential operators. An operator

L

annihilates a function

f(x)

L[f(x)] = 0

Common Annihilators (where

D = d/dx

$f(x) = x^k$ is annihilated by $D^{k+1}$ .
$f(x) = e^{ax}$ is annihilated by $(D - a)$ .
$f(x) = \sin(\beta x)$ or $\cos(\beta x)$ is annihilated by $(D^2 + \beta^2)$ .
$f(x) = e^{\alpha x}\sin(\beta x)$ or $e^{\alpha x}\cos(\beta x)$ is annihilated by $(D^2 - 2\alpha D + (\alpha^2 + \beta^2))$ .

Method: Apply the annihilator to both sides of the non-homogeneous equation

L_{orig}[y] = f(x)

to get

L_{ann}[L_{orig}[y]] = 0

, a higher-order homogeneous equation. Solve this new homogeneous equation for its general solution. Discard the terms that belong to the original

y_c

; the remaining terms form the structure of

y_p

. Then, substitute this

y_p

back into the original non-homogeneous DE to solve for the undetermined coefficients.

Method 2: Variation of Parameters

This is a general method that works for any

f(x)

, even if it is not a finite sum of polynomials, exponentials, or sinusoids (e.g.,

f(x) = \sec(x), \tan(x), \ln(x)

Formula for Variation of Parameters

For a second-order linear DE

y'' + P(x)y' + Q(x)y = f(x)

with independent complementary solutions

y_1

and

y_2

y_p = u_1(x)y_1(x) + u_2(x)y_2(x)

The functions

u_1

and

u_2

are calculated by:

u_1 = \int \frac{-y_2 f(x)}{W} \, dx \quad \text{and} \quad u_2 = \int \frac{y_1 f(x)}{W} \, dx

where

W = W(y_1, y_2) = y_1 y_2' - y_2 y_1'

is the Wronskian of the fundamental solution set.

Wronskian Properties

The Wronskian, denoted as

W(y_1, y_2, \dots, y_n)

, determines whether a set of solutions to a homogeneous differential equation is linearly independent:

If $W \neq 0$ at any point in the interval, the functions are linearly independent.
If $W = 0$ everywhere in the interval, and the functions are solutions to a regular linear ODE, they are linearly dependent.
Abel's Identity states that $W(x) = W(x_0)e^{-\int_{x_0}^x P(t)dt}$ . This means the Wronskian is either never zero or always zero on an interval where $P(x)$ is continuous.

Key Takeaways

General Solution Form: Always write the solution as $y = y_c + y_p$ . Find $y_c$ first.
Superposition Principle: If $f(x) = f_1(x) + f_2(x)$ , solve for $y_{p1}$ and $y_{p2}$ separately and sum them.
Undetermined Coefficients: Fast and effective for polynomial, exponential, and sinusoidal forcing functions, including their sums and products. The assumed form of $y_p$ mirrors $f(x)$ .
The Multiplication Rule: Crucial step. If your guessed $y_p$ contains terms from $y_c$ , multiply the guess by $x$ (or $x^2$ ) until no duplication exists.
Variation of Parameters: A universally applicable method using the Wronskian $W$ . It works for any continuous $f(x)$ , but relies on integrating potentially complex rational functions.

PreviousHigher-Order Homogeneous DEs - Examples & Applications

Quiz Me

NextHigher-Order Non-Homogeneous DEs - Examples & Applications