<<

11/24/2018 - Wikipedia

Generalized eigenvector

In linear , a generalized eigenvector of an n × n is a vector which satisfies certain criteria which are more relaxed than those for an (ordinary) eigenvector.[1]

Let be an n-dimensional ; let be a linear in L(V), the of all linear maps from into itself; and let be the of with respect to some ordered .

There may not always exist a full set of n linearly independent eigenvectors of that form a complete basis for . That is, the matrix may not be diagonalizable.[2][3] This happens when the algebraic multiplicity of at least one eigenvalue is greater than its geometric multiplicity (the nullity of the matrix , or the of its nullspace). In this case, is called a defective eigenvalue and is called a .[4]

A generalized eigenvector corresponding to , together with the matrix generate a Jordan chain of linearly independent generalized eigenvectors which form a basis for an subspace of .[5][6][7]

Using generalized eigenvectors, a set of linearly independent eigenvectors of can be extended, if necessary, to a complete basis for .[8] This basis can be used to determine an "almost matrix" in , similar to , which is useful in computing certain matrix functions of .[9] The matrix is also useful in solving the system of linear differential where need not be diagonalizable.[10][11]

Contents

Overview and definition Examples Example 1 Example 2 Jordan chains Computation of generalized eigenvectors Example 3 Generalized Jordan normal form Example 4 Example 5 Applications Matrix functions Differential equations Notes References https://en.wikipedia.org/wiki/Generalized_eigenvector 1/15 11/24/2018 Generalized eigenvector - Wikipedia Overview and definition

There are several equivalent ways to define an ordinary eigenvector.[12][13][14][15][16][17][18][19] For our purposes, an eigenvector associated with an eigenvalue of an × matrix is a nonzero vector for which , where is the × and is the zero vector of length .[20] That is, is in the of the transformation . If has linearly independent eigenvectors, then is similar to a . That is, there exists an such that is diagonalizable through the similarity transformation .[21][22] The matrix is called a spectral matrix for . The matrix is called a modal matrix for .[23] Diagonalizable matrices are of particular interest since matrix functions of them can be computed easily.[24]

On the other hand, if does not have linearly independent eigenvectors associated with it, then is not diagonalizable.[25][26]

Definition: A vector is a generalized eigenvector of m of the matrix and corresponding to the eigenvalue if

but

[27]

Clearly, a generalized eigenvector of rank 1 is an ordinary eigenvector.[28] Every × matrix has linearly independent generalized eigenvectors associated with it and can be shown to be similar to an "almost diagonal" matrix in Jordan normal form.[29] That is, there exists an invertible matrix such that .[30] The matrix in this case is called a generalized modal matrix for .[31] If is an eigenvalue of algebraic multiplicity , then will have linearly independent generalized eigenvectors corresponding to .[32] These results, in turn, provide a straightforward method for computing certain matrix functions of .[33]

Note: For an matrix over a to be expressed in Jordan normal form, all eigenvalues of must be in . That is, the must factor completely into linear factors. For example, if has real-valued elements, then it may be necessary for the eigenvalues and the components of the eigenvectors to have complex values.[34][35][36]

The set spanned by all generalized eigenvectors for a given , forms the generalized eigenspace for .[37]

Examples

Here are some examples to illustrate the concept of generalized eigenvectors. Some of the details will be described later.

Example 1 This example is simple but clearly illustrates the point. This type of matrix is used frequently in textbooks.[38][39][40] Suppose

Then there is only one eigenvalue, , and its algebraic multiplicity is m = 2.

https://en.wikipedia.org/wiki/Generalized_eigenvector 2/15 11/24/2018 Generalized eigenvector - Wikipedia Notice that this matrix is in Jordan normal form but is not diagonal. Hence, this matrix is not diagonalizable. Since there is one superdiagonal entry, there will be one generalized eigenvector of rank greater than 1 (or one could note that the vector space is of dimension 2, so there can be at most one generalized eigenvector of rank greater than 1). Alternatively, one could compute the dimension of the nullspace of to be p = 1, and thus there are m – p = 1 generalized eigenvectors of rank greater than 1.

The ordinary eigenvector is computed as usual (see the eigenvector page for examples). Using this eigenvector, we compute the generalized eigenvector by solving

Writing out the values:

This simplifies to

The has no restrictions. The generalized eigenvector of rank 2 is then , where a can have any value. The choice of a = 0 is usually the simplest.

Note that

so that is a generalized eigenvector,

so that is an ordinary eigenvector, and that and are linearly independent and hence constitute a basis for the vector space .

Example 2 This example is more complex than Example 1. Unfortunately, it is a little difficult to construct an interesting example of low order.[41] The matrix

https://en.wikipedia.org/wiki/Generalized_eigenvector 3/15 11/24/2018 Generalized eigenvector - Wikipedia has eigenvalues and with algebraic multiplicities and , but geometric multiplicities and .

