Difference between revisions of "Method of Undetermined Coefficients"

From Department of Mathematics at UTSA
Jump to navigation Jump to search
Line 158: Line 158:
  
 
:<math>y= t^2 - 2 t + 2 + c_1 e^{-t}</math>
 
:<math>y= t^2 - 2 t + 2 + c_1 e^{-t}</math>
 +
 +
==Higher Order==
 +
<p>Recall from <a href="/the-method-of-undetermined-coefficients">The Method of Undetermined Coefficients</a> page that if we had a second order linear nonhomogenous differential equation whose coefficients were constant, that is <math>a \frac{d^2y}{dt^2} + b \frac{dy}{dt} + cy = g(t)</math>, then to solve this differential equation, all we need to do is solve the corresponding second order linear homogenous differential equation <math>a \frac{d^2y}{dt^2} + b \frac{dy}{dt} + cy = 0</math>, and then find a partial solution by assuming the form of the particular solution. More precisely, if <math>g(t)</math> was a polynomial, exponential, or sine/cosine function (or a combination of these), then we could assume a form for the particular solutions (see the linked page above for more details) and solve for the coefficients of this form to obtain a particular solution.</p>
 +
<p>We will now look at this method for higher order linear nonhomogenous differential equations. Consider the following <math>n^{\mathrm{th}}</math> order linear nonhomogenous differential equation with coefficients <math>a_0, a_1, ..., a_n \in \mathbb{R}</math>:</p>
 +
 +
<math>\begin{align} \quad a_0 \frac{d^ny}{dt^n} + a_1 \frac{d^{n-1}y}{dt^{n-1}} + ... + a_{n-1}\frac{dy}{dt} + a_n y = g(t) \end{align}</math>
 +
<p>Suppose that <math>g(t)</math> is of a form containing a polynomial, exponential function, or a sine/cosine function (like with when we were dealing with the method of undetermined coefficients for second order linear nonhomogenous differential equations). We can then find a particular solution <math>Y(t)</math> if we extend the assumed forms of combinations of polynomials, exponential functions, or sine/cosine functions, and multiplying by powers of <math>t</math> to ensure our particular solution does not contain part of a solution to the corresponding higher order linear homogenous differential equation. We can then solve for these constants by plugging our assumed form <math>Y(t)</math> into our differential equation above.</p>
 +
<p>Providing a list of possible forms is trivial, so we will instead look at some examples of apply this method.</p>
 +
<h2 id="toc1">Example 1</h2>
 +
<p><strong>Find the general solution to the differential equation <math>\frac{d^3y}{dt^3} + \frac{d^2y}{dt^2} + \frac{dy}{dt} + y = e^{-t} + 4t</math>.</strong></p>
 +
<p>We will need to first solve the corresponding third order linear homogenous differential equation <math>\frac{d^3y}{dt^3} + \frac{d^2y}{dt^2} + \frac{dy}{dt} + y = 0</math>. This characteristic equation to this differential equation is:</p>
 +
 +
<math>\begin{align} \quad r^3 + r^2 + r + 1 = 0 \end{align}</math>
 +
<p>We can (by trial and error) see that <math>r_1 = -1</math> is a solution to this characteristic equation. Applying long division with the factor <math>(r + 1)</math> and we have that our characteristic equation can be written as:</p>
 +
 +
<math>\begin{align} \quad (r + 1)(r^2 + 1) = 0 \end{align}</math>
 +
<p>Therefore we can see that <math>r_2 = i</math> and <math>r_3 = -i</math>. Therefore the general solution to our third order linear homogenous differential equation is:</p>
 +
 +
<math>\begin{align} \quad y_h(t) = C_1e^{-t} + C_2\cos(t) + C_3\sin(t) \end{align}</math>
 +
<p>Now we note that <math>g(t) = e^{-t} + 4t</math> has an exponential term and a cosine term, so we expect the form of our particular solution to be <math>Y(t) = Ae^{-t} ...</math>. We now compute the first, second, and third derivatives of <math>Y</math>.</p>
 +
  
 
==Licensing==
 
==Licensing==
 
Content obtained and/or adapted from:
 
Content obtained and/or adapted from:
 
