Skip to main content

## SectionMRMatrix Representations

We have seen that linear transformations whose domain and codomain are vector spaces of columns vectors have a close relationship with matrices (Theorem MBLT, Theorem MLTCV). In this section, we will extend the relationship between matrices and linear transformations to the setting of linear transformations between abstract vector spaces.

### SubsectionMRMatrix Representations

This is a fundamental definition.

#### DefinitionMR.Matrix Representation.

Suppose that $$\ltdefn{T}{U}{V}$$ is a linear transformation, $$B=\set{\vectorlist{u}{n}}$$ is a basis for $$U$$ of size $$n\text{,}$$ and $$C$$ is a basis for $$V$$ of size $$m\text{.}$$ Then the matrix representation of $$T$$ relative to $$B$$ and $$C$$ is the $$m\times n$$ matrix,

\begin{equation*} \matrixrep{T}{B}{C}=\left[ \left.\vectrep{C}{\lteval{T}{\vect{u}_1}}\right| \left.\vectrep{C}{\lteval{T}{\vect{u}_2}}\right| \left.\vectrep{C}{\lteval{T}{\vect{u}_3}}\right| \ldots \left|\vectrep{C}{\lteval{T}{\vect{u}_n}}\right. \right]\text{.} \end{equation*}

Consider the linear transformation, $$\ltdefn{S}{P_3}{M_{22}}\text{,}$$ given by

\begin{align*} \lteval{S}{a+bx+cx^2+dx^3}= \begin{bmatrix} 3a+7b-2c-5d & 8a+14b-2c-11d\\ -4a-8b+2c+6d & 12a+22b-4c-17d \end{bmatrix}\text{.} \end{align*}

First, we build a representation relative to the bases,

\begin{align*} B&=\set{1+2x+x^2-x^3,\,1+3x+x^2+x^3,\,-1-2x+2x^3,\,2+3x+2x^2-5x^3}\\ C&=\set{ \begin{bmatrix}1 & 1 \\ 1 & 2\end{bmatrix},\, \begin{bmatrix}2 & 3 \\ 2 & 5\end{bmatrix},\, \begin{bmatrix}-1 & -1 \\ 0 & -2\end{bmatrix},\, \begin{bmatrix}-1 & -4 \\ -2 & -4\end{bmatrix} }\text{.} \end{align*}

We evaluate $$S$$ with each element of the basis for the domain, $$B\text{,}$$ and coordinatize the result relative to the vectors in the basis for the codomain, $$C\text{.}$$ Notice here how we take elements of vector spaces and decompose them into linear combinations of basis elements as the key step in constructing coordinatizations of vectors. There is a system of equations involved almost every time, but we will omit these details since this should be a routine exercise at this stage. We have

\begin{align*} &\vectrep{C}{\lteval{S}{1+2x+x^2-x^3}} =\vectrep{C}{\begin{bmatrix}20 & 45 \\ -24 & 69\end{bmatrix}}\\ &\quad=\vectrep{C}{ (-90)\begin{bmatrix}1 & 1 \\ 1 & 2\end{bmatrix}+ 37\begin{bmatrix}2 & 3 \\ 2 & 5\end{bmatrix}+ (-40)\begin{bmatrix}-1 & -1 \\ 0 & -2\end{bmatrix}+ 4\begin{bmatrix}-1 & -4 \\ -2 & -4\end{bmatrix} } =\colvector{-90\\37\\-40\\4}\\ &\vectrep{C}{\lteval{S}{1+3x+x^2+x^3}} =\vectrep{C}{\begin{bmatrix}17 & 37 \\ -20 & 57\end{bmatrix}}\\ &\quad=\vectrep{C}{ (-72)\begin{bmatrix}1 & 1 \\ 1 & 2\end{bmatrix}+ 29\begin{bmatrix}2 & 3 \\ 2 & 5\end{bmatrix}+ (-34)\begin{bmatrix}-1 & -1 \\ 0 & -2\end{bmatrix}+ 3\begin{bmatrix}-1 & -4 \\ -2 & -4\end{bmatrix} } =\colvector{-72\\29\\-34\\3}\\ &\vectrep{C}{\lteval{S}{-1-2x+2x^3}} =\vectrep{C}{\begin{bmatrix}-27 & -58 \\ 32 & -90\end{bmatrix}}\\ &\quad=\vectrep{C}{ 114\begin{bmatrix}1 & 1 \\ 1 & 2\end{bmatrix}+ (-46)\begin{bmatrix}2 & 3 \\ 2 & 5\end{bmatrix}+ 54\begin{bmatrix}-1 & -1 \\ 0 & -2\end{bmatrix}+ (-5)\begin{bmatrix}-1 & -4 \\ -2 & -4\end{bmatrix} } =\colvector{114\\-46\\54\\-5}\\ &\vectrep{C}{\lteval{S}{2+3x+2x^2-5x^3}} =\vectrep{C}{\begin{bmatrix}48 & 109 \\ -58 & 167\end{bmatrix}}\\ &\quad=\vectrep{C}{ -220\begin{bmatrix}1 & 1 \\ 1 & 2\end{bmatrix}+ 91\begin{bmatrix}2 & 3 \\ 2 & 5\end{bmatrix}+ (-96)\begin{bmatrix}-1 & -1 \\ 0 & -2\end{bmatrix}+ 10\begin{bmatrix}-1 & -4 \\ -2 & -4\end{bmatrix} } =\colvector{-220\\91\\-96\\10}\text{.} \end{align*}

Thus, employing Definition MR

\begin{equation*} \matrixrep{S}{B}{C}= \begin{bmatrix} -90 & -72 & 114 & -220 \\ 37 & 29 & -46 & 91 \\ -40 & -34 & 54 & -96 \\ 4 & 3 & -5 & 10 \end{bmatrix}\text{.} \end{equation*}

Often we use “nice” bases to build matrix representations and the work involved is much easier. Suppose we take bases

\begin{align*} D&=\set{1,\,x,\,x^2,\,x^3} & E&=\set{ \begin{bmatrix}1&0\\0&0\end{bmatrix},\, \begin{bmatrix}0&1\\0&0\end{bmatrix},\, \begin{bmatrix}0&0\\1&0\end{bmatrix},\, \begin{bmatrix}0&0\\0&1\end{bmatrix}}\text{.} \end{align*}

The evaluation of $$S$$ at the elements of $$D$$ is easy and coordinatization relative to $$E$$ can be done on sight,

\begin{align*} \vectrep{E}{\lteval{S}{1}} &=\vectrep{E}{\begin{bmatrix}3 & 8 \\ -4 & 12\end{bmatrix}}\\ &=\vectrep{E}{ 3\begin{bmatrix}1&0\\0&0\end{bmatrix}+ 8\begin{bmatrix}0&1\\0&0\end{bmatrix}+ (-4)\begin{bmatrix}0&0\\1&0\end{bmatrix}+ 12\begin{bmatrix}0&0\\0&1\end{bmatrix} } =\colvector{3 \\ 8 \\ -4 \\ 12}\\ \vectrep{E}{\lteval{S}{x}} &=\vectrep{E}{\begin{bmatrix}7 & 14 \\ -8 & 22\end{bmatrix}}\\ &=\vectrep{E}{ 7\begin{bmatrix}1&0\\0&0\end{bmatrix}+ 14\begin{bmatrix}0&1\\0&0\end{bmatrix}+ (-8)\begin{bmatrix}0&0\\1&0\end{bmatrix}+ 22\begin{bmatrix}0&0\\0&1\end{bmatrix} } =\colvector{7 \\ 14 \\ -8 \\ 22}\\ \vectrep{E}{\lteval{S}{x^2}} &=\vectrep{E}{\begin{bmatrix}-2 & -2 \\ 2 & -4\end{bmatrix}}\\ &=\vectrep{E}{ (-2)\begin{bmatrix}1&0\\0&0\end{bmatrix}+ (-2)\begin{bmatrix}0&1\\0&0\end{bmatrix}+ 2\begin{bmatrix}0&0\\1&0\end{bmatrix}+ (-4)\begin{bmatrix}0&0\\0&1\end{bmatrix} } =\colvector{-2 \\ -2 \\ 2 \\ -4}\\ \vectrep{E}{\lteval{S}{x^3}} &=\vectrep{E}{\begin{bmatrix}-5 & -11 \\ 6 & -17\end{bmatrix}}\\ &=\vectrep{E}{ (-5)\begin{bmatrix}1&0\\0&0\end{bmatrix}+ (-11)\begin{bmatrix}0&1\\0&0\end{bmatrix}+ 6\begin{bmatrix}0&0\\1&0\end{bmatrix}+ (-17)\begin{bmatrix}0&0\\0&1\end{bmatrix} }\\ &=\colvector{-5 \\ -11 \\ 6 \\ -17}\text{.} \end{align*}

So the matrix representation of $$S$$ relative to $$D$$ and $$E$$ is

\begin{equation*} \matrixrep{S}{D}{E}= \begin{bmatrix} 3 & 7 & -2 & -5 \\ 8 & 14 & -2 & -11 \\ -4 & -8 & 2 & 6 \\ 12 & 22 & -4 & -17 \end{bmatrix}\text{.} \end{equation*}

One more time, but now let us use bases

\begin{align*} F&=\set{1+x-x^2+2x^3,\,-1+2x+2x^3,\,2+x-2x^2+3x^3,\,1+x+2x^3}\\ G&=\set{ \begin{bmatrix}1&1\\-1&2\end{bmatrix},\, \begin{bmatrix}-1&2\\0&2\end{bmatrix},\, \begin{bmatrix}2&1\\-2&3\end{bmatrix},\, \begin{bmatrix}1&1\\0&2\end{bmatrix} } \end{align*}

and evaluate $$S$$ with the elements of $$F\text{,}$$ then coordinatize the results relative to $$G\text{,}$$

\begin{align*} \vectrep{G}{\lteval{S}{1+x-x^2+2x^3}} &=\vectrep{G}{\begin{bmatrix}2 & 2 \\ -2 & 4\end{bmatrix}} =\vectrep{G}{ 2\begin{bmatrix}1&1\\-1&2\end{bmatrix} } =\colvector{2\\0\\0\\0}\\ \vectrep{G}{\lteval{S}{-1+2x+2x^3}} &=\vectrep{G}{\begin{bmatrix}1 & -2 \\ 0 & -2\end{bmatrix}} =\vectrep{G}{ (-1)\begin{bmatrix}-1&2\\0&2\end{bmatrix} } =\colvector{0\\-1\\0\\0}\\ \vectrep{G}{\lteval{S}{2+x-2x^2+3x^3}} &=\vectrep{G}{\begin{bmatrix}2 & 1 \\ -2 & 3\end{bmatrix}} =\vectrep{G}{ \begin{bmatrix}2&1\\-2&3\end{bmatrix} } =\colvector{0\\0\\1\\0}\\ \vectrep{G}{\lteval{S}{1+x+2x^3}} &=\vectrep{G}{\begin{bmatrix}0 & 0 \\ 0 & 0\end{bmatrix}} =\vectrep{G}{ 0\begin{bmatrix}1&1\\0&2\end{bmatrix} } =\colvector{0\\0\\0\\0}\text{.} \end{align*}

So we arrive at an especially economical matrix representation,

\begin{equation*} \matrixrep{S}{F}{G}= \begin{bmatrix} 2 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 \end{bmatrix}\text{.} \end{equation*}

We may choose to use whatever terms we want when we make a definition. Some are arbitrary, while others make sense, but only in light of subsequent theorems. Matrix representation is in the latter category. We begin with a linear transformation and produce a matrix. So what? Here is the theorem that justifies the term matrix representation.

Let $$B=\set{\vectorlist{u}{n}}$$ be the basis of $$U\text{.}$$ Since $$\vect{u}\in U\text{,}$$ there are scalars $$\scalarlist{a}{n}$$ such that

\begin{equation*} \vect{u}=\lincombo{a}{u}{n}\text{.} \end{equation*}

Then,

\begin{align*} &\matrixrep{T}{B}{C}\vectrep{B}{\vect{u}}\\ &= \left[ \left.\vectrep{C}{\lteval{T}{\vect{u}_1}}\right| \left.\vectrep{C}{\lteval{T}{\vect{u}_2}}\right| \left.\vectrep{C}{\lteval{T}{\vect{u}_3}}\right| \ldots \left|\vectrep{C}{\lteval{T}{\vect{u}_n}}\right. \right] \vectrep{B}{\vect{u}}&& \knowl{./knowl/definition-MR.html}{\text{Definition MR}}\\ &= \left[ \left.\vectrep{C}{\lteval{T}{\vect{u}_1}}\right| \left.\vectrep{C}{\lteval{T}{\vect{u}_2}}\right| \left.\vectrep{C}{\lteval{T}{\vect{u}_3}}\right| \ldots \left|\vectrep{C}{\lteval{T}{\vect{u}_n}}\right. \right] \colvector{a_1\\a_2\\a_3\\\vdots\\a_n}&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\\ &= a_1\vectrep{C}{\lteval{T}{\vect{u}_1}}+ a_2\vectrep{C}{\lteval{T}{\vect{u}_2}}+ \cdots+ a_n\vectrep{C}{\lteval{T}{\vect{u}_n}}&& \knowl{./knowl/definition-MVP.html}{\text{Definition MVP}}\\ &= \vectrep{C}{ a_1\lteval{T}{\vect{u}_1}+ a_2\lteval{T}{\vect{u}_2}+ a_3\lteval{T}{\vect{u}_3}+\cdots+ a_n\lteval{T}{\vect{u}_n} }&& \knowl{./knowl/theorem-LTLC.html}{\text{Theorem LTLC}}\\ &= \vectrep{C}{ \lteval{T}{ \lincombo{a}{u}{n} }}&& \knowl{./knowl/theorem-LTLC.html}{\text{Theorem LTLC}}\\ &= \vectrep{C}{\lteval{T}{\vect{u}}}\text{.} \end{align*}

The alternative conclusion is obtained as

\begin{align*} \lteval{T}{\vect{u}} &=\lteval{I_V}{\lteval{T}{\vect{u}}}&& \knowl{./knowl/definition-IDLT.html}{\text{Definition IDLT}}\\ &=\lteval{\left(\compose{\vectrepinvname{C}}{\vectrepname{C}}\right)}{\lteval{T}{\vect{u}}}&& \knowl{./knowl/definition-IVLT.html}{\text{Definition IVLT}}\\ &=\vectrepinv{C}{\vectrep{C}{\lteval{T}{\vect{u}}}}&& \knowl{./knowl/definition-LTC.html}{\text{Definition LTC}}\\ &=\vectrepinv{C}{\matrixrep{T}{B}{C}\left(\vectrep{B}{\vect{u}}\right)}\text{.} \end{align*}

This theorem says that we can apply $$T$$ to $$\vect{u}$$ and coordinatize the result relative to $$C$$ in $$V\text{,}$$ or we can first coordinatize $$\vect{u}$$ relative to $$B$$ in $$U\text{,}$$ then multiply by the matrix representation. Either way, the result is the same. So the effect of a linear transformation can always be accomplished by a matrix-vector product (Definition MVP). That is important enough to say again. The effect of a linear transformation is a matrix-vector product.

