Exercise 18.  In the value of an annuity due model   [Graphics:Images/ZTransformIntroModHome_gr_1262.gif]  

use the parameters   [Graphics:Images/ZTransformIntroModHome_gr_1263.gif].

18 (a).  Use the z-transform and tables to find the solution.    Hint.  Get  [Graphics:Images/ZTransformIntroModHome_gr_1264.gif].  

18 (b).  Use residues to find the solution.

Solution 18.

See text and/or instructor's solution manual.

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

Solution 18 (a).

        Substitute   [Graphics:../Images/ZTransformIntroModHome_gr_1266.gif]   into   [Graphics:../Images/ZTransformIntroModHome_gr_1267.gif],  

and get the difference equation  

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

Take the z-transform of both sides and use the assumed initial condition  [Graphics:../Images/ZTransformIntroModHome_gr_1269.gif]:  

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

                    [Graphics:../Images/ZTransformIntroModHome_gr_1271.gif].  
        
Solve for  [Graphics:../Images/ZTransformIntroModHome_gr_1272.gif]  and get  

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

then use Table 9.1 to find the inverse z-transform   

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

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

 

We are done.   

 

Solution 18 (b).

        Using residues we get

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

                    and

                    [Graphics:../Images/ZTransformIntroModHome_gr_1276.gif]  
                    
Thus

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

Therefore,

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

 

We are done.   

Aside.  We can let Mathematica double check our work.

 

Take the z-transform of both sides.

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

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


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

[Graphics:../Images/ZTransformIntroModHome_gr_1282.gif]
[Graphics:../Images/ZTransformIntroModHome_gr_1283.gif]
[Graphics:../Images/ZTransformIntroModHome_gr_1284.gif]

Find the inverse z-transform.

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

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


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

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


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

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


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

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


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

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


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

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

          The Maple commands are similar  

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

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


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

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


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

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


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

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

 

Aside.  We can use Mathematica's InverseZTransform subroutine.

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

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


[Graphics:../Images/ZTransformIntroModHome_gr_1307.gif]
[Graphics:../Images/ZTransformIntroModHome_gr_1288.gif] [Graphics:../Images/ZTransformIntroModHome_gr_1308.gif] [Graphics:../Images/ZTransformIntroModHome_gr_1309.gif]

 

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

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

          The Maple command is similar  

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

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

 

We are really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

 

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

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

          The Maple command is similar  

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

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

Aside.  The solution can be written as  

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

 

We are really really done.   

Aside.  We can graph the solution.

 

          [Graphics:../Images/ZTransformIntroModHome_gr_1319.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1320.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1321.gif]

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

 

 

We are really really really done.   

The Details for the Partial Fractions.   

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

It is natural to try the command:

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

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

But this is not the desired form for using Table 9.1 of z-transforms.

 

Method (i).   Use the following algebra steps  

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

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

Then use the standard procedure for expanding in partial fractions   

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

Then make the substitution   [Graphics:../Images/ZTransformIntroModHome_gr_1341.gif]   in   [Graphics:../Images/ZTransformIntroModHome_gr_1342.gif]   and get  

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


Therefore, the desired form is  

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

Aside.   The Mathematica commands for Method (iii)  are

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

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


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

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


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

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


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

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


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

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

Now use the substitution

                    [Graphics:../Images/ZTransformIntroModHome_gr_1355.gif].     
                    
Therefore, the desired form is  

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

 

  

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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