3.2. Matrix Operations

3.2.1. Introduction

In Chapter 2 matrices were introduced to represent systems of linear equations. The coefficients of a linear system were put into the coefficient matrix \(A\), and a system as a whole could be squeezed into the augmented matrix. In Section 3.1 we used matrices to construct linear transformations. In this chapter we will study matrices as entities on their own, though every now and then we will keep in mind their role in the two contexts just mentioned.

3.2.2. Sum, Scalar Multiple and Transpose

In this section we will define the sum and the product of two matrices, and the transpose of a matrix. Recall that an \(m\times n\) matrix has \(m\) (horizontal) rows of size \(n\) or, equivalently, \(n\) (vertical) columns of size \(m\).

Definition 3.2.1 (Equality of matrices)

Two matrices are said to have the same size if they have the same number of rows and the same number of columns.

Two matrices \(A\) and \(B\) are equal if they have the same size, say \(m\) rows and \(n\) columns, and all the corresponding entries are equal, i.e.

\[ a_{ij} = b_{ij},\,\quad \text{for} \quad i = 1,\ldots,m, \quad j = 1,\ldots,n. \]

Definition 3.2.2

A zero matrix \(O\) is a matrix with all entries equal to 0. If the context requires clarity as to its size it may be denoted by \(O_{mn}\).

Definition 3.2.3 (Scalar multiplication)

If \(A\) is an \(m\times n\) matrix and \(c\) is a scalar, then \(cA\) is the \(m \times n\) matrix that is the result of multiplying each entry of \(A\) by \(c\):

\[\begin{split} c \left[\begin{array}{cccc} a_{11} & a_{12}& \ldots& a_{1n} \\ a_{21} & a_{22}& \ldots& a_{2n} \\ \vdots & \vdots& \cdots& \vdots \\ a_{m1} & a_{m2}& \ldots& a_{mn} \end{array} \right]= \left[\begin{array}{cccc} ca_{11} & ca_{12}& \ldots& ca_{1n} \\ ca_{21} & ca_{22}& \ldots& ca_{2n} \\ \vdots & \vdots& \cdots& \vdots \\ ca_{m1} & ca_{m2}& \ldots& ca_{mn} \end{array} \right]. \end{split}\]

We then say that \(cA\) is a scalar multiple of \(A\), or simply a multiple of \(A\).

Definition 3.2.4 (The sum of two matrices)

If \(A\) and \(B\) are two \(m\times n\) matrices then the sum \(A+B\) is the \(m\times n\) matrix of which the entry on the position \((i,j)\) is the sum of the corresponding entries of \(A\) and \(B\):

\[\begin{split} \left[\begin{array}{cccc} a_{11} & a_{12}& \ldots& a_{1n} \\ a_{21} & a_{22}& \ldots& a_{2n} \\ \vdots & \vdots& \cdots& \vdots \\ a_{m1} & a_{m2}& \ldots& a_{mn} \end{array} \right] + \left[\begin{array}{cccc} b_{11} & b_{12}& \ldots& b_{1n} \\ b_{21} & b_{22}& \ldots& b_{2n} \\ \vdots & \vdots& \cdots& \vdots \\ b_{m1} & b_{m2}& \ldots& b_{mn} \end{array} \right]= \end{split}\]

\[\begin{split} = \left[\begin{array}{cccc} a_{11}+b_{11} & a_{12}+b_{12}& \ldots& a_{1n}+b_{1n} \\ a_{21}+b_{21} & a_{22}+b_{22}& \ldots& a_{2n}+b_{2n} \\ \vdots & \vdots& \cdots& \vdots \\ a_{m1}+b_{m1} & a_{m2}+b_{m2}& \ldots& a_{mn}+b_{mn} \end{array} \right]. \end{split}\]

If \(A\) and \(B\) are not of the same size their sum is not defined.

Example 3.2.1

\[\begin{split} \begin{bmatrix} 1 & 3 \\ 5 & 2 \\ 6 & -4 \end{bmatrix} + \begin{bmatrix} 3 & 2 \\ 4 & -5 \\ 2 & 5 \end{bmatrix} = \begin{bmatrix} 4 & 5 \\ 9 & -3 \\ 8 & 1 \end{bmatrix}, \end{split}\]

\[\begin{split} \begin{bmatrix} 1 & 3 \\ 5 & 2 \\ 6 & -4 \end{bmatrix} + \begin{bmatrix} 0 & 0 \\ 0 & 0 \\ 0 & 0 \end{bmatrix} = \begin{bmatrix} 0 & 0 \\ 0 & 0 \\ 0 & 0 \end{bmatrix} + \begin{bmatrix} 1 & 3 \\ 5 & 2 \\ 6 & -4 \end{bmatrix} = \begin{bmatrix} 1 & 3 \\ 5 & 2 \\ 6 & -4 \end{bmatrix}, \end{split}\]

\[\begin{split} \begin{array}{lcl} \begin{bmatrix} 1 & 3 & 5 \\ 2 & 4 & 1 \end{bmatrix} + (-1)\begin{bmatrix} 1 & 3 & 5 \\ 2 & 4 & 1 \end{bmatrix} &=& \begin{bmatrix} 1 & 3 & 5 \\ 2 & 4 & 1 \end{bmatrix} + \begin{bmatrix} -1 & -3 &-5 \\ -2 & -4 & -1 \end{bmatrix} \\ &=& \begin{bmatrix} 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix}. \end{array} \end{split}\]

The multiple \((-1)A\) is also written as \(-A\). An obvious property, illustrated in the third example, is:

\[ A + (-A) = O, \]

where \(O\) is the zero matrix.

Example 3.2.2

\[\begin{split} \begin{bmatrix} 1 & 3 \\ 5 & 2 \\ 6 & -4 \end{bmatrix} + \begin{bmatrix} 1 & 3 & 5 \\ 2 & 4 & 1 \end{bmatrix} \end{split}\]

is not defined. This is because the matrices do not have the same size.

Remark 3.2.1

The two definitions of sum and scalar multiple are called componentwise definitions. They are completely analogous to the definitions of the scalar multiples of a vector and the sum of two vectors. Hence it is not surprising that they obey exactly the same rules, as is summarized in the next proposition. (cf. Section Vectors.)

Proposition 3.2.1

Suppose \(A, B\) and \(C\) are \(m\times n\) matrices and let \(c_{1},c_{2}\) be two real numbers. Then we have:

• \(A+O_{mn}=A=O_{mn}+A\)

• \((A+B)+C=A+(B+C)\)

• \(A+B=B+A\)

• \(A+(-A)=O\)

• \(1A=A\)

• \(c_{1}(A+B)=c_{1}A+c_{1}B\)

• \((c_{1}+c_{2})A=c_{1}A+c_{2}A\)

• \(c_{1}(c_{2}A)=(c_{1}c_{2})A\).

An operator of which the usefulness is not immediately clear, but which fits well in this section with matrix operations, is the following:

Definition 3.2.5

The transpose of an \(m \times n\) matrix \(A\) with entries \(a_{ij}\) is the \(n \times m\) matrix \(B\) with entries \(b_{ij}\) defined by

\( b_{ij} = a_{ji}\).   It is denoted by \(B = A^T\).

