Unit 1 · 19 min read

Matrices and Determinants

What the syllabus asks for 5
  • Matrices and operations
  • Determinants
  • Inverse of a matrix
  • The determinant as an area
  • Eight uses of matrices in computing

Matrices store tables of numbers and describe linear change: a rotation of an image, a system of equations, a network's connections, a neural-network layer. This unit builds the algebra of matrices, determinants and inverses, the foundation of the whole subject.

What this unit gives you
  • Matrix types and operations (sum, product, transpose) with their laws.
  • Determinants: expansion, properties, shortcuts.
  • Inverse by adjoint and by row operations; properties of the inverse.

1.1 Matrices and operations

An m \times n matrix has m rows and n columns; a_{ij} is the entry in row i, column j.

Type Description
Square m = n
Diagonal non-zero entries only on the main diagonal
Identity I_{n} diagonal of 1s; AI = I A = A
Zero all entries 0
Symmetric A^{T} = A; skew-symmetric: A^{T} = - A
Upper/lower triangular zeros below/above the diagonal
Orthogonal A^{T} A = I

Sum and scalar multiple are entry by entry (same size). The transpose A^{T} swaps rows and columns; (A B)^{T} = B^{T} A^{T}.

Product. If A is m \times n and B is n \times p, AB is m \times p with (A B)_{ij} = \sum_{k} a_{ik} b_{kj} (row of A times column of B). The inner sizes must agree.

Worked example

A = \begin{bmatrix} 2 & 1 \\ 4 & 3 \end{bmatrix} ,\ B = \begin{bmatrix} 1 & 0 \\ 2 & 5 \end{bmatrix}.

A B = \begin{bmatrix} 2 \cdot 1 + 1 \cdot 2 & 2 \cdot 0 + 1 \cdot 5 \\ 4 \cdot 1 + 3 \cdot 2 & 4 \cdot 0 + 3 \cdot 5 \end{bmatrix} = \begin{bmatrix} 4 & 5 \\ 10 & 15 \end{bmatrix}.

B A = \begin{bmatrix} 1 \cdot 2 + 0 \cdot 4 & 1 \cdot 1 + 0 \cdot 3 \\ 2 \cdot 2 + 5 \cdot 4 & 2 \cdot 1 + 5 \cdot 3 \end{bmatrix} = \begin{bmatrix} 2 & 1 \\ 24 & 17 \end{bmatrix}.

So A B \neq B A: matrix multiplication is not commutative. It is associative, (AB)C = A (B C), and distributive, A (B + C) = A B + A C.

Warning signs: A B = 0 does not imply A = 0 or B = 0 ;\ (A + B)^{2} = A^{2} + A B + B A + B^{2} (not 2AB). Powers: A^{n} for square A; for a diagonal matrix, raise each diagonal entry.

Cost: the product of two n \times n matrices takes about n^{3} multiplications, which is why algorithms (Strassen) and GPUs matter.

1.2 Determinants

The determinant \det A (or |A|) is a number defined for square matrices.

  • 2 \times 2: \det \begin{bmatrix} a & b \\ c & d \end{bmatrix} = a d - b c.
  • 3 \times 3: cofactor expansion along a row or column: \det A = \sum_{j} a_{ij} C_{ij}, where C_{ij} = (- 1)^{i+j} M_{ij} and the minor M_{ij} is the determinant after deleting row i and column j. Sign pattern: [+ − +; − + −; + − +].
Worked example

\det \begin{bmatrix} 2 & 0 & 1 \\ 1 & 3 & 2 \\ 1 & 1 & 1 \end{bmatrix} expanded along row 1 = 2·(3·1 − 2·1) − 0 \cdot (\ldots) + 1·(1·1 − 3·1) = 2 (1) + (- 2) = 0. A zero determinant means the matrix is singular (no inverse).

\det \begin{bmatrix} 1 & 2 & 3 \\ 0 & 1 & 4 \\ 5 & 6 & 0 \end{bmatrix} = 1 (0 - 24) - 2 (0 - 20) + 3 (0 - 5) = - 24 + 40 - 15 = 1.

Properties (the shortcuts that make hand computation practical):

  1. \det A^{T} = \det A ;\ \det (A B) = \det A \cdot \det B ;\ \det (k A) = k^{n} \det A for n \times n.
  2. Swapping two rows changes the sign. A row of zeros, or two equal (or proportional) rows, gives \det = 0.
  3. Multiplying one row by k multiplies the determinant by k.
  4. Adding a multiple of one row to another changes nothing (use this to create zeros).
  5. For triangular matrices, det = product of the diagonal.
