7 Vector Spaces

7.1 Vector Spaces

A vector space is a nonempty set VV of vectors, and is closed under vector addition and scalar multiplication. Specifically, these operations must follow these 10 axioms, for all vectors 𝒖\mathbf{u}, 𝒗\mathbf{v}, 𝒘\mathbf{w} and scalars cc, dd:
Addition axioms:

  1. 𝒖+𝒗\mathbf{u}+\mathbf{v}, must be in VV.

  2. 𝒖+𝒗=𝒗+𝒖\mathbf{u}+\mathbf{v}=\mathbf{v}+\mathbf{u}.

  3. (𝒖+𝒗)+𝒘=𝒖+(𝒗+𝒘)(\mathbf{u}+\mathbf{v})+\mathbf{w}=\mathbf{u}+(\mathbf{v}+\mathbf{w}).

  4. There is a zero vector 𝟎\mathbf{0} in VV such that 𝒖+𝟎=𝒖\mathbf{u}+\mathbf{0}=\mathbf{u}.

  5. For each 𝒖\mathbf{u} in VV, there is a vector −𝒖-\mathbf{u} in VV such that 𝒖+(−𝒖)=𝟎\mathbf{u}+(-\mathbf{u})=\mathbf{0}.

Scalar multiplication axioms:

  1. The scalar multiple of uu by cc, denoted by c𝒖c\mathbf{u}, must be in VV.

  2. c(𝒖+𝒗)=c𝒖+c𝒗c(\mathbf{u}+\mathbf{v})=c\mathbf{u}+c\mathbf{v}.

  3. c(𝒖+𝒗)=c𝒖+c𝒗c(\mathbf{u}+\mathbf{v})=c\mathbf{u}+c\mathbf{v}.(c+d)𝒖=c𝒖+d𝒖(c+d)\mathbf{u}=c\mathbf{u}+d\mathbf{u}.

  4. c(d𝒖)=(cd)𝒖c(d\mathbf{u})=(cd)\mathbf{u}.

  5. 1𝒖=𝒖1\mathbf{u}=\mathbf{u}.

Following from the above axioms,

An element of a vector space can be thought of as a single "point" or vector in that space.
Examples of vector spaces

7.2 Subspaces

A subspace is a subset of a vector space that is itself a vector space, i.e. it satisfies the 3 criteria:

  1. A subspace must contain the zero vector of the parent space.

  2. For any two vectors 𝒖\mathbf{u} and 𝒗\mathbf{v} that are in the subspace, 𝒖+𝒗\mathbf{u}+\mathbf{v} must also be in the subspace.

  3. For any vector 𝒖\mathbf{u} in in the subspace and any scalar cc, c𝒖c\mathbf{u} must also be in the subspace.

Examples of subspaces:

7.3 Subspace spanned by a set

Span: The span of a set of vectors 𝒗1,⋅⋅⋅,𝒗p\mathbf{v}_{1},\cdot\cdot\cdot,\mathbf{v}_{p}, denoted as Span{𝒗1,⋅⋅⋅,𝒗p}Span\{\mathbf{v}_{1},\cdot\cdot\cdot,\mathbf{v}_{p}\}, is the set of all possible vectors that can be created using linear combinations of 𝒗1,⋅⋅⋅,𝒗p\mathbf{v}_{1},\cdot\cdot\cdot,\mathbf{v}_{p}.

7.3.0.1 Theorem

If vectors 𝒗1,⋅⋅⋅,𝒗p\mathbf{v}_{1},\cdot\cdot\cdot,\mathbf{v}_{p} exist in a vector space VV, then Span{𝒗1,⋅⋅⋅,𝒗p}Span\{\mathbf{v}_{1},\cdot\cdot\cdot,\mathbf{v}_{p}\} is a valid subspace of VV.
Terminology:

7.3.0.2 Example

Given 𝒗1\mathbf{v}_1 and 𝒗2\mathbf{v}_2 in a vector space VV, let H=Span{𝒗1,𝒗2}H = Span\{\mathbf{v}_1, \mathbf{v}_2\}. Show that HH is a subspace of VV.
Prove the 3 subspace properties hold for HH:

  1. Contains the Zero Vector: Set the scalar weights to zero, i.e. 0𝒗𝟏+0𝒗𝟐=𝟎0\mathbf{v_1} + 0\mathbf{v_2} = \mathbf{0}.

  2. Closed under Vector Addition: Let 𝒖\mathbf{u} and 𝒘\mathbf{w} be any two vectors in HH. By definition, they can be written as linear combinations: 𝒖=c1𝒗1+c2𝒗2\mathbf{u} = c_1\mathbf{v}_1 + c_2\mathbf{v}_2 and 𝒘=d1𝒗1+d2𝒗2\mathbf{w} = d_1\mathbf{v}_1 + d_2\mathbf{v}_2. Adding them together gives: 𝒖+𝒘=(c1+d1)𝒗1+(c2+d2)𝒗2\mathbf{u} + \mathbf{w} = (c_1 + d_1)\mathbf{v}_1 + (c_2 + d_2)\mathbf{v}_2

  3. Closed under Scalar Multiplication: Let 𝒖=c1𝒗1+c2𝒗2\mathbf{u} = c_1\mathbf{v}_1 + c_2\mathbf{v}_2 be in HH, and let kk be any real scalar. Multiplying 𝒖\mathbf{u} by kk gives: k𝒖=(kc1)𝒗1+(kc2)𝒗2k\mathbf{u} = (kc_1)\mathbf{v}_1 + (kc_2)\mathbf{v}_2

Hence, HH qualifies as a subspace.

7.3.0.3 Example

Let HH be the set of all vectors of the form (a−3b,b−a,a,b)(a-3b, b-a, a, b) where aa and bb are real numbers. Show that HH is a subspace of ℝ4\mathbb{R}^4.
Expressing HH as linear combinations of 𝒗𝟏\mathbf{v_1} and 𝒗𝟐\mathbf{v_2}, H=[a−3b−a+bab]=a[1−110]+b[−3101]=a𝒗𝟏+b𝒗𝟐H = \begin{bmatrix} a - 3b \\ -a + b \\ a \\ b \end{bmatrix} = a \begin{bmatrix} 1 \\ -1 \\ 1 \\ 0 \end{bmatrix} + b \begin{bmatrix} -3 \\ 1 \\ 0 \\ 1 \end{bmatrix} = a\mathbf{v_1} + b\mathbf{v_2} Since 𝒗𝟏\mathbf{v_1}, 𝒗𝟐\mathbf{v_2} are in ℝ4\mathbb{R}^4 and H=Span{𝒗1,𝒗2}H = Span\{\mathbf{v}_1, \mathbf{v}_2\}, HH is a valid subspace of ℝ4\mathbb{R}^4.

7.4 Null space of matrix

The null space of an m×nm \times n matrix AA is the set of all vectors 𝒙\mathbf{x} in ℝn\mathbb{R}^n that satisfy the homogenous equation A𝒙=𝟎A\mathbf{x}=\mathbf{0}. 𝑵(A)={𝒙:𝒙∈ℝn,A𝒙=𝟎}\mathbf{N}(A) = \{\mathbf{x}: \mathbf{x} \in \mathbb{R}^n,\: A\mathbf{x} = \mathbf{0}\} The null space consists of all vectors 𝒙\mathbf{x} in ℝn\mathbb{R}^n that are mapped to the zero vector in ℝm\mathbb{R}^m by the linear transformation 𝒙↦A𝒙\mathbf{x} \mapsto A\mathbf{x}.

7.4.0.1 Theorem

