3 System of Linear Equations

3.1 Linear Equations

A linear equation in nn variables is a1x1+a2x2+โ‹ฏ+anxn=ba_1x_1 + a_2x_2 + \cdots + a_nx_n = b The coefficients a1,โ€ฆ,ana_1,\ldots,a_n cannot all be 0, and variables x1,โ€ฆ,xnx_1,\ldots,x_n are raised only to the first power.

A linear system consisting of mm linear equations in nn variables x1,x2,โ‹ฏ,xnx_1, x_2, \cdots, x_n, is represented as a11x1+a12x2+โ‹ฏ+a1nxn=b1a21x1+a22x2+โ‹ฏ+a2nxn=b2โ‹ฎโ‹ฎโ‹ฎam1x1+am2x2+โ‹ฏ+amnxn=bm\begin{alignat*} {4} a_{11}x_1 &+ a_{12}x_2 &+ \cdots &+ a_{1n}x_n &&= b_1 \\ a_{21}x_1 &+ a_{22}x_2 &+ \cdots &+ a_{2n}x_n &&= b_2 \\ \vdots & \quad & & \quad \vdots && \quad \vdots \\ a_{m1}x_1 &+ a_{m2}x_2 &+ \cdots &+ a_{mn}x_n &&= b_m \end{alignat*}

The solution of a linear system in nn unknowns is a sequence of nn numbers s1,s2,โ€ฆ,sns_1, s_2,\ldots,s_n.

Any system of linear equations behave in one of these ways:

  1. The system has no solution

  2. The system has a unique solution

  3. The system has infintely many solutions

If the linear system has at least one solution, it is consistent.

3.2 Homogenous Linear System

If b1=b2=โ€ฆ=bn=0b_1=b_2=\ldots=b_n=0, it is a homogenous linear system. Every homogenous system has at least one solution, known as the trivial solution, x1=x2=โ€ฆ=xn=0x_1 = x_2 = \ldots = x_n = 0 Hence, every homogenous system is consistent, and

  1. has only the trivial solution, or

  2. has infintely many solutions in addition to the trivial solution

3.3 Representation as Matrix

Consider a linear system of two equations: 5x+y=32xโˆ’y=4\begin{align*} 5x + y &= 3 \\ 2x - y &= 4 \end{align*}

Coefficient matrix: [512โˆ’1]\begin{bmatrix} 5 & 1 \\ 2 & -1 \end{bmatrix}

Augmented matrix: [5132โˆ’14]\begin{bmatrix} 5 & 1 & 3 \\ 2 & -1 & 4 \end{bmatrix}

3.4 The two fundamental questions about any linear system

  1. Is the system consistent?

  2. If a solution exists, is it unique?

The following concepts are tools to answer these questions.

3.5 Row Echelon and Reduced Row Echelon Form

Row equivalent matrices can be transformed to each other by one of these elementary row operations:

  1. Swap: Swapping two rows

  2. Scale: Multiply a row by a non-zero constant

  3. Pivot: Add a multiple of one row to another row

3.5.1 Row Echelon Form

Properties:

  1. All rows having only zero entries are at the bottom.

  2. The leftmost non-zero entry (i.e. leading entry) of every non-zero row (i.e. the pivot) is to the right of the leading entry of every row above.

A matrix can have multiple REF forms, depending on the steps taken. Matrix in row echelon form (and not in reduced row echelon form): [1a0a1a2a3002a4a50001a600000]\begin{bmatrix} 1 & a_0 & a_1 & a_2 & a_3 \\ 0 & 0 & 2 & a_4 & a_5 \\ 0 & 0 & 0 & 1 & a_6 \\ 0 & 0 & 0 & 0 & 0 \end{bmatrix}

3.5.2 Reduced Row Echelon Form

Properties:

  1. The matrix is in row echelon form.

  2. The leading entry (or pivot) of every non-zero row is 1.

  3. Every column containing a leading 1 has zeroes in its other entries.

Every matrix has a unique RREF matrix, regardless of the steps taken.

Example matrix in reduced row echelon form: [100a1a2010a3a4001a5a600000]\begin{bmatrix} 1 & 0 & 0 & a_1 & a_2 \\ 0 & 1 & 0 & a_3 & a_4 \\ 0 & 0 & 1 & a_5 & a_6 \\ 0 & 0 & 0 & 0 & 0 \end{bmatrix}

3.5.3 Gaussian and Gauss-Jordan elimination

Consider this linear of system equations: 2x+yโˆ’z=8โˆ’3xโˆ’y+2z=โˆ’11โˆ’2x+y+2z=โˆ’3\begin{alignat*} {4} 2x &+ y &&- z &&= & 8 \\ -3x &- y &&+ 2z &&= &-11 \\ -2x &+ y &&+ 2z &&= &-3 \end{alignat*}

Row operations Augmented matrix
[21โˆ’18โˆ’3โˆ’12โˆ’11โˆ’212โˆ’3]\left[ \begin{array}{ccc|c} 2 & 1 & -1 & 8 \\ -3 & -1 & 2 & -11 \\ -2 & 1 & 2 & -3 \end{array} \right]
L2+32L1โ†’L2L3+L1โ†’L3\begin{aligned} L_2 + \tfrac{3}{2}L_1 &\rightarrow L_2 \\ L_3 + L_1 &\rightarrow L_3 \end{aligned} [21โˆ’180121210215]\left[ \begin{array}{ccc|c} 2 & 1 & -1 & 8 \\ 0 & \tfrac{1}{2} & \tfrac{1}{2} & 1 \\ 0 & 2 & 1 & 5 \end{array} \right]
L3+โˆ’4L2โ†’L3L_3 + -4L_2 \rightarrow L_3 [21โˆ’1801212100โˆ’11]\left[ \begin{array}{ccc|c} 2 & 1 & -1 & 8 \\ 0 & \frac{1}{2} & \frac{1}{2} & 1 \\ 0 & 0 & -1 & 1 \end{array} \right]

