Ten problems on eigenvalues and eigenvectors worked by hand: 2×2 and 3×3 matrices, the characteristic polynomial, a rotation with complex eigenvalues, a shear with too few eigenvectors, diagonalisation and powers, orthogonal eigenvectors of symmetric matrices, power iteration and singular values, with worked solutions and the classic sign and null-space mistakes.
Before you start
Every eigenvalue problem on this page comes down to two computations: a determinant that gives the eigenvalues, and a small linear system for each one that gives the eigenvectors. The errors come from the same few places: a sign in ,A−λI, an eigenvector read off a row instead of solved for, a repeated eigenvalue assumed to bring two eigenvectors. These ten problems start with 2×2 matrices, meet a rotation and a shear that break the easy pattern, and end with the three uses machine learning makes of eigenvectors: powers of a matrix, power iteration and singular values.
A is a real n×n matrix and I the n×n identity. A number λ is an eigenvalue of ,A, with eigenvector ,v, when Av=λv and .v=0.
Av=λv is the same as .(A−λI)v=0. That has a nonzero solution exactly when A−λI is singular, that is when .det(A−λI)=0.det(A−λI) is a polynomial of degree n in ,λ, the characteristic polynomial, and its roots are the eigenvalues.
An eigenvector is any nonzero solution of ,(A−λI)v=0, so the answers below give one representative each; any nonzero multiple is also correct.
Vectors are columns; (1,−1)⊤ is the column with entries 1 and ,−1, written on one line to save space.
For 2×2 matrices: ,det[acbd]=ad−bc, and when that is not ,0,.[acbd]−1=ad−bc1[d−c−ba].,trA, the trace, is the sum of the diagonal entries.
P is a matrix whose columns are eigenvectors of .A.D (Problem 7) and Λ (Problem 8) are both diagonal matrices of eigenvalues, listed in the same order as the columns of ;P;diag(a,b) is the diagonal matrix with a and b on the diagonal. Λ is used when the eigenvectors are orthonormal, and then the eigenvector matrix is written :Q: an orthogonal matrix, meaning ,Q⊤Q=I, so .Q−1=Q⊤.
i is the imaginary unit, .i2=−1.σ1,σ2 are singular values (Problem 10).
Problems
·
Find the eigenvalues and eigenvectors of .A=[2112].
·
Find the characteristic polynomial, eigenvalues and eigenvectors of .A=[4213].
··
Find the eigenvalues of the rotation .R=[01−10]. Does it have a real eigenvector? What does that mean geometrically?
··
Find the eigenvalues and eigenvectors of the shear .A=[1011]. Can it be diagonalised?
··
For any 2×2 matrix, show that the trace is the sum of the eigenvalues and the determinant is their product. Verify on Problem 2.
···
Find the eigenvalues and eigenvectors of .A=210121012.
··
Diagonalise the matrix of Problem 2 as A=PDP−1 and use it to compute .A5.
··
Show that the eigenvectors of Problem 1 are orthogonal, and write A=QΛQ⊤ with Q orthogonal. Why does this always happen for symmetric matrices?
···
Apply A=[2112] repeatedly to .x0=(1,0)⊤. Find Akx0 in closed form and say what direction it approaches and how fast.
···
For ,A=[3405], find the eigenvalues of A⊤A and hence the singular values of .A. Check them against .detA.
,Q=21[111−1],,Λ=diag(3,1),;A=QΛQ⊤; for symmetric ,A,,λ1v1⊤v2=(Av1)⊤v2=v1⊤Av2=λ2v1⊤v2, so distinct eigenvalues force v1⊤v2=0
;Akx0=23k(1,1)⊤+21(1,−1)⊤; the direction approaches ,(1,1)/2, the error shrinking like (1/3)k — the ratio of the two eigenvalues
A⊤A=[25202025] with eigenvalues 45 and ;5; singular values 35 and ;5; their product 15=∣detA∣
Worked solutions
Problem 1
Find the eigenvalues and eigenvectors of .A=[2112].
,A−λI=[2−λ112−λ], so .det(A−λI)=(2−λ)2−1.Subtract λ from each diagonal entry, then .ad−bc.
,(2−λ)2=1, so 2−λ=±1 and λ=1 or .λ=3.The eigenvalues are the roots of the characteristic polynomial; this one is a square minus ,1, so take square roots rather than expanding.
:λ=3:,A−3I=[−111−1], and with ,v=(x,y)⊤,(A−3I)v=0 says ,−x+y=0, so .v=(1,1)⊤.The eigenvectors for λ=3 are the null space of .A−3I. Its two rows are multiples of each other, as they must be when the determinant is ,0, so one equation remains; choose .x=1.
:λ=1:,A−I=[1111], and (A−I)v=0 says ,x+y=0, so .v=(1,−1)⊤.Null space of ;A−I; the equation x+y=0 fixes v up to scale.
λ=3 with ;v=(1,1)⊤;λ=1 with v=(1,−1)⊤Check: A(1,1)⊤=(3,3)⊤ and .A(1,−1)⊤=(1,−1)⊤.
Problem 2
Find the characteristic polynomial, eigenvalues and eigenvectors of .A=[4213].
.det(A−λI)=(4−λ)(3−λ)−1⋅2.ad−bc for .A−λI=[4−λ213−λ].
.(4−λ)(3−λ)−2=λ2−7λ+12−2=λ2−7λ+10.Expand the product and collect terms.
.λ2−7λ+10=(λ−5)(λ−2).Two numbers with product 10 and sum .7.
:λ=5:A−5I=[−121−2] gives −x+y=0 for ,v=(x,y)⊤, so .v=(1,1)⊤.Null space of ;A−5I; the second row is −2 times the first, so it adds nothing.
:λ=2:A−2I=[2211] gives ,2x+y=0, so .v=(1,−2)⊤.Null space of ,A−2I, with .x=1.
;λ2−7λ+10=(λ−5)(λ−2);,λ=5,;v=(1,1)⊤;,λ=2,v=(1,−2)⊤Check: A(1,1)⊤=(5,5)⊤ and .A(1,−2)⊤=(2,−4)⊤. The eigenvectors are not perpendicular; A is not symmetric, and Problem 8 shows why that matters.
Problem 3
Find the eigenvalues of the rotation .R=[01−10]. Does it have a real eigenvector? What does that mean geometrically?
λ2+1=0 has no real root; .λ=±i.λ2≥0 for every real ,λ, so ;λ2+1≥1; the two roots are complex.
No real eigenvector.A real eigenvector v would make Rv=λv with both sides real, so λ would be real; there is no real eigenvalue. The eigenvectors are complex: λ=i has ,v=(1,−i)⊤, since .R(1,−i)⊤=(i,1)⊤=i(1,−i)⊤.
,R(x,y)⊤=(−y,x)⊤, which is v=(x,y)⊤ turned 90° anticlockwise.v⋅Rv=−xy+yx=0 and ,∥Rv∥=∥v∥, so for every nonzero real ,v,Rv is a nonzero vector perpendicular to ,v, off the line through .v. An eigenvector is exactly a vector that stays on its own line.
,λ2+1=0,;λ=±i; no real eigenvector: a rotation by 90° sends every real direction somewhere elseSteps 2 to 4.
Problem 4
Find the eigenvalues and eigenvectors of the shear .A=[1011]. Can it be diagonalised?
.det(A−λI)=(1−λ)2−1⋅0=(λ−1)2.A−λI is upper triangular, so its determinant is the product of the diagonal entries.
,λ=1, a double root.The only root, counted twice: its algebraic multiplicity, the number of times it is a root, is .2.
A−I=[0010] has one nonzero row, so .rank(A−I)=1.The rank decides how many independent eigenvectors λ=1 has: the null space of a 2×2 matrix has dimension .2−rank.
For ,v=(x,y)⊤,(A−I)v=0 says ,y=0, so every eigenvector is ,(x,0)⊤, a multiple of .(1,0)⊤.The null space has dimension :2−1=1: the geometric multiplicity, the number of independent eigenvectors, is .1.
Not diagonalisable.A=PDP−1 needs two independent eigenvectors for the columns of an invertible ,P, and there is only one direction. Directly: D would have to be ,diag(1,1)=I, and then .PIP−1=I=A.
,(λ−1)2,λ=1 (double); A−I=[0010] has rank ,1, so the only eigenvectors are multiples of ;(1,0)⊤; not diagonalisableSteps 1 to 5.
Problem 5
For any 2×2 matrix, show that the trace is the sum of the eigenvalues and the determinant is their product. Verify on Problem 2.
Let .A=[acbd]. Then .det(A−λI)=(a−λ)(d−λ)−bc.ad−bc for .A−λI.
.(a−λ)(d−λ)−bc=λ2−(a+d)λ+(ad−bc).Expand and collect powers of .λ.
.det(A−λI)=(λ−λ1)(λ−λ2)=λ2−(λ1+λ2)λ+λ1λ2.λ1 and λ2 are the two roots, complex or repeated if need be, and a quadratic whose λ2 coefficient is 1 is the product of (λ−root) over its roots.
a+d=λ1+λ2 and .ad−bc=λ1λ2.Steps 2 and 3 are the same polynomial, so their coefficients of λ and their constant terms agree.
,trA=λ1+λ2,;detA=λ1λ2; Problem 2: ,7=5+2,10=5⋅2For Problem 2, trA=4+3=7 and .detA=4⋅3−1⋅2=10. This is a quick check on any eigenvalue computation.
Problem 6
Find the eigenvalues and eigenvectors of .A=210121012.
Write ,t=2−λ, so the diagonal entries of A−λI are all t and the off-diagonal entries are those of .A. Each eigenvector is written .v=(x,y,z)⊤.
.det(A−λI)=tdet[t11t]−1⋅det[101t]+0.Cofactor expansion along the first row, :(t,1,0): each entry times the determinant left after deleting its row and column, with signs .+,−,+.
.t(t2−1)−t=t3−2t=t(t2−2).The two 2×2 determinants are t2−1 and .t−0.
.det(A−λI)=(2−λ)((2−λ)2−2).Substitute t=2−λ back.
,λ=2, or ,2−λ=±2, that is .λ=2∓2.Set each factor to .0. Check with the trace: .2+(2−2)+(2+2)=6=trA.
λ=2 ():t=0): the rows of A−2I give ,y=0,,x+z=0,,y=0, so .v=(1,0,−1)⊤.The null space of ,A−2I, whose rows are ,(0,1,0),,(1,0,1),.(0,1,0).
λ=2+2 ():t=−2): the rows give ,−2x+y=0,,x−2y+z=0,.y−2z=0.The null space of ,A−(2+2)I, whose diagonal entries are .−2.
y=2x and ,y=2z, so ;z=x; with ,x=1,.v=(1,2,1)⊤.The first and third equations; the middle one then holds, ,1−2+1=0, as it must because the determinant is .0.
λ=2−2 ():t=2): the same equations with 2 in place of −2 give ,y=−2x=−2z, so .v=(1,−2,1)⊤.Only the sign of t changed, so only the sign of y changes.
λ=2 with ;(1,0,−1)⊤;λ=2±2 with (1,±2,1)⊤Check: .A(1,2,1)⊤=(2+2,2+22,2+2)⊤=(2+2)(1,2,1)⊤.A is symmetric, and the three eigenvectors are mutually perpendicular (Problem 8).
Problem 7
Diagonalise the matrix of Problem 2 as A=PDP−1 and use it to compute .A5.
,P=[111−2],.D=diag(5,2).The columns of P are the eigenvectors from Problem 2, and D lists their eigenvalues in the same order. Then AP=PD column by column, which is Av=λv twice.
,detP=1⋅(−2)−1⋅1=−3, so .P−1=−31[−2−1−11]=[2/31/31/3−1/3].The 2×2 inverse: swap the diagonal, negate the off-diagonal, divide by the determinant. detP=0 because the two eigenvectors are independent.
PD=[552−4] and .PDP−1=[4213]=A.Check the diagonalisation by multiplying it out: the first row of PDP−1 is .(10/3+2/3,5/3−2/3)=(4,1).
.A5=PDP−1PDP−1⋯PDP−1=PD5P−1.Each inner P−1P is ,I, so the five copies collapse. Powers of a diagonal matrix are powers of its entries: .D5=diag(3125,32).
.PD5=[3125312532−64].Multiplying by a diagonal matrix on the right scales the columns of P by 55 and .25.
.PD5P−1=31[6250+326250−643125−323125+64]=[2094206210311063].Row times column with .P−1=31[211−1].
,P=[111−2],,D=diag(5,2),;P−1=−31[−2−1−11];A5=PD5P−1=[2094206210311063]Check: (1,1)⊤ is an eigenvector with eigenvalue ,5, so ,A5(1,1)⊤=3125(1,1)⊤, and both rows do sum to .3125.
Problem 8
Show that the eigenvectors of Problem 1 are orthogonal, and write A=QΛQ⊤ with Q orthogonal. Why does this always happen for symmetric matrices?
.(1,1)⊤⋅(1,−1)⊤=1−1=0.Orthogonal means a zero dot product.
,Q=21[111−1],.Λ=diag(3,1).Divide each eigenvector by its length 2 so the columns are unit vectors; they are still eigenvectors. Λ lists the eigenvalues in column order.
.Q⊤Q=21[2002]=I.The columns are orthogonal unit vectors, so Q is orthogonal and .Q−1=Q⊤.
.QΛQ⊤=21[331−1][111−1]=21[4224]=A.A=QΛQ−1 as in Problem 7, with Q−1 replaced by ;Q⊤; multiplying it out confirms it.
For symmetric A with ,Av1=λ1v1,:Av2=λ2v2:.λ1v1⊤v2=(Av1)⊤v2=v1⊤A⊤v2=v1⊤Av2=λ2v1⊤v2.Move A across the dot product: ,(Av1)⊤=v1⊤A⊤, and A⊤=A is where symmetry is used.
,(λ1−λ2)v1⊤v2=0, so v1⊤v2=0 when .λ1=λ2.Subtract the two ends of step 5; a product is 0 only if a factor is. A repeated eigenvalue of a symmetric matrix has as many independent eigenvectors as its multiplicity (the spectral theorem), and they can be chosen orthogonal, so A=QΛQ⊤ holds for every real symmetric matrix.
,Q=21[111−1],,Λ=diag(3,1),;A=QΛQ⊤; for symmetric ,A,,λ1v1⊤v2=(Av1)⊤v2=v1⊤Av2=λ2v1⊤v2, so distinct eigenvalues force v1⊤v2=0Steps 1 to 6.
Problem 9
Apply A=[2112] repeatedly to .x0=(1,0)⊤. Find Akx0 in closed form and say what direction it approaches and how fast.
.x0=21(1,1)⊤+21(1,−1)⊤.Write x0 in the eigenvectors from Problem 1, where A acts by scaling; the coefficients solve ,c1+c2=1,.c1−c2=0.
Ak(1,1)⊤=3k(1,1)⊤ and .Ak(1,−1)⊤=1k(1,−1)⊤=(1,−1)⊤.Each application of A multiplies an eigenvector by its eigenvalue, so k applications multiply it by the k-th power.
.Akx0=23k(1,1)⊤+21(1,−1)⊤.Ak is linear, so apply step 2 to each term of step 1. Check: k=1 gives ,(2,1)⊤=Ax0, and k=2 gives .(5,4)⊤.
.Akx0=23k((1,1)⊤+3−k(1,−1)⊤).Factor out the growing term; the positive scale factor does not change the direction.
,tanθk=3−k, where θk is the angle between Akx0 and .(1,1)⊤.In step 4 the two eigenvectors are perpendicular and have the same length ,2, so the ratio of their coefficients, ,3−k, is the tangent of that angle, and it goes to .0.
;Akx0=23k(1,1)⊤+21(1,−1)⊤; the direction approaches ,(1,1)/2, the error shrinking like (1/3)k — the ratio of the two eigenvaluesThis is power iteration: repeated multiplication (with a rescaling each step to keep the length finite, since it grows like )3k) finds the eigenvector of the largest eigenvalue, at a rate set by .λ2/λ1=1/3.
Problem 10
For ,A=[3405], find the eigenvalues of A⊤A and hence the singular values of .A. Check them against .detA.
.A⊤A=[3045][3405]=[9+16202025]=[25202025].Row times column; the result is symmetric, as A⊤A always is.
,det(A⊤A−λI)=(25−λ)2−400, so 25−λ=±20 and λ=45 or .5.The same shape as Problem 1: a square minus a constant. The eigenvectors are (1,1)⊤ and (1,−1)⊤ again.
σ1=45=35 and .σ2=5.The singular values are the square roots of the eigenvalues of :A⊤A: if A=UΣV⊤ with ,U,V orthogonal and ,Σ=diag(σ1,σ2), then ,A⊤A=VΣ2V⊤, an orthogonal diagonalisation as in Problem 8 with .Λ=Σ2. The eigenvalues are never negative, since .v⊤A⊤Av=∥Av∥2≥0.
,σ1σ2=35⋅5=15, and .detA=3⋅5−0⋅4=15.det(A⊤A)=(detA)2 is the product of its eigenvalues σ12σ22 (Problem 5), so ;σ1σ2=∣detA∣; here .45⋅5=225=152.
A⊤A=[25202025] with eigenvalues 45 and ;5; singular values 35 and ;5; their product 15=∣detA∣The eigenvalues of A itself are 3 and ,5, its diagonal entries since it is triangular: same product, but not the singular values. They coincide, up to sign, only for special matrices such as symmetric ones.
Where this goes wrong
1. Wrong sign in the characteristic matrix
Av=λv has λ on the right, and moving it across is where the sign slips.
Av=λv for A=[2112]Right so far: the eigen-equation of Problem 1.
“Collect everything on one side: A and λI together.”The step that causes the mistake: the two terms are gathered without tracking the sign.
,det(A+λI)=(2+λ)2−1=0, so λ=−1 or λ=−3Av=λv rearranges to ,Av−λv=0, that is .(A−λI)v=0. With the plus sign every root comes out negated: −1 and −3 instead of 1 and .3. The trace catches it: the eigenvalues must sum to trA=4 (Problem 5), and these sum to .−4.
2. Reading the eigenvector off a row of A − λI
The rows of A−λI are the first vectors on the page after the eigenvalue is found, and it is tempting to take one as the answer.
λ=3 in Problem 1, and A−3I=[−111−1]Right so far.
“The eigenvector comes from ,A−3I, so read it off the first row.”The step that causes the mistake: the row is taken as the vector instead of as an equation the vector must satisfy.
v=(−1,1)⊤The eigenvector is in the null space of ,A−3I, which is orthogonal to its rows: the row (−1,1) says ,−x+y=0, and the solution (1,1)⊤ is perpendicular to it. Here the wrong vector differs from the right one only in the sign of one entry, which is why the habit survives, but :A(−1,1)⊤=(−1,1)⊤: it is the eigenvector for ,λ=1, not .3. For a 2×2 row ,(a,b), the correct move is to swap the entries and negate one, giving .(b,−a).
3. Assuming a double eigenvalue gives two eigenvectors
Two eigenvalues usually bring two eigenvectors, and a double eigenvalue looks like two eigenvalues.
The shear A=[1011] has det(A−λI)=(λ−1)2Right so far: Problem 4.
“A 2×2 matrix has two eigenvalues, so it has two eigenvectors.”The count that causes the mistake: it counts roots, not independent solutions.
“λ=1 twice, so R2 is spanned by eigenvectors and .A=PDP−1.”The algebraic multiplicity is 2 but the geometric multiplicity is :1:A−I has rank ,1, so its null space holds only multiples of (1,0)⊤ (Problem 4). The rank computation is the test; without it the claim would force .A=PIP−1=I.
4. P and P⁻¹ swapped
P and P−1 appear on opposite sides of ,D, and which side is which is easy to forget.
,P=[111−2],D=diag(5,2) for the matrix of Problem 2Right so far: Problem 7.
“A is D with P on one side and P−1 on the other.”The half-memory that causes the mistake: it does not say which side.
A5=P−1D5PA=PDP−1 because :AP=PD:P holds the eigenvectors as columns and multiplies on the left. The swapped product is not ;A5; here it is ,[2094103120621063], the transpose of the right answer — only because this P happens to be symmetric; in general the swapped product is just a different matrix. The eigenvector check catches it: A5(1,1)⊤ must be ,3125(1,1)⊤, but the first row of this matrix sums to .4156.
5. Complex eigenvalues read as no eigenvalues
In school algebra a negative discriminant means “no solution”, and the habit carries over.
The rotation R=[01−10] has det(R−λI)=λ2+1Right so far: Problem 3.
“λ2+1=0 has no solution.”The habit that causes the mistake: no real solution is read as no solution.
“The rotation has no eigenvalues.”The characteristic polynomial of an n×n matrix always has n roots over ,C, counted with multiplicity; here they are .±i. What the rotation lacks is a real eigenvector, because no real direction stays on its own line under a 90° turn.