The null space of an m×nm \times n matrix AA is a subspace of ℝn\mathbb{R}^n. Equivalently, the set of all solutions to the system of mm homogeneous linear equations in nn unknowns is a subspace of ℝn\mathbb{R}^n.
Proof. We show that 𝑵(A)\mathbf{N}(A) satisfies the 3 properties of a subspace.

  1. The zero vector is in 𝑵(A)\mathbf{N}(A). A(𝟎)=𝟎A(\mathbf{0})=\mathbf{0} Any matrix multiplied by the zero vector results in the zero vector, so the zero vector is in 𝑵(A)\mathbf{N}(A).

  2. 𝑵(A)\mathbf{N}(A) is closed under vector addition.
    To show to that if 𝒖\mathbf{u} and 𝒗\mathbf{v} are in 𝑵(A)\mathbf{N}(A), then 𝒖+𝒗\mathbf{u} + \mathbf{v} is also in 𝑵(A)\mathbf{N}(A). We have A𝒖=𝟎,A𝒗=𝟎A\mathbf{u}=\mathbf{0},\quad A\mathbf{v}=\mathbf{0} Using the distributive property of matrix multiplication, A(𝒖+𝒗)=A𝒖+A𝒗=𝟎+𝟎=𝟎A(\mathbf{u} + \mathbf{v}) = A\mathbf{u} + A\mathbf{v} = \mathbf{0} + \mathbf{0} = \mathbf{0} Hence 𝒖+𝒗∈𝑵(A)\mathbf{u} + \mathbf{v} \in \mathbf{N}(A).

  3. N(A)N(A) is closed under scalar multiplication.
    To show that if a vector uu is in N(A)N(A), then multiplying it by any real scalar cc results in a vector c𝒖c\mathbf{u} that is also in N(A)N(A). We have A𝒖=𝟎A\mathbf{u} = \mathbf{0} Using properties of matrix scalar multiplication, A(c𝒖)=c(A𝒖)=c(𝟎)=𝟎A(c\mathbf{u}) = c(A\mathbf{u}) = c(\mathbf{0}) = \mathbf{0} Hence c𝒖∈𝑵(A)c\mathbf{u} \in \mathbf{N}(A).

7.4.0.2 Example

For this matrix AA and vector 𝒖\mathbf{u}, A=[1−3−2−591],𝒖=[53−2]A = \begin{bmatrix} 1 & -3 & -2 \\ -5 & 9 & 1 \end{bmatrix}, \quad \mathbf{u} = \begin{bmatrix} 5 \\ 3 \\ -2 \end{bmatrix} A𝒖=[1−3−2−591][53−2]=[5−9+4−25+27−2]=[00]A \mathbf{u} = \begin{bmatrix} 1 & -3 & -2 \\ -5 & 9 & 1 \end{bmatrix} \begin{bmatrix} 5 \\ 3 \\ -2 \end{bmatrix} = \begin{bmatrix} 5 - 9 + 4 \\ -25 + 27 - 2 \end{bmatrix} = \begin{bmatrix} 0 \\ 0 \end{bmatrix} So 𝒖∈𝑵(A)\mathbf{u} \in \mathbf{N}(A).

7.4.0.3 Example

Find the null space of A=[1236]A = \begin{bmatrix} 1 & 2 \\ 3 & 6 \end{bmatrix}.
Performing row operations on augmented matrix to reach RREF, [12|036|0]⟹[12|000|0]\begin{bmatrix} 1 & 2 & \big| & 0 \\ 3 & 6 & \big| & 0 \end{bmatrix} \implies \begin{bmatrix} 1 & 2 & \big| & 0 \\ 0 & 0 & \big| & 0 \end{bmatrix} x1+2x2=0⟹x1=−2x2⟹𝒙=[x1x2]=[−2x2x2]=x2[−21]x_1 + 2x_2 = 0 \implies x_1 = -2x_2 \implies \mathbf{x} = \begin{bmatrix} x_1 \\ x_2 \end{bmatrix} = \begin{bmatrix} -2x_2 \\ x_2 \end{bmatrix} = x_2 \begin{bmatrix} -2 \\ 1 \end{bmatrix} So 𝑵(A)\mathbf{N}(A) spans the line defined by all multiples of the vector [−21]\begin{bmatrix} -2 \\ 1 \end{bmatrix}.

7.4.0.4 Example

Find a spanning set for the null space of the 3×53 \times 5 matrix A=[−36−11−71−223−12−458−4]A = \begin{bmatrix} -3 & 6 & -1 & 1 & -7 \\ 1 & -2 & 2 & 3 & -1 \\ 2 & -4 & 5 & 8 & -4 \end{bmatrix} Performing elimination on augmented matrix [A0][A \quad 0] to reach RREF, we get [1−20−13|00012−2|000000|0]\begin{bmatrix} 1 & -2 & 0 & -1 & 3 & \big| & 0 \\ 0 & 0 & 1 & 2 & -2 & \big| & 0 \\ 0 & 0 & 0 & 0 & 0 & \big| & 0 \end{bmatrix} Expressing the pivots in terms of free variables, x1=2x2+x4−3x5,x3=−2x4+2x5x_1 = 2x_2 + x_4 - 3x_5, \quad x_3 = -2x_4 + 2x_5 Finally expressing the solution vector as a linear combination where the weights are the free variables, [x1x2x3x4x5]=[2x2+x4−3x5x2−2x4+x5x4x5]=x2[21000]+x4[10−210]+x5[−30201]=x2𝒖+x4𝒗+x5𝒘\begin{bmatrix} x_1 \\ x_2 \\ x_3 \\ x_4 \\ x_5 \end{bmatrix} = \begin{bmatrix} 2x_2 + x_4 - 3x_5 \\ x_2 \\ -2x_4 + x_5 \\ x_4 \\ x_5 \end{bmatrix} = x_2 \begin{bmatrix} 2 \\ 1 \\ 0 \\ 0 \\ 0 \end{bmatrix} + x_4 \begin{bmatrix} 1 \\ 0 \\ -2 \\ 1 \\ 0 \end{bmatrix} + x_5 \begin{bmatrix} -3 \\ 0 \\ 2 \\ 0 \\ 1 \end{bmatrix} = x_2\mathbf{u} + x_4\mathbf{v} + x_5\mathbf{w} Every linear combination of 𝒖\mathbf{u}, 𝒗\mathbf{v}, 𝒘\mathbf{w} is an element of 𝑵(A)\mathbf{N}(A). Thus {𝒖,𝒗,𝒘}\{\mathbf{u}, \mathbf{v}, \mathbf{w}\} form the spanning set for 𝑵(A)\mathbf{N}(A).

7.5 The Column Space of a Matrix

The column space of an m×nm \times n matrix AA is the set of all possible linear combinations of the columns of AA. If A=[𝒂1⋯𝒂n]A = [\mathbf{a}_1 \cdots \mathbf{a}_n], then 𝑪(A)=Span{𝒂1,…,𝒂n}\mathbf{C}(A) = Span\{\mathbf{a}_1, \dots, \mathbf{a}_n\} Alternatively, it can be defined as the set of all vectors 𝒃\mathbf{b} that can be produced by multiplying AA with some vector 𝒙\mathbf{x} in ℝn\mathbb{R}^n: 𝑪(A)={𝒃:𝒃=A𝒙 for some 𝒙 in ℝn}\mathbf{C}(A) = \{\mathbf{b} : \mathbf{b} = A\mathbf{x} \text{ for some } \mathbf{x} \text{ in } \mathbb{R}^n\}.