The alternative conclusion of this result might be even more striking. It says that to effect a linear transformation ($$T$$) of a vector ($$\vect{u}$$), coordinatize the input (with $$\vectrepname{B}$$), do a matrix-vector product (with $$\matrixrep{T}{B}{C}$$), and un-coordinatize the result (with $$\vectrepinvname{C}$$). So, absent some bookkeeping about vector representations, a linear transformation is a matrix. To adjust the diagram, we “reverse” the arrow on the right, which means inverting the vector representation $$\vectrepname{C}$$ on $$V\text{.}$$ Now we can go directly across the top of the diagram, computing the linear transformation between the abstract vector spaces. Or, we can around the other three sides, using vector representation, a matrix-vector product, followed by un-coordinatization.

Here is an example to illustrate how the “action” of a linear transformation can be effected by matrix multiplication.

In Example OLTTR we found three representations of the linear transformation $$S\text{.}$$ In this example, we will compute a single output of $$S$$ in four different ways. First “normally,” then three times over using Theorem FTMR.

Choose $$p(x)=3-x+2x^2-5x^3\text{,}$$ for no particular reason. Then the straightforward application of $$S$$ to $$p(x)$$ yields

\begin{align*} \lteval{S}{p(x)} &=\lteval{S}{3-x+2x^2-5x^3}\\ &= \begin{bmatrix} 3(3)+7(-1)-2(2)-5(-5) & 8(3)+14(-1)-2(2)-11(-5)\\ -4(3)-8(-1)+2(2)+6(-5) & 12(3)+22(-1)-4(2)-17(-5) \end{bmatrix}\\ &= \begin{bmatrix} 23 & 61 \\ -30 & 91 \end{bmatrix}\text{.} \end{align*}

Now use the representation of $$S$$ relative to the bases $$B$$ and $$C$$ and Theorem FTMR. Note that we will employ the following linear combination in moving from the second line to the third,

\begin{align*} 3-x+2x^2-5x^3 &= 48(1+2x+x^2-x^3)+(-20)(1+3x+x^2+x^3)+\\ &\quad\quad (-1)(-1-2x+2x^3)+(-13)(2+3x+2x^2-5x^3)\text{.} \end{align*}

We have

\begin{align*} \lteval{S}{p(x)} &=\vectrepinv{C}{\matrixrep{S}{B}{C}\vectrep{B}{p(x)}}\\ &=\vectrepinv{C}{\matrixrep{S}{B}{C}\vectrep{B}{3-x+2x^2-5x^3}}\\ &=\vectrepinv{C}{\matrixrep{S}{B}{C}\colvector{48\\-20\\-1\\-13}}\\ &=\vectrepinv{C}{ \begin{bmatrix} -90 & -72 & 114 & -220 \\ 37 & 29 & -46 & 91 \\ -40 & -34 & 54 & -96 \\ 4 & 3 & -5 & 10 \end{bmatrix} \colvector{48\\-20\\-1\\-13}}\\ &=\vectrepinv{C}{\colvector{-134\\59\\-46\\7}}\\ &= (-134)\begin{bmatrix}1 & 1 \\ 1 & 2\end{bmatrix}+ 59\begin{bmatrix}2 & 3 \\ 2 & 5\end{bmatrix}+ (-46)\begin{bmatrix}-1 & -1 \\ 0 & -2\end{bmatrix}+ 7\begin{bmatrix}-1 & -4 \\ -2 & -4\end{bmatrix}\\ &= \begin{bmatrix} 23 & 61 \\ -30 & 91 \end{bmatrix}\text{.} \end{align*}

Again, but now with “nice” bases like $$D$$ and $$E\text{,}$$ and the computations are more transparent. We have

\begin{align*} \lteval{S}{p(x)} &=\vectrepinv{E}{\matrixrep{S}{D}{E}\vectrep{D}{p(x)}}\\ &=\vectrepinv{E}{\matrixrep{S}{D}{E}\vectrep{D}{3-x+2x^2-5x^3}}\\ &=\vectrepinv{E}{\matrixrep{S}{D}{E}\vectrep{D}{3(1)+(-1)(x)+2(x^2)+(-5)(x^3)}}\\ &=\vectrepinv{E}{\matrixrep{S}{D}{E}\colvector{3\\-1\\2\\-5}}\\ &=\vectrepinv{E}{ \begin{bmatrix} 3 & 7 & -2 & -5 \\ 8 & 14 & -2 & -11 \\ -4 & -8 & 2 & 6 \\ 12 & 22 & -4 & -17 \end{bmatrix} \colvector{3\\-1\\2\\-5} }\\ &=\vectrepinv{E}{\colvector{23 \\ 61 \\ -30 \\ 91}}\\ &= 23\begin{bmatrix}1&0\\0&0\end{bmatrix}+ 61\begin{bmatrix}0&1\\0&0\end{bmatrix}+ (-30)\begin{bmatrix}0&0\\1&0\end{bmatrix}+ 91\begin{bmatrix}0&0\\0&1\end{bmatrix}\\ &= \begin{bmatrix} 23 & 61 \\ -30 & 91 \end{bmatrix}\text{.} \end{align*}

OK, last time, now with the bases $$F$$ and $$G\text{.}$$ The coordinatizations will take some work this time, but the matrix-vector product (Definition MVP) (which is the actual action of the linear transformation) will be especially easy, given the diagonal nature of the matrix representation, $$\matrixrep{S}{F}{G}\text{.}$$ Here we go,

\begin{align*} \lteval{S}{p(x)} &=\vectrepinv{G}{\matrixrep{S}{F}{G}\vectrep{F}{p(x)}}\\ &=\vectrepinv{G}{\matrixrep{S}{F}{G}\vectrep{F}{3-x+2x^2-5x^3}}\\ &=\vectrepinvname{G}\left(\matrixrep{S}{F}{G}\vectrepname{F}\left( 32(1+x-x^2+2x^3) -7(-1+2x+2x^3)\right.\right.\\ &\left.\left.\quad\quad\quad\quad\quad\quad\quad -17(2+x-2x^2+3x^3) -2(1+x+2x^3) \right)\right)\\ &=\vectrepinv{G}{\matrixrep{S}{F}{G}\colvector{32\\-7\\-17\\-2}}\\ &=\vectrepinv{G}{ \begin{bmatrix} 2 & 0 & 0 & 0 \\ 0 & -1 & 0 & 0 \\ 0 & 0 & 1 & 0 \\ 0 & 0 & 0 & 0 \end{bmatrix} \colvector{32\\-7\\-17\\-2} }\\ &=\vectrepinv{G}{\colvector{64 \\ 7 \\ -17 \\ 0}}\\ &= 64\begin{bmatrix}1&1\\-1&2\end{bmatrix}+ 7\begin{bmatrix}-1&2\\0&2\end{bmatrix}+ (-17)\begin{bmatrix}2&1\\-2&3\end{bmatrix}+ 0\begin{bmatrix}1&1\\0&2\end{bmatrix}\\ &= \begin{bmatrix} 23 & 61 \\ -30 & 91 \end{bmatrix}\text{.} \end{align*}

This example is not meant to necessarily illustrate that any one of these four computations is simpler than the others. Instead, it is meant to illustrate the many different ways we can arrive at the same result, with the last three all employing a matrix representation to effect the linear transformation.

We will use Theorem FTMR frequently in the next few sections. A typical application will feel like the linear transformation $$T$$ “commutes” with a vector representation, $$\vectrepname{C}\text{,}$$ and as it does the transformation morphs into a matrix, $$\matrixrep{T}{B}{C}\text{,}$$ while the vector representation changes to a new basis, $$\vectrepname{B}\text{.}$$ Or vice-versa.

#### SageMR.Matrix Representations.

It is very easy to create matrix representations of linear transformations in Sage. Our examples in the text have used abstract vector spaces to make it a little easier to keep track of the domain and codomain. With Sage we will have to consistently use variants of QQ^n, but we will use non-obvious bases to make it nontrivial and interesting. Here are the components of our first example, one linear function, two vector spaces, four bases.

x1, x2, x3, x4 = var('x1, x2, x3, x4')
outputs = [3*x1 + 7*x2 + x3 - 4*x4,
2*x1 + 5*x2 + x3 - 3*x4,
-x1 - 2*x2      +   x4]
T_symbolic(x1, x2, x3, x4) = outputs
U = QQ^4
V = QQ^3
b0 = vector(QQ, [ 1, 1, -1, 0])
b1 = vector(QQ, [-1, 0, -2, 7])
b2 = vector(QQ, [ 0, 1, -2, 4])
b3 = vector(QQ, [-2, 0, -1, 6])
B = [b0, b1, b2, b3]
c0 = vector(QQ, [ 1,  6, -6])
c1 = vector(QQ, [ 0,  1, -1])
c2 = vector(QQ, [-2, -3,  4])
C = [c0, c1, c2]
d0 = vector(QQ, [ 1, -3,  2, -1])
d1 = vector(QQ, [ 0,  1,  0,  1])
d2 = vector(QQ, [-1,  2, -1, -1])
d3 = vector(QQ, [ 2, -8,  4, -3])
D = [d0, d1, d2, d3]
e0 = vector(QQ, [ 0,  1, -3])
e1 = vector(QQ, [-1,  2, -1])
e2 = vector(QQ, [ 2, -4,  3])
E = [e0, e1, e2]


We create alternate versions of the domain and codomain by providing them with bases other than the standard ones. Then we build the linear transformation and ask for its .matrix(). We will use numerals to distinguish our two examples.

U1 = U.subspace_with_basis(B)
V1 = V.subspace_with_basis(C)
T1 = linear_transformation(U1, V1, T_symbolic)
T1.matrix(side='right')


Now, we do it again, but with the bases Dand E.

U2 = U.subspace_with_basis(D)
V2 = V.subspace_with_basis(E)
T2 = linear_transformation(U2, V2, T_symbolic)
T2.matrix(side='right')


So once the pieces are in place, it is quite easy to obtain a matrix representation. You might experiment with examples from the text, by first mentally converting elements of abstract vector spaces into column vectors via vector representations using simple bases for the abstract spaces.

The linear transformations T1and T2are not different as functions, despite the fact that Sage treats them as unequal (since they have different matrix representations). We can check that the two functions behave identically, first with random testing. Repeated execution of the following compute cell should always produce True.

u = random_vector(QQ, 4)
T1(u) == T2(u)


A better approach would be to see if T1and T2agree on a basis for the domain, and to avoid any favoritism, we will use the basis of Uitself. Convince yourself that this is a proper application of Theorem LTDB to demonstrate equality.

all([T1(u) == T2(u) for u in U.basis()])


Or we can just ask if they are equal functions (a method that is implemented using the previous idea about agreeing on a basis).

T1.is_equal_function(T2)


We can also demonstrate Theorem FTMR — we will use the representation for T1. We need a couple of vector representation linear transformations.

rhoB = linear_transformation(U1, U, U.basis())
rhoC = linear_transformation(V1, V, V.basis())


Now we experiment with a “random” element of the domain. Study the third computation carefully, as it is the Sage version of Theorem FTMR. You might compute some of the pieces of this expression to build up the final expression, and you might duplicate the experiment with a different input and/or with the T2version.

T_symbolic(-3, 4, -5, 2)

u = vector(QQ, [-3, 4, -5, 2])
T1(u)

(rhoC^-1)(T1.matrix(side='right')*rhoB(u))


One more time, but we will replace the vector representation linear transformations with Sage conveniences. Recognize that the .linear_combination_of_basis()method is the inverse of vector representation (it is “un-coordinatization”).

output_coord = T1.matrix(side='right')*U1.coordinate_vector(u)
V1.linear_combination_of_basis(output_coord)


We have concentrated here on the first version of the conclusion of Theorem FTMR. You could experiment with the second version in a similar manner. Extra credit: what changes do you need to make in any of these experiments if you remove the side='right'option?

### SubsectionNRFONew Representations from Old

In Subsection LT.NLTFO we built new linear transformations from other linear transformations. Sums, scalar multiples and compositions. These new linear transformations will have matrix representations as well. How do the new matrix representations relate to the old matrix representations? Here are the three theorems.

Let $$\vect{x}$$ be any vector in $$\complex{n}\text{.}$$ Define $$\vect{u}\in U$$ by $$\vect{u}=\vectrepinv{B}{\vect{x}}\text{,}$$ so $$\vect{x}=\vectrep{B}{\vect{u}}\text{.}$$ Then,

\begin{align*} \matrixrep{T+S}{B}{C}\vect{x}&= \matrixrep{T+S}{B}{C}\vectrep{B}{\vect{u}}&&\text{Substitution}\\ &=\vectrep{C}{\lteval{\left(T+S\right)}{\vect{u}}}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrep{C}{\lteval{T}{\vect{u}}+\lteval{S}{\vect{u}}}&& \knowl{./knowl/definition-LTA.html}{\text{Definition LTA}}\\ &=\vectrep{C}{\lteval{T}{\vect{u}}}+\vectrep{C}{\lteval{S}{\vect{u}}}&& \knowl{./knowl/definition-LT.html}{\text{Definition LT}}\\ &=\matrixrep{T}{B}{C}\left(\vectrep{B}{\vect{u}}\right)+\matrixrep{S}{B}{C}\left(\vectrep{B}{\vect{u}}\right)&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\left(\matrixrep{T}{B}{C}+\matrixrep{S}{B}{C}\right)\vectrep{B}{\vect{u}}&& \knowl{./knowl/theorem-MMDAA.html}{\text{Theorem MMDAA}}\\ &=\left(\matrixrep{T}{B}{C}+\matrixrep{S}{B}{C}\right)\vect{x}&&\text{Substitution}\text{.} \end{align*}

