Gilbert Strang第八课时:求解Ax=b:可解性和解的结构

本文是Gilbert Strang的线性代数导论课程笔记。课程地址: http://v.163.com/special/opencourse/daishu.html  
第八课时:求解Ax=b:可解性和解的结构
本课时的目标是Ax=b,可能有解,也可能无解,需要通过需要消元才知道,有解的话是唯一解还是很多解。

继续用上课时的例子。
注意到,方程组中,第三行是第一行和第二行的和。如果方程组有解,b1 b2 b3需要满足什么条件?必须满足b3=b1+b2。消元告诉我们,这是必须的。换句话说,左侧行的线性组合得到0,那么右侧常量线性组合也比为0。

第一阶段消元,增广矩阵Augmented matrix=[A b]:(哈哈,上个教授的图)
第二阶段消元,可以看到最后一个方程是什么,常量0=b3-b2-b1,验证了打头说的必须满足的有解的条件:b3=b1+b2
我们需要求解的是Ax=b中的x,现在b我们可以根据b3=b1+b2来选定一组值,假设b=[1 5 6]。现在只有两个方程,4个未知数。

Ax=b可解性Solvability:有解时右侧向量b须满足的条件
1)有解,仅当b属于A的列空间时成立,即,b必须是A的各列的线性组合,
2)行的线性组合如果得到零行,那么b中元素的同样组合必然也是零。
这两种描述是等价的!他们同样是描述方程组有解的条件。

求解Ax=b的解
1)特解,将所有的自由量设置为0,然后解出主变量得到特解Xp
2)零空间中的任意x,Xn
因此Ax=b的所有解为特解加上零空间中任意向量。  
对于方程组某解,其与零空间内任意向量之和仍为解,因为零空间右侧得到的是0。
如此,就可以得到上式方程组的所有解:
 
把所有这些解在四维空间中都画出来,想象一下,Xp是一个非原点的点,Xn是一个穿过原点的平面,那么Xp+Xn是两者的组合,是一个不经过原点的经过Xp的二维平面,注意它不是子空间。

秩r与Ax=b的解的关系
秩r的m×n的矩阵A始终有:r<=m, r<=n
1)列满秩r=n,各列线性无关
每一列都有主元,0个自由变量,此时零空间N(A)只有零向量,因为没有自由变量能够赋值,列的线性组合无法产生0列(回顾下第六课时和第三课时,其中3中讲到:如果存在非0向量X,使的AX=0,即X对A的列向量的线性组合可以得到0向量(有一列在线性组合中不起作用),那么A是不可逆的。)。Ax=b的全部解:0个或一个解,如果有解,即是唯一解特解Xp
以上例子 A通过消元变成行最简阶梯式,很容易看出来A的两列线性无关,所以R中两个主元。同时我也一眼能看出来有两个行是多余的,肯定R下面会有两个0行,因为行向量是2维的,因为前两行是线性无关的,2维平面中有两个向量线性无关,那该平面的所有向量都可以由这两个向量线性组合得到,所有会出现两个0行
上式特解为(1,1),即A的两列的和的线性组合
2)行满秩r=m,各行线性无关
此时消元会得到每一行都有一个主元,自由变量n-r(n-m)个,此时对任意的b,Ax=b都有解。
3)r=m=n,满秩方阵,行,列线性无关
零空间只包含0向量,此时对于任意的b,Ax=b都有解。由r=n知道有唯一解。


总结:矩阵A为m×n的矩阵,Ax=b的解的情况
r=m=n    R=I   有唯一解       b是A列向量的线性组合         
r=n<m             有0解(无解)或唯一解     b如果恰好是A的列的线性组合则有唯一解
r=m<n             有无穷个解         特解+零空间
r<m,r<n           有0解或无穷解      如果b的行和A的行向量之间有相同的组合关系,那有无穷解,否则有0解

猜你喜欢

转载自blog.csdn.net/neptune4751/article/details/78921713