7.5.0.1 Theorem

The column space of an m×nm \times n matrix AA is a valid subspace of ℝm\mathbb{R}^m.
Proof. Each of 𝒂1,…,𝒂n\mathbf{a}_1, \dots, \mathbf{a}_n has mm entries, so they are vectors in ℝm\mathbb{R}^m. By earlier theorem, 𝑪(A)\mathbf{C}(A) is the span of some vectors in ℝm\mathbb{R}^m, so 𝑪(A)\mathbf{C}(A) is a subspace of ℝm\mathbb{R}^m.

Connection to consistency Considering the equation A𝒙=𝒃A\mathbf{x}=\mathbf{b}, if it has a solution, it means that 𝒃\mathbf{b} can be expressed as linear combination of columns of AA, with the entries of 𝒙\mathbf{x} as the weights.

For a m×nm \times n matrix AA,

Comparison: Null Space vs. Column Space
Feature Null Space 𝑵(A)\mathbf{N}(A) Column Space 𝑪(A)\mathbf{C}(A)
Space 𝑵(A)\mathbf{N}(A) is a subspace of ℝn\mathbb{R}^n. 𝑪(A)\mathbf{C}(A) is a subspace of ℝm\mathbb{R}^m.
Definition Style 𝑵(A)\mathbf{N}(A) is implicitly defined; that is, you are given only a condition (A𝒙=𝟎A\mathbf{x}=\mathbf{0}) that vectors in 𝑵(A)\mathbf{N}(A) must satisfy. 𝑪(A)\mathbf{C}(A) is explicitly defined; that is, you are told how to build vectors in 𝑪(A)\mathbf{C}(A).
Finding Vectors It takes time to find vectors in 𝑵(A)\mathbf{N}(A). Row operations on [A𝟎][A \quad \mathbf{0}] are required. It is easy to find vectors in 𝑪(A)\mathbf{C}(A). The columns of AA are displayed; others are formed from them.
Relation to Matrix AA There is no obvious relation between 𝑵(A)\mathbf{N}(A) and the entries in AA. There is an obvious relation between 𝑪(A)\mathbf{C}(A) and the entries in AA, since each column of AA is in 𝑪(A)\mathbf{C}(A).
Checking a Specific Vector vv A typical vector 𝒗\mathbf{v} in 𝑵(A)\mathbf{N}(A) has the property that A𝒗=0A\mathbf{v}=0. Given a specific vector 𝒗\mathbf{v}, it is easy to tell if 𝒗\mathbf{v} is in 𝑵(A)\mathbf{N}(A). Just compute A𝒗A\mathbf{v}. A typical vector 𝒗\mathbf{v} in 𝑪(A)\mathbf{C}(A) has the property that the equation A𝒙=𝒗A\mathbf{x}=\mathbf{v} is consistent. Given a specific vector 𝒗\mathbf{v}, it may take time to tell if 𝒗\mathbf{v} is in 𝑪(A)\mathbf{C}(A). Row operations on [A𝒗][A \quad \mathbf{v}] are required.
Mapping Conditions 𝑵(A)={𝟎}\mathbf{N}(A) = \{\mathbf{0}\} if and only if the equation A𝒙=𝟎A\mathbf{x}=\mathbf{0} has only the trivial solution. This happens if and only if the linear transformation 𝒙↦A𝒙\mathbf{x} \mapsto A\mathbf{x} is one-to-one. 𝑪(A)=ℝm\mathbf{C}(A) = \mathbb{R}^m if and only if the equation A𝒙=𝒃A\mathbf{x}=\mathbf{b} has a solution for every 𝒃\mathbf{b} in ℝm\mathbb{R}^m. This happens if and only if the linear transformation 𝒙↦A𝒙\mathbf{x} \mapsto A\mathbf{x} maps ℝn\mathbb{R}^n onto ℝm\mathbb{R}^m.

7.5.0.2 Example

A=[24−21−2−57337−86]A = \begin{bmatrix} 2 & 4 & -2 & 1 \\ -2 & -5 & 7 & 3 \\ 3 & 7 & -8 & 6 \end{bmatrix}

  1. If the column space of A is a subspace of ℝk\mathbb{R}^k, what is kk?

    • k=3k = 3. Because AA has 3 rows, each column vector contains exactly 3 entries. Therefore, the columns are vectors in ℝ3\mathbb{R}^3, making C(A)C(A) a subspace of ℝ3\mathbb{R}^3.

  2. If the null space of A is a subspace of ℝk\mathbb{R}^k, what is kk?

    • k=4k = 4. The null space N(A)N(A) consists of all vectors 𝒙\mathbf{x} such that A𝒙=𝟎A\mathbf{x} = \mathbf{0}. Because AA has 4 columns, the vector xx must have 4 entries for the matrix multiplication to be defined. Therefore, 𝒙\mathbf{x} is a vector in ℝ4\mathbb{R}^4, making N(A)N(A) a subspace of ℝ4\mathbb{R}^4.

  3. Find a nonzero vector in C(A)C(A) and a nonzero vector in N(A)N(A).

    1. Select any column directly from matrix A. Choosing the first column gives the vector [2−23]\begin{bmatrix} 2 \\ -2 \\ 3 \end{bmatrix}.

    2. We must find a vector 𝒙\mathbf{x} that satisfies A𝒙=𝟎A\mathbf{x} = \mathbf{0}. We set up the augmented matrix [A0][A \quad 0] and perform row reduction: [24−210−2−573037−860]∼[12−10.500−15400008.50]\begin{bmatrix} 2 & 4 & -2 & 1 & 0 \\ -2 & -5 & 7 & 3 & 0 \\ 3 & 7 & -8 & 6 & 0 \end{bmatrix} \sim \begin{bmatrix} 1 & 2 & -1 & 0.5 & 0 \\ 0 & -1 & 5 & 4 & 0 \\ 0 & 0 & 0 & 8.5 & 0 \end{bmatrix} We get x4=0x_4 = 0, x2=5x3x_2 = 5x_3, x1=−9x3x_1 = -9x_3. Hence 𝒙=[−9510]\mathbf{x} = \begin{bmatrix} -9 \\ 5 \\ 1 \\ 0 \end{bmatrix} is a valid nonzero vector in N(A)N(A).

  4. Given 𝒖=[3−2−10]\mathbf{u} = \begin{bmatrix} 3 \\ -2 \\ -1 \\ 0 \end{bmatrix} and 𝒗=[3−13]\mathbf{v} = \begin{bmatrix} 3 \\ -1 \\ 3 \end{bmatrix},

    1. Determine if 𝒖\mathbf{u} is in 𝑵(A)\mathbf{N}(A). Could 𝒖\mathbf{u} be in 𝑪(A)\mathbf{C}(A)?
      A𝒖=[24−21−2−57337−86][3−2−10]=[6−8+2+0−6+10−7+09−14+8+0]=[0−33]A\mathbf{u} = \begin{bmatrix} 2 & 4 & -2 & 1 \\ -2 & -5 & 7 & 3 \\ 3 & 7 & -8 & 6 \end{bmatrix} \begin{bmatrix} 3 \\ -2 \\ -1 \\ 0 \end{bmatrix} = \begin{bmatrix} 6 - 8 + 2 + 0 \\ -6 + 10 - 7 + 0 \\ 9 - 14 + 8 + 0 \end{bmatrix} = \begin{bmatrix} 0 \\ -3 \\ 3 \end{bmatrix} Since A𝒖≠0A\mathbf{u} \neq 0, the vector 𝒖\mathbf{u} is not in 𝑵(A)\mathbf{N}(A).
      The vector 𝒖\mathbf{u} has 4 entries (ℝ4\mathbb{R}^4), but 𝑪(A)\mathbf{C}(A) is strictly a subspace of ℝ3\mathbb{R}^3, so 𝒖\mathbf{u} cannot be in 𝑪(A)\mathbf{C}(A).

    2. Determine if 𝒗\mathbf{v} is in 𝑪(A)\mathbf{C}(A). Could 𝒗\mathbf{v} be in 𝑵(A)\mathbf{N}(A)?
      Set up the augmented matrix [A𝒗][A \quad \mathbf{v}] and row reduce it: [24−213−2−573−137−863]∼…∼[12−10.51.50−15420008.50.5]\begin{bmatrix} 2 & 4 & -2 & 1 & 3 \\ -2 & -5 & 7 & 3 & -1 \\ 3 & 7 & -8 & 6 & 3 \end{bmatrix} \sim \dots \sim \begin{bmatrix} 1 & 2 & -1 & 0.5 & 1.5 \\ 0 & -1 & 5 & 4 & 2 \\ 0 & 0 & 0 & 8.5 & 0.5 \end{bmatrix} The system is consistent, so 𝒗\mathbf{v} is in 𝑪(A)\mathbf{C}(A).
      The vector 𝒗\mathbf{v} has 3 entries (ℝ3\mathbb{R}^3), but 𝑵(A)\mathbf{N}(A) is strictly a subspace of ℝ4\mathbb{R}^4, so 𝒗\mathbf{v} cannot be in 𝑵(A)\mathbf{N}(A).