Since the matrices $$\matrixrep{T+S}{B}{C}$$ and $$\matrixrep{T}{B}{C}+\matrixrep{S}{B}{C}$$ have equal matrix-vector products for every vector in $$\complex{n}\text{,}$$ by Theorem EMMVP they are equal matrices. (Now would be a good time to double-back and study the proof of Theorem EMMVP. You did promise to come back to this theorem sometime, didn't you?)

Let $$\vect{x}$$ be any vector in $$\complex{n}\text{.}$$ Define $$\vect{u}\in U$$ by $$\vect{u}=\vectrepinv{B}{\vect{x}}\text{,}$$ so $$\vect{x}=\vectrep{B}{\vect{u}}\text{.}$$ Then,

\begin{align*} \matrixrep{\alpha T}{B}{C}\vect{x}&= \matrixrep{\alpha T}{B}{C}\vectrep{B}{\vect{u}}&&\text{Substitution}\\ &=\vectrep{C}{\lteval{\left(\alpha T\right)}{\vect{u}}}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrep{C}{\alpha\lteval{T}{\vect{u}}}&& \knowl{./knowl/definition-LTSM.html}{\text{Definition LTSM}}\\ &=\alpha\vectrep{C}{\lteval{T}{\vect{u}}}&& \knowl{./knowl/definition-LT.html}{\text{Definition LT}}\\ &=\alpha\left(\matrixrep{T}{B}{C}\vectrep{B}{\vect{u}}\right)&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\left(\alpha\matrixrep{T}{B}{C}\right)\vectrep{B}{\vect{u}}&& \knowl{./knowl/theorem-MMSMM.html}{\text{Theorem MMSMM}}\\ &=\left(\alpha\matrixrep{T}{B}{C}\right)\vect{x}&&\text{Substitution}\text{.} \end{align*}

Since the matrices $$\matrixrep{\alpha T}{B}{C}$$ and $$\alpha\matrixrep{T}{B}{C}$$ have equal matrix-vector products for every vector in $$\complex{n}\text{,}$$ by Theorem EMMVP they are equal matrices.

The vector space of all linear transformations from $$U$$ to $$V$$ is now isomorphic to the vector space of all $$m\times n$$ matrices.

Let $$\vect{x}$$ be any vector in $$\complex{n}\text{.}$$ Define $$\vect{u}\in U$$ by $$\vect{u}=\vectrepinv{B}{\vect{x}}\text{,}$$ so $$\vect{x}=\vectrep{B}{\vect{u}}\text{.}$$ Then,

\begin{align*} \matrixrep{\compose{S}{T}}{B}{D}\vect{x}&= \matrixrep{\compose{S}{T}}{B}{D}\vectrep{B}{\vect{u}}&&\text{Substitution}\\ &=\vectrep{D}{\lteval{\left(\compose{S}{T}\right)}{\vect{u}}}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrep{D}{\lteval{S}{\lteval{T}{\vect{u}}}}&& \knowl{./knowl/definition-LTC.html}{\text{Definition LTC}}\\ &=\matrixrep{S}{C}{D}\vectrep{C}{\lteval{T}{\vect{u}}}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\matrixrep{S}{C}{D}\left(\matrixrep{T}{B}{C}\vectrep{B}{\vect{u}}\right)&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\left(\matrixrep{S}{C}{D}\matrixrep{T}{B}{C}\right)\vectrep{B}{\vect{u}}&& \knowl{./knowl/theorem-MMA.html}{\text{Theorem MMA}}\\ &=\left(\matrixrep{S}{C}{D}\matrixrep{T}{B}{C}\right)\vect{x}&&\text{Substitution}\text{.} \end{align*}

Since the matrices $$\matrixrep{\compose{S}{T}}{B}{D}$$ and $$\matrixrep{S}{C}{D}\matrixrep{T}{B}{C}$$ have equal matrix-vector products for every vector in $$\complex{n}\text{,}$$ by Theorem EMMVP they are equal matrices.

This is the second great surprise of introductory linear algebra. Matrices are linear transformations (functions, really), and matrix multiplication is function composition! We can form the composition of two linear transformations, then form the matrix representation of the result. Or we can form the matrix representation of each linear transformation separately, then multiply the two representations together via Definition MM. In either case, we arrive at the same result.

Consider the two linear transformations,

\begin{align*} \ltdefn{T}{\complex{2}}{P_2}&\quad \lteval{T}{\colvector{a\\b}}=(-a + 3b)+(2a + 4b)x+(a - 2b)x^2\\ \ltdefn{S}{P_2}{M_{22}}&\quad\ \lteval{S}{a+bx+cx^2}= \begin{bmatrix} 2a + b + 2c & a + 4b - c\\ -a + 3c & 3a + b + 2 c\end{bmatrix} \end{align*}

and bases for $$\complex{2}\text{,}$$ $$P_2$$ and $$M_{22}$$ (respectively),

\begin{align*} B&=\set{\colvector{3\\1},\,\colvector{2\\1}}\\ C&=\set{1-2x+x^2,\,-1+3x,\,2x+3x^2}\\ D&=\set{ \begin{bmatrix} 1 & -2\\ 1 & -1 \end{bmatrix},\, \begin{bmatrix} 1 & -1\\ 1 & 2\end{bmatrix},\, \begin{bmatrix} -1 & 2 \\ 0 & 0 \end{bmatrix},\, \begin{bmatrix} 2 & -3 \\ 2 & 2\end{bmatrix} }\text{.} \end{align*}

Begin by computing the new linear transformation that is the composition of $$T$$ and $$S$$ (Definition LTC, Theorem CLTLT), $$\ltdefn{\left(\compose{S}{T}\right)}{\complex{2}}{M_{22}}\text{,}$$

\begin{align*} &\lteval{\left(\compose{S}{T}\right)}{\colvector{a\\b}}=\lteval{S}{\lteval{T}{\colvector{a\\b}}}\\ &\quad\quad=\lteval{S}{(-a + 3b)+(2a + 4b)x+(a - 2b)x^2}\\ &\quad\quad= \begin{bmatrix} 2(-a + 3b) + (2a + 4b) + 2(a - 2b) & (-a + 3b) + 4(2a + 4b) - (a - 2b)\\ -(-a + 3b) + 3(a - 2b) & 3(-a + 3b) + (2a + 4b) + 2(a - 2b) \end{bmatrix}\\ &\quad\quad= \begin{bmatrix} 2a + 6b & 6a + 21b\\ 4a - 9b &a + 9b \end{bmatrix}\text{.} \end{align*}

Now compute the matrix representations (Definition MR) for each of these three linear transformations ($$T\text{,}$$ $$S\text{,}$$ $$\compose{S}{T}$$), relative to the appropriate bases. First for $$T\text{,}$$

\begin{align*} \vectrep{C}{\lteval{T}{\colvector{3\\1}}}&=\vectrep{C}{10x+x^2}\\ &=\vectrep{C}{28(1-2x+x^2)+28(-1+3x)+(-9)(2x+3x^2)}\\ &=\colvector{28\\28\\-9}\\ \vectrep{C}{\lteval{T}{\colvector{2\\1}}}&=\vectrep{C}{1+8x}\\ &=\vectrep{C}{33(1-2x+x^2)+32(-1+3x)+(-11)(2x+3x^2)}\\ &=\colvector{33\\32\\-11}\text{.} \end{align*}

So we have the matrix representation of $$T\text{,}$$

\begin{equation*} \matrixrep{T}{B}{C}= \begin{bmatrix} 28 & 33 \\ 28 & 32 \\ -9 & -11 \end{bmatrix}\text{.} \end{equation*}

Now, a representation of $$S\text{,}$$

\begin{align*} &\vectrep{D}{\lteval{S}{1-2x+x^2}}=\vectrep{D}{\begin{bmatrix}2 & -8 \\ 2 & 3 \end{bmatrix}}\\ &\quad\quad=\vectrep{D}{ (-11)\begin{bmatrix} 1 & -2\\ 1 & -1 \end{bmatrix}+ (-21)\begin{bmatrix} 1 & -1\\ 1 & 2\end{bmatrix}+ 0\begin{bmatrix} -1 & 2 \\ 0 & 0 \end{bmatrix}+ (17)\begin{bmatrix} 2 & -3 \\ 2 & 2\end{bmatrix} }\\ &\quad\quad=\colvector{-11\\-21\\0\\17}\\ &\vectrep{D}{\lteval{S}{-1+3x}}=\vectrep{D}{\begin{bmatrix}1 & 11 \\ 1 & 0\end{bmatrix}}\\ &\quad\quad=\vectrep{D}{ 26\begin{bmatrix} 1 & -2\\ 1 & -1 \end{bmatrix}+ 51\begin{bmatrix} 1 & -1\\ 1 & 2\end{bmatrix}+ 0\begin{bmatrix} -1 & 2 \\ 0 & 0 \end{bmatrix}+ (-38)\begin{bmatrix} 2 & -3 \\ 2 & 2\end{bmatrix} }\\ &\quad\quad=\colvector{26\\51\\0\\-38}\\ &\vectrep{D}{\lteval{S}{2x+3x^2}}=\vectrep{D}{\begin{bmatrix}8 & 5 \\ 9 & 8\end{bmatrix}}\\ &\quad\quad=\vectrep{D}{ 34\begin{bmatrix} 1 & -2\\ 1 & -1 \end{bmatrix}+ 67\begin{bmatrix} 1 & -1\\ 1 & 2\end{bmatrix}+ 1\begin{bmatrix} -1 & 2 \\ 0 & 0 \end{bmatrix}+ (-46)\begin{bmatrix} 2 & -3 \\ 2 & 2\end{bmatrix} }\\ &\quad\quad=\colvector{34\\67\\1\\-46}\text{.} \end{align*}

So we have the matrix representation of $$S\text{,}$$

\begin{equation*} \matrixrep{S}{C}{D}= \begin{bmatrix} -11 & 26 & 34 \\ -21 & 51 & 67 \\ 0 & 0 & 1 \\ 17 & -38 & -46 \end{bmatrix}\text{.} \end{equation*}

Finally, a representation of $$\compose{S}{T}\text{,}$$

\begin{align*} &\vectrep{D}{\lteval{\left(\compose{S}{T}\right)}{\colvector{3\\1}}}=\vectrep{D}{\begin{bmatrix}12 & 39\\3 &12\end{bmatrix}}\\ &\quad\quad=\vectrep{D}{ 114\begin{bmatrix} 1 & -2\\ 1 & -1 \end{bmatrix}+ 237\begin{bmatrix} 1 & -1\\ 1 & 2\end{bmatrix}+ (-9)\begin{bmatrix} -1 & 2 \\ 0 & 0 \end{bmatrix}+ (-174)\begin{bmatrix} 2 & -3 \\ 2 & 2\end{bmatrix} }\\ &\quad\quad=\colvector{114\\237\\-9\\-174}\\ &\vectrep{D}{\lteval{\left(\compose{S}{T}\right)}{\colvector{2\\1}}}=\vectrep{D}{\begin{bmatrix}10 & 33\\-1 & 11\end{bmatrix}}\\ &\quad\quad=\vectrep{D}{ 95\begin{bmatrix} 1 & -2\\ 1 & -1 \end{bmatrix}+ 202\begin{bmatrix} 1 & -1\\ 1 & 2\end{bmatrix}+ (-11)\begin{bmatrix} -1 & 2 \\ 0 & 0 \end{bmatrix}+ (-149)\begin{bmatrix} 2 & -3 \\ 2 & 2\end{bmatrix} }\\ &\quad\quad=\colvector{95\\202\\-11\\-149}\text{.} \end{align*}

So we have the matrix representation of $$\compose{S}{T}\text{,}$$

\begin{equation*} \matrixrep{\compose{S}{T}}{B}{D}= \begin{bmatrix} 114 & 95 \\ 237 & 202 \\ -9 & -11 \\ -174 & -149 \end{bmatrix}\text{.} \end{equation*}

Now, we are all set to verify the conclusion of Theorem MRCLT,

\begin{align*} \matrixrep{S}{C}{D}\matrixrep{T}{B}{C} &= \begin{bmatrix} -11 & 26 & 34 \\ -21 & 51 & 67 \\ 0 & 0 & 1 \\ 17 & -38 & -46 \end{bmatrix} \begin{bmatrix} 28 & 33 \\ 28 & 32 \\ -9 & -11 \end{bmatrix}\\ &= \begin{bmatrix} 114 & 95 \\ 237 & 202 \\ -9 & -11 \\ -174 & -149 \end{bmatrix}\\ &=\matrixrep{\compose{S}{T}}{B}{D}\text{.} \end{align*}

We have intentionally used nonstandard bases. If you were to choose “nice” bases for the three vector spaces, then the result of the theorem might be rather transparent. But this would still be a worthwhile exercise — give it a go.

A diagram, similar to ones we have seen earlier, might make the importance of this theorem clearer.

One of our goals in the first part of this book is to make the definition of matrix multiplication (Definition MVP, Definition MM) seem as natural as possible. However, many of us are brought up with an entry-by-entry description of matrix multiplication (Theorem EMP) as the definition of matrix multiplication, and then theorems about columns of matrices and linear combinations follow from that definition. With this unmotivated definition, the realization that matrix multiplication is function composition is quite remarkable. It is an interesting exercise to begin with the question, “What is the matrix representation of the composition of two linear transformations?” and then, without using any theorems about matrix multiplication, finally arrive at the entry-by-entry description of matrix multiplication. Try it yourself (Exercise MR.T80).

#### SageSUTH3.Sage Under The Hood, Round 3.

We have seen that Sage is able to add two linear transformations, or multiply a single linear transformation by a scalar. Under the hood, this is accomplished by simply adding matrix representations, or multiplying a matrix by a scalar, according to Theorem MRSLT and Theorem MRMLT (respectively). Theorem MRCLT says linear transformation composition is matrix representation multiplication. Is it still a mystery why we use the symbol *for linear transformation composition in Sage?

We could do several examples here, but you should now be able to construct them yourselves. We will do just one, linear transformation composition is matrix representation multiplication.

x1, x2, x3, x4 = var('x1, x2, x3, x4')
outputs = [-5*x1 - 2*x2 +   x3,
4*x1 - 3*x2 - 3*x3,
4*x1 - 6*x2 - 4*x3,
5*x1 + 3*x2       ]
T_symbolic(x1, x2, x3) = outputs
outputs = [-3*x1 - x2 + x3 + 2*x4,
7*x1 + x2 + x3 -   x4]
S_symbolic(x1, x2, x3, x4) = outputs
b0 = vector(QQ, [-1, -2,  2])
b1 = vector(QQ, [ 1,  1,  0])
b2 = vector(QQ, [ 0,  3, -5])
U = (QQ^3).subspace_with_basis([b0, b1, b2])
c0 = vector(QQ, [ 0,  0,  2,  1])
c1 = vector(QQ, [ 2, -3, -1, -6])
c2 = vector(QQ, [-2,  3,  2,  7])
c3 = vector(QQ, [ 1, -2, -4, -6])
V = (QQ^4).subspace_with_basis([c0, c1, c2, c3])
d0 = vector(QQ, [3, 4])
d1 = vector(QQ, [2, 3])
W = (QQ^2).subspace_with_basis([d0, d1])
T = linear_transformation(U, V, T_symbolic)
S = linear_transformation(V, W, S_symbolic)
(S*T).matrix('right')

S.matrix(side='right')*T.matrix(side='right')


Extra credit: what changes do you need to make if you dropped the side='right'option on these three matrix representations?

### SubsectionPMRProperties of Matrix Representations

It will not be a surprise to discover that the kernel and range of a linear transformation are closely related to the null space and column space of the transformation's matrix representation. Perhaps this idea has been bouncing around in your head already, even before seeing the definition of a matrix representation. However, with a formal definition of a matrix representation (Definition MR), and a fundamental theorem to go with it (Theorem FTMR) we can be formal about the relationship, using the idea of isomorphic vector spaces (Definition IVS). Here are the twin theorems.

To establish that two vector spaces are isomorphic, we must find an isomorphism between them, an invertible linear transformation (Definition IVS). The kernel of the linear transformation $$T\text{,}$$ $$\krn{T}\text{,}$$ is a subspace of $$U\text{,}$$ while the null space of the matrix representation, $$\nsp{\matrixrep{T}{B}{C}}$$ is a subspace of $$\complex{n}\text{.}$$ The function $$\vectrepname{B}$$ is defined as a function from $$U$$ to $$\complex{n}\text{,}$$ but we can just as well employ the definition of $$\vectrepname{B}$$ as a function from $$\krn{T}$$ to $$\nsp{\matrixrep{T}{B}{C}}\text{.}$$

We must first insure that if we choose an input for $$\vectrepname{B}$$ from $$\krn{T}$$ that then the output will be an element of $$\nsp{\matrixrep{T}{B}{C}}\text{.}$$ So suppose that $$\vect{u}\in\krn{T}\text{.}$$ Then

\begin{align*} \matrixrep{T}{B}{C}\vectrep{B}{\vect{u}} &=\vectrep{C}{\lteval{T}{\vect{u}}}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrep{C}{\zerovector}&& \knowl{./knowl/definition-KLT.html}{\text{Definition KLT}}\\ &=\zerovector&& \knowl{./knowl/theorem-LTTZZ.html}{\text{Theorem LTTZZ}}\text{.} \end{align*}

This says that $$\vectrep{B}{\vect{u}}\in\nsp{\matrixrep{T}{B}{C}}\text{,}$$ as desired.

The restriction in the size of the domain and codomain $$\vectrepname{B}$$ will not affect the fact that $$\vectrepname{B}$$ is a linear transformation (Theorem VRLT), nor will it affect the fact that $$\vectrepname{B}$$ is injective (Theorem VRI). Something must be done though to verify that $$\vectrepname{B}$$ is surjective. To this end, appeal to the definition of surjective (Definition SLT), and suppose that we have an element of the codomain, $$\vect{x}\in\nsp{\matrixrep{T}{B}{C}}\subseteq\complex{n}$$ and we wish to find an element of the domain with $$\vect{x}$$ as its image. We now show that the desired element of the domain is $$\vect{u}=\vectrepinv{B}{\vect{x}}\text{.}$$ First, verify that $$\vect{u}\in\krn{T}\text{,}$$

\begin{align*} \lteval{T}{\vect{u}} &=\lteval{T}{\vectrepinv{B}{\vect{x}}}\\ &=\vectrepinv{C}{\matrixrep{T}{B}{C}\left(\vectrep{B}{\vectrepinv{B}{\vect{x}}}\right)}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrepinv{C}{\matrixrep{T}{B}{C}\left(\lteval{I_{\complex{n}}}{\vect{x}}\right)}&& \knowl{./knowl/definition-IVLT.html}{\text{Definition IVLT}}\\ &=\vectrepinv{C}{\matrixrep{T}{B}{C}\vect{x}}&& \knowl{./knowl/definition-IDLT.html}{\text{Definition IDLT}}\\ &=\vectrepinv{C}{\zerovector_{\complex{n}}}&& \knowl{./knowl/definition-KLT.html}{\text{Definition KLT}}\\ &=\zerovector_{V}&& \knowl{./knowl/theorem-LTTZZ.html}{\text{Theorem LTTZZ}}\\ \end{align*}

Second, verify that the proposed isomorphism, $$\vectrepname{B}\text{,}$$ takes $$\vect{u}$$ to $$\vect{x}\text{,}$$

\begin{align*} \vectrep{B}{\vect{u}}&=\vectrep{B}{\vectrepinv{B}{\vect{x}}}&&\text{Substitution}\\ &=\lteval{I_{\complex{n}}}{\vect{x}}&& \knowl{./knowl/definition-IVLT.html}{\text{Definition IVLT}}\\ &=\vect{x}&& \knowl{./knowl/definition-IDLT.html}{\text{Definition IDLT}}\text{.} \end{align*}

With $$\vectrepname{B}$$ demonstrated to be an injective and surjective linear transformation from $$\krn{T}$$ to $$\nsp{\matrixrep{T}{B}{C}}\text{,}$$ Theorem ILTIS tells us $$\vectrepname{B}$$ is invertible, and so by Definition IVS, we say $$\krn{T}$$ and $$\nsp{\matrixrep{T}{B}{C}}$$ are isomorphic.

Consider the kernel of the linear transformation, $$\ltdefn{T}{M_{22}}{P_2}\text{,}$$ given by

\begin{equation*} \lteval{T}{\begin{bmatrix}a&b\\c&d\end{bmatrix}}= (2a-b+c-5d)+(a+4b+5c+2d)x+(3a-2b+c-8d)x^2\text{.} \end{equation*}

We will begin with a matrix representation of $$T$$ relative to the bases for $$M_{22}$$ and $$P_2$$ (respectively),

\begin{align*} B&=\set{ \begin{bmatrix}1 & 2 \\ -1 & -1\end{bmatrix},\, \begin{bmatrix}1 & 3 \\ -1 & -4\end{bmatrix},\, \begin{bmatrix}1 & 2 \\ 0 & -2\end{bmatrix},\, \begin{bmatrix}2 & 5 \\ -2 & -4\end{bmatrix} }\\ C&=\set{1+x+x^2,\,2+3x,\,-1-2x^2}\text{.} \end{align*}

Then,

\begin{align*} \vectrep{C}{\lteval{T}{\begin{bmatrix}1 & 2 \\ -1 & -1\end{bmatrix}}} &=\vectrep{C}{4+2x+6x^2}\\ &=\vectrep{C}{2(1+x+x^2)+0(2+3x)+(-2)(-1-2x^2)}\\ &=\colvector{2\\0\\-2}\\ \vectrep{C}{\lteval{T}{\begin{bmatrix}1 & 3 \\ -1 & -4\end{bmatrix}}} &=\vectrep{C}{18+28x^2}\\ &=\vectrep{C}{(-24)(1+x+x^2)+8(2+3x)+(-26)(-1-2x^2)}\\ &=\colvector{-24\\8\\-26}\\ \vectrep{C}{\lteval{T}{\begin{bmatrix}1 & 2 \\ 0 & -2\end{bmatrix}}} &=\vectrep{C}{10+5x+15x^2}\\ &=\vectrep{C}{5(1+x+x^2)+0(2+3x)+(-5)(-1-2x^2)}\\ &=\colvector{5\\0\\-5}\\ \vectrep{C}{\lteval{T}{\begin{bmatrix}2 & 5 \\ -2 & -4\end{bmatrix}}} &=\vectrep{C}{17+4x+26x^2}\\ &=\vectrep{C}{(-8)(1+x+x^2)+(4)(2+3x)+(-17)(-1-2x^2)}\\ &=\colvector{-8\\4\\-17}\text{.} \end{align*}

So the matrix representation of $$T$$ (relative to $$B$$ and $$C$$) is

\begin{equation*} \matrixrep{T}{B}{C} = \begin{bmatrix} 2 & -24 & 5 & -8 \\ 0 & 8 & 0 & 4 \\ -2 & -26 & -5 & -17 \end{bmatrix}\text{.} \end{equation*}

We know from Theorem KNSI that the kernel of the linear transformation $$T$$ is isomorphic to the null space of the matrix representation $$\matrixrep{T}{B}{C}$$ and by studying the proof of Theorem KNSI we learn that $$\vectrepname{B}$$ is an isomorphism between these null spaces. Rather than trying to compute the kernel of $$T$$ using definitions and techniques from Chapter LT we will instead analyze the null space of $$\matrixrep{T}{B}{C}$$ using techniques from way back in Chapter V. First row-reduce $$\matrixrep{T}{B}{C}\text{,}$$

\begin{equation*} \begin{bmatrix} 2 & -24 & 5 & -8 \\ 0 & 8 & 0 & 4 \\ -2 & -26 & -5 & -17 \end{bmatrix} \rref \begin{bmatrix} \leading{1} & 0 & \frac{5}{2} & 2 \\ 0 & \leading{1} & 0 & \frac{1}{2} \\ 0 & 0 & 0 & 0 \end{bmatrix}\text{.} \end{equation*}

So, by Theorem BNS, a basis for $$\nsp{\matrixrep{T}{B}{C}}$$ is

\begin{equation*} \set{\colvector{-\frac{5}{2}\\0\\1\\0},\,\colvector{-2\\-\frac{1}{2}\\0\\1}}\text{.} \end{equation*}

We can now convert this basis of $$\nsp{\matrixrep{T}{B}{C}}$$ into a basis of $$\krn{T}$$ by applying $$\vectrepinvname{B}$$ to each element of the basis,

\begin{align*} \vectrepinv{B}{\colvector{-\frac{5}{2}\\0\\1\\0}}&= (-\frac{5}{2})\begin{bmatrix}1 & 2 \\ -1 & -1\end{bmatrix}+ 0\begin{bmatrix}1 & 3 \\ -1 & -4\end{bmatrix}+ 1\begin{bmatrix}1 & 2 \\ 0 & -2\end{bmatrix}+ 0\begin{bmatrix}2 & 5 \\ -2 & -4\end{bmatrix}\\ &=\begin{bmatrix}-\frac{3}{2} & -3 \\ \frac{5}{2} & \frac{1}{2}\end{bmatrix}\\ \vectrepinv{B}{\colvector{-2\\-\frac{1}{2}\\0\\1}}&= (-2)\begin{bmatrix}1 & 2 \\ -1 & -1\end{bmatrix}+ (-\frac{1}{2})\begin{bmatrix}1 & 3 \\ -1 & -4\end{bmatrix}+ 0\begin{bmatrix}1 & 2 \\ 0 & -2\end{bmatrix}+ 1\begin{bmatrix}2 & 5 \\ -2 & -4\end{bmatrix}\\ &=\begin{bmatrix}-\frac{1}{2} & -\frac{1}{2} \\ \frac{1}{2} &0\end{bmatrix}\text{.} \end{align*}

So the set

\begin{equation*} \set{ \begin{bmatrix}-\frac{3}{2} & -3 \\ \frac{5}{2} & \frac{1}{2}\end{bmatrix},\, \begin{bmatrix}-\frac{1}{2} & -\frac{1}{2} \\ \frac{1}{2} &0\end{bmatrix} } \end{equation*}

is a basis for $$\krn{T}\text{.}$$ Just for fun, you might evaluate $$T$$ with each of these two basis vectors and verify that the output is the zero polynomial (Exercise MR.C10).

An entirely similar result applies to the range of a linear transformation and the column space of a matrix representation of the linear transformation.

To establish that two vector spaces are isomorphic, we must find an isomorphism between them, an invertible linear transformation (Definition IVS). The range of the linear transformation $$T\text{,}$$ $$\rng{T}\text{,}$$ is a subspace of $$V\text{,}$$ while the column space of the matrix representation, $$\csp{\matrixrep{T}{B}{C}}$$ is a subspace of $$\complex{m}\text{.}$$ The function $$\vectrepname{C}$$ is defined as a function from $$V$$ to $$\complex{m}\text{,}$$ but we can just as well employ the definition of $$\vectrepname{C}$$ as a function from $$\rng{T}$$ to $$\csp{\matrixrep{T}{B}{C}}\text{.}$$

We must first insure that if we choose an input for $$\vectrepname{C}$$ from $$\rng{T}$$ that then the output will be an element of $$\csp{\matrixrep{T}{B}{C}}\text{.}$$ So suppose that $$\vect{v}\in\rng{T}\text{.}$$ Then there is a vector $$\vect{u}\in U\text{,}$$ such that $$\lteval{T}{\vect{u}}=\vect{v}\text{.}$$ Consider

\begin{align*} \matrixrep{T}{B}{C}\vectrep{B}{\vect{u}} &=\vectrep{C}{\lteval{T}{\vect{u}}}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrep{C}{\vect{v}}&& \knowl{./knowl/definition-RLT.html}{\text{Definition RLT}}\text{.} \end{align*}

This says that $$\vectrep{C}{\vect{v}}\in\csp{\matrixrep{T}{B}{C}}\text{,}$$ as desired.

The restriction in the size of the domain and codomain will not affect the fact that $$\vectrepname{C}$$ is a linear transformation (Theorem VRLT), nor will it affect the fact that $$\vectrepname{C}$$ is injective (Theorem VRI). Something must be done though to verify that $$\vectrepname{C}$$ is surjective. This all gets a bit confusing, since the domain of our isomorphism is the range of the linear transformation, so think about your objects as you go. To establish that $$\vectrepname{C}$$ is surjective, appeal to the definition of a surjective linear transformation (Definition SLT), and suppose that we have an element of the codomain, $$\vect{y}\in\csp{\matrixrep{T}{B}{C}}\subseteq\complex{m}$$ and we wish to find an element of the domain with $$\vect{y}$$ as its image. Since $$\vect{y}\in\csp{\matrixrep{T}{B}{C}}\text{,}$$ there exists a vector, $$\vect{x}\in\complex{n}$$ with $$\matrixrep{T}{B}{C}\vect{x}=\vect{y}\text{.}$$

We now show that the desired element of the domain is $$\vect{v}=\vectrepinv{C}{\vect{y}}\text{.}$$ First, verify that $$\vect{v}\in\rng{T}$$ by applying $$T$$ to $$\vect{u}=\vectrepinv{B}{\vect{x}}\text{,}$$

\begin{align*} \lteval{T}{\vect{u}}&=\lteval{T}{\vectrepinv{B}{\vect{x}}}\\ &=\vectrepinv{C}{\matrixrep{T}{B}{C}\left(\vectrep{B}{\vectrepinv{B}{\vect{x}}}\right)}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrepinv{C}{\matrixrep{T}{B}{C}\left(\lteval{I_{\complex{n}}}{\vect{x}}\right)}&& \knowl{./knowl/definition-IVLT.html}{\text{Definition IVLT}}\\ &=\vectrepinv{C}{\matrixrep{T}{B}{C}\vect{x}}&& \knowl{./knowl/definition-IDLT.html}{\text{Definition IDLT}}\\ &=\vectrepinv{C}{\vect{y}}&& \knowl{./knowl/definition-CSM.html}{\text{Definition CSM}}\\ &=\vect{v}&&\text{Substitution}\\ \end{align*}

Second, verify that the proposed isomorphism, $$\vectrepname{C}\text{,}$$ takes $$\vect{v}$$ to $$\vect{y}\text{,}$$

\begin{align*} \vectrep{C}{\vect{v}}&=\vectrep{C}{\vectrepinv{C}{\vect{y}}}&&\text{Substitution}\\ &=\lteval{I_{\complex{m}}}{\vect{y}}&& \knowl{./knowl/definition-IVLT.html}{\text{Definition IVLT}}\\ &=\vect{y}&& \knowl{./knowl/definition-IDLT.html}{\text{Definition IDLT}}\text{.} \end{align*}

With $$\vectrepname{C}$$ demonstrated to be an injective and surjective linear transformation from $$\rng{T}$$ to $$\csp{\matrixrep{T}{B}{C}}\text{,}$$ Theorem ILTIS tells us $$\vectrepname{C}$$ is invertible, and so by Definition IVS, we say $$\rng{T}$$ and $$\csp{\matrixrep{T}{B}{C}}$$ are isomorphic.

In this example, we will recycle the linear transformation $$T$$ and the bases $$B$$ and $$C$$ of Example KVMR but now we will compute the range of $$\ltdefn{T}{M_{22}}{P_2}\text{,}$$ given by

\begin{equation*} \lteval{T}{\begin{bmatrix}a&b\\c&d\end{bmatrix}}= (2a-b+c-5d)+(a+4b+5c+2d)x+(3a-2b+c-8d)x^2\text{.} \end{equation*}

With bases $$B$$ and $$C\text{,}$$

\begin{align*} B&=\set{ \begin{bmatrix}1 & 2 \\ -1 & -1\end{bmatrix},\, \begin{bmatrix}1 & 3 \\ -1 & -4\end{bmatrix},\, \begin{bmatrix}1 & 2 \\ 0 & -2\end{bmatrix},\, \begin{bmatrix}2 & 5 \\ -2 & -4\end{bmatrix} }\\ C&=\set{1+x+x^2,\,2+3x,\,-1-2x^2} \end{align*}

we obtain the matrix representation

\begin{equation*} \matrixrep{T}{B}{C} = \begin{bmatrix} 2 & -24 & 5 & -8 \\ 0 & 8 & 0 & 4 \\ -2 & -26 & -5 & -17 \end{bmatrix}\text{.} \end{equation*}

We know from Theorem RCSI that the range of the linear transformation $$T$$ is isomorphic to the column space of the matrix representation $$\matrixrep{T}{B}{C}$$ and by studying the proof of Theorem RCSI we learn that $$\vectrepname{C}$$ is an isomorphism between these subspaces. Notice that since the range is a subspace of the codomain, we will employ $$\vectrepname{C}$$ as the isomorphism, rather than $$\vectrepname{B}\text{,}$$ which was the correct choice for an isomorphism between the null spaces of Example KVMR.

Rather than trying to compute the range of $$T$$ using definitions and techniques from Chapter LT we will instead analyze the column space of $$\matrixrep{T}{B}{C}$$ using techniques from way back in Chapter M. First row-reduce $$\transpose{\left(\matrixrep{T}{B}{C}\right)}\text{,}$$

\begin{equation*} \begin{bmatrix} 2 & 0 & -2 \\ -24 & 8 & -26 \\ 5 & 0 & -5 \\ -8 & 4 & -17 \end{bmatrix} \rref \begin{bmatrix} \leading{1} & 0 & -1 \\ 0 & \leading{1} & -\frac{25}{4} \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix}\text{.} \end{equation*}