Worked example: row reduction for a determinant

\det \begin{bmatrix} 2 & 4 & 6 \\ 1 & 3 & 5 \\ 1 & 1 & 2 \end{bmatrix}: R_{1} \to R_{1} - 2 R_{2} gives [0 ,\ - 2 ,\ - 4] ;\ R_{3} \to R_{3} - R_{2} gives [0, −2, −3].

Now expand along column 1: its only non-zero entry is the 1 in position (2, 1), whose sign is −, so \det = - 1 \cdot \det \begin{bmatrix} - 2 & - 4 \\ - 2 & - 3 \end{bmatrix} = - 1 \cdot (6 - 8) = 2.

(Direct expansion gives the same: 2 (1) - 4 (- 3) + 6 (- 2) = 2.)

Geometry: |det A| is the factor by which the transformation x ↦ Ax scales areas (2 \times 2) or volumes (3 \times 3); a negative sign means orientation is flipped; \det = 0 means the transformation collapses space to a lower dimension.

The area of the triangle with vertices (x_{1} ,\ y_{1}) ,\ (x_{2} ,\ y_{2}) ,\ (x_{3} ,\ y_{3}) is \frac{1}{2} |det \begin{bmatrix} x_{1} & y_{1} & 1 \\ x_{2} & y_{2} & 1 \\ x_{3} & y_{3} & 1 \end{bmatrix} \mid; three points are collinear iff that determinant is 0.

1.3 Inverse of a matrix

A square A is invertible if there is A^{-1} with A A^{-1} = A^{-1} A = I. This happens iff \det A \neq 0. Properties: (A^{-1})^{-1} = A ;\ (A B)^{-1} = B^{-1} A^{-1} ;\ (A^{T})^{-1} = (A^{-1})^{T} ;\ \det A^{-1} = \frac{1}{\det A}.

By the adjoint: A^{-1} = \frac{\mathrm{adj} (A)}{\det A}, where adj(A) is the transpose of the cofactor matrix. For 2 \times 2: \begin{bmatrix} a & b \\ c & d \end{bmatrix}⁻¹ = (1/(ad − bc)) \begin{bmatrix} d & - b \\ - c & a \end{bmatrix}.

By row operations (Gauss–Jordan): write [A | I] and row-reduce until the left side is I; the right side is then A^{-1}. This is how computers do it, in about n^{3} operations, and it also tells you if A is singular.

Worked example

A = \begin{bmatrix} 2 & 1 \\ 4 & 3 \end{bmatrix}: \det = 2, so A^{-1} = \frac{1}{2} \begin{bmatrix} 3 & - 1 \\ - 4 & 2 \end{bmatrix} = \begin{bmatrix} 1.5 & - 0.5 \\ - 2 & 1 \end{bmatrix}.

Check row 1 of A times column 1 of A^{-1}: 2 (1.5) + 1 (- 2) = 1 ✓.

N = \begin{bmatrix} 1 & 2 & 3 \\ 0 & 1 & 4 \\ 5 & 6 & 0 \end{bmatrix} ,\ \det N = 1.

Cofactors give N^{-1} = \begin{bmatrix} - 24 & 18 & 5 \\ 20 & - 15 & - 4 \\ - 5 & 4 & 1 \end{bmatrix}; check the first entry of N N^{-1}: 1 (- 24) + 2 (20) + 3 (- 5) = 1 ✓.

Solving A X = B by X = A^{-1} B is elegant but in practice Gaussian elimination (Unit 2) is faster and more stable. Elementary matrices perform one row operation each; any invertible matrix is a product of them.

1.4 The determinant as an area

A 2 \times 2 determinant has a picture. Take the matrix A with columns (3, 1) and (1, 2). It sends the unit square to a parallelogram whose sides are those two columns (Figure 1.1).

The area of that parallelogram is |det A| = |3·2 − 1·1| = 5, so A stretches every area by a factor of 5. The sign tells whether the picture is flipped.

If the columns point along the same line, as in (2, 1) and (4, 2), the parallelogram collapses to a line segment, the area is 0 ,\ \det = 0, and the matrix has no inverse: a squashed picture cannot be un-squashed.

diagram
Figure 1.1. Determinant as area. Left: the matrix with columns (3,1) and (1,2) turns the unit square (light) into a parallelogram of area \det = 5. Right: columns (2,1) and (4,2) lie on one line, the parallelogram collapses, \det = 0 and the matrix is singular.