Example 3.2.3

\[\begin{split} \begin{bmatrix} 1 & 3 \\ 5 & 2 \\ 6 & 4 \end{bmatrix}^T = \begin{bmatrix} 1 & 5 & 6 \\ 3 & 2 & 4 \end{bmatrix} \quad \text{and} \quad \begin{bmatrix} -1 & 2 & -4 & 0\end{bmatrix}^T = \begin{bmatrix} -1 \\ 2 \\-4 \\ 0\end{bmatrix}. \end{split}\]

The following rules involving the three operators defined so far in this section are easy to prove:

Proposition 3.2.2

Let \(A\) and \(B\) be \(m\times n\) matrices and \(c\) a scalar. Then we have

• \((cA)^T = c A^T\)

• \((A+B)^T = A^T + B^T\)

• \((A^T)^T = A\).

Proof. We will prove the second statement and leave the other two to the diligent reader. See Exercise 3.2.1.

So, suppose \(A\) and \(B\) are two \(m \times n\) matrices. Then \(A+B\) is an \(m \times n\) matrix too, hence \((A+B)^T\) is an \(n \times m\) matrix. The matrix \(A^T + B^T\) on the right-hand side of the equation is the sum of two \(n \times m\) matrices, which is again an \(n \times m\) matrix. So the matrices on both sides of the equation have the same size.

Next we have to show that they have equal entries on the corresponding positions. If we put

\[ E = (A+B)^T \quad \text{and}\quad F = A^T + B^T \]

we see that

\[ e_{ij} = \text{ entry of } (A+B) \,\text{ on position }(j,i) \]

and

\[\begin{split} \begin{array}{rl} f_{ij} &= \text{ entry of } A^T \,\text{ on position }(i,j)\,+ \text{ entry of } B^T \,\text{ on position }(i,j) \\ &= \text{ entry of } A \,\text{ on position }(j,i)\,+ \text{ entry of } B \,\text{ on position }(j,i)\\ &= \text{ entry of } (A+B) \,\text{ on position }(j,i)\\ &= \,\,\,\,e_{ij}, \end{array} \end{split}\]

so we are done.

If you are lost in the forest of indices, have a look at Example 3.2.4.

Example 3.2.4

We check property (ii) for two general \(3\times 4\) matrices \(A\) and \(B\) on the position \((2,3)\). Let

\[\begin{split} A = \begin{bmatrix} a_{11}& a_{12} & a_{13} & a_{14} \\ a_{21}& a_{22} & a_{23} & a_{24} \\ a_{31} & \fbox{$a_{32}$} & a_{33} & a_{34} \end{bmatrix} \quad \text{and} \quad B = \begin{bmatrix} b_{11}& b_{12} & b_{13} & b_{14} \\ b_{21}& b_{22} & b_{23} & b_{24} \\ b_{31} & \fbox{$b_{32}$} & b_{33} & b_{34} \end{bmatrix}. \end{split}\]

Then

\[\begin{split} E = (A+B)^T = \begin{bmatrix} a_{11}+b_{11}& a_{12}+b_{12} & a_{13}+b_{13} & a_{14} +b_{14}\\ a_{21}+b_{21}& a_{22}+b_{22} & a_{23}+b_{23} & a_{24}+b_{24} \\ a_{31}+b_{31} & \fbox{$a_{32}+b_{32}$} & a_{33}+b_{33} & a_{34}+b_{34} \end{bmatrix}^T \end{split}\]

so

\[\begin{split} E = \begin{bmatrix} a_{11}+b_{11}& a_{21}+b_{21} & a_{31}+b_{31} \\ a_{12}+b_{12}& a_{22}+b_{22} & \fbox{$a_{32}+b_{32}$} \\ a_{13}+b_{13}& a_{23}+b_{23} & a_{33}+b_{33} \\ a_{14} +b_{14}& a_{24}+b_{24} & a_{34}+b_{34} \end{bmatrix}, \end{split}\]

and on position \((2,3)\) we have \(a_{32}+b_{32}\).

On the other hand

\[\begin{split} F = A^T + B^T = \begin{bmatrix} a_{11}& a_{21} & a_{31} \\ a_{12}& a_{22} & \fbox{$a_{32}$} \\ a_{13}& a_{23} & a_{33}\\ a_{14}& a_{24} & a_{34}\end{bmatrix} + \begin{bmatrix} b_{11}& b_{21} & b_{31} \\ b_{12}& b_{22} & \fbox{$b_{32}$} \\ b_{13}& b_{23} & b_{33}\\ b_{14}& b_{24} & b_{34} \end{bmatrix}, \end{split}\]

with on position \((2,3)\) the value \(a_{32}\) + \(b_{32}\).

Exercise 3.2.1

Prove statements (i) and (iii) of Proposition 3.2.2.

Solution to Exercise 3.2.1 (click to show)

Suppose \(A = \left[\begin{array}{cccc} a_{11} & a_{12}& \ldots& a_{1n} \\ a_{21} & a_{22}& \ldots& a_{2n} \\ \vdots & \vdots& \cdots& \vdots \\ a_{m1} & a_{m2}& \ldots& a_{mn} \end{array} \right]\) is an arbitrary \(m \times n\) matrix. Then

\[\begin{split} \begin{array}{rcl} (cA)^T &=& \left[\begin{array}{cccc} ca_{11} & ca_{12}& \ldots& ca_{1n} \\ ca_{21} & ca_{22}& \ldots& ca_{2n} \\ \vdots & \vdots& \cdots& \vdots \\ ca_{m1} & ca_{m2}& \ldots& ca_{mn} \end{array}^T \right] = \left[\begin{array}{cccc} ca_{11} & ca_{21}& \ldots& ca_{m1} \\ ca_{12} & ca_{22}& \ldots& ca_{m2} \\ \vdots & \vdots& \cdots& \vdots \\ ca_{1n} & ca_{2n}& \ldots& ca_{mn} \end{array} \right] \\ &=& c \left[\begin{array}{cccc} a_{11} & a_{21}& \ldots& a_{m1} \\ a_{12} & a_{22}& \ldots& a_{m2} \\ \vdots & \vdots& \cdots& \vdots \\ a_{1n} & a_{2n}& \ldots& a_{mn} \end{array} \right] = c\,A^T \end{array} \end{split}\]

Example 3.2.5

We will solve the equation \(A + 2X^T + B = C\) for \(X\), where

\[\begin{split} A = \begin{bmatrix} 1 & 1 & 2 \\ 3 & 1 & 0 \end{bmatrix}, \quad B = \begin{bmatrix} 2 & 0 & 3 \\ 2 & 3 & 4 \end{bmatrix}, \text{ and} \quad C = \begin{bmatrix} 7 & 5 & 1 \\ 1 & 4 & 2 \end{bmatrix}. \end{split}\]

We will extricate \(X\) step by step:

\[ A + 2X^T + B = C \,\, \iff \,\,2X^T = C-A-B \,\, \iff \,\, X^T = \tfrac12(C-A-B). \]

Next we transpose both terms to find

