Homogeneous linear differential equation with constant coefficients: Difference between revisions

From Calculus
(Created page with "==Definition== A '''homogeneous linear differential equation with constant coefficients''', which can also be thought of as a linear differential equation that is simulta...")
 
 
(3 intermediate revisions by the same user not shown)
Line 5: Line 5:
<math>y^{(k)} + p_{k-1}y^{(k-1)} + \dots + p_1y' + p_0y = 0</math>
<math>y^{(k)} + p_{k-1}y^{(k-1)} + \dots + p_1y' + p_0y = 0</math>


where <math>p_0,p_1,\dots,p_{k-1}</math> are all constants (i.e., real numbers).
where <math>p_0,p_1,\dots,p_{k-1}</math> are all constants (i.e., real numbers). Here <math>y</math> is the dependent variable and <math>x</math> (which does not appear explicitly above) is the independent variable with respect to which the differentiations occur.


===Solution method===
===Solution method===
Line 11: Line 11:
Consider the following polynomial:
Consider the following polynomial:


<math>t^k + p_{k-1}t^{k-1} + \dots + p_1t +p_0 = 0</math>
<math>t^k + p_{k-1}t^{k-1} + \dots + p_1t +p_0</math>


This polynomial is called the ''characteristic polynomial'' of the differential equation. We consider various cases:
This polynomial is called the ''characteristic polynomial'' of the differential equation. We consider various cases:
Line 24: Line 24:
| The polynomial splits completely over the complex numbers into distinct linear factors, but some of the roots are not real || For any real root <math>\alpha</math>, use <math>e^{\alpha x}</math> as a basis vector. Non-real roots occur in complex conjugate pairs. For a pair <math>a \pm ib</math>, choose the vectors <math>e^{ax}\cos(bx)</math> and <math>e^{ax}\sin(bx)</math>. Combining, we get a basis of <math>k</math> vectors.
| The polynomial splits completely over the complex numbers into distinct linear factors, but some of the roots are not real || For any real root <math>\alpha</math>, use <math>e^{\alpha x}</math> as a basis vector. Non-real roots occur in complex conjugate pairs. For a pair <math>a \pm ib</math>, choose the vectors <math>e^{ax}\cos(bx)</math> and <math>e^{ax}\sin(bx)</math>. Combining, we get a basis of <math>k</math> vectors.
|-
|-
| The general case || {{fillin}} -- combine ideas of two preceding generalizations.
| The general case || For a real root <math>\alpha</math> of multiplicity <math>s</math>, the <math>s</math> basis vectors are <math>x^de^{\alpha x}, 0 \le d < s</math>. For a pair of complex conjugates <math>a \pm ib</math> of multiplicity <math>s</math>, the <math>2s</math> basis vectors are <math>x^de^{ax}\cos(bx), 0 \le d < s</math> and <math>x^de^{ax}\sin(bx), 0 \le d < s</math>.
|}
|}
==Examples==
===Examples with distinct real roots===
Consider the differential equation:
<math>y'' - 3y' + 2y = 0</math>
The characteristic polynomial is:
<math>t^2 - 3t + 2</math>
We want to find the roots of this polynomial. The roots are <math>\alpha_1 = 1, \alpha_2 = 2</math>. Thus, the general solution is:
<math>y = C_1e^x + C_2e^{2x}, \qquad C_1,C_2 \in \R</math>
Here is another example:
<math>y''' + 2y'' -2y' - 4y = 0</math>
The characteristic polynomial is:
<math>t^3 + 2t^2 - 2t - 4</math>
The roots are <math>\alpha_1 = -2, \alpha_2 = -\sqrt{2}, \alpha_3 = \sqrt{2}</math>. The general solution is thus:
<math>y = C_1e^{-2x} + C_2e^{-\sqrt{2} x} + C_3e^{\sqrt{2} x}, \qquad C_1,C_2,C_3 \in \R</math>
===Examples with repeated real roots===
Consider the differential equation:
<math>y'' = 6y' - 9y</math>
First, we bring everything to one side, to get:
<math>y'' - 6y' + 9y = 0</math>
The characteristic polynomial is:
<math>t^2 - 6t + 9</math>
This factors as <math>(t - 3)^2</math>, so it has a root 3 with multiplicity 2. The general solution is thus:
<math>y = (C_1 + C_2x)e^{3x}, \qquad C_1,C_2 \in \R</math>
===Examples with real and complex roots without repetitions===
Consider the differential equation:
<math>y''' = y</math>
First, we bring everything to one side, to get:
<math>y''' - y = 0</math>
The characteristic polynomial is:
<math>t^3 - 1</math>
The only real root is 1, corresponding to a solution basis element <math>e^x</math>. The complex roots are the roots of <math>t^2 + t + 1</math>, and they are complex conjugates:
<math>\frac{-1 \pm \sqrt{3}i}{2}</math>
The basis solutions for this pair are thus <math>e^{-x/2}\cos(\sqrt{3} x/2)</math> and <math>e^{-x/2}\sin(\sqrt{3} x/2)</math>. The general solution is thus:
<math>y = C_1e^x + C_2e^{-x/2}\cos(\sqrt{3} x/2) + C_3e^{-x/2}\sin(\sqrt{3}x/2)</math>

Latest revision as of 00:03, 7 July 2012

Definition

A homogeneous linear differential equation with constant coefficients, which can also be thought of as a linear differential equation that is simultaneously an autonomous differential equation, is a differential equation of the form:

y(k)+pk1y(k1)++p1y+p0y=0

where p0,p1,,pk1 are all constants (i.e., real numbers). Here y is the dependent variable and x (which does not appear explicitly above) is the independent variable with respect to which the differentiations occur.

Solution method

Consider the following polynomial:

tk+pk1tk1++p1t+p0

This polynomial is called the characteristic polynomial of the differential equation. We consider various cases:

Case Solution in that case
The polynomial has pairwise distinct real roots α1,α2,,αk The solution space has basis eα1x,eα2x,,eαkx. In other words, the general solution is C1eα1x+C2eα2x++Ckeαkx where C1,C2,,Ck are freely varying real parameters.
The polynomial splits completely into linear factors over the reals, but with possible repetitions. α1 occurs s1 times, α2 occurs s2 times, and so on till αr, which occurs sr times. We have s1+s2++sr=k. The solution space has basis all functions of the form xdeαix where 0d<si with d an integer. Thus, for each i, there are si basis vectors corresponding to αi. We get a total of k basis vectors.
The polynomial splits completely over the complex numbers into distinct linear factors, but some of the roots are not real For any real root α, use eαx as a basis vector. Non-real roots occur in complex conjugate pairs. For a pair a±ib, choose the vectors eaxcos(bx) and eaxsin(bx). Combining, we get a basis of k vectors.
The general case For a real root α of multiplicity s, the s basis vectors are xdeαx,0d<s. For a pair of complex conjugates a±ib of multiplicity s, the 2s basis vectors are xdeaxcos(bx),0d<s and xdeaxsin(bx),0d<s.

Examples

Examples with distinct real roots

Consider the differential equation:

y3y+2y=0

The characteristic polynomial is:

t23t+2

We want to find the roots of this polynomial. The roots are α1=1,α2=2. Thus, the general solution is:

y=C1ex+C2e2x,C1,C2R

Here is another example:

y+2y2y4y=0

The characteristic polynomial is:

t3+2t22t4

The roots are α1=2,α2=2,α3=2. The general solution is thus:

y=C1e2x+C2e2x+C3e2x,C1,C2,C3R

Examples with repeated real roots

Consider the differential equation:

y=6y9y

First, we bring everything to one side, to get:

y6y+9y=0

The characteristic polynomial is:

t26t+9

This factors as (t3)2, so it has a root 3 with multiplicity 2. The general solution is thus:

y=(C1+C2x)e3x,C1,C2R

Examples with real and complex roots without repetitions

Consider the differential equation:

y=y

First, we bring everything to one side, to get:

yy=0

The characteristic polynomial is:

t31

The only real root is 1, corresponding to a solution basis element ex. The complex roots are the roots of t2+t+1, and they are complex conjugates:

1±3i2

The basis solutions for this pair are thus ex/2cos(3x/2) and ex/2sin(3x/2). The general solution is thus:

y=C1ex+C2ex/2cos(3x/2)+C3ex/2sin(3x/2)