Now employ Theorem CSRST and Theorem BRS (there are other methods we could choose here to compute the column space, such as Theorem BCS) to obtain the basis for $$\csp{\matrixrep{T}{B}{C}}\text{,}$$

\begin{equation*} \set{ \colvector{1\\0\\-1},\, \colvector{0\\1\\-\frac{25}{4}} }\text{.} \end{equation*}

We can now convert this basis of $$\csp{\matrixrep{T}{B}{C}}$$ into a basis of $$\rng{T}$$ by applying $$\vectrepinvname{C}$$ to each element of the basis,

\begin{align*} \vectrepinv{C}{\colvector{1\\0\\-1}}&= (1+x+x^2)-(-1-2x^2)=2+x+3x^2\\ \vectrepinv{C}{\colvector{0\\1\\-\frac{25}{4}}}&= (2+3x)-\frac{25}{4}(-1-2x^2)=\frac{33}{4}+3x+\frac{31}{2}x^2\text{.} \end{align*}

So the set

\begin{equation*} \set{ 2+3x+3x^2,\, \frac{33}{4}+3x+\frac{31}{2}x^2 } \end{equation*}

is a basis for $$\rng{T}\text{.}$$

Theorem KNSI and Theorem RCSI can be viewed as further formal evidence for the The Coordinatization Principle, though they are not direct consequences.