We obtain the the row echelon form using Gaussian elimination. Note that top-down approach uses the rows above to introduce zeroes below the pivots.

This is also the forward phase of Gauss-Jordan elimination.

Row operations Augmented matrix
L1โˆ’L3โ†’L1L2+12L3โ†’L2\begin{aligned} L_1 - L_3 &\rightarrow L_1 \\ L_2 + \frac{1}{2}L_3 &\rightarrow L_2 \end{aligned} [210701203200โˆ’11]\left[ \begin{array}{ccc|c} 2 & 1 & 0 & 7 \\ 0 & \frac{1}{2} & 0 & \frac{3}{2} \\ 0 & 0 & -1 & 1 \end{array} \right]
2L2โ†’L2โˆ’L3โ†’L3\begin{aligned} 2L_2 &\rightarrow L_2 \\ -L_3 &\rightarrow L_3 \end{aligned} [21070103001โˆ’1]\left[ \begin{array}{ccc|c} 2 & 1 & 0 & 7 \\ 0 & 1 & 0 & 3 \\ 0 & 0 & 1 & -1 \end{array} \right]
L1โˆ’L2โ†’L112L1โ†’L1\begin{aligned} L_1 - L_2 &\rightarrow L_1 \\ \frac{1}{2}L_1 &\rightarrow L_1 \end{aligned} [10020103001โˆ’1]\left[ \begin{array}{ccc|c} 1 & 0 & 0 & 2 \\ 0 & 1 & 0 & 3 \\ 0 & 0 & 1 & -1 \end{array} \right]

We obtain the RREF using Gauss-Jordan elimination. Note that bottom-up approach uses the rows below to introduce zeroes above pivots, and to transform the pivots to 1.

This is the backward phase of Gauss-Jordan elimination.

Note that of the EROs, only Scale and Pivot were used.

3.6 Answering the fundamental questions

Using row reduction to solve a linear system:

  1. Write the augmented matrix of the linear system.

  2. Use row reduction algorithms to obtain row echelon form of the augmented matrix. Here we can decide if the system in consistent. If it is, go to the next step.

  3. Continue row reduction algorithm to obtain reduced row echelon form.

  4. Write the system of equations corresponding to this matrix.

  5. Rewrite each nonzero equation so that the leading variable is expressed in terms of any free variables.

3.6.1 Example (Not consistent system)

Determine if this linear system is consistent: x2โˆ’4x3=82x1โˆ’3x2+2x3=14x1โˆ’8x2+12x3=1\begin{alignat*} {4} & & x_2 & {}-{} & 4x_3 &= 8 \\ 2x_1 & {}-{} & 3x_2 & {}+{} & 2x_3 &= 1 \\ 4x_1 & {}-{} & 8x_2 & {}+{} & 12x_3 &= 1 \end{alignat*}

Row operations Augmented matrix
Initial Matrix [01โˆ’482โˆ’3214โˆ’8121]\left[ \begin{array}{ccc|c} 0 & 1 & -4 & 8 \\ 2 & -3 & 2 & 1 \\ 4 & -8 & 12 & 1 \end{array} \right]
R1โ†”R2R_1 \leftrightarrow R_2 [2โˆ’32101โˆ’484โˆ’8121]\left[ \begin{array}{ccc|c} 2 & -3 & 2 & 1 \\ 0 & 1 & -4 & 8 \\ 4 & -8 & 12 & 1 \end{array} \right]
R3โˆ’2R1โ†’R3R_3 - 2R_1 \rightarrow R_3 [2โˆ’32101โˆ’480โˆ’28โˆ’1]\left[ \begin{array}{ccc|c} 2 & -3 & 2 & 1 \\ 0 & 1 & -4 & 8 \\ 0 & -2 & 8 & -1 \end{array} \right]
R3+2R2โ†’R3R_3 + 2R_2 \rightarrow R_3 [2โˆ’32101โˆ’4800015]\left[ \begin{array}{ccc|c} 2 & -3 & 2 & 1 \\ 0 & 1 & -4 & 8 \\ 0 & 0 & 0 & 15 \end{array} \right]

The last equation is 0x1+0x2+0x3=150x_1 + 0x_2 + 0x_3 = 15 i.e. 0=150=15. Therefore the system is inconsistent.

3.6.2 Example (Consistent system, infinite solution)

Determine the solution to this linear system: x1+3x2โˆ’2x3+2x5=02x1+6x2โˆ’5x3โˆ’2x4+4x5โˆ’3x6=โˆ’15x3+10x4+15x6=52x1+6x2+8x4+4x5+18x6=6\begin{alignat*} {7} x_1 & {}+{} & 3x_2 & {}-{} & 2x_3 & & & {}+{} & 2x_5 & & &= 0 \\ 2x_1 & {}+{} & 6x_2 & {}-{} & 5x_3 & {}-{} & 2x_4 & {}+{} & 4x_5 & {}-{} & 3x_6 &= -1 \\ & & & & 5x_3 & {}+{} & 10x_4 & & & {}+{} & 15x_6 &= 5 \\ 2x_1 & {}+{} & 6x_2 & & & {}+{} & 8x_4 & {}+{} & 4x_5 & {}+{} & 18x_6 &= 6 \end{alignat*}