7.6 Kernel and Range of a Linear Transformation

Defining a Linear Transformation A linear transformation TT from a vector space VV to a vector space WW is a rule that assigns every vector 𝒙\mathbf{x} in VV to a unique vector T(𝒙)T(\mathbf{x}) in WW. It must sastify:

  1. Addition: T(𝒖+𝒗)=T(𝒖)+T(𝒗)T(\mathbf{u}+\mathbf{v}) = T(\mathbf{u}) + T(\mathbf{v}) for all vectors 𝒖\mathbf{u} and 𝒗\mathbf{v} in VV.

  2. Scalar Multiplication: T(c𝒖)=cT(𝒖)T(c\mathbf{u}) = cT(\mathbf{u}) for all vectors 𝒖\mathbf{u} in VV and all scalars cc.

The kernel (or null space) of TT is the set of all 𝒖\mathbf{u} in VV such that T(𝒖)=𝟎T(\mathbf{u}) = \mathbf{0} (zero vector in WW).
The range of TT is the set of vectors in WW of the form T(𝒙)T(\mathbf{x}) for some 𝒙\mathbf{x} in VV.
Defining linear transformation as matrix multiplication of the form T(𝒙)=A𝒙T(\mathbf{x}) = A\mathbf{x}, the kernel is the null space 𝑵(A)\mathbf{N}(A) and the range is the column space 𝑪(A)\mathbf{C}(A).

7.7 Bases

7.7.0.1 Recall Linear Independence

An indexed set of vectors {𝒗1,...,𝒗p}\{\mathbf{v}_{1},...,\mathbf{v}_{p}\} in a vector space VV is considered linearly independent if the vector equation x1𝒗1+x2𝒗2+⋅⋅⋅+xp𝒗p=𝟎x_{1}\mathbf{v}_{1}+x_{2}\mathbf{v}_{2}+\cdot\cdot\cdot+x_{p}\mathbf{v}_{p}=\mathbf{0} has only the trivial solution.
Conversely, the set is linearly dependent if there are specific weights c1,...,cpc_{1},...,c_{p} (where not all of them are zero) that make the equation c1𝒗1+c2𝒗2+⋅⋅⋅+cp𝒗p=𝟎c_{1}\mathbf{v}_{1}+c_{2}\mathbf{v}_{2}+\cdot\cdot\cdot+c_{p}\mathbf{v}_{p}=\mathbf{0} true.
Some quick checks for independence:

7.7.0.2 Theorem

An indexed set {𝒗1,…,𝒗p}\{\mathbf{v}_{1},\ldots,\mathbf{v}_{p}\} consisting of two or more vectors (where 𝒗1≠𝟎\mathbf{v}_{1}\ne\mathbf{0}) is linearly dependent if and only if some vector 𝒗j\mathbf{v}_{j} in the set is a linear combination of the vectors that precede it (𝒗1,…,𝒗j−1\mathbf{v}_{1},\ldots,\mathbf{v}_{j-1}).

7.7.0.3 Basis

A basis is a set of vectors that acts as a foundational coordinate system for a vector space, and is the smallest spanning set for some space.
Let HH be a subspace of a vector space VV. An indexed set of vectors ℬ={𝒃1,…,𝒃p}\mathcal{B}=\{\mathbf{b}_{1},\ldots,\mathbf{b}_{p}\} in VV is defined as a basis for HH if it satisfies two strict conditions:

  1. Linear Independence: The set ℬ\mathcal{B} must be linearly independent. (There are no redundant vectors.)

  2. Spanning: The subspace that is spanned by ℬ\mathcal{B} must exactly coincide with HH, meaning H=Span{𝒃1,...,𝒃p}H=Span\{\mathbf{b}_{1},...,\mathbf{b}_{p}\}.

Every vector space is a subspace of itself, so the basis of VV is simply a linearly independent set that completely spans VV.
Some common bases:

7.7.0.4 Example

Let 𝒗1=[30−6]\mathbf{v}_{1}=\left[\begin{matrix}3\\ 0\\ -6\end{matrix}\right], 𝒗2=[−417]\mathbf{v}_{2}=\left[\begin{matrix}-4\\ 1\\ 7\end{matrix}\right], and 𝒗3=[−215]\mathbf{v}_{3}=\left[\begin{matrix}-2\\ 1\\ 5\end{matrix}\right]. Determine if {𝒗1,𝒗2,𝒗3}\{\mathbf{v}_{1},\mathbf{v}_{2},\mathbf{v}_{3}\} is a basis for ℝ3\mathbb{R}^{3}. A=[3−4−2011−675]⟹[3−4−2011002]A = \begin{bmatrix} 3 & -4 & -2 \\ 0 & 1 & 1 \\ -6 & 7 & 5 \end{bmatrix} \implies \begin{bmatrix} 3 & -4 & -2 \\ 0 & 1 & 1 \\ 0 & 0 & 2 \end{bmatrix} There are 3 pivots, hence the columns are linearly independent, and span ℝ3\mathbb{R}^3, so the set is a basis for ℝ3\mathbb{R}^3.

7.8 Spanning Set Theorem

A set of vectors may span a subspace, but it may contain redundant vectors. How do we trim a spanning set into only the necessary vectors?

7.8.0.1 Theorem

The Spanning Set Theorem Let S={𝒗1,…,𝒗p}S=\{\mathbf{v}_{1},\ldots,\mathbf{v}_{p}\} be a set of vectors in a vector space VV, and let H=Span{𝒗1,...,𝒗p}H=Span\{\mathbf{v}_{1},...,\mathbf{v}_{p}\}.

  1. Removing redundant vectors: If one vector 𝒗k\mathbf{v}_k in SS is a linear combination of the others, then removing vkv_{k} still yields a spanning set for HH.

  2. Extracting the basis: If H≠{𝟎}H\ne\{\mathbf{0}\}, some subset of SS will form a basis for HH.