1.5 Eight uses of matrices in computing

  1. Computer graphics. Every rotation, scaling and translation of an image or game object is a matrix multiplied by a coordinate vector. A phone screen applies thousands of such products every frame. Understanding 2 \times 2 and 3 \times 3 matrices is the first step to 3D graphics.
  2. Digital images. A grey-scale photograph is a matrix of brightness numbers. Blurring, sharpening and edge detection are matrix operations on those numbers. Compression methods such as JPEG use matrix transforms to discard detail the eye cannot see.
  3. Search ranking. Google's PageRank treats the web as a huge matrix of links and finds its dominant eigenvector (Unit 4). The same idea ranks pages, researchers and social-media accounts. It is a matrix problem with billions of rows.
  4. Machine learning. Every layer of a neural network multiplies a weight matrix by an input vector. Training adjusts the matrix entries to reduce error. Fast matrix multiplication on graphics cards is why modern AI is possible.
  5. Networks and routing. An adjacency matrix records which towns, routers or villages are connected. Powers of this matrix count the number of routes of a given length. Planners use it for road, fibre and power-grid design.
  6. Economics and planning. Input-output tables show how the output of each industry feeds the others. Solving a linear system with the inverse matrix tells a planner how much each sector must produce to meet demand. The Nepal Rastra Bank and planning bodies use such tables.
  7. Cryptography. The Hill cipher encrypts blocks of letters by multiplying with a key matrix; decryption needs the inverse matrix modulo 26. It teaches why a determinant of zero makes a key useless.
  8. Statistics and data tables. A census table (districts by age group) is a matrix. Totals, averages and correlations come from matrix products. Spreadsheets and databases are everyday matrix tools.

photograph
(a) Arthur Cayley (matrix algebra, 1858)
diagram
(b) AB and BA differ; A times its inverse is I
Plate 1.1. Order matters. (a) Arthur Cayley, who in 1858 published the first paper on the algebra of matrices, including the rule for multiplying them. (b) Take A = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} and the swap matrix B = \begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix}. Then A B = \begin{bmatrix} 2 & 1 \\ 4 & 3 \end{bmatrix} (the columns of A are swapped) but B A = \begin{bmatrix} 3 & 4 \\ 1 & 2 \end{bmatrix} (the rows of A are swapped), so AB is not BA. The determinant of A is 1 x 4 - 2 x 3 = - 2, which is not zero, so A has the inverse A^{- 1} = (1 / - 2) \begin{bmatrix} 4 & - 2 \\ - 3 & 1 \end{bmatrix} = \begin{bmatrix} - 2 & 1 \\ \frac{3}{2} & - \frac{1}{2} \end{bmatrix}, and the figure checks that A times A^{- 1} is the identity. Also \det (A B) = \det (A) \det (B) = (- 2) (- 1) = 2. All entries are computed with exact fractions.Photo credits: (a) Herbert Beraud, Public domain.

photograph
(a) A dot-matrix LED display
diagram
(b) Letter A as a matrix, its transpose and a turn
Plate 1.2. A picture is a matrix. (a) A small dot-matrix LED display: every character is a grid of lit and unlit dots, which is a matrix of 1s and 0s. (b) The letter A as a 7 *x* 5 matrix with 18 lit dots, and its transpose, a 5 *x* 7 matrix whose rows are the columns of the original. Turning a picture a quarter turn clockwise is the transpose followed by reversing each row; the program checks that the transpose of the transpose is the original, that the 18 lit dots are all still there, and that the rotation matrix \begin{bmatrix} 0 & - 1 \\ 1 & 0 \end{bmatrix} has determinant +1 (a pure turn) while the mirror matrix \begin{bmatrix} - 1 & 0 \\ 0 & 1 \end{bmatrix} has determinant -1 (a reflection). This is how phones rotate photographs.Photo credits: (a) LHF, CC BY-SA 3.0.

Formula and fact sheet

  • (A B)_{ij} = \sum a_{ik} b_{kj} ;\ A B \neq B A ;\ (A B)^{T} = B^{T} A^{T} ;\ (A B)^{-1} = B^{-1} A^{-1}.
  • \det \begin{bmatrix} a & b \\ c & d \end{bmatrix} = a d - b c. \det (A B) = \det A \det B ;\ \det A^{T} = \det A ;\ \det (k A) = k^{n} \det A.
  • Row swap flips the sign; adding a row multiple changes nothing; equal rows give 0.
  • A^{-1} = \frac{\mathrm{adj} (A)}{\det A}; exists iff \det A \neq 0. 2 \times 2 inverse formula above.
  • Triangle area = \frac{1}{2} |det \begin{bmatrix} x_{1} & y_{1} & 1 \\ x_{2} & y_{2} & 1 \\ x_{3} & y_{3} & 1 \end{bmatrix} \mid.