Row operations Augmented matrix
Initial System [13โˆ’2020026โˆ’5โˆ’24โˆ’3โˆ’100510015526084186]\left[ \begin{array}{cccccc|c} 1 & 3 & -2 & 0 & 2 & 0 & 0 \\ 2 & 6 & -5 & -2 & 4 & -3 & -1 \\ 0 & 0 & 5 & 10 & 0 & 15 & 5 \\ 2 & 6 & 0 & 8 & 4 & 18 & 6 \end{array} \right]
R2โˆ’2R1โ†’R2R4โˆ’2R1โ†’R4\begin{aligned} R_2 - 2R_1 &\rightarrow R_2 \\ R_4 - 2R_1 &\rightarrow R_4 \end{aligned} [13โˆ’2020000โˆ’1โˆ’20โˆ’3โˆ’100510015500480186]\left[ \begin{array}{cccccc|c} 1 & 3 & -2 & 0 & 2 & 0 & 0 \\ 0 & 0 & -1 & -2 & 0 & -3 & -1 \\ 0 & 0 & 5 & 10 & 0 & 15 & 5 \\ 0 & 0 & 4 & 8 & 0 & 18 & 6 \end{array} \right]
โˆ’R2โ†’R2-R_2 \rightarrow R_2 [13โˆ’20200001203100510015500480186]\left[ \begin{array}{cccccc|c} 1 & 3 & -2 & 0 & 2 & 0 & 0 \\ 0 & 0 & 1 & 2 & 0 & 3 & 1 \\ 0 & 0 & 5 & 10 & 0 & 15 & 5 \\ 0 & 0 & 4 & 8 & 0 & 18 & 6 \end{array} \right]
R3โˆ’5R2โ†’R3R4โˆ’4R2โ†’R4\begin{aligned} R_3 - 5R_2 &\rightarrow R_3 \\ R_4 - 4R_2 &\rightarrow R_4 \end{aligned} [13โˆ’20200001203100000000000062]\left[ \begin{array}{cccccc|c} 1 & 3 & -2 & 0 & 2 & 0 & 0 \\ 0 & 0 & 1 & 2 & 0 & 3 & 1 \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 6 & 2 \end{array} \right]
R3โ†”R416R3โ†’R3\begin{aligned} R_3 &\leftrightarrow R_4 \\ \tfrac{1}{6}R_3 &\rightarrow R_3 \end{aligned} [13โˆ’202000012031000001130000000]\left[ \begin{array}{cccccc|c} 1 & 3 & -2 & 0 & 2 & 0 & 0 \\ 0 & 0 & 1 & 2 & 0 & 3 & 1 \\ 0 & 0 & 0 & 0 & 0 & 1 & \frac{1}{3} \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 \end{array} \right]
Note: Matrix is now in Row Echelon Form (REF)
R2โˆ’3R3โ†’R2R_2 - 3R_3 \rightarrow R_2 [13โˆ’202000012000000001130000000]\left[ \begin{array}{cccccc|c} 1 & 3 & -2 & 0 & 2 & 0 & 0 \\ 0 & 0 & 1 & 2 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 1 & \frac{1}{3} \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 \end{array} \right]
R1+2R2โ†’R1R_1 + 2R_2 \rightarrow R_1 [13042000012000000001130000000]\left[ \begin{array}{cccccc|c} 1 & 3 & 0 & 4 & 2 & 0 & 0 \\ 0 & 0 & 1 & 2 & 0 & 0 & 0 \\ 0 & 0 & 0 & 0 & 0 & 1 & \frac{1}{3} \\ 0 & 0 & 0 & 0 & 0 & 0 & 0 \end{array} \right]
Note: Matrix is now in Reduced Row Echelon Form (RREF)

We first note that the system is consistent.

Also, not every variable column has a pivot (in either REF or RREF). Hence, the system has infinitely many solutions.

From the RREF, we obtain this system of equations: x1+3x2+4x4+2x5=0x3+2x4=0x6=13\begin{align*} x_1 + 3x_2 + 4x_4 + 2x_5 &= 0 \\ x_3 + 2x_4 &= 0 \\ x_6 &= \frac{1}{3} \end{align*}

Solving for the leading variables in terms of the free variables, x1=โˆ’3x2โˆ’4x4โˆ’2x5x3=โˆ’2x4x6=13\begin{align*} x_1 &= -3x_2 - 4x_4 - 2x_5 \\ x_3 &= -2x_4 \\ x_6 &= \frac{1}{3} \end{align*}

The general solution is expressed parametrically, x1=โˆ’3rโˆ’4sโˆ’2tx2=rx3=โˆ’2sx4=sx5=tx6=13\begin{align*} x_1 &= -3r - 4s - 2t \\ x_2 &= r \\ x_3 &= -2s \\ x_4 &= s \\ x_5 &= t \\ x_6 &= \frac{1}{3} \end{align*}

Alternative: Gaussian elimination with back-substitution for larger linear systems (instead of Gauss-Jordan elimination).

From the REF, x1+3x2โˆ’2x3+2x5=0x3+2x4+3x6=0x6=13\begin{align*} x_1 + 3x_2 - 2x_3 + 2x_5 &= 0 \\ x_3 + 2x_4 + 3x_6 &= 0 \\ x_6 &= \frac{1}{3} \end{align*}

Solving for the leading variables, x1=โˆ’3x2+2x3โˆ’2x5x3=1โˆ’2x4โˆ’3x6x6=13\begin{align*} x_1 &= -3x_2 + 2x_3 - 2x_5 \\ x_3 &= 1 - 2x_4 - 3x_6 \\ x_6 &= \frac{1}{3} \end{align*}

Back-substitute the bottom equations into those above it: x1=โˆ’3x2โˆ’4x4โˆ’2x5x3=โˆ’2x4x6=13\begin{align*} x_1 &= -3x_2 - 4x_4 - 2x_5 \\ x_3 &= -2x_4 \\ x_6 &= \frac{1}{3} \end{align*}

Assigning arbitary parameters r,s,tr,s,t to the free variables, we would get the same general solution.

3.6.3 Example (Consistent system, unique solution)

Determine the solution to this linear system: 2x1+x2โˆ’x3=8โˆ’3x1โˆ’x2+2x3=โˆ’11โˆ’2x1+x2+2x3=โˆ’3\begin{align*} 2x_1 + x_2 - x_3 &= 8 \\ -3x_1 - x_2 + 2x_3 &= -11 \\ -2x_1 + x_2 + 2x_3 &= -3 \end{align*} Row echelon form: [21โˆ’1801212100โˆ’11]\left[ \begin{array}{ccc|c} 2 & 1 & -1 & 8 \\ 0 & \frac{1}{2} & \frac{1}{2} & 1 \\ 0 & 0 & -1 & 1 \end{array} \right] First, note that the system is consistent.