Note There are infinitely many different possible bases for a vector space. But, every possible basis for some vector space will always have the same number of vectors. For example, In ℝ2\mathbb{R}^2, every valid basis must contain exactly 2 vectors. In ℝ3\mathbb{R}^3, every valid basis must contain exactly 3 vectors. If a subspace HH has a basis of 4 vectors, then every other basis you can possibly find for HH will also have exactly 4 vectors.

7.8.0.2 Example

Let 𝒗1=[02−1]\mathbf{v}_{1}=\begin{bmatrix}0\\ 2\\ -1\end{bmatrix}, 𝒗2=[220]\mathbf{v}_{2}=\begin{bmatrix}2\\ 2\\ 0\end{bmatrix} and 𝒗3=[616−5]\mathbf{v}_{3}=\begin{bmatrix}6\\ 16\\ -5\end{bmatrix}, and let H=Span{𝒗1,𝒗2,𝒗3}H=Span\{\mathbf{v}_{1},\mathbf{v}_{2},\mathbf{v}_{3}\}. Note that 𝒗3=5𝒗1+3𝒗2\mathbf{v}_{3}=5\mathbf{v}_{1}+3\mathbf{v}_{2}.

  1. Span{𝐯1,𝐯2,𝐯3}=Span{𝐯1,𝐯2}Span\{\textbf{v}_{1},\textbf{v}_{2},\textbf{v}_{3}\}=Span\{\textbf{v}_{1},\textbf{v}_{2}\}.
    For any vector 𝒙\mathbf{x} in HH, 𝐱=c1𝐯1+c2𝐯2+c3𝐯3=c1𝐯1+c2𝐯2+c3(5𝐯1+3𝐯2)=(c1+5c3)𝐯1+(c2+3c3)𝐯2\textbf{x} = c_{1}\textbf{v}_{1} + c_{2}\textbf{v}_{2} + c_{3}\textbf{v}_{3} = c_{1}\textbf{v}_{1} + c_{2}\textbf{v}_{2} + c_{3}(5\textbf{v}_{1} + 3\textbf{v}_{2}) = (c_{1} + 5c_{3})\textbf{v}_{1} + (c_{2} + 3c_{3})\textbf{v}_{2}

  2. Find a basis for the subspace HH.
    We know that {𝒗1,𝒗2}\{\mathbf{v}_1,\mathbf{v}_2\} spans HH. Checking for linear indpendence, there is no scalar cc such that c[02−1]=[220]c\begin{bmatrix}0\\ 2\\ -1\end{bmatrix} = \begin{bmatrix}2\\ 2\\ 0\end{bmatrix}. Hence {𝒗1,𝒗2}\{\mathbf{v}_1,\mathbf{v}_2\} forms a basis for HH.

Note that HH is a two-dimensional subspace, in a 3-dimensional parent vector space.

7.9 Bases for 𝑵(A)\mathbf{N}(A) and 𝑪(A)\mathbf{C}(A)

7.9.0.1 Finding a basis for the null space 𝑵(A)\mathbf{N}(A)

We want the exact set of vectors that describe every possible solution to the homogeneous equation A𝐱=𝟎A\textbf{x} = \textbf{0}.

  1. Set up the augmented matrix [A𝟎][A \quad \textbf{0}].

  2. Row reduce to RREF to identify the pivots and free variables.

  3. Express basic variables in terms of free variables.

  4. Express the generic solution vector 𝒙\mathbf{x} in Parametric Vector Form, and the vectors attached to the free variables form the basis for 𝑵(A)\mathbf{N}(A).

7.9.0.2 Example

Imagine after row reducing A𝐱=𝟎A\textbf{x} = \textbf{0}, we get the equations: x1=2x2−5x4x_1 = 2x_2 - 5x_4 and x3=3x4x_3 = 3x_4. Here, x2x_2 and x4x_4 are free variables. 𝐱=[x1x2x3x4]=[2x2−5x4x23x4x4]=x2[2100]+x4[−5031]=x2𝐮+x4𝐯\textbf{x} = \begin{bmatrix} x_1 \\ x_2 \\ x_3 \\ x_4 \end{bmatrix} = \begin{bmatrix} 2x_2 - 5x_4 \\ x_2 \\ 3x_4 \\ x_4 \end{bmatrix} = x_2 \begin{bmatrix} 2 \\ 1 \\ 0 \\ 0 \end{bmatrix} + x_4 \begin{bmatrix} -5 \\ 0 \\ 3 \\ 1 \end{bmatrix} = x_2\textbf{u} + x_4\textbf{v} Note the 1s and 0s in the second and fourth position of 𝒖\mathbf{u} and 𝒗\mathbf{v}, guaranteeing they are independent.
{𝒖,𝒗}\{\mathbf{u}, \mathbf{v}\} is independent and spans the solution space, so they form the basis for 𝑵(A)\mathbf{N}(A).

7.9.0.3 Finding a basis for the column space 𝑪(A)\mathbf{C}(A)

𝑪(A)\mathbf{C}(A) is the span of all the columns of matrix AA, so they are definitely a spanning set. But all the columns may not form a basis, because some columns may be linearly dependent (redundant).

  1. Row reduce the original matrix AA to REF or RREF.

  2. Identify the pivot columns.

  3. Extract the corresponding columns from the ORIGINAL matrix to get the basis for AA. (Row operations alter the column space.)

Row operations alter the columns themselves and hence the column space, so we must take the columns from original matrix as the basis. However, row operations preserve the relationships between the columns. Non-pivot columns are linear combinations of the original pivot columns. By the Spanning Set Theorem, we can ignore the non-pivot columns. The original pivot columns are linearly indepedent and span the space, forming the basis for 𝑪(A)\mathbf{C}(A).

7.9.0.4 Example

Find the basis for the column space for A=[13322695−1−330]A = \begin{bmatrix} 1 & 3 & 3 & 2 \\ 2 & 6 & 9 & 5 \\ -1 & -3 & 3 & 0 \end{bmatrix} Row reducing to echelon form matrix BB, B=[133200310000]B = \begin{bmatrix} 1 & 3 & 3 & 2 \\ 0 & 0 & 3 & 1 \\ 0 & 0 & 0 & 0 \end{bmatrix} The pivots are located in Column 1 and Column 3. Note that Column 2 is just 3 times Column 1, and Column 4 is a combination of Columns 1 and 3. Extracting pivot columns from AA, 𝐚1=[12−1],𝐚3=[393]\textbf{a}_{1} = \begin{bmatrix} 1 \\ 2 \\ -1 \end{bmatrix}, \quad \textbf{a}_{3} = \begin{bmatrix} 3 \\ 9 \\ 3 \end{bmatrix} The basis for C(A)C(A) is the set {𝐚1,𝐚3}\{\textbf{a}_{1}, \textbf{a}_{3}\}. Note that the basis is a 2-dimensional subspace of ℝ3\mathbb{R}^3.

7.10 Coordinate Systems

