Second Order Difference EquationsSecond Order Difference Equations

Exercise 6.  Solve the nonhomogeneous difference equations.  

Exercise 6 (c).  [Graphics:Images/ZTransformDEModHome_gr_1210.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_1211.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_1212.gif].  

Solution 6 (c).

See text and/or instructor's solution manual.

Answer.   [Graphics:../Images/ZTransformDEModHome_gr_1408.gif].  

Remark.   The preferred method to use involves Z-transforms, limits and residues.  

Solution.   Take the z-transform of both sides and use the initial conditions   [Graphics:../Images/ZTransformDEModHome_gr_1409.gif]:  

                    [Graphics:../Images/ZTransformDEModHome_gr_1410.gif],  

and then get    

                    [Graphics:../Images/ZTransformDEModHome_gr_1411.gif].  

New solve for  [Graphics:../Images/ZTransformDEModHome_gr_1412.gif]  and obtain:    

                    [Graphics:../Images/ZTransformDEModHome_gr_1413.gif].  

Using Tables.   Using Table 9.1 and the formula   [Graphics:../Images/ZTransformDEModHome_gr_1414.gif]   we get:  

                    [Graphics:../Images/ZTransformDEModHome_gr_1415.gif]  

Remark.  The details for the partial fraction expansion are at the bottom of the page.

Therefore,  

                    [Graphics:../Images/ZTransformDEModHome_gr_1416.gif].  

We are done.   

Using Residues.   Calculate residues of  [Graphics:../Images/ZTransformDEModHome_gr_1417.gif]  at the poles   [Graphics:../Images/ZTransformDEModHome_gr_1418.gif].  

                    [Graphics:../Images/ZTransformDEModHome_gr_1419.gif]  

                    and

                    [Graphics:../Images/ZTransformDEModHome_gr_1420.gif]  

At the conjugate pole we can use the computation

                    [Graphics:../Images/ZTransformDEModHome_gr_1421.gif].  

Thus,  

                    [Graphics:../Images/ZTransformDEModHome_gr_1422.gif]  

Therefore,  

                    [Graphics:../Images/ZTransformDEModHome_gr_1423.gif].  

We are done.   

Aside.  We can let Mathematica double check our work.

[Graphics:../Images/ZTransformDEModHome_gr_1424.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1425.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1426.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1427.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1428.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1429.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1430.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1431.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1432.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1433.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1434.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1435.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1436.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1437.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1438.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1439.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1440.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1441.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1442.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1443.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1444.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1445.gif]

Aside.  The Maple commands are similar  

[Graphics:../Images/ZTransformDEModHome_gr_1446.gif]  

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1447.gif]

  
[Graphics:../Images/ZTransformDEModHome_gr_1448.gif]  

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1449.gif]
                                                                                                                                                                                          
                                                            
[Graphics:../Images/ZTransformDEModHome_gr_1450.gif]  

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1451.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1452.gif]  

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1453.gif]

We are really done.   

Aside.  The solution can be written in the alternative form:  

[Graphics:../Images/ZTransformDEModHome_gr_1454.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1455.gif]

We are really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

[Graphics:../Images/ZTransformDEModHome_gr_1456.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1457.gif]

Aside.  The Maple command is similar  

[Graphics:../Images/ZTransformDEModHome_gr_1458.gif]  

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1459.gif]

We are really really really done.   

 

Aside.  We can graph some of the terms in the sequence.

 

               [Graphics:../Images/ZTransformDEModHome_gr_1460.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1461.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1462.gif]

          [Graphics:../Images/ZTransformDEModHome_gr_1463.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1464.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1465.gif]

                    The sequence   [Graphics:../Images/ZTransformDEModHome_gr_1466.gif].  

 

We are really really really really done.   

The Details for the Partial Fractions.   

Aside.  How can we expand   [Graphics:../Images/ZTransformDEModHome_gr_1467.gif]   into the proper partial fractions?

It is natural to use the standard partial fraction expansion and the command:

[Graphics:../Images/ZTransformDEModHome_gr_1468.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1469.gif]

However, as we have seen, this will produce a solution involving the  [Graphics:../Images/ZTransformDEModHome_gr_1470.gif]  functions.   

This can be overcome if we use a special partial fraction expansion that is easier to use with Table 9.1.

Method (i).   Use the following algebra steps  

                    [Graphics:../Images/ZTransformDEModHome_gr_1471.gif]  

Now we have the desired partial fraction form:

                    [Graphics:../Images/ZTransformDEModHome_gr_1472.gif].

Method (ii).   Find the linear combination of   [Graphics:../Images/ZTransformDEModHome_gr_1473.gif],  

                    [Graphics:../Images/ZTransformDEModHome_gr_1474.gif].  

Equate the numerators   [Graphics:../Images/ZTransformDEModHome_gr_1475.gif],  

and solve the linear system  

                    [Graphics:../Images/ZTransformDEModHome_gr_1476.gif]  

and get   [Graphics:../Images/ZTransformDEModHome_gr_1477.gif].   

Therefore, the desired partial fraction form is  

                    [Graphics:../Images/ZTransformDEModHome_gr_1478.gif].  

Aside.   The Mathematica commands for Method (ii)  are

[Graphics:../Images/ZTransformDEModHome_gr_1479.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1480.gif]


[Graphics:../Images/ZTransformDEModHome_gr_1481.gif]

[Graphics:../Images/ZTransformDEModHome_gr_1482.gif]

Method (iii).   The substitution   [Graphics:../Images/ZTransformDEModHome_gr_1483.gif]   does not apply when there are complex roots.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

(c) 2008 John H. Mathews, Russell W. Howell