Mathematics / Algebra Basic & Higher Algebra 100% Free Open Access
Chapter 7 • Theory & Derivations

Elementary Operations, RREF & Block Matrices

Elementary Row Operations, Echelon Uniqueness, Rank-Nullity Theorem, Rouché-Capelli Consistency & Schur Complements

§7.1 Elementary Row Operations & Elementary Matrices

1. The Three Elementary Row Operations

For a matrix $A \in M_{m \times n}(\mathbb{F})$, the elementary row operations are:

  1. Type I (Row Interchange): $R_i \leftrightarrow R_j$ (swap rows $i$ and $j$).
  2. Type II (Row Scaling): $R_i \to c R_i$ with $c \ne 0$ (multiply row $i$ by a non-zero scalar).
  3. Type III (Row Addition): $R_i \to R_i + c R_j$ with $i \ne j$ (add a scalar multiple of row $j$ to row $i$).
Each operation is strictly reversible by an elementary operation of the identical type!

2. Elementary Matrices

An elementary matrix $E$ is obtained by applying a single elementary row operation to the identity matrix $I_m$.
Fundamental Principle: Performing a row operation on $A$ is algebraically identical to pre-multiplying $A$ by the corresponding elementary matrix: $$\mathbf{A \xrightarrow{\text{Row Op}} B \iff B = E A}$$ Since each elementary matrix is invertible, two matrices $A$ and $B$ are row equivalent ($A \sim B$) if and only if there exist elementary matrices $E_1, \dots, E_k$ such that: $$B = E_k \cdots E_2 E_1 A = P A, \quad \text{where } P \text{ is invertible}$$

§7.2 Row Echelon Form (REF) & Reduced Row Echelon Form (RREF)

1. Row Echelon Form (REF)

A matrix is in Row Echelon Form if:

  • All rows consisting entirely of zeros are at the bottom.
  • The leading entry (first non-zero entry from the left, called the pivot) of each non-zero row is strictly to the right of the leading entry of the row above it.
  • All entries in a column below a leading pivot are zero.

2. Reduced Row Echelon Form (RREF)

A matrix is in Reduced Row Echelon Form (RREF) if it satisfies REF and additionally:

  1. Every leading pivot entry is equal to $1$.
  2. Each leading pivot $1$ is the sole non-zero entry in its column (all entries above and below the pivot are zero).
Theorem (Uniqueness of RREF): Every matrix $A \in M_{m \times n}$ is row equivalent to a uniquely determined reduced row echelon matrix $\text{rref}(A)$.

§7.3 The Rank of a Matrix & The Rank-Nullity Theorem

1. Row Rank, Column Rank & Matrix Rank

  • The row space $\text{Row}(A) \subseteq \mathbb{F}^n$ is the subspace spanned by the row vectors of $A$. Its dimension is the row rank.
  • The column space $\text{Col}(A) \subseteq \mathbb{F}^m$ is the subspace spanned by the column vectors of $A$. Its dimension is the column rank.
Fundamental Rank Theorem: For any matrix $A \in M_{m \times n}$: $$\mathbf{\text{row rank}(A) = \text{column rank}(A) = \text{rank}(A)}$$ The rank equals the number of non-zero rows (or pivot columns) in $\text{rref}(A)$.

2. The Rank-Nullity Theorem

The nullspace (kernel) of $A$ is $\text{Null}(A) = \{X \in \mathbb{F}^n \mid AX = 0\}$. Its dimension is the nullity of $A$, which equals the number of free variables (non-pivot columns) in $\text{rref}(A)$.
Theorem (Rank-Nullity): For any $m \times n$ matrix $A$: $$\mathbf{\text{rank}(A) + \text{nullity}(A) = n \quad (\text{number of columns})}$$

§7.4 Systems of Linear Equations & The Rouché-Capelli Theorem

1. The Augmented Matrix