\[\begin{split} X = \tfrac12(C-A-B)^T = \frac12\begin{bmatrix} 4 & 4 & -4 \\ -4 & 0 & -2 \end{bmatrix}^T = \begin{bmatrix} 2 & -2 \\ 2 & 0 \\ -2 & -1 \end{bmatrix} \end{split}\]

3.2.3. Grasple Exercises

Grasple Exercise 3.2.1

To compute the sum of two matrices

Show/Hide Content

Grasple Exercise 3.2.2

To compute \(c_1A + c_2B\).

Show/Hide Content

Grasple Exercise 3.2.3

To compute \(c_1A + c_2B\).

Show/Hide Content

Grasple Exercise 3.2.4

To solve equations involving sum and transpose

Show/Hide Content

Grasple Exercise 3.2.5

True/False questions involving sum and transpose

Show/Hide Content

3.2.4. The product of two matrices

Next we turn our attention to the most important matrix operation, namely the product \(AB\) of two matrices. In the previous chapter we have already seen the special case where \(B\) is a matrix of just one column, i.e.,

\[\begin{split} B = \mathbf{x} = \begin{bmatrix}x_1 \\ x_2 \\ \vdots \\ x_n \end{bmatrix}, \end{split}\]

a vector in \(\mathbb{R}^n\), which we can identify with an \(n \times 1\) matrix. We want of course the definition of the general matrix product to be consistent with this.

Definition 3.2.6

The product of an \(m\times n\) matrix \(A\) and an \(n\times p\) matrix \(B = [\,{\vect{b}_1}\quad {\vect{b}_2}\quad \ldots \quad {\mathbf{b}_p}]\) is defined by

\[ AB = [\,A\mathbf{b}_1\quad A\mathbf{b}_2\quad \ldots \quad A\mathbf{b}_p]. \]

So we have

\[ j\text{-th column of } AB = A\text{ times $j$-th column of } B, \quad j = 1,2,\ldots,p. \]

Note that this makes \(AB\) an \(m \times p\) matrix.

If the number of columns of \(A\) is not equal to the number of rows of \(B\) the product \(AB\) is not defined.

Example 3.2.6

\[\begin{split} \begin{bmatrix} 1 & -3 \\ -1 & 2 \\ 3& -2 \end{bmatrix} \begin{bmatrix} 2 & 1 & 1\\ 3 & 0 & 2 \end{bmatrix} = \begin{bmatrix} -7 & 1 & -5 \\ 4 & -1 & 3 \\ 0 & 3 &-1 \end{bmatrix}. \end{split}\]

For instance, the third column is computed as

\[\begin{split} \begin{bmatrix} 1 & -3 \\ -1 & 2 \\ 3& -2 \end{bmatrix} \begin{bmatrix} 1\\ 2 \end{bmatrix} = 1\begin{bmatrix} 1 \\ -1 \\ 3\end{bmatrix} + 2\begin{bmatrix} -3 \\ 2 \\ -2 \end{bmatrix} \,\,=\, \begin{bmatrix} -5 \\ 3 \\ -1\end{bmatrix}. \end{split}\]

Proposition 3.2.3

The product of the \(m\times n\) matrix \(A\) and the \(n\times p\) matrix \(B\) is the \(m\times p\) matrix \(C\) for which the entry on the position \((i,j)\) is given by

\[\begin{split} c_{ij} = a_{i1}b_{1j} + a_{i2}b_{2j} + \ldots + a_{in}b_{nj} = \begin{bmatrix}a_{i1} & a_{i2} & \ldots & a_{in} \end{bmatrix} \begin{bmatrix} b_{1j} \\ b_{2j} \\ \vdots \\ b_{nj}\end{bmatrix}. \end{split}\]

This is sometimes called the row-column expansion of the product.

Proof. We already saw this row-column expansion in Section 2.4.

The following scheme nicely visualizes the row-column expansion

\[\begin{split} \begin{array}{ccc} & \begin{bmatrix} b_{11} & b_{12}& \ldots& {\color{red}b_{1j}} & \ldots& b_{1p} \\ b_{21} & b_{22}& \ldots& {\color{red}b_{2j}} & \ldots& b_{2p} \\ \vdots & \vdots& \ldots& & \ldots& \vdots \\ b_{n1} & b_{n2}& \ldots& {\color{red}b_{nj}} & \ldots& b_{np} \end{bmatrix} \\ \begin{bmatrix} a_{11} & a_{12}& \ldots& \ldots& a_{1n} \\ a_{21} & a_{22}& \ldots& \ldots& a_{2n} \\ \cdots & \ddots& \ldots& \ldots& \vdots \\ {\color{red}a_{i1}} & {\color{red}a_{i2}}& {\color{red}\cdots}& \ldots& {\color{red}a_{in}} \\ \vdots & \vdots& \ldots& \ldots& \vdots \\ a_{m1} & a_{m2}& \ldots& \ldots& a_{mn} \end{bmatrix} \!\! & \! \begin{bmatrix} c_{11} & c_{12}& \ldots& c_{1j} &\ldots& c_{1p} \\ c_{21} & c_{22}& \ldots& c_{2j} &\ldots& c_{2p} \\ \vdots & \vdots& \ldots& & \ldots& \vdots \\ c_{i1} & c_{i2}& \cdots&{\color{red}c_{ij}} &\ldots& c_{ip} \\ \vdots & \vdots& \ldots& &\ldots& \vdots \\ c_{m1} & c_{m2}& \ldots& c_{mn} &\ldots& c_{mp} \end{bmatrix} \end{array}\end{split}\]

Example 3.2.7

Let us consider the same matrix product

\[\begin{split} \begin{bmatrix} 1 & -3 \\ -1 & 2 \\ 3& -2 \end{bmatrix} \begin{bmatrix} 2 & 1 & 1\\ 3 & 0 & 2 \end{bmatrix} = \begin{bmatrix} -7 & 1 & -5 \\ 4 & -1 & 3 \\ 0 & 3 &-1 \end{bmatrix}. \end{split}\]

The \(-5\) on position \((1,3)\) and the \(3\) on position \((3,2)\) in the product come from

\[\begin{split} -5 = \begin{bmatrix} 1 & -3 \end{bmatrix} \begin{bmatrix} 1\\ 2 \end{bmatrix} \quad \text{and} \quad 3 = \begin{bmatrix} 3 & -2 \end{bmatrix} \begin{bmatrix} 1\\ 0 \end{bmatrix}. \end{split}\]

Exercise 3.2.2

Explain why the product

\[\begin{split} \begin{bmatrix} 1 & -3 \\ -1 & 2 \\ 3& -2 \end{bmatrix}\begin{bmatrix} 1 & -3 \\ -1 & 2 \\ 3& -2 \end{bmatrix} \end{split}\]

is not defined.

Remark 3.2.2

The product of a matrix \(A\) with itself is only defined if \(A\) is an \(n \times n\) matrix. In that case we use the obvious notation

\[ A^2 = A\cdot A. \]

Example 3.2.8

\[\begin{split} \begin{bmatrix} a_{11}& a_{12} & a_{13} \\ a_{21}& a_{22} & a_{23} \\ a_{31}& a_{32} & a_{33} \\ a_{41} & a_{42} & a_{43} \end{bmatrix} \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0\\ 0 & 0 & 1 \end{bmatrix} = \begin{bmatrix} a_{11}& a_{12} & a_{13} \\ a_{21}& a_{22} & a_{23} \\ a_{31}& a_{32} & a_{33} \\ a_{41} & a_{42} & a_{43} \end{bmatrix}. \end{split}\]

