查看“︁线性代数（英文维基教科书）/Gauss' Method 高斯消元法”︁的源代码

<noinclude>{{navigation|Book=线性代数（英文维基教科书）
|current=Gauss' Method 高斯消元法
|previous=Solving Linear Systems 求解线性方程组
|next=Describing the Solution Set 解集的表示}}</noinclude>

{{TextBox|1=
;Definition 1.1{{anchor|linear equation}}:
A '''linear equation''' in variables
<math>x_1,x_2,\ldots,x_n</math> has the form

未知数<math>x_1,x_2,\ldots,x_n</math>的一个线性方程组，形如

:<math>a_1x_1+a_2x_2+a_3x_3+\cdots+a_nx_n=d</math>

where the numbers <math>a_1,\dots,a_n\in\R</math> are the equation's
'''coefficients''' and
<math>d\in\R</math> is the '''constant'''.
An <math>n</math>-tuple <math>(s_1,s_2,\dots,s_n)\in\R^n </math> is a 
'''solution''' 
of, or '''satisfies''', that equation if substituting the numbers
<math>s_1,\dots,s_n</math> for the variables gives a true statement:
<math>a_1s_1+a_2s_2+\ldots+a_ns_n=d</math>.

其中<math>a_1,\dots,a_n\in\R</math>是方程的系数，<math>d\in\R</math>是一个常数项。
当一个<math>n</math>元数组<math>(s_1,s_2,\dots,s_n)\in\R^n </math>满足：<math>a_1s_1+a_2s_2+\ldots+a_ns_n=d</math>成立
那么我们称这个数组为为方程的一个解。

A '''system of linear equations'''

一个'''线性方程组'''

:<math>
\begin{cases} 
\begin{array}{*{4}{rc}r}
a_{1,1}x_1 &+ &a_{1,2}x_2 &+ &\cdots &+ &a_{1,n}x_n &=  &d_1  \\
a_{2,1}x_1 &+ &a_{2,2}x_2 &+ &\cdots &+ &a_{2,n}x_n &=  &d_2  \\
&  &            &   &       &  &           &\vdots   \\
a_{m,1}x_1 &+ &a_{m,2}x_2 &+ &\cdots &+ &a_{m,n}x_n &=  &d_m
\end{array}
\end{cases}
</math>

has the solution
<math>(s_1,s_2,\ldots,s_n)</math> if that <math>n</math>-tuple is a solution of all of the equations in the system.

当有一个<math>n</math>元数组<math>(s_1,s_2,\dots,s_n)\in\R^n </math>，是方程组中每一个方程的解，那么我们称
这个<math>n</math>元数组为方程组的解。

}}

{{TextBox|1=
;Example 1.2:
The ordered pair <math>(-1,5)</math> is a solution of this system.

有序数对<math>(-1,5)</math>是这个方程组的解。
:<math>
\begin{cases} 
\begin{array}{*{2}{rc}r}
4x_1 &+ &2x_2 &= &6  \\
-x_1 &+ &x_2  &= &6
\end{array}
\end{cases} 
</math>

In contrast, <math> (5,-1) </math> is not a solution.

相反，<math> (5,-1) </math>就不是这个方程组的解。
}}

