Exercise 9 (a).  Solve   [Graphics:Images/ZTransformDEModHome_gr_1992.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_1993.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_1994.gif].  

Solution 9 (a).

See text and/or instructor's solution manual.

Answer.   [Graphics:../Images/ZTransformDEModHome_gr_1998.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_1999.gif]:  

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

and then get    

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

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

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

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

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

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

Therefore,  

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

We are done.   

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

and   
[Graphics:../Images/ZTransformDEModHome_gr_2007.gif]   and   [Graphics:../Images/ZTransformDEModHome_gr_2008.gif]   with   [Graphics:../Images/ZTransformDEModHome_gr_2009.gif]   we get:  

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

        Since   [Graphics:../Images/ZTransformDEModHome_gr_2011.gif] and   [Graphics:../Images/ZTransformDEModHome_gr_2012.gif]we can write  

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

                    [Graphics:../Images/ZTransformDEModHome_gr_2014.gif]   for   [Graphics:../Images/ZTransformDEModHome_gr_2015.gif].  

Therefore, the solution has the following form:

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

We are done.   

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

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

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


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

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

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

 

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

[Graphics:../Images/ZTransformDEModHome_gr_2024.gif]
[Graphics:../Images/ZTransformDEModHome_gr_2025.gif]

We are really done.   

Using Residues.   Calculate residue of  [Graphics:../Images/ZTransformDEModHome_gr_2026.gif]  at the pole   [Graphics:../Images/ZTransformDEModHome_gr_2027.gif].  

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

Therefore,  

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

We are really really done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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


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

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

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

                                                            [Graphics:../Images/ZTransformDEModHome_gr_2047.gif]
                                                                            
                                                            
[Graphics:../Images/ZTransformDEModHome_gr_2048.gif]  

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


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

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

We are really really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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

Aside.  The Maple command is similar  

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

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

We are really really really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_2056.gif]     [Graphics:../Images/ZTransformDEModHome_gr_2057.gif]     [Graphics:../Images/ZTransformDEModHome_gr_2058.gif]

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

 

We are really really really really really done.   

The Details for the Partial Fractions.   

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

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

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

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

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

Now we have the desired partial fraction form:

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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