This example illustrates the existence of a unit element with respect to the multiplication. To identify it we first introduce some more terminology.

Definition 3.2.7

An \(n\times n \) matrix \(A\) is called a square matrix. So it is a matrix where the number of columns is equal to the number of rows.

For a square matrix \(A\) we call the elements \(a_{ii}\) the diagonal elements. Together the diagonal elements form the (main) diagonal of \(A\).

A square matrix where all non-diagonal elements are equal to 0 is called a diagonal matrix.

Remark 3.2.3

The other diagonal of a square matrix, the one from bottom left to top right, plays a minor role. For this reason we don’t reserve a name for it. By ‘diagonal’ we will always mean: main diagonal.

Example 3.2.9

Consider the matrices

\[\begin{split} A = \begin{bmatrix} 2 & 2 \\ 3 & 3 \end{bmatrix}, \quad B = \begin{bmatrix} 2 & 0 & 0 \\ 0 & 3 & 0 \\ 0 & 0 & 6 \end{bmatrix}, \quad C = \begin{bmatrix} 1 & 0 \\ 0 & 1 \\ 0 & 0 \end{bmatrix}. \end{split}\]

The matrices \(A\) and \(B\) are square, and only \(B\) is a diagonal matrix.

Exercise 3.2.3

Is the following statement true or false?

The \(n \times n\) zero matrix \(O_{nn}\) is a diagonal matrix.

Exercise 3.2.4

Suppose \(A = \begin{bmatrix} \mathbf{a}_1 & \mathbf{a}_2 & \ldots & \mathbf{a}_n \end{bmatrix} \) is an \(m\times n\) matrix and
\(B= \begin{bmatrix} \mathbf{b}_1 & \mathbf{b}_2 & \ldots & \mathbf{b}_p \end{bmatrix} \) an \(n\times p\) matrix. Show that

\[\begin{split} A^TB = %\begin{bmatrix} \vect{a_1}^T \vect{b_1} & \vect{a_1}^T \vect{b_2} & \ldots & \vect{a_1}^T \vect{b_p} \\ % \vect{a_2}^T \vect{b_1} & \vect{a_2}^T \vect{b_2} & \ldots & \vect{a_2}^T \vect{b_p} \\ % \vdots & \vdots & & \vdots \\ % \vect{a_n}^T \vect{b_1} & \vect{a_n}^T \vect{b_2} & \ldots & \vect{a_n}^T \vect{b_p} \\ %\end{bmatrix} % = \begin{bmatrix} \mathbf{a}_1\ip \mathbf{b}_1 & \mathbf{a}_1\ip\mathbf{b}_2 & \ldots & \mathbf{a}_1\ip \mathbf{b}_p \\ \mathbf{a}_2\ip \mathbf{b}_1 & \mathbf{a}_2\ip\mathbf{b}_2 & \ldots & \mathbf{a}_2\ip \mathbf{b}_p \\ \vdots & \vdots & & \vdots \\ \mathbf{a}_n\ip \mathbf{b}_1 & \mathbf{a}_n\ip\mathbf{b}_2 & \ldots & \mathbf{a}_n\ip \mathbf{b}_p \\ \\ \end{bmatrix},\end{split}\]

where \(\mathbf{a}\ip\mathbf{b}\) is the dot product of the vectors \(\mathbf{a}\) and \(\mathbf{b}\).

Exercise 3.2.5

The special case in the previous exercise where \(A = B\) will become very important when we will look at orthogonal projections. For now, show that: the columns of a matrix \(A\) are orthogonal if and only if the matrix \(A^TA\) is a diagonal matrix.

Definition 3.2.8

The identity matrix \(I_n\) is the \(n \times n\) diagonal matrix with 1’s on the diagonal. If the size is irrelevant or clear from the context, we denote it simply by \(I\).

Exercise 3.2.6

Let

\[\begin{split} I = I_4 = \begin{bmatrix}1 & 0 & 0 & 0 \\ 0 & 1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 1 \end{bmatrix} \quad \text{and} \quad A = \begin{bmatrix} a_{11}& a_{12} & a_{13} \\ a_{21}& a_{22} & a_{23} \\ a_{31}& a_{32} & a_{33} \\ a_{41} & a_{42} & a_{43} \end{bmatrix}. \end{split}\]

Show that \(IA = A\).

The definition of the product of two matrices and the earlier definition of the product of a matrix and a vector (Definition 2.4.1) immediately imply that the columns of the product of two matrices are linear combinations of the columns of the first matrix.
As is often the case in linear algebra things can be looked at from a different perspective. From Proposition 3.2.3 it follows that the elements \(c_{i1},c_{i2},\ldots,c_{in}\) of the \(i\)-th row of the product \(C = AB\) as far as \(A\) is concerned only depend on the elements \(a_{ik}\) of its \(i\)-th row. The following proposition explains in which way.

Proposition 3.2.4

The \(i\)-th row of the product \(AB\) is the linear combination of the rows of the second matrix, \(B\), with the entries of the \(i\)-th row of \(A\) as coefficients.

Proof. The indicated linear combination yields:

\[ a_{i1} \begin{bmatrix}b_{11} & b_{12} & \ldots &b_{1p} \end{bmatrix} + a_{i2} \begin{bmatrix}b_{21} & b_{22} & \,\, \ldots \,\, &b_{2p} \end{bmatrix} + \ldots + a_{in} \begin{bmatrix}b_{n1} & b_{n2} & \ldots &b_{np} \end{bmatrix} \]

\[ = \begin{bmatrix} (a_{i1}b_{11} + a_{i2}b_{21}+ \ldots +a_{in}b_{n1}) & \quad\ldots\quad & (a_{i1}b_{1p} + a_{i2} b_{2p} + \ldots + a_{in}b_{np}) \end{bmatrix}. \]

This is a row vector with on the \(j\)-th position the number

\[ (a_{i1}b_{1j} + a_{i2} b_{2j} + \ldots + a_{in}b_{nj}), \]

and that is precisely the entry \(c_{ij}\) of the matrix \(C = AB\).

Interestingly this opens the way to describe the row operations of Chapter 2 via matrix multiplication. The following example illustrates this for the three basic row operations.

Example 3.2.10

The following multiplication adds the first row of the matrix

\[\begin{split} A = \begin{bmatrix} a_{11}& a_{12} & a_{13} \\ a_{21}& a_{22} & a_{23} \\ a_{31}& a_{32} & a_{33} \end{bmatrix} \end{split}\]

four times to the second row:

\[\begin{split} \begin{bmatrix} 1 & 0 & 0 \\ 4 & 1 & 0 \\ 0 & 0 & 1\end{bmatrix} \begin{bmatrix} a_{11}& a_{12} & a_{13} \\ a_{21}& a_{22} & a_{23} \\ a_{31}& a_{32} & a_{33} \\ \end{bmatrix} = \begin{bmatrix} a_{11}& a_{12} & a_{13} \\ 4a_{11}+a_{21}&4a_{12} +a_{22}& 4a_{13}+a_{23} \\ a_{31}& a_{32} & a_{33} \\ \end{bmatrix}. \end{split}\]