Let ℬ={𝒃1,…,𝒃n}\mathcal{B} = \{\mathbf{b}_{1}, \dots, \mathbf{b}_{n}\} be a basis for a vector space VV. Because a basis spans the space and is linearly independent, for each vector 𝒙\mathbf{x} in VV, there exists a unique set of scalars c1,…,cnc_{1}, \dots, c_{n} such that: 𝒙=c1𝒃1+…+cn𝒃n\mathbf{x} = c_{1}\mathbf{b}_{1} + \dots + c_{n}\mathbf{b}_{n} These weights c1,…,cnc_{1}, \dots, c_{n} are the coordinates of 𝒙\mathbf{x} relative to the basis ℬ\mathcal{B} (or the ℬ\mathcal{B}-coordinates of 𝒙\mathbf{x}).
When we stack these weights into a column vector, it is called the coordinate vector of 𝒙\mathbf{x} relative to ℬ\mathcal{B}, denoted as [𝒙]ℬ[\mathbf{x}]_{\mathcal{B}}. [𝒙]ℬ=[c1⋮cn][\mathbf{x}]_{\mathcal{B}} = \begin{bmatrix} c_{1} \\ \vdots \\ c_{n} \end{bmatrix} If we place the basis vectors into the columns of a matrix, we get Pℬ=[𝒃1⋯𝒃n]P_{\mathcal{B}} = \begin{bmatrix} \mathbf{b}_{1} & \cdots & \mathbf{b}_{n} \end{bmatrix}. This matrix changes the ℬ\mathcal{B}-coordinates of a vector 𝒙\mathbf{x} into the standard coordinates of 𝒙\mathbf{x}, by this equation: 𝒙=Pℬ[𝒙]ℬ\mathbf{x} = P_{\mathcal{B}}[\mathbf{x}]_{\mathcal{B}} Because the columns of PℬP_{\mathcal{B}} form a basis for ℝn\mathbb{R}^{n}, they are linearly independent, which means PℬP_{\mathcal{B}} is an invertible matrix. Hence we also have Pℬ−1x=[x]ℬP_{\mathcal{B}}^{-1}x = [x]_{\mathcal{B}} The mapping x↦[x]ℬx \mapsto [x]_{\mathcal{B}} is a one-to-one linear transformation from VV onto ℝn\mathbb{R}^n.
Since this mapping preserves vector addition and scalar multiplication, [c1𝒖1+…+cp𝒖p]ℬ=c1[𝒖1]ℬ+…+cp[𝒖p]ℬ[c_{1}\mathbf{u}_{1} + \dots + c_{p}\mathbf{u}_{p}]_{\mathcal{B}} = c_{1}[\mathbf{u}_{1}]_{\mathcal{B}} + \dots + c_{p}[\mathbf{u}_{p}]_{\mathcal{B}} i.e. the ℬ\mathcal{B}-coordinate vector of a linear combination of vectors is the exact same linear combination of their coordinate vectors.

7.10.0.1 Example

Given a basis ℬ={b1,b2}\mathcal{B} = \{b_1, b_2\} for ℝ2\mathbb{R}^2, where b1=[10]b_1 = \begin{bmatrix} 1 \\ 0 \end{bmatrix} and b2=[12]b_2 = \begin{bmatrix} 1 \\ 2 \end{bmatrix}.
A vector 𝒙\mathbf{x} has the coordinate vector [𝒙]ℬ=[−23][\mathbf{x}]_{\mathcal{B}} = \begin{bmatrix} -2 \\ 3 \end{bmatrix}. Find the standard vector 𝒙\mathbf{x}. x=−2[10]+3[12]=[−20]+[36]=[16]x = -2\begin{bmatrix} 1 \\ 0 \end{bmatrix} + 3\begin{bmatrix} 1 \\ 2 \end{bmatrix} = \begin{bmatrix} -2 \\ 0 \end{bmatrix} + \begin{bmatrix} 3 \\ 6 \end{bmatrix} = \begin{bmatrix} 1 \\ 6 \end{bmatrix}

7.10.0.2 Example

Given a basis ℬ={𝒃1,𝒃2}\mathcal{B} = \{\mathbf{b}_1, \mathbf{b}_2\}, 𝒃1=[21]\mathbf{b}_1 = \begin{bmatrix} 2 \\ 1 \end{bmatrix}, 𝒃2=[−11]\mathbf{b}_2 = \begin{bmatrix} -1 \\ 1 \end{bmatrix}.
We have a vector 𝒙\mathbf{x} with standard vector coordinates [45]\begin{bmatrix} 4 \\ 5 \end{bmatrix}.
Find the coordinate vector [𝒙]ℬ[\mathbf{x}]_{\mathcal{B}} of 𝒙\mathbf{x} relative to the basis ℬ\mathcal{B}, i.e. find c1c_1 and c2c_2 such that c1𝒃1+c2𝒃2=𝒙c_1\mathbf{b}_1 + c_2\mathbf{b}_2 = \mathbf{x}. [2−111][c1c2]=[45]\begin{bmatrix} 2 & -1 \\ 1 & 1 \end{bmatrix}\begin{bmatrix} c_1 \\ c_2 \end{bmatrix} = \begin{bmatrix} 4 \\ 5 \end{bmatrix} Solving, we get [x]ℬ=[32][x]_{\mathcal{B}} = \begin{bmatrix} 3 \\ 2 \end{bmatrix}.

7.11 Dimension of a Vector Space

For a vector space VV spanned by a basis, the dimension of VV (or dim⁡V\dim V), is simply the number of vectors in that basis.
dim⁡{𝟎}=0\dim \{\mathbf{0}\} = 0, because there are no linearly independent vectors to form a basis.
dim⁡ℝn=n\dim \mathbb{R}^{n} = n, because the standard basis for ℝn\mathbb{R}^{n} consists of the columns of the n×nn \times n identity matrix (e1,e2,…,ene_{1}, e_{2}, \dots, e_{n}).
dim⁡ℙn=n+1\dim \mathbb{P}_{n} = n + 1, because the standard basis is {1,t,…,tn}\{1, t, \ldots, t^n\}.
In ℝ3\mathbb{R}^3, a 0-D space is a point, a 1-D space is a line, a 2-D space is a plane, and a 3-D space is the entire ℝ3\mathbb{R}^3.

7.11.0.1 Example

Find the dimension of the subspace H={[a−3b+6c5a+4db−2c−d5d]:a,b,c,d in ℝ}H = \left\{ \begin{bmatrix} a - 3b + 6c \\ 5a + 4d \\ b - 2c - d \\ 5d \end{bmatrix} : a, b, c, d \text{ in } \mathbb{R} \right\}. [a−3b+6c5a+4db−2c−d5d]=a[1500]+b[−3010]+c[60−20]+d[04−15]\begin{bmatrix} a - 3b + 6c \\ 5a + 4d \\ b - 2c - d \\ 5d \end{bmatrix} = a \begin{bmatrix} 1 \\ 5 \\ 0 \\ 0 \end{bmatrix} + b \begin{bmatrix} -3 \\ 0 \\ 1 \\ 0 \end{bmatrix} + c \begin{bmatrix} 6 \\ 0 \\ -2 \\ 0 \end{bmatrix} + d \begin{bmatrix} 0 \\ 4 \\ -1 \\ 5 \end{bmatrix} We can see that 𝒗3=−2𝒗2\mathbf{v}_3 = -2\mathbf{v}_2, hence the basis is {𝒗1,𝒗2,𝒗4}\{\mathbf{v}_{1}, \mathbf{v}_{2}, \mathbf{v}_{4}\}. Therefore, dim⁡H=3\dim H = 3.

7.12 Dimensions of 𝑪(A)\mathbf{C}(A) and 𝑵(A)\mathbf{N}(A)