* [https://en.wikipedia.org/wiki/Method_of_undetermined_coefficients Method of undetermined coefficients, Wikipedia] under a CC BY-SA license
 
* [https://en.wikipedia.org/wiki/Method_of_undetermined_coefficients Method of undetermined coefficients, Wikipedia] under a CC BY-SA license

Revision as of 14:49, 26 October 2021

In mathematics, the method of undetermined coefficients is an approach to finding a particular solution to certain nonhomogeneous ordinary differential equations and recurrence relations. It is closely related to the annihilator method, but instead of using a particular kind of differential operator (the annihilator) in order to find the best possible form of the particular solution, a "guess" is made as to the appropriate form, which is then tested by differentiating the resulting equation. For complex equations, the annihilator method or variation of parameters is less time-consuming to perform.

Undetermined coefficients is not as general a method as variation of parameters, since it only works for differential equations that follow certain forms.[1]

Description of the method

Consider a linear non-homogeneous ordinary differential equation of the form

where denotes the i-th derivative of , and denotes a function of .

The method of undetermined coefficients provides a straightforward method of obtaining the solution to this ODE when two criteria are met:[2]

  1. are constants.
  2. g(x) is a constant, a polynomial function, exponential function , sine or cosine functions or , or finite sums and products of these functions (, constants).

The method consists of finding the general homogeneous solution for the complementary linear homogeneous differential equation

and a particular integral of the linear non-homogeneous ordinary differential equation based on . Then the general solution to the linear non-homogeneous ordinary differential equation would be

[3]

If consists of the sum of two functions and we say that is the solution based on and the solution based on . Then, using a superposition principle, we can say that the particular integral is[3]

Typical forms of the particular integral

In order to find the particular integral, we need to 'guess' its form, with some coefficients left as variables to be solved for. This takes the form of the first derivative of the complementary function. Below is a table of some typical functions and the solution to guess for them.

Function of x Form for y

If a term in the above particular integral for y appears in the homogeneous solution, it is necessary to multiply by a sufficiently large power of x in order to make the solution independent. If the function of x is a sum of terms in the above table, the particular integral can be guessed using a sum of the corresponding terms for y.[1]

Examples

Example 1

Find a particular integral of the equation

The right side t cos t has the form

with n = 2, α = 0, and β = 1.

Since α + = i is a simple root of the characteristic equation

we should try a particular integral of the form

Substituting yp into the differential equation, we have the identity

Comparing both sides, we have

which has the solution

We then have a particular integral

Example 2

Consider the following linear nonhomogeneous differential equation:

This is like the first example above, except that the nonhomogeneous part () is not linearly independent to the general solution of the homogeneous part (); as a result, we have to multiply our guess by a sufficiently large power of x to make it linearly independent.

Here our guess becomes:

By substituting this function and its derivative into the differential equation, one can solve for A:

So, the general solution to this differential equation is:

Example 3

Find the general solution of the equation:

is a polynomial of degree 2, so we look for a solution using the same form,

Plugging this particular function into the original equation yields,

which gives:

Solving for constants we get:

To solve for the general solution,

where is the homogeneous solution , therefore, the general solution is:

Higher Order

Recall from <a href="/the-method-of-undetermined-coefficients">The Method of Undetermined Coefficients</a> page that if we had a second order linear nonhomogenous differential equation whose coefficients were constant, that is , then to solve this differential equation, all we need to do is solve the corresponding second order linear homogenous differential equation , and then find a partial solution by assuming the form of the particular solution. More precisely, if was a polynomial, exponential, or sine/cosine function (or a combination of these), then we could assume a form for the particular solutions (see the linked page above for more details) and solve for the coefficients of this form to obtain a particular solution.

We will now look at this method for higher order linear nonhomogenous differential equations. Consider the following order linear nonhomogenous differential equation with coefficients :

Suppose that is of a form containing a polynomial, exponential function, or a sine/cosine function (like with when we were dealing with the method of undetermined coefficients for second order linear nonhomogenous differential equations). We can then find a particular solution if we extend the assumed forms of combinations of polynomials, exponential functions, or sine/cosine functions, and multiplying by powers of to ensure our particular solution does not contain part of a solution to the corresponding higher order linear homogenous differential equation. We can then solve for these constants by plugging our assumed form into our differential equation above.

Providing a list of possible forms is trivial, so we will instead look at some examples of apply this method.

Example 1

Find the general solution to the differential equation .

We will need to first solve the corresponding third order linear homogenous differential equation . This characteristic equation to this differential equation is:

We can (by trial and error) see that is a solution to this characteristic equation. Applying long division with the factor and we have that our characteristic equation can be written as:

Therefore we can see that and . Therefore the general solution to our third order linear homogenous differential equation is:

Now we note that has an exponential term and a cosine term, so we expect the form of our particular solution to be . We now compute the first, second, and third derivatives of .


Licensing

Content obtained and/or adapted from:

  • 1.0 1.1 Ralph P. Grimaldi (2000). "Nonhomogeneous Recurrence Relations". Section 3.3.3 of Handbook of Discrete and Combinatorial Mathematics. Kenneth H. Rosen, ed. CRC Press. Template:Isbn.
  • Template:Cite book
  • 3.0 3.1 Template:Cite book