Since every variable column contains a pivot, there is a unique solution. (There are no free variables)

Reduced row echelon form: [10020103001โˆ’1]\left[ \begin{array}{ccc|c} 1 & 0 & 0 & 2 \\ 0 & 1 & 0 & 3 \\ 0 & 0 & 1 & -1 \end{array} \right] Note that the RREF of every linear system with a unique solution is the identity matrix.

The solution to the linear system can be read off directly from the RREF: x1=1,x2=3,x3=โˆ’1x_1=1,\:x_2=3,\:x_3=-1

3.7 Linear combination of vectors

Given vectors ๐’—1,๐’—2,โ‹ฏ,๐’—p\mathbf{v}_1, \mathbf{v}_2, \cdots, \mathbf{v}_p in โ„n\mathbb{R}^n and given scalars c1,c2,โ€ฆ,cpc_1, c_2, \dots, c_p, the vector ๐’š\mathbf{y} defined by ๐’š=c1๐’—1+โ€ฆ+cp๐’—p\mathbf{y} = c_1\mathbf{v}_1 + \dots + c_p\mathbf{v}_p is called a linear combination of ๐’—1,๐’—2,โ€ฆ,๐’—p\mathbf{v}_1, \mathbf{v}_2, \dots, \mathbf{v}_p with weights c1,c2,โ€ฆ,cpc_1, c_2, \dots, c_p.

3.7.1 Example

Let ๐’‚1=[1โˆ’2โˆ’5]\mathbf{a}_1 = \begin{bmatrix} 1 \\ -2 \\ -5 \end{bmatrix}, ๐’‚2=[256]\mathbf{a}_2 = \begin{bmatrix} 2 \\ 5 \\ 6 \end{bmatrix}, and ๐’ƒ=[74โˆ’3]\mathbf{b} = \begin{bmatrix} 7 \\ 4 \\ -3 \end{bmatrix}. Determine whether ๐’ƒ\mathbf{b} can be written as a linear combination of ๐’‚1\mathbf{a}_1 and ๐’‚2\mathbf{a}_2.

The vector equation x1๐’‚1+x2๐’‚2=๐’ƒx_1\mathbf{a}_1 + x_2\mathbf{a}_2 = \mathbf{b} is represented by the augmented matrix [๐’‚1๐’‚2โˆฃ๐’ƒ][ \mathbf{a}_1 \ \mathbf{a}_2 \mid \mathbf{b} ]: [127โˆ’254โˆ’56โˆ’3]\left[ \begin{array}{cc|c} 1 & 2 & 7 \\ -2 & 5 & 4 \\ -5 & 6 & -3 \end{array} \right]

After row reduction, we obtain the Reduced Row Echelon Form (RREF): [103012000]\left[ \begin{array}{cc|c} 1 & 0 & 3 \\ 0 & 1 & 2 \\ 0 & 0 & 0 \end{array} \right]

From the RREF, we can read the unique solution for the weights: x1=3,x2=2x_1 = 3, \quad x_2 = 2

Thus, ๐’ƒ\mathbf{b} is a linear combination of ๐’‚1\mathbf{a}_1 and ๐’‚2\mathbf{a}_2: ๐’ƒ=3๐’‚1+2๐’‚2\mathbf{b} = 3\mathbf{a}_1 + 2\mathbf{a}_2

3.8 Span

If ๐’—1,๐’—2,โ€ฆ,๐’—p\mathbf{v}_1,\mathbf{v}_2,\ldots,\mathbf{v}_p are vectors in โ„n\mathbb{R}^n, then Span{๐’—1,๐’—2,โ€ฆ,๐’—p}\{\mathbf{v}_1,\mathbf{v}_2,\ldots,\mathbf{v}_p\} is the set of all vectors that can be written in the form c1๐’—1+c2๐’—2+โ€ฆ+cp๐’—pc_1\mathbf{v}_1 + c_2\mathbf{v}_2 + \ldots + c_p\mathbf{v}_p where c1,c2,โ€ฆ,cpc_1,c_2,\ldots,c_p are scalars.

It is also called the subset of โ„n\mathbb{R}^n spanned by ๐’—1,๐’—2,โ€ฆ,๐’—n\mathbf{v}_1,\mathbf{v}_2,\ldots,\mathbf{v}_n.

3.8.1 Geometric intepretation

Spans as lines, planes, and real space

3.8.2 Example

Let ๐’‚1=[1โˆ’23]\mathbf{a}_1 = \begin{bmatrix} 1 \\ -2 \\ 3 \end{bmatrix}, ๐’‚2=[5โˆ’13โˆ’3]\mathbf{a}_2 = \begin{bmatrix} 5 \\ -13 \\ -3 \end{bmatrix}, and ๐’ƒ=[โˆ’381]\mathbf{b} = \begin{bmatrix} -3 \\ 8 \\ 1 \end{bmatrix}. Determine if ๐’ƒ\mathbf{b} is in the Span{๐’‚1,๐’‚2}\text{Span}\{\mathbf{a}_1, \mathbf{a}_2\}.

The augmented matrix is: [15โˆ’3โˆ’2โˆ’1383โˆ’31]\left[ \begin{array}{cc|c} 1 & 5 & -3 \\ -2 & -13 & 8 \\ 3 & -3 & 1 \end{array} \right]

Performing row reduction to obtain row echelon form, [15โˆ’30โˆ’3200โˆ’2]\left[ \begin{array}{cc|c} 1 & 5 & -3 \\ 0 & -3 & 2 \\ 0 & 0 & -2 \end{array} \right]

