Exercise 7.  Solve the nonhomogeneous difference equations.  

Exercise 7 (a).  [Graphics:Images/ZTransformDEModHome_gr_1484.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_1485.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_1486.gif].  

Solution 7 (a).

See text and/or instructor's solution manual.

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

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

and then get    

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

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

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

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

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

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

Therefore,  

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

We are done.   

Alternative Solution Using Tables.   Using ordinary partial fractions and Table 9.1  

and the formula   
[Graphics:../Images/ZTransformDEModHome_gr_1502.gif]   we get:  

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

        Since   [Graphics:../Images/ZTransformDEModHome_gr_1504.gif] and   [Graphics:../Images/ZTransformDEModHome_gr_1505.gif]we can write  

                    [Graphics:../Images/ZTransformDEModHome_gr_1506.gif]   and   

                    [Graphics:../Images/ZTransformDEModHome_gr_1507.gif]   for   [Graphics:../Images/ZTransformDEModHome_gr_1508.gif].  

Therefore, the solution has the following form:

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

We are done.   

Aside.  The commands for the ordinary partial fraction expansion are:  

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

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


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

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

        Here   [Graphics:../Images/ZTransformDEModHome_gr_1514.gif]   and   [Graphics:../Images/ZTransformDEModHome_gr_1515.gif],   and we can corroborate this solution.

 

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

[Graphics:../Images/ZTransformDEModHome_gr_1517.gif]
[Graphics:../Images/ZTransformDEModHome_gr_1518.gif]

We are really done.   

Using Residues.   Calculate the residues  [Graphics:../Images/ZTransformDEModHome_gr_1519.gif]  at the poles   [Graphics:../Images/ZTransformDEModHome_gr_1520.gif].  

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

                    and

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

                    and

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

Therefore,  

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

Therefore,  

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

We are really really done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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


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

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


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

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1553.gif]
                                                                                                      
                                                            
[Graphics:../Images/ZTransformDEModHome_gr_1554.gif]  

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


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

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

We are really really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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

Aside.  The Maple command is similar  

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

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

We are really really really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_1562.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1563.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1564.gif]

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

 

We are really really really really really done.   

The Details for the Partial Fractions.   

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

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

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

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

However, as we have seen, this will produce a solution involving the  [Graphics:../Images/ZTransformDEModHome_gr_1569.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_1570.gif]  

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

Method (iii).  (For distinct real roots)   First make the substitution   [Graphics:../Images/ZTransformDEModHome_gr_1581.gif]   in   [Graphics:../Images/ZTransformDEModHome_gr_1582.gif]   and get  

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

Then use the standard procedure for expanding in partial fractions   

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

Then make the substitution   [Graphics:../Images/ZTransformDEModHome_gr_1585.gif]   in   [Graphics:../Images/ZTransformDEModHome_gr_1586.gif]   and get  

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (iii)  are

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

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


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

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


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

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


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

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

Now we have the desired partial fraction form:

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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