Linear Differential Equations with Constant Coefficients

Welcome to the study of linear differential equations with constant coefficients. These equations form a fundamental cornerstone in the field of mathematics and have widespread applications in physics, engineering, economics, and biology. Understanding their properties and methods of solution is crucial for anyone looking to model and analyze dynamic systems.

1. Introduction to Linear Differential Equations with Constant Coefficients

A linear differential equation is an equation that is linear in its dependent variable and its derivatives. When the coefficients of these derivatives are constants, we call it a linear differential equation with constant coefficients.

The general form of a linear differential equation of order 'n' is:

an(x) y(n) + an-1(x) y(n-1) + ... + a1(x) y' + a0(x) y = g(x)

Where y(k) denotes the k-th derivative of y with respect to x, and ai(x) are functions of x.

However, when the coefficients are constants, i.e., an, an-1, ..., a0 are constants, the equation simplifies to:

an y(n) + an-1 y(n-1) + ... + a1 y' + a0 y = g(x)

Here, an ≠ 0.

If g(x) = 0, the equation is called a homogeneous linear differential equation with constant coefficients:

an y(n) + an-1 y(n-1) + ... + a1 y' + a0 y = 0

If g(x) ≠ 0, it is a non-homogeneous equation.

In this unit, we will primarily focus on second-order equations, which are of the form:

ay'' + by' + cy = g(x)

where 'a', 'b', and 'c' are constants, and a ≠ 0.

2. Existence of Solutions

For linear differential equations, the existence and uniqueness of solutions are guaranteed under certain conditions.

Consider the second-order linear differential equation:

y'' + p(x) y' + q(x) y = f(x)

where p(x), q(x), and f(x) are continuous functions on an interval I.

The Existence and Uniqueness Theorem states that if p(x), q(x), and f(x) are continuous on an open interval I containing a point x0, then there exists a unique solution y(x) to the initial value problem (IVP):

y'' + p(x) y' + q(x) y = f(x)

satisfying the initial conditions:

y(x0) = y0 and y'(x0) = y'0

for any given initial values y0 and y'0. The unique solution is defined on the entire interval I.

For linear differential equations with constant coefficients, the coefficients p(x) and q(x) are constants. Since constants are continuous everywhere, the conditions for existence and uniqueness are always met for any interval. This means that for any given initial conditions, a unique solution always exists for these equations.

3. The Wronskian

The Wronskian is a determinant that helps us determine whether a set of solutions to a homogeneous linear differential equation is linearly independent.

Consider a second-order homogeneous linear differential equation:

ay'' + by' + cy = 0

Suppose we have two solutions, y1(x) and y2(x), to this equation.

The Wronskian of y1 and y2, denoted by W(y1, y2)(x), is defined as the determinant:

W(y1, y2)(x) = | y1(x) y2(x) |

| y'1(x) y'2(x) |

This expands to:

W(y1, y2)(x) = y1(x) y'2(x) - y2(x) y'1(x)

For a third-order equation with solutions y1, y2, y3, the Wronskian is:

W(y1, y2, y3)(x) = | y1 y2 y3 |

| y'1 y'2 y'3 |

| y''1 y''2 y''3 |

The Wronskian is a powerful tool for checking linear independence.

4. Independence of Solutions

Linear independence is a crucial concept when constructing the general solution to a homogeneous linear differential equation. A set of solutions {y1, y2, ..., yn} to an n-th order homogeneous linear differential equation is said to be linearly independent if the only solution to the equation:

c1 y1(x) + c2 y2(x) + ... + cn yn(x) = 0

for all x in the interval of interest is c1 = c2 = ... = cn = 0.

The relationship between the Wronskian and linear independence is as follows:

Theorem: If y1(x) and y2(x) are two solutions to the second-order homogeneous linear differential equation y'' + p(x) y' + q(x) y = 0 on an interval I, where p(x) and q(x) are continuous on I, then the following are equivalent:

1. y1 and y2 are linearly independent on I.
2. W(y1, y2)(x) ≠ 0 for all x in I.
3. W(y1, y2)(x) ≠ 0 for at least one x in I.

In simpler terms, if the Wronskian is non-zero at even one point in the interval, the solutions are linearly independent across the entire interval. For linear differential equations with constant coefficients, the derivatives of exponential functions (which are the typical solutions) are also exponential functions, and their Wronskians are generally non-zero unless the roots of the characteristic equation are repeated.

If y1(x) and y2(x) are linearly independent solutions to the homogeneous equation ay'' + by' + cy = 0, then the general solution is given by:

y(x) = c1 y1(x) + c2 y2(x)

where c1 and c2 are arbitrary constants. These linearly independent solutions y1 and y2 are often referred to as a fundamental set of solutions.

5. Solving Homogeneous Linear Differential Equations with Constant Coefficients

To find solutions to ay'' + by' + cy = 0, we assume a solution of the form y = erx.

Taking derivatives:

y' = r erx

y'' = r2 erx

Substituting into the differential equation:

a(r2 erx) + b(r erx) + c(erx) = 0

Since erx is never zero, we can divide by it:

ar2 + br + c = 0

This is called the characteristic equation (or auxiliary equation). The roots of this quadratic equation determine the form of the solutions.

Case 1: Distinct Real Roots

If the characteristic equation has two distinct real roots, r1 and r2 (i.e., b2 - 4ac > 0), then two linearly independent solutions are y1 = er1x and y2 = er2x.

The general solution is:

y(x) = c1 er1x + c2 er2x

