上一讲介绍了列空间和零空间,接下来我们开始讲解如何找出这些空间中的向量。
1. 求解 A X = 0 AX=0 AX=0
1.1 秩
A
=
[
1
2
2
2
2
4
6
8
3
6
8
10
]
(1)
A=\begin{bmatrix} 1&2&2&2\\ 2&4&6&8\\ 3&6&8&10\\ \end{bmatrix} \tag{1}
A=
1232462682810
(1)
求解
A
X
=
0
AX=0
AX=0,同样用到消元法,此过程不会改变解,右侧向量也始终为
0
0
0。
A
=
[
1
2
2
2
2
4
6
8
3
6
8
10
]
→
[
1
2
2
2
0
0
2
4
0
0
2
4
]
A=\begin{bmatrix} \boldsymbol{1}&2&2&2\\ 2&4&6&8\\ 3&6&8&10 \end{bmatrix}\to \begin{bmatrix} \boldsymbol{1}&2&2&2\\ 0&0&2&4\\ 0&0&2&4 \end{bmatrix}
A=
1232462682810
→
100200222244
此时主元2的位置都为0,意味着列1和列2是线性相关的,那么我们跳过列2,从列3继续消元。
A → [ 1 2 2 2 0 0 2 4 0 0 0 0 ] = U (2) A\to \begin{bmatrix} \boldsymbol{1}&2&2&2\\ 0&0&\boldsymbol{2}&4\\ 0&0&0&0 \end{bmatrix}=U\tag{2} A→ 100200220240 =U(2)
消元成上三角矩阵后,主元位置以阶梯的形式呈现,我们称这种形式为阶梯形式。
并且其主元只有
2
2
2个,该数字称为矩阵的秩(rank)。
由此,我们给出矩阵秩的定义:矩阵 A A A的秩(rank),即主元的个数。
由于消元过程不会影响解和零空间,因此,求解 A X = 0 AX=0 AX=0变为求解 U X = 0 UX=0 UX=0。
对于(2)式中的上三角矩阵
U
=
[
1
2
2
2
0
0
2
4
0
0
0
0
]
U=\begin{bmatrix} \boldsymbol{1}&2&2&2\\ 0&0&\boldsymbol{2}&4\\ 0&0&0&0 \end{bmatrix}
U=
100200220240
将列1和列3称为主列(主元列),也就是主元所在的列,
列2和列4称为自由列,即自由变量所在的列。
因此, x 2 x_2 x2和 x 4 x_4 x4可以是任意数,我们只需解出 x 1 x_1 x1和 x 3 x_3 x3:
假设 X = [ x 1 1 x 3 0 ] X = \begin{bmatrix} x_1 \\ 1 \\ x_3 \\ 0 \end{bmatrix} X= x11x30 ,将 U U U 写成方程组形式: { x 1 + 2 x 2 + 2 x 3 + 2 x 4 = 0 2 x 3 + 4 x 4 = 0 \begin{cases} x_1 + 2x_2 + 2x_3 + 2x_4 = 0 \\ 2x_3 + 4x_4 = 0 \end{cases} {x1+2x2+2x3+2x4=02x3+4x4=0
回代 x 2 = 1 x_2 = 1 x2=1, x 4 = 0 x_4 = 0 x4=0,得到 { x 1 + 2 x 3 + 2 = 0 2 x 3 = 0 \begin{cases} x_1 + 2x_3 + 2 = 0 \\ 2x_3 = 0 \end{cases} {x1+2x3+2=02x3=0
解得 X 1 = [ − 2 1 0 0 ] X_1 = \begin{bmatrix} -2 \\ 1 \\ 0 \\ 0 \end{bmatrix} X1= −2100 ,此即为 A X = 0 AX = 0 AX=0的一个特解。
同样地,将 X = [ x 1 0 x 3 1 ] X= \begin{bmatrix} x_1\\ 0\\ x_3\\ 1\\ \end{bmatrix} X= x10x31 ,回代后,得到 X 2 = [ 2 0 − 2 1 ] X_2= \begin{bmatrix} 2\\0\\-2\\1 \end{bmatrix} X2= 20−21
这样我们就找到了该零空间中另一个向量。
将
X
1
,
X
2
X_1,X_2
X1,X2称为特解,那么整个零空间便是这些特解的线性组合,
即
X
=
C
1
[
−
2
1
0
0
]
+
C
2
[
2
0
−
2
1
]
X=C_1 \begin{bmatrix} -2\\1\\0\\0\\ \end{bmatrix} + C_2 \begin{bmatrix} 2\\0\\-2\\1 \end{bmatrix}
X=C1
−2100
+C2
20−21
。
通过以上求解过程,我们有:
对于 m × n m\times n m×n的矩阵 A A A,其有 n n n个未知数,若主元有 r r r个,即秩为 r r r,对应地自由变量有 ( n − r ) (n-r) (n−r)个,那么求解时有 r r r个方程起作用。
1.2 行简化
对于阶梯形式的矩阵 U U U,我们还可以对其进一步简化成简化行阶梯形式,即行最简形矩阵 R R R。
U
=
[
1
2
2
2
0
0
2
4
0
0
0
0
]
→
[
1
2
0
−
2
0
0
2
4
0
0
0
0
]
→
[
1
2
0
−
2
0
0
1
2
0
0
0
0
]
=
R
U=\begin{bmatrix} \boldsymbol{1}&2&2&2\\ 0&0&\boldsymbol{2}&4\\ 0&0&0&0 \end{bmatrix}\to \begin{bmatrix} \boldsymbol{1}&2&0& - 2\\ 0&0&\boldsymbol{2}&4\\ 0&0&0&0 \end{bmatrix}\to \begin{bmatrix} \boldsymbol{1}&2&0& - 2\\ 0&0&\boldsymbol{1}&2\\ 0&0&0&0 \end{bmatrix}=R
U=
100200220240
→
100200020−240
→
100200010−220
=R
对(2)式中的U,我们继续向上进行消元,使得主元均变为1,就得到行最简形矩阵
R
R
R。
求解
R
X
=
0
RX=0
RX=0就是
A
X
=
0
AX=0
AX=0的最简化形式。
将主元列和自由列分开,得到:
主元列
I
=
[
1
0
0
1
]
,自由列
F
=
[
2
−
2
0
2
]
主元列 I=\begin{bmatrix}1&0\\0&1\end{bmatrix},自由列 F = \begin{bmatrix}2& - 2\\0&2\end{bmatrix}
主元列I=[1001],自由列F=[20−22]
对于解 X = C 1 [ − 2 1 0 0 ] + C 2 [ 2 0 − 2 1 ] X = C_1\begin{bmatrix}-2\\1\\0\\0\end{bmatrix}+C_2\begin{bmatrix}2\\0\\-2\\1\end{bmatrix} X=C1 −2100 +C2 20−21 ,因为对于给定的自由变量值,要使等式 A X = 0 AX=0 AX=0成立,即移到等式另一侧后是相反的,所以有如下关系:
主变量解 [ − 2 2 0 − 2 ] ,自由变量值 [ 1 0 0 1 ] 主变量解 \begin{bmatrix}-2&2\\0& - 2\end{bmatrix},自由变量值 \begin{bmatrix}1&0\\0&1\end{bmatrix} 主变量解[−202−2],自由变量值[1001]。
如何直接求解 A X = 0 AX = 0 AX=0:
构造零空间矩阵 N N N,其各列由特解组成,使得 R N = 0 RN = 0 RN=0。
若 R = [ I F 0 0 ] R=\begin{bmatrix}I&F\\0&0\end{bmatrix} R=[I0F0],则 N = [ − F I ] N=\begin{bmatrix}-F\\I\end{bmatrix} N=[−FI],即 R R R的主元部分变为 N N N的自由部分,自由部分取反变为 N N N的主元部分。
对此的解释是:
因为
R
X
=
0
RX = 0
RX=0,其中
R
=
[
I
F
]
,
X
=
[
x
主元
x
自由
]
R=\begin{bmatrix}I&F\end{bmatrix},X=\begin{bmatrix}x_{主元}\\x_{自由}\end{bmatrix}
R=[IF],X=[x主元x自由],即
[
I
F
]
[
x
主元
x
自由
]
=
0
\begin{bmatrix}I&F\end{bmatrix}\begin{bmatrix}x_{主元}\\x_{自由}\end{bmatrix}=0
[IF][x主元x自由]=0。
可看作主元列与 F F F乘以自由变量的和为 0 0 0,即 x 主元 + F x 自由 = 0 x_{主元}+F x_{自由}=0 x主元+Fx自由=0,那么 x 主元 = − F x 自由 x_{主元}=-F x_{自由} x主元=−Fx自由。
这时给自由变量分配单位矩阵 I I I,当自由变量部分对应单位矩阵时, x 主元 = − F x_{主元}=-F x主元=−F。
下面给出另一示例:
A
=
[
1
2
3
2
4
6
2
6
8
2
8
10
]
→
[
1
2
3
0
0
0
0
2
2
0
4
4
]
→
[
1
2
3
0
2
2
0
0
0
0
0
0
]
=
U
A=\begin{bmatrix} 1&2&3\\ 2&4&6\\ 2&6&8\\ 2&8&10 \end{bmatrix}\to \begin{bmatrix} \boldsymbol{1}&2&3\\ 0&0&0\\ 0&2&2\\ 0&4&4 \end{bmatrix}\to \begin{bmatrix} \boldsymbol{1}&2&3\\ 0&\boldsymbol{2}&2\\ 0&0&0\\ 0&0&0 \end{bmatrix}=U
A=
1222246836810
→
100020243024
→
100022003200
=U
主元列
2
2
2 个,自由列
1
1
1个(即自由变量
1
1
1个)。
写作方程组:
{
x
1
+
2
x
2
+
3
x
3
=
0
2
x
2
+
2
x
3
=
0
\begin{cases} x_1 + 2x_2+3x_3=0\\ 2x_2+2x_3=0 \end{cases}
{x1+2x2+3x3=02x2+2x3=0
解出
X
1
=
[
−
1
−
1
1
]
X_1=\begin{bmatrix}-1\\-1\\1\end{bmatrix}
X1=
−1−11
,通解
X
=
C
[
−
1
−
1
1
]
X = C\begin{bmatrix}-1\\-1\\1\end{bmatrix}
X=C
−1−11
。
上述为一般的方法。
若继续求
R
R
R:
R
=
[
1
0
1
0
1
1
0
0
0
0
0
0
]
R=\begin{bmatrix} \boldsymbol{1}&0&1\\ 0&\boldsymbol{1}&1\\ 0&0&0\\ 0&0&0 \end{bmatrix}
R=
100001001100
分开主元部分和自由部分:
主元部分
[
1
0
0
1
]
,自由部分
[
1
1
]
主元部分\begin{bmatrix}1&0\\0&1\end{bmatrix},自由部分 \begin{bmatrix}1\\1\end{bmatrix}
主元部分[1001],自由部分[11]。
若自由变量分配单位矩阵,主变量就是自由部分取反,得
N
=
[
−
1
−
1
1
]
N=\begin{bmatrix}-1\\-1\\1\end{bmatrix}
N=
−1−11
其中矩阵 N N N各列由特解组成,对其特解线性组合即得到通解。
总结:
对于求解 A X = 0 AX=0 AX=0,我们只需将 A A A化简为 R R R,主变量为 R R R中自由部分取反,自由变量分配单位矩阵,就能得到特解矩阵 N N N, N N N的各列即 A X = 0 AX=0 AX=0的特解,这些特解的线性组合就是 X X X。