The third row corresponds to the equation 0x1+0x2=โˆ’20x_1 + 0x_2 = -2, which is a contradiction.

Therefore, the system is inconsistent, and ๐’ƒ\mathbf{b} is not in the Span{๐’‚1,๐’‚2}\text{Span}\{\mathbf{a}_1, \mathbf{a}_2\}.

3.9 The Matrix Equation A๐’™=๐’ƒA\mathbf{x} = \mathbf{b}

Recall matrix multiplication as linear combinations of columns of AA.

image

If AA is an mร—nm\times n matrix, with columns a1,a2,โ€ฆ,ana_1, a_2,\ldots,a_n, and if ๐’™\mathbf{x} is in โ„n\mathbb{R}^n, then A๐’™A\mathbf{x} is the linear combination of the columns of A using the corresponding entries in ๐’™\mathbf{x} as weights, A๐’™=[๐’‚1๐’‚1โ‹ฏ๐’‚n][x1โ‹ฎxn]=x1๐’‚1+x1๐’‚1+โ€ฆ+xn๐’‚nA\mathbf{x} = \begin{bmatrix} \mathbf{a}_1 & \mathbf{a}_1 & \cdots & \mathbf{a}_n \\ \end{bmatrix} \begin{bmatrix} x_1 \\ \vdots \\ x_n \end{bmatrix} = x_1\mathbf{a}_1 + x_1\mathbf{a}_1 + \ldots + x_n\mathbf{a}_n

Further, if ๐’ƒ\mathbf{b} is in โ„m\mathbb{R}^m, then A๐’™=๐’ƒA\mathbf{x} = \mathbf{b} has the same solution as x1๐’‚1+x1๐’‚1+โ€ฆ+xn๐’‚n=๐’ƒx_1\mathbf{a}_1 + x_1\mathbf{a}_1 + \ldots + x_n\mathbf{a}_n = \mathbf{b} i.e. the same solution to the linear system represented by this augmented matrix [๐’‚1๐’‚2โ‹ฏ๐’‚n๐’ƒ]\begin{bmatrix} \mathbf{a}_1 & \mathbf{a}_2 & \cdots & \mathbf{a}_n & \mathbf{b} \end{bmatrix}

3.9.1 Example

Let A=[134โˆ’42โˆ’6โˆ’3โˆ’2โˆ’7]A = \begin{bmatrix} 1 & 3 & 4 \\ -4 & 2 & -6 \\ -3 & -2 & -7 \end{bmatrix} and ๐’ƒ=[b1b2b3]\mathbf{b} = \begin{bmatrix} b_1 \\ b_2 \\ b_3 \end{bmatrix}. Is the equation A๐’™=๐’ƒA\mathbf{x} = \mathbf{b} consistent for all possible values of b1,b2,b3b_1, b_2, b_3?

The augmented matrix is [134b1โˆ’42โˆ’6b2โˆ’3โˆ’2โˆ’7b3]\left[ \begin{array}{ccc|c} 1 & 3 & 4 & b_1 \\ -4 & 2 & -6 & b_2 \\ -3 & -2 & -7 & b_3 \end{array} \right]

Performing row reduction, the row echelon form is [134b1014104b1+b2000b1โˆ’12b2+b3]\left[ \begin{array}{ccc|c} 1 & 3 & 4 & b_1 \\ 0 & 14 & 10 & 4b_1 + b_2 \\ 0 & 0 & 0 & b_1 - \frac{1}{2}b_2 + b_3 \end{array} \right] For the system to be consistent, we need b1โˆ’12b2+b3=0b_1 - \frac{1}{2}b_2 + b_3 = 0 The columns of A lie on a plane passing through the origin in โ„3\mathbb{R}^3, and the equation of the plane is 2xโˆ’y+2x=02x-y+2x=0.

3.9.2 Spanning โ„m\mathbb{R}^m

If Span{๐’—1,๐’—2,โ€ฆ,๐’—p\mathbf{v}_1,\mathbf{v}_2,\ldots,\mathbf{v}_p}=โ„m\mathbb{R}^m, every vector in โ„m\mathbb{R}^m is a linear combination of ๐’—1,โ€ฆ,๐’—p\mathbf{v}_1,\ldots,\mathbf{v}_p.

These statements are logically equivalent:

  1. A๐’™=๐’ƒA\mathbf{x} = \mathbf{b} has a solution for each ๐’ƒ\mathbf{b} in โ„m\mathbb{R}^m.
    i.e. Any vector in โ„m\mathbb{R}^m can be found through AA.

  2. Each ๐’ƒ\mathbf{b} in ๐‘นm\mathbf{R}^m is a linear combination of the columns of AA.

  3. The columns of AA span โ„m\mathbb{R}^m.

  4. AA has a pivot position in every row.
    If not, the last row will have all zero entries. Then ๐’ƒ\mathbf{b} can be constructed with a 1 in the last row, making the linear system inconsistent.
    **Not to be confused with a pivot in every variable column, which shows uniqueness of solution for the linear system.

Some properties, for ๐’–\mathbf{u} and ๐’—\mathbf{v} in โ„m\mathbb{R}^m, and cc a scalar:

  1. A(๐’–+๐’—)=A๐’–+A๐’—A(\mathbf{u} + \mathbf{v}) = A\mathbf{u} + A\mathbf{v}

  2. A(c๐’–)=c(A๐’–)A(c\mathbf{u}) = c(A\mathbf{u})

3.9.3 Note on Pivots

A pivot in every column โ‡”\iff there are no free variables โ‡”\iff the linear system A๐’™=๐’ƒA\mathbf{x}=\mathbf{b} has at most one solution. โ‡”\iff Columns of AA are linearly independent. The system is either inconsistent or has a unique solution.

[125012000]\left[ \begin{array}{cc|c} \boxed{1} & 2 & 5 \\ 0 & \boxed{1} & 2 \\ 0 & 0 & 0 \end{array} \right]

