Difference between revisions of "Reduction of the Order"

From Department of Mathematics at UTSA
Jump to navigation Jump to search
Line 1: Line 1:
'''Reduction of order''' is a technique in [[mathematics]] for solving second-order linear [[Ordinary differential equation|ordinary]] [[differential equation]]s.  It is employed when one solution <math>y_1(x)</math> is known and a second [[Linear independence|linearly independent]] solution <math>y_2(x)</math> is desired. The method also applies to ''n''-th order equations. In this case the [[ansatz]] will yield an (''n''−1)-th order equation for <math>v</math>.
+
'''Reduction of order''' is a technique in mathematics for solving second-order linear ordinary differential equations.  It is employed when one solution <math>y_1(x)</math> is known and a second linearly independent solution <math>y_2(x)</math> is desired. The method also applies to ''n''-th order equations. In this case the ansatz will yield an (''n''−1)-th order equation for <math>v</math>.
  
 
== Second-order linear ordinary differential equations==
 
== Second-order linear ordinary differential equations==
Line 9: Line 9:
 
:<math> a y''(x) + b y'(x) + c y(x) = 0,</math>
 
:<math> a y''(x) + b y'(x) + c y(x) = 0,</math>
  
where <math>a, b, c</math> are real non-zero coefficients. Two linearly independent solutions for this ODE can be straightforwardly found using [[characteristic equation (calculus)|characteristic equations]] except for the case when the [[discriminant]], <math>b^2 - 4 a c</math>, vanishes. In this case,
+
where <math>a, b, c</math> are real non-zero coefficients. Two linearly independent solutions for this ODE can be straightforwardly found using characteristic equations except for the case when the discriminant, <math>b^2 - 4 a c</math>, vanishes. In this case,
  
 
:<math> a y''(x) + b y'(x) + \frac{b^2}{4a} y(x) = 0,</math>
 
:<math> a y''(x) + b y'(x) + \frac{b^2}{4a} y(x) = 0,</math>
Line 39: Line 39:
 
:<math> a y_1 v'' = 0.</math>
 
:<math> a y_1 v'' = 0.</math>
  
Since <math>a</math> is assumed non-zero and <math>y_1(x)</math> is an [[exponential function]] (and thus always non-zero), we have
+
Since <math>a</math> is assumed non-zero and <math>y_1(x)</math> is an exponential function (and thus always non-zero), we have
  
 
:<math> v'' = 0.</math>
 
:<math> v'' = 0.</math>
Line 55: Line 55:
 
:<math> y_2(x) = x y_1(x) = x e^{-\frac{b}{2 a} x}.</math>
 
:<math> y_2(x) = x y_1(x) = x e^{-\frac{b}{2 a} x}.</math>
  
