The approach to solving linear constant coefficient difference equations is to find the general form of all possible solutions to the equation and then apply a number of conditions to find the appropriate solution. The two main types of problems are initial value problems, which involve constraints on the solution at several consecutive points, and boundary value problems, which involve constraints on the solution at nonconsecutive points.
The number of initial conditions needed for an th order difference equation, which is the order of the highest order difference or the largest delay parameter of the output in the equation, is , and a unique solution is always guaranteed if these are supplied. Boundary value probelms can be slightly more complicated and will not necessarily have a unique solution or even a solution at all for a given set of conditions. Thus, this section will focus exclusively on initial value problems.
Consider some linear constant coefficient difference equation given by , in which is a difference operator of the form
where is the first difference operator
Let and be two functions such that and . By the linearity of , note that . Thus, the form of the general solution to any linear constant coefficient ordinary differential equation is the sum of a homogeneous solution to the equation and a particular solution that is specific to the forcing function .
We wish to determine the forms of the homogeneous and nonhomogeneous solutions in full generality in order to avoid incorrectly restricting the form of the solution before applying any conditions. Otherwise, a valid set of initial or boundary conditions might appear to have no corresponding solution trajectory. The following sections discuss how to accomplish this for linear constant coefficient difference equations.
In order to find the homogeneous solution to a difference equation described by the recurrence relation , consider the difference equation . We know that the solutions have the form for some complex constants . Since for a solution it follows that
so it also follows that
Therefore, the solution exponential are the roots of the above polynomial, called the characteristic polynomial.
For equations of order two or more, there will be several roots. If all of the roots are distinct, then the general form of the homogeneous solution is simply
If a root has multiplicity that is greater than one, the repeated solutions must be multiplied by each power of from 0 to one less than the root multipicity (in order to ensure linearly independent solutions). For instance, if had a multiplicity of 2 and had multiplicity 3, the homogeneous solution would be
Recall that the Fibonacci sequence describes a (very unrealistic) model of what happens when a pair rabbits get left alone in a black box... The assumptions are that a pair of rabits never die and produce a pair of offspring every month starting on their second month of life. This system is defined by the recursion relation for the number of rabit pairs at month
with the initial conditions and .
Note that the forcing function is zero, so only the homogenous solution is needed. It is easy to see that the characteristic polynomial is , so there are two roots with multiplicity one. These are and . Thus, the solution is of the form
Using the initial conditions, we determine that
Hence, the Fibonacci sequence is given by
Finding the particular solution is a slightly more complicated task than finding the homogeneous solution. It can be found through convolution of the input with the unit impulse response once the unit impulse response is known. Finding the particular solution ot a differential equation is discussed further in the chapter concerning the z-transform, which greatly simplifies the procedure for solving linear constant coefficient differential equations using frequency domain tools.
Consider the following difference equation describing a system with feedback
In order to find the homogeneous solution, consider the difference equation
It is easy to see that the characteristic polynomial is , so is the only root. Thus the homogeneous solution is of the form
In order to find the particular solution, consider the output for the unit impulse case
By inspection, it is clear that the impulse response is . Hence, the particular solution for a given is
Therefore, the general solution is
Initial conditions and a specific input can further tailor this solution to a specific situation.
Linear constant coefficient difference equations are useful for modeling a wide variety of discrete time systems. The approach to solving them is to find the general form of all possible solutions to the equation and then apply a number of conditions to find the appropriate solution. This is done by finding the homogeneous solution to the difference equation that does not depend on the forcing function input and a particular solution to the difference equation that does depend on the forcing function input.