Chapter summary

  • Matrix multiplication combines rows with columns, is associative but not commutative.
  • The determinant measures scaling and invertibility; row operations compute it efficiently.
  • A matrix is invertible exactly when its determinant is non-zero; inverses come from the adjoint or Gauss–Jordan.

Practice

  1. For A = \begin{bmatrix} 1 & 2 \\ 3 & 4 \end{bmatrix} and B = \begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix} find AB, BA, A^{2} and A^{T}.
  2. Find \det \begin{bmatrix} 3 & 1 & 2 \\ 0 & 4 & 1 \\ 2 & 0 & 5 \end{bmatrix}.
  3. Use properties to show \det \begin{bmatrix} 1 & a & b + c \\ 1 & b & c + a \\ 1 & c & a + b \end{bmatrix} = 0.
  4. Find the inverse of \begin{bmatrix} 3 & 5 \\ 1 & 2 \end{bmatrix} and verify.
  5. Find the inverse of \begin{bmatrix} 1 & 2 & 3 \\ 0 & 1 & 4 \\ 5 & 6 & 0 \end{bmatrix} by Gauss–Jordan (or check the given N^{-1}).
  6. If A is 3 \times 3 with \det A = 4, find \det (2 A) ,\ \det (A^{-1}) ,\ \det (A^{2}) and \det (A^{T}).
  7. Are the points (1, 2), (3, 8), (5, 14) collinear? Use a determinant.
  8. Find x for which \begin{bmatrix} x & 2 \\ 3 & x + 1 \end{bmatrix} is singular.
  9. If A = \begin{bmatrix} 1 & 1 \\ 0 & 1 \end{bmatrix}, find A^{n}.
  10. Show that the inverse of an orthogonal matrix is its transpose, and verify for the rotation by 90 ^{\circ} ,\ \begin{bmatrix} 0 & - 1 \\ 1 & 0 \end{bmatrix}.
Answers and worked solutions

1. A B = \begin{bmatrix} 2 & 1 \\ 4 & 3 \end{bmatrix} (columns swapped); B A = \begin{bmatrix} 3 & 4 \\ 1 & 2 \end{bmatrix} (rows swapped); A^{2} = \begin{bmatrix} 7 & 10 \\ 15 & 22 \end{bmatrix} ;\ A^{T} = \begin{bmatrix} 1 & 3 \\ 2 & 4 \end{bmatrix}.

2. 3 ( 4 \cdot 5 − 1·0) − 1 ( 0 \cdot 5 − 1·2) + 2 ( 0 \cdot 0 − 4·2) = 60 + 2 - 16 = 46.

3. Replace C_{3} by C_{3} + C_{2} (the determinant does not change): every entry of the new third column is a + b + c, so it is (a + b + c) times the first column [1 ,\ 1 ,\ 1]^{T}. Two proportional columns give det = 0.

4. \det = 6 - 5 = 1; inverse = \begin{bmatrix} 2 & - 5 \\ - 1 & 3 \end{bmatrix}.

Check \begin{bmatrix} 3 & 5 \\ 1 & 2 \end{bmatrix} \cdot \begin{bmatrix} 2 & - 5 \\ - 1 & 3 \end{bmatrix} = \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix}.

5. \begin{bmatrix} - 24 & 18 & 5 \\ 20 & - 15 & - 4 \\ - 5 & 4 & 1 \end{bmatrix} (see §1.3; multiply to confirm it gives I).

6. det 2A = 2^{3} \cdot 4 = 32; \det A^{-1} = \frac{1}{4}; \det A^{2} = 16; \det A^{T} = 4.

7. \det \begin{bmatrix} 1 & 2 & 1 \\ 3 & 8 & 1 \\ 5 & 14 & 1 \end{bmatrix} = 1 (8 - 14) - 2 (3 - 5) + 1 (42 - 40) = - 6 + 4 + 2 = 0, so collinear (they lie on y = 3 x - 1).

8. x (x + 1) - 6 = 0 \Rightarrow x^{2} + x - 6 = 0 \Rightarrow (x + 3) (x - 2) = 0 ⇒ x = - 3 or x = 2.

