Exercise 4.  Solve the homogeneous difference equations.  

Exercise 4 (a).  [Graphics:Images/ZTransformDEModHome_gr_770.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_771.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_772.gif].  

Solution 4 (a).

See text and/or instructor's solution manual.

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

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

Solution.   

Method 1.
  The characteristic equation  

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

has roots [Graphics:../Images/ZTransformDEModHome_gr_781.gif] and [Graphics:../Images/ZTransformDEModHome_gr_782.gif].

The general solution is  

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

Solve the linear system  

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

and get  [Graphics:../Images/ZTransformDEModHome_gr_785.gif]  and  [Graphics:../Images/ZTransformDEModHome_gr_786.gif].

Therefore,  

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

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

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

and then get  

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

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

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

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

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

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

Therefore,  

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

We are done.   

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

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

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

        Since   [Graphics:../Images/ZTransformDEModHome_gr_798.gif] and   [Graphics:../Images/ZTransformDEModHome_gr_799.gif]we can write  

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

                    [Graphics:../Images/ZTransformDEModHome_gr_801.gif]   for   [Graphics:../Images/ZTransformDEModHome_gr_802.gif].  

Therefore, the solution has the following form:

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

We are done.   

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

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

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


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

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

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

 

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

[Graphics:../Images/ZTransformDEModHome_gr_811.gif]
[Graphics:../Images/ZTransformDEModHome_gr_812.gif]

We are really done.   

Using Residues.   Calculate the residues  [Graphics:../Images/ZTransformDEModHome_gr_813.gif]  at the poles   [Graphics:../Images/ZTransformDEModHome_gr_814.gif].  

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

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

Therefore,  

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

We are really really done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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

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

                                                            [Graphics:../Images/ZTransformDEModHome_gr_839.gif]
                                                            
                                                            
[Graphics:../Images/ZTransformDEModHome_gr_840.gif]  

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


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

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

We are really really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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

Aside.  The Maple command is similar  

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

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

We are really really really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_848.gif]     [Graphics:../Images/ZTransformDEModHome_gr_849.gif]     [Graphics:../Images/ZTransformDEModHome_gr_850.gif]

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

 

We are really really really really really done.   

The Details for the Partial Fractions.   

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

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

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

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

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

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

Then use the standard procedure for expanding in partial fractions   

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

Then make the substitution   [Graphics:../Images/ZTransformDEModHome_gr_871.gif]   in   [Graphics:../Images/ZTransformDEModHome_gr_872.gif]   and get  

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


Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (iii)  are

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

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


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

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


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

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


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

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

Now we have the desired partial fraction form:

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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