Unique Solution.

[125012001]\left[ \begin{array}{cc|c} \boxed{1} & 2 & 5 \\ 0 & \boxed{1} & 2 \\ 0 & 0 & 1 \end{array} \right]

No Solution (0=10=1)
Pivots in every column.

A pivot in every row โ‡”\iff the linear system A๐’™=๐’ƒA\mathbf{x}=\mathbf{b} has at least one solution, for every ๐’ƒ\mathbf{b}. This also means the columns of AA span โ„m\mathbb{R}^m, where mm is the number of rows.

mm pivots mean mm dimensions (mm linearly independent rows) to span โ„m\mathbb{R}^m.

[102b1013b2]\left[ \begin{array}{ccc|c} \boxed{1} & 0 & 2 & b_1 \\ 0 & \boxed{1} & 3 & b_2 \end{array} \right]

Wide matrices with pivot in every row are always consistent and have infinitely many solutions (at least one free variable left over).

[10b101b200b3]\left[ \begin{array}{cc|c} \boxed{1} & 0 & b_1 \\ 0 & \boxed{1} & b_2 \\ 0 & 0 & b_3 \end{array} \right]

Tall matrices cannot have a pivot in every row. Row of zeros leads to "No Solution" if b3โ‰ 0b_3 \neq 0.
Pivots on every row.

There is a pivot in every row and column โ‡”\iff AA is a square matrix โ‡”\iff The system A๐’™=๐’ƒA\mathbf{x}=\mathbf{b} has a unique solution for any ๐’ƒ\mathbf{b}.

[1005010โˆ’20014]\left[ \begin{array}{ccc|c} \boxed{1} & 0 & 0 & 5 \\ 0 & \boxed{1} & 0 & -2 \\ 0 & 0 & \boxed{1} & 4 \end{array} \right]

Pivots in every row and column lead to a unique solution.

Hence, wide matrix can never have unique solution, tall matrix can never have solution for every b.

3.10 Solution Sets of Linear Systems

Parametric vector form of the solution set of A๐’™=๐’ƒA\mathbf{x}=\mathbf{b} is ๐’™=๐’‘+t๐’—,tโˆˆโ„\mathbf{x} = \mathbf{p} + t\mathbf{v},\quad t\in\mathbb{R}

If ๐’ƒ=๐ŸŽ\mathbf{b}=\mathbf{0}, it is a homogenous system of equations.

Parametric vector form for the solution to a homogenous system is ๐’™=t๐’—\mathbf{x} = t\mathbf{v}

๐’‘\mathbf{p} is then one particular solution to A๐’™=๐’ƒA\mathbf{x}=\mathbf{b} (corresponding to t=0).

๐’‘\mathbf{p} basically translates t๐’—t\mathbf{v} to pass through ๐’‘\mathbf{p} instead of the origin.

3.10.1 Example (homogenous system)

Determine if this homogenous system has a nontrivial solution: 3x1+5x2โˆ’4x3=0โˆ’3x1โˆ’2x2+4x3=06x1+x2โˆ’8x3=0\begin{aligned} 3x_1 + 5x_2 - 4x_3 &= 0 \\ -3x_1 - 2x_2 + 4x_3 &= 0 \\ 6x_1 + x_2 - 8x_3 &= 0 \end{aligned} The augmented matrix is [35โˆ’40โˆ’3โˆ’24061โˆ’80]\left[ \begin{array}{ccc|c} 3 & 5 & -4 & 0 \\ -3 & -2 & 4 & 0 \\ 6 & 1 & -8 & 0 \end{array} \right] Performing row reduction operations, the row echelon form is [35โˆ’4003000000]\left[ \begin{array}{ccc|c} 3 & 5 & -4 & 0 \\ 0 & 3 & 0 & 0 \\ 0 & 0 & 0 & 0 \end{array} \right] We get 3x1+5x2โˆ’4x3=03x2=0\begin{aligned} 3x_1 + 5x_2 - 4x_3 &= 0 \\ 3x_2 &= 0 \end{aligned} Solving, 3x1=4x3โŸนx1=43x33x_1 = 4x_3 \implies x_1 = \frac{4}{3}x_3 Since x3x_3 is a free variable, we can express the solution vector ๐’™\mathbf{x} as: ๐’™=[x1x2x3]=[43x30x3]=x3[4301]\mathbf{x} = \begin{bmatrix} x_1 \\ x_2 \\ x_3 \end{bmatrix} = \begin{bmatrix} \frac{4}{3}x_3 \\ 0 \\ x_3 \end{bmatrix} = x_3 \begin{bmatrix} \frac{4}{3} \\ 0 \\ 1 \end{bmatrix} Every solution to A๐’™=๐ŸŽA\mathbf{x}=\mathbf{0} in this case is a scalar multiple of ๐’—=[4301]\mathbf{v}=\begin{bmatrix} \frac{4}{3} \\ 0 \\ 1 \end{bmatrix}

3.10.2 Example (non-homogenous system)