The generalized eigenspaces of are calculated below. is the ordinary eigenvector associated with . is a generalized eigenvector associated with . is the ordinary eigenvector associated with . and are generalized eigenvectors associated with .

This results in a basis for each of the generalized eigenspaces of . Together the two chains of generalized eigenvectors span the space of all 5-dimensional column vectors.

https://en.wikipedia.org/wiki/Generalized_eigenvector 4/15 11/24/2018 Generalized eigenvector - Wikipedia

An "almost diagonal" matrix in Jordan normal form, similar to is obtained as follows:

where is a generalized modal matrix for , the columns of are a canonical basis for , and .[42]

Jordan chains

Definition: Let be a generalized eigenvector of rank m corresponding to the matrix and the eigenvalue . The chain generated by is a set of vectors given by

(1)

Thus, in general,

(2) https://en.wikipedia.org/wiki/Generalized_eigenvector 5/15 11/24/2018 Generalized eigenvector - Wikipedia

The vector , given by (2), is a generalized eigenvector of rank j corresponding to the eigenvalue . A chain is a linearly independent set of vectors.[43]

Canonical basis

Definition: A set of n linearly independent generalized eigenvectors is a canonical basis if it is composed entirely of Jordan chains.

Thus, once we have determined that a generalized eigenvector of rank m is in a canonical basis, it follows that the m − 1 vectors that are in the Jordan chain generated by are also in the canonical basis.[44]

Let be an eigenvalue of of algebraic multiplicity . First, find the ranks (matrix ranks) of the matrices . The integer is determined to be the first integer for which has rank (n being the number of rows or columns of , that is, is n × n).

Now define

The variable designates the number of linearly independent generalized eigenvectors of rank k corresponding to the eigenvalue that will appear in a canonical basis for . Note that

.[45]

Computation of generalized eigenvectors

In the preceding sections we have seen techniques for obtaining the n linearly independent generalized eigenvectors of a canonical basis for the vector space associated with an n × n matrix . These techniques can be combined into a procedure:

Solve the characteristic of for eigenvalues and their algebraic multiplicities ; For each

Determine ; Determine ; Determine for ; Determine each Jordan chain for ;

Example 3 The matrix

https://en.wikipedia.org/wiki/Generalized_eigenvector 6/15 11/24/2018 Generalized eigenvector - Wikipedia

has an eigenvalue of algebraic multiplicity and an eigenvalue of algebraic multiplicity . We also have n = 4. For we have .

The first integer for which has rank is .

We now define

Consequently, there will be three linearly independent generalized eigenvectors; one each of ranks 3, 2 and 1. Since corresponds to a single chain of three linearly independent generalized eigenvectors, we know that there is a generalized eigenvector of rank 3 corresponding to such that

(3)

but

(4)

https://en.wikipedia.org/wiki/Generalized_eigenvector 7/15 11/24/2018 Generalized eigenvector - Wikipedia Equations (3) and (4) represent linear systems that can be solved for . Let

Then

and

Thus, in order to satisfy the conditions (3) and (4), we must have and . No restrictions are placed on and . By choosing , we obtain

as a generalized eigenvector of rank 3 corresponding to . Note that it is possible to obtain infinitely many other generalized eigenvectors of rank 3 by choosing different values of , and , with . Our first choice, however, is the simplest.[46]

Now using equations (1), we obtain and as generalized eigenvectors of rank 2 and 1 respectively, where

and https://en.wikipedia.org/wiki/Generalized_eigenvector 8/15 11/24/2018 Generalized eigenvector - Wikipedia

The simple eigenvalue can be dealt with using standard techniques and has an ordinary eigenvector

A canonical basis for is

and are generalized eigenvectors associated with . is the ordinary eigenvector associated with .

It should be noted that this is a fairly simple example. In general, the numbers of linearly independent generalized eigenvectors of rank k will not always be equal. That is, there may be several chains of different lengths corresponding to a particular eigenvalue.[47]

Generalized modal matrix

Let be an n × n matrix. A generalized modal matrix for is an n × n matrix whose columns, considered as vectors, form a canonical basis for and appear in according to the following rules:

All Jordan chains consisting of one vector (that is, one vector in length) appear in the first columns of . All vectors of one chain appear together in adjacent columns of . Each chain appears in in order of increasing rank (that is, the generalized eigenvector of rank 1 appears before the generalized eigenvector of rank 2 of the same chain, which appears before the generalized eigenvector of rank 3 of the same chain, etc.).[48]

Jordan normal form

Let be an n-dimensional vector space; let be a in L(V), the set of all linear maps from into itself; and let be the matrix representation of with respect to some ordered basis. It can be shown that if the characteristic polynomial of factors into linear factors, so that has the form

