Exercise 14 (a).  Solve   [Graphics:Images/ZTransformDEModHome_gr_2900.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_2901.gif].   
                            Hint.   Get  [Graphics:Images/ZTransformDEModHome_gr_2902.gif].  

Solution 14 (a).

See text and/or instructor's solution manual.

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

Alternative Answer.   [Graphics:../Images/ZTransformDEModHome_gr_2907.gif].  

Remark.   The preferred method to use involves Z-transforms, limits and residues.  

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

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

and then get    

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

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

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

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

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

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

Therefore,  

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

We are done.   

Alternative Solution Using Tables.   Using ordinary partial fractions and formula   [Graphics:../Images/ZTransformDEModHome_gr_2916.gif]   from Table 9.1

and the substitution  [Graphics:../Images/ZTransformDEModHome_gr_2917.gif]  and    [Graphics:../Images/ZTransformDEModHome_gr_2918.gif]  and get

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

Aside.  We can verify the substitution.

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

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

We are done.   

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

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

At the conjugate pole we can use the computation

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

Thus,  

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

Therefore,  

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

We are done.   

Aside.  We can let Mathematica double check our work.

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

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

Remark.  The following two computations were done using Mathematica 7.

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

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

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

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

 

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

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


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

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


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

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

Now use limits.

First, find the limit    [Graphics:../Images/ZTransformDEModHome_gr_2940.gif]

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

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


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

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


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

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


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

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

Second, find the limit    [Graphics:../Images/ZTransformDEModHome_gr_2949.gif]

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

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


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

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


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

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


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

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

Then  

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

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

Aside.  The Maple commands are similar  

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

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

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

                                                            [Graphics:../Images/ZTransformDEModHome_gr_2963.gif]
                                                            
                                                            
[Graphics:../Images/ZTransformDEModHome_gr_2964.gif]  

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


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

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

We are really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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

Aside.  The Maple command is similar  

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

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

For curiosities sake, we can compute some of the terms in the sequence with the various formulae.

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

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


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

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


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

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


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

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

We leave it for the reader to prove that the four forms of the answer are equivalent.

We are really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_2980.gif]     [Graphics:../Images/ZTransformDEModHome_gr_2981.gif]     [Graphics:../Images/ZTransformDEModHome_gr_2982.gif]

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

 

We are really really really done.   

The Details for the Partial Fractions.   

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

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

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

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

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

Now we have the desired partial fraction form:

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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