### Section FS  Four Subsets

From A First Course in Linear Algebra
Version 2.99
http://linear.ups.edu/

There are four natural subsets associated with a matrix. We have met three already: the null space, the column space and the row space. In this section we will introduce a fourth, the left null space. The objective of this section is to describe one procedure that will allow us to find linearly independent sets that span each of these four sets of column vectors. Along the way, we will make a connection with the inverse of a matrix, so Theorem FS will tie together most all of this chapter (and the entire course so far).

#### Subsection LNS: Left Null Space

Definition LNS
Left Null Space
Suppose A is an m × n matrix. Then the left null space is defined as ℒ\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left ({A}^{t}\right ) ⊆ {ℂ}^{m}.

(This definition contains Notation LNS.)

The left null space will not feature prominently in the sequel, but we can explain its name and connect it to row operations. Suppose y ∈ℒ\kern -1.95872pt \left (A\right ). Then by Definition LNS, {A}^{t}y = 0. We can then write

\eqalignno{ {0}^{t} & ={ \left ({A}^{t}y\right )}^{t} & &\text{@(a href="#definition.LNS")Definition LNS@(/a)} & & & & \cr & = {y}^{t}{\left ({A}^{t}\right )}^{t} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.MMT")Theorem MMT@(/a)} & & & & \cr & = {y}^{t}A & &\text{@(a href="fcla-jsmath-2.99li30.html#theorem.TT")Theorem TT@(/a)} & & & & }

The product {y}^{t}A can be viewed as the components of y acting as the scalars in a linear combination of the rows of A. And the result is a “row vector”, {0}^{t} that is totally zeros. When we apply a sequence of row operations to a matrix, each row of the resulting matrix is some linear combination of the rows. These observations tell us that the vectors in the left null space are scalars that record a sequence of row operations that result in a row of zeros in the row-reduced version of the matrix. We will see this idea more explicitly in the course of proving Theorem FS.

Example LNS
Left null space
We will find the left null space of

 A = \left [\array{ 1 &−3&1\cr −2 & 1 &1 \cr 1 & 5 &1\cr 9 &−4 &0 } \right ]

We transpose A and row-reduce,

 { A}^{t} = \left [\array{ 1 &−2&1& 9\cr −3 & 1 &5 &−4 \cr 1 & 1 &1& 0 } \right ]\mathop{\longrightarrow}\limits_{}^{\text{RREF}}\left [\array{ \text{1}&0&0& 2\cr 0&\text{1 } &0 &−3 \cr 0&0&\text{1}& 1 } \right ]

Applying Definition LNS and Theorem BNS we have

 ℒ\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left ({A}^{t}\right ) = \left \langle \left \{\left [\array{ −2\cr 3 \cr −1\cr 1 } \right ]\right \}\right \rangle

If you row-reduce A you will discover one zero row in the reduced row-echelon form. This zero row is created by a sequence of row operations, which in total amounts to a linear combination, with scalars {a}_{1} = −2, {a}_{2} = 3, {a}_{3} = −1 and {a}_{4} = 1, on the rows of A and which results in the zero vector (check this!). So the components of the vector describing the left null space of A provide a relation of linear dependence on the rows of A.

#### Subsection CCS: Computing Column Spaces

We have three ways to build the column space of a matrix. First, we can use just the definition, Definition CSM, and express the column space as a span of the columns of the matrix. A second approach gives us the column space as the span of some of the columns of the matrix, but this set is linearly independent (Theorem BCS). Finally, we can transpose the matrix, row-reduce the transpose, kick out zero rows, and transpose the remaining rows back into column vectors. Theorem CSRST and Theorem BRS tell us that the resulting vectors are linearly independent and their span is the column space of the original matrix.

We will now demonstrate a fourth method by way of a rather complicated example. Study this example carefully, but realize that its main purpose is to motivate a theorem that simplifies much of the apparent complexity. So other than an instructive exercise or two, the procedure we are about to describe will not be a usual approach to computing a column space.

Example CSANS
Column space as null space
Lets find the column space of the matrix A below with a new approach.

 A = \left [\array{ 10 & 0 & 3 & 8 & 7\cr −16 &−1 &−4 &−10 &−13 \cr −6 & 1 &−3& −6 & −6\cr 0 & 2 &−2 & −3 & −2 \cr 3 & 0 & 1 & 2 & 3\cr −1 &−1 & 1 & 1 & 0 } \right ]

By Theorem CSCS we know that the column vector b is in the column space of A if and only if the linear system ℒS\kern -1.95872pt \left (A,\kern 1.95872pt b\right ) is consistent. So let’s try to solve this system in full generality, using a vector of variables for the vector of constants. In other words, which vectors b lead to consistent systems? Begin by forming the augmented matrix \left [\left .A\kern 1.95872pt \right \vert \kern 1.95872pt b\right ] with a general version of b,

 \left [\left .A\kern 1.95872pt \right \vert \kern 1.95872pt b\right ] = \left [\array{ 10 & 0 & 3 & 8 & 7 &{b}_{1} \cr −16&−1&−4&−10&−13&{b}_{2} \cr −6 & 1 &−3& −6 & −6 &{b}_{3} \cr 0 & 2 &−2& −3 & −2 &{b}_{4} \cr 3 & 0 & 1 & 2 & 3 &{b}_{5} \cr −1 &−1& 1 & 1 & 0 &{b}_{6} } \right ]

To identify solutions we will row-reduce this matrix and bring it to reduced row-echelon form. Despite the presence of variables in the last column, there is nothing to stop us from doing this. Except our numerical routines on calculators can’t be used, and even some of the symbolic algebra routines do some unexpected maneuvers with this computation. So do it by hand. Yes, it is a bit of work. But worth it. We’ll still be here when you get back. Notice along the way that the row operations are exactly the same ones you would do if you were just row-reducing the coefficient matrix alone, say in connection with a homogeneous system of equations. The column with the {b}_{i} acts as a sort of bookkeeping device. There are many different possibilities for the result, depending on what order you choose to perform the row operations, but shortly we’ll all be on the same page. Here’s one possibility (you can find this same result by doing additional row operations with the fifth and sixth rows to remove any occurrences of {b}_{1} and {b}_{2} from the first four rows of your result):

\eqalignno{ \left [\array{ \text{1}&0&0&0& 2 & {b}_{3} − {b}_{4} + 2{b}_{5} − {b}_{6} \cr 0&\text{1}&0&0&−3&−2{b}_{3} + 3{b}_{4} − 3{b}_{5} + 3{b}_{6} \cr 0&0&\text{1}&0& 1 & {b}_{3} + {b}_{4} + 3{b}_{5} + 3{b}_{6} \cr 0&0&0&\text{1}&−2& −2{b}_{3} + {b}_{4} − 4{b}_{5} \cr 0&0&0&0& 0 &{b}_{1} + 3{b}_{3} − {b}_{4} + 3{b}_{5} + {b}_{6} \cr 0&0&0&0& 0 & {b}_{2} − 2{b}_{3} + {b}_{4} + {b}_{5} − {b}_{6} } \right ] & & }

Our goal is to identify those vectors b which make ℒS\kern -1.95872pt \left (A,\kern 1.95872pt b\right ) consistent. By Theorem RCLS we know that the consistent systems are precisely those without a leading 1 in the last column. Are the expressions in the last column of rows 5 and 6 equal to zero, or are they leading 1’s? The answer is: maybe. It depends on b. With a nonzero value for either of these expressions, we would scale the row and produce a leading 1. So we get a consistent system, and b is in the column space, if and only if these two expressions are both simultaneously zero. In other words, members of the column space of A are exactly those vectors b that satisfy

\eqalignno{ {b}_{1} + 3{b}_{3} − {b}_{4} + 3{b}_{5} + {b}_{6} & = 0 & & \cr {b}_{2} − 2{b}_{3} + {b}_{4} + {b}_{5} − {b}_{6} & = 0 & & }

Hmmm. Looks suspiciously like a homogeneous system of two equations with six variables. If you’ve been playing along (and we hope you have) then you may have a slightly different system, but you should have just two equations. Form the coefficient matrix and row-reduce (notice that the system above has a coefficient matrix that is already in reduced row-echelon form). We should all be together now with the same matrix,

 L = \left [\array{ \text{1}&0& 3 &−1&3& 1\cr 0&\text{1 } &−2 & 1 &1 &−1 } \right ]

So, C\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left (L\right ) and we can apply Theorem BNS to obtain a linearly independent set to use in a span construction,

 C\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left (L\right ) = \left \langle \left \{\left [\array{ −3\cr 2 \cr 1\cr 0 \cr 0\cr 0 } \right ],\kern 1.95872pt \left [\array{ 1\cr −1 \cr 0\cr 1 \cr 0\cr 0 } \right ],\kern 1.95872pt \left [\array{ −3\cr −1 \cr 0\cr 0 \cr 1\cr 0 } \right ],\kern 1.95872pt \left [\array{ −1\cr 1 \cr 0\cr 0 \cr 0\cr 1 } \right ]\right \}\right \rangle

Whew! As a postscript to this central example, you may wish to convince yourself that the four vectors above really are elements of the column space? Do they create consistent systems with A as coefficient matrix? Can you recognize the constant vector in your description of these solution sets?

OK, that was so much fun, let’s do it again. But simpler this time. And we’ll all get the same results all the way through. Doing row operations by hand with variables can be a bit error prone, so let’s see if we can improve the process some. Rather than row-reduce a column vector b full of variables, let’s write b = {I}_{6}b and we will row-reduce the matrix {I}_{6} and when we finish row-reducing, then we will compute the matrix-vector product. You should first convince yourself that we can operate like this (this is the subject of a future homework exercise). Rather than augmenting A with b, we will instead augment it with {I}_{6} (does this feel familiar?),

 M = \left [\array{ 10 & 0 & 3 & 8 & 7 &1&0&0&0&0&0\cr −16 &−1 &−4 &−10 &−13 &0 &1 &0 &0 &0 &0 \cr −6 & 1 &−3& −6 & −6 &0&0&1&0&0&0\cr 0 & 2 &−2 & −3 & −2 &0 &0 &0 &1 &0 &0 \cr 3 & 0 & 1 & 2 & 3 &0&0&0&0&1&0\cr −1 &−1 & 1 & 1 & 0 &0 &0 &0 &0 &0 &1 } \right ]

We want to row-reduce the left-hand side of this matrix, but we will apply the same row operations to the right-hand side as well. And once we get the left-hand side in reduced row-echelon form, we will continue on to put leading 1’s in the final two rows, as well as clearing out the columns containing those two additional leading 1’s. It is these additional row operations that will ensure that we all get to the same place, since the reduced row-echelon form is unique (Theorem RREFU),

 N = \left [\array{ 1&0&0&0& 2 &0&0& 1 &−1& 2 &−1\cr 0&1 &0 &0 &−3 &0 &0 &−2 & 3 &−3 & 3 \cr 0&0&1&0& 1 &0&0& 1 & 1 & 3 & 3\cr 0&0 &0 &1 &−2 &0 &0 &−2 & 1 &−4 & 0 \cr 0&0&0&0& 0 &1&0& 3 &−1& 3 & 1\cr 0&0 &0 &0 & 0 &0 &1 &−2 & 1 & 1 &−1 } \right ]

We are after the final six columns of this matrix, which we will multiply by b

 J = \left [\array{ 0&0& 1 &−1& 2 &−1\cr 0&0 &−2 & 3 &−3 & 3 \cr 0&0& 1 & 1 & 3 & 3\cr 0&0 &−2 & 1 &−4 & 0 \cr 1&0& 3 &−1& 3 & 1\cr 0&1 &−2 & 1 & 1 &−1 } \right ]

so

 Jb = \left [\array{ 0&0& 1 &−1& 2 &−1\cr 0&0 &−2 & 3 &−3 & 3 \cr 0&0& 1 & 1 & 3 & 3\cr 0&0 &−2 & 1 &−4 & 0 \cr 1&0& 3 &−1& 3 & 1\cr 0&1 &−2 & 1 & 1 &−1 } \right ]\left [\array{ {b}_{1} \cr {b}_{2} \cr {b}_{3} \cr {b}_{4} \cr {b}_{5} \cr {b}_{6} } \right ] = \left [\array{ {b}_{3} − {b}_{4} + 2{b}_{5} − {b}_{6} \cr −2{b}_{3} + 3{b}_{4} − 3{b}_{5} + 3{b}_{6} \cr {b}_{3} + {b}_{4} + 3{b}_{5} + 3{b}_{6} \cr −2{b}_{3} + {b}_{4} − 4{b}_{5} \cr {b}_{1} + 3{b}_{3} − {b}_{4} + 3{b}_{5} + {b}_{6} \cr {b}_{2} − 2{b}_{3} + {b}_{4} + {b}_{5} − {b}_{6}\cr } \right ]

So by applying the same row operations that row-reduce A to the identity matrix (which we could do with a calculator once {I}_{6} is placed alongside of A), we can then arrive at the result of row-reducing a column of symbols where the vector of constants usually resides. Since the row-reduced version of A has two zero rows, for a consistent system we require that

\eqalignno{ {b}_{1} + 3{b}_{3} − {b}_{4} + 3{b}_{5} + {b}_{6} & = 0 & & \cr {b}_{2} − 2{b}_{3} + {b}_{4} + {b}_{5} − {b}_{6} & = 0 & & }

Now we are exactly back where we were on the first go-round. Notice that we obtain the matrix L as simply the last two rows and last six columns of N.

This example motivates the remainder of this section, so it is worth careful study. You might attempt to mimic the second approach with the coefficient matrices of Archetype I and Archetype J. We will see shortly that the matrix L contains more information about A than just the column space.

#### Subsection EEF: Extended echelon form

The final matrix that we row-reduced in Example CSANS should look familiar in most respects to the procedure we used to compute the inverse of a nonsingular matrix, Theorem CINM. We will now generalize that procedure to matrices that are not necessarily nonsingular, or even square. First a definition.

Definition EEF
Extended Echelon Form
Suppose A is an m × n matrix. Extend A on its right side with the addition of an m × m identity matrix to form an m × (n + m) matrix M. Use row operations to bring M to reduced row-echelon form and call the result N. N is the extended reduced row-echelon form of A, and we will standardize on names for five submatrices (B, C, J, K, L) of N.

Let B denote the m × n matrix formed from the first n columns of N and let J denote the m × m matrix formed from the last m columns of N. Suppose that B has r nonzero rows. Further partition N by letting C denote the r × n matrix formed from all of the non-zero rows of B. Let K be the r × m matrix formed from the first r rows of J, while L will be the (m − r) × m matrix formed from the bottom m − r rows of J. Pictorially,

 M = [A\vert {I}_{m}]\mathop{\longrightarrow}\limits_{}^{\text{RREF}}N = [B\vert J] = \left [\array{ C&K ̲ ̲ \cr 0 &L } \right ]

Example SEEF
Submatrices of extended echelon form
We illustrate Definition EEF with the matrix A,

 A = \left [\array{ 1 &−1&−2& 7 & 1 & 6\cr −6 & 2 &−4 &−18 &−3 &−26 \cr 4 &−1& 4 & 10 & 2 & 17\cr 3 &−1 & 2 & 9 & 1 & 12 } \right ]

Augmenting with the 4 × 4 identity matrix, M=

 \left [\array{ 1 &−1&−2& 7 & 1 & 6 &1&0&0&0\cr −6 & 2 &−4 &−18 &−3 &−26 &0 &1 &0 &0 \cr 4 &−1& 4 & 10 & 2 & 17 &0&0&1&0\cr 3 &−1 & 2 & 9 & 1 & 12 &0 &0 &0 &1 } \right ]

and row-reducing, we obtain

 N = \left [\array{ \text{1}&0&2& 1 &0& 3 &0& 1 &1& 1\cr 0&\text{1 } &4 &−6 &0 &−1 &0 & 2 &3 & 0 \cr 0&0&0& 0 &\text{1}& 2 &0&−1&0&−2\cr 0&0 &0 & 0 &0 & 0 &\text{1 } & 2 &2 & 1 } \right ]

So we then obtain

\eqalignno{ B & = \left [\array{ \text{1}&0&2& 1 &0& 3\cr 0&\text{1 } &4 &−6 &0 &−1 \cr 0&0&0& 0 &\text{1}& 2\cr 0&0 &0 & 0 &0 & 0 } \right ] & & \cr C & = \left [\array{ \text{1}&0&2& 1 &0& 3\cr 0&\text{1 } &4 &−6 &0 &−1 \cr 0&0&0& 0 &\text{1}& 2 } \right ] & & \cr J & = \left [\array{ 0& 1 &1& 1\cr 0& 2 &3 & 0 \cr 0&−1&0&−2\cr \text{1 }& 2 &2 & 1 } \right ] & & \cr K & = \left [\array{ 0& 1 &1& 1\cr 0& 2 &3 & 0 \cr 0&−1&0&−2 } \right ] & & \cr L & = \left [\array{ \text{1}&2&2&1 } \right ] & & }

You can observe (or verify) the properties of the following theorem with this example.

Theorem PEEF
Properties of Extended Echelon Form
Suppose that A is an m × n matrix and that N is its extended echelon form. Then

1. J is nonsingular.
2. B = JA.
3. If x ∈ {ℂ}^{n} and y ∈ {ℂ}^{m}, then Ax = y if and only if Bx = Jy.
4. C is in reduced row-echelon form, has no zero rows and has r pivot columns.
5. L is in reduced row-echelon form, has no zero rows and has m − r pivot columns.

Proof   J is the result of applying a sequence of row operations to {I}_{m}, as such J and {I}_{m} are row-equivalent. ℒS\kern -1.95872pt \left ({I}_{m},\kern 1.95872pt 0\right ) has only the zero solution, since {I}_{m} is nonsingular (Theorem NMRRI). Thus, ℒS\kern -1.95872pt \left (J,\kern 1.95872pt 0\right ) also has only the zero solution (Theorem REMES, Definition ESYS) and J is therefore nonsingular (Definition NSM).

To prove the second part of this conclusion, first convince yourself that row operations and the matrix-vector are commutative operations. By this we mean the following. Suppose that F is an m × n matrix that is row-equivalent to the matrix G. Apply to the column vector Fw the same sequence of row operations that converts F to G. Then the result is Gw. So we can do row operations on the matrix, then do a matrix-vector product, or do a matrix-vector product and then do row operations on a column vector, and the result will be the same either way. Since matrix multiplication is defined by a collection of matrix-vector products (Definition MM), if we apply to the matrix product FH the same sequence of row operations that converts F to G then the result will equal GH. Now apply these observations to A.

Write A{I}_{n} = {I}_{m}A and apply the row operations that convert M to N. A is converted to B, while {I}_{m} is converted to J, so we have B{I}_{n} = JA. Simplifying the left side gives the desired conclusion.

For the third conclusion, we now establish the two equivalences

\eqalignno{ Ax & = y & &\kern 3.26288pt \mathrel{⇔}\kern 3.26288pt &JAx & = Jy & &\kern 3.26288pt \mathrel{⇔}\kern 3.26288pt &Bx & = Jy & & & & & & & & & & }

The forward direction of the first equivalence is accomplished by multiplying both sides of the matrix equality by J, while the backward direction is accomplished by multiplying by the inverse of J (which we know exists by Theorem NI since J is nonsingular). The second equivalence is obtained simply by the substitutions given by JA = B.

The first r rows of N are in reduced row-echelon form, since any contiguous collection of rows taken from a matrix in reduced row-echelon form will form a matrix that is again in reduced row-echelon form. Since the matrix C is formed by removing the last n entries of each these rows, the remainder is still in reduced row-echelon form. By its construction, C has no zero rows. C has r rows and each contains a leading 1, so there are r pivot columns in C.

The final m − r rows of N are in reduced row-echelon form, since any contiguous collection of rows taken from a matrix in reduced row-echelon form will form a matrix that is again in reduced row-echelon form. Since the matrix L is formed by removing the first n entries of each these rows, and these entries are all zero (they form the zero rows of B), the remainder is still in reduced row-echelon form. L is the final m − r rows of the nonsingular matrix J, so none of these rows can be totally zero, or J would not row-reduce to the identity matrix. L has m − r rows and each contains a leading 1, so there are m − r pivot columns in L.

Notice that in the case where A is a nonsingular matrix we know that the reduced row-echelon form of A is the identity matrix (Theorem NMRRI), so B = {I}_{n}. Then the second conclusion above says JA = B = {I}_{n}, so J is the inverse of A. Thus this theorem generalizes Theorem CINM, though the result is a “left-inverse” of A rather than a “right-inverse.”

The third conclusion of Theorem PEEF is the most telling. It says that x is a solution to the linear system ℒS\kern -1.95872pt \left (A,\kern 1.95872pt y\right ) if and only if x is a solution to the linear system ℒS\kern -1.95872pt \left (B,\kern 1.95872pt Jy\right ). Or said differently, if we row-reduce the augmented matrix \left [\left .A\kern 1.95872pt \right \vert \kern 1.95872pt y\right ] we will get the augmented matrix \left [\left .B\kern 1.95872pt \right \vert \kern 1.95872pt Jy\right ]. The matrix J tracks the cumulative effect of the row operations that converts A to reduced row-echelon form, here effectively applying them to the vector of constants in a system of equations having A as a coefficient matrix. When A row-reduces to a matrix with zero rows, then Jy should also have zero entries in the same rows if the system is to be consistent.

#### Subsection FS: Four Subsets

With all the preliminaries in place we can state our main result for this section. In essence this result will allow us to say that we can find linearly independent sets to use in span constructions for all four subsets (null space, column space, row space, left null space) by analyzing only the extended echelon form of the matrix, and specifically, just the two submatrices C and L, which will be ripe for analysis since they are already in reduced row-echelon form (Theorem PEEF).

Theorem FS
Four Subsets
Suppose A is an m × n matrix with extended echelon form N. Suppose the reduced row-echelon form of A has r nonzero rows. Then C is the submatrix of N formed from the first r rows and the first n columns and L is the submatrix of N formed from the last m columns and the last m − r rows. Then

1. The null space of A is the null space of C, N\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left (C\right ).
2. The row space of A is the row space of C, ℛ\kern -1.95872pt \left (A\right ) = ℛ\kern -1.95872pt \left (C\right ).
3. The column space of A is the null space of L, C\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left (L\right ).
4. The left null space of A is the row space of L, ℒ\kern -1.95872pt \left (A\right ) = ℛ\kern -1.95872pt \left (L\right ).

Proof   First, N\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left (B\right ) since B is row-equivalent to A (Theorem REMES). The zero rows of B represent equations that are always true in the homogeneous system ℒS\kern -1.95872pt \left (B,\kern 1.95872pt 0\right ), so the removal of these equations will not change the solution set. Thus, in turn, N\kern -1.95872pt \left (B\right ) = N\kern -1.95872pt \left (C\right ).

Second, ℛ\kern -1.95872pt \left (A\right ) = ℛ\kern -1.95872pt \left (B\right ) since B is row-equivalent to A (Theorem REMRS). The zero rows of B contribute nothing to the span that is the row space of B, so the removal of these rows will not change the row space. Thus, in turn, ℛ\kern -1.95872pt \left (B\right ) = ℛ\kern -1.95872pt \left (C\right ).

Third, we prove the set equality C\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left (L\right ) with Definition SE. Begin by showing that C\kern -1.95872pt \left (A\right ) ⊆N\kern -1.95872pt \left (L\right ). Choose y ∈C\kern -1.95872pt \left (A\right ) ⊆ {ℂ}^{m}. Then there exists a vector x ∈ {ℂ}^{n} such that Ax = y (Theorem CSCS). Then for 1 ≤ k ≤ m − r,

\eqalignno{ {\left [Ly\right ]}_{k} & ={ \left [Jy\right ]}_{r+k} & &\text{$L$ a submatrix of $J$} & & & & \cr & ={ \left [Bx\right ]}_{r+k} & &\text{@(a href="#theorem.PEEF")Theorem PEEF@(/a)} & & & & \cr & ={ \left [Ox\right ]}_{k} & &\text{Zero matrix a submatrix of $B$} & & & & \cr & ={ \left [0\right ]}_{k} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.MMZM")Theorem MMZM@(/a)} & & & & }

So, for all 1 ≤ k ≤ m − r, {\left [Ly\right ]}_{k} ={ \left [0\right ]}_{k}. So by Definition CVE we have Ly = 0 and thus y ∈N\kern -1.95872pt \left (L\right ).

Now, show that N\kern -1.95872pt \left (L\right ) ⊆C\kern -1.95872pt \left (A\right ). Choose y ∈N\kern -1.95872pt \left (L\right ) ⊆ {ℂ}^{m}. Form the vector Ky ∈ {ℂ}^{r}. The linear system ℒS\kern -1.95872pt \left (C,\kern 1.95872pt Ky\right ) is consistent since C is in reduced row-echelon form and has no zero rows (Theorem PEEF). Let x ∈ {ℂ}^{n} denote a solution to ℒS\kern -1.95872pt \left (C,\kern 1.95872pt Ky\right ).

Then for 1 ≤ j ≤ r,

\eqalignno{ {\left [Bx\right ]}_{j} & ={ \left [Cx\right ]}_{j} & &\text{$C$ a submatrix of $B$} & & & & \cr & ={ \left [Ky\right ]}_{j} & &\text{$x$ a solution to $ℒS\kern -1.95872pt \left (C,\kern 1.95872pt Ky\right )$} & & & & \cr & ={ \left [Jy\right ]}_{j} & &\text{$K$ a submatrix of $J$} & & & & \cr & & & & }

And for r + 1 ≤ k ≤ m,

\eqalignno{ {\left [Bx\right ]}_{k} & ={ \left [Ox\right ]}_{k−r} & &\text{Zero matrix a submatrix of $B$} & & & & \cr & ={ \left [0\right ]}_{k−r} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.MMZM")Theorem MMZM@(/a)} & & & & \cr & ={ \left [Ly\right ]}_{k−r} & &\text{$y$ in $N\kern -1.95872pt \left (L\right )$} & & & & \cr & ={ \left [Jy\right ]}_{k} & &\text{$L$ a submatrix of $J$} & & & & \cr & & & & }

So for all 1 ≤ i ≤ m, {\left [Bx\right ]}_{i} ={ \left [Jy\right ]}_{i} and by Definition CVE we have Bx = Jy. From Theorem PEEF we know then that Ax = y, and therefore y ∈C\kern -1.95872pt \left (A\right ) (Theorem CSCS). By Definition SE we now have C\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left (L\right ).

Fourth, we prove the set equality ℒ\kern -1.95872pt \left (A\right ) = ℛ\kern -1.95872pt \left (L\right ) with Definition SE. Begin by showing that ℛ\kern -1.95872pt \left (L\right ) ⊆ℒ\kern -1.95872pt \left (A\right ). Choose y ∈ℛ\kern -1.95872pt \left (L\right ) ⊆ {ℂ}^{m}. Then there exists a vector w ∈ {ℂ}^{m−r} such that y = {L}^{t}w (Definition RSM, Theorem CSCS). Then for 1 ≤ i ≤ n,

\eqalignno{ {\left [{A}^{t}y\right ]}_{ i} & ={ \mathop{∑ }}_{k=1}^{m}{\left [{A}^{t}\right ]}_{ ik}{\left [y\right ]}_{k} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.EMP")Theorem EMP@(/a)} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{m}{\left [{A}^{t}\right ]}_{ ik}{\left [{L}^{t}w\right ]}_{ k} & &\text{Definition of $w$} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{m}{\left [{A}^{t}\right ]}_{ ik}{ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [{L}^{t}\right ]}_{ kℓ}{\left [w\right ]}_{ℓ} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.EMP")Theorem EMP@(/a)} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{m}{ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [{A}^{t}\right ]}_{ ik}{\left [{L}^{t}\right ]}_{ kℓ}{\left [w\right ]}_{ℓ} & &\text{@(a href="fcla-jsmath-2.99li69.html#property.DCN")Property DCN@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m−r}{ \mathop{∑ }}_{k=1}^{m}{\left [{A}^{t}\right ]}_{ ik}{\left [{L}^{t}\right ]}_{ kℓ}{\left [w\right ]}_{ℓ} & &\text{@(a href="fcla-jsmath-2.99li69.html#property.CACN")Property CACN@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m−r}\left ({\mathop{∑ }}_{k=1}^{m}{\left [{A}^{t}\right ]}_{ ik}{\left [{L}^{t}\right ]}_{ kℓ}\right ){\left [w\right ]}_{ℓ} & &\text{@(a href="fcla-jsmath-2.99li69.html#property.DCN")Property DCN@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m−r}\left ({\mathop{∑ }}_{k=1}^{m}{\left [{A}^{t}\right ]}_{ ik}{\left [{J}^{t}\right ]}_{ k,r+ℓ}\right ){\left [w\right ]}_{ℓ} & &\text{$L$ a submatrix of $J$} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [{A}^{t}{J}^{t}\right ]}_{ i,r+ℓ}{\left [w\right ]}_{ℓ} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.EMP")Theorem EMP@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [{\left (JA\right )}^{t}\right ]}_{ i,r+ℓ}{\left [w\right ]}_{ℓ} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.MMT")Theorem MMT@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [{B}^{t}\right ]}_{ i,r+ℓ}{\left [w\right ]}_{ℓ} & &\text{@(a href="#theorem.PEEF")Theorem PEEF@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m−r}0{\left [w\right ]}_{ ℓ} & &\text{Zero rows in $B$} & & & & \cr & = 0 & &\text{@(a href="fcla-jsmath-2.99li69.html#property.ZCN")Property ZCN@(/a)} & & & & \cr & ={ \left [0\right ]}_{i} & &\text{@(a href="fcla-jsmath-2.99li18.html#definition.ZCV")Definition ZCV@(/a)} & & & & }

