Exercise 1.  Solve the homogeneous difference equations.  

Exercise 1 (a).  [Graphics:Images/ZTransformDEModHome_gr_1.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_2.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_3.gif].  

Solution 1 (a).

See text and/or instructor's solution manual.

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

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

Solution.   

Method 1.
  The characteristic equation [Graphics:../Images/ZTransformDEModHome_gr_11.gif] has roots  [Graphics:../Images/ZTransformDEModHome_gr_12.gif] and [Graphics:../Images/ZTransformDEModHome_gr_13.gif].  

The general solution is  

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

Solve the linear system  

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

and get  [Graphics:../Images/ZTransformDEModHome_gr_16.gif]  and  [Graphics:../Images/ZTransformDEModHome_gr_17.gif].

Therefore,  

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

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

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

and then get  

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

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

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

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

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

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

Therefore,  

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

We are done.   

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

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

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

        Since   [Graphics:../Images/ZTransformDEModHome_gr_29.gif] and   [Graphics:../Images/ZTransformDEModHome_gr_30.gif]we can write  

                    [Graphics:../Images/ZTransformDEModHome_gr_31.gif]   and   [Graphics:../Images/ZTransformDEModHome_gr_32.gif]  for  [Graphics:../Images/ZTransformDEModHome_gr_33.gif].  

Therefore, the solution has the following form:

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

We are done.   

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

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

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


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

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

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

 

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

[Graphics:../Images/ZTransformDEModHome_gr_42.gif]
[Graphics:../Images/ZTransformDEModHome_gr_43.gif]

We are really done.   

Using Residues.   Calculate residues of  [Graphics:../Images/ZTransformDEModHome_gr_44.gif]  at the poles  [Graphics:../Images/ZTransformDEModHome_gr_45.gif]  and  [Graphics:../Images/ZTransformDEModHome_gr_46.gif].

                    [Graphics:../Images/ZTransformDEModHome_gr_47.gif]  
and
                    [Graphics:../Images/ZTransformDEModHome_gr_48.gif]  

Therefore,  

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

We are really really done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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


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

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


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

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

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

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


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

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

We are really really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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

Aside.  The Maple command is similar  

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

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

We are really really really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_82.gif]     [Graphics:../Images/ZTransformDEModHome_gr_83.gif]     [Graphics:../Images/ZTransformDEModHome_gr_84.gif]

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

 

We are really really really really really done.   

The Details for the Partial Fractions.   

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

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

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

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

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

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

Then use the standard procedure for expanding in partial fractions   

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

Then make the substitution   [Graphics:../Images/ZTransformDEModHome_gr_105.gif]   in   [Graphics:../Images/ZTransformDEModHome_gr_106.gif]   and get  

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


Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (iii)  are

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

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


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

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


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

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


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

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

Now we have the desired partial fraction form:

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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