Annihilator method

From Wikipedia, the free encyclopedia

In mathematics, the annihilator method is a procedure used to find a particular solution to certain types of inhomogeneous ordinary differential equations. It is equivalent to the method of undetermined coefficients, and the two names are sometimes used interchangeably, although they can describe different techniques. The phrase undetermined coefficients can also be used to refer to the step in the annihilator method in which the coefficients are calculated.

The annihilator method is used as follows. Given the ODE P(D)y = f(x), find another differential operator A(D) such that A(D)f(x) = 0. This operator is called the annihilator, thus giving the method its name. Applying A(D) to both sides of the ODE gives a homogeneous ODE \big(A(D)P(D)\big)y = 0 for which we find a solution basis \{y_1,\ldots,y_n\} as before. Then the original inhomogeneous ODE is used to construct a system of equations restricting the coefficients of the linear combinations to satisfy the ODE.

This method is not as general as variation of parameters in the sense that an annihilator does not always exist.

[edit] Example

Given y'' − 4y' + 5y = sin(kx), P(D) = D2 − 4D + 5. The simplest annihilator of sin(kx) is A(D) = D2 + k2. The zeros of A(z)P(z) are {2 + i,2 − i,ik, − ik}, so the solution basis of A(D)P(D) is {y1,y2,y3,y4} = {e(2 + i)x,e(2 − i)x,eikx,e ikx}.

Setting y = c1y1 + c2y2 + c3y3 + c4y4 we find

sin(kx) = P(D)y
= P(D)(c1y1 + c2y + c3y3 + c4y4)
= c1P(D)y1 + c2P(D)y2 + c3P(D)y3 + c4P(D)y4
= 0 + 0 + c3( − k2 − 4ik + 5)y3 + c4( − k2 + 4ik + 5)y4
= c3( − k2 − 4ik + 5)(cos(kx) + isin(kx)) + c4( − k2 + 4ik + 5)(cos(kx) − isin(kx))

giving the system

i = (k2 + 4ik − 5)c3 + ( − k2 + 4ik + 5)c4
0 = (k2 + 4ik − 5)c3 + (k2 − 4ik − 5)c4

which has solutions

c_3=\frac i{2(k^2+4ik-5)}, c_4=\frac i{2(-k^2+4ik+5)}

giving the solution set

y\, =c_1y_1+c_2y_2+\frac i{2(k^2+4ik-5)}y_3+\frac i{2(-k^2+4ik+5)}y_4
=c_1y_1+c_2y_2+\frac{4k\cos(kx)-(k^2-5)\sin(kx)}
{(k^2+4ik-5)(k^2-4ik-5)}
=c_1y_1+c_2y_2
+\frac{4k\cos(kx)+(5-k^2)\sin(kx)}{k^4+6k^2+25}.