9. A^{2} = \begin{bmatrix} 1 & 2 \\ 0 & 1 \end{bmatrix} ,\ A^{3} = \begin{bmatrix} 1 & 3 \\ 0 & 1 \end{bmatrix}, so A^{n} = \begin{bmatrix} 1 & n \\ 0 & 1 \end{bmatrix} (induction: multiply by A).

10. A^{T} A = I means A^{T} is a left inverse of A, and for square matrices a left inverse is the inverse.

For R = \begin{bmatrix} 0 & - 1 \\ 1 & 0 \end{bmatrix} ,\ R^{T} R = \begin{bmatrix} 0 & 1 \\ - 1 & 0 \end{bmatrix} \cdot \begin{bmatrix} 0 & - 1 \\ 1 & 0 \end{bmatrix} = \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix} ✓.

Model answers: medium (5 marks) and long (10 marks) questions

Q1 (5 marks). Find the inverse of N = \begin{bmatrix} 1 & 2 & 3 \\ 0 & 1 & 4 \\ 5 & 6 & 0 \end{bmatrix} by the adjoint method.

Introduction. For a square matrix, A^{-1} = adj(A) ÷ \det (A). The inverse exists only when \det (A) \neq 0.

Step 1: determinant. Expanding along the first row: \det N = 1 ( 1 \cdot 0 − 4·6) − 2 ( 0 \cdot 0 − 4·5) + 3 ( 0 \cdot 6 − 1·5) = - 24 + 40 - 15 = 1. Since \det N = 1 \neq 0 ,\ N is invertible.

Step 2: cofactors. C11 = −24, C12 = 20, C13 = −5; C21 = 18, C22 = −15, C23 = 4; C31 = 5, C32 = −4, C33 = 1.

Step 3: adjoint and inverse. adj N is the transpose of the cofactor matrix, so N^{-1} = \begin{bmatrix} - 24 & 18 & 5 \\ 20 & - 15 & - 4 \\ - 5 & 4 & 1 \end{bmatrix} (the determinant is 1, so no division is needed).

Check. Row 1 of N times column 1 of N^{-1} is 1 (- 24) + 2 (20) + 3 (- 5) = 1, and row 1 times column 2 is 18 - 30 + 12 = 0, as required for the identity.

Conclusion N^{-1} is the matrix above.

Q2 (10 marks). Explain the operations on matrices and the properties of the determinant, with examples and the geometric meaning.

Introduction. A matrix is a rectangular array of numbers; operations on matrices model transformations, systems of equations and data tables. The determinant is a single number that records whether a square matrix can be inverted.

Operations.

  1. Addition and scalar multiplication work entry by entry and need equal sizes.
  2. Multiplication. (A B)_{ij} is row i of A times column j of B; the inner sizes must match.
  3. Non-commutative. For A = \begin{bmatrix} 2 & 1 \\ 4 & 3 \end{bmatrix} and B = \begin{bmatrix} 1 & 0 \\ 2 & 5 \end{bmatrix} ,\ A B = \begin{bmatrix} 4 & 5 \\ 10 & 15 \end{bmatrix} but B A = \begin{bmatrix} 2 & 1 \\ 24 & 17 \end{bmatrix}, so A B \neq B A in general.
  4. Transpose. Rows become columns, and (A B)^{T} = B^{T} A^{T}.
  5. Inverse. A^{-1} exists exactly when \det A \neq 0; for the 2 \times 2 case \begin{bmatrix} a & b \\ c & d \end{bmatrix}⁻¹ = (1 ÷ (ad − bc)) \begin{bmatrix} d & - b \\ - c & a \end{bmatrix}.

Determinant properties.

  1. Swapping two rows changes the sign; a row of zeros or two equal rows gives 0.
  2. Multiplying a row by k multiplies the determinant by k, and \det (A B) = \det A \cdot \det B.
  3. Geometrically |det A| is the factor by which areas (2 \times 2) or volumes (3 \times 3) are multiplied (Figure 1.1 ) ;\ \det = 0 means the transformation squashes space and cannot be undone.
Property Rule
AB versus BA usually different
(A B)^{T} B^{T} A^{T}
\det (A B) \det A \times \det B
A invertible \det A \neq 0
2 \times 2 inverse swap diagonal, negate off-diagonal, divide by det

Example. For A = \begin{bmatrix} 3 & 1 \\ 1 & 2 \end{bmatrix} ,\ \det = 5 and A^{-1} = \left( \frac{1}{5} \right) \begin{bmatrix} 2 & - 1 \\ - 1 & 3 \end{bmatrix}.

Conclusion Matrices combine data and transformations; the determinant tells at a glance whether a transformation is reversible.