ConceptComplete

Generalized Eigenvectors and Defective Matrices

When a matrix AA has repeated eigenvalues with insufficient eigenvectors, generalized eigenvectors and Jordan chains provide the tools to construct a complete solution.


Defective Matrices

Definition5.4Defective matrix

A matrix AA is defective if it does not have nn linearly independent eigenvectors. This occurs when the geometric multiplicity of some eigenvalue is strictly less than its algebraic multiplicity. For such matrices, Jordan chains replace the missing eigenvectors.

Definition5.5Generalized eigenvector

A generalized eigenvector of rank kk for eigenvalue λ\lambda is a nonzero vector v\mathbf{v} satisfying (AλI)kv=0(A - \lambda I)^k \mathbf{v} = \mathbf{0} but (AλI)k1v0(A - \lambda I)^{k-1}\mathbf{v} \neq \mathbf{0}. A Jordan chain of length kk is v1,v2,,vk\mathbf{v}_1, \mathbf{v}_2, \ldots, \mathbf{v}_k where (AλI)v1=0(A-\lambda I)\mathbf{v}_1 = \mathbf{0} and (AλI)vj+1=vj(A-\lambda I)\mathbf{v}_{j+1} = \mathbf{v}_j.


Solutions from Jordan Chains

Theorem5.5Solutions from generalized eigenvectors

If λ\lambda has algebraic multiplicity mm and geometric multiplicity p<mp < m, the mm linearly independent solutions corresponding to λ\lambda involve terms of the form tkeλtvt^k e^{\lambda t}\mathbf{v} where v\mathbf{v} are generalized eigenvectors. Specifically, for a Jordan chain v1,,vk\mathbf{v}_1, \ldots, \mathbf{v}_k:

xj(t)=eλti=0j1tii!vji,j=1,,k.\mathbf{x}_j(t) = e^{\lambda t}\sum_{i=0}^{j-1}\frac{t^i}{i!}\mathbf{v}_{j-i}, \quad j = 1, \ldots, k.

ExampleDefective $3 \times 3$ system

Solve x=Ax\mathbf{x}' = A\mathbf{x} where A=(210021002)A = \begin{pmatrix}2&1&0\\0&2&1\\0&0&2\end{pmatrix} (one eigenvalue λ=2\lambda=2, multiplicity 3).

Only one eigenvector: v1=(1,0,0)T\mathbf{v}_1 = (1,0,0)^T. Generalized: (A2I)v2=v1(A-2I)\mathbf{v}_2 = \mathbf{v}_1 gives v2=(0,1,0)T\mathbf{v}_2 = (0,1,0)^T. Then (A2I)v3=v2(A-2I)\mathbf{v}_3 = \mathbf{v}_2 gives v3=(0,0,1)T\mathbf{v}_3 = (0,0,1)^T.

x(t)=e2t[c1(100)+c2((010)+t(100))+c3((001)+t(010)+t22(100))].\mathbf{x}(t) = e^{2t}\left[c_1\begin{pmatrix}1\\0\\0\end{pmatrix} + c_2\left(\begin{pmatrix}0\\1\\0\end{pmatrix}+t\begin{pmatrix}1\\0\\0\end{pmatrix}\right) + c_3\left(\begin{pmatrix}0\\0\\1\end{pmatrix}+t\begin{pmatrix}0\\1\\0\end{pmatrix}+\frac{t^2}{2}\begin{pmatrix}1\\0\\0\end{pmatrix}\right)\right].


The Jordan Normal Form

RemarkJordan blocks

Every square matrix is similar to a Jordan normal form J=diag(J1,J2,,Jr)J = \text{diag}(J_1, J_2, \ldots, J_r) where each Jordan block is

Jk=(λk10λk110λk).J_k = \begin{pmatrix}\lambda_k & 1 & & 0 \\ & \lambda_k & 1 & \\ & & \ddots & 1 \\ 0 & & & \lambda_k\end{pmatrix}.

The matrix exponential of a Jordan block is eJkt=eλkt(1tt2/201t01)e^{J_k t} = e^{\lambda_k t}\begin{pmatrix}1 & t & t^2/2 & \cdots \\ 0 & 1 & t & \cdots \\ & & \ddots & \\ 0 & & & 1\end{pmatrix}.

Example$2 \times 2$ Jordan block example

For A=(3115)A = \begin{pmatrix}3 & 1\\-1 & 5\end{pmatrix}, eigenvalue λ=4\lambda = 4 (multiplicity 2, geometric multiplicity 1).

Eigenvector: v1=(1,1)T\mathbf{v}_1 = (1, 1)^T. Generalized: (A4I)v2=v1    v2=(0,1)T(A-4I)\mathbf{v}_2 = \mathbf{v}_1 \implies \mathbf{v}_2 = (0, 1)^T.

x(t)=c1e4t(11)+c2e4t[(01)+t(11)].\mathbf{x}(t) = c_1 e^{4t}\begin{pmatrix}1\\1\end{pmatrix} + c_2 e^{4t}\left[\begin{pmatrix}0\\1\end{pmatrix} + t\begin{pmatrix}1\\1\end{pmatrix}\right].