Figure KRI is meant to suggest Theorem KNSI and Theorem RCSI, in addition to their proofs (and so carry the same notation as the statements of these two theorems). The dashed lines indicate a subspace relationship, with the smaller vector space lower down in the diagram. The central square is highly reminiscent of Figure FTMR. Each of the four vector representations is an isomorphism, so the inverse linear transformation could be depicted with an arrow pointing in the other direction. The four vector spaces across the bottom are familiar from the earliest days of the course, while the four vector spaces across the top are completely abstract. The vector representations that are restrictions (far left and far right) are the functions shown to be invertible representations as the key technique in the proofs of Theorem KNSI and Theorem RCSI. So this diagram could be helpful as you study those two proofs.

#### SageLTR.Linear Transformation Restrictions.

Theorem KNSI and Theorem RCSI have two of the most subtle proofs we have seen so far. The conclusion that two vector spaces are isomorphic is established by actually constructing an isomorphism between the vector spaces. To build the isomorphism, we begin with a familar object, a vector representation linear transformation, but the hard work is showing that we can “restrict” the domain and codomain of this function and still arrive at a legitimate (invertible) linear transformation. In an effort to make the proofs more concrete, we will walk through a nontrivial example for Theorem KNSI, and you might try to do the same for Theorem RCSI. (An understanding of this subsection is not needed for the remainder — its second purpose is to demonstrate some of the powerful tools Sage provides.)

Here are the pieces. We build a linear transformation with two different representations, one with respect to standard bases, the other with respect to less-obvious bases.

x1, x2, x3, x4, x5 = var('x1, x2, x3, x4, x5')
outputs = [ x1 - x2 - 5*x3 +   x4 +   x5,
x1      - 2*x3 -   x4 -   x5,
- x2 - 3*x3 + 2*x4 + 2*x5,
-x1 + x2 + 5*x3 -   x4 -   x5]
T_symbolic(x1, x2, x3, x4, x5) = outputs
b0 = vector(QQ, [-1,  6,  5,  5,  1])
b1 = vector(QQ, [-1,  5,  4,  4,  1])
b2 = vector(QQ, [-2,  4,  3,  2,  5])
b3 = vector(QQ, [ 1, -1,  0,  1, -5])
b4 = vector(QQ, [ 3, -7, -6, -5, -4])
U = (QQ^5).subspace_with_basis([b0, b1, b2, b3, b4])
c0 = vector(QQ, [1, 1, 1, -3])
c1 = vector(QQ, [-2, 3, -6, -7])
c2 = vector(QQ, [0, -1, 1, 2])
c3 = vector(QQ, [-1, 3, -4, -7])
V = (QQ^4).subspace_with_basis([c0, c1, c2, c3])
T_plain = linear_transformation(QQ^5, QQ^4, T_symbolic)
T_fancy = linear_transformation(   U,    V, T_symbolic)


Now we compute the kernel of the linear transformation using the “plain” version, and the null space of a matrix representation coming from the “fancy” version.

K = T_plain.kernel()
K

MK = T_fancy.matrix(side='right').right_kernel()
MK


So quite obviously, the kernel of the linear transformation is quite different looking from the null space of the matrix representation. Though it is no accident that they have the same dimension. Now we build the necessary vector representation, and use two Sage commands to “restrict” the function to a smaller domain (the kernel of the linear transformation) and a smaller codomain (the null space of the matrix representation relative to nonstandard bases).

rhoB = linear_transformation(U, QQ^5, (QQ^5).basis())
rho_restrict = rhoB.restrict_domain(K).restrict_codomain(MK)
rho_restrict


The first success is that the restriction was even created. Sage would recognize if the original linear transformation ever carried an input from the restricted domain to an output that was not contained in the proposed codomain, and would have raised an error in that event. Phew! Guaranteeing this success was the first big step in the proof of Theorem KNSI. Notice that the matrix representation of the restriction is a $$3\times 3$$ matrix, since the restriction runs between a domain and codomain that each have dimension 3. These two vector spaces (the domain and codomain of the restriction) have dimension 3 but still contain vectors with 5 entries in their un-coordinatized versions.

The next two steps of the proof show that the restriction is injective (easy in the proof) and surjective (hard in the proof). In Sage, here is the second success,

rho_restrict.is_injective()

rho_restrict.is_surjective()


Verified as invertible, rho_restrictqualifies as an isomorphism between the linear transformation kernel, K, and the matrix representation null space, MK. Only an example, but still very nice. Your turn — can you create a verfication of Theorem RCSI (for this example, or some other nontrivial example you might create yourself)?

### SubsectionIVLTInvertible Linear Transformations

We have seen, both in theorems and in examples, that questions about linear transformations are often equivalent to questions about matrices. It is the matrix representation of a linear transformation that makes this idea precise. Here is our final theorem that solidifies this connection.

##### (⇒)

Suppose $$T$$ is invertible, so the inverse linear transformation $$\ltdefn{\ltinverse{T}}{V}{U}$$ exists (Definition IVLT). Both linear transformations have matrix representations relative to the bases of $$U$$ and $$V\text{,}$$ namely $$\matrixrep{T}{B}{C}$$ and $$\matrixrep{\ltinverse{T}}{C}{B}$$ (Definition MR).

Then

\begin{align*} \matrixrep{\ltinverse{T}}{C}{B}\matrixrep{T}{B}{C} &=\matrixrep{\compose{\ltinverse{T}}{T}}{B}{B}&& \knowl{./knowl/theorem-MRCLT.html}{\text{Theorem MRCLT}}\\ &=\matrixrep{I_U}{B}{B}&& \knowl{./knowl/definition-IVLT.html}{\text{Definition IVLT}}\\ &=\left[ \left.\vectrep{B}{\lteval{I_U}{\vect{u}_1}}\right| \left.\vectrep{B}{\lteval{I_U}{\vect{u}_2}}\right| \ldots \left|\vectrep{B}{\lteval{I_U}{\vect{u}_n}}\right. \right]&& \knowl{./knowl/definition-MR.html}{\text{Definition MR}}\\ &=\left[ \left.\vectrep{B}{\vect{u}_1}\right| \left.\vectrep{B}{\vect{u}_2}\right| \left.\vectrep{B}{\vect{u}_3}\right| \ldots \left|\vectrep{B}{\vect{u}_n}\right. \right]&& \knowl{./knowl/definition-IDLT.html}{\text{Definition IDLT}}\\ &=\matrixcolumns{e}{n}&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\\ &=I_n&& \knowl{./knowl/definition-IM.html}{\text{Definition IM}}\text{.} \end{align*}

And

\begin{align*} \matrixrep{T}{B}{C}\matrixrep{\ltinverse{T}}{C}{B}&= \matrixrep{\compose{T}{\ltinverse{T}}}{C}{C}&& \knowl{./knowl/theorem-MRCLT.html}{\text{Theorem MRCLT}}\\ &=\matrixrep{I_V}{C}{C}&& \knowl{./knowl/definition-IVLT.html}{\text{Definition IVLT}}\\ &=\left[ \left.\vectrep{C}{\lteval{I_V}{\vect{v}_1}}\right| \left.\vectrep{C}{\lteval{I_V}{\vect{v}_2}}\right| \ldots \left|\vectrep{C}{\lteval{I_V}{\vect{v}_n}}\right. \right]&& \knowl{./knowl/definition-MR.html}{\text{Definition MR}}\\ &=\left[ \left.\vectrep{C}{\vect{v}_1}\right| \left.\vectrep{C}{\vect{v}_2}\right| \left.\vectrep{C}{\vect{v}_3}\right| \ldots \left|\vectrep{C}{\vect{v}_n}\right. \right]&& \knowl{./knowl/definition-IDLT.html}{\text{Definition IDLT}}\\ &=\matrixcolumns{e}{n}&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\\ &=I_n&& \knowl{./knowl/definition-IM.html}{\text{Definition IM}}\text{.} \end{align*}

These two equations show that $$\matrixrep{T}{B}{C}$$ and $$\matrixrep{\ltinverse{T}}{C}{B}$$ are inverse matrices (Definition MI) and establish that when $$T$$ is invertible, then $$\matrixrep{\ltinverse{T}}{C}{B}=\inverse{\left(\matrixrep{T}{B}{C}\right)}\text{.}$$

##### (⇐)

Suppose now that $$\matrixrep{T}{B}{C}$$ is an invertible matrix and hence nonsingular (Theorem NI). We compute the nullity of $$T\text{,}$$

\begin{align*} \nullity{T} &=\dimension{\krn{T}}&& \knowl{./knowl/definition-KLT.html}{\text{Definition KLT}}\\ &=\dimension{\nsp{\matrixrep{T}{B}{C}}}&& \knowl{./knowl/theorem-KNSI.html}{\text{Theorem KNSI}}\\ &=\nullity{\matrixrep{T}{B}{C}}&& \knowl{./knowl/definition-NOM.html}{\text{Definition NOM}}\\ &=0&& \knowl{./knowl/theorem-RNNM.html}{\text{Theorem RNNM}}\text{.} \end{align*}

So the kernel of $$T$$ is trivial, and by Theorem KILT, $$T$$ is injective.

We now compute the rank of $$T\text{,}$$

\begin{align*} \rank{T} &=\dimension{\rng{T}}&& \knowl{./knowl/definition-RLT.html}{\text{Definition RLT}}\\ &=\dimension{\csp{\matrixrep{T}{B}{C}}}&& \knowl{./knowl/theorem-RCSI.html}{\text{Theorem RCSI}}\\ &=\rank{\matrixrep{T}{B}{C}}&& \knowl{./knowl/definition-ROM.html}{\text{Definition ROM}}\\ &=\dimension{V}&& \knowl{./knowl/theorem-RNNM.html}{\text{Theorem RNNM}}\text{.} \end{align*}

Since the dimension of the range of $$T$$ equals the dimension of the codomain $$V\text{,}$$ by Theorem EDYES, $$\rng{T}=V\text{.}$$ Which says that $$T$$ is surjective by Theorem RSLT.

Because $$T$$ is both injective and surjective, by Theorem ILTIS, $$T$$ is invertible.

By now, the connections between matrices and linear transformations should be starting to become more transparent, and you may have already recognized the invertibility of a matrix as being tantamount to the invertibility of the associated matrix representation. The next example shows how to apply this theorem to the problem of actually building a formula for the inverse of an invertible linear transformation.

Consider the linear transformation

\begin{equation*} \ltdefn{R}{P_3}{M_{22}},\quad \lteval{R}{a+bx+cx^2+dx^3}= \begin{bmatrix} a +b - c + 2d & 2a + 3b - 2c + 3d\\ a + b + 2 d & -a + b + 2c - 5d \end{bmatrix}\text{.} \end{equation*}

If we wish to quickly find a formula for the inverse of $$R$$ (presuming it exists), then choosing “nice” bases will work best. So build a matrix representation of $$R$$ relative to the bases $$B$$ and $$C\text{,}$$

\begin{align*} B&=\set{1,\,x,\,x^2,\,x^3}\\ C&=\set{ \begin{bmatrix} 1 & 0 \\ 0 & 0 \end{bmatrix},\, \begin{bmatrix} 0 & 1 \\ 0 & 0 \end{bmatrix},\, \begin{bmatrix} 0 & 0 \\ 1 & 0 \end{bmatrix},\, \begin{bmatrix} 0 & 0 \\ 0 & 1 \end{bmatrix} }\text{.} \end{align*}

Then,

\begin{align*} \vectrep{C}{\lteval{R}{1}} &=\vectrep{C}{\begin{bmatrix} 1 & 2 \\ 1 & -1 \end{bmatrix}} =\colvector{1\\2\\1\\-1}\\ \vectrep{C}{\lteval{R}{x}} &=\vectrep{C}{\begin{bmatrix} 1 & 3 \\ 1 & 1 \end{bmatrix}} =\colvector{1\\3\\1\\1}\\ \vectrep{C}{\lteval{R}{x^2}} &=\vectrep{C}{\begin{bmatrix} -1 & -2 \\ 0 & 2 \end{bmatrix}} =\colvector{-1\\-2\\0\\2}\\ \vectrep{C}{\lteval{R}{x^3}} &=\vectrep{C}{\begin{bmatrix} 2 & 3 \\ 2 & -5 \end{bmatrix}} =\colvector{2\\3\\2\\-5}\text{.} \end{align*}

So a representation of $$R$$ is