https://en.wikipedia.org/wiki/Generalized_eigenvector 9/15 11/24/2018 Generalized eigenvector - Wikipedia where are the distinct eigenvalues of , then each is the algebraic multiplicity of its corresponding eigenvalue and is similar to a matrix in Jordan normal form, where each appears consecutive times on the diagonal, and the entry directly above each (that is, on the superdiagonal) is either 0 or 1: the entry above the first occurrence of each is always 0; all other entries on the superdiagonal are 1. All other entries (that is, off the diagonal and superdiagonal) are 0. The matrix is as close as one can come to a diagonalization of . If is diagonalizable, then all entries above the diagonal are zero.[49] Note that some textbooks have the ones on the subdiagonal, that is, immediately below the main diagonal instead of on the superdiagonal. The eigenvalues are still on the main diagonal.[50][51]

Every n × n matrix is similar to a matrix in Jordan normal form, obtained through the similarity transformation , where is a generalized modal matrix for .[52] (See Note above.)

Example 4 An example of a matrix in Jordan normal Find a matrix in Jordan normal form that is similar to form. The grey blocks are called Jordan blocks.

Solution: The characteristic equation of is , hence, is an eigenvalue of algebraic multiplicity three. Following the procedures of the previous sections, we find that

and

Thus, and , which implies that a canonical basis for will contain one linearly independent generalized eigenvector of rank 2 and two linearly independent generalized eigenvectors of rank 1, or equivalently, one chain of two vectors and one chain of one vector . Designating , we find that

and

https://en.wikipedia.org/wiki/Generalized_eigenvector 10/15 11/24/2018 Generalized eigenvector - Wikipedia where is a generalized modal matrix for , the columns of are a canonical basis for , and .[53] Note that since generalized eigenvectors themselves are not unique, and since some of the columns of both and may be interchanged, it follows that both and are not unique.[54]

Example 5 In Example 3, we found a canonical basis of linearly independent generalized eigenvectors for a matrix . A generalized modal matrix for is

A matrix in Jordan normal form, similar to is

so that .

Applications

Matrix functions Three of the most fundamental operations which can be performed on square matrices are matrix , multiplication by a scalar, and .[55] These are exactly those operations necessary for defining a polynomial of an n × n matrix .[56] If we recall from basic that many functions can be written as a Maclaurin series, then we can define more general functions of matrices quite easily.[57] If is diagonalizable, that is

with

then https://en.wikipedia.org/wiki/Generalized_eigenvector 11/15 11/24/2018 Generalized eigenvector - Wikipedia

and the evaluation of the Maclaurin series for functions of is greatly simplified.[58] For example, to obtain any power k of , we need only compute , premultiply by , and postmultiply the result by .[59]