Here the third row is scaled with a factor 5:

\[\begin{split} \begin{bmatrix} 1 & 0 & 0 \\ 0 & 1 & 0 \\ 0 & 0 & 5\end{bmatrix} \begin{bmatrix} a_{11}& a_{12} & a_{13} \\ a_{21}& a_{22} & a_{23} \\ a_{31} & a_{32} & a_{33} \end{bmatrix} = \begin{bmatrix} a_{11}& a_{12} & a_{13} \\ a_{21}& a_{22} & a_{23} \\ 5a_{31}& 5a_{32} & 5a_{33} \\ \end{bmatrix}. \end{split}\]

And with the following multiplication the first and third row of \(A\) are swapped:

\[\begin{split} \begin{bmatrix} 0 & 0 & 1 \\ 0 & 1 & 0 \\ 1 & 0 & 0\end{bmatrix} \begin{bmatrix} a_{11}& a_{12} & a_{13} \\ a_{21}& a_{22} & a_{23} \\ a_{31}& a_{32} & a_{33} \end{bmatrix} = \begin{bmatrix}a_{31}& a_{32} & a_{33} \\ a_{21}& a_{22} & a_{23} \\ a_{11}& a_{12} & a_{13} \end{bmatrix}. \end{split}\]

For future reference we give these matrices a name:

Definition 3.2.9

The matrices \(E\) that perform one single row operation (row replacement, row scaling, or row exchange) via \(A \mapsto EA\) are called elementary matrices.

Exercise 3.2.7

Describe in words which row operations are the effect of pre-multiplying a \(4\times n\) matrix \(A\) with the following elementary matrices:

\[\begin{split} E_1 = \begin{bmatrix} 1 & 0 & 0 & 0\\ 0 & 1 & 0 & -1\\ 0 & 0 & 1 & 0\\ 0 & 0 & 0 & 1 \end{bmatrix}, \quad \quad E_2 = \begin{bmatrix} 1 & 0 & 0 & 0\\ 0 & 0 & 0 & 1\\ 0 & 0 & 1 & 0\\ 0 & 1 & 0 & 0 \end{bmatrix}. \end{split}\]

Example 3.2.11

The following product may at first sight seem a bit odd, but it is exactly according to the definition:

\[\begin{split} \begin{bmatrix} 1 \\-2\\3\\4 \end{bmatrix}\begin{bmatrix} 2&4&0& -1 \end{bmatrix} = \begin{bmatrix} 2 & 4 & 0 & -1 \\ -4 & -8 & 0 & 2 \\ 6 & 12 & 0 & -3 \\ 8 & 16 & 0 & -4 \end{bmatrix}. \end{split}\]

The column-row product in the last example is the building block for yet another way to look at the matrix product. The next exercise explains how.

Exercise 3.2.8

Column-row expansion of the product.

Denote the columns of the \(m\times n\) matrix \(A\) by \(A_{(1)}, \ldots, A_{(n)}\), and the rows of the \(n\times p\) matrix \(B\) by \(B^{(1)}, \ldots, B^{(p)}\), so

\[\begin{split} A_{(j)} = \begin{bmatrix} a_{1j} \\ \vdots \\ a_{mj}\end{bmatrix} \quad \text{and} \quad B^{(i)} = \begin{bmatrix} b_{i1} & b_{i2} & \ldots & b_{ip}\end{bmatrix}. \end{split}\]

Show that

\[ AB = A_{(1)} B^{(1)} + A_{(2)} B^{(2)} + \ldots + A_{(n)} B^{(n)}, \]

i.e., \(AB\) is the sum of \(n\) column-row products (like in Example 3.2.11).

3.2.5. Properties of the matrix product

Now let us have a look which of the rules of the products of numbers also hold for products of matrices. And which do not.

Proposition 3.2.5

For all \(m \times n\) matrices \(A,A_1,A_2\), all \(n \times p\) matrices \(B,B_1,B_2\), all \(p \times q\) matrices \(C\) and all real numbers \(c\) the following are true:

1. \(A(B_1+B_2) = AB_1 + AB_2\) and \((A_1+A_2)B = A_1B+A_2B\);

2. \(A(cB) = c(AB) = (cA)B\);

3. \(AI_n = A\) and \(I_mA = A\) (the identity matrix \(I\) acts as a unit element);

4. \(A(BC) = (AB)C\).

Example 3.2.12

As an illustration of rule iv. we compute the two triple products for the three matrices

\[\begin{split} A = \begin{bmatrix} 3 & 1 \\ 2 & 1 \\ 0 & 5 \end{bmatrix}, \quad B = \begin{bmatrix} 1 & 2 \\ 3 & 0 \end{bmatrix}, \quad C = \begin{bmatrix} 1 & 2 & 0 \\ 2 & 1 & 2 \end{bmatrix}. \end{split}\]

On the one hand

\[\begin{split} A(BC) = \begin{bmatrix} 3 & 1 \\ 2 & 1 \\ 0 & 5 \end{bmatrix} \begin{bmatrix} 5 & 4 & 4\\ 3 & 6 & 0 \end{bmatrix} = \begin{bmatrix} 18 & 18 & 12\\ 13 & 14 & 8 \\ 15 & 30 & 0 \end{bmatrix}, \end{split}\]

and on the other hand

\[\begin{split} (AB)C = \begin{bmatrix} 6 & 6 \\ 5 & 4 \\ 15 & 0 \end{bmatrix} \begin{bmatrix} 1 & 2 & 0 \\ 2 & 1 & 2 \end{bmatrix} = \begin{bmatrix} 18 & 18 & 12 \\ 13 & 14 & 8 \\ 15 & 30 & 0 \end{bmatrix}. \end{split}\]

So the products are indeed equal. But it is not immediately clear how. For instance, the value 14 on position (2,2) comes about in two ways

\[ \text{via } A(BC)\!: \,14 = 2\cdot4 + 1\cdot 6, \quad \,\, \text{via } (AB)C\!: \,14 = 5\cdot2 + 4\cdot1. \]

We need a good perspective to give a proof of the general case.

Proof. (of Proposition 3.2.5)
Rules i. and ii. are checked in a straightforward way. See Exercise 3.2.9.

3. We saw instances of this property already in Example 3.2.8 and Exercise 3.2.6. For the general case, one way to show validity of the first statement is to note that the \(j\)-th column of \(AI_n\) is \(A\mathbf{e}_j\) where \(\mathbf{e}_j\) is the \(j\)-th column of the identity matrix \(I_n\). This gives the linear combination

   
   
   \[ A\mathbf{e}_j = 0\mathbf{a}_1 + 0\mathbf{a}_2 + \ldots + 1\mathbf{a}_j +\dots + 0\mathbf{a}_n = \mathbf{a}_j \]

   which shows that the \(j\)-th column of \(AI_n\) is equal to the \(j\)-th column of \(A\). And this holds for any column.

   The identity \(\quad I_mA = A\quad\) is shown in an analogous way, working row by row.

