Exercise 14.  Solve the difference equation   [Graphics:Images/ZTransformIntroModHome_gr_887.gif]   with initial condition   [Graphics:Images/ZTransformIntroModHome_gr_888.gif].

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

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

Solution 14.

See text and/or instructor's solution manual.

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

Solution 14 (a).  

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

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

                    [Graphics:../Images/ZTransformIntroModHome_gr_893.gif].  
        
Solve for  [Graphics:../Images/ZTransformIntroModHome_gr_894.gif]  and get  

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

then use Table 9.1 to find the inverse z-transform   

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

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

 

We are done.   

 

Solution 14 (b).  

        Using residues we get

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

                    and

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

Thus

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

Therefore,

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

 

We are done.   

Aside.  We can let Mathematica double check our work.

 

Take the z-transform of both sides.

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

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


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

[Graphics:../Images/ZTransformIntroModHome_gr_904.gif]
[Graphics:../Images/ZTransformIntroModHome_gr_905.gif]
[Graphics:../Images/ZTransformIntroModHome_gr_906.gif]

Find the inverse z-transform.

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

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


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

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


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

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


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

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


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

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


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

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

          The Maple commands are similar  

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

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


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

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


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

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


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

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

 

Aside.  We can use Mathematica's InverseZTransform subroutine.

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

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


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

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

 

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

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

 

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

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

          The Maple command is similar  

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

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

 

We are really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

 

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

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

          The Maple command is similar  

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

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

 

We are really really done.   

Aside.  We can graph the solution.

 

          [Graphics:../Images/ZTransformIntroModHome_gr_941.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_942.gif]     [Graphics:../Images/ZTransformIntroModHome_gr_943.gif]

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

 

 

We are really really really done.   

The Details for the Partial Fractions.   

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

It is natural to try the command:

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

[Graphics:../Images/ZTransformIntroModHome_gr_947.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_948.gif]  

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

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

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

and solve the linear system  

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

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

Therefore, the desired form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

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

Then use the standard procedure for expanding in partial fractions   

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

Then make the substitution   [Graphics:../Images/ZTransformIntroModHome_gr_963.gif]   in   [Graphics:../Images/ZTransformIntroModHome_gr_964.gif]   and get  

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


Therefore, the desired form is  

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

Aside.   The Mathematica commands for Method (iii)  are

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

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


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

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


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

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


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

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


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

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

Now we have the desired form:

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

 

  

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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