Using generalized eigenvectors, we can obtain the Jordan normal form for and these results can be generalized to a straightforward method for computing functions of nondiagonalizable matrices.[60] (See Matrix function#Jordan decomposition.)

Differential equations Consider the problem of solving the system of linear ordinary differential equations

(5)

where

and

If the matrix is a diagonal matrix so that for , then the system (5) reduces to a system of n equations which take the form

(6)

In this case, the general solution is given by

https://en.wikipedia.org/wiki/Generalized_eigenvector 12/15 11/24/2018 Generalized eigenvector - Wikipedia

In the general case, we try to diagonalize and reduce the system (5) to a system like (6) as follows. If is diagonalizable, we have , where is a modal matrix for . Substituting , equation (5) takes the form , or

(7)

where

(8)

The solution of (7) is

The solution of (5) is then obtained using the (8).[61]

On the other hand, if is not diagonalizable, we choose to be a generalized modal matrix for , such that is the Jordan normal form of . The system has the form

(9)

where the are the eigenvalues from the main diagonal of and the are the ones and zeros from the superdiagonal of . The system (9) is often more easily solved than (5). We may solve the last equation in (9) for , obtaining . We then substitute this solution for into the next to last equation in (9) and solve for . Continuing this procedure, we work through (9) from the last equation to the first, solving the entire system for . The solution is then obtained using the relation (8).[62] https://en.wikipedia.org/wiki/Generalized_eigenvector 13/15 11/24/2018 Generalized eigenvector - Wikipedia Notes

1. Bronson (1970, p. 189) 32. Bronson (1970, p. 196) 2. Beauregard & Fraleigh (1973, p. 310) 33. Bronson (1970, pp. 189,209–215) 3. Nering (1970, p. 118) 34. Golub & Van Loan (1996, p. 316) 4. Golub & Van Loan (1996, p. 316) 35. Herstein (1964, p. 259) 5. Beauregard & Fraleigh (1973, p. 319) 36. Nering (1970, p. 118) 6. Bronson (1970, pp. 194–195) 37. Nering (1970, p. 118) 7. Golub & Van Loan (1996, p. 311) 38. Nering (1970, p. 118) 8. Bronson (1970, p. 196) 39. Herstein (1964, p. 261) 9. Bronson (1970, p. 189) 40. Beauregard & Fraleigh (1973, p. 310) 10. Beauregard & Fraleigh (1973, pp. 316–318) 41. Nering (1970, pp. 122,123) 11. Nering (1970, p. 118) 42. Bronson (1970, pp. 189–209) 12. Anton (1987, pp. 301–302) 43. Bronson (1970, pp. 194–195) 13. Beauregard & Fraleigh (1973, p. 266) 44. Bronson (1970, pp. 196,197) 14. Burden & Faires (1993, p. 401) 45. Bronson (1970, pp. 197,198) 15. Golub & Van Loan (1996, pp. 310–311) 46. Bronson (1970, pp. 190–191) 16. Harper (1976, p. 58) 47. Bronson (1970, pp. 197–198) 17. Herstein (1964, p. 225) 48. Bronson (1970, p. 205) 18. Kreyszig (1972, pp. 273,684) 49. Beauregard & Fraleigh (1973, p. 311) 19. Nering (1970, p. 104) 50. Cullen (1966, p. 114) 20. Burden & Faires (1993, p. 401) 51. Franklin (1968, p. 122) 21. Beauregard & Fraleigh (1973, pp. 270–274) 52. Bronson (1970, p. 207) 22. Bronson (1970, pp. 179–183) 53. Bronson (1970, pp. 208) 23. Bronson (1970, p. 181) 54. Bronson (1970, p. 206) 24. Bronson (1970, p. 179) 55. Beauregard & Fraleigh (1973, pp. 57–61) 25. Beauregard & Fraleigh (1973, pp. 270–274) 56. Bronson (1970, p. 104) 26. Bronson (1970, pp. 179–183) 57. Bronson (1970, p. 105) 27. Bronson (1970, p. 189) 58. Bronson (1970, p. 184) 28. Bronson (1970, pp. 190,202) 59. Bronson (1970, p. 185) 29. Bronson (1970, pp. 189,203) 60. Bronson (1970, pp. 209–218) 30. Bronson (1970, pp. 206–207) 61. Beauregard & Fraleigh (1973, pp. 274–275) 31. Bronson (1970, p. 205) 62. Beauregard & Fraleigh (1973, p. 317)

References

Anton, Howard (1987), Elementary (5th ed.), New York: Wiley, ISBN 0-471-84819-0 Axler, Sheldon (1997). Linear Algebra Done Right (2nd ed.). Springer. ISBN 978-0-387-98258-8. https://en.wikipedia.org/wiki/Generalized_eigenvector 14/15 11/24/2018 Generalized eigenvector - Wikipedia Beauregard, Raymond A.; Fraleigh, John B. (1973), A First Course In Linear Algebra: with Optional Introduction to Groups, Rings, and Fields, Boston: Houghton Mifflin Co., ISBN 0- 395-14017-X Bronson, Richard (1970), Matrix Methods: An Introduction, New York: Academic Press, LCCN 70097490 (https://lccn.loc.gov/70097490) Burden, Richard L.; Faires, J. Douglas (1993), (5th ed.), Boston: Prindle, Weber and Schmidt, ISBN 0-534-93219-3 Cullen, Charles G. (1966), Matrices and Linear Transformations, Reading: Addison-Wesley, LCCN 66021267 (https://lccn.loc.gov/66021267) Franklin, Joel N. (1968), Matrix Theory, Englewood Cliffs: Prentice-Hall, LCCN 68016345 (https://lccn.loc.gov/68016345) Golub, Gene H.; Van Loan, Charles F. (1996), Matrix Computations (3rd ed.), Baltimore: Johns Hopkins University Press, ISBN 0-8018-5414-8 Harper, Charlie (1976), Introduction to , New Jersey: Prentice-Hall, ISBN 0-13-487538-9 Herstein, I. N. (1964), Topics In Algebra, Waltham: Blaisdell Publishing Company, ISBN 978-1114541016 Kreyszig, Erwin (1972), Advanced Mathematics (3rd ed.), New York: Wiley, ISBN 0-471-50728-8 Nering, Evar D. (1970), Linear Algebra and Matrix Theory (2nd ed.), New York: Wiley, LCCN 76091646 (https://lccn.loc.gov/76091646)

Retrieved from "https://en.wikipedia.org/w/index.php?title=Generalized_eigenvector&oldid=868856580"

This page was last edited on 14 November 2018, at 22:03 (UTC).

Text is available under the Creative Commons Attribution-ShareAlike License; additional terms may apply. By using this site, you agree to the Terms of Use and Privacy Policy. Wikipedia® is a registered trademark of the Wikimedia Foundation, Inc., a non-profit organization.

https://en.wikipedia.org/wiki/Generalized_eigenvector 15/15