\begin{equation*} \matrixrep{R}{B}{C}= \begin{bmatrix} 1 & 1 & -1 & 2 \\ 2 & 3 & -2 & 3 \\ 1 & 1 & 0 & 2 \\ -1 & 1 & 2 & -5 \end{bmatrix}\text{.} \end{equation*}

The matrix $$\matrixrep{R}{B}{C}$$ is invertible (as you can check) so we know for sure that $$R$$ is invertible by Theorem IMR. Furthermore,

\begin{equation*} \matrixrep{\ltinverse{R}}{C}{B} =\inverse{\left(\matrixrep{R}{B}{C}\right)} = \inverse{ \begin{bmatrix} 1 & 1 & -1 & 2 \\ 2 & 3 & -2 & 3 \\ 1 & 1 & 0 & 2 \\ -1 & 1 & 2 & -5 \end{bmatrix} } = \begin{bmatrix} 20 & -7 & -2 & 3 \\ -8 & 3 & 1 & -1 \\ -1 & 0 & 1 & 0 \\ -6 & 2 & 1 & -1 \end{bmatrix}\text{.} \end{equation*}

We can use this representation of the inverse linear transformation, in concert with Theorem FTMR, to determine an explicit formula for the inverse itself,

\begin{align*} \lteval{\ltinverse{R}}{\begin{bmatrix}a&b\\c&d\end{bmatrix}} &=\vectrepinv{B}{\matrixrep{\ltinverse{R}}{C}{B}\vectrep{C}{\begin{bmatrix}a&b\\c&d\end{bmatrix}}}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrepinv{B}{\inverse{\left(\matrixrep{R}{B}{C}\right)}\vectrep{C}{\begin{bmatrix}a&b\\c&d\end{bmatrix}}}&& \knowl{./knowl/theorem-IMR.html}{\text{Theorem IMR}}\\ &=\vectrepinv{B}{\inverse{\left(\matrixrep{R}{B}{C}\right)}\colvector{a\\b\\c\\d}}&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\\ &=\vectrepinv{B}{ \begin{bmatrix} 20 & -7 & -2 & 3 \\ -8 & 3 & 1 & -1 \\ -1 & 0 & 1 & 0 \\ -6 & 2 & 1 & -1 \end{bmatrix} \colvector{a\\b\\c\\d}}&& \knowl{./knowl/definition-MI.html}{\text{Definition MI}}\\ &=\vectrepinv{B}{ \colvector{ 20a - 7b - 2c + 3d\\ -8a + 3b + c -d\\ -a + c\\ -6a + 2b + c - d}}&& \knowl{./knowl/definition-MVP.html}{\text{Definition MVP}}\\ &=(20a - 7b - 2c + 3d)+(-8a + 3b + c -d)x\\ &\quad\quad +(-a + c)x^2+(-6a + 2b + c - d)x^3&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\text{.} \end{align*}

You might look back at Example AIVLT, where we first witnessed the inverse of a linear transformation and recognize that the inverse ($$S$$) was built from using the method of Example ILTVR with a matrix representation of $$T\text{.}$$

Choose bases $$B=C=\set{\vectorlist{e}{n}}$$ consisting of the standard unit vectors as a basis of $$\complex{n}$$ (Theorem SUVB) and build a matrix representation of $$T$$ relative to $$B$$ and $$C\text{.}$$ Then

\begin{align*} \vectrep{C}{\lteval{T}{\vect{e}_i}} &=\vectrep{C}{A\vect{e}_i}\\ &=\vectrep{C}{\vect{A}_i}\\ &=\vect{A}_i\text{.} \end{align*}

So then the matrix representation of $$T\text{,}$$ relative to $$B$$ and $$C\text{,}$$ is simply $$\matrixrep{T}{B}{C}=A\text{.}$$ With this observation, the proof becomes a specialization of Theorem IMR,

\begin{gather*} T\text{ is invertible} \iff\matrixrep{T}{B}{C}\text{ is invertible} \iff A\text{ is invertible}\text{.} \end{gather*}

This theorem may seem gratuitous. Why state such a special case of Theorem IMR? Because it adds another condition to our NMEx series of theorems, and in some ways it is the most fundamental expression of what it means for a matrix to be nonsingular — the associated linear transformation is invertible. This is our final update.

By Theorem IMILT, the new addition to this list is equivalent to the statement that $$A$$ is invertible, so we can expand Theorem NME8.

#### SageNME9.Nonsingular Matrix Equivalences, Round 9.

Our final fact about nonsingular matrices expresses the correspondence between invertible matrices and invertible linear transformations. As a Sage demonstration, we will begin with an invertible linear transformation and examine two matrix representations. We will create the linear transformation with nonstandard bases and then compute its representation relative to standard bases.

x1, x2, x3, x4 = var('x1, x2, x3, x4')
outputs = [ x1        - 2*x3 - 4*x4,
x2 -   x3 - 5*x4,
-x1 - 2*x2 + 2*x3 + 7*x4,
-   x2        +   x4]
T_symbolic(x1, x2, x3, x4) = outputs
b0 = vector(QQ, [ 1, -2, -1,  8])
b1 = vector(QQ, [ 0,  1,  0, -2])
b2 = vector(QQ, [-1, -2,  2, -5])
b3 = vector(QQ, [-1, -3,  2, -2])
U = (QQ^4).subspace_with_basis([b0, b1, b2, b3])
c0 = vector(QQ, [ 3, -1,  4, -8])
c1 = vector(QQ, [ 1,  0,  1, -1])
c2 = vector(QQ, [ 0,  2, -1,  6])
c3 = vector(QQ, [-1,  2, -2,  8])
V = (QQ^4).subspace_with_basis([c0, c1, c2, c3])
T = linear_transformation(U, V, T_symbolic)
T

T.is_invertible()

T.matrix(side='right').is_invertible()

(T^-1).matrix(side='right') == (T.matrix(side='right'))^-1


We can convert Tto a new representation using standard bases for QQ^4by computing images of the standard basis.

images = [T(u) for u in (QQ^4).basis()]
T_standard = linear_transformation(QQ^4, QQ^4, images)
T_standard

T_standard.matrix(side='right').is_invertible()


Understand that any matrix representation of T_symbolicwill have an invertible matrix representation, no matter which bases are used. If you look at the matrix representation of T_standardand the definition of T_symbolicthe construction of this example will be transparent, especially if you know the random matrix constructor,

A = random_matrix(QQ, 4, algorithm='unimodular', upper_bound=9)
A                                                      # random


### Reading QuestionsMRReading Questions

#### 1.

Why does Theorem FTMR deserve the moniker “fundamental”?

#### 2.

Find the matrix representation, $$\matrixrep{T}{B}{C}$$ of the linear transformation

\begin{equation*} \ltdefn{T}{\complex{2}}{\complex{2}},\quad\lteval{T}{\colvector{x_1\\x_2}}=\colvector{2x_1-x_2\\3x_1+2x_2} \end{equation*}

relative to the bases

\begin{align*} B&=\set{\colvector{2\\3},\,\colvector{-1\\2}}& C&=\set{\colvector{1\\0},\,\colvector{1\\1}}\text{.} \end{align*}

#### 3.

What is the second “surprise,” and why is it surprising?

### ExercisesMRExercises

#### C10.

Example KVMR concludes with a basis for the kernel of the linear transformation $$T\text{.}$$ Compute the value of $$T$$ for each of these two basis vectors. Did you get what you expected?

#### C20.

Compute the matrix representation of $$T$$ relative to the bases $$B$$ and $$C\text{,}$$

\begin{gather*} \ltdefn{T}{P_3}{\complex{3}},\quad \lteval{T}{a+bx+cx^2+dx^3}= \colvector{2a-3b+4c-2d\\a+b-c+d\\3a+2c-3d}\\ B=\set{1,\,x,\,x^2,\,x^3}\quad\quad C=\set{\colvector{1\\0\\0},\,\colvector{1\\1\\0},\,\colvector{1\\1\\1}}\text{.} \end{gather*}
Solution.

Apply Definition MR,

\begin{align*} \vectrep{C}{\lteval{T}{1}}&=\vectrep{C}{\colvector{2\\1\\3}} =\vectrep{C}{1\colvector{1\\0\\0}+(-2)\colvector{1\\1\\0}+3\colvector{1\\1\\1}}=\colvector{1\\-2\\3}\\ \vectrep{C}{\lteval{T}{x}}&=\vectrep{C}{\colvector{-3\\1\\0}} =\vectrep{C}{(-4)\colvector{1\\0\\0}+1\colvector{1\\1\\0}+0\colvector{1\\1\\1}}=\colvector{-4\\1\\0}\\ \vectrep{C}{\lteval{T}{x^2}}&=\vectrep{C}{\colvector{4\\-1\\2}} =\vectrep{C}{5\colvector{1\\0\\0}+(-3)\colvector{1\\1\\0}+2\colvector{1\\1\\1}}=\colvector{5\\-3\\2}\\ \vectrep{C}{\lteval{T}{x^3}}&=\vectrep{C}{\colvector{-2\\1\\-3}} =\vectrep{C}{(-3)\colvector{1\\0\\0}+4\colvector{1\\1\\0}+(-3)\colvector{1\\1\\1}}=\colvector{-3\\4\\-3}\text{.} \end{align*}

These four vectors are the columns of the matrix representation,

\begin{equation*} \matrixrep{T}{B}{C}= \begin{bmatrix} 1 & -4 & 5 & -3\\ -2 & 1 & -3 & 4\\ 3 & 0 & 2 & -3 \end{bmatrix}\text{.} \end{equation*}

#### C21.

Find a matrix representation of the linear transformation $$T$$ relative to the bases $$B$$ and $$C\text{,}$$

\begin{align*} &\ltdefn{T}{P_2}{\complex{2}},\quad \lteval{T}{p(x)}=\colvector{p(1)\\p(3)}\\ &B=\set{2-5x+x^2,\,1+x-x^2,\,x^2}\\ &C=\set{\colvector{3\\4},\,\colvector{2\\3}}\text{.} \end{align*}
Solution.

Applying Definition MR,

\begin{align*} \vectrep{C}{\lteval{T}{2-5x+x^2}}&=\vectrep{C}{\colvector{-2\\-4}} =\vectrep{C}{2\colvector{3\\4}+(-4)\colvector{2\\3}}=\colvector{2\\-4}\\ \vectrep{C}{\lteval{T}{1+x-x^2}}&=\vectrep{C}{\colvector{1\\-5}} =\vectrep{C}{13\colvector{3\\4}+(-19)\colvector{2\\3}}=\colvector{13\\-19}\\ \vectrep{C}{\lteval{T}{x^2}}&=\vectrep{C}{\colvector{1\\9}} =\vectrep{C}{(-15)\colvector{3\\4}+23\colvector{2\\3}}=\colvector{-15\\23}\text{.} \end{align*}

So the resulting matrix representation is

\begin{align*} \matrixrep{T}{B}{C} &= \begin{bmatrix} 2 & 13 & -15 \\ -4 & -19 & 23 \end{bmatrix}\text{.} \end{align*}

#### C22.

Let $$S_{22}$$ be the vector space of $$2\times 2$$ symmetric matrices. Build the matrix representation of the linear transformation $$\ltdefn{T}{P_2}{S_{22}}$$ relative to the bases $$B$$ and $$C$$ and then use this matrix representation to compute $$\lteval{T}{3+5x-2x^2}\text{,}$$

\begin{align*} B&=\set{1,\,1+x,\,1+x+x^2} & C&=\set{ \begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix},\, \begin{bmatrix}0 & 1 \\ 1 & 0\end{bmatrix},\, \begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} }\\ \lteval{T}{a+bx+cx^2}&= \begin{bmatrix} 2a-b+c & a+3b-c \\ a+3b-c & a-c \end{bmatrix}\text{.} \end{align*}
Solution.

Input to $$T$$ the vectors of the basis $$B$$ and coordinatize the outputs relative to $$C\text{,}$$

\begin{align*} \vectrep{C}{\lteval{T}{1}}&= \vectrep{C}{ \begin{bmatrix}2 & 1 \\ 1 & 1\end{bmatrix} } = \vectrep{C}{ 2\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ 1\begin{bmatrix}0 & 1 \\ 1 & 0\end{bmatrix}+ 1\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{2\\1\\1}\\ \vectrep{C}{\lteval{T}{1+x}}&= \vectrep{C}{ \begin{bmatrix}1 & 4 \\ 4 & 1\end{bmatrix} } = \vectrep{C}{ 1\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ 4\begin{bmatrix}0 & 1 \\ 1 & 0\end{bmatrix}+ 1\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{1\\4\\1}\\ \vectrep{C}{\lteval{T}{1+x+x^2}}&= \vectrep{C}{ \begin{bmatrix}2 & 3 \\ 3 & 0\end{bmatrix} } = \vectrep{C}{ 2\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ 3\begin{bmatrix}0 & 1 \\ 1 & 0\end{bmatrix}+ 0\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{2\\3\\0}\text{.} \end{align*}

Applying Definition MR we have the matrix representation

\begin{equation*} \matrixrep{T}{B}{C}= \begin{bmatrix}2 & 1 & 2 \\ 1 & 4 & 3 \\ 1 & 1 & 0 \end{bmatrix}\text{.} \end{equation*}

To compute $$\lteval{T}{3+5x-2x^2}$$ employ Theorem FTMR,

\begin{align*} \lteval{T}{3+5x-2x^2}&= \vectrepinv{C}{\matrixrep{T}{B}{C}\vectrep{B}{3+5x-2x^2}}\\ &=\vectrepinv{C}{\matrixrep{T}{B}{C}\vectrep{B}{(-2)(1)+7(1+x)+(-2)(1+x+x^2)}}\\ &=\vectrepinv{C}{ \begin{bmatrix}2&1&2\\1&4&3\\1&1&0\end{bmatrix} \colvector{-2\\7\\-2} }\\ &=\vectrepinv{C}{\colvector{-1\\20\\5}}\\ &= (-1)\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ 20\begin{bmatrix}0 & 1 \\ 1 & 0\end{bmatrix}+ 5\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix}\\ &= \begin{bmatrix}-1 & 20 \\ 20 & 5\end{bmatrix}\text{.} \end{align*}

You can, of course, check your answer by evaluating $$\lteval{T}{3+5x-2x^2}$$ directly.

#### C25.

Use a matrix representation to determine if the linear transformation $$\ltdefn{T}{P_3}{M_{22}}$$ is surjective.

\begin{equation*} \lteval{T}{a+bx+cx^2+dx^3}= \begin{bmatrix} -a+4b+c+2d & 4a-b+6c-d\\ a+5b-2c+2d & a+2c+5d \end{bmatrix} \end{equation*}
Solution.

Choose bases $$B$$ and $$C$$ for the matrix representation,

\begin{align*} B&=\set{1,\,x,\,x^2,\,x^3} & C&=\set{ \begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix},\, \begin{bmatrix}0 & 1 \\ 0 & 0\end{bmatrix},\, \begin{bmatrix}0 & 0 \\ 1 & 0\end{bmatrix},\, \begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} }\text{.} \end{align*}

Input to $$T$$ the vectors of the basis $$B$$ and coordinatize the outputs relative to $$C\text{,}$$

