Link Search Menu Expand Document

Equations for solvability of inhomogeneous systems

Gaussian Elimination and row echelon form provide a simple and efficient means to check if a system $(A \mid b)$ has a solution or not.

Perform Gaussian Elimination to reduce $A$ into row echelon form and then see if we have a row of the form $(\mathbf{0} \mid c)$ with $c \neq 0$. If so, there is no solution. If not, there will be a solution.

One useful observation: if our system is of the form $A \mathbf{x} = \mathbf{0}$ to start with, then we will always have $0$ entries in the augmented portion of the matrix at any step of Gaussian Elimination. Thus, systems with $\mathbf{b} = \mathbf{0}$ will always have a solution.

In retrospect, this is pretty easy see. We know $\mathbf{x} = \mathbf{0}$ solves $A \mathbf{x} = \mathbf{0}$.

For such systems, we know that there either

  • there is a unique solution or
  • there are infinitely many solutions.

Homogeneous vs inhomogeneous systems

Given what we learned above, it makes sense to give such systems a name.

A system of the form \(A\mathbf{x} = \mathbf{0}\)
is called a homogeneous linear system.

A system \(A\mathbf{x} = \mathbf{b}\) with $\mathbf{b} \neq \mathbf{0}$ is called inhomogeneous.

Restating our observations, a homogeneous system always has at least one solution, though it may be infinitely many.

On the other hand, inhomogeneous systems may or may not have solutions in general.

What can we say about the set \(\mathcal R(A) := \lbrace \mathbf{b} \mid A\mathbf{x} = \mathbf{b} \text { is solvable} \rbrace \ ?\) This set is called the range of $A$. It, itself, is the solution set of a different linear system!

Ranges as null spaces

Let’s see how this works in an example. For general $\mathbf{b}$, let’s use Gaussian Elimination to simplify the system: \(5 x_1 + 3x_2 + x_3 = b_1 \\ x_1 - x_2 + x_3 = b_2 \\ x_2 + x_3 = b_3 \\ 3x_1 - x_3 = b_4\)

The augemented matrix is \(\begin{pmatrix} 5 & 3 & 1 & \bigg| & b_1 \\ 1 & -1 & 1 & \bigg| & b_2 \\ 0 & 1 & 1 & \bigg| & b_3 \\ 3 & 0 & -1 & \bigg| & b_4 \end{pmatrix}\)

If we convert this to row echelon form, we get \(\begin{pmatrix} 1 & -1 & 1 & \bigg| & b_2 \\ 0 & 1 & 1 & \bigg| & b_3 \\ 0 & 0 & -7 & \bigg| & b_4-3b_2-3b_3 \\ 0 & 0 & 0 & \bigg| & b_1 + b_2/7 - 20b_3/7 - 12b_4/7 \end{pmatrix}\) Thus, we have a solution to $A \mathbf{x} = \mathbf{b}$ if and only if \(b_1 + \frac{1}{7}b_2 - \frac{20}{7}b_3 - \frac{12}{7} b_4 = 0\) or in terms of the new notation \(\mathcal R(A) = \left\lbrace \mathbf{b} \mid b_1 + \frac{1}{7}b_2 - \frac{20}{7}b_3 - \frac{12}{7} b_4 = 0 \right\rbrace\)

If we set \(B = \begin{pmatrix} 1 & 1/7 & -20/7 & -12/7 \end{pmatrix}\) then \(\mathcal R(A) = \mathcal Z(B).\) In other words, the set of $\mathbf{b}$ making $A\mathbf{x} = \mathbf{b}$ solvable is itself, in fact, the solution set of another system $B \mathbf{b} = \mathbf{0}$.

We know how to give a parametric presentation of $\mathcal Z(B)$; we just eliminate the variables \(b_1 = -\frac{1}{7}b_2 + \frac{20}{7}b_3 + \frac{12}{7} b_4.\)

