Exercise 7.  Solve the nonhomogeneous difference equations.  

Exercise 7 (c).  [Graphics:Images/ZTransformDEModHome_gr_1490.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_1491.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_1492.gif].  

Solution 7 (c).

See text and/or instructor's solution manual.

Answer.   [Graphics:../Images/ZTransformDEModHome_gr_1688.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_1689.gif]:  

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

and then get    

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

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

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

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

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

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

Therefore,  

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

We are done.   

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

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

                    and

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

At the conjugate pole we can use the computation

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

Thus,  

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

Therefore,  

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

We are done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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


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

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1729.gif]
                                                            
                                                            
[Graphics:../Images/ZTransformDEModHome_gr_1730.gif]  

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


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

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

We are really done.   

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

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

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

We are really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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

Aside.  The Maple command is similar  

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

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

We are really really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_1740.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1741.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1742.gif]

  

          [Graphics:../Images/ZTransformDEModHome_gr_1743.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1744.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1745.gif]

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

 

We are really really really really done.    

The Details for the Partial Fractions.   

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

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

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

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

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

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_1751.gif]  

Now we have the desired partial fraction form:

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

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

                    [Graphics:../Images/ZTransformDEModHome_gr_1754.gif][Graphics:../Images/ZTransformDEModHome_gr_1755.gif].  

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 



This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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