4. First we observe that both triple products yield \(m \times q\) matrices. Then the identity can be proved ‘column by column’, as the previous one.

   We are done if we can show that

   
   
   \[\begin{split} \begin{array}{rcl} k\text{-th column of }A(BC) &=& k\text{-th column of }(AB)C \\ &=& (AB)( k\text{-th column of }C) = (AB)\mathbf{c}_k, \end{array} \end{split}\]

   for \( k = 1,2,\ldots q \).

   Now recall that (by definition)

   
   
   \[ k\text{-th column of }BC = B\vect{c}_k, \]

   so

   \[ k\text{-th column of }A(BC) = A\,(B\vect{c}_k) \]

   Making extensive use of the rule

   
   
   \[ A(c_1\mathbf{x} + c_2\mathbf{y}) = c_1A\mathbf{x} + c_2A\mathbf{y} \]

   we find

   \[\begin{split} \begin{array}{ccl} A\,(B\mathbf{c_k}) & = & A \,(c_{1k}\mathbf{b}_1 +c_{2k}\mathbf{b}_2 + \ldots + c_{pk}\mathbf{b}_p)\\ & = & c_{1k}(A\mathbf{b}_1) +c_{2k}(A\mathbf{b}_2) + \ldots + c_{pk}(A\mathbf{b}_p)\\ & = & \begin{bmatrix} A\mathbf{b}_1 & A\mathbf{b}_2 & \ldots & A\mathbf{b}_p \end{bmatrix} \begin{bmatrix} c_{1k} \\ \vdots \\ c_{pk} \end{bmatrix} \\ & = & (AB)\mathbf{c}_k. \end{array} \end{split}\]

Exercise 3.2.9

Prove rules i. and ii. of Proposition 3.2.5.

Recall that matrices are equal when they have the same size and the entries on corresponding positions are equal (which may be checked column by column or row by row).

Remark 3.2.4

The proof of Proposition 3.2.5 iv. can be seen in another light. In the Section Linear Transformations we saw that an \(m\times n\) matrix \(A\) defines a transformation \(T\) from \(\mathbb{R}^n\) to \(\mathbb{R}^m\), namely

\[ \text{for } \mathbf{x} \in \mathbb{R}^n: \quad \mathbf{x} \mapsto T(\mathbf{x}) = A\mathbf{x}. \]

The definition of the product of two matrices then precisely matches the composition of two of such transformations: if \(A\) is an \(m\times n\) matrix and \(B\) is an \(n\times p\) matrix

\[ \mathbf{x}\in\mathbb{R}^p \,\,\stackrel{B}{\longrightarrow}\,\, \vect{y}_1 = B\vect{x}\in\mathbb{R}^n \,\, \stackrel{A}{\longrightarrow} \,\,\, \vect{y}_2 = A(B\mathbf{x}) \in \mathbb{R}^m\, \]

and

\[ \mathbf{x}\in\mathbb{R}^p\,\,\,\stackrel{AB}{\longrightarrow} \,\,\,\vect{y}_3 \,=\, (AB)\mathbf{x} \in \mathbb{R}^m \]

yield the same vector:

\[ \vect{y}_2 = \vect{y}_3. \]

So far so good: matrix multiplication behaves as multiplication of numbers. However, in two important respects the concepts deviate. First of all, commutativity no longer holds.

Example 3.2.13

For the matrices

\[\begin{split} A = \begin{bmatrix} 2 & 2 & 1\\ 3 & 3 & 0 \end{bmatrix} \quad \text{and} \quad B = \begin{bmatrix} 1 & 3 \\ 3 & 1 \\ 4 & 0 \end{bmatrix} \end{split}\]

it is clear that

\[ AB \neq BA,\]

simply because the two products are not of the same size: \(AB\) is a \(2\times 2\) matrix, \(BA\) a \(3\times3\) matrix.

The following example illustrates that \(AB = BA\) is not even guaranteed for two \(n\times n\) matrices \(A\) and \(B\):

\[\begin{split} \begin{bmatrix} 1 & 3 \\ 2 & 1 \end{bmatrix} \begin{bmatrix} 0 & 1 \\ 1 & 2 \end{bmatrix} = \begin{bmatrix} 3 & 7 \\ 1 & 4 \end{bmatrix} \neq \begin{bmatrix} 2 & 1 \\ 5 & 5 \end{bmatrix} = \begin{bmatrix} 0 & 1 \\ 1 & 2 \end{bmatrix} \begin{bmatrix} 1 & 3 \\ 2 & 1 \end{bmatrix}. \end{split}\]

The fact that \(AB \neq BA\) can be understood by thinking about the composition of the two transformations corresponding to \(A\) and \(B\). (See Section Linear Transformations.)

The following two exercises shed some light on the non-commutativity.

Example 3.2.14

Consider the two matrices

\[\begin{split} A = \begin{bmatrix} 2 & 0 \\ 0 & 1 \end{bmatrix} \quad \text{and} \quad B = \begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix} \end{split}\]

and the corresponding linear transformations

\[ T_A: \mathbb{R}^2 \to \mathbb{R}^2, \quad \mathbf{x} \mapsto T_A(\mathbf{x}) = A \mathbf{x} \]

and

\[ T_B: \mathbb{R}^2 \to \mathbb{R}^2, \quad \mathbf{x} \mapsto T_B(\mathbf{x}) = B \mathbf{x}. \]

We get

\[\begin{split} \mathbf{x} = \begin{bmatrix} x_1\\ x_2 \end{bmatrix} \quad \mapsto \quad A\mathbf{x} = \begin{bmatrix} 2 & 0 \\ 0 & 1 \end{bmatrix} \begin{bmatrix} x_1\\ x_2 \end{bmatrix} \,\, = \,\, \begin{bmatrix} 2x_1\\ x_2 \end{bmatrix} \end{split}\]

and likewise

\[\begin{split} T_B(\mathbf{x}) = \begin{bmatrix} 0 & 1 \\ 1 & 0 \end{bmatrix} \begin{bmatrix} x_1\\ x_2 \end{bmatrix} = \begin{bmatrix} x_2\\ x_1 \end{bmatrix}. \end{split}\]

Fig. 3.2.1 \( \begin{bmatrix} 2 & 0 \\ 0 & 1 \end{bmatrix}\begin{bmatrix} 0 & 1 \\ 1&0 \end{bmatrix} \neq \begin{bmatrix} 0 & 1 \\ 1&0 \end{bmatrix}\begin{bmatrix} 2 & 0 \\ 0 & 1 \end{bmatrix}\)

Note that \(T_A\) is a transformation that ‘stretches’ horizontally, and \(T_B\) is a reflection. Figure 3.2.1 visualizes the transformations corresponding to \(AB\) and \(BA\). When we apply the transformations one after another, the order in which we do this is important.

Exercise 3.2.10

Another way to understand why \(AB\neq BA\) is the following.

Recall that the matrices

\[\begin{split} E_1 = \begin{bmatrix} 1 & 0 \\ 2 & 1 \end{bmatrix} \quad \text{and} \quad E_2 = \begin{bmatrix} 3 & 0 \\ 0 & 1 \end{bmatrix} \end{split}\]

perform row operations, when multiplied with a \(2 \times n\) matrix \(A\).

