Theorem 5.1.λ is an eigenvalue of A if and only if det(A−λI)=0.
The polynomial p(λ)=det(A−λI) is called the characteristic polynomial of A. Its roots (in the algebraic closure of F) are the eigenvalues of A.
If p(λ)=(lambda−lambda1)m1(lambda−lambda2)m2⋯(lambda−lambdak)mk With lambda1,…,lambdak distinct, then mi is the algebraic multiplicity of lambdai.
Proposition 5.2. For each eigenvalue λ, 1≤dim(Eλ)≤mλ (geometric multiplicity does not exceed algebraic multiplicity).
Definition.A is diagonalisable if there exists an invertible matrix P and a diagonal Matrix D such that
A=PDP−1
Theorem 5.3.A∈Mn×n(F) is diagonalisable (over F) if and only if A has n linearly independent Eigenvectors (over F). Equivalently, the sum of the geometric multiplicities equals n.
Corollary 5.4. If A has n distinct eigenvalues, then A is diagonalisable.
Proof. Eigenvectors corresponding to distinct eigenvalues are linearly independent. With n distinct Eigenvalues, we obtain n linearly independent eigenvectors, which form a basis of Fn. ■
Theorem 5.5 (Cayley—Hamilton). Every square matrix satisfies its own characteristic polynomial: If p(λ)=det(λI−A) Then p(A)=0 (the zero matrix).
Proof sketch. Let p(λ)=λn+cn−1λn−1+⋯+c1λ+c0. By the adjugate formula (Theorem 3.5), (λI−A)⋅adj(λI−A)=p(λ)⋅I. Each entry of adj(λI−A) is a polynomial in λ of degree at most n−1 So we can write adj(λI−A)=Bn−1λn−1+⋯+B1λ+B0 for Matrices Bi. Multiplying out and comparing coefficients of λk:
Bn−1=I,Bn−2−ABn−1=cn−1I,…,−AB0=c0I
Multiplying the k-th equation on the left by Ak and summing over k:
When a matrix is not diagonalisable, the Jordan normal form provides the next-best canonical Representation.
Theorem 5.6. Let A∈Mn×n(C). Then A is similar to a block-diagonal Matrix
J=J1⋱Jk
Where each Jordan block has the form
Ji=λi1λi⋱⋱1λi
The Jordan form is unique up to permutation of the blocks.
Intuition. Each Jordan block corresponds to one eigenvalue. The size of the block equals the Number of steps in the chain v,(A−λI)v,(A−λI)2v,… Of generalised eigenvectors. A diagonalisable matrix has all Jordan blocks of size 1×1.
Problem. Find the Jordan normal form of
A=(3013)
Solution
The characteristic polynomial is det(A−λI)=(3−λ)2 So λ=3 is the Only eigenvalue with algebraic multiplicity 2.
A−3I=(0010)Which has rank 1, so the geometric Multiplicity is dim(ker(A−3I))=2−1=1.
Since the geometric multiplicity (1) is less than the algebraic multiplicity (2), A is not Diagonalisable. The Jordan form has one block of size 2:
Theorem 5.7 (Spectral Theorem). If A∈Mn×n(R) is symmetric (A=AT), then:
All eigenvalues of A are real.
A has n linearly independent orthonormal eigenvectors.
A is orthogonally diagonalisable: A=QDQT where Q is orthogonal (QTQ=I).
Proof. We prove (1) and then (2) and (3) by induction on n.
(1) Let λ∈C be an eigenvalue with eigenvector v∈Cnv=0. Then
vTAv=vT(λv)=λvTv
Since A=AT and A has real entries, A=A=AT So
vTAv=(Av)Tv=(Av)Tv=(λv)Tv=λvTv
Therefore (λ−λ)vTv=0. Since vTv>0 We have λ=λ So λ∈R.
(2) and (3) By induction. For n=1 the result is trivial. Assume it holds for (n−1)×(n−1) Symmetric matrices. Since all eigenvalues are real, A has a real eigenvalue λ1 with real Eigenvector v1. Normalise: q1=v1/∥v1∥.
Let W=q1⊥={w∈Rn:q1Tw=0}. For any w∈W:
q1T(Aw)=(Aq1)Tw=(λ1q1)Tw=λ1⋅0=0
So Aw∈W. Therefore A restricts to a symmetric linear map A∣W:W→W on an (n−1)-dimensional space. By the inductive hypothesis, W has an orthonormal basis {q2,…,qn} of eigenvectors of A∣W.
Then {q1,q2,…,qn} is an orthonormal eigenbasis for Rn And A=QDQT with Q=[q1∣⋯∣qn]. ■
An eigenvector of a linear transformation A is a direction that is preserved by A: applying A to an eigenvector only stretches or compresses it, without rotating it. The eigenvalue λ measures how much the eigenvector is stretched: if λ>1, the eigenvector is elongated; if 0<λ<1, it is compressed; if λ<0, it is reflected and then scaled.
Geometrically, the eigenvectors of a matrix reveal the “natural axes” of the transformation. For a 2×2 matrix, the eigenvectors define the directions along which the transformation acts as simple scaling. In the eigenbasis, the matrix is diagonal: the transformation is just independent stretching along each axis.
This is why diagonalisation is powerful. Computing A100 directly requires 100 matrix multiplications, but D100 is trivial: just raise each diagonal entry to the 100th power. The change-of-basis matrix P handles the translation between the standard basis and the eigenbasis.
Physical examples:
In mechanics, the eigenvectors of the inertia tensor are the principal axes of rotation.
In vibration analysis, the eigenvectors of the stiffness matrix are the normal modes, and the eigenvalues are the squared natural frequencies.
In Google’s PageRank, the eigenvector corresponding to eigenvalue 1 of the web matrix gives the importance ranking of all pages.
Why eigenvalues matter for stability: In the ODE x′=Ax, the solution is x(t)=eAtx0. If A is diagonalisable, this becomes x(t)=PeDtP−1x0. The behavior is governed by eλit: if all Re(λi)<0, the system decays to zero (stable); if any Re(λi)>0, the system grows without bound (unstable).
Since A is upper triangular, the eigenvalues are the diagonal entries: λ1=2, λ2=3 (with algebraic multiplicity 2).
For λ1=2: Solve (A−2I)v=0.
A−2I=000110001→000100010
Free variable: x1=t. Eigenvector: v1=(1,0,0)T.
For λ2=3: Solve (A−3I)v=0.
A−3I=−100100000
Free variables: x2=s, x3=t. Then x1=s. Eigenvectors: s(1,1,0)T+t(0,0,1)T.
The geometric multiplicity of λ2=3 is 2, equal to its algebraic multiplicity. Therefore A is diagonalisable with
P=100110001,D=200030003
Key observation: When the geometric multiplicity equals the algebraic multiplicity for every eigenvalue, the matrix is diagonalisable. When they differ (as in the Jordan form example in Section 5.5), the matrix is not diagonalisable, and the Jordan normal form is the best alternative. ■
Problem. Use the Cayley—Hamilton theorem to compute A10 for the same matrix A above.
Solution
The characteristic polynomial is p(λ)=λ2−7λ+10 So by Cayley—Hamilton, A2=7A−10I.
To find A10Divide λ10 by p(λ):
λ10=q(λ)(λ2−7λ+10)+r(λ)
Where r(λ)=aλ+b has degree less than 2. Then A10=r(A)=aA+bI.
To find a and bEvaluate at the eigenvalues:
λ10λ=5=510=9765625=5a+b
λ10λ=2=210=1024=2a+b
Subtracting: 3a=9765625−1024=9764601 So a=3254867.