The column space of AA is the span of all its columns. To find a basis for 𝑪(A)\mathbf{C}(A), we row reduce AA to find the pivot columns, and dim⁡𝑪(A)=\dim \mathbf{C}(A) = the number of pivot columns in AA.
The null space of AA is the set of all vectors 𝒙\mathbf{x} that satisfy the homogeneous equation Ax=0Ax = 0. To find a basis for 𝑵(A)\mathbf{N}(A), we row reduce AA, express the basic variables in terms of the free variables, and express the solution in parametric vector form: x=x2𝒗1+𝒙4𝒗2+…x = x_{2}\mathbf{v}_{1} + \mathbf{x}_{4}\mathbf{v}_{2} + \dots Notice that every free variable in the system produces exactly one basis vector (𝒗1,𝒗2\mathbf{v}_{1}, \mathbf{v}_{2}, etc.) in the parametric vector form. Therefore, dim⁡𝑵(A)=\dim \mathbf{N}(A) = to the number of free variables in the equation A𝒙=𝟎A\mathbf{x} = \mathbf{0}.

7.12.0.1 Example

Find the dimension of 𝑪(A)\mathbf{C}(A) and 𝑵(A)\mathbf{N}(A) for A=[−36−11−71−223−12−458−4]A = \begin{bmatrix} -3 & 6 & -1 & 1 & -7 \\ 1 & -2 & 2 & 3 & -1 \\ 2 & -4 & 5 & 8 & -4 \end{bmatrix} The echelon form of the augmented matrix [A𝟎][A \quad \mathbf{0}] is [1−223−100012−20000000]\begin{bmatrix} 1 & -2 & 2 & 3 & -1 & 0 \\ 0 & 0 & 1 & 2 & -2 & 0 \\ 0 & 0 & 0 & 0 & 0 & 0 \end{bmatrix} There are 3 free variables, so dim⁡𝑵(A)=3\dim \mathbf{N}(A) = 3 and there are 2 pivot columns, so dim⁡𝑪(A)=2\dim \mathbf{C}(A) = 2.

7.13 Row Space

The row space of a matrix AA is the set of all linear combinations of the rows (or row vectors) of AA.
The row space is formally 𝑪(AT)\mathbf{C}(A^T), since the rows of a matrix AA are simply the columns of its transpose, ATA^T.
For an m×nm \times n matrix, each row vector contains exactly nn entries. Consequently, the row space C(AT)C(A^T) is a subspace living in ℝn\mathbb{R}^n. (Contrast this with the column space, C(A)C(A), which lives in ℝm\mathbb{R}^m).

7.13.0.1 Theorem

If two matrices AA and BB are row equivalent (can be transformed into the other using EROs), then their row spaces are identical.
Row operations alter the row vectors, but they do not change the span of those rows. Therefore, if we reduce a matrix AA to an echelon form matrix BB, the non-zero rows of BB form a basis for the row space of AA (and also BB).

7.13.0.2 Example

Find bases for the row space, column space, and null space of matrix AA: A=[−2−580−1713−515311−197117−135−3]A = \begin{bmatrix} -2 & -5 & 8 & 0 & -17 \\ 1 & 3 & -5 & 1 & 5 \\ 3 & 11 & -19 & 7 & 1 \\ 1 & 7 & -13 & 5 & -3 \end{bmatrix} First, reduce AA to echelon form matrix BB: A∼B=[13−51501−22−7000−42000000]A \sim B = \begin{bmatrix} 1 & 3 & -5 & 1 & 5 \\ 0 & 1 & -2 & 2 & -7 \\ 0 & 0 & 0 & -4 & 20 \\ 0 & 0 & 0 & 0 & 0 \end{bmatrix}

  1. Basis for 𝑪(AT)={(1,3,−5,1,5),(0,1,−2,2,−7),(0,0,0,−4,20)}\mathbf{C}(A^T) = \{(1, 3, -5, 1, 5),\: (0, 1, -2, 2, -7),\: (0, 0, 0, -4, 20)\}
    Note that we use the rows of BB, not AA.

  2. Basis for 𝑪(A)={[−2131],[−53117],[0175]}\mathbf{C}(A) = \left\{ \begin{bmatrix} -2 \\ 1 \\ 3 \\ 1 \end{bmatrix},\: \begin{bmatrix} -5 \\ 3 \\ 11 \\ 7 \end{bmatrix},\: \begin{bmatrix} 0 \\ 1 \\ 7 \\ 5 \end{bmatrix} \right\}
    Note that we identify the pivot columns using BB, but use the columns of AA as the basis, since EROs do change the column space.

We reduce BB further to RREF matrix CC: B∼C=[1010101−2030001−500000]B \sim C = \begin{bmatrix} 1 & 0 & 1 & 0 & 1 \\ 0 & 1 & -2 & 0 & 3 \\ 0 & 0 & 0 & 1 & -5 \\ 0 & 0 & 0 & 0 & 0 \end{bmatrix} We extract the equations and express the basic variables in terms of the free variables 𝒙3\mathbf{x}_3 and 𝒙5\mathbf{x}_5, then decompose the general solution into parametric vector form.
Hence basis for 𝑵(A)={[−12100],[−1−3051]}\mathbf{N}(A) = \left\{ \begin{bmatrix} -1 \\ 2 \\ 1 \\ 0 \\ 0 \end{bmatrix}, \begin{bmatrix} -1 \\ -3 \\ 0 \\ 5 \\ 1 \end{bmatrix} \right\}

7.14 Rank

The rank of matrix AA is the dimension of 𝑪(A)\mathbf{C}(A).
Recall that the equation Ax=𝒃Ax = \mathbf{b} is only consistent when 𝒃\mathbf{b} is in the column space of AA. Therefore, the rank of AA represents the exact dimension of the set of all vectors 𝒃\mathbf{b} that make the system consistent.
Implication: If you have an n×nn \times n matrix AA with a rank of pp, the transformation 𝒙↦A𝒙\mathbf{x} \mapsto A\mathbf{x} takes vectors from an nn-dimensional space and maps them into a pp-dimensional space.

For an m×nm \times n matrix AA, dim⁡𝑪(A)=dim⁡𝑪(AT)\dim \mathbf{C}(A) = \dim \mathbf{C}(A^T). A matrix has the exact same number of linearly independent rows as it does linearly independent columns. (because they are both defined by the number of pivots.)

7.14.0.1 Example

Considering this matrix AA and row reducing it to REF form: A=[−2−580−1713−515311−197117−135−3]∼B=[13−51501−22−7000−42000000]A = \begin{bmatrix} -2 & -5 & 8 & 0 & -17 \\ 1 & 3 & -5 & 1 & 5 \\ 3 & 11 & -19 & 7 & 1 \\ 1 & 7 & -13 & 5 & -3 \end{bmatrix} \sim B = \begin{bmatrix} 1 & 3 & -5 & 1 & 5 \\ 0 & 1 & -2 & 2 & -7 \\ 0 & 0 & 0 & -4 & 20 \\ 0 & 0 & 0 & 0 & 0 \end{bmatrix} The rank of AA is 3 since there are 3 pivot columns. The basis for the row space of AA is formed by the rows in BB that contain those exact same pivots.

7.14.0.2 Theorem

The Rank Theorem: If matrix AA has nn columns, then rank of A+dim⁡𝑵(A)=n\text{rank of } A + \dim \mathbf{N}(A) = n (Because rank of AA is dim⁡𝑪(A)\dim \mathbf{C}(A), i.e. the number of pivot columns, and dim⁡𝑵(A)\dim \mathbf{N}(A) is the number of free variables, i.e. the number of non-pivot columns.)
This links the solutions of A𝒙=𝟎A\mathbf{x}=\mathbf{0}, i.e. 𝑵(A)\mathbf{N}(A), to the vectors 𝒃\mathbf{b} that make A𝒙=𝒃A\mathbf{x}=\mathbf{b} consistent, i.e. 𝑪(A)\mathbf{C}(A).

7.14.0.3 Example

