4.8 行满秩方程

4.8 行满秩方程

高斯约当消元法还可以应用于行满秩矩阵，行满秩矩阵进行高斯约当消元法，最终矩阵变为单位矩阵和自由矩阵。例如方程

$$\begin{gathered} 2x+4y+6z=12 \\ 4x+9y+13z=36 \end{gathered}$$

系数矩阵为
$$A=\left[\begin{array}{ccc}2& 4& 6\\ 4& 9& 13\end{array}\right]$$
是行满秩矩阵。

方程有 $3$ 个未知数，但是只有 $2$ 个方程，有无穷多解。

增广矩阵进行高斯消元法

$$\left[\begin{array}{cccc}2& 4& 6& 12\\ 4& 9& 13& 36\end{array}\right]\Rightarrow \left[\begin{array}{cccc}2& 4& 6& 12\\ 0& 1& 1& 12\end{array}\right]$$

方程变为
$$\begin{gathered} 2x+4y+6z=12 \\ y+z=12 \end{gathered}$$

未知数 $z$ 是自由变量，可取任意实数，方程解可以表示为
$$\begin{gathered} x=-18-z \\ y=12-z \\ z=0＋z \end{gathered}$$

写成向量形式

$$\left[\begin{array}{c}x\\ y\\ z\end{array}\right]=\left[\begin{array}{c}-18\\ 12\\ 0\end{array}\right]+z\left[\begin{array}{c}-1\\ -1\\ 1\end{array}\right]$$

$z$ 取任意实数。

这是只有一个自由变量的情况，如果有多个呢？举例说明。
$$\begin{gathered} 2x+4y+6z+w=12 \\ 4x+9y+13z+4w=36 \end{gathered}$$

系数矩阵为
$$A=\left[\begin{array}{cccc}2& 4& 6& 1\\ 4& 9& 13& 4\end{array}\right]$$
是行满秩矩阵。

方程有 $4$ 个未知数，但是只有 $2$ 个方程，有无穷多解。

增广矩阵进行高斯消元法

$$\left[\begin{array}{ccccc}2& 4& 6& 1& 12\\ 4& 9& 13& 4& 36\end{array}\right]\Rightarrow \left[\begin{array}{ccccc}2& 4& 6& 1& 12\\ 0& 1& 1& 2& 12\end{array}\right]$$

方程变为
$$\begin{gathered} 2x+4y+6z+w=12 \\ y+z+2w=12 \end{gathered}$$

未知数 $z,w$ 是自由变量，可取任意实数，方程解可以表示为
$$\begin{gathered} x=-18-z-w \\ y=12-z-2w \\ z=0＋z+0w \\ w=0+0z+w \end{gathered}$$

写成向量形式

$$\left[\begin{array}{c}x\\ y\\ z\\ w\end{array}\right]=\left[\begin{array}{c}-18\\ 12\\ 0\\ 0\end{array}\right]+z\left[\begin{array}{c}-1\\ -1\\ 1\\ 0\end{array}\right]+w\left[\begin{array}{c}-1\\ -2\\ 0\\ 1\end{array}\right]$$

$z,w$ 取任意实数。

有几点需要说明，第 $1$ 点，自由变量不一定非取 $z,w$ ，有多种取法，比如可取 $y,w$ 。

方程为

$$\begin{gathered} 2x+4y+6z+w=12 \\ y+z+2w=12 \end{gathered}$$

未知数 $y,w$ 是自由变量，可取任意实数，方程解可以表示为
$$\begin{gathered} x=-60+2y+11w \\ z=12-y-2w \\ y=0+y+0w \\ w=0+0y+w \end{gathered}$$

写成向量形式

$$\left[\begin{array}{c}x\\ z\\ y\\ w\end{array}\right]=\left[\begin{array}{c}-60\\ 12\\ 0\\ 0\end{array}\right]+y\left[\begin{array}{c}2\\ -1\\ 1\\ 0\end{array}\right]+w\left[\begin{array}{c}11\\ -2\\ 0\\ 1\end{array}\right]$$

$y,w$ 取任意实数。

这两组解的形式虽然差别很大，但其实都是方程的解。

第 $2$ 点，自由变量虽然有多种取法，但自由变量的数量是固定的，为 $n-m$ ，是矩阵 $A$ 零空间的维度。

第 $3$ 点，向量 $\left[\begin{array}{c}2\\ -1\\ 1\\ 0\end{array}\right]$ 和向量 $\left[\begin{array}{c}11\\ -2\\ 0\\ 1\end{array}\right]$ 是方程 $A\mathbf{x}=0$ 解，称为零解，它们线性组合构成矩阵 $A$ 零空间。一般来说，每个零解中自由变量组均取单位向量 ${\mathbf{e}}_i$ ，即只有一个自由变量取 $1$ ，其它均为 $0$ ，这样计算比较方便。

