Exercise 4.  Solve the homogeneous difference equations.  

Exercise 4 (c).  [Graphics:Images/ZTransformDEModHome_gr_776.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_777.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_778.gif].  

Solution 4 (c).

See text and/or instructor's solution manual.

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

Remark.   It is our goal to learn Z-transforms and use residues.  

Solution.   

Method 1.
  The characteristic equation  

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

has complex roots [Graphics:../Images/ZTransformDEModHome_gr_977.gif].   

The general solution is  

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

Solve the linear system  

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

and get  [Graphics:../Images/ZTransformDEModHome_gr_980.gif]  and  [Graphics:../Images/ZTransformDEModHome_gr_981.gif].  

Therefore,  

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

Method 2.  Take the z-transform of both sides and use the initial conditions   [Graphics:../Images/ZTransformDEModHome_gr_983.gif]:  

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

and then get    

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

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

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

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

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

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

Therefore,  

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

We are done.   

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

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

At the conjugate pole we can use the computation

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

Thus,  

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

Therefore,  

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

We are done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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


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

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


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

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


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

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

We are really done.   

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

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

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

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

We are really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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


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

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

Aside.  The Maple command is similar  

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

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

We are really really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_1032.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1033.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1034.gif]

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

 

We are really really really really done.    

The Details for the Partial Fractions.   

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

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

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

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

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

Now we have the desired partial fraction form:

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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