Also Eigenvectors · Eigenvektor · Eigenspace · Eigenraum · Eigenbasis
An eigenvector of a square matrix A is a vector x ≠ 0 whose direction A does not change, Ax = λx; the eigenvectors for λ together with 0 form the eigenspace ker(A − λI).
Only the direction is determined: if x is an eigenvector for λ, so is cx for every c=0, since A(cx)=cAx=λ(cx) (MML p. 105). An eigenvector is therefore usually given up to scaling, e.g. normalized to ∥x∥=1. Two vectors that point in the same or the opposite direction are called collinear.
Eigenspaces are subspaces:Eλ is the solution set of the homogeneous linear system(A−λI)x=0, so it is the kernel of A−λI. Every linear combination of eigenvectors for the same λ is again in Eλ.
Multiplicities:1≤dimEλ≤ algebraic multiplicity of λ (MML p. 108).
A and A⊤ have the same eigenvalues, but in general not the same eigenvectors.
ProofProof (beyond the lecture)
Induction on k; k=1 holds since x1=0. Suppose ∑i=1kcixi=0. Applying A gives ∑iciλixi=0, and multiplying the first equation by λk gives ∑iciλkxi=0. Subtracting: ∑i=1k−1ci(λi−λk)xi=0. By induction ci(λi−λk)=0, and since λi=λk, ci=0 for i<k. Then ckxk=0 gives ck=0.
A matrix with fewer than n linearly independent eigenvectors is defective (MML Def. 4.13). It has an eigenvalue whose geometric multiplicity is smaller than its algebraic multiplicity.
ProofWhy
AP=[Ap1,…,Apn] and PD=[λ1p1,…,λnpn], so AP=PD iff Api=λipi for all i. P is invertible iff its columns form a basis, and then AP=PD⟺A=PDP−1 (MML pp. 116–117).
Read A=PDP−1 from right to left: P−1 changes to coordinates in the eigenbasis, D scales each coordinate by its eigenvalue, and P changes back. Powers become easy: Ak=PDkP−1.
Eigenvectors of a symmetric matrix for different eigenvalues are automatically orthogonal. Within one eigenspace of dimension >1 an orthonormal basis can be chosen with Gram–Schmidt (MML Example 4.8). This is the decomposition behind the eigenvalue test for positive semidefinite matrices.
Power iteration (MML Example 4.9, PageRank): if A has an eigenvalue λ1 with ∣λ1∣ strictly larger than all others, then the directions of x,Ax,A2x,… converge to the eigenvector of λ1 for almost every start vector x. See Intuition for why.
A linear map usually turns a vector as well as stretching it. Left: the unit vectors x (dots) and their images Ax (arrows) for A=(4123). Almost every arrow leaves the line through its dot. Right: the angle by which A turns x, over the direction of x. Its zeros are the eigenvectors p1∝(2,1)⊤ and p2∝(1,−1)⊤, the only directions that come out on their own line, scaled by λ=5 and λ=2. Directions, not vectors, are the point: the whole line through p1 is E5.
The eigenvectors give A its natural coordinate grid. A cell spanned by p1,p2 is mapped to the cell spanned by 5p1,2p2: the same directions, only stretched. In these coordinates A is the diagonal matrix D=diag(5,2), which is what A=PDP−1 says. In the standard grid the same map shears the unit square (dashed) into a slanted parallelogram. A defective matrix (e.g. the shear in Eigenvalue > Intuition) has too few eigen directions to build such a grid.
Write a start vector in the eigenbasis, x0=c1p1+c2p2. Each application of A multiplies the p1-part by 5 and the p2-part by 2, so Akx0=c15kp1+c22kp2. The p1-part wins: the direction of Akx0 turns toward p1, and the angle shrinks by about λ2/λ1=2/5 per step (right, log scale). This is power iteration, and the reason why the eigenvector of the largest eigenvalue dominates, e.g. in PageRank.