\begin{align*} \vectrep{C}{\lteval{T}{1}}&= \vectrep{C}{\begin{bmatrix}-1 & 4 \\ 1 & 1\end{bmatrix}} = \vectrep{C}{ (-1)\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ 4\begin{bmatrix}0 & 1 \\ 0 & 0\end{bmatrix}+ 1\begin{bmatrix}0 & 0 \\ 1 & 0\end{bmatrix}+ 1\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{-1\\4\\1\\1}\\ \vectrep{C}{\lteval{T}{x}}&= \vectrep{C}{\begin{bmatrix}4 & -1 \\ 5 & 0\end{bmatrix}} = \vectrep{C}{ 4\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ (-1)\begin{bmatrix}0 & 1 \\ 0 & 0\end{bmatrix}+ 5\begin{bmatrix}0 & 0 \\ 1 & 0\end{bmatrix}+ 0\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{4\\-1\\5\\0}\\ \vectrep{C}{\lteval{T}{x^2}}&= \vectrep{C}{\begin{bmatrix}1 & 6 \\ -2 & 2\end{bmatrix}} = \vectrep{C}{ 1\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ 6\begin{bmatrix}0 & 1 \\ 0 & 0\end{bmatrix}+ (-2)\begin{bmatrix}0 & 0 \\ 1 & 0\end{bmatrix}+ 2\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{1\\6\\-2\\2}\\ \vectrep{C}{\lteval{T}{x^3}}&= \vectrep{C}{\begin{bmatrix}2 & -1 \\ 2 & 5\end{bmatrix}} = \vectrep{C}{ 2\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ (-1)\begin{bmatrix}0 & 1 \\ 0 & 0\end{bmatrix}+ 2\begin{bmatrix}0 & 0 \\ 1 & 0\end{bmatrix}+ 5\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{2\\-1\\2\\5}\text{.} \end{align*}

Applying Definition MR we have the matrix representation

\begin{equation*} \matrixrep{T}{B}{C}= \begin{bmatrix} -1 & 4 & 1 & 2 \\ 4 & -1 & 6 & -1 \\ 1 & 5 & -2 & 2 \\ 1 & 0 & 2 & 5 \end{bmatrix}\text{.} \end{equation*}

Properties of this matrix representation will translate to properties of the linear transformation. The matrix representation is nonsingular since it row-reduces to the identity matrix (Theorem NMRRI) and therefore has a column space equal to $$\complex{4}$$ (Theorem CNMB). The column space of the matrix representation is isomorphic to the range of the linear transformation (Theorem RCSI). So the range of $$T$$ has dimension 4, equal to the dimension of the codomain $$M_{22}\text{.}$$ By Theorem ROSLT, $$T$$ is surjective.

#### C30.

Find bases for the kernel and range of the linear transformation $$S$$ below.

\begin{equation*} \ltdefn{S}{M_{22}}{P_2},\quad\lteval{S}{\begin{bmatrix}a&b\\c&d\end{bmatrix}}= (a+2b+5c-4d)+(3a-b+8c+2d)x+(a+b+4c-2d)x^2 \end{equation*}
Solution.

These subspaces will be easiest to construct by analyzing a matrix representation of $$S\text{.}$$ Since we can use any matrix representation, we might as well use natural bases that allow us to construct the matrix representation quickly and easily,

\begin{align*} B&=\set{ \begin{bmatrix}1&0\\0&0\end{bmatrix},\, \begin{bmatrix}0&1\\0&0\end{bmatrix},\, \begin{bmatrix}0&0\\1&0\end{bmatrix},\, \begin{bmatrix}0&0\\0&1\end{bmatrix} } & C&=\set{1,\,x,\,x^2} \end{align*}

then we can practically build the matrix representation on sight,

\begin{equation*} \matrixrep{S}{B}{C}= \begin{bmatrix} 1 & 2 & 5 & -4\\ 3 & -1 & 8 & 2\\ 1 & 1 & 4 & -2 \end{bmatrix}\text{.} \end{equation*}

The first step is to find bases for the null space and column space of the matrix representation. Row-reducing the matrix representation we find,

\begin{equation*} \begin{bmatrix} \leading{1} & 0 & 3 & 0\\ 0 & \leading{1} & 1 & -2\\ 0 & 0 & 0 & 0 \end{bmatrix}\text{.} \end{equation*}

So by Theorem BNS and Theorem BCS, we have

\begin{align*} \nsp{\matrixrep{S}{B}{C}}&=\spn{\set{\colvector{-3\\-1\\1\\0},\,\colvector{0\\2\\0\\1}}} & \csp{\matrixrep{S}{B}{C}}&=\spn{\set{\colvector{1\\3\\1},\,\colvector{2\\-1\\1}}}\text{.} \end{align*}

Now, the proofs of Theorem KNSI and Theorem RCSI tell us that we can apply $$\ltinverse{\vectrepname{B}}$$ and $$\vectrepinvname{C}$$ (respectively) to “un-coordinatize” and get bases for the kernel and range of the linear transformation $$S$$ itself,

\begin{align*} \krn{S}&=\spn{\set{ \begin{bmatrix} -3&-1\\1&0 \end{bmatrix},\, \begin{bmatrix} 0&2\\0&1 \end{bmatrix} }} & \rng{S}&=\spn{\set{1+3x+x^2,\,2-x+x^2}}\text{.} \end{align*}

#### C40.

Let $$S_{22}$$ be the set of $$2\times 2$$ symmetric matrices. Verify that the linear transformation $$R$$ is invertible and find $$\ltinverse{R}\text{.}$$

\begin{equation*} \ltdefn{R}{S_{22}}{P_2},\quad\lteval{R}{\begin{bmatrix}a&b\\b&c\end{bmatrix}}= (a-b)+(2a-3b-2c)x+(a-b+c)x^2 \end{equation*}
Solution.

The analysis of $$R$$ will be easiest if we analyze a matrix representation of $$R\text{.}$$ Since we can use any matrix representation, we might as well use natural bases that allow us to construct the matrix representation quickly and easily,

\begin{align*} B&=\set{ \begin{bmatrix}1&0\\0&0\end{bmatrix},\, \begin{bmatrix}0&1\\1&0\end{bmatrix},\, \begin{bmatrix}0&0\\0&1\end{bmatrix} } & C&=\set{1,\,x,\,x^2} \end{align*}

then we can practically build the matrix representation on sight,

\begin{equation*} \matrixrep{R}{B}{C}= \begin{bmatrix} 1 & -1 & 0\\ 2 & -3 & -2\\ 1 & -1 & 1 \end{bmatrix}\text{.} \end{equation*}

This matrix representation is invertible (it has a nonzero determinant of $$-1\text{,}$$ Theorem SMZD, Theorem NI) so Theorem IMR tells us that the linear transformation $$R$$ is also invertible. To find a formula for $$\ltinverse{R}$$ we compute,

\begin{align*} \lteval{\ltinverse{R}}{a+bx+cx^2} &=\vectrepinv{B}{\matrixrep{\ltinverse{R}}{C}{B}\vectrep{C}{a+bx+cx^2}}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrepinv{B}{\inverse{\left(\matrixrep{R}{B}{C}\right)}\vectrep{C}{a+bx+cx^2}}&& \knowl{./knowl/theorem-IMR.html}{\text{Theorem IMR}}\\ &=\vectrepinv{B}{\inverse{\left(\matrixrep{R}{B}{C}\right)}\colvector{a\\b\\c}}&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\\ &=\vectrepinv{B}{ \begin{bmatrix}5&-1&-2\\4&-1&-2\\-1&0&1\end{bmatrix} \colvector{a\\b\\c}}&& \knowl{./knowl/definition-MI.html}{\text{Definition MI}}\\ &=\vectrepinv{B}{\colvector{5a-b-2c\\4a-b-2c\\-a+c}}&& \knowl{./knowl/definition-MVP.html}{\text{Definition MVP}}\\ &=\begin{bmatrix}5a-b-2c&4a-b-2c\\4a-b-2c&-a+c\end{bmatrix}&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\text{.} \end{align*}

#### C41.

Prove that the linear transformation $$S$$ is invertible. Then find a formula for the inverse linear transformation, $$\ltinverse{S}\text{,}$$ by employing a matrix inverse.

\begin{equation*} \ltdefn{S}{P_1}{M_{12}},\quad \lteval{S}{a+bx}= \begin{bmatrix} 3a+b & 2a+b \end{bmatrix} \end{equation*}
Solution.

First, build a matrix representation of $$S$$ (Definition MR). We are free to choose whatever bases we wish, so we should choose ones that are easy to work with, such as

\begin{align*} B&=\set{1,\,x}\\ C&=\set{\begin{bmatrix}1 & 0\end{bmatrix},\,\begin{bmatrix}0 & 1\end{bmatrix}}\text{.} \end{align*}

The resulting matrix representation is then

\begin{equation*} \matrixrep{T}{B}{C}= \begin{bmatrix} 3 & 1\\ 2 & 1 \end{bmatrix}\text{.} \end{equation*}

This matrix is invertible since it has a nonzero determinant, so by Theorem IMR the linear transformation $$S$$ is invertible. We can use the matrix inverse and Theorem IMR to find a formula for the inverse linear transformation,

\begin{align*} \lteval{\ltinverse{S}}{\begin{bmatrix}a&b\end{bmatrix}} &=\vectrepinv{B}{\matrixrep{\ltinverse{S}}{C}{B}\vectrep{C}{\begin{bmatrix}a&b\end{bmatrix}}}&& \knowl{./knowl/theorem-FTMR.html}{\text{Theorem FTMR}}\\ &=\vectrepinv{B}{\inverse{\left(\matrixrep{S}{B}{C}\right)}\vectrep{C}{\begin{bmatrix}a&b\end{bmatrix}}}&& \knowl{./knowl/theorem-IMR.html}{\text{Theorem IMR}}\\ &=\vectrepinv{B}{\inverse{\left(\matrixrep{S}{B}{C}\right)}\colvector{a\\b}}&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\\ &=\vectrepinv{B}{ \inverse{\left( \begin{bmatrix} 3 & 1\\ 2 & 1 \end{bmatrix} \right)} \colvector{a\\b}}\\ &=\vectrepinv{B}{ \begin{bmatrix} 1 & -1\\ -2 & 3 \end{bmatrix} \colvector{a\\b}}&& \knowl{./knowl/definition-MI.html}{\text{Definition MI}}\\ &=\vectrepinv{B}{ \colvector{ a-b\\ -2a+3b }}&& \knowl{./knowl/definition-MVP.html}{\text{Definition MVP}}\\ &=(a-b)+(-2a+3b)x&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\text{.} \end{align*}

#### C42.

The linear transformation $$\ltdefn{R}{M_{12}}{M_{21}}$$ is invertible. Use a matrix representation to determine a formula for the inverse linear transformation $$\ltdefn{\ltinverse{R}}{M_{21}}{M_{12}}\text{.}$$

\begin{equation*} \lteval{R}{\begin{bmatrix}a & b\end{bmatrix}} = \begin{bmatrix} a+3b\\ 4a+11b \end{bmatrix} \end{equation*}
Solution.

Choose bases $$B$$ and $$C$$ for $$M_{12}$$ and $$M_{21}$$ (respectively),

\begin{align*} B=\set{ \begin{bmatrix}1 & 0\end{bmatrix},\, \begin{bmatrix}0 & 1\end{bmatrix} } C=\set{ \begin{bmatrix}1 \\ 0\end{bmatrix},\, \begin{bmatrix}0 \\ 1\end{bmatrix} }\text{.} \end{align*}

The resulting matrix representation is

\begin{equation*} \matrixrep{R}{B}{C} = \begin{bmatrix} 1 & 3\\ 4 & 11 \end{bmatrix}\text{.} \end{equation*}

This matrix is invertible (its determinant is nonzero, Theorem SMZD), so by Theorem IMR, we can compute the matrix representation of $$\ltinverse{R}$$ with a matrix inverse (Theorem TTMI),

\begin{equation*} \matrixrep{\ltinverse{R}}{C}{B} =\inverse{\begin{bmatrix}1 & 3\\4 & 11\end{bmatrix}} =\begin{bmatrix}-11 & 3\\4 & -1\end{bmatrix}\text{.} \end{equation*}

To obtain a general formula for $$\ltinverse{R}\text{,}$$ use Theorem FTMR,

\begin{align*} \lteval{\ltinverse{R}}{\begin{bmatrix}x \\ y\end{bmatrix}} &=\vectrepinv{B}{\matrixrep{\ltinverse{R}}{C}{B}\vectrep{C}{\begin{bmatrix}x\\y\end{bmatrix}}}\\ &=\vectrepinv{B}{\begin{bmatrix}-11&3\\4&-1\end{bmatrix}\colvector{x\\y}}\\ &=\vectrepinv{B}{\colvector{-11x+3y\\4x-y}}\\ &=\begin{bmatrix}-11x+3y&4x-y\end{bmatrix}\text{.} \end{align*}

#### C50.

Use a matrix representation to find a basis for the range of the linear transformation $$L\text{.}$$

\begin{equation*} \ltdefn{L}{M_{22}}{P_2},\quad \lteval{L}{\begin{bmatrix}a&b\\c&d\end{bmatrix}}= (a+2b+4c+d)+(3a+c-2d)x+(-a+b+3c+3d)x^2 \end{equation*}
Solution.

As usual, build any matrix representation of $$L\text{,}$$ most likely using “nice” bases, such as

\begin{align*} B&=\set{ \begin{bmatrix} 1 & 0 \\ 0 & 0 \end{bmatrix},\, \begin{bmatrix} 0 & 1 \\ 0 & 0 \end{bmatrix},\, \begin{bmatrix} 0 & 0 \\ 1 & 0 \end{bmatrix},\, \begin{bmatrix} 0 & 0 \\ 0 & 1 \end{bmatrix} }\\ C&=\set{1,\,x,\,x^2}\text{.} \end{align*}

Then the matrix representation (Definition MR) is,

\begin{equation*} \matrixrep{L}{B}{C}= \begin{bmatrix} 1 & 2 & 4 & 1 \\ 3 & 0 & 1 & -2 \\ -1 & 1 & 3 & 3 \end{bmatrix}\text{.} \end{equation*}

Theorem RCSI tells us that we can compute the column space of the matrix representation, then use the isomorphism $$\vectrepinvname{C}$$ to convert the column space of the matrix representation into the range of the linear transformation. So we first analyze the matrix representation,

\begin{equation*} \begin{bmatrix} 1 & 2 & 4 & 1 \\ 3 & 0 & 1 & -2 \\ -1 & 1 & 3 & 3 \end{bmatrix} \rref \begin{bmatrix} \leading{1} & 0 & 0 & -1 \\ 0 & \leading{1} & 0 & -1 \\ 0 & 0 & \leading{1} & 1 \end{bmatrix}\text{.} \end{equation*}

With three nonzero rows in the reduced row-echelon form of the matrix, we know the column space has dimension 3. Since $$P_2$$ has dimension 3 (Theorem DP), the range must be all of $$P_2\text{.}$$ So any basis of $$P_2$$ would suffice as a basis for the range. For instance, $$C$$ itself would be a correct answer.

A more laborious approach would be to use Theorem BCS and choose the first three columns of the matrix representation as a basis for the range of the matrix representation. These could then be “un-coordinatized” with $$\vectrepinvname{C}$$ to yield a (“not nice”) basis for $$P_2\text{.}$$

#### C51.

Use a matrix representation to find a basis for the kernel of the linear transformation $$L\text{.}$$

\begin{equation*} \ltdefn{L}{M_{22}}{P_2},\quad \lteval{L}{\begin{bmatrix}a&b\\c&d\end{bmatrix}}= (a+2b+4c+d)+(3a+c-2d)x+(-a+b+3c+3d)x^2 \end{equation*}

#### C52.

Find a basis for the kernel of the linear transformation $$\ltdefn{T}{P_2}{M_{22}}\text{.}$$

\begin{equation*} \lteval{T}{a+bx+cx^2}= \begin{bmatrix} a+2b-2c & 2a+2b \\ -a+b-4c & 3a+2b+2c \end{bmatrix} \end{equation*}
Solution.

Choose bases $$B$$ and $$C$$ for the matrix representation,

\begin{align*} B&=\set{1,\,x,\,x^2} & C&=\set{ \begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix},\, \begin{bmatrix}0 & 1 \\ 0 & 0\end{bmatrix},\, \begin{bmatrix}0 & 0 \\ 1 & 0\end{bmatrix},\, \begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} }\text{.} \end{align*}

Input to $$T$$ the vectors of the basis $$B$$ and coordinatize the outputs relative to $$C\text{,}$$

\begin{align*} \vectrep{C}{\lteval{T}{1}}&= \vectrep{C}{\begin{bmatrix}1 & 2 \\ -1 & 3\end{bmatrix}} = \vectrep{C}{ 1\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ 2\begin{bmatrix}0 & 1 \\ 0 & 0\end{bmatrix}+ (-1)\begin{bmatrix}0 & 0 \\ 1 & 0\end{bmatrix}+ 3\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{1\\2\\-1\\3}\\ \vectrep{C}{\lteval{T}{x}}&= \vectrep{C}{\begin{bmatrix}2 & 2 \\ 1 & 2 \end{bmatrix}} = \vectrep{C}{ 2\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ 2\begin{bmatrix}0 & 1 \\ 0 & 0\end{bmatrix}+ 1\begin{bmatrix}0 & 0 \\ 1 & 0\end{bmatrix}+ 2\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{2\\2\\1\\2}\\ \vectrep{C}{\lteval{T}{x^2}}&= \vectrep{C}{\begin{bmatrix}-2 & 0 \\ -4 & 2\end{bmatrix}} = \vectrep{C}{ (-2)\begin{bmatrix}1 & 0 \\ 0 & 0\end{bmatrix}+ 0\begin{bmatrix}0 & 1 \\ 0 & 0\end{bmatrix}+ (-4)\begin{bmatrix}0 & 0 \\ 1 & 0\end{bmatrix}+ 2\begin{bmatrix}0 & 0 \\ 0 & 1\end{bmatrix} } = \colvector{-2\\0\\-4\\2}\text{.} \end{align*}

Applying Definition MR we have the matrix representation

\begin{equation*} \matrixrep{T}{B}{C}= \begin{bmatrix} 1 & 2 & -2 \\ 2 & 2 & 0 \\ -1 & 1 & -4 \\ 3 & 2 & 2 \end{bmatrix}\text{.} \end{equation*}

The null space of the matrix representation is isomorphic (via $$\vectrepname{B}$$) to the kernel of the linear transformation (Theorem KNSI). So we compute the null space of the matrix representation by first row-reducing the matrix to,

\begin{equation*} \begin{bmatrix} \leading{1} & 0 & 2 \\ 0 & \leading{1} & -2 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix}\text{.} \end{equation*}

Employing Theorem BNS we have

\begin{equation*} \nsp{\matrixrep{T}{B}{C}}=\spn{\set{\colvector{-2\\2\\1}}}\text{.} \end{equation*}

We only need to uncoordinatize this one basis vector to get a basis for $$\krn{T}\text{,}$$

\begin{equation*} \krn{T} =\spn{\set{\vectrepinv{B}{\colvector{-2\\2\\1}}}} =\spn{\set{-2+2x+x^2}}\text{.} \end{equation*}

#### M20.

The linear transformation $$D$$ performs differentiation on polynomials. Use a matrix representation of $$D$$ to find the rank and nullity of $$D\text{.}$$

\begin{equation*} \ltdefn{D}{P_n}{P_n},\quad \lteval{D}{p(x)}=p^\prime(x) \end{equation*}
Solution.

Build a matrix representation (Definition MR) with the set

\begin{equation*} B=\set{1,\,x,\,x^2,\,\dots,\,x^n} \end{equation*}

employed as a basis of both the domain and codomain. Then

\begin{align*} \vectrep{B}{\lteval{D}{1}}&=\vectrep{B}{0} =\colvector{0\\0\\0\\\vdots\\0\\0} & \vectrep{B}{\lteval{D}{x}}&=\vectrep{B}{1} =\colvector{1\\0\\0\\\vdots\\0\\0}\\ \vectrep{B}{\lteval{D}{x^2}}&=\vectrep{B}{2x} =\colvector{0\\2\\0\\\vdots\\0\\0} & \vectrep{B}{\lteval{D}{x^3}}&=\vectrep{B}{3x^2} =\colvector{0\\0\\3\\\vdots\\0\\0}\\ &\vdots\\ \vectrep{B}{\lteval{D}{x^n}}&=\vectrep{B}{nx^{n-1}} =\colvector{0\\0\\0\\\vdots\\n\\0} \end{align*}

and the resulting matrix representation is

\begin{equation*} \matrixrep{D}{B}{B}= \begin{bmatrix} 0 & 1 & 0 & 0 & \dots & 0 & 0 \\ 0 & 0 & 2 & 0 & \dots & 0 & 0 \\ 0 & 0 & 0 & 3 & \dots & 0 & 0 \\ & \vdots &&& \ddots & & \vdots \\ 0 & 0 & 0 & 0 & \dots & 0 & n \\ 0 & 0 & 0 & 0 & \dots & 0 & 0 \end{bmatrix}\text{.} \end{equation*}

This $$(n+1)\times(n+1)$$ matrix is very close to being in reduced row-echelon form. Multiply row $$i$$ by $$\frac{1}{i}\text{,}$$ for $$1\leq i\leq n\text{,}$$ to convert it to reduced row-echelon form. From this we can see that matrix representation $$\matrixrep{D}{B}{B}$$ has rank $$n$$ and nullity $$1\text{.}$$ Applying Theorem RCSI and Theorem KNSI tells us that the linear transformation $$D$$ will have the same values for the rank and nullity, as well.

#### M60.

Suppose $$U$$ and $$V$$ are vector spaces and define a function $$\ltdefn{Z}{U}{V}$$ by $$\lteval{Z}{\vect{u}}=\zerovector_{V}$$ for every $$\vect{u}\in U\text{.}$$ Then Exercise IVLT.M60 asks you to formulate the theorem: $$Z$$ is invertible if and only if $$U=\set{\zerovector_U}$$ and $$V=\set{\zerovector_V}\text{.}$$ What would a matrix representation of $$Z$$ look like in this case? How does Theorem IMR read in this case?

#### T20.

Construct a new solution to Exercise B.T50 along the following outline. From the $$n\times n$$ matrix $$A\text{,}$$ construct the linear transformation $$\ltdefn{T}{\complex{n}}{\complex{n}}\text{,}$$ $$\lteval{T}{\vect{x}}=A\vect{x}\text{.}$$ Use Theorem NI, Theorem IMILT and Theorem ILTIS to translate between the nonsingularity of $$A$$ and the surjectivity/injectivity of $$T\text{.}$$ Then apply Theorem ILTB and Theorem SLTB to connect these properties with bases.

Solution.

Given the nonsingular $$n\times n$$ matrix $$A\text{,}$$ create the linear transformation $$\ltdefn{T}{\complex{n}}{\complex{n}}$$ defined by $$\lteval{T}{\vect{x}}=A\vect{x}\text{.}$$ Then

\begin{align*} A\text{ nonsingular} &\iff A\text{ invertible}&& \knowl{./knowl/theorem-NI.html}{\text{Theorem NI}}\\ &\iff T\text{ invertible}&& \knowl{./knowl/theorem-IMILT.html}{\text{Theorem IMILT}}\\ &\iff T\text{ injective and surjective}&& \knowl{./knowl/theorem-ILTIS.html}{\text{Theorem ILTIS}}\\ &\iff C\text{ linearly independent, and}&& \knowl{./knowl/theorem-ILTB.html}{\text{Theorem ILTB}}\\ &\quad\quad\quad\quad C\text{ spans }\complex{n}&& \knowl{./knowl/theorem-SLTB.html}{\text{Theorem SLTB}}\\ &\iff C\text{ basis for }\complex{n}&& \knowl{./knowl/definition-B.html}{\text{Definition B}}\text{.} \end{align*}

#### T40.

Theorem VSLT defines the vector space $$\vslt{U}{V}$$ containing all linear transformations with domain $$U$$ and codomain $$V\text{.}$$ Suppose $$\dimension{U}=n$$ and $$\dimension{V}=m\text{.}$$ Prove that $$\vslt{U}{V}$$ is isomorphic to $$M_{mn}\text{,}$$ the vector space of all $$m\times n$$ matrices (Example VSM). (Hint: we could have suggested this exercise in Chapter LT, but have postponed it to this section. Why?)

#### T41.

Theorem VSLT defines the vector space $$\vslt{U}{V}$$ containing all linear transformations with domain $$U$$ and codomain $$V\text{.}$$ Determine a basis for $$\vslt{U}{V}\text{.}$$ (Hint: study Exercise MR.T40 first.)

#### T60.

Create an entirely different proof of Theorem IMILT that relies on Definition IVLT to establish the invertibility of $$T\text{,}$$ and that relies on Definition MI to establish the invertibility of $$A\text{.}$$

#### T80.

Suppose that $$\ltdefn{T}{U}{V}$$ and $$\ltdefn{S}{V}{W}$$ are linear transformations, and that $$B\text{,}$$ $$C$$ and $$D$$ are bases for $$U\text{,}$$ $$V\text{,}$$ and $$W\text{.}$$ Using only Definition MR define matrix representations for $$T$$ and $$S\text{.}$$ Using these two definitions, and Definition MR, derive a matrix representation for the composition $$\compose{S}{T}$$ in terms of the entries of the matrices $$\matrixrep{T}{B}{C}$$ and $$\matrixrep{S}{C}{D}\text{.}$$ Explain how you would use this result to motivate a definition for matrix multiplication that is strikingly similar to Theorem EMP.

Solution.

Suppose that $$B=\set{\vectorlist{u}{m}}\text{,}$$ $$C=\set{\vectorlist{v}{n}}$$ and $$D=\set{\vectorlist{w}{p}}\text{.}$$ For convenience, set $$M=\matrixrep{T}{B}{C}\text{,}$$ $$m_{ij}=\matrixentry{M}{ij}\text{,}$$ $$1\leq i\leq n\text{,}$$ $$1\leq j\leq m\text{,}$$ and similarly, set $$N=\matrixrep{S}{C}{D}\text{,}$$ $$n_{ij}=\matrixentry{N}{ij}\text{,}$$ $$1\leq i\leq p\text{,}$$ $$1\leq j\leq n\text{.}$$ We want to learn about the matrix representation of $$\ltdefn{\compose{S}{T}}{V}{W}$$ relative to $$B$$ and $$D\text{.}$$ (Zach Jarvis helped with this proof.)

We will examine a single (generic) entry of this representation.

\begin{align*} \matrixentry{\matrixrep{\compose{S}{T}}{B}{D}}{ij} &=\vectorentry{\vectrep{D}{\lteval{\left(\compose{S}{T}\right)}{\vect{u}_j}}}{i}&& \knowl{./knowl/definition-MR.html}{\text{Definition MR}}\\ &=\vectorentry{\vectrep{D}{\lteval{S}{\lteval{T}{\vect{u}_j}}}}{i}&& \knowl{./knowl/definition-LTC.html}{\text{Definition LTC}}\\ &=\vectorentry{\vectrep{D}{\lteval{S}{\sum_{k=1}^{n}m_{kj}\vect{v}_k}}}{i}&& \knowl{./knowl/definition-MR.html}{\text{Definition MR}}\\ &=\vectorentry{\vectrep{D}{\sum_{k=1}^{n}m_{kj}\lteval{S}{\vect{v}_k}}}{i}&& \knowl{./knowl/theorem-LTLC.html}{\text{Theorem LTLC}}\\ &=\vectorentry{\vectrep{D}{\sum_{k=1}^{n}m_{kj}\sum_{\ell=1}^{p}n_{\ell k}\vect{w}_\ell}}{i}&& \knowl{./knowl/definition-MR.html}{\text{Definition MR}}\\ &=\vectorentry{\vectrep{D}{\sum_{k=1}^{n}\sum_{\ell=1}^{p}m_{kj}n_{\ell k}\vect{w}_\ell}}{i}&& \knowl{./knowl/property-DVA.html}{\text{Property DVA}}\\ &=\vectorentry{\vectrep{D}{\sum_{\ell=1}^{p}\sum_{k=1}^{n}m_{kj}n_{\ell k}\vect{w}_\ell}}{i}&& \knowl{./knowl/property-C.html}{\text{Property C}}\\ &=\vectorentry{\vectrep{D}{\sum_{\ell=1}^{p}\left(\sum_{k=1}^{n}m_{kj}n_{\ell k}\right)\vect{w}_\ell}}{i}&& \knowl{./knowl/property-DSA.html}{\text{Property DSA}}\\ &=\vectorentry{ \colvector{ \sum_{k=1}^{n}m_{kj}n_{1k}\\ \sum_{k=1}^{n}m_{kj}n_{2k}\\ \sum_{k=1}^{n}m_{kj}n_{3k}\\ \vdots\\ \sum_{k=1}^{n}m_{kj}n_{pk}} }{i}&& \knowl{./knowl/definition-VR.html}{\text{Definition VR}}\\ &=\sum_{k=1}^{n}m_{kj}n_{ik}&& \knowl{./knowl/definition-CV.html}{\text{Definition CV}}\\ &=\sum_{k=1}^{n}n_{ik}m_{kj}&& \knowl{./knowl/property-CMCN.html}{\text{Property CMCN}}\\ &=\sum_{k=1}^{n}\matrixentry{\matrixrep{S}{C}{D}}{ik}\matrixentry{\matrixrep{T}{B}{C}}{kj}&& \knowl{./knowl/property-CMCN.html}{\text{Property CMCN}}\text{.} \end{align*}

This formula for the entry of a matrix should remind you of Theorem EMP. However, while the theorem presumed we knew how to multiply matrices, the solution before us never uses any understanding of matrix products. It uses the definitions of vector and matrix representations, properties of linear transformations and vector spaces. So if we began a course by first discussing vector space, and then linear transformations between vector spaces, we could carry matrix representations into a motivation for a definition of matrix multiplication that is grounded in function composition. That is worth saying again — a definition of matrix representations of linear transformations results in a matrix product being the representation of a composition of linear transformations.

This exercise is meant to explain why many authors take the formula in Theorem EMP as their definition of matrix multiplication, and why it is a natural choice when the proper motivation is in place. If we first defined matrix multiplication in the style of Theorem EMP, then the above argument, followed by a simple application of the definition of matrix equality (Definition ME), would yield Theorem MRCLT.