Example: Solve y'' - 5y' + 6y = 0. The characteristic equation is r2 - 5r + 6 = 0. Factoring gives (r - 2)(r - 3) = 0. The roots are r1 = 2 and r2 = 3. These are distinct real roots. The general solution is y(x) = c1 e2x + c2 e3x.

Case 2: Repeated Real Roots

If the characteristic equation has one repeated real root, r (i.e., b2 - 4ac = 0), then one solution is y1 = erx. To find a second linearly independent solution, we use the method of reduction of order, which leads to y2 = x erx.

The general solution is:

y(x) = c1 erx + c2 x erx = (c1 + c2 x) erx

Example: Solve y'' + 4y' + 4y = 0. The characteristic equation is r2 + 4r + 4 = 0. Factoring gives (r + 2)2 = 0. The root is r = -2 (repeated). The general solution is y(x) = c1 e-2x + c2 x e-2x = (c1 + c2 x) e-2x.

Case 3: Complex Conjugate Roots

If the characteristic equation has complex conjugate roots, r = α ± iβ, where β ≠ 0 (i.e., b2 - 4ac < 0), then the solutions are of the form e(α + iβ)x and e(α - iβ)x.

Using Euler's formula, eiθ = cos(θ) + i sin(θ), we can rewrite these solutions in terms of real-valued functions:

e(α + iβ)x = eαx eiβx = eαx (cos(βx) + i sin(βx))

e(α - iβ)x = eαx e-iβx = eαx (cos(βx) - i sin(βx))

Linear combinations of these two solutions give two real, linearly independent solutions: y1 = eαx cos(βx) and y2 = eαx sin(βx).

The general solution is:

y(x) = c1 eαx cos(βx) + c2 eαx sin(βx) = eαx (c1 cos(βx) + c2 sin(βx))

Example: Solve y'' + 2y' + 5y = 0. The characteristic equation is r2 + 2r + 5 = 0. Using the quadratic formula r = [-b ± sqrt(b2 - 4ac)] / 2a: r = [-2 ± sqrt(22 - 4 * 1 * 5)] / (2 * 1) r = [-2 ± sqrt(4 - 20)] / 2 r = [-2 ± sqrt(-16)] / 2 r = [-2 ± 4i] / 2 r = -1 ± 2i Here, α = -1 and β = 2. The general solution is y(x) = e-x (c1 cos(2x) + c2 sin(2x)).

6. Initial Value Problems (IVPs) for Second-Order Equations

An initial value problem for a second-order linear differential equation requires finding a particular solution that satisfies not only the differential equation but also specific values for the function and its first derivative at a given point.

Consider the equation:

ay'' + by' + cy = g(x)

with initial conditions:

y(x0) = y0

y'(x0) = y'0

To solve such a problem, we first find the general solution to the differential equation.

If the equation is homogeneous (g(x) = 0), the general solution is y(x) = c1 y1(x) + c2 y2(x), where y1 and y2 form a fundamental set of solutions.

We then use the initial conditions to find the specific values of the constants c1 and c2.

Steps:

1. Find the general solution y(x) = c1 y1(x) + c2 y2(x).
2. Differentiate the general solution to find y'(x).
3. Substitute x0 and y0 into the general solution: y0 = c1 y1(x0) + c2 y2(x0).
4. Substitute x0 and y'0 into the derivative of the general solution: y'0 = c1 y'1(x0) + c2 y'2(x0).
5. Solve the system of two linear equations for c1 and c2.
6. Substitute the values of c1 and c2 back into the general solution to obtain the particular solution.

Example: Solve the IVP y'' + 4y = 0, with y(0) = 1 and y'(0) = -2.

1. Find the general solution: The characteristic equation is r2 + 4 = 0. The roots are r2 = -4, so r = ±2i. This is a complex case with α = 0 and β = 2. The general solution is y(x) = e0x (c1 cos(2x) + c2 sin(2x)), which simplifies to y(x) = c1 cos(2x) + c2 sin(2x).

2. Differentiate the general solution: y'(x) = -2c1 sin(2x) + 2c2 cos(2x).

3. Apply the first initial condition y(0) = 1: 1 = c1 cos(0) + c2 sin(0) 1 = c1 * 1 + c2 * 0 c1 = 1.

4. Apply the second initial condition y'(0) = -2: -2 = -2c1 sin(0) + 2c2 cos(0) -2 = -2c1 * 0 + 2c2 * 1 -2 = 2c2 c2 = -1.

5. Substitute c1 and c2 back into the general solution: y(x) = 1 * cos(2x) + (-1) * sin(2x) y(x) = cos(2x) - sin(2x).

This is the unique solution to the IVP.

7. Non-homogeneous Equations (Brief Overview)

While the focus of this topic is on homogeneous equations and the fundamental properties, it's important to note that the general solution to a non-homogeneous equation ay'' + by' + cy = g(x) is given by:

y(x) = yc(x) + yp(x)

where yc(x) is the complementary solution (the general solution to the associated homogeneous equation) and yp(x) is a particular solution to the non-homogeneous equation. Methods like Undetermined Coefficients and Variation of Parameters are used to find yp(x).

8. Higher-Order Equations

The concepts discussed – the characteristic equation, its roots, and the Wronskian – extend to linear differential equations of order higher than two. For an n-th order equation, the characteristic equation will be a polynomial of degree n, and finding its roots will yield n solutions (considering multiplicity and complex conjugates). The Wronskian will be an n x n determinant. The general solution will involve n arbitrary constants, and an IVP will require n initial conditions (y(x₀), y'(x₀), ..., y^(n-1)(x₀)).

Summary of Key Concepts