Jordan Normal Form Article Index for
Jordan
Limousines in
Jordan
Website Links For
Jordan
 

Information About

Jordan Normal Form





MOTIVATION

Consider the situation of Matrix Diagonalization . A square matrix is diagonalizable if the sum of the dimensions of the eigenspaces is the number of rows or columns of the matrix. Let us examine the following matrix:

:A=\begin{pmatrix}
322 & -323 & -323 & 322 \
325 & -326 & -325 & 326 \
-259 & 261 & 261 & -260 \
-237 & 237 & 238 & -237 \end{pmatrix}.

We have eigenvalues of ''A'' being only λ = 5, 5, 5, 5. Now, the Dimension of the Kernel of ''A'' − 5''I'' is 1, so ''A'' is not diagonalizable. However, we can construct the Jordan form of this matrix. Since the above dimension is 1, we know that the Jordan form is comprised of only one Jordan Block , that is, the Jordan form of ''A'' is

:J=J_4(5)=\begin{pmatrix}
5 & 1 & 0 & 0 \
0 & 5 & 1 & 0 \
0 & 0 & 5 & 1 \
0 & 0 & 0 & 5 \end{pmatrix}.

Observe that ''J'' can be written as 5''I'' + ''N'', where ''N'' is a Nilpotent Matrix . Since we have now ''A'' similar to such a simple matrix, we can perform calculations involving ''A'' by using the Jordan form, which can ease the calculations in many cases. For example calculating powers of matrices is significantly easier by using the Jordan form.


GENERAL CASE

It is not possible to make all such matrices ''M'' diagonal, even when ''K'' is , ''N'' is Nilpotent , and ''D'' commutes with ''N''.

The way the normal form is usually written is explicitly as the Direct Sum of Block Square Matrices , known as ''Jordan blocks''. Jordan blocks are of the form λ''I'' + ''N'' = ''J''''n''(λ), where λ is one of the eigenvalues of ''M'', ''n'' the number of rows or columns of the Jordan block, and ''N'' is a special nilpotent matrix defined as ''N''''ij'' = δ''i'',''j'' − 1 (where δ is the Kronecker Delta ). This form is valid whenever ''K'' contains the eigenvalues of ''M''. That is, one typical Jordan block looks like

:\begin{pmatrix}
\lambda & 1 & 0 & 0 & \cdots & 0 \
0 & \lambda & 1 & 0 & \cdots & 0 \
0 & 0 & \lambda & 1 & \cdots & 0 \
dots & dots & dots & \ddots & \ddots & dots \
0 & 0 & 0 & 0 & \lambda & 1 \
0 & 0 & 0 & 0 & 0 & \lambda \\end{pmatrix}.

If one knows the dimensions of the kernels of (''M''-λ''I'')''k'' for 1 ≤ ''k'' ≤ ''m'', where ''m'' is the algebraic multiplicity of the eigenvalue λ, one can determine the Jordan form that exists for ''M''. We may view the underlying vector space ''V'' as a ''K'' {Link without Title} - Module by regarding the action of ''x'' on ''V'' as application of ''M'' and extending by ''K''-linearity. Then the polynomials (''x'' − λ)''k'' are called the elementary divisors of ''M'', and the Jordan canonical form is concerned with representing ''M'' in terms of blocks associated to these elementary divisors.

Calculating the invertible transition matrix ''P'' such that ''P''−1''MP'' = ''J'' can be done by considering eigenvectors.

The proof of the Jordan normal form is usually carried out as an application to the Ring ''K'' {Link without Title} of the structure theorem for finitely-generated Modules over Principal Ideal Domain s, of which it is a corollary.


ALGORITHMS AND METHODS

Let us examine the methods of determining the transition matrix by example.


Example 1

Consider the calculation of the transition matrix for the matrix above. Recall

:A=\begin{pmatrix}
322 & -323 & -323 & 322 \
325 & -326 & -325 & 326 \
-259 & 261 & 261 & -260 \
-237 & 237 & 238 & -237 \end{pmatrix}.

We concern ourselves with obtaining Generalized Eigenvector s, that is, solutions to

: (A-\lambda I)^k\mathbf{v} = \mathbf{0}

which will allow us to calculate "chains" of vectors, whose elements form the columns of the transition matrix.

For ''A'' above, we know there is only one Jordan block (see above), so we firstly obtain one generalized eigenvector - since (''A'' − 5''I'')4 is the zero matrix, ker(''A'' − 5''I'')4 is the entire space, so we can pick one of the standard basis vectors for the space, ''v'' = (1,0,0,0)''T'', since none of the standard basis vectors is an eigenvector of (''A'' − 5''I'')3, (''A'' − 5''I'')2, or ''A'' − 5''I''. Then, forming the chain

: \left\{(A-5I)^3\mathbf{v}, (A-5I)^2\mathbf{v}, (A-5I)\mathbf{v}, \mathbf{v} ight\}
: =\left\{
\begin{pmatrix} 5922 \ 4230 \ -3572 \ -5170 \end{pmatrix},
\begin{pmatrix} 2857 \ 2363 \ -1962 \ -2392 \end{pmatrix},
\begin{pmatrix} 317 \ 325 \ -259 \ -237 \end{pmatrix},
\begin{pmatrix} 1 \ 0 \0 \ 0 \end{pmatrix} ight\}

so, we can form the transition matrix as

: P=\begin{pmatrix}
5922 & 2857 & 317 & 1 \
4230 & 2363 & 325 & 0 \
-3572 & -1962 & -259 & 0 \
-5170 & -2392 & -237 & 0 \
\end{pmatrix}.


Example 2

Say we have

:B =
\begin{pmatrix}
5 & 4 & 2 & 1 \
0 & 1 & -1 & -1 \
-1 & -1 & 3 & 0 \
1 & 1 & -1 & 2 \
\end{pmatrix}.

The eigenvalues of ''B'' are 4, 4, 2 and 1. Now, we have

: \mathrm{dim}\ \ker{B-I} = 1, \mathrm{dim}\ \ker{B-2I} = 1,
\mathrm{dim}\ \ker{B-4I} = 1, \mathrm{dim}\ \ker({B-4I})^2 = 2

so we can say that the Jordan form of the matrix is the Direct Sum

:J=J_1(1)\oplus J_1(2)\oplus J_2(4)

since vectors in ker ''B'' − 4''I'' are also in ker (''B'' − 4''I'')2.

We have that

: \ker{(B-4I)}^4 = \mathrm{span}\ \left\{\begin{pmatrix} 0 \ 0 \ -1 \ 1\end{pmatrix}, \begin{pmatrix} 1 \ 0 \ -1 \ 1\end{pmatrix} ight\}

but we pick a vector in the span that is not in any of the kernels of (''B'' − 4''I'')3, (''B'' − 4''I'')2, or ''B'' − 4''I'', so choose ''v'' = (0,0,−1,1)''T'' since (1,0,−1,1) is in the kernel of ''B'' − 4''I''.

Now, there are three chains, {(''B'' − 4''I'')''v'', ''v''}, {'''''w'''''}, and {'''''x'''''}, where '''''w''''' = (1,−1,0,1)''T'' is the basis vector of the 1-dimensional kernel of ''B''-2''I'' and likewise '''''x''''' = ( -1,1,0,0) is the basis vector of the 1-dimensional kernel of ''B'' − ''I''. Form the transition matrix from these chain vectors as follows:

: P=\begin{pmatrix}
-1 & 0 & 1 & -1\
0 & 0 & -1 & 1\
1 & -1 & 0 & 0\