Finally, we can prove that the second solution <math>y_2(x)</math> found via this method is linearly independent of the first solution by calculating the [[Wronskian]]
+
Finally, we can prove that the second solution <math>y_2(x)</math> found via this method is linearly independent of the first solution by calculating the Wronskian
  
 
:<math>W(y_1,y_2)(x) = \begin{vmatrix} y_1 & x y_1 \\ y_1' & y_1 + x y_1' \end{vmatrix} = y_1 ( y_1 + x y_1' ) - x y_1 y_1' = y_1^{2} + x y_1 y_1' - x y_1 y_1' = y_1^{2} = e^{-\frac{b}{a}x} \neq 0.</math>
 
:<math>W(y_1,y_2)(x) = \begin{vmatrix} y_1 & x y_1 \\ y_1' & y_1 + x y_1' \end{vmatrix} = y_1 ( y_1 + x y_1' ) - x y_1 y_1' = y_1^{2} + x y_1 y_1' - x y_1 y_1' = y_1^{2} = e^{-\frac{b}{a}x} \neq 0.</math>
Line 91: Line 91:
 
:<math>v''+\left(\frac{2y_1'(t)}{y_1(t)}+p(t)\right)\,v'=\frac{r(t)}{y_1(t)}.</math>
 
:<math>v''+\left(\frac{2y_1'(t)}{y_1(t)}+p(t)\right)\,v'=\frac{r(t)}{y_1(t)}.</math>
  
The [[integrating factor]] is <math>\mu(t)=e^{\int(\frac{2y_1'(t)}{y_1(t)}+p(t))dt}=y_1^2(t)e^{\int p(t) dt}</math>.
+
The integrating factor is <math>\mu(t)=e^{\int(\frac{2y_1'(t)}{y_1(t)}+p(t))dt}=y_1^2(t)e^{\int p(t) dt}</math>.
  
 
Multiplying the differential equation by the integrating factor <math>\mu(t)</math>, the equation for <math>v(t)</math> can be reduced to
 
Multiplying the differential equation by the integrating factor <math>\mu(t)</math>, the equation for <math>v(t)</math> can be reduced to
Line 102: Line 102:
  
 
==Resources==
 
==Resources==
* [https://en.wikipedia.org/wiki/Reduction_of_order Reduction of order], Wikipedia
 
 
* [https://tutorial.math.lamar.edu/Classes/DE/ReductionofOrder.aspx Reduction of the Order Notes]. Produced by Paul Dawkins, Lamar University
 
* [https://tutorial.math.lamar.edu/Classes/DE/ReductionofOrder.aspx Reduction of the Order Notes]. Produced by Paul Dawkins, Lamar University
 +
 +
==References==
 +
# W. E. Boyce and R. C. DiPrima, Elementary Differential Equations and Boundary Value Problems (8th edition), John Wiley & Sons, Inc., 2005. ISBN 0-471-43338-1.
 +
# Teschl, Gerald (2012). Ordinary Differential Equations and Dynamical Systems. Providence: American Mathematical Society. ISBN 978-0-8218-8328-0.
 +
# Eric W. Weisstein, Second-Order Ordinary Differential Equation Second Solution, From MathWorld—A Wolfram Web Resource.

Revision as of 15:00, 20 October 2021

Reduction of order is a technique in mathematics for solving second-order linear ordinary differential equations. It is employed when one solution is known and a second linearly independent solution is desired. The method also applies to n-th order equations. In this case the ansatz will yield an (n−1)-th order equation for .

Second-order linear ordinary differential equations

An example

Consider the general, homogeneous, second-order linear constant coefficient ordinary differential equation. (ODE)

where are real non-zero coefficients. Two linearly independent solutions for this ODE can be straightforwardly found using characteristic equations except for the case when the discriminant, , vanishes. In this case,

from which only one solution,

can be found using its characteristic equation.

The method of reduction of order is used to obtain a second linearly independent solution to this differential equation using our one known solution. To find a second solution we take as a guess

where is an unknown function to be determined. Since must satisfy the original ODE, we substitute it back in to get

Rearranging this equation in terms of the derivatives of we get

Since we know that is a solution to the original problem, the coefficient of the last term is equal to zero. Furthermore, substituting into the second term's coefficient yields (for that coefficient)

Therefore, we are left with

Since is assumed non-zero and is an exponential function (and thus always non-zero), we have

This can be integrated twice to yield

where are constants of integration. We now can write our second solution as

Since the second term in is a scalar multiple of the first solution (and thus linearly dependent) we can drop that term, yielding a final solution of

Finally, we can prove that the second solution found via this method is linearly independent of the first solution by calculating the Wronskian

Thus is the second linearly independent solution we were looking for.

General method

Given the general non-homogeneous linear differential equation

and a single solution of the homogeneous equation [], let us try a solution of the full non-homogeneous equation in the form:

where is an arbitrary function. Thus

and

If these are substituted for , , and in the differential equation, then

Since is a solution of the original homogeneous differential equation, , so we can reduce to

which is a first-order differential equation for (reduction of order). Divide by , obtaining

The integrating factor is .

Multiplying the differential equation by the integrating factor , the equation for can be reduced to

After integrating the last equation, is found, containing one constant of integration. Then, integrate to find the full solution of the original non-homogeneous second-order equation, exhibiting two constants of integration as it should:


Resources

References

  1. W. E. Boyce and R. C. DiPrima, Elementary Differential Equations and Boundary Value Problems (8th edition), John Wiley & Sons, Inc., 2005. ISBN 0-471-43338-1.
  2. Teschl, Gerald (2012). Ordinary Differential Equations and Dynamical Systems. Providence: American Mathematical Society. ISBN 978-0-8218-8328-0.
  3. Eric W. Weisstein, Second-Order Ordinary Differential Equation Second Solution, From MathWorld—A Wolfram Web Resource.