本讲的主要内容:
-
Ax=b 的求解过程
- 讨论
Ax=b 各种情况是否存在解
求解过程
这一部分使用的例子,这里我直接写成矩阵形式:
Ax=b⇔⎝⎛1232462682810⎠⎞x=⎝⎛b1b2b3⎠⎞
接下来使用写成增广矩阵的形式
(A,b) 进行消元:
⎝⎛1232462682810b1b2b3⎠⎞⇒⎝⎛100200220240b1b2−2b1b3−b2−b1⎠⎞
至此,我们可以很清楚的看到这个方程可解的条件:
b3−b2−b1=0
对于方程
Ax=b 方程可解的条件,或者说可解性(Solvability) 即(其实在之前列空间部分讲到过):
-
b 存在于
A 的列空间中,也就是
b∈C(A)
或者说:
-
A 的行向量的线性组合为0的时候,右侧向量的组合也必须为0,也就是同步
那么,接下来来写出方程的完整解:
分为两个步骤:
- 1.找到方程的特殊解(particular solution)
- 2.找到方程的零空间的任意解(也就是基础解系)
为方便,取
⎝⎛b1b2b3⎠⎞=⎝⎛156⎠⎞代入方程中,得到消元后增广矩阵:
⎝⎛100200220240130⎠⎞
按照求解的算法:
Step1: 求出特定解,也就是将所有的自由变量(free variables)赋值为0,接出来得到:
xparticular=⎝⎜⎜⎛−203/20⎠⎟⎟⎞
Step2: 求出零空间的任意一个解,也就是求解对应的
Ax=0 的解,具体做法按照之前的算法,将所有的自由变量依次赋值为1,上述方程有两个自由变量,最终结果可以记作:
xnullspace=c1⎝⎜⎜⎛−2100⎠⎟⎟⎞+c2⎝⎜⎜⎛20−20⎠⎟⎟⎞,这个结果又叫做基础解系
Step3: 将上述两部分加起来即可得到完整解,也就是:
xcompelete=xparticular+xnullspace=⎝⎜⎜⎛−203/20⎠⎟⎟⎞+c1⎝⎜⎜⎛−2100⎠⎟⎟⎞+c2⎝⎜⎜⎛20−20⎠⎟⎟⎞
这个结果即为方程的所有解,这个结果证明如下:
{Axparticular=bAxnullspace=0⇒A(xparticular+xnullspace)=b
如果从图像中理解这个结果,也就是在四维空间中的一张平面(不太直观):我们可以看到基础解系确定的是这个四维空间的一个过原点的平面,也就是
Ax=0 的情况,将这个平面移动到过
xparticular 这个点,此时,这个平面就表示结果。
讨论不同情况方程的可解性
总结对于一个
m×n 的矩阵(秩为
r )所有的解的可能情况:
- 1.
m=n=r: 唯一解,没有自由变量,所以零空间只有零向量,所以解就等于特殊解,最后的形式为
- 2.
r=m<n: 有无穷多解,有自由变量,零空间的存在无穷多向量。
- 3.
r=n<m: 有唯一解或者只有0解,没有自由变量,所以零空间只有一个零向量,只存在特殊解,特殊解可能存在或者不存在。
- 4.
r<m,r<n: 无解或者有无穷多解
这几个结论,其实只要画一下矩阵的大致形状就很好理解了。
以上~