The matrix
has eigenvalues and corresponding eigenvectors
A diagonal matrix D {\displaystyle D} , similar to A {\displaystyle A} is
One possible choice for an invertible matrix M {\displaystyle M} such that D = M − 1 A M , {\displaystyle D=M^{-1}AM,} is
Note that since eigenvectors themselves are not unique, and since the columns of both M {\displaystyle M} and D {\displaystyle D} may be interchanged, it follows that both M {\displaystyle M} and D {\displaystyle D} are not unique.4
Let A {\displaystyle A} be an n × n matrix. A generalized modal matrix M {\displaystyle M} for A {\displaystyle A} is an n × n matrix whose columns, considered as vectors, form a canonical basis for A {\displaystyle A} and appear in M {\displaystyle M} according to the following rules:
One can show that
where J {\displaystyle J} is a matrix in Jordan normal form. By premultiplying by M − 1 {\displaystyle M^{-1}} , we obtain
Note that when computing these matrices, equation (1) is the easiest of the two equations to verify, since it does not require inverting a matrix.6
This example illustrates a generalized modal matrix with four Jordan chains. Unfortunately, it is a little difficult to construct an interesting example of low order.7 The matrix
has a single eigenvalue λ 1 = 1 {\displaystyle \lambda _{1}=1} with algebraic multiplicity μ 1 = 7 {\displaystyle \mu _{1}=7} . A canonical basis for A {\displaystyle A} will consist of one linearly independent generalized eigenvector of rank 3 (generalized eigenvector rank; see generalized eigenvector), two of rank 2 and four of rank 1; or equivalently, one chain of three vectors { x 3 , x 2 , x 1 } {\displaystyle \left\{\mathbf {x} _{3},\mathbf {x} _{2},\mathbf {x} _{1}\right\}} , one chain of two vectors { y 2 , y 1 } {\displaystyle \left\{\mathbf {y} _{2},\mathbf {y} _{1}\right\}} , and two chains of one vector { z 1 } {\displaystyle \left\{\mathbf {z} _{1}\right\}} , { w 1 } {\displaystyle \left\{\mathbf {w} _{1}\right\}} .
An "almost diagonal" matrix J {\displaystyle J} in Jordan normal form, similar to A {\displaystyle A} is obtained as follows:
where M {\displaystyle M} is a generalized modal matrix for A {\displaystyle A} , the columns of M {\displaystyle M} are a canonical basis for A {\displaystyle A} , and A M = M J {\displaystyle AM=MJ} .8 Note that since generalized eigenvectors themselves are not unique, and since some of the columns of both M {\displaystyle M} and J {\displaystyle J} may be interchanged, it follows that both M {\displaystyle M} and J {\displaystyle J} are not unique.9
Bronson (1970, pp. 179–183) - Bronson, Richard (1970), Matrix Methods: An Introduction, New York: Academic Press, LCCN 70097490 https://lccn.loc.gov/70097490 ↩
Bronson (1970, p. 181) - Bronson, Richard (1970), Matrix Methods: An Introduction, New York: Academic Press, LCCN 70097490 https://lccn.loc.gov/70097490 ↩
Beauregard & Fraleigh (1973, pp. 271, 272) - 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 https://archive.org/details/firstcourseinlin0000beau ↩
Bronson (1970, p. 205) - Bronson, Richard (1970), Matrix Methods: An Introduction, New York: Academic Press, LCCN 70097490 https://lccn.loc.gov/70097490 ↩
Bronson (1970, pp. 206–207) - Bronson, Richard (1970), Matrix Methods: An Introduction, New York: Academic Press, LCCN 70097490 https://lccn.loc.gov/70097490 ↩
Nering (1970, pp. 122, 123) - Nering, Evar D. (1970), Linear Algebra and Matrix Theory (2nd ed.), New York: Wiley, LCCN 76091646 https://lccn.loc.gov/76091646 ↩
Bronson (1970, pp. 208, 209) - Bronson, Richard (1970), Matrix Methods: An Introduction, New York: Academic Press, LCCN 70097490 https://lccn.loc.gov/70097490 ↩
Bronson (1970, p. 206) - Bronson, Richard (1970), Matrix Methods: An Introduction, New York: Academic Press, LCCN 70097490 https://lccn.loc.gov/70097490 ↩