• Describe in words the row operations corresponding to \(E_1\) and \(E_2\).

• Describe in words the combined row operations corresponding to \(E_1E_2\) and \(E_2E_1\). Can you explain why \(E_1E_2 \neq E_2E_1\)?

• Compute \(E_1E_2\) and \(E_2E_1\) to double check the last non-identity.

The second major difference between the product of numbers and the product of matrices: for two (e.g. real) numbers \(a\) and \(b\) it is known that

\[ \text{if} \quad a \neq 0 \quad \text{and} \quad b \neq 0\quad \text{then} \quad ab \neq 0, \]

or, equivalently,

\[ ab = 0 \quad \Rightarrow \quad a = 0 \,\,\text{ or } \,\,b = 0. \]

As the following example shows, things are different in the realm of matrices.

Example 3.2.15

\[\begin{split} \begin{bmatrix} 1 & 2 \\ 2 & 4 \end{bmatrix} \begin{bmatrix} 2 & 6 \\ -1 & -3 \end{bmatrix} = \begin{bmatrix} 0 & 0 \\ 0 & 0 \end{bmatrix}. \end{split}\]

So the product of two nonzero matrices may be the zero matrix.

The following example shows that things are even ‘worse’:

Example 3.2.16

\[\begin{split} \begin{bmatrix} 1 & -3 & 2 \\ 1 & -3 & 2 \\ 1 & -3 & 2 \end{bmatrix} \begin{bmatrix}1 & -3 & 2 \\ 1 & -3 & 2 \\ 1 & -3 & 2 \end{bmatrix} = \begin{bmatrix} 0 & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix}, \end{split}\]

which shows that we cannot even conclude from \(A\cdot A = O\) that \(A\) itself must be the zero matrix.

And here is another example of a nonzero matrix whose square is the zero matrix. In this case it can be seen geometrically what is going on.

Example 3.2.17

For the matrix \(A = \begin{bmatrix} 0 & 1 \\ 0 & 0 \end{bmatrix}\) we again have that

\[\begin{split} A^2 = \begin{bmatrix} 0 & 1 \\ 0 & 0 \end{bmatrix}^2 = \begin{bmatrix} 0 & 0 \\ 0 & 0 \end{bmatrix}. \end{split}\]

It holds that

\[\begin{split} A = A_2A_1 = \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix}\begin{bmatrix} 0 & 0 \\ 0 & 1 \end{bmatrix} \end{split}\]

Now consider the transformations corresponding to the matrices \(A_1\) and \(A_2\).

\[\begin{split} T_1(\vect{x}) = A_1\vect{x} = \begin{bmatrix} 0 & 0 \\ 0 & 1 \end{bmatrix}\begin{bmatrix} x_1 \\ x_2 \end{bmatrix} = \begin{bmatrix} 0\\ x_2 \end{bmatrix} \end{split}\]

is the projection onto the \(x_2\)-axis, and

\[\begin{split} T_2(\vect{x}) = A_2\vect{x} = \begin{bmatrix} 0 & 1 \\ -1 & 0 \end{bmatrix}\begin{bmatrix} x_1 \\ x_2 \end{bmatrix} = \begin{bmatrix} x_2\\ -x_1 \end{bmatrix} \end{split}\]

is the clockwise rotation about an angle \(\frac12\pi\).

Now let us see, step by step, what is the effect of the transformation \(T_2T_1T_2T_1\), corresponding to \(A^2\).

An arbitrary vector \(\vect{x}\) is sent to a vector \(T_1(\vect{x})\) on the \(x_2\) axis by \(T_1\).

The rotation sends this to a vector \(T_2(T_1(\vect{x}))\) on the \(x_1\)-axis. Projecting onto the \(x_2\)-axis again, will bring this to \(T_1(T_2(T_1(\vect{x}))) = \vect{0}\).
Lastly, applying \(T_2\) leaves the zero vector where it is.

So \(A^2\vect{x} = T_2(T_1(T_2(T_1(\vect{x})))) = \vect{0}\), for each vector \(\vect{x}\).

See Figure 3.2.2.

Fig. 3.2.2 Visualisation of \(\vect{x} \mapsto A^2\vect{x}\).

Remark 3.2.5

The next list gives six situations where matrix multiplication acts differently than multiplication of numbers.
In fact, all statements can be related to one of the first two.

1. In general, \(AB = BA\) does not hold for two \(n\times n\) matrices \(A\) and \(B\).

2. In general, from \(AB = O\) it does not follow that either \(A =O\) or \(B = O\).

3. In general, \((A+B)(A+B) = A^2 + 2AB + B^2\) does not hold for two \(n\times n\) matrices \(A\) and \(B\).

4. In general, \((A+B)(A-B) = A^2 - B^2\) does not hold for two \(n\times n\) matrices \(A\) and \(B\).

5. In general, from \(AB = AC\) and \(A \neq O\) it does not follow that \(B = C\).

6. In general, from \(A^2 = I\) it does not follow that either \(A = I\) or \(A = -I\).

For each statement counterexamples can be given, as we already did for the first two. To get more insight in what is really going on, we can also try to find out how the third till the sixth statement relate to the first two statements.

For instance, the third statement is closely related to the first. Let us check where ‘things go wrong’.

\[\begin{split} \begin{array}{cl} (A+B)(A+B)& = A(A+B) +B(A+B)\\ & = A^2 + AB + BA + B^2 \end{array} \end{split}\]

The last expression is only equal to

\[ A^2 + 2AB + B^2 \]

if

\[ AB + BA = 2AB \]

And that is only the case if

\[ BA = AB. \]

So any pair of two matrices \(A\) and \(B\) with

\[ AB \neq BA \]

provides a counterexample where

\[ (A+B)(A+B) \neq A^2 + 2AB + B^2. \]

Likewise, v. follows from ii. Namely,

\[ AB = AC \iff AB - AC = O \iff A(B-C) = O. \]

According to ii. from the last equation we cannot deduce that

\[ \text{either } A = O \quad \text{or}\quad B-C = O. \]

We can create a counterexample by taking for \(A\) and \(B\) nonzero matrices for which

\[ AB = O, \]

and we let \(C\) be the zero matrix. Then \(B \neq C\), whereas

\[ AB = AC = O \quad \text{and (by assumption)} \quad A \neq O. \]

Statement vi. also relates to ii. Namely,

\[ A^2 = I \quad \iff \quad A^2 - I = (A+I)(A-I) = O. \]

From the last equality we cannot conclude that one of the factors \((A+I)\) or \((A-I)\) must be the zero matrix. In this case we do not get a counterexample for free. You are asked to construct counterexamples in Exercise 3.2.11.

Exercise 3.2.11

1. Give a \(2 \times 2\) matrix \(A \neq \pm\, I\) for which \(A^2 = I\).

2. Give a \(2 \times 2\) matrix \(A\) not containing any zeros, for which \(A^2 = I\).

3. Give a \(2 \times 2\) matrix \(B\) for which \(B^2 = -I\).

The following property connects the two operations matrix transposition and matrix multiplication.

Proposition 3.2.6

If \(A\) is an \(m\times n\) matrix and \(B\) an \(n\times p\) matrix, then

\[ (AB)^T = B^TA^T. \]

