Exercise 6.  Solve the nonhomogeneous difference equations.  

Exercise 6 (b).  [Graphics:Images/ZTransformDEModHome_gr_1207.gif]   with   [Graphics:Images/ZTransformDEModHome_gr_1208.gif].   
                          Hint.  Get  [Graphics:Images/ZTransformDEModHome_gr_1209.gif].  

Solution 6 (b).

See text and/or instructor's solution manual.

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

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

and then get    

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

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

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

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

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

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

Therefore,  

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

We are done.   

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

and   
[Graphics:../Images/ZTransformDEModHome_gr_1327.gif]   and   [Graphics:../Images/ZTransformDEModHome_gr_1328.gif]   we get:  

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

        Since   [Graphics:../Images/ZTransformDEModHome_gr_1330.gif] and   [Graphics:../Images/ZTransformDEModHome_gr_1331.gif]we can write  

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

                    [Graphics:../Images/ZTransformDEModHome_gr_1333.gif][Graphics:../Images/ZTransformDEModHome_gr_1334.gif]   for   [Graphics:../Images/ZTransformDEModHome_gr_1335.gif].  

Therefore, the solution has the following form:

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

We are done.   

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

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

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


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

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

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

 

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

[Graphics:../Images/ZTransformDEModHome_gr_1344.gif]
[Graphics:../Images/ZTransformDEModHome_gr_1345.gif]

We are really done.   

Using Residues.   Calculate the residues  [Graphics:../Images/ZTransformDEModHome_gr_1346.gif]  at the poles   [Graphics:../Images/ZTransformDEModHome_gr_1347.gif].  

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

                    and  

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

Thus,  

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

Therefore,  

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

We are really really done.   

Aside.  We can let Mathematica double check our work.

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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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


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

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

Aside.  The Maple commands are similar  

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

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


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

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


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

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


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

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


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

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

We are really really really done.   

Aside.  We can use Mathematica's Rsolve subroutine.

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

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

Aside.  The Maple command is similar  

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

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

We are really really really really done.   

 

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

 

          [Graphics:../Images/ZTransformDEModHome_gr_1388.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1389.gif]     [Graphics:../Images/ZTransformDEModHome_gr_1390.gif]

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

 

We are really really really really really done.   

The Details for the Partial Fractions.   

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

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

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

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

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

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

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

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

and solve the linear system  

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

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

Therefore, the desired partial fraction form is  

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

Aside.   The Mathematica commands for Method (ii)  are

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

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


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

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

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

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

This solution is complements of the authors.

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

 

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