第 $4$ 点，向量 $\left[\begin{array}{c}-60\\ 12\\ 0\\ 0\end{array}\right]$ 是方程 $A\mathbf{x}=\mathbf{b}$ 解，称为特解。一般来说，特解中自由变量均取 $0$ ，这样计算比较方便，当然也可取任意值。

所以行满秩方程解的结构为：特解 $+$ 零解的线性组合。

第 $5$ 点，非自由变量对应的列向量，是矩阵列空间的基。本例中，$\left[\begin{array}{cc}2& 4\\ 4& 9\end{array}\right]$ 这两个列向量是列空间的基。由于自由变量有多种取法，故基也有多种。

前面都是非自由变量对应的列向量在矩阵的最前列，即最前面的上三角阵是可逆的，但有时不可逆。例如方程

$$\begin{gathered} 2x+y+3z+2w=5 \\ 4x+2y+7z+5w=12 \end{gathered}$$

增广矩阵进行高斯消元法

$$\left[\begin{array}{ccccc}2& 1& 3& 2& 5\\ 4& 2& 7& 5& 12\end{array}\right]\Rightarrow \left[\begin{array}{ccccc}2& 1& 3& 2& 5\\ 0& 0& 1& 1& 2\end{array}\right]$$

前面两列是矩阵
$$\left[\begin{array}{cc}2& 1\\ 0& 0\end{array}\right]$$

对角元素有 $0$ ，上三角阵不可逆。对应方程为

$$\begin{gathered} 2x+y+3z+2w=5 \\ 0x+0y+z+w=2 \end{gathered}$$

那怎么选非自由变量呢？我们把非自由变量称为主元。每行选一个，从左到右选首个非零元素对应的变量为主元。第一行首个非零元素是 $2$ ，对于主元是 $x$ ，第二行首个非零元素是 $1$ ，对于主元是 $z$ ；则自由变量就是 $y,w$ ，可取任意实数，方程解可以表示为
$$\begin{gathered} x=-1-y+w \\ z=2-0y-w \\ y=0+y+0w \\ w=0+0y+w \end{gathered}$$

写成向量形式

$$\left[\begin{array}{c}x\\ z\\ y\\ w\end{array}\right]=\left[\begin{array}{c}-1\\ 2\\ 0\\ 0\end{array}\right]+y\left[\begin{array}{c}-1\\ 0\\ 1\\ 0\end{array}\right]+w\left[\begin{array}{c}1\\ -1\\ 0\\ 1\end{array}\right]$$

$y,w$ 取任意实数。

矩阵 $A$ 经过高斯消元变换为 $\left[\begin{array}{cccc}2& 1& 3& 2\\ 0& 0& 1& 1\end{array}\right]$ ，但我们希望前两列向量是上三角阵，这样更有利于理论分析和对比。这可通过列对调操作实现，即对调 $2,3$ 列向量，则前两列向量是上三角阵。列对调操作类似行对调操作，可以通过置换矩阵实现，差别在于：行对调是左乘，列对调是右乘，即 $P_{ij}A$ 对调矩阵 $i,j$ 两行， $A{P}_{ij}$ 对调矩阵 $i,j$ 两列。

总结如下，行满秩矩阵 $A_{mn}$，高斯消元法变换为 $\left[\begin{array}{c}{U}_{mm},F_{m,n-m}\end{array}\right]$ ，$U_{mm}$ 是 $m$ 阶上三角阵，其对角元素是矩阵 $A$ 的主元，$F_{m,n-m}$ 是自由矩阵。用矩阵乘法表示为 $L_{mm}A=\left[\begin{array}{c}{U}_{mm},F_{m,n-m}\end{array}\right]$ ，$L_{mm}$ 是 $m$ 阶单位下三角阵。上三角阵对应的变量是主元，自由矩阵对应的变量是自由变量。方程解结构为：特解加零解。如有必要，需进行列对调。矩阵乘法表示，即对任意行满秩矩阵 $A$ ，存在可逆矩阵 $P,Q$ ，使 $PAQ=\left[\begin{array}{c}{U}_{mm},F_{m,n-m}\end{array}\right]$ 成立，进一步对 $U_{mm}$ 进行高斯约当消元，则可表示为，存在可逆矩阵 $P,Q$ ，使 $PAQ=\left[\begin{array}{c}{E}_{mm},F_{m,n-m}\end{array}\right]$ 成立。

前面章节介绍了，行满秩矩阵的列向量组是相关组，其极大无关组是 $m$ 维空间的基。利用高斯消元法可以找到行满秩矩阵的极大无关组，主元对应的列向量即是极大无关组。