Since {\left [{A}^{t}y\right ]}_{ i} ={ \left [0\right ]}_{i} for 1 ≤ i ≤ n, Definition CVE implies that {A}^{t}y = 0. This means that y ∈N\kern -1.95872pt \left ({A}^{t}\right ).

Now, show that ℒ\kern -1.95872pt \left (A\right ) ⊆ℛ\kern -1.95872pt \left (L\right ). Choose y ∈ℒ\kern -1.95872pt \left (A\right ) ⊆ {ℂ}^{m}. The matrix J is nonsingular (Theorem PEEF), so {J}^{t} is also nonsingular (Theorem MIT) and therefore the linear system ℒS\kern -1.95872pt \left ({J}^{t},\kern 1.95872pt y\right ) has a unique solution. Denote this solution as x ∈ {ℂ}^{m}. We will need to work with two “halves” of x, which we will denote as z and w with formal definitions given by

\eqalignno{ {\left [z\right ]}_{j} & ={ \left [x\right ]}_{i} & &1 ≤ j ≤ r, & & &{\left [w\right ]}_{k} & ={ \left [x\right ]}_{r+k} & &1 ≤ k ≤ m − r & & & & & & & & & & }

Now, for 1 ≤ j ≤ r,

\eqalignno{ {\left [{C}^{t}z\right ]}_{ j} & ={ \mathop{∑ }}_{k=1}^{r}{\left [{C}^{t}\right ]}_{ jk}{\left [z\right ]}_{k} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.EMP")Theorem EMP@(/a)} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{r}{\left [{C}^{t}\right ]}_{ jk}{\left [z\right ]}_{k} +{ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [O\right ]}_{ jℓ}{\left [w\right ]}_{ℓ} & &\text{@(a href="fcla-jsmath-2.99li30.html#definition.ZM")Definition ZM@(/a)} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{r}{\left [{B}^{t}\right ]}_{ jk}{\left [z\right ]}_{k} +{ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [{B}^{t}\right ]}_{ j,r+ℓ}{\left [w\right ]}_{ℓ} & &\text{$C$, $O$ submatrices of $B$} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{r}{\left [{B}^{t}\right ]}_{ jk}{\left [x\right ]}_{k} +{ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [{B}^{t}\right ]}_{ j,r+ℓ}{\left [x\right ]}_{r+ℓ} & &\text{Definitions of $z$ and $w$} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{r}{\left [{B}^{t}\right ]}_{ jk}{\left [x\right ]}_{k} +{ \mathop{∑ }}_{k=r+1}^{m}{\left [{B}^{t}\right ]}_{ jk}{\left [x\right ]}_{k} & &\text{Re-index second sum} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{m}{\left [{B}^{t}\right ]}_{ jk}{\left [x\right ]}_{k} & &\text{Combine sums} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{m}{\left [{\left (JA\right )}^{t}\right ]}_{ jk}{\left [x\right ]}_{k} & &\text{@(a href="#theorem.PEEF")Theorem PEEF@(/a)} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{m}{\left [{A}^{t}{J}^{t}\right ]}_{ jk}{\left [x\right ]}_{k} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.MMT")Theorem MMT@(/a)} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{m}{ \mathop{∑ }}_{ℓ=1}^{m}{\left [{A}^{t}\right ]}_{ jℓ}{\left [{J}^{t}\right ]}_{ ℓk}{\left [x\right ]}_{k} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.EMP")Theorem EMP@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m}{ \mathop{∑ }}_{k=1}^{m}{\left [{A}^{t}\right ]}_{ jℓ}{\left [{J}^{t}\right ]}_{ ℓk}{\left [x\right ]}_{k} & &\text{@(a href="fcla-jsmath-2.99li69.html#property.CACN")Property CACN@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m}{\left [{A}^{t}\right ]}_{ jℓ}\left ({\mathop{∑ }}_{k=1}^{m}{\left [{J}^{t}\right ]}_{ ℓk}{\left [x\right ]}_{k}\right ) & &\text{@(a href="fcla-jsmath-2.99li69.html#property.DCN")Property DCN@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m}{\left [{A}^{t}\right ]}_{ jℓ}{\left [{J}^{t}x\right ]}_{ ℓ} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.EMP")Theorem EMP@(/a)} & & & & \cr & ={ \mathop{∑ }}_{ℓ=1}^{m}{\left [{A}^{t}\right ]}_{ jℓ}{\left [y\right ]}_{ℓ} & &\text{Definition of $x$} & & & & \cr & ={ \left [{A}^{t}y\right ]}_{ j} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.EMP")Theorem EMP@(/a)} & & & & \cr & ={ \left [0\right ]}_{j} & &\text{$y ∈ℒ\kern -1.95872pt \left (A\right )$} & & & & }

So, by Definition CVE, {C}^{t}z = 0 and the vector z gives us a linear combination of the columns of {C}^{t} that equals the zero vector. In other words, z gives a relation of linear dependence on the rows of C. However, the rows of C are a linearly independent set by Theorem BRS. According to Definition LICV we must conclude that the entries of z are all zero, i.e. z = 0.

Now, for 1 ≤ i ≤ m, we have

\eqalignno{ {\left [y\right ]}_{i} & ={ \left [{J}^{t}x\right ]}_{ i} & &\text{Definition of $x$} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{m}{\left [{J}^{t}\right ]}_{ ik}{\left [x\right ]}_{k} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.EMP")Theorem EMP@(/a)} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{r}{\left [{J}^{t}\right ]}_{ ik}{\left [x\right ]}_{k} +{ \mathop{∑ }}_{k=r+1}^{m}{\left [{J}^{t}\right ]}_{ ik}{\left [x\right ]}_{k} & &\text{Break apart sum} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{r}{\left [{J}^{t}\right ]}_{ ik}{\left [z\right ]}_{k} +{ \mathop{∑ }}_{k=r+1}^{m}{\left [{J}^{t}\right ]}_{ ik}{\left [w\right ]}_{k−r} & &\text{Definition of $z$ and $w$} & & & & \cr & ={ \mathop{∑ }}_{k=1}^{r}{\left [{J}^{t}\right ]}_{ ik}0 +{ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [{J}^{t}\right ]}_{ i,r+ℓ}{\left [w\right ]}_{ℓ} & &\text{$z = 0$, re-index} & & & & \cr & = 0 +{ \mathop{∑ }}_{ℓ=1}^{m−r}{\left [{L}^{t}\right ]}_{ i,ℓ}{\left [w\right ]}_{ℓ} & &\text{$L$ a submatrix of $J$} & & & & \cr & ={ \left [{L}^{t}w\right ]}_{ i} & &\text{@(a href="fcla-jsmath-2.99li31.html#theorem.EMP")Theorem EMP@(/a)} & & & & \cr & & & & }

So by Definition CVE, y = {L}^{t}w. The existence of w implies that y ∈ℛ\kern -1.95872pt \left (L\right ), and therefore ℒ\kern -1.95872pt \left (A\right ) ⊆ℛ\kern -1.95872pt \left (L\right ). So by Definition SE we have ℒ\kern -1.95872pt \left (A\right ) = ℛ\kern -1.95872pt \left (L\right ).

The first two conclusions of this theorem are nearly trivial. But they set up a pattern of results for C that is reflected in the latter two conclusions about L. In total, they tell us that we can compute all four subsets just by finding null spaces and row spaces. This theorem does not tell us exactly how to compute these subsets, but instead simply expresses them as null spaces and row spaces of matrices in reduced row-echelon form without any zero rows (C and L). A linearly independent set that spans the null space of a matrix in reduced row-echelon form can be found easily with Theorem BNS. It is an even easier matter to find a linearly independent set that spans the row space of a matrix in reduced row-echelon form with Theorem BRS, especially when there are no zero rows present. So an application of Theorem FS is typically followed by two applications each of Theorem BNS and Theorem BRS.

The situation when r = m deserves comment, since now the matrix L has no rows. What is C\kern -1.95872pt \left (A\right ) when we try to apply Theorem FS and encounter N\kern -1.95872pt \left (L\right )? One interpretation of this situation is that L is the coefficient matrix of a homogeneous system that has no equations. How hard is it to find a solution vector to this system? Some thought will convince you that any proposed vector will qualify as a solution, since it makes all of the equations true. So every possible vector is in the null space of L and therefore C\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left (L\right ) = {ℂ}^{m}. OK, perhaps this sounds like some twisted argument from Alice in Wonderland. Let us try another argument that might solidly convince you of this logic.

If r = m, when we row-reduce the augmented matrix of ℒS\kern -1.95872pt \left (A,\kern 1.95872pt b\right ) the result will have no zero rows, and all the leading 1’s will occur in first n columns, so by Theorem RCLS the system will be consistent. By Theorem CSCS, b ∈C\kern -1.95872pt \left (A\right ). Since b was arbitrary, every possible vector is in the column space of A, so we again have C\kern -1.95872pt \left (A\right ) = {ℂ}^{m}. The situation when a matrix has r = m is known by the term full rank, and in the case of a square matrix coincides with nonsingularity (see Exercise FS.M50).

The properties of the matrix L described by this theorem can be explained informally as follows. A column vector y ∈ {ℂ}^{m} is in the column space of A if the linear system ℒS\kern -1.95872pt \left (A,\kern 1.95872pt y\right ) is consistent (Theorem CSCS). By Theorem RCLS, the reduced row-echelon form of the augmented matrix \left [\left .A\kern 1.95872pt \right \vert \kern 1.95872pt y\right ] of a consistent system will have zeros in the bottom m − r locations of the last column. By Theorem PEEF this final column is the vector Jy and so should then have zeros in the final m − r locations. But since L comprises the final m − r rows of J, this condition is expressed by saying y ∈N\kern -1.95872pt \left (L\right ).

Additionally, the rows of J are the scalars in linear combinations of the rows of A that create the rows of B. That is, the rows of J record the net effect of the sequence of row operations that takes A to its reduced row-echelon form, B. This can be seen in the equation JA = B (Theorem PEEF). As such, the rows of L are scalars for linear combinations of the rows of A that yield zero rows. But such linear combinations are precisely the elements of the left null space. So any element of the row space of L is also an element of the left null space of A. We will now illustrate Theorem FS with a few examples.

Example FS1
Four subsets, #1
In Example SEEF we found the five relevant submatrices of the matrix

 A = \left [\array{ 1 &−1&−2& 7 & 1 & 6\cr −6 & 2 &−4 &−18 &−3 &−26 \cr 4 &−1& 4 & 10 & 2 & 17\cr 3 &−1 & 2 & 9 & 1 & 12 } \right ]

To apply Theorem FS we only need C and L,

\eqalignno{ C & = \left [\array{ \text{1}&0&2& 1 &0& 3\cr 0&\text{1 } &4 &−6 &0 &−1 \cr 0&0&0& 0 &\text{1}& 2 } \right ] &L & = \left [\array{ \text{1}&2&2&1 } \right ] & & & & }

Then we use Theorem FS to obtain

\eqalignno{ N\kern -1.95872pt \left (A\right ) & = N\kern -1.95872pt \left (C\right ) = \left \langle \left \{\left [\array{ −2\cr −4 \cr 1\cr 0 \cr 0\cr 0 } \right ],\kern 1.95872pt \left [\array{ −1\cr 6 \cr 0\cr 1 \cr 0\cr 0 } \right ],\kern 1.95872pt \left [\array{ −3\cr 1 \cr 0\cr 0 \cr −2\cr 1 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li26.html#theorem.BNS")Theorem BNS@(/a)} & & & & \cr ℛ\kern -1.95872pt \left (A\right ) & = ℛ\kern -1.95872pt \left (C\right ) = \left \langle \left \{\left [\array{ 1\cr 0 \cr 2\cr 1 \cr 0\cr 3 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr 4\cr −6 \cr 0\cr −1 } \right ],\kern 1.95872pt \left [\array{ 0\cr 0 \cr 0\cr 0 \cr 1\cr 2 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li34.html#theorem.BRS")Theorem BRS@(/a)} & & & & \cr C\kern -1.95872pt \left (A\right ) & = N\kern -1.95872pt \left (L\right ) = \left \langle \left \{\left [\array{ −2\cr 1 \cr 0\cr 0 } \right ],\kern 1.95872pt \left [\array{ −2\cr 0 \cr 1\cr 0 } \right ],\kern 1.95872pt \left [\array{ −1\cr 0 \cr 0\cr 1 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li26.html#theorem.BNS")Theorem BNS@(/a)} & & & & \cr ℒ\kern -1.95872pt \left (A\right ) & = ℛ\kern -1.95872pt \left (L\right ) = \left \langle \left \{\left [\array{ 1\cr 2 \cr 2\cr 1 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li34.html#theorem.BRS")Theorem BRS@(/a)} & & & & }

Boom!

Example FS2
Four subsets, #2
Now lets return to the matrix A that we used to motivate this section in Example CSANS,

 A = \left [\array{ 10 & 0 & 3 & 8 & 7\cr −16 &−1 &−4 &−10 &−13 \cr −6 & 1 &−3& −6 & −6\cr 0 & 2 &−2 & −3 & −2 \cr 3 & 0 & 1 & 2 & 3\cr −1 &−1 & 1 & 1 & 0 } \right ]

We form the matrix M by adjoining the 6 × 6 identity matrix {I}_{6},

 M = \left [\array{ 10 & 0 & 3 & 8 & 7 &1&0&0&0&0&0\cr −16 &−1 &−4 &−10 &−13 &0 &1 &0 &0 &0 &0 \cr −6 & 1 &−3& −6 & −6 &0&0&1&0&0&0\cr 0 & 2 &−2 & −3 & −2 &0 &0 &0 &1 &0 &0 \cr 3 & 0 & 1 & 2 & 3 &0&0&0&0&1&0\cr −1 &−1 & 1 & 1 & 0 &0 &0 &0 &0 &0 &1 } \right ]

and row-reduce to obtain N

 N = \left [\array{ \text{1}&0&0&0& 2 &0&0& 1 &−1& 2 &−1\cr 0&\text{1 } &0 &0 &−3 &0 &0 &−2 & 3 &−3 & 3 \cr 0&0&\text{1}&0& 1 &0&0& 1 & 1 & 3 & 3\cr 0&0 &0 &\text{1 } &−2 &0 &0 &−2 & 1 &−4 & 0 \cr 0&0&0&0& 0 &\text{1}&0& 3 &−1& 3 & 1\cr 0&0 &0 &0 & 0 &0 &\text{1 } &−2 & 1 & 1 &−1 } \right ]

To find the four subsets for A, we only need identify the 4 × 5 matrix C and the 2 × 6 matrix L,

\eqalignno{ C & = \left [\array{ \text{1}&0&0&0& 2\cr 0&\text{1 } &0 &0 &−3 \cr 0&0&\text{1}&0& 1\cr 0&0 &0 &\text{1 } &−2 } \right ] &L & = \left [\array{ \text{1}&0& 3 &−1&3& 1\cr 0&\text{1 } &−2 & 1 &1 &−1 } \right ] & & & & }

Then we apply Theorem FS,

\eqalignno{ N\kern -1.95872pt \left (A\right ) & = N\kern -1.95872pt \left (C\right ) = \left \langle \left \{\left [\array{ −2\cr 3 \cr −1\cr 2 \cr 1 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li26.html#theorem.BNS")Theorem BNS@(/a)} & & & & \cr ℛ\kern -1.95872pt \left (A\right ) & = ℛ\kern -1.95872pt \left (C\right ) = \left \langle \left \{\left [\array{ 1\cr 0 \cr 0\cr 0 \cr 2 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr 0\cr 0 \cr −3 } \right ],\kern 1.95872pt \left [\array{ 0\cr 0 \cr 1\cr 0 \cr 1 } \right ],\kern 1.95872pt \left [\array{ 0\cr 0 \cr 0\cr 1 \cr −2 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li34.html#theorem.BRS")Theorem BRS@(/a)} & & & & \cr C\kern -1.95872pt \left (A\right ) & = N\kern -1.95872pt \left (L\right ) = \left \langle \left \{\left [\array{ −3\cr 2 \cr 1\cr 0 \cr 0\cr 0 } \right ],\kern 1.95872pt \left [\array{ 1\cr −1 \cr 0\cr 1 \cr 0\cr 0 } \right ],\kern 1.95872pt \left [\array{ −3\cr −1 \cr 0\cr 0 \cr 1\cr 0 } \right ],\kern 1.95872pt \left [\array{ −1\cr 1 \cr 0\cr 0 \cr 0\cr 1 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li26.html#theorem.BNS")Theorem BNS@(/a)} & & & & \cr ℒ\kern -1.95872pt \left (A\right ) & = ℛ\kern -1.95872pt \left (L\right ) = \left \langle \left \{\left [\array{ 1\cr 0 \cr 3\cr −1 \cr 3\cr 1 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr −2\cr 1 \cr 1\cr −1 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li34.html#theorem.BRS")Theorem BRS@(/a)} & & & & }

The next example is just a bit different since the matrix has more rows than columns, and a trivial null space.

Example FSAG
Four subsets, Archetype G
Archetype G and Archetype H are both systems of m = 5 equations in n = 2 variables. They have identical coefficient matrices, which we will denote here as the matrix G,

 G = \left [\array{ 2 & 3\cr −1 & 4 \cr 3 &10\cr 3 &−1 \cr 6 & 9 } \right ]

Adjoin the 5 × 5 identity matrix, {I}_{5}, to form

 M = \left [\array{ 2 & 3 &1&0&0&0&0\cr −1 & 4 &0 &1 &0 &0 &0 \cr 3 &10&0&0&1&0&0\cr 3 &−1 &0 &0 &0 &1 &0 \cr 6 & 9 &0&0&0&0&1 } \right ]

This row-reduces to

 N = \left [\array{ \text{1}&0&0&0&0& {3\over _ 11} & {1\over _ 33} \cr 0&\text{1}&0&0&0&−{2\over _ 11}& {1\over _ 11} \cr 0&0&\text{1}&0&0& 0 &−{1\over 3} \cr 0&0&0&\text{1}&0& 1 &−{1\over 3} \cr 0&0&0&0&\text{1}& 1 &−1 } \right ]

The first n = 2 columns contain r = 2 leading 1’s, so we obtain C as the 2 × 2 identity matrix and extract L from the final m − r = 3 rows in the final m = 5 columns.

\eqalignno{ C & = \left [\array{ \text{1}&0\cr 0&\text{1} } \right ] &L & = \left [\array{ \text{1}&0&0&0&−{1\over 3} \cr 0&\text{1}&0&1&−{1\over 3} \cr 0&0&\text{1}&1&−1 } \right ] & & & & }

Then we apply Theorem FS,

\eqalignno{ N\kern -1.95872pt \left (G\right ) = N\kern -1.95872pt \left (C\right ) & = \left \langle ∅\right \rangle = \left \{0\right \} & &\text{@(a href="fcla-jsmath-2.99li26.html#theorem.BNS")Theorem BNS@(/a)} & & & & \cr ℛ\kern -1.95872pt \left (G\right ) = ℛ\kern -1.95872pt \left (C\right ) & = \left \langle \left \{\left [\array{ 1\cr 0 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 } \right ]\right \}\right \rangle = {ℂ}^{2} & &\text{@(a href="fcla-jsmath-2.99li34.html#theorem.BRS")Theorem BRS@(/a)} & & & & \cr C\kern -1.95872pt \left (G\right ) = N\kern -1.95872pt \left (L\right ) & = \left \langle \left \{\left [\array{ 0\cr −1 \cr −1\cr 1 \cr 0 } \right ],\kern 1.95872pt \left [\array{ {1\over 3} \cr {1\over 3} \cr 1\cr 0 \cr 1 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li26.html#theorem.BNS")Theorem BNS@(/a)} & & & & \cr & = \left \langle \left \{\left [\array{ 0\cr −1 \cr −1\cr 1 \cr 0 } \right ],\kern 1.95872pt \left [\array{ 1\cr 1 \cr 3\cr 0 \cr 3 } \right ]\right \}\right \rangle & & & & \cr ℒ\kern -1.95872pt \left (G\right ) = ℛ\kern -1.95872pt \left (L\right ) & = \left \langle \left \{\left [\array{ 1\cr 0 \cr 0\cr 0 \cr −{1\over 3} } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr 0\cr 1 \cr −{1\over 3} } \right ],\kern 1.95872pt \left [\array{ 0\cr 0 \cr 1\cr 1 \cr −1 } \right ]\right \}\right \rangle & &\text{@(a href="fcla-jsmath-2.99li34.html#theorem.BRS")Theorem BRS@(/a)} & & & & \cr & = \left \langle \left \{\left [\array{ 3\cr 0 \cr 0\cr 0 \cr −1 } \right ],\kern 1.95872pt \left [\array{ 0\cr 3 \cr 0\cr 3\ \cr −1 } \right ],\kern 1.95872pt \left [\array{ 0\cr 0 \cr 1\cr 1 \cr −1 } \right ]\right \}\right \rangle & & & & }

As mentioned earlier, Archetype G is consistent, while Archetype H is inconsistent. See if you can write the two different vectors of constants from these two archetypes as linear combinations of the two vectors in C\kern -1.95872pt \left (G\right ). How about the two columns of G, can you write each individually as a linear combination of the two vectors in C\kern -1.95872pt \left (G\right )? They must be in the column space of G also. Are your answers unique? Do you notice anything about the scalars that appear in the linear combinations you are forming?

Example COV and Example CSROI each describes the column space of the coefficient matrix from Archetype I as the span of a set of r = 3 linearly independent vectors. It is no accident that these two different sets both have the same size. If we (you?) were to calculate the column space of this matrix using the null space of the matrix L from Theorem FS then we would again find a set of 3 linearly independent vectors that span the range. More on this later.

So we have three different methods to obtain a description of the column space of a matrix as the span of a linearly independent set. Theorem BCS is sometimes useful since the vectors it specifies are equal to actual columns of the matrix. Theorem BRS and Theorem CSRST combine to create vectors with lots of zeros, and strategically placed 1’s near the top of the vector. Theorem FS and the matrix L from the extended echelon form gives us a third method, which tends to create vectors with lots of zeros, and strategically placed 1’s near the bottom of the vector. If we don’t care about linear independence we can also appeal to Definition CSM and simply express the column space as the span of all the columns of the matrix, giving us a fourth description.

With Theorem CSRST and Definition RSM, we can compute column spaces with theorems about row spaces, and we can compute row spaces with theorems about row spaces, but in each case we must transpose the matrix first. At this point you may be overwhelmed by all the possibilities for computing column and row spaces. Diagram CSRST is meant to help. For both the column space and row space, it suggests four techniques. One is to appeal to the definition, another yields a span of a linearly independent set, and a third uses Theorem FS. A fourth suggests transposing the matrix and the dashed line implies that then the companion set of techniques can be applied. This can lead to a bit of silliness, since if you were to follow the dashed lines twice you would transpose the matrix twice, and by Theorem TT would accomplish nothing productive.

Although we have many ways to describe a column space, notice that one tempting strategy will usually fail. It is not possible to simply row-reduce a matrix directly and then use the columns of the row-reduced matrix as a set whose span equals the column space. In other words, row operations do not preserve column spaces (however row operations do preserve row spaces, Theorem REMRS). See Exercise CRS.M21.

1. Find a nontrivial element of the left null space of A.
 A = \left [\array{ 2 & 1 &−3& 4\cr −1 &−1 & 2 &−1 \cr 0 &−1& 1 & 2} \right ]
2. Find the matrices C and L in the extended echelon form of A.
 A = \left [\array{ −9& 5 &−3\cr 2 &−1 & 1 \cr −5& 3 &−1} \right ]
3. Why is Theorem FS a great conclusion to Chapter M?

#### Subsection EXC: Exercises

C20 Example FSAG concludes with several questions. Perform the analysis suggested by these questions.
Contributed by Robert Beezer

C25 Given the matrix A below, use the extended echelon form of A to answer each part of this problem. In each part, find a linearly independent set of vectors, S, so that the span of S, \left \langle S\right \rangle , equals the specified set of vectors.

 A = \left [\array{ −5& 3 &−1\cr −1 & 1 & 1 \cr −8& 5 &−1\cr 3 &−2 & 0 } \right ]

(a) The row space of A, ℛ\kern -1.95872pt \left (A\right ).
(b) The column space of A, C\kern -1.95872pt \left (A\right ).
(c) The null space of A, N\kern -1.95872pt \left (A\right ).
(d) The left null space of A, ℒ\kern -1.95872pt \left (A\right ).

Contributed by Robert Beezer Solution [831]

C26 For the matrix D below use the extended echelon form to find
(a) a linearly independent set whose span is the column space of D.
(b) a linearly independent set whose span is the left null space of D.

\eqalignno{ D & = \left [\array{ −7&−11&−19&−15\cr 6 & 10 & 18 & 14 \cr 3 & 5 & 9 & 7\cr −1 & −2 & −4 & −3 } \right ] & & }

Contributed by Robert Beezer Solution [833]

C41 The following archetypes are systems of equations. For each system, write the vector of constants as a linear combination of the vectors in the span construction for the column space provided by Theorem FS and Theorem BNS (these vectors are listed for each of these archetypes).
Archetype A
Archetype B
Archetype C
Archetype D
Archetype E
Archetype F
Archetype G
Archetype H
Archetype I
Archetype J

Contributed by Robert Beezer

C43 The following archetypes are either matrices or systems of equations with coefficient matrices. For each matrix, compute the extended echelon form N and identify the matrices C and L. Using Theorem FS, Theorem BNS and Theorem BRS express the null space, the row space, the column space and left null space of each coefficient matrix as a span of a linearly independent set.
Archetype A
Archetype B
Archetype C
Archetype D/Archetype E
Archetype F
Archetype G/Archetype H
Archetype I
Archetype J
Archetype K
Archetype L

Contributed by Robert Beezer

C60 For the matrix B below, find sets of vectors whose span equals the column space of B (C\kern -1.95872pt \left (B\right )) and which individually meet the following extra requirements.
(a) The set illustrates the definition of the column space.
(b) The set is linearly independent and the members of the set are columns of B.
(c) The set is linearly independent with a “nice pattern of zeros and ones” at the top of each vector.
(d) The set is linearly independent with a “nice pattern of zeros and ones” at the bottom of each vector.

 B = \left [\array{ 2 &3&1& 1\cr 1 &1 &0 & 1 \cr −1&2&3&−4 } \right ]

Contributed by Robert Beezer Solution [835]

C61 Let A be the matrix below, and find the indicated sets with the requested properties.

 A = \left [\array{ 2 &−1& 5 &−3\cr −5 & 3 &−12 & 7 \cr 1 & 1 & 4 &−3 } \right ]

(a) A linearly independent set S so that C\kern -1.95872pt \left (A\right ) = \left \langle S\right \rangle and S is composed of columns of A.
(b) A linearly independent set S so that C\kern -1.95872pt \left (A\right ) = \left \langle S\right \rangle and the vectors in S have a nice pattern of zeros and ones at the top of the vectors.
(c) A linearly independent set S so that C\kern -1.95872pt \left (A\right ) = \left \langle S\right \rangle and the vectors in S have a nice pattern of zeros and ones at the bottom of the vectors.
(d) A linearly independent set S so that ℛ\kern -1.95872pt \left (A\right ) = \left \langle S\right \rangle .
Contributed by Robert Beezer Solution [838]

M50 Suppose that A is a nonsingular matrix. Extend the four conclusions of Theorem FS in this special case and discuss connections with previous results (such as Theorem NME4).
Contributed by Robert Beezer

M51 Suppose that A is a singular matrix. Extend the four conclusions of Theorem FS in this special case and discuss connections with previous results (such as Theorem NME4).
Contributed by Robert Beezer

#### Subsection SOL: Solutions

C25 Contributed by Robert Beezer Statement [826]
Add a 4 × 4 identity matrix to the right of A to form the matrix M and then row-reduce to the matrix N,

 M = \left [\array{ −5& 3 &−1&1&0&0&0\cr −1 & 1 & 1 &0 &1 &0 &0 \cr −8& 5 &−1&0&0&1&0\cr 3 &−2 & 0 &0 &0 &0 &1 } \right ]\mathop{\longrightarrow}\limits_{}^{\text{RREF}}\left [\array{ \text{1}&0&2&0&0&−2&−5\cr 0&\text{1 } &3 &0 &0 &−3 &−8 \cr 0&0&0&\text{1}&0&−1&−1\cr 0&0 &0 &0 &\text{1 } & 1 & 3 } \right ] = N

To apply Theorem FS in each of these four parts, we need the two matrices,

\eqalignno{ C & = \left [\array{ \text{1}&0&2\cr 0&\text{1 } &3} \right ] &L & = \left [\array{ \text{1}&0&−1&−1\cr 0&\text{1 } & 1 & 3 } \right ] & & & & }

(a)

\eqalignno{ ℛ\kern -1.95872pt \left (A\right ) & = ℛ\kern -1.95872pt \left (C\right ) & &\text{@(a href="#theorem.FS")Theorem FS@(/a)} & & & & \cr & = \left \langle \left [\array{ 1\cr 0 \cr 2 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr 3 } \right ]\right \rangle & &\text{@(a href="fcla-jsmath-2.99li34.html#theorem.BRS")Theorem BRS@(/a)} & & & & }

(b)

\eqalignno{ C\kern -1.95872pt \left (A\right ) & = N\kern -1.95872pt \left (L\right ) & &\text{@(a href="#theorem.FS")Theorem FS@(/a)} & & & & \cr & = \left \langle \left [\array{ 1\cr −1 \cr 1\cr 0 } \right ],\kern 1.95872pt \left [\array{ 1\cr −3 \cr 0\cr 1 } \right ]\right \rangle & &\text{@(a href="fcla-jsmath-2.99li26.html#theorem.BNS")Theorem BNS@(/a)} & & & & }

(c)

\eqalignno{ N\kern -1.95872pt \left (A\right ) & = N\kern -1.95872pt \left (C\right ) & &\text{@(a href="#theorem.FS")Theorem FS@(/a)} & & & & \cr & = \left \langle \left [\array{ −2\cr −3 \cr 1 } \right ]\right \rangle & &\text{@(a href="fcla-jsmath-2.99li26.html#theorem.BNS")Theorem BNS@(/a)} & & & & }

(d)

\eqalignno{ ℒ\kern -1.95872pt \left (A\right ) & = ℛ\kern -1.95872pt \left (L\right ) & &\text{@(a href="#theorem.FS")Theorem FS@(/a)} & & & & \cr & = \left \langle \left [\array{ 1\cr 0 \cr −1\cr −1 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr 1\cr 3 } \right ]\right \rangle & &\text{@(a href="fcla-jsmath-2.99li34.html#theorem.BRS")Theorem BRS@(/a)} & & & & }

C26 Contributed by Robert Beezer Statement [826]
For both parts, we need the extended echelon form of the matrix.

\eqalignno{ \left [\array{ −7&−11&−19&−15&1&0&0&0\cr 6 & 10 & 18 & 14 &0 &1 &0 &0 \cr 3 & 5 & 9 & 7 &0&0&1&0\cr −1 & −2 & −4 & −3 &0 &0 &0 &1 } \right ]\mathop{\longrightarrow}\limits_{}^{\text{RREF}}\left [\array{ \text{1}&0&−2&−1&0&0& 2 & 5\cr 0&\text{1 } & 3 & 2 &0 &0 &−1 &−3 \cr 0&0& 0 & 0 &\text{1}&0& 3 & 2\cr 0&0 & 0 & 0 &0 &\text{1 } &−2 & 0 } \right ]&& }

From this matrix we extract the last two rows, in the last four columns to form the matrix L,

\eqalignno{ L = \left [\array{ \text{1}&0& 3 &2\cr 0&\text{1 } &−2 &0 } \right ] & & }

(a) By Theorem FS and Theorem BNS we have

\eqalignno{ C\kern -1.95872pt \left (D\right ) = N\kern -1.95872pt \left (L\right ) = \left \langle \left \{\left [\array{ −3\cr 2 \cr 1\cr 0 } \right ],\kern 1.95872pt \left [\array{ −2\cr 0 \cr 0\cr 1 } \right ]\right \}\right \rangle & & }

(b) By Theorem FS and Theorem BRS we have

\eqalignno{ ℒ\kern -1.95872pt \left (D\right ) = ℛ\kern -1.95872pt \left (L\right ) = \left \langle \left \{\left [\array{ 1\cr 0 \cr 3\cr 2 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr −2\cr 0 } \right ]\right \}\right \rangle & & }

C60 Contributed by Robert Beezer Statement [828]
(a) The definition of the column space is the span of the set of columns (Definition CSM). So the desired set is just the four columns of B,

 S = \left \{\left [\array{ 2\cr 1 \cr −1 } \right ],\kern 1.95872pt \left [\array{ 3\cr 1 \cr 2 } \right ],\kern 1.95872pt \left [\array{ 1\cr 0 \cr 3 } \right ],\kern 1.95872pt \left [\array{ 1\cr 1 \cr −4 } \right ]\right \}

(b) Theorem BCS suggests row-reducing the matrix and using the columns of B that correspond to the pivot columns.

 B\mathop{\longrightarrow}\limits_{}^{\text{RREF}}\left [\array{ \text{1}&0&−1& 2\cr 0&\text{1 } & 1 &−1 \cr 0&0& 0 & 0 } \right ]

So the pivot columns are numbered by elements of D = \left \{1,\kern 1.95872pt 2\right \}, so the requested set is

 S = \left \{\left [\array{ 2\cr 1 \cr −1 } \right ],\kern 1.95872pt \left [\array{ 3\cr 1 \cr 2 } \right ]\right \}

(c) We can find this set by row-reducing the transpose of B, deleting the zero rows, and using the nonzero rows as column vectors in the set. This is an application of Theorem CSRST followed by Theorem BRS.

 { B}^{t}\mathop{\longrightarrow}\limits_{}^{\text{RREF}}\left [\array{ \text{1}&0& 3\cr 0&\text{1 } &−7 \cr 0&0& 0\cr 0&0 & 0 } \right ]

So the requested set is

 S = \left \{\left [\array{ 1\cr 0 \cr 3 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr −7 } \right ]\right \}

(d) With the column space expressed as a null space, the vectors obtained via Theorem BNS will be of the desired shape. So we first proceed with Theorem FS and create the extended echelon form,

 \left [\left .B\kern 1.95872pt \right \vert \kern 1.95872pt {I}_{3}\right ] = \left [\array{ 2 &3&1& 1 &1&0&0\cr 1 &1 &0 & 1 &0 &1 &0 \cr −1&2&3&−4&0&0&1 } \right ]\mathop{\longrightarrow}\limits_{}^{\text{RREF}}\left [\array{ \text{1}&0&−1& 2 &0& {2\over 3} &{−1\over 3} \cr 0&\text{1}& 1 &−1&0& {1\over 3} & {1\over 3} \cr 0&0& 0 & 0 &\text{1}&{−7\over 3} &{−1\over 3} } \right ]

So, employing Theorem FS, we have C\kern -1.95872pt \left (B\right ) = N\kern -1.95872pt \left (L\right ), where

 L = \left [\array{ \text{1}&{−7\over 3} &{−1\over 3} } \right ]

We can find the desired set of vectors from Theorem BNS as

 S = \left \{\left [\array{ {7\over 3} \cr 1\cr 0 } \right ],\kern 1.95872pt \left [\array{ {1\over 3} \cr 0\cr 1 } \right ]\right \}

C61 Contributed by Robert Beezer Statement [829]
(a) First find a matrix B that is row-equivalent to A and in reduced row-echelon form

 B = \left [\array{ \text{1}&0&3&−2\cr 0&\text{1 } &1 &−1 \cr 0&0&0& 0 } \right ]

By Theorem BCS we can choose the columns of A that correspond to dependent variables (D = \left \{1, 2\right \}) as the elements of S and obtain the desired properties. So

 S = \left \{\left [\array{ 2\cr −5 \cr 1 } \right ],\kern 1.95872pt \left [\array{ −1\cr 3 \cr 1 } \right ]\right \}

(b) We can write the column space of A as the row space of the transpose (Theorem CSRST). So we row-reduce the transpose of A to obtain the row-equivalent matrix C in reduced row-echelon form

 C = \left [\array{ 1&0&8\cr 0&1 &3 \cr 0&0&0\cr 0&0 &0} \right ]

The nonzero rows (written as columns) will be a linearly independent set that spans the row space of {A}^{t}, by Theorem BRS, and the zeros and ones will be at the top of the vectors,

 S = \left \{\left [\array{ 1\cr 0 \cr 8 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr 3 } \right ]\right \}

(c) In preparation for Theorem FS, augment A with the 3 × 3 identity matrix {I}_{3} and row-reduce to obtain the extended echelon form,

 \left [\array{ 1&0&3&−2&0&−{1\over 8}& {3\over 8} \cr 0&1&1&−1&0& {1\over 8} & {5\over 8} \cr 0&0&0& 0 &1& {3\over 8} &−{1\over 8} } \right ]

Then since the first four columns of row 3 are all zeros, we extract

 L = \left [\array{ \text{1}&{3\over 8}&−{1\over 8} } \right ]

Theorem FS says that C\kern -1.95872pt \left (A\right ) = N\kern -1.95872pt \left (L\right ). We can then use Theorem BNS to construct the desired set S, based on the free variables with indices in F = \left \{2, 3\right \} for the homogeneous system ℒS\kern -1.95872pt \left (L,\kern 1.95872pt 0\right ), so

 S = \left \{\left [\array{ −{3\over 8} \cr 1\cr 0 } \right ],\kern 1.95872pt \left [\array{ {1\over 8} \cr 0\cr 1 } \right ]\right \}

Notice that the zeros and ones are at the bottom of the vectors.
(d) This is a straightforward application of Theorem BRS. Use the row-reduced matrix B from part (a), grab the nonzero rows, and write them as column vectors,

 S = \left \{\left [\array{ 1\cr 0 \cr 3\cr −2 } \right ],\kern 1.95872pt \left [\array{ 0\cr 1 \cr 1\cr −1 } \right ]\right \}