Exercise 15.  Solve the difference equation   [Graphics:Images/ZTransformIntroModHome_gr_978.gif]   with initial condition   [Graphics:Images/ZTransformIntroModHome_gr_979.gif].

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

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

Solution 15.

See text and/or instructor's solution manual.

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

Solution 15 (a).  

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

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

                    [Graphics:../Images/ZTransformIntroModHome_gr_984.gif].  
        
Solve for  [Graphics:../Images/ZTransformIntroModHome_gr_985.gif]  and get  

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

then use Table 9.1 to find the inverse z-transform   

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

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

 

We are done.   

 

Solution 15 (b).  

        Using residues we get

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

                    and

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

Thus

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

Therefore,

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

 

We are done.   

Aside.  We can let Mathematica double check our work.

 

Take the z-transform of both sides.

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

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


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

[Graphics:../Images/ZTransformIntroModHome_gr_995.gif]
[Graphics:../Images/ZTransformIntroModHome_gr_996.gif]
[Graphics:../Images/ZTransformIntroModHome_gr_997.gif]


[Graphics:../Images/ZTransformIntroModHome_gr_998.gif]
[Graphics:../Images/ZTransformIntroModHome_gr_999.gif]
[Graphics:../Images/ZTransformIntroModHome_gr_1000.gif]

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

Find the inverse z-transform.

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

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


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

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


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

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


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

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


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

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


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

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

          The Maple commands are similar  

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

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


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

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


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

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


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

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

 

Aside.  We can use Mathematica's InverseZTransform subroutine.

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

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


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

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

 

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

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

 

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

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

          The Maple command is similar  

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

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

 

We are really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

 

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

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

          The Maple command is similar  

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

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

 

We are really really done.   

Aside.  We can graph the solution.

 

          [Graphics:../Images/ZTransformIntroModHome_gr_1036.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1037.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_1038.gif]

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

 

 

We are really really done.   

The Details for the Partial Fractions.   

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

It is natural to try the command:

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

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

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

 

Method (i).   Use the following algebra steps  

                      

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

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

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

and solve the linear system  

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

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

Therefore, the desired form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

 

  

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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