Exercise 5 (b).  Lucas numbers.  Solve  [Graphics:Images/ZTransformDEModHome_gr_1056.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_1057.gif].  
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_1058.gif].  

Solution 5 (b).

See text and/or instructor's solution manual.

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

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

Solution.   

Method 1.
  The characteristic equation  

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

has roots [Graphics:../Images/ZTransformDEModHome_gr_1133.gif] and [Graphics:../Images/ZTransformDEModHome_gr_1134.gif].  

The general solution is  

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

Solve the linear system  

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

and get  [Graphics:../Images/ZTransformDEModHome_gr_1137.gif]  and  [Graphics:../Images/ZTransformDEModHome_gr_1138.gif].

Therefore,  

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

Aside.  A few terms in the sequence are

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

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

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

and then get    

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

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

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

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

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

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

Therefore,  

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

We are done.   

Using Residues.   Calculate the residues  [Graphics:../Images/ZTransformDEModHome_gr_1149.gif]  at the poles   [Graphics:../Images/ZTransformDEModHome_gr_1150.gif].  

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

and

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

Thus,  

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

Therefore,  

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

We are done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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

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

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1176.gif]
                                                            
                                                            
[Graphics:../Images/ZTransformDEModHome_gr_1177.gif]  

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


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

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

We are really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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

Aside.  The Maple command is similar  

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

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

Remark. The Lucas numbers start with the subscript  [Graphics:../Images/ZTransformDEModHome_gr_1185.gif]  and are

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

We are really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_1187.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1188.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1189.gif]

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

 

We are really really really done.   

The Details for the Partial Fractions.   

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

Method (i).   Use the following algebra steps  

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

Therefore, the desired partial fraction form is  

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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