So $b_2,b_3,b_4$ are free variables and $b_1$ is determined by from these by this equation.

This is a general phenomenon.

Theorem. Given a matrix $A$, there exists another matrix $B$ so that \(\mathcal R(A) = \mathcal Z(B).\) Moreover, one can find such a $B$ by applying Gaussian Elimination to $(A \mid b)$ and take $B$ to be the coefficient matrix for the system given by all the zero rows.

Proof. We have built up enough technology to show this quickly.

If $(A^\prime \mid b^\prime)$ is the resulting matrix in row echelon form, then

  • each $b_i^\prime$ is a linear combination of the original $b_j$
  • and $b_i^\prime = 0$ must be satisfied for each zero row $\mathbf{R}_i^{A^\prime}$.

Null spaces as ranges

The previous Theorem says that any range can be written as null space of a (possibly) different matrix. The converse is also true. Any null space can be written as a range. Let’s see why.

Recall that using Gaussian Elimination to achieve reduced row echelon form allows us to identify a set of free variables and bound variables for the system $(A \mid 0)$.

For simplicity of notation let’s permute the columns of the matrix so that $x_1, \ldots, x_r$ are the bound variables and $x_{r+1}, \ldots, x_n$ are the free variables.

If $A^\prime$ is the reduced row echelon form of $A$, then let $C$ be the matrix obtained by deleting the zero rows and the columns containing the pivots for each non-zero.

The parametric solution to the system is obtained by as \(\begin{aligned} x_1 & = -C_{1(r+1)} x_{r+1} - \cdots - C_{1n} x_n \\ x_2 & = -C_{2(r+1)} x_{r+1} - \cdots - C_{2n} x_n \\ & \vdots \\ x_r & = -C_{r(r+1)} x_{r+1} - \cdots - C_{rn} x_n \end{aligned}\)

Consider the $n \times (n-r)$-matrix obtained by stacking the $(n-r) \times (n-r)$ identity matrix underneath $-C$. So \(B = \begin{pmatrix} -C \\ I_{n-r} \end{pmatrix}\)

Let’s see that $B$ is the matrix we seek. We want to check that $\mathcal R(B) = \mathcal Z(A)$.

Let’s suggestively write the input variables to $B$ as $y_{n-r},\ldots,y_n$. Then, using the definition of matrix multiplication, writing out $B \mathbf{y}$ gives \(\begin{aligned} x_1 & = -C_{1(r+1)} y_{r+1} - \cdots - C_{1n} y_n \\ x_2 & = -C_{2(r+1)} y_{r+1} - \cdots - C_{2n} y_n \\ & \vdots \\ x_r & = -C_{r(r+1)} y_{r+1} - \cdots - C_{rn} y_n \\ x_{r+1} & = y_{r+1} \\ x_{r+2} & = y_{r+2} \\ & \vdots \\ x_n & = y_n \end{aligned}\) which reduces to the equations above if we eliminate $y_i$ using $x_i = y_i$.

We have established that every null space can be written as range.

Theorem. For any matrix $A$, there exists another matrix $B$ with \(\mathcal Z(A) = \mathcal R(B).\)

Let’s see this in an example. \(A = \begin{pmatrix} 1 & -1 & 0 & 0 \\ 0 & 2 & -1 & 1 \\ 1 & 1 & 0 & 1 \end{pmatrix}\) Converting $A$ to reduced row echelon form gives \(A = \begin{pmatrix} 1 & 0 & 0 & 1/2 \\ 0 & 1 & 0 & 1/2 \\ 0 & 0 & 1 & 0 \end{pmatrix}\) We don’t have rearrange any of the columns of $A$. The only free variable is $x_4$ and \(C = \begin{pmatrix} 1/2 \\ 1/2 \\ 0 \end{pmatrix}\) so \(B = \begin{pmatrix} -1/2 \\ -1/2 \\ 0 \\ 1 \end{pmatrix}\)