There is a homogeneous system of 40 equations in 42 variables. Someone has found two linearly independent solutions that span N(A)N(A). Can we be certain that an associated nonhomogeneous system (with the same coefficients) will definitely have a solution?
We have a 40×4240 \times 42 matrix AA. The two solutions form a basis for 𝑵(A)\mathbf{N}(A), so dim⁡𝑵(A)=2\dim \mathbf{N}(A) = 2.
By the Rank Theorem, the rank of AA or dim⁡𝑪(A)=42−2=40\dim \mathbf{C}(A) = 42 - 2 = 40.
The column space is a subspace of ℝ40\mathbb{R}^{40} (vectors with 40 rows). Since we know its dimension is 40, we know that 𝑪(A)\mathbf{C}(A) spans the entirety of ℝ40\mathbb{R}^{40}. Hence, every possible vector 𝒃∈ℝ40\mathbf{b} \in \mathbb{R}^{40} is in 𝑪(A)\mathbf{C}(A) and every non-homogeneous equation A𝒙=𝒃A\mathbf{x} = \mathbf{b} is guaranteed to have a solution.

7.14.0.4 Expanding the Invertible Matrix Theorem

For an n×nn \times n matrix AA, the following statements are logically equivalent to the statement "AA is invertible":

  1. The columns of AA form a basis for ℝn\mathbb{R}^{n}.

  2. 𝑪(A)=ℝn\mathbf{C}(A) = \mathbb{R}^{n}.

  3. dim⁡𝑪(A)=n\dim \mathbf{C}(A) = n.

  4. Rank of A=nA = n.

  5. 𝑵(A)={𝟎}\mathbf{N}(A) = \{\mathbf{0}\}.

  6. dim⁡𝑵(A)=0\dim \mathbf{N}(A) = 0.

7.15 Change of bases

Earlier, we saw how to translate a vector’s coordinates between a custom basis (ℬ\mathcal{B}) and the standard basis (ℰ\mathcal{E}) using the matrix PℬP_{\mathcal{B}}.
We can bypass the standard basis entirely if we want translate coordinates in some basis ℬ\mathcal{B} to another basis 𝒞\mathcal{C}. (i.e. we are looking at the same point in space, but changing how the grid lines are drawn.)
We can build a change of coordinates matrix from ℬ\mathcal{B} to 𝒞\mathcal{C}, i.e. P𝒞←ℬP_{\mathcal{C}\leftarrow\mathcal{B}}.
How to build it: We need to know how to express the ℬ\mathcal{B} basis vectors (𝒃1,…,𝒃n\mathbf{b}_1, \dots, \mathbf{b}_n) in terms of 𝒞\mathcal{C}-coordinates. Then,
P𝒞←ℬ=[[𝒃1]𝒞[𝒃2]𝒞…[𝒃n]𝒞]P_{\mathcal{C}\leftarrow\mathcal{B}} = \begin{bmatrix} [\mathbf{b}_1]_{\mathcal{C}} & [\mathbf{b}_2]_{\mathcal{C}} & \dots & [\mathbf{b}_n]_{\mathcal{C}} \end{bmatrix} With this matrix, we can translate any vector 𝒙\mathbf{x} by matrix multiplication: [𝒙]𝒞=P𝒞←ℬ⋅[𝒙]ℬ[\mathbf{x}]_{\mathcal{C}} = P_{\mathcal{C}\leftarrow\mathcal{B}} \cdot [\mathbf{x}]_{\mathcal{B}}

7.15.0.1 Example

We have a vector space with two different bases, ℬ={b1,b2}\mathcal{B} = \{b_1, b_2\} and 𝒞={c1,c2}\mathcal{C} = \{c_1, c_2\}, and are given a vector 𝒙\mathbf{x} whose coordinates in the ℬ\mathcal{B} system are [x]ℬ=[31][x]_{\mathcal{B}} = \begin{bmatrix} 3 \\ 1 \end{bmatrix}.
Also, [𝒃1]𝒞=[41][\mathbf{b}_1]_{\mathcal{C}} = \begin{bmatrix} 4 \\ 1 \end{bmatrix} and [𝒃2]𝒞=[−61][\mathbf{b}_2]_{\mathcal{C}} = \begin{bmatrix} -6 \\ 1 \end{bmatrix}. Find [𝒙]𝒞[\mathbf{x}]_{\mathcal{C}}.
Building the change of coordinates of matrix, P𝒞←ℬ=[4−611]P_{\mathcal{C}\leftarrow\mathcal{B}} = \begin{bmatrix} 4 & -6 \\ 1 & 1 \end{bmatrix} Performing matrix multiplication, [x]𝒞=[4−611][31]=[64][x]_{\mathcal{C}} = \begin{bmatrix} 4 & -6 \\ 1 & 1 \end{bmatrix} \begin{bmatrix} 3 \\ 1 \end{bmatrix} = \begin{bmatrix} 6 \\ 4 \end{bmatrix} Alternatively, we know that [𝒙]ℬ=[31]⟹𝒙=3𝒃1+1𝒃2[\mathbf{x}]_{\mathcal{B}} = \begin{bmatrix} 3 \\ 1 \end{bmatrix} \implies \mathbf{x} = 3\mathbf{b}_1 + 1\mathbf{b}_2.
Applying the 𝒞\mathcal{C}-coordinate mapping to the whole equation, [𝒙]𝒞=[3𝒃1+𝒃2]𝒞=3[𝒃1]𝒞+[𝒃2]𝒞\begin{align*} [\mathbf{x}]_{\mathcal{C}} &= [3\mathbf{b}_1 + \mathbf{b}_2]_{\mathcal{C}} \\ &= 3[\mathbf{b}_1]_{\mathcal{C}} + [\mathbf{b}_2]_{\mathcal{C}} \\ \end{align*}

Invertibility of PP: Because both ℬ\mathcal{B} and 𝒞\mathcal{C} are valid bases, their vectors are linearly independent. This means that the change of coordinates matrix P𝒞←ℬP_{\mathcal{C}\leftarrow\mathcal{B}} is always invertible. Pℬ←𝒞=(P𝒞←ℬ)−1P_{\mathcal{B}\leftarrow\mathcal{C}} = (P_{\mathcal{C}\leftarrow\mathcal{B}})^{-1}

Alternative to calculating PP: If we are not directly given the 𝒞\mathcal{C}-coordinates of the ℬ\mathcal{B} basis vectors (like [b1]𝒞[b_1]_{\mathcal{C}}), then we need both bases expressed in a common coordinate system (likely the standard ℝn\mathbb{R}^n vectors).
We set up an augmented matrix with the target basis 𝒞\mathcal{C} and the starting basis 𝒞\mathcal{C}: [𝒞ℬ]\left[ \begin{array}{c|c} \mathcal{C} & \mathcal{B} \end{array} \right] Then perform row reduction until we get II on the left: [IP𝒞←ℬ]\left[ \begin{array}{c|c} I & P_{\mathcal{C}\leftarrow\mathcal{B}} \end{array} \right]

7.15.0.2 Example

We have

Set up the augmented matrix: [11212354]∼[101−10112]\left[ \begin{array}{cc|cc} 1 & 1 & 2 & 1 \\ 2 & 3 & 5 & 4 \end{array} \right] \sim \left[ \begin{array}{cc|cc} 1 & 0 & 1 & -1 \\ 0 & 1 & 1 & 2 \end{array} \right] Hence P𝒞←ℬ=[1−112]P_{\mathcal{C}\leftarrow\mathcal{B}} = \begin{bmatrix} 1 & -1 \\ 1 & 2 \end{bmatrix}