Exercise 8 (b).  Solve   [Graphics:Images/ZTransformDEModHome_gr_1770.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_1771.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_1772.gif].  

Solution 8 (b).

See text and/or instructor's solution manual.

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

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

and then get    

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

Solve for  

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

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

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

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

Therefore,  

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

We are done.   

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

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

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

        Since   [Graphics:../Images/ZTransformDEModHome_gr_1890.gif] and   [Graphics:../Images/ZTransformDEModHome_gr_1891.gif]we can write  

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

                    [Graphics:../Images/ZTransformDEModHome_gr_1893.gif]   for   [Graphics:../Images/ZTransformDEModHome_gr_1894.gif].  

Therefore, the solution has the following form:

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

We are done.   

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

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

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


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

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

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

 

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

[Graphics:../Images/ZTransformDEModHome_gr_1903.gif]
[Graphics:../Images/ZTransformDEModHome_gr_1904.gif]

We are really done.   

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

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

                    and  

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

                    and  

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

Thus,  

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

Therefore,  

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

We are really really done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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


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

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1937.gif]
                                                            
  
[Graphics:../Images/ZTransformDEModHome_gr_1938.gif]  

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

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

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


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

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

We are really really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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


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

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

Aside.  The Maple command is similar  

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

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

We are really really really really done.   

 

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

 

           [Graphics:../Images/ZTransformDEModHome_gr_1950.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1951.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1952.gif]

  

          [Graphics:../Images/ZTransformDEModHome_gr_1953.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1954.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1955.gif]

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

 

We are really really really really really done.   

The Details for the Partial Fractions.   

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

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

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

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

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

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

Then use the standard procedure for expanding in partial fractions   

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

Then make the substitution   [Graphics:../Images/ZTransformDEModHome_gr_1976.gif]   in   [Graphics:../Images/ZTransformDEModHome_gr_1977.gif]   and get  

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (iii)  are

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

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


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

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


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

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


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

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

Use the substitutions

                    [Graphics:../Images/ZTransformDEModHome_gr_1988.gif],  
                    
                    [Graphics:../Images/ZTransformDEModHome_gr_1989.gif],  
                    
                    [Graphics:../Images/ZTransformDEModHome_gr_1990.gif].  

Now we have the desired partial fraction form:

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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