Study Vault Search notes, aliases and full text
Graph
ConceptLinear Algebra used inMSIAI-written · not yet proofread

Eigenvector

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).

Know these firstEigenvalueLinear Map

Properties#

  • Only the direction is determined: if xx is an eigenvector for λ\lambda, so is cxcx for every c≠0c \neq 0, since A(cx)=cAx=λ(cx)A(cx) = cAx = \lambda(cx) (MML p. 105). An eigenvector is therefore usually given up to scaling, e.g. normalized to ∥x∥=1\|x\| = 1. Two vectors that point in the same or the opposite direction are called collinear.
  • Eigenspaces are subspaces: EλE_\lambda is the solution set of the homogeneous linear system (A−λI)x=0(A - \lambda I)x = 0, so it is the kernel of A−λIA - \lambda I. Every linear combination of eigenvectors for the same λ\lambda is again in EλE_\lambda.
  • Multiplicities: 1≤dim⁡Eλ≤1 \le \dim E_\lambda \le algebraic multiplicity of λ\lambda (MML p. 108).
  • AA and A⊤A^\top have the same eigenvalues, but in general not the same eigenvectors.
ProofProof (beyond the lecture)

Induction on kk; k=1k = 1 holds since x1≠0x_1 \neq 0. Suppose ∑i=1kcixi=0\sum_{i=1}^k c_i x_i = 0. Applying AA gives ∑iciλixi=0\sum_i c_i\lambda_i x_i = 0, and multiplying the first equation by λk\lambda_k gives ∑iciλkxi=0\sum_i c_i\lambda_k x_i = 0. Subtracting: ∑i=1k−1ci(λi−λk)xi=0\sum_{i=1}^{k-1} c_i(\lambda_i - \lambda_k)x_i = 0. By induction ci(λi−λk)=0c_i(\lambda_i - \lambda_k) = 0, and since λi≠λk\lambda_i \neq \lambda_k, ci=0c_i = 0 for i<ki < k. Then ckxk=0c_k x_k = 0 gives ck=0c_k = 0.

A matrix with fewer than nn 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]AP = [Ap_1, \dots, Ap_n] and PD=[λ1p1,…,λnpn]PD = [\lambda_1 p_1, \dots, \lambda_n p_n], so AP=PDAP = PD iff Api=λipiAp_i = \lambda_i p_i for all ii. PP is invertible iff its columns form a basis, and then AP=PD  ⟺  A=PDP−1AP = PD \iff A = PDP^{-1} (MML pp. 116–117).

Read A=PDP−1A = PDP^{-1} from right to left: P−1P^{-1} changes to coordinates in the eigenbasis, DD scales each coordinate by its eigenvalue, and PP changes back. Powers become easy: Ak=PDkP−1A^k = PD^kP^{-1}.

Eigenvectors of a symmetric matrix for different eigenvalues are automatically orthogonal. Within one eigenspace of dimension >1> 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 AA has an eigenvalue λ1\lambda_1 with ∣λ1∣|\lambda_1| strictly larger than all others, then the directions of x,Ax,A2x,…x, Ax, A^2x, \dots converge to the eigenvector of λ1\lambda_1 for almost every start vector xx. See Intuition for why.

Example#

Intuition#

la eigenvector directions
A linear map usually turns a vector as well as stretching it. Left: the unit vectors xx (dots) and their images AxAx (arrows) for A=(4213)A = \begin{pmatrix} 4 & 2 \\ 1 & 3 \end{pmatrix}. Almost every arrow leaves the line through its dot. Right: the angle by which AA turns xx, over the direction of xx. Its zeros are the eigenvectors p1∝(2,1)⊤p_1 \propto (2, 1)^\top and p2∝(1,−1)⊤p_2 \propto (1, -1)^\top, the only directions that come out on their own line, scaled by λ=5\lambda = 5 and λ=2\lambda = 2. Directions, not vectors, are the point: the whole line through p1p_1 is E5E_5.

la eigenvector basis
The eigenvectors give AA its natural coordinate grid. A cell spanned by p1,p2p_1, p_2 is mapped to the cell spanned by 5p1,2p25p_1, 2p_2: the same directions, only stretched. In these coordinates AA is the diagonal matrix D=diag⁡(5,2)D = \operatorname{diag}(5, 2), which is what A=PDP−1A = 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.

la eigenvector power
Write a start vector in the eigenbasis, x0=c1p1+c2p2x_0 = c_1p_1 + c_2p_2. Each application of AA multiplies the p1p_1-part by 55 and the p2p_2-part by 22, so Akx0=c15kp1+c22kp2A^kx_0 = c_15^kp_1 + c_22^kp_2. The p1p_1-part wins: the direction of Akx0A^kx_0 turns toward p1p_1, and the angle shrinks by about λ2/λ1=2/5\lambda_2/\lambda_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.

Sources#

  • MML p. 105 (Def. 4.6, Def. 4.7 collinearity, non-uniqueness of eigenvectors)
  • MML p. 106 (Def. 4.10 eigenspace, Eλ=ker⁡(A−λI)E_\lambda = \ker(A - \lambda I)), MML pp. 107–108 (Examples 4.5, 4.6, Def. 4.11 geometric multiplicity)
  • MML pp. 111–112 (Thm. 4.12, Def. 4.13 defective, Thm. 4.15 spectral theorem, Example 4.8), MML pp. 113–114 (Example 4.9 PageRank)
  • MML pp. 116–117 (Def. 4.19 diagonalizable, Thm. 4.20 eigendecomposition, Fig. 4.7)
  • Use in MSI: spectral decomposition behind the eigenvalue test for PSD matrices

Search

Type In
↑ ↓ move ↵ open ⇥ next section Ctrl↵ results page Prefix # to search text only