Describe all solutions of this non-homogenous system: 3x1+5x2โˆ’4x3=7โˆ’3x1โˆ’2x2+4x3=โˆ’16x1+x2โˆ’8x3=4\begin{aligned} 3x_1 + 5x_2 - 4x_3 &= 7 \\ -3x_1 - 2x_2 + 4x_3 &= -1 \\ 6x_1 + x_2 - 8x_3 &= 4 \end{aligned} The augmented matrix is [35โˆ’47โˆ’3โˆ’24โˆ’161โˆ’8โˆ’4]\left[ \begin{array}{ccc|c} 3 & 5 & -4 & 7 \\ -3 & -2 & 4 & -1 \\ 6 & 1 & -8 & -4 \end{array} \right] The reduced row echelon form is [10โˆ’43โˆ’101020000]\left[ \begin{array}{ccc|c} 1 & 0 & -\frac{4}{3} & -1 \\ 0 & 1 & 0 & 2 \\ 0 & 0 & 0 & 0 \end{array} \right] We get {x1โˆ’43x3=โˆ’1x2=2x3 is free\left\{ \begin{aligned} x_1 - \frac{4}{3}x_3 &= -1 \\ x_2 &= 2 \\ x_3 &\text{ is free} \end{aligned} \right. The solution set is ๐’™=[x1x2x3]=[โˆ’120]+t[4301],tโˆˆโ„\mathbf{x} = \begin{bmatrix} x_1 \\ x_2 \\ x_3 \end{bmatrix} = \begin{bmatrix} -1 \\ 2 \\ 0 \end{bmatrix} + t \begin{bmatrix} \frac{4}{3} \\ 0 \\ 1 \end{bmatrix},\quad t\in\mathbb{R}

3.11 Linear Independence

A indexed set of vectors {๐’—1,โ€ฆ,๐’—p}\{\mathbf{v}_1,\ldots,\mathbf{v}_p\} in โ„n\mathbb{R}^n is linearly independent iff the vector equation x1๐’—1+x2๐’—2+โ‹ฏ+xp๐’—p=๐ŸŽx_1\mathbf{v}_1 + x_2\mathbf{v}_2 + \cdots + x_p\mathbf{v}_p = \mathbf{0} has only the trivial solution.

If โˆƒ weights c1,โ€ฆ,cp\exists \text{ weights } c_1,\ldots,c_p that are not all zero, s.t. c1๐’—1+c2๐’—2+โ‹ฏ+cp๐’—p=๐ŸŽc_1\mathbf{v}_1 + c_2\mathbf{v}_2 +\cdots+ c_p\mathbf{v}_p = \mathbf{0} then the set of vectors is linearly dependent.

In other words, if S={๐’—1,โ‹ฏ,๐’—p}S = \{\mathbf{v}_1,\cdots,\mathbf{v}_p\} is linearly dependent, and ๐’—1โ‰ 0\mathbf{v}_1 \neq 0, then some ๐’—j\mathbf{v}_j, j>1j>1, is a linear combination of the preceding vectors ๐’—1,โ‹ฏ,๐’—jโˆ’1\mathbf{v}_1,\cdots,\mathbf{v}_{j-1}.

In a matrix, if all the columns have a pivot, then the columns are linearly independent.

3.11.1 Example

Consider the set of vectors in โ„3\mathbb{R}^3: ๐’—1=[110],๐’—2=[101],๐’—3=[312],๐’—4=[011]\mathbf{v}_1 = \begin{bmatrix} 1 \\ 1 \\ 0 \end{bmatrix}, \quad \mathbf{v}_2 = \begin{bmatrix} 1 \\ 0 \\ 1 \end{bmatrix}, \quad \mathbf{v}_3 = \begin{bmatrix} 3 \\ 1 \\ 2 \end{bmatrix}, \quad \mathbf{v}_4 = \begin{bmatrix} 0 \\ 1 \\ 1 \end{bmatrix}

First construct the matrix A=[๐’—1๐’—2๐’—3๐’—4]A = [\mathbf{v}_1 \ \mathbf{v}_2 \ \mathbf{v}_3 \ \mathbf{v}_4] and reduce it to REF to locate the pivots: Aโˆผ[๐Ÿ1300๐Ÿ2โˆ’1000๐Ÿ]A \sim \begin{bmatrix} \mathbf{1} & 1 & 3 & 0 \\ 0 & \mathbf{1} & 2 & -1 \\ 0 & 0 & 0 & \mathbf{1} \end{bmatrix} Column 1 and Column 2 have pivots, Column 3 is the first column without a pivot. This indicates that ๐’—3\mathbf{v}_3 is the first vector in the sequence that is linearly dependent on those preceding it.

We can extract the vectors up to ๐’—3\mathbf{v}_3 into an augmented matrix [๐’—1๐’—2โˆฃ๐’—3][\mathbf{v}_1 \ \mathbf{v}_2 \mid \mathbf{v}_3]: [๐Ÿ130๐Ÿ2000]\left[ \begin{array}{cc|c} \mathbf{1} & 1 & 3 \\ 0 & \mathbf{1} & 2 \\ 0 & 0 & 0 \end{array} \right]

Since ๐’—3\mathbf{v}_3 corresponds to a non-pivot column in the original set, it is redundant. Specifically: ๐’—3=1๐’—1+2๐’—2\mathbf{v}_3 = 1\mathbf{v}_1 + 2\mathbf{v}_2

Takeaway: If there are columns with no pivot, the set of vectors is linearly dependent.

3.11.2 Wide matrices

Naturally, columns of a wide matrix are linearly dependent. np[***************]n \begin{array}{c} \scriptstyle p \\ \left[ \begin{array}{ccccc} * & * & * & * & * \\ * & * & * & * & * \\ * & * & * & * & * \end{array} \right] \end{array} If p>np>n, the matrix cannot have a pivot in every column โ‡’\Rightarrow there must be a free variable โ‡’\Rightarrow A๐’™=๐ŸŽA\mathbf{x} = \mathbf{0} has non trivial solution โ‡’\Rightarrow columns of AA are linerly dependent.

3.11.3 Set with zero vector

If S={๐’—1,โ€ฆ,๐’—p}S = \{\mathbf{v}_1,\ldots,\mathbf{v}_p\} in โ„n\mathbb{R}^n contains the zero vector, then SS is linearly dependent.

WLOG, suppose ๐’—1=๐ŸŽ\mathbf{v}_1 = \mathbf{0}, then 1๐’—1+0๐’—2+โ‹ฏ+0๐’—p=๐ŸŽ1\mathbf{v}_1 + 0\mathbf{v}_2 +\cdots +0\mathbf{v}_p = \mathbf{0} shows that SS is linearly dependent.

3.12 Linear Transformations

A transformation T:โ„nโ†’โ„mT : \mathbb{R}^n \to \mathbb{R}^m is a rule that assigns to each vector ๐’™\mathbf{x} in โ„n\mathbb{R}^n a vector T(๐’™)T(\mathbf{x}) in โ„m\mathbb{R}^m.

A transformation TT is linear if it satisfies:

  1. T(๐’–+๐’—)=T(๐’–)+T(๐’—)T(\mathbf{u} + \mathbf{v}) = T(\mathbf{u}) + T(\mathbf{v})

  2. T(c๐’–)=cT(๐’–)T(c\mathbf{u}) = cT(\mathbf{u})

Generalization: T(c1๐’—1+โ€ฆ+cp๐’—p)=c1T(๐’—1)+โ€ฆ+cpT(๐’—p)T(c_1 \mathbf{v}_1 + \dots + c_p \mathbf{v}_p) = c_1 T(\mathbf{v}_1) + \dots + c_p T(\mathbf{v}_p)

Example: T(๐’™)=r๐’™T(\mathbf{x}) = r\mathbf{x}, T:โ„2โ†’โ„2T : \mathbb{R}^2 \to \mathbb{R}^2 is called a contraction when 0โ‰คrโ‰ค10 \leq r \leq 1 and a dilation when r>1r > 1.

3.12.1 Finding the Standard Matrix

Let T:โ„nโ†’โ„mT:\mathbb{R}^n\rightarrow\mathbb{R}^m be a linear transformation. Then there exists a unique matrix AA such that T(๐’™)=A๐’™ for all ๐’™ in โ„nT(\mathbf{x}) = A\mathbf{x} \text{ for all } \mathbf{x} \text{ in } \mathbb{R}^n AA is also known as the standard matrix for the linear transformation TT.

Here, AA is the mร—nm\times n matrix whose jthj^{\text{th}} column is the vector T(ej)T(e_j), where eje_j is the jthj^{\text{th}} column of the identity matrix in โ„n\mathbb{R}^n, i.e. A=[T(๐’†1)T(๐’†2)โ€ฆT(๐’†n)]A = \begin{bmatrix} T(\mathbf{e}_1) & T(\mathbf{e}_2) & \dots & T(\mathbf{e}_n) \end{bmatrix}

For T:โ„nโ†’โ„mT:\mathbb{R}^n\rightarrow\mathbb{R}^m, T(๐’™)=A๐’™T(\mathbf{x})=A\mathbf{x}, A is an mร—n matrix.A \text{ is an } m\times n \text{ matrix.}

Example:

Let I2=[1001]I_2 = \begin{bmatrix} 1 & 0 \\ 0 & 1 \end{bmatrix}, so ๐’†1=[10]\mathbf{e}_1 = \begin{bmatrix} 1 \\ 0 \end{bmatrix} and ๐’†2=[01]\mathbf{e}_2 = \begin{bmatrix} 0 \\ 1 \end{bmatrix}. Suppose TT is a linear transformation from โ„2\mathbb{R}^2 to โ„3\mathbb{R}^3 such that T(๐’†1)=[5โˆ’72]T(\mathbf{e}_1) = \begin{bmatrix} 5 \\ -7 \\ 2 \end{bmatrix} and T(๐’†2)=[โˆ’380]T(\mathbf{e}_2) = \begin{bmatrix} -3 \\ 8 \\ 0 \end{bmatrix}. Find the formula for the image of an arbitrary ๐’™\mathbf{x} in โ„2\mathbb{R}^2.

There exists a unique matrix AA such that T(๐’™)=A๐’™T(\mathbf{x}) = A\mathbf{x} for all ๐’™\mathbf{x} in โ„n\mathbb{R}^n. The columns of AA are the images of the basis vectors ๐’†1\mathbf{e}_1 and ๐’†2\mathbf{e}_2: A=[T(๐’†1)T(๐’†2)]=[5โˆ’3โˆ’7820]A = [T(\mathbf{e}_1) \quad T(\mathbf{e}_2)] = \begin{bmatrix} 5 & -3 \\ -7 & 8 \\ 2 & 0 \end{bmatrix} For an arbitrary ๐’™=[x1x2]\mathbf{x} = \begin{bmatrix} x_1 \\ x_2 \end{bmatrix} in โ„2\mathbb{R}^2, the formula for the image is: T(๐’™)=[5โˆ’3โˆ’7820][x1x2]=[5x1โˆ’3x2โˆ’7x1+8x22x1]T(\mathbf{x}) = \begin{bmatrix} 5 & -3 \\ -7 & 8 \\ 2 & 0 \end{bmatrix} \begin{bmatrix} x_1 \\ x_2 \end{bmatrix} = \begin{bmatrix} 5x_1 - 3x_2 \\ -7x_1 + 8x_2 \\ 2x_1 \end{bmatrix}

Example:

Let T:โ„2โ†’โ„2T : \mathbb{R}^2 \to \mathbb{R}^2 be the transformation that rotates each point in โ„2\mathbb{R}^2 about the origin through an angle ฯ•\phi, with counterclockwise rotation for a positive angle. Find the standard matrix AA for this transformation.

To find the standard matrix AA, we determine the images of the standard basis vectors ๐’†1=[10]\mathbf{e}_1 = \begin{bmatrix} 1 \\ 0 \end{bmatrix} and ๐’†2=[01]\mathbf{e}_2 = \begin{bmatrix} 0 \\ 1 \end{bmatrix} under the rotation ฯ•\phi:

The standard matrix AA is formed by these columns: A=[T(๐’†1)T(๐’†2)]=[cos⁡ฯ•โˆ’sin⁡ฯ•sin⁡ฯ•cos⁡ฯ•]A = [T(\mathbf{e}_1) \quad T(\mathbf{e}_2)] = \begin{bmatrix} \cos \phi & -\sin \phi \\ \sin \phi & \cos \phi \end{bmatrix}