A system of $m$ linear equations in $n$ variables $AX = B$ is represented by the augmented matrix $[A \mid B] \in M_{m \times (n+1)}$. Applying Gauss-Jordan elimination transforms $[A \mid B]$ into $[R \mid B']$ in RREF without altering the solution set.

2. The Rouché–Capelli Consistency Theorem

Theorem: The linear system $AX = B$ is consistent (possesses at least one solution) if and only if the rank of the coefficient matrix equals the rank of the augmented matrix: $$\mathbf{\text{rank}(A) = \text{rank}([A \mid B])}$$ Classification of Solution Sets:

  • Inconsistent (No Solutions): $\text{rank}(A) < \text{rank}([A \mid B])$. This occurs if and only if $\text{rref}([A \mid B])$ contains a row of the form $[0, 0, \dots, 0 \mid 1]$.
  • Unique Solution: $\text{rank}(A) = \text{rank}([A \mid B]) = n$ (every column has a pivot, nullity = 0).
  • Infinitely Many Solutions: $\text{rank}(A) = \text{rank}([A \mid B]) = r < n$. The general solution depends on $k = n - r$ arbitrary parameters (free variables).

§7.5 Gauss-Jordan Inversion & Block Matrices

1. Gauss-Jordan Inversion Algorithm

To invert an $n \times n$ matrix $A$, form the partitioned augmented matrix $[A \mid I_n]$. Apply elementary row operations to reduce $A$ to $I_n$: $$\mathbf{[A \mid I_n] \xrightarrow{\text{Gauss-Jordan}} [I_n \mid A^{-1}]}$$ If $\text{rref}(A)$ has fewer than $n$ pivots, $A$ is singular and has no inverse.

2. Block Matrices and the Schur Complement

Let $M = \begin{pmatrix} A & B \\ C & D \end{pmatrix}$ be a partitioned block matrix with $A$ invertible.
The Schur complement of $A$ in $M$ is defined by: $$\mathbf{S \equiv D - C A^{-1} B}$$ We factor $M$ via block Gaussian elimination: $$\begin{pmatrix} A & B \\ C & D \end{pmatrix} = \begin{pmatrix} I & 0 \\ C A^{-1} & I \end{pmatrix} \begin{pmatrix} A & 0 \\ 0 & S \end{pmatrix} \begin{pmatrix} I & A^{-1} B \\ 0 & I \end{pmatrix}$$ Consequently: $$\det(M) = \det(A) \det(S) = \det(A) \det(D - C A^{-1} B)$$ If $S$ is also invertible, the explicit block inverse is: $$\mathbf{M^{-1} = \begin{pmatrix} A^{-1} + A^{-1} B S^{-1} C A^{-1} & -A^{-1} B S^{-1} \\ -S^{-1} C A^{-1} & S^{-1} \end{pmatrix}}$$

TIERED UNIVERSITY HONORS PROBLEMS

Step-by-Step Solved Examination Problems

Comprehensive analytical derivations, multi-tier solutions (Foundational, Intermediate Exam, and Honors/Proof Challenge) with complete line-by-line verification.

Foundational Mechanics Example 7.1: Computing RREF and Matrix Rank

Find the Reduced Row Echelon Form (RREF) and determine the rank of the $3 \times 4$ matrix $A = \begin{pmatrix} 1 & 2 & -1 & 3 \\ 2 & 4 & 1 & 9 \\ 3 & 6 & 2 & 14 \end{pmatrix}$.

Step 1: Eliminate Entries Below Pivot 1
$$\begin{pmatrix} 1 & 2 & -1 & 3 \\ 2 & 4 & 1 & 9 \\ 3 & 6 & 2 & 14 \end{pmatrix} \xrightarrow{\substack{R_2 \to R_2 - 2R_1 \\ R_3 \to R_3 - 3R_1}} \begin{pmatrix} 1 & 2 & -1 & 3 \\ 0 & 0 & 3 & 3 \\ 0 & 0 & 5 & 5 \end{pmatrix}$$

Create zeros in column 1 below row 1.

Step 2: Normalize Pivot 2
$$\xrightarrow{R_2 \to \frac{1}{3}R_2} \begin{pmatrix} 1 & 2 & -1 & 3 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 5 & 5 \end{pmatrix}$$

Scale row 2 so pivot in column 3 becomes 1.

Step 3: Eliminate Entries Above and Below Pivot 2
$$\xrightarrow{\substack{R_1 \to R_1 + R_2 \\ R_3 \to R_3 - 5R_2}} \begin{pmatrix} 1 & 2 & 0 & 4 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 0 \end{pmatrix}$$

Eliminate above and below pivot column 3. Row 3 vanishes completely.

Step 4: Conclude Rank and Pivot Positions
$$\text{Pivots are at } (1, 1) \text{ and } (2, 3). \quad \text{Number of non-zero rows} = 2 \implies \text{rank}(A) = 2$$

There are 2 pivot columns (1 and 3) and 2 free columns (2 and 4).

Final Answer & Physical Insight

\mathbf{\text{rref}(A) = \begin{pmatrix} 1 & 2 & 0 & 4 \\ 0 & 0 & 1 & 1 \\ 0 & 0 & 0 & 0 \end{pmatrix}}; \qquad \mathbf{\text{rank}(A) = 2}.

Intermediate University Exam Example 7.2: Gauss-Jordan Matrix Inversion & Parametric System

Use Gauss-Jordan elimination on $[A \mid I_3]$ to compute the inverse of $A = \begin{pmatrix} 1 & 1 & 2 \\ 2 & 1 & 1 \\ 1 & 2 & 1 \end{pmatrix}$.

Step 1: Set Up Augmented Matrix [A | I₃]
$$\left(\begin{array}{ccc|ccc} 1 & 1 & 2 & 1 & 0 & 0 \\ 2 & 1 & 1 & 0 & 1 & 0 \\ 1 & 2 & 1 & 0 & 0 & 1 \end{array}\right) \xrightarrow{\substack{R_2 \to R_2 - 2R_1 \\ R_3 \to R_3 - R_1}} \left(\begin{array}{ccc|ccc} 1 & 1 & 2 & 1 & 0 & 0 \\ 0 & -1 & -3 & -2 & 1 & 0 \\ 0 & 1 & -1 & -1 & 0 & 1 \end{array}\right)$$

Eliminate entries below first pivot.

Step 2: Pivot on Column 2
$$\xrightarrow{R_2 \to -R_2} \left(\begin{array}{ccc|ccc} 1 & 1 & 2 & 1 & 0 & 0 \\ 0 & 1 & 3 & 2 & -1 & 0 \\ 0 & 1 & -1 & -1 & 0 & 1 \end{array}\right) \xrightarrow{\substack{R_1 \to R_1 - R_2 \\ R_3 \to R_3 - R_2}} \left(\begin{array}{ccc|ccc} 1 & 0 & -1 & -1 & 1 & 0 \\ 0 & 1 & 3 & 2 & -1 & 0 \\ 0 & 0 & -4 & -3 & 1 & 1 \end{array}\right)$$

Clear column 2 above and below row 2.

Step 3: Pivot on Column 3 and Clean Columns Above
$$\xrightarrow{R_3 \to -\frac{1}{4}R_3} \left(\begin{array}{ccc|ccc} 1 & 0 & -1 & -1 & 1 & 0 \\ 0 & 1 & 3 & 2 & -1 & 0 \\ 0 & 0 & 1 & 3/4 & -1/4 & -1/4 \end{array}\right) \xrightarrow{\substack{R_1 \to R_1 + R_3 \\ R_2 \to R_2 - 3R_3}} \left(\begin{array}{ccc|ccc} 1 & 0 & 0 & -1/4 & 3/4 & -1/4 \\ 0 & 1 & 0 & -1/4 & -1/4 & 3/4 \\ 0 & 0 & 1 & 3/4 & -1/4 & -1/4 \end{array}\right)$$

Normalize pivot 3 and eliminate above.

Final Answer & Physical Insight

\mathbf{A^{-1} = \frac{1}{4}\begin{pmatrix} -1 & 3 & -1 \\ -1 & -1 & 3 \\ 3 & -1 & -1 \end{pmatrix}}.

Honors / Proof Challenge Example 7.3: Rouché-Capelli Parameter-Dependent System Analysis

Analyze the consistency and determine the complete solution set of the system for all values of $\lambda \in \mathbb{R}$: $\begin{cases} x + y + z = 1 \\ x + 2y + 4z = \lambda \\ x + 4y + 10z = \lambda^2 \end{cases}$.

Step 1: Set Up Augmented Matrix and Perform Row Reduction
$$\left(\begin{array}{ccc|c} 1 & 1 & 1 & 1 \\ 1 & 2 & 4 & \lambda \\ 1 & 4 & 10 & \lambda^2 \end{array}\right) \xrightarrow{\substack{R_2 \to R_2 - R_1 \\ R_3 \to R_3 - R_1}} \left(\begin{array}{ccc|c} 1 & 1 & 1 & 1 \\ 0 & 1 & 3 & \lambda - 1 \\ 0 & 3 & 9 & \lambda^2 - 1 \end{array}\right) \\ \xrightarrow{R_3 \to R_3 - 3R_2} \left(\begin{array}{ccc|c} 1 & 1 & 1 & 1 \\ 0 & 1 & 3 & \lambda - 1 \\ 0 & 0 & 0 & (\lambda^2 - 1) - 3(\lambda - 1) \end{array}\right)$$

Reduce coefficient matrix toward upper triangular form.

Step 2: Analyze the Inconsistency Entry
$$R_3 \text{ entry: } \lambda^2 - 1 - 3\lambda + 3 = \lambda^2 - 3\lambda + 2 = (\lambda - 1)(\lambda - 2) \\ \text{Rank condition: } \text{rank}(A) = 2 \text{ always.} \\ \text{If } (\lambda - 1)(\lambda - 2) \ne 0 \implies \text{rank}([A \mid B]) = 3 > 2 \implies \text{Inconsistent (No solutions)!}$$

If $\lambda \notin \{1, 2\}$, the third equation is $0 = \text{non-zero}$, so no solution exists.

Step 3: Solve for Consistent Values λ = 1 and λ = 2
$$\text{Case 1: } \lambda = 1: \quad \left(\begin{array}{ccc|c} 1 & 1 & 1 & 1 \\ 0 & 1 & 3 & 0 \\ 0 & 0 & 0 & 0 \end{array}\right) \xrightarrow{R_1 \to R_1 - R_2} \left(\begin{array}{ccc|c} 1 & 0 & -2 & 1 \\ 0 & 1 & 3 & 0 \\ 0 & 0 & 0 & 0 \end{array}\right) \\ \text{Let } z = t \implies x = 1 + 2t, \; y = -3t, \; z = t \\ \text{Case 2: } \lambda = 2: \quad \left(\begin{array}{ccc|c} 1 & 1 & 1 & 1 \\ 0 & 1 & 3 & 1 \\ 0 & 0 & 0 & 0 \end{array}\right) \xrightarrow{R_1 \to R_1 - R_2} \left(\begin{array}{ccc|c} 1 & 0 & -2 & 0 \\ 0 & 1 & 3 & 1 \\ 0 & 0 & 0 & 0 \end{array}\right) \\ \text{Let } z = t \implies x = 2t, \; y = 1 - 3t, \; z = t$$

For $\lambda \in \{1, 2\}$, $\text{rank}(A) = \text{rank}([A \mid B]) = 2 < 3$, yielding 1 free variable $t$.

Final Answer & Physical Insight

\begin{cases} \lambda \ne 1, 2: & \text{Inconsistent (No solution)} \\ \lambda = 1: & (x, y, z) = (1 + 2t, -3t, t), \; t \in \mathbb{R} \\ \lambda = 2: & (x, y, z) = (2t, 1 - 3t, t), \; t \in \mathbb{R} \end{cases}