Exercise 6.  Solve the nonhomogeneous difference equations.  

Exercise 6 (a).  [Graphics:Images/ZTransformDEModHome_gr_1204.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_1205.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_1206.gif].  

Solution 6 (a).

See text and/or instructor's solution manual.

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

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

and get  

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

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

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

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

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

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

Therefore,  

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

We are done.   

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

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

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

        Since   [Graphics:../Images/ZTransformDEModHome_gr_1224.gif] and   [Graphics:../Images/ZTransformDEModHome_gr_1225.gif]we can write  

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

                    [Graphics:../Images/ZTransformDEModHome_gr_1227.gif]   for   [Graphics:../Images/ZTransformDEModHome_gr_1228.gif].  

Therefore, the solution has the following form:

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

We are done.   

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

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

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


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

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

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

 

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

[Graphics:../Images/ZTransformDEModHome_gr_1237.gif]
[Graphics:../Images/ZTransformDEModHome_gr_1238.gif]

We are really done.   

Using Residues.   Calculate the residues  [Graphics:../Images/ZTransformDEModHome_gr_1239.gif]  at the poles   [Graphics:../Images/ZTransformDEModHome_gr_1240.gif].  

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

                    and

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

                    and

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

Thus,  

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

Therefore,  

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

We are really really done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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


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

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


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

                                                            [Graphics:../Images/ZTransformDEModHome_gr_1273.gif]
                                                
                                                            
[Graphics:../Images/ZTransformDEModHome_gr_1274.gif]  

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


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

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

We are really really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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

Aside.  The Maple command is similar  

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

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

We are really really really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_1282.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1283.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1284.gif]

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

 

We are really really really done.   

The Details for the Partial Fractions.   

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

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

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

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

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

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

Then use the standard procedure for expanding in partial fractions   

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

Then make the substitution   [Graphics:../Images/ZTransformDEModHome_gr_1305.gif]   in   [Graphics:../Images/ZTransformDEModHome_gr_1306.gif]   and get  

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


Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (iii)  are

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

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


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

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


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

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


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

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

Now we have the desired partial fraction form:

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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