Finding the set of all solutions is '''solving'''{{anchor|system of linear equations!solving}} the system. No guesswork or good fortune is needed to solve a linear system. There is an algorithm that always works. The next example introduces that algorithm, called '''Gauss' method'''{{anchor|Gauss' method}}. It transforms the system, step by step, into one with a form that is easily solved.

所有满足方程组的解的集合就是这个方程组的解。解一个线性方程不需要猜测或好运。有一种方法法总是有效的。在下一个例子，我们将介绍该算法，它被称为“高斯消元法”。它将方程组一步一步地转换成一个易于求解的形式。

{{TextBox|1=
;Example 1.3:
To solve this system

解方程组

:<math>
\begin{cases} 
\begin{array}{*{3}{rc}r}
&   &      &   &3x_3  &=  &9  \\
x_1 &+  &5x_2  &-  &2x_3  &=  &2  \\
\frac{1}{3}x_1 &+  &2x_2  &   &      &=  &3  
\end{array}
\end{cases} 
</math>

we repeatedly transform it until it is in a form that is easy to solve.

我们反复变换它，直到它变成一种易于求解的形式。

:<math>\begin{array}{rcl}
\quad
&\xrightarrow[]{ \text{swap row 1 with row 3   将第1行与第3行交换} }
& \begin{array}{*{3}{rc}r}
\frac{1}{3}x_1 &+  &2x_2  &   &      &=  &3  \\
x_1 &+  &5x_2  &-  &2x_3  &=  &2  \\
&   &      &   &3x_3  &=  &9  
\end{array}                                    \\
&\xrightarrow[]{ \text{multiply row 1 by 3   第1行乘以3} }
&\begin{array}{*{3}{rc}r}
x_1 &+  &6x_2  &   &      &=  &9  \\
x_1 &+  &5x_2  &-  &2x_3  &=  &2  \\
&   &      &   &3x_3  &=  &9  
\end{array}                                          \\
&\xrightarrow[]{ \text{add }-1\text{ times row 1 to row 2   第一行乘以-1加到第二行} }
&\begin{array}{*{3}{rc}r}
x_1 &+  &6x_2  &   &      &=  &9  \\
&   &-x_2  &-  &2x_3  &=  &-7 \\
&   &      &   &3x_3  &=  &9  
\end{array}
\end{array}
</math>
The third step is the only nontrivial one. We've mentally multiplied both sides of the first row by <math>-1</math>, mentally added that to the old second row, and written the result in as the new second row.

第三步是特殊的一步。我们把第一行的两边都乘以<math>-1</math>，把它加到原来的第二行，然后把结果写成新的第二行。

Now we can find the value of each variable.
The bottom equation shows that <math>x_3=3</math>. Substituting <math>3</math> for <math>x_3</math> in the middle equation shows that <math>x_2=1</math>. Substituting those two into the top equation gives that <math>x_1=3</math> and so the system has a unique solution: the solution set is <math>\{\,(3,1,3)\,\}</math>.

现在我们可以找到每个变量的值。
下面的等式显示<math>x_3=3</math>。将中间方程中的<math>x_3</math>替换为<math>x_2=1</math>。把这两个代入最上面的方程，得到<math>x_1=3</math>，因此方程组有一个唯一的解：解集是<math>\{\,(3,1,3)\,\}</math>。
}}

Most of this subsection and the next one consists of examples of solving linear systems by Gauss' method. We will use it throughout this book. It is fast and easy. But, before we get to those examples, we will first show that this method is also safe in that it never loses solutions or picks up extraneous solutions.

{{TextBox|1=
;Theorem 1.4 (Gauss' method) {{anchor|th:GaussMethod}}:<!--\label{th:GaussMethod}-->
If a linear system is changed to another by one of these operations
<ol type=1 start=1>
<li>
an equation is swapped with another
<li>
an equation has both sides multiplied by a nonzero constant
<li>
an equation is replaced by the sum of itself and a multiple of another
</ol>
then the two systems have the same set of solutions.
}}

Each of those three operations has a restriction. Multiplying a row by <math>0</math> is not allowed because obviously that can change the solution set of the system. Similarly, adding a multiple of a row to itself is not allowed because adding <math>-1</math> times the row to itself has the effect of multiplying the row by <math>0</math>. Finally, swapping a row with itself is disallowed to make some results in the fourth chapter easier to state and remember (and besides, self-swapping doesn't accomplish anything).

{{TextBox|1=
;Proof:
We will cover the equation swap operation here and save the other two
cases for [[#ex:ProveGaussMethod|Problem 14]]<!--\ref{ex:ProveGaussMethod}-->.

Consider this swap of row <math>i</math> with row <math>j</math>.

:<math>
\begin{array}{*{4}{rc}r}
a_{1,1}x_1  &+  &a_{1,2}x_2 &+  &\cdots  &&a_{1,n}x_n  &=  &d_1  \\
&   &           &   &        &   &            &\vdots   \\
a_{i,1}x_1  &+  &a_{i,2}x_2 &+  &\cdots  &&a_{i,n}x_n  &=  &d_i  \\
&   &           &   &        &   &            &\vdots   \\
a_{j,1}x_1  &+  &a_{j,2}x_2 &+  &\cdots  &&a_{j,n}x_n  &=  &d_j  \\
&   &           &   &        &   &            &\vdots   \\
a_{m,1}x_1  &+  &a_{m,2}x_2 &+  &\cdots  &&a_{m,n}x_n  &=  &d_m  
\end{array}
\xrightarrow[]{}
\begin{array}{*{4}{rc}r}
a_{1,1}x_1  &+  &a_{1,2}x_2 &+  &\cdots  &&a_{1,n}x_n  &=  &d_1  \\
&   &           &   &        &   &            &\vdots   \\
a_{j,1}x_1  &+  &a_{j,2}x_2 &+  &\cdots  &&a_{j,n}x_n  &=  &d_j  \\
&   &           &   &        &   &            &\vdots   \\
a_{i,1}x_1  &+  &a_{i,2}x_2 &+  &\cdots  &&a_{i,n}x_n  &=  &d_i  \\
&   &           &   &        &   &            &\vdots   \\
a_{m,1}x_1  &+  &a_{m,2}x_2 &+  &\cdots  &&a_{m,n}x_n  &=  &d_m  
\end{array}
</math>


The <math>n</math>-tuple <math>(s_1,\ldots,s_n)</math> satisfies the system before the swap if and only if substituting the values, the <math>s</math>'s, for the variables, the <math>x</math>'s, gives true statements:
<math>a_{1,1}s_1+a_{1,2}s_2+\cdots+a_{1,n}s_n=d_1</math>
and ... 
<math>a_{i,1}s_1+a_{i,2}s_2+\cdots+a_{i,n}s_n=d_i</math>
and ...  <math>a_{j,1}s_1+a_{j,2}s_2+\cdots+a_{j,n}s_n=d_j</math>
and ...  <math>a_{m,1}s_1+a_{m,2}s_2+\cdots+a_{m,n}s_n=d_m</math>.

In a requirement consisting of statements and-ed together we can rearrange the order of the statements, so that this requirement is met if and only if
<math>a_{1,1}s_1+a_{1,2}s_2+\cdots+a_{1,n}s_n=d_1</math>
and ...  <math>a_{j,1}s_1+a_{j,2}s_2+\cdots+a_{j,n}s_n=d_j</math>
and ...  <math>a_{i,1}s_1+a_{i,2}s_2+\cdots+a_{i,n}s_n=d_i</math>
and ...  <math>a_{m,1}s_1+a_{m,2}s_2+\cdots+a_{m,n}s_n=d_m</math>.
This is exactly the requirement that <math>(s_1,\ldots,s_n)</math> solves the system after the row swap.
}}

{{TextBox|1=
;Definition 1.5:{{anchor|elementary operations}}
The three operations from [[#th:GaussMethod|Theorem 1.4]]<!--\ref{th:GaussMethod}-->
are the '''elementary reduction operations''',
or '''row operations''', or '''Gaussian operations'''.
They are '''swapping''', '''multiplying by a scalar''' or '''rescaling''', and '''pivoting'''.
}}

When writing out the calculations, we will abbreviate "row <math>i</math>" by "<math>\rho_i</math>". For instance, we will denote a pivot operation by <math>k\rho_i+\rho_j</math>, with the row that is changed written second. We will also, to save writing, often list pivot steps together when they use the same <math>\rho_i</math>.

{{TextBox|1=
;Example 1.6:
A typical use of Gauss' method is to solve this system.

:<math>
\begin{array}{*{3}{rc}r}
x  &+  &y  &   &   &=  &0  \\
2x  &-  &y  &+  &3z &=  &3  \\
x  &-  &2y &-  &z  &=  &3  
\end{array}
</math>

The first transformation of the system involves using the first row to eliminate the <math>x</math> in the second row and the <math>x</math> in the third. To get rid of the second row's <math>2x</math>, we multiply the entire first row by <math>-2</math>, 
add that to the second row, and write the result in as the new second row. To get rid of the third row's <math>x</math>, we multiply the first row by <math>-1</math>, add that to the third row, and write the result in as the new third row.

:<math>\begin{array}{rcl}
&\xrightarrow[-\rho_1+\rho_3]{-2\rho_1+\rho_2}
&\begin{array}{*{3}{rc}r}
x  &+  &y  &   &   &=  &0  \\
&   &-3y&+  &3z &=  &3  \\
&   &-3y&-  &z  &=  &3  
\end{array}   
\end{array}
</math>
(Note that the two <math>\rho_1</math> steps <math>-2\rho_1+\rho_2</math> and <math>-\rho_1+\rho_3</math> are written as one operation.) In this second system, the last two equations involve only two unknowns. To finish we transform the second system into a third system, where the last equation involves only one unknown. This transformation uses the second row to eliminate <math>y</math> from the third row.

:<math>\begin{array}{rcl}
&\xrightarrow[]{-\rho_2+\rho_3}
&\begin{array}{*{3}{rc}r}
x  &+  &y  &   &   &=  &0  \\
&   &-3y&+  &3z &=  &3  \\
&   &   &   &-4z&=  &0
\end{array}
\end{array}
</math>
Now we are set up for the solution. The third row shows that <math>z=0</math>. {{anchor|back substitution}}Substitute that back into the second row to get <math>y=-1</math>, and then substitute back into the first row to get <math> x=1 </math>.
}}

{{TextBox|1=
;Example 1.7{{anchor|statics problem}}:
For the Physics problem from the start of this
chapter, Gauss' method gives this.

:<math>\begin{array}{rcl}
\begin{array}{*{2}{rc}r}
40h  &+  &15c  &=  &100      \\
-50h &+  &25c  &=  &50         
\end{array}
&\xrightarrow[]{5/4\rho_1 +\rho_2}
&\begin{array}{*{2}{rc}r}
40h  &+  &15c       &=  &100      \\
&   &(175/4)c  &=  &175 
\end{array}
\end{array}
</math>
So <math> c=4 </math>, and back-substitution gives that <math> h=1 </math>.
(The Chemistry problem is solved later.)
}}

{{TextBox|1=
;Example 1.8:
The reduction

:<math>\begin{array}{rcl}
\begin{array}{*{3}{rc}r}
x  &+  &y  &+  &z  &=  &9  \\
2x  &+  &4y &-  &3z &=  &1  \\
3x  &+  &6y &-  &5z &=  &0  
\end{array}
&\xrightarrow[-3\rho_1 +\rho_3]{-2\rho_1 +\rho_2}
&\begin{array}{*{3}{rc}r}
x  &+  &y  &+  &z  &=  &9  \\
&   &2y &-  &5z &=  &-17\\
&   &3y &-  &8z&=  &-27
\end{array}                                    \\
&\xrightarrow[]{-(3/2)\rho_2+\rho_3}
&\begin{array}{*{3}{rc}r}
x  &+  &y  &+  &z            &=  &9  \\
&   &2y &-  &5z           &=  &-17\\
&   &   &   &-(1/2)z      &=  &-(3/2) 
\end{array}
\end{array}
</math>
shows that <math> z=3 </math>, <math> y=-1 </math>, and <math> x=7 </math>.
}}

As these examples illustrate, Gauss' method uses the elementary reduction
operations to set up back-substitution.

{{TextBox|1=
;Definition 1.9{{anchor|echelon form}}:
In each row, the first variable with a nonzero coefficient is the row's
'''leading variable'''. 
A system is in '''echelon form'''
if each leading variable
is to the right of the leading variable in the row above it
(except for the leading variable in the first row).
}}

{{TextBox|1=
;Example 1.10: 
The only operation needed in the examples above is pivoting.
Here is a linear system that requires the operation of swapping equations. 
After the first pivot

:<math>\begin{array}{rcl}
\begin{array}{*{4}{rc}r}
x  &-  &y  &   &   &   &   &=  &0  \\
2x  &-  &2y &+  &z  &+  &2w &=  &4  \\
&   &y  &   &   &+  &w  &=  &0  \\
&   &   &   &2z &+  &w  &=  &5  
\end{array}
&\xrightarrow[]{-2\rho_1 +\rho_2}
&\begin{array}{*{4}{rc}r}
x  &-  &y  &   &   &   &   &=  &0  \\
&   &   &   &z  &+  &2w &=  &4  \\
&   &y  &   &   &+  &w  &=  &0  \\
&   &   &   &2z &+  &w  &=  &5  
\end{array}    
\end{array}
</math>
the second equation has no leading <math>y</math>.
To get one, 
we look lower down in the system for a row that has a leading <math>y</math> and 
swap it in.

:<math>\begin{array}{rcl}
&\xrightarrow[]{\rho_2 \leftrightarrow\rho_3}
&\begin{array}{*{4}{rc}r}
x  &-  &y  &   &   &   &   &=  &0  \\
&   &y  &   &   &+  &w  &=  &0  \\
&   &   &   &z  &+  &2w &=  &4  \\
&   &   &   &2z &+  &w  &=  &5  
\end{array}    
\end{array}
</math>
(Had there been more than one row below the second with a leading <math>y</math>
then we could have swapped in any one.)
The rest of Gauss' method goes as before.

:<math>\begin{array}{rcl}
&\xrightarrow[]{-2\rho_3 +\rho_4}
&\begin{array}{*{4}{rc}r}
x  &-  &y  &   &   &   &   &=  &0  \\
&   &y  &   &   &+  &w  &=  &0  \\
&   &   &   &z  &+  &2w &=  &4  \\
&   &   &   &   &   &-3w&=  &-3 
\end{array}
\end{array}
</math>
Back-substitution gives <math> w=1 </math>, <math> z=2 </math>, <math> y=-1 </math>, and <math> x=-1 </math>.
}}

Strictly speaking, the operation of rescaling rows is 
not needed to solve linear systems.  
We have included it because we will use it 
later in this chapter as part of a variation on Gauss' method, 
the Gauss-Jordan method.

All of the systems seen so far have the same number of equations as unknowns.
All of them have a solution, and for all of them there is only one solution.
We finish this subsection by seeing for contrast
some other things that can happen.

{{TextBox|1=
;Example 1.11{{anchor|ex:MoreEqsThanUnks}}: <!--\label{ex:MoreEqsThanUnks}-->
Linear systems need not have the same number of equations as unknowns.
This system 

:<math>
\begin{array}{*{2}{rc}r}
x  &+  &3y  &=  &1  \\
2x  &+  &y   &=  &-3 \\
2x  &+  &2y  &=  &-2 
\end{array}  
</math>

has more equations than variables.
Gauss' method helps us understand this system also, since this

:<math>\begin{array}{rcl}
&\xrightarrow[-2\rho_1 +\rho_3]{-2\rho_1 +\rho_2}
&\begin{array}{*{2}{rc}r}
x  &+  &3y  &=  &1  \\
&   &-5y &=  &-5 \\
&   &-4y &=  &-4 
\end{array}
\end{array}
</math>
shows that one of the equations is redundant.
Echelon form

:<math>\begin{array}{rcl}
&\xrightarrow[]{-(4/5)\rho_2 +\rho_3}
&\begin{array}{*{2}{rc}r}
x  &+  &3y  &=  &1  \\
&   &-5y &=  &-5 \\
&   &0   &=  &0  
\end{array}
\end{array}
</math>
gives <math> y=1 </math> and <math> x=-2 </math>.
The "<math> 0=0 </math>" is derived from the redundancy.
}}

That example's system has more equations than variables.
Gauss' method is also useful on systems with more variables than equations.
Many examples are in the next subsection.

Another way that linear systems can differ from the examples shown earlier 
is that some linear systems do not have a unique solution.
This can happen in two ways.

The first is that it can fail to have any solution at all.

{{TextBox|1=
;Example 1.12{{anchor|ex:MoreEqsThanUnksInconsis}}: <!--\label{ex:MoreEqsThanUnksInconsis}-->
Contrast the system in the last example with this one.

:<math>\begin{array}{rcl}
\begin{array}{*{2}{rc}r}
x  &+  &3y  &=  &1  \\
2x  &+  &y   &=  &-3 \\
2x  &+  &2y  &=  &0  
\end{array}
&\xrightarrow[-2\rho_1 +\rho_3]{-2\rho_1 +\rho_2}
&\begin{array}{*{2}{rc}r}
x  &+  &3y  &=  &1  \\
&   &-5y &=  &-5 \\
&   &-4y &=  &-2
\end{array}
\end{array}
</math>
Here the system is inconsistent: no pair of numbers satisfies 
all of the equations simultaneously.
Echelon form makes this inconsistency obvious.

:<math>\begin{array}{rcl}
&\xrightarrow[]{-(4/5)\rho_2 +\rho_3}
&\begin{array}{*{2}{rc}r}
x  &+  &3y  &=  &1  \\
&   &-5y &=  &-5 \\
&   &0   &=  &2 
\end{array}
\end{array}
</math>
The solution set is empty.
}}

{{TextBox|1=
;Example 1.13:
The prior system has more equations than unknowns, but
that is not what causes the inconsistency&mdash; 
[[#ex:MoreEqsThanUnks|Example 1.11]]<!--\ref{ex:MoreEqsThanUnks}-->
has more equations than unknowns and yet is consistent.
Nor is having more equations than unknowns sufficient for
inconsistency, as is illustrated by this inconsistent system with the 
same number of equations as unknowns.

:<math>\begin{array}{rcl}
\begin{array}{*{2}{rc}r}
x  &+  &2y  &=  &8  \\
2x  &+  &4y  &=  &8  
\end{array}
&\xrightarrow[]{-2\rho_1 + \rho_2}
&\begin{array}{*{2}{rc}r}
x  &+  &2y  &=  &8  \\
&   &0   &=  &-8
\end{array}
\end{array}
</math>
}}

The other way that a linear system can fail to have a unique solution
is to have many solutions.

{{TextBox|1=
;Example 1.14:
In this system

:<math>
\begin{array}{*{2}{rc}r}
x  &+  &y   &=  &4  \\
2x  &+  &2y  &=  &8  
\end{array}
</math>

any pair of numbers satisfying the first equation automatically
satisfies the second.
The solution set <math> \{ (x,y)\,\big|\, x+y=4 \} </math> is infinite; some
of its members are <math>(0,4)</math>, <math>(-1,5)</math>, and <math>(2.5,1.5)</math>.
The result of applying Gauss' method here contrasts with the prior example
because we do not get a contradictory equation.

:<math>\begin{array}{rcl}
&\xrightarrow[]{-2\rho_1 + \rho_2}
&\begin{array}{*{2}{rc}r}
x  &+  &y   &=  &4  \\
&   &0   &=  &0  
\end{array}
\end{array}
</math>
}}

Don't be fooled by the "<math> 0=0 </math>" equation in that example.
It is not the signal that a system has many solutions.

{{TextBox|1=
;Example 1.15{{anchor|ex:NoZerosInfManySols}}: <!--\label{ex:NoZerosInfManySols}-->
The absence of a "<math> 0=0 </math>" does not keep a system from having
many different solutions.
This system is in echelon form

:<math>
\begin{array}{*{3}{rc}r}
x  &+  &y  &+  &z  &=  &0  \\
&   &y  &+  &z  &=  &0  
\end{array}
</math>

has no "<math>0=0</math>", and yet has infinitely many solutions.
(For instance, each of these is a solution: <math>(0,1,-1)</math>, 
<math>(0,1/2,-1/2)</math>, <math>(0,0,0)</math>, and <math>(0,-\pi,\pi)</math>.
There are infinitely many solutions because
any triple whose first component is <math>0</math> and whose second component is the 
negative of the third is a solution.)

Nor does the presence of a "<math> 0=0 </math>" mean that the system must have 
many solutions.
[[#ex:MoreEqsThanUnks|Example 1.11]]<!--\ref{ex:MoreEqsThanUnks}--> shows that.
So does this system, which does not have 
many solutions&mdash; in fact it has none&mdash; despite that 
when it is brought to echelon form it has a "<math>0=0</math>" row.

:<math>\begin{array}{rcl}
\begin{array}{*{3}{rc}r}
2x  &   &   &-   &2z  &=  &6  \\
&   &y  &+   &z   &=  &1  \\
2x  &+  &y  &-   &z   &=  &7  \\
&   &3y &+   &3z  &=  &0  
\end{array}
&\xrightarrow[]{-\rho_1 +\rho_3}
&\begin{array}{*{3}{rc}r}
2x  &   &   &-   &2z  &=  &6  \\
&   &y  &+   &z   &=  &1  \\
&   &y  &+   &z   &=  &1  \\
&   &3y &+   &3z  &=  &0  
\end{array}     \\
&\xrightarrow[-3\rho_2 +\rho_4]{-\rho_2 +\rho_3}
&\begin{array}{*{3}{rc}r}
2x  &   &   &-   &2z  &=  &6  \\
&   &y  &+   &z   &=  &1  \\
&   &   &    &0   &=  &0  \\
&   &   &    &0   &=  &-3 
\end{array}
\end{array}
</math>
}}

We will finish this subsection with a summary of
what we've seen so far about Gauss' method.

Gauss' method uses the three row operations to 
set a system up for back substitution.
If any step shows a contradictory equation then we can stop with the
conclusion that the system has no solutions.
If we reach echelon form without a contradictory equation,
and each variable is a leading variable in its
row, then the system has a unique solution and we find it by
back substitution.
Finally, if we reach echelon form without a contradictory equation,
and there is not a unique solution
(at least one variable is not a leading variable) then the system has many
solutions.

The next subsection deals with the third case&mdash; 
we will see how to describe the
solution set of a system with many solutions.

== Exercises ==
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 1: 
Use Gauss' method to find the unique solution for each system.
<ol type=1 start=1>
<li> 
<math>\begin{array}{*{2}{rc}r}
2x  &+  &3y  &=  &13  \\
x   &-  &y   &=  &-1
\end{array}</math>
<br><p><br><li> 
<math>\begin{array}{*{3}{rc}r}
x   &  &  &-  &z  &=  &0  \\
3x  &+ &y &   &   &=  &1  \\
-x  &+ &y &+  &z  &=  &4
\end{array}</math>
</ol>
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 2:  
Use Gauss' method to solve each system
or conclude "many solutions" or "no solutions".
<ol type=1 start=1>
<li> <math>
\begin{array}{*{2}{rc}r}
2x  &+  &2y  &=  &5  \\
x  &-  &4y  &=  &0  
\end{array}
</math>
<br><p><br>
<li> <math>
\begin{array}{*{2}{rc}r}
-x  &+  &y   &=  &1  \\
x  &+  &y   &=  &2  
\end{array}
</math> 
<br><p><br>
<li>  <math>
\begin{array}{*{3}{rc}r}
x  &-  &3y  &+  &z  &=  &1  \\
x  &+  &y   &+  &2z &=  &14 
\end{array}
</math> 
<br><p><br>
<li>  <math>
\begin{array}{*{2}{rc}r}
-x  &-  &y   &=  &1  \\
-3x  &-  &3y  &=  &2  
\end{array}
</math> 
<br><p><br>
<li>  <math>
\begin{array}{*{3}{rc}r}
&   &4y  &+  &z  &=  &20 \\
2x  &-  &2y  &+  &z  &=  &0  \\
x  &   &    &+  &z  &=  &5  \\
x  &+  &y   &-  &z  &=  &10 
\end{array}
</math>
<br><p><br>
<li> <math> \begin{array}{*{4}{rc}r}
2x  &   &   &+  &z  &+  &w  &=  &5  \\
&   &y  &   &   &-  &w  &=  &-1 \\
3x  &   &   &-  &z  &-  &w  &=  &0  \\
4x  &+  &y  &+  &2z &+  &w  &=  &9  
\end{array}
</math>
</ol>
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 3: 
There are methods for solving linear systems other 
than Gauss' method.
One often taught in high school is to solve one of the 
equations for a variable, then substitute the resulting expression into
other equations.
That step is repeated until there is an equation with only one
variable.
From that, the first number in the solution is derived, and then 
back-substitution can be done.
This method takes longer than Gauss' method, since it involves
more arithmetic operations, and is also more
likely to lead to errors.
To illustrate how it can lead to wrong conclusions, we will use the system 

:<math>
\begin{array}{*{2}{rc}r}
x  &+  &3y  &=  &1  \\
2x  &+  &y   &=  &-3 \\
2x  &+  &2y  &=  &0  
\end{array}
</math>

from [[#ex:MoreEqsThanUnksInconsis|Example 1.12]]<!--\ref{ex:MoreEqsThanUnksInconsis}-->.
<ol type=1 start=1>
<li> Solve the first equation for <math>x</math> and 
substitute that expression into the second equation.
Find the resulting <math>y</math>.
<li> Again solve the first equation for <math>x</math>, 
but this time substitute that expression into the third equation.
Find this <math>y</math>.
</ol>
What extra step must a user of this method take to avoid 
erroneously concluding a system has a solution?
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 4: 
For which values of <math> k </math> are
there no solutions, many solutions, or a unique solution
to this system?

:<math>
\begin{array}{*{2}{rc}r}
x  &-  &y  &=  &1  \\
3x  &-  &3y &=  &k  
\end{array}
</math>
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 5: 
This system is not linear, in some sense,

:<math>
\begin{array}{*{3}{rc}r}
2\sin\alpha  &-  &\cos\beta  &+  &3\tan\gamma  &=  &3  \\
4\sin\alpha  &+  &2\cos\beta &-  &2\tan\gamma  &=  &10  \\
6\sin\alpha  &-  &3\cos\beta &+  &\tan\gamma   &=  &9  
\end{array}
</math>

and yet we can nonetheless apply Gauss' method.
Do so. Does the system have a solution?
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 6: 
What conditions must the constants, the <math>b</math>'s,
satisfy so that each of these systems has a solution?
''Hint.'' 
Apply Gauss' method and see what happens to the right side {{harv|Anton|1987}}.

<ol type=1 start=1>
<li>  <math>
\begin{array}{*{2}{rc}r}
x  &-  &3y  &=  &b_1 \\
3x  &+  &y   &=  &b_2 \\
x  &+  &7y  &=  &b_3 \\
2x  &+  &4y  &=  &b_4 
\end{array}   </math>
<br><p><br>
<li> <math>
\begin{array}{*{3}{rc}r}
x_1  &+  &2x_2  &+  &3x_3  &=  &b_1  \\
2x_1  &+  &5x_2  &+  &3x_3  &=  &b_2  \\
x_1  &   &      &+  &8x_3  &=  &b_3  
\end{array}   </math>
</ol>
}}
{{TextBox|1=
;Problem 7: 
True or false: a system with more unknowns than equations
has at least one solution.
(As always, to say "true" you must prove it, while to say 
"false" you must produce a counterexample.)
}}
{{TextBox|1=
;Problem 8{{anchor|chemistry problem}}: 
Must any Chemistry problem like
the one that starts this subsection&mdash;
a balance the reaction problem&mdash; have infinitely many solutions?
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 9: 
Find the coefficients
<math> a </math>, <math> b </math>, and <math> c </math> so that the graph of <math> f(x)=ax^2+bx+c </math> 
passes through the points <math> (1,2) </math>, <math> (-1,6) </math>, and <math> (2,3) </math>.
}}
{{TextBox|1=
;Problem 10: 
Gauss' method works by combining the equations in a system to make new
equations.
<ol type=1 start=1>
<li> Can the equation <math> 3x-2y=5 </math> be derived, by a sequence of
Gaussian reduction steps, from the equations in this system?

:<math>
\begin{array}{*{2}{rc}r}
x  &+  &y  &=  &1  \\
4x  &-  &y  &=  &6
\end{array}
</math>

<li> Can the equation <math> 5x-3y=2 </math> be derived, by a sequence of
Gaussian reduction steps, from the equations in this system?

:<math>
\begin{array}{*{2}{rc}r}
2x  &+  &2y &=  &5  \\
3x  &+  &y  &=  &4
\end{array}
</math>

<li> Can the equation <math> 6x-9y+5z=-2 </math> be derived, 
by a sequence of
Gaussian reduction steps, from the equations in the system?

:<math>
\begin{array}{*{3}{rc}r}
2x  &+  &y  &-  &z  &=  &4  \\
6x  &-  &3y &+  &z  &=  &5
\end{array}
</math>

</ol>
}}
{{TextBox|1=
;Problem 11: 
Prove that, where <math> a,b,\ldots,e </math> are real numbers
and <math> a\neq 0 </math>, if

:<math>
ax+by=c
</math>

has the same solution set as

:<math>
ax+dy=e
</math>

then they are the same equation.
What if <math> a=0 </math>?
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 12: 
Show that if <math> ad-bc\neq 0 </math> then

:<math>
\begin{array}{*{2}{rc}r}
ax  &+  &by  &=  &j  \\
cx  &+  &dy  &=  &k  
\end{array}
</math>

has a unique solution.
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 13: 
In the system

:<math>
\begin{array}{*{2}{rc}r}
ax  &+  &by  &=  &c  \\
dx  &+  &ey  &=  &f  
\end{array}
</math>

each of the equations describes a line in the <math> xy </math>-plane.
By geometrical reasoning, show that there are three possibilities:
there is a unique solution, there is no solution, 
and there are infinitely many solutions.
}}
{{TextBox|1=
;Problem 14 {{anchor|ex:ProveGaussMethod}}: <!--\label{ex:ProveGaussMethod}-->
Finish the proof of [[#th:GaussMethod|Theorem 1.4]]<!--\ref{th:GaussMethod}-->.
}}
{{TextBox|1=
;Problem 15: 
Is there a two-unknowns
linear system whose solution set is all of <math> \mathbb{R}^2 </math>?
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;Problem 16: 
Are any of the operations used in Gauss' method
redundant?
That is, can any of the operations be synthesized from the others?
}}
{{TextBox|1=
;Problem 17: 
Prove that each operation of Gauss' method is reversible.
That is, show that if two systems are related by a row operation
<math>S_1\rightarrow S_2</math> then there is a row operation to go back
<math>S_2\rightarrow S_1</math>.
}}
{{TextBox|1=
;? Problem 18:  
A box holding pennies, nickels and dimes contains
thirteen coins with a total value of <math> 83 </math> cents.
How many coins of each type are in the box?
{{harv|Anton|1987}}
}}
{{TextBox|1=
;? Problem 19: 
Four positive integers are given.
Select any three of the integers, find their arithmetic average,
and add this result to the fourth integer.
Thus the numbers 29, 23, 21, and 17 are obtained.
One of the original integers is:
<ol type=1 start=1>
<li> 19
<li> 21
<li> 23
<li> 29
<li> 17
</ol>
{{harv|Salkind|1975|loc=1955 problem 38}}
}}
{{Linear Algebra/Book 2/Recommended}}
{{TextBox|1=
;? Problem 20:  
Laugh at this:  <math> \mbox{AHAHA}+\mbox{TEHE}=\mbox{TEHAW} </math>.
It resulted from substituting a code letter for each digit of a simple
example in addition, and it is required to identify the letters
and prove the solution unique {{harv|Ransom|Gupta|1935}}.
}}
{{TextBox|1=
;? Problem 21: 
The Wohascum County Board of Commissioners, which has 20 members, 
recently had to elect a President.
There were three candidates (<math>A</math>, <math>B</math>, and <math>C</math>); on each ballot
the three
candidates were to be listed in order of preference, with no abstentions.
It was found that 11 members, a majority, preferred <math>A</math> over <math>B</math>
(thus the other 9 preferred <math>B</math> over <math>A</math>).
Similarly, it was found that 12 members preferred <math>C</math> over <math>A</math>.
Given these results, it was suggested that <math>B</math> should withdraw, to enable
a runoff election between <math>A</math> and <math>C</math>.
However, <math>B</math> protested, and it was then found that 14 members preferred
<math>B</math> over <math>C</math>!
The Board has not yet recovered from the resulting confusion.
Given that every possible order of <math>A</math>, <math>B</math>, <math>C</math> appeared on at least 
one ballot, how many members voted for <math>B</math> as their first choice {{harv|Gilbert|Krusemeyer|Larson|1993|loc=Problem number 2}}?

}}
{{TextBox|1=
;? Problem 22:   
"This system of <math>n</math> linear equations with
<math>n</math> unknowns," said the Great Mathematician, "has a curious
property."

"Good heavens!" said the Poor Nut,  "What is it?"

"Note," said the Great Mathematician, "that the constants are in
arithmetic progression."

"It's all so clear when you explain it!" said the Poor Nut.
"Do you mean like <math> 6x+9y=12 </math> and <math> 15x+18y=21 </math>?"

"Quite so," said the Great Mathematician, pulling out his bassoon.
"Indeed, even larger systems can be solved regardless of their progression.
Can you find their solution?"

"Good heavens!" cried the Poor Nut, "I am baffled."

Are you? {{harv|Dudley|Lebow|Rothman|1963}}
}}
<noinclude>
[[/Solutions/]]



{{navigation|Book=线性代数（英文维基教科书）
|current=Gauss' Method 高斯消元法
|previous=Solving Linear Systems 求解线性方程组
|next=Describing the Solution Set 解集的表示}}