Before we present the proof, we consider a typical example.

Example 3.2.18

We verify the rule for the two matrices

\[\begin{split} A = \begin{bmatrix} 2 & 1 & -1 \\ 1 & -1 & 3 \end{bmatrix} \quad\text{and}\quad B = \begin{bmatrix} 1 & -3 & 0\\ 4 & 2 & -1 \\ 5 & 2 & 1\end{bmatrix}. \end{split}\]

We compute:

\[\begin{split} AB = \begin{bmatrix} 2 & 1 & -1 \\ 1 & -1 & 3 \end{bmatrix} \begin{bmatrix} 1 & -3 & 0 \\ 4 & 2 & -1\\ 5 & 2 & 1 \end{bmatrix} = \begin{bmatrix} 1 & -6 & -2 \\ 12 & 1 & 4\end{bmatrix} \end{split}\]

and

\[\begin{split} B^TA^T = \begin{bmatrix} 1 & 4 & 5 \\ -3 & 2 & 2 \\ 0 & -1 & 1 \end{bmatrix} \begin{bmatrix} 2& 1 \\ 1 & -1 \\ -1 & 3 \end{bmatrix} = \begin{bmatrix} 1 & 12 \\ -6 & 1 \\ -2 & 4 \end{bmatrix}, \end{split}\]

so that indeed

\[\begin{split} B^TA^T = \begin{bmatrix} 1 & 12 \\ -6 & 1 \\ -2 & 4 \end{bmatrix} = \begin{bmatrix} 1 & -6 & -2 \\ 12 & 1 & 4\end{bmatrix}^T = (AB)^T. \end{split}\]

Careful inspection learns that for the two matrix products exactly the same sums and products of numbers have to be computed. For instance, in both products \(12\) is the sum of products

\[ 12 = 1\cdot1 +4\cdot(-1) +5\cdot3 = 1\cdot1 +(-1)\cdot4 +3\cdot5. \]

As Example 3.2.18 illustrates the rule is not restricted to square matrices \(A\) and \(B\). The proof for general matrices \(A\) and \(B\) for which the product \(AB\) is well defined is as follows

Proof. To show that

\[ (AB)^T = B^TA^T \]

we have to show that the matrices have the same size, and are equal entry by entry. First, we see that \(AB\) is an \(m \times p\) matrix, so \((AB)^T\) is a \(p \times m\) matrix, and \(B^TA^T\), being the product of a \(p \times n\) matrix with an \(n \times m\), is also a \(p \times m\) matrix.

Second, the \((i,j)\) entry of \((AB)^T\) is the \((j,i)\) entry of \(AB\), which is the (row-column) product of the \(j\)-th row of \(A\) and the \(i\)-th column of \(B\):

\[\begin{split} [(AB)^T]_{ij} = \begin{bmatrix} a_{j1} & a_{j2} & \ldots & a_{jn} \end{bmatrix}\begin{bmatrix} b_{1i} \\ b_{2i} \\ \vdots \\ b_{ni} \end{bmatrix}. \end{split}\]

The \((i,j)\) entry of \(B^TA^T\) is the product of the \(i\)-th row of \(B^T\) and the \(j\)-th column of \(A^T\).
Now the \(i\)-th row of \(B^T\) is the \(i\)-th column of \(B\) written as a row, and the \(j\)-th column of \(A^T\) is the \(j\)-th row of \(A\) written as a column:

\[\begin{split} [B^TA^T]_{ij} = \begin{bmatrix} b_{1i} & b_{2i} & \ldots & b_{ni} \end{bmatrix}\begin{bmatrix} a_{j1} \\ a_{j2} \\ \vdots \\ a_{jn} \end{bmatrix}. \end{split}\]

Both row-column products end up as the same value

\[ a_{j1}b_{1i} + a_{j2}b_{2i} + \ldots + a_{jn}b_{ni} = b_{1i}a_{j1} + b_{2i}a_{j2} + \ldots + b_{ni}a_{jn}. \]

Remark 3.2.6

We already defined \(A^2\) for a square matrix \(A\). We can extend this to higher powers of \(A\) in an obvious way:

\[ A^3 = A(A^2),\quad A^4 = A(A^3), \,\,.\,.\,.\,. \]

Since

\[ A(A^2) = A(AA) = (AA)A, \]

we can do without the parentheses.

For the same reason

\[ A^kA^{\ell} = A^{k+\ell}, \quad \text{for integers} \quad k,\ell \geq 1. \]

If we define

\[ A^0 = I, \]

then

\[ A^kA^{\ell} = A^{k+\ell} \quad \text{holds for all integers} \quad \quad k,\ell \geq 0. \]

And what can we say about \(A^{-1}\)?

We will dedicate Section 3.4 to this topic.

3.2.6. Grasple Exercises

Grasple Exercise 3.2.6

To compute a product \(AB\)

Show/Hide Content

Grasple Exercise 3.2.7

To compute a product \(AB\)

Show/Hide Content

Grasple Exercise 3.2.8

To compute a product \(AB\)

Show/Hide Content

Grasple Exercise 3.2.9

To compute a product \(AB\)

Show/Hide Content

Grasple Exercise 3.2.10

To compute several matrix products.

Show/Hide Content

Grasple Exercise 3.2.11

To compute \(\vect{u}^T\vect{v}\) and \(\vect{u}\vect{v}^T\).

Show/Hide Content

Grasple Exercise 3.2.12

To find \(k\) for which \(AB=BA\).

Show/Hide Content

Grasple Exercise 3.2.13

To find \(k\) for which \(AB=BA\).

Show/Hide Content

Grasple Exercise 3.2.14

To find two products \(AD_1\) and \(D_2A\).

Show/Hide Content

Grasple Exercise 3.2.15

To find a high power of a special matrix

Show/Hide Content

Grasple Exercise 3.2.16

To find a high power of a special matrix.

Show/Hide Content

The remaining exercises are less of a compuational character.

Grasple Exercise 3.2.17

Is the zero matrix a diagonal matrix?

Show/Hide Content

Grasple Exercise 3.2.18

To explain why a certain product does not exist.

Show/Hide Content

Grasple Exercise 3.2.19

To find a matrix \(A \neq \pm I\) for which \(A^2 = I\).

Show/Hide Content

Grasple Exercise 3.2.20

To show that \((cA)^T = cA^T\).

Show/Hide Content

Grasple Exercise 3.2.21

To show: \(A^TA = D \iff A\) has orthogonal columns.

Show/Hide Content

Grasple Exercise 3.2.22

To find the size of \(C\) if \(AC = B\).

Show/Hide Content

Grasple Exercise 3.2.23

Number of columns of \(C\) is \(AC=B\).

Show/Hide Content

Grasple Exercise 3.2.24

To find the number of rows of \(B\) if \(BC\) is an \(m\times n\) matrix.

Show/Hide Content

Grasple Exercise 3.2.25

Finding \(E\) such that \(EA = M\)   (or \(AE = M\)).

Show/Hide Content

Grasple Exercise 3.2.26

A bit like the previous one ‘by inspection’.

Show/Hide Content

Grasple Exercise 3.2.27

Two True/False questions about products and symmmetric matrices.

Show/Hide Content