Exercise 20.  Consider the system of difference equations

                    [Graphics:Images/ZTransformIntroModHome_gr_1468.gif]

with the initial conditions   [Graphics:Images/ZTransformIntroModHome_gr_1469.gif].  

20 (b).  Use z-transforms and residues to construct the solution in part (a).

Solution 20 (b).

Method (i)   Take the z-transform of both sides and use the initial conditions   [Graphics:../Images/ZTransformIntroModHome_gr_1493.gif].  

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

                    and    

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

then get  

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

                    and    

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

which can be written as    

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

                    and    

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

Solve for [Graphics:../Images/ZTransformIntroModHome_gr_1500.gif] and [Graphics:../Images/ZTransformIntroModHome_gr_1501.gif] and get

                    [Graphics:../Images/ZTransformIntroModHome_gr_1502.gif]   
                    
Then get

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

                    and    

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

Use Table 9.1 to find the inverse z-transform   

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

                    and    

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

 

We are done.   

 

Method (ii)   Using residues we obtain  [Graphics:../Images/ZTransformIntroModHome_gr_1507.gif]  

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

                    and

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

Thus,

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

Using residues we obtain  [Graphics:../Images/ZTransformIntroModHome_gr_1511.gif]  

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

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

Thus,

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

Therefore,   

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

                    and   

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

 

We are done.   

Aside.  We can let Mathematica double check our work.

 

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

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

Using limits and residues we obtain  [Graphics:../Images/ZTransformIntroModHome_gr_1519.gif]  

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

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


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

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


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

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

Using limits and residues we obtain  [Graphics:../Images/ZTransformIntroModHome_gr_1526.gif]  

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

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


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

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


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

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

          The Maple commands are similar  


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

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


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

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

Using limits and residues we obtain  [Graphics:../Images/ZTransformIntroModHome_gr_1537.gif]  

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

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


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

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

Using limits and residues we obtain  [Graphics:../Images/ZTransformIntroModHome_gr_1542.gif]  

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

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


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

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

 

Aside.  We can use Mathematica's InverseZTransform subroutine.

The following results were obtained using Mathematica 7.0

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

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


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

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


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

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


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

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

 

          The Maple commands are similar.  However, the inverse Z-transforms are complicated.

Hence, we will verify that the Z-transforms of   [Graphics:../Images/ZTransformIntroModHome_gr_1555.gif]   and   [Graphics:../Images/ZTransformIntroModHome_gr_1556.gif],  

are   
[Graphics:../Images/ZTransformIntroModHome_gr_1557.gif]   and   [Graphics:../Images/ZTransformIntroModHome_gr_1558.gif],   respectively.


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

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


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

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

 

We are really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

The following results were obtained using Mathematica 5.0

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

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

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

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

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

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

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

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

 

          The Maple command is similar 

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

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

 

We are really really really done.   

Aside.  We can graph the solution.

 

          [Graphics:../Images/ZTransformIntroModHome_gr_1573.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1574.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1575.gif]

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

          [Graphics:../Images/ZTransformIntroModHome_gr_1577.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1578.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1579.gif]

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

          [Graphics:../Images/ZTransformIntroModHome_gr_1581.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1582.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1583.gif]

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

 

  

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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