Matrices provide a compact algebraic framework for representing and manipulating systems of linear
Equations, geometric transformations, and — at a more advanced level — quantum mechanical states and
Data structures. This topic develops the algebra of matrices and their interpretation as linear
Transformations of the plane and space.
Board Coverage
Board
Paper
Notes
AQA
Paper 1
3D transformations, eigenvalues and eigenvectors
Edexcel
FP1
2D transformations, 3×3 matrices, determinants
OCR (A)
Paper 1
2D transformations, 3×3 matrices
CIE
P1
2D and basic 3D transformations, inverses, determinants
1. Matrix Notation and Operations
Definition. An m×nmatrix is a rectangular array of numbers arranged in m rows and
n columns. We write A=(aij) where aij is the entry in row iColumn j. The
set of all m×n matrices with real entries is denoted Mm×n(R).
A matrix with a single row is a row vector, and a matrix with a single column is a column
Vector.
1.1 Matrix Addition and Scalar Multiplication
For matrices A,B∈Mm×n(R) and scalar k∈R:
(A+B)ij(kA)ij=aij+bij=k⋅aij
Matrix addition is commutative (A+B=B+A) and
Associative.
1.2 Matrix Multiplication
Definition. If A∈Mm×p(R) and
B∈Mp×n(R)The product AB∈Mm×n(R)
is defined by:
(AB)ij=k=1∑paikbkj
Matrix multiplication is associative but not commutative : AB=BA.
The n×n identity matrix In satisfies
AIn=InA=A for any
A∈Mm×n(R).
:::caution warning “apply B first, then A.” When composing transformations,
The rightmost matrix is applied first.
:::
2. Determinants
2.1 The 2×2 Determinant
Definition. For A=(acbd):
det(A)=ad−bc
2.2 The 3×3 Determinant
Definition. For
A=a11a21a31a12a22a32a13a23a33
The determinant is computed by cofactor expansion along any row or column:
Intuition. The determinant measures how a matrix scales area (in 2D) or volume (in 3D).
Composing two transformations multiplies their area/volume scaling factors, which is why
Determinants multiply.
3. Inverse Matrices
Definition. The inverse of a square matrix AWritten A−1Is the Unique
matrix satisfying:
AA−1=A−1A=I
An inverse exists if and only if det(A)=0. Such a matrix is called non-singular.
3.1 Inverse of a 2×2 Matrix
For A=(acbd) with
det(A)=ad−bc=0:
A−1=ad−bc1(d−c−ba)
3.2 The Adjugate Method (3×3)
For a 3×3 matrix A:
A−1=L◆B◆1◆RB◆◆LB◆det(A)◆RB◆adj(A)
Where the adjugate (or adjoint) matrix adj(A) is the transpose of the
cofactor matrix.
Definition. The cofactorCij of entry aij is (−1)i+j times the determinant of
The submatrix obtained by deleting row i and column j. The cofactor matrix has entries
CijAnd adj(A)=(Cij)T.
Proof that the inverse is unique
Suppose B and C are both inverses of A. Then:
B=BI=B(AC)=(BA)C=IC=C□
:::tip To verify your inverse, always check that AA−1=I. This
Catches sign errors and arithmetic mistakes immediately.
:::
A system of n linear equations in n unknowns can be written as
Ax=b where A is the coefficient matrix, x is
The column vector of unknowns, and b is the column vector of constants.
If A is non-singular, the unique solution is:
x=A−1b
4.1 Geometric Interpretation (2D)
For a 2×2 system:
det(A)=0: the two lines intersect at a unique point.
det(A)=0 and the equations are consistent: the lines are coincident (infinitely many
solutions).
det(A)=0 and the equations are inconsistent: the lines are parallel (no solutions).
A 2×2 matrix T represents a linear transformation of R2: the point
(x,y) is mapped to (x′,y′) where
(x′y′)=T(xy).
5.1 Reflection
Reflection in the x-axis:(100−1)
Reflection in the y-axis:(−1001)
Reflection in the line y=x:(0110)
Reflection in the line y=tanθx:
R=(cos2θsin2θsin2θ−cos2θ)
Proof of the reflection matrix
The reflection of a vector in a line through the origin making angle θ with the x-axis can
Be decomposed: first rotate by −θ to align the mirror with the x-axis, reflect in the
x-axis, then rotate back by θ.
Rotation anticlockwise by angle θ about the origin:
Rθ=(cosθsinθ−sinθcosθ)
Note: det(Rθ)=cos2θ+sin2θ=1.
5.3 Enlargement
Enlargement by scale factor k about the origin:(k00k)
det=k2Confirming the area is scaled by k2.
5.4 Shear
Shear parallel to the x-axis by factor k:(10k1)
Shear parallel to the y-axis by factor k:(1k01)
Note: det=1 for shears, so area is preserved.
5.5 Combining Transformations
If transformation A is followed by transformation BThe combined Transformation
is represented by BA (rightmost applied first).
:::caution Order matters. A rotation followed by a reflection generally produces a different result
From a reflection followed by a rotation. The combined matrix is BA (not AB)
When A is applied first.
:::
:::tip To find the matrix of a combined transformation, multiply the matrices in reverse order of
Application. If the question says “reflect then rotate,” compute
Rrot×Rref.
:::
Worked Example: Combined transformation
Find the matrix representing a rotation of 90∘ anticlockwise about the origin followed by a
Reflection in the line y=x.
Rotation by 90∘: R=(01−10)
Reflection in y=x: S=(0110)
Combined (reflection applied after rotation):
SR=(0110)(01−10)=(100−1)
This is a reflection in the x-axis. det=−1Consistent with an orientation-reversing
Transformation.
6. Matrix Transformations in 3D
:::info 3D transformations are required by AQA and appear on CIE P1. Edexcel and OCR focus primarily
On 2D but may include basic 3×3 determinant and inverse calculations.
:::
6.1 Rotations in 3D
Rotation about the x-axis by angle θ:
Rx=1000cosθsinθ0−sinθcosθ
Rotation about the y-axis by angle θ:
Ry=cosθ0−sinθ010sinθ0cosθ
Rotation about the z-axis by angle θ:
Rz=cosθsinθ0−sinθcosθ0001
Each has det=1 and represents a rigid motion preserving distances and orientation.
6.2 Reflections in 3D
Reflection in the plane x=0 (the yz-plane):−100010001
Reflection in the plane y=0 (the xz-plane):1000−10001
Reflection in the plane z=0 (the xy-plane):10001000−1
Each has det=−1Confirming orientation reversal.
7. Eigenvalues and Eigenvectors
:::info info require this topic at A Level.
:::
Definition. Let A be an n×n matrix. A scalar λ is an eigenvalue of
A if there exists a non-zero vector v such that:
Av=λv
The vector v is called an eigenvector corresponding to λ.
7.1 The Characteristic Equation
Av=λv⟺(A−λI)v=0.
For a non-trivial solution (v=0), we require
det(A−λI)=0.
det(A−λI)=0
This is the characteristic equation of A. Its roots are the eigenvalues.
7.2 Finding Eigenvectors
For each eigenvalue λSolve (A−λI)v=0 by Row
reduction.
7.3 Diagonalisation
Definition. A matrix A is diagonalisable if there exists an invertible matrix
P and a diagonal matrix D such that:
A=PDP−1
The columns of P are the eigenvectors of AAnd the diagonal entries of
D are the corresponding eigenvalues.
A matrix is diagonalisable if and only if it has n linearly independent eigenvectors (always true
For n distinct eigenvalues).
Intuition. Diagonalisation changes to a coordinate system where the transformation acts
Independently on each axis (stretching by eigenvalues). In this basis, the matrix takes its simplest
Possible form.
Worked Example: Eigenvalues and eigenvectors
Find the eigenvalues and eigenvectors of
A=(4213).
Characteristic equation:
det(4−λ213−λ)=0
(4−λ)(3−λ)−2=λ2−7λ+10=(λ−5)(λ−2)=0
Eigenvalues: λ1=5, λ2=2.
For λ1=5:(−121−2)(xy)=(00)⟹−x+y=0⟹y=x.
Eigenvector: v1=(11) (or any non-zero scalar
Multiple).
Problem 4. The triangle with vertices (0, 0)$$(2, 0)$$(0, 1) is transformed by the matrix
T=(3012). Find the coordinates of the vertices of
The image, and verify that the area scales by ∣det(T)∣.
Hint
Apply T to each vertex. The original triangle has area 1.
Problem 6. Find the single 2×2 matrix that represents an enlargement by scale factor 2
About the origin followed by a rotation of 90∘ anticlockwise. Show that this is equivalent to
A single rotation of 90∘ combined with an enlargement by factor 2.
Hint
The enlargement matrix is 2I and the rotation is Rπ/2. Since
2I commutes with all matrices, the order doesn’t matter.
Answer
Enlargement by 2: E=(2002).
Rotation by 90∘: R=(01−10).
Enlargement then rotation:
RE=(01−10)(2002)=(02−20).
Rotation then enlargement:
ER=(2002)(01−10)=(02−20).
Both give the same result:
(02−20)=2(01−10).
This is a rotation by 90∘ combined with an enlargement by factor 2, and the order is
Irrelevant because scalar multiples of the identity commute with all matrices.
Problem 7. AQA only. Find the eigenvalues and corresponding eigenvectors of
A=(5142)And write down a matrix P And
diagonal matrix D such that A=PDP−1.
Problem 9. Find the 3×3 matrix representing a rotation of 90∘ anticlockwise
About the z-axis. Verify that this matrix has determinant 1 and that it maps
100 to 010.
The x-axis is correctly rotated to the y-axis by a 90∘ anticlockwise rotation about z.
Problem 10. Prove that if A and B are non-singular n×n matrices,
Then (AB)−1=B−1A−1.
Hint
Show that B−1A−1 satisfies the definition of the inverse of AB
And invoke uniqueness.
Answer
We need to show that (AB)(B−1A−1)=I and
(B−1A−1)(AB)=I.
(AB)(B−1A−1)=A(BB−1)A−1=AIA−1=AA−1=I
(B−1A−1)(AB)=B−1(A−1A)B=B−1IB=B−1B=I
Since the inverse is unique, (AB)−1=B−1A−1. □
Intuition. The order reverses, just like putting on and taking off socks and shoes. To undo “A
Then B,” you must undo B first, then undo A.
:::
8. Advanced Worked Examples
Example 8.1: Diagonalisation of a 3×3 matrix
Problem. Diagonalise A=200121003.
Solution. Find eigenvalues: det(A−λI)=0.
det2−λ0012−λ1003−λ=(2−λ)2(3−λ)=0
λ1=2 (repeated), λ2=3.
For λ=2: (A−2I)v=0 gives
000101001v=0So v2=0 and
v3=0With v1 free. Only one eigenvector: (1,0,0). Since the geometric multiplicity (1) is
Less than the algebraic multiplicity (2), A is not diagonalisable.
Example 8.2: Finding An using Cayley—Hamilton
Problem. For A=(1023)Find A5.
Solution. By Cayley—Hamilton, A2−4A+3I=OSo
A2=4A−3I.
Problem.A=(3012). Find all invariant points and
Invariant lines of the transformation x↦Ax.
Solution.Invariant points:Ax=x⟹(A−I)x=0.
(2011)x=0⟹x2=0,2x1=0
Only the origin (0,0) is an invariant point.
Invariant lines through the origin: These are the eigenspaces. Eigenvalues:
(3−λ)(2−λ)=0So λ=3 and λ=2.
For λ=3: (001−1)v=0⟹v2=0.
Line: y=0 (the x-axis).
For λ=2:
(1010)v=0⟹v1+v2=0. Line:
y=−x.
Example 8.4: Determinant as a scaling factor
Problem. The matrix T=(2−113) represents a
Transformation. A triangle has vertices (0,0), (1,0), (0,1). Find the area of its image.
Solution.det(T)=6−(−1)=7.
Original area =21. Image area
=∣det(T)∣×original area=7×21=3.5.
Example 8.5: Commutator and non-commuting matrices
Problem. For A=(0010) and
B=(0100)Compute the commutator
[A,B]=AB−BA.
Solution.AB=(0010)(0100)=(1000).
BA=(0100)(0010)=(0001).
[A,B]=(100−1)
Since [A,B]=O, A and B do not commute.
Example 8.6: Matrix representing successive transformations
Problem. Transformation R is a reflection in the line y=x3. Transformation S is a
Rotation by 90° anticlockwise about the origin. Find the matrix representing RS.
Solution. The line y=x3 makes angle 60° with the x-axis.
det(RS)=−43−41=−1 and tr(RS)=0Confirming this is a
Reflection.
Example 8.7: Finding the inverse of a 3×3 matrix
Problem. Find A−1 where
A=1010112−10.
Solution.det(A)=1(0+1)−0+2(0−1)=1−2=−1=0.
A−1=−111−1−12−2−1−211=−111−2212−1−1
9. Common Pitfalls
Pitfall
Correct Approach
Assuming all matrices are diagonalisable
Check geometric multiplicity equals algebraic multiplicity for each eigenvalue
Forgetting that (AB)−1=B−1A−1
The order reverses
Computing det(A+B)=det(A)+det(B)
, det(A+B)=det(A)+det(B)
Mixing up row and column operations
Column operations change the determinant differently from row operations
10. Additional Exam-Style Questions
Question 8
Find the eigenvalues and eigenvectors of A=(42−11). Hence
Write down P and D such that
P−1AP=D.
Solution
det(A−λI)=(4−λ)(1−λ)+2=λ2−5λ+6=0.
λ=2,3.
λ=2:
(22−1−1)v=0⟹2v1−v2=0⟹v=(1,2).
λ=3:
(12−1−2)v=0⟹v1=v2⟹v=(1,1).
P=(1211)D=(2003).
Question 9
Prove that if λ is an eigenvalue of A with eigenvector vThen
λ2 is an eigenvalue of A2 with the same eigenvector.
Solution
Av=λv.
A2v=A(Av)=A(λv)=λ(Av)=λ(λv)=λ2v.
Therefore λ2 is an eigenvalue of A2 with eigenvector v.
■
Question 10
The transformation represented by M=(acbd) maps the unit
Square to a parallelogram of area 6. Given a+d=5 and ad−bc=6Find the eigenvalues of
M.
Solution
The characteristic equation:
λ2−(a+d)λ+det(M)=λ2−5λ+6=0.
(λ−2)(λ−3)=0.
λ=2 and λ=3
11. Connections to Other Topics
11.1 Matrices and complex numbers
Complex eigenvalues lead to rotation-scaling transformations. See
Complex Numbers.
11.2 Matrices and vectors
The cross product can be written as a matrix multiplication. See
Vectors in 3D.
11.3 Matrices and further algebra
Cayley—Hamilton connects matrices to polynomial algebra. See
Further Algebra.
12. Key Results Summary
| Result | Formula/Condition |
| -------------------- | --------------------------------------------------------------------------- | ---------------- | ------------------------------------ |
| Invertibility | A is invertible ⟺det(A)=0 |
| (AB)−1 | B−1A−1 |
| (AB)T | BTAT |
| det(AB) | det(A)det(B) |
| Trace of product | tr(AB)=tr(BA) |
| Cayley—Hamilton | A satisfies its own characteristic equation |
| Diagonalisability | All eigenvalues must have geometric multiplicity = algebraic multiplicity |
| Area scaling | ∣det(T)∣× original area = image area |
13. Further Exam-Style Questions
Question 11
Find the matrix representing a stretch of scale factor 3 parallel to the y-axis followed by a
Reflection in the x-axis.
Solution
Stretch: S=(1003). Reflection:
R=(100−1).
Combined:
RS=(100−1)(1003)=(100−3).
(100−3)
Question 12
Prove thatdet(AT)=det(A) for any square matrix A.
Solution
The determinant can be computed by cofactor expansion along any row or column. Expanding
det(A) along row i and det(AT) along column i (which is row i of
A) gives the same expression, since the cofactors are the same.
Formally: this follows from the permutation definition of the determinant. ■
14. Advanced Topics
14.1 Eigenvalues and the characteristic polynomial
The characteristic polynomial of A is p(λ)=det(A−λI).
Properties:
The sum of eigenvalues equals the trace: ∑λi=tr(A)
The product of eigenvalues equals the determinant: ∏λi=det(A)
A is invertible iff no eigenvalue is zero
14.2 Jordan normal form
Not all matrices are diagonalisable. The Jordan normal form is a generalisation where the diagonal
Matrix D may have 1s on the superdiagonal (Jordan blocks).
For example, if A has a repeated eigenvalue λ with only one eigenvector:
P−1AP=(λ01λ)
14.3 Orthogonal diagonalisation
A real symmetric matrix A can always be diagonalised by an orthogonal matrix:
A=QDQT where QT=Q−1.
The Spectral Theorem states that A=∑i=1nλiqiqiT
Where λi are eigenvalues and qi are orthonormal eigenvectors.
14.4 Matrix norms
The Frobenius norm:
∥A∥F=◆LB◆∑i,jaij2◆RB◆=◆LB◆tr(ATA)◆RB◆.
The spectral norm: ∥A∥2=σmax (largest singular value).
15. Further Exam-Style Questions
Question 14
Find the eigenvalues and eigenvectors of
A=111111111.
Prove that similar matrices have the same eigenvalues.
Solution
If B=P−1APThen:
det(B−λI)=det(P−1AP−λI)=det(P−1(A−λI)P)
=det(P−1)det(A−λI)det(P)=det(A−λI).
Since the characteristic polynomials are identical, the eigenvalues are the same. ■
16. Further Advanced Topics
16.1 LU decomposition
Any square matrix A can be decomposed as A=LU where
L is lower triangular and U is upper triangular. This is used for efficient
Numerical solution of systems Ax=b.
16.2 The Cayley—Hamilton theorem — applications
Since A satisfies p(A)=O where p is the characteristic
Polynomial:
A−1 can be computed from An terms
An for large n can be reduced using the recurrence
16.3 Singular Value Decomposition (SVD)
Any m×n matrix A can be written as
A=UΣVT where U and V are
Orthogonal and Σ is diagonal with non-negative singular values.
16.4 Positive definite matrices
A symmetric matrix A is positive definite if xTAx>0 for
All x=0.
Equivalent conditions: all eigenvalues positive, all leading principal minors positive.
17. Further Exam-Style Questions
Question 16
Find det(A) where A=147258360.
Solution
det(A)=1(0−48)−2(0−42)+3(32−35)=−48+84−9=27.
Question 17
Prove that the trace of a matrix equals the sum of its eigenvalues (with multiplicity).
Solution
The characteristic polynomial is
p(λ)=(−1)n[λn−(trA)λn−1+⋯+(−1)ndet(A)].
By Vieta’s formulae, the coefficient of λn−1 equals
−(λ1+λ2+⋯+λn).