第二课时、矩阵的消元(Elimination)回代(back substitution)法
本课时主要讲解了矩阵的消元回代法,来求解方程组。该方法时高斯提出的,也叫做高斯消元法,主要是求解方程组的解,复杂矩阵一般不采用此方法,下面对高斯消元法进行简单的了解和学习。
一.消元法(Elimination)的两种情况
(1)主元存在,可直接消元
例子如下:
上述式子转换成矩阵形式:
其中是系数矩阵,是未知数,是方程组右侧值。
对进行消元,过程如下:
首先找到主元(主元不能为零),从上边矩阵可以知道,第一个主元为1,在第一行,我们保持第一行不变,消去第二行的主元下边也就是3这个数消为0.
上述就是消元的过程,第一个主元是第一行的第一个数;第二个主元是第二行的第二个数。执行完消元以后的矩阵如上所示。通过消元之后,就直接可以求出方程组的解了,其实消元的目的就是求解线性方程组,通过把不同行的未知数不断的消去。最后求解出方程组的解。但是如果第一个主元是零的话怎么办呢?
(2) 主元(pivot)不存在
1. (例子1):方程组无解的情况,消元失败
上述左边是一个方程组,右边是通过常规方法得到的方程组的解,我们可以看到方程组是无解的。在几何图像中我们也可以看出,下图是方程组的几何图。
图上可以看出,从方程组的行来看,这是两条平行线,永远是没有交点的,也就是说他们永远没有解。而从右边的列向量来看,他们是两个相对方向的向量,永远无法组合得到向量b。关于方程组的几何意义,我已经在上一节方程组的几何解释中说明了(行图像和列图像)https://blog.csdn.net/LLLLLLLLLLLLLLYYYYYY/article/details/103513854。上述的情况,我们是无法进行消元处理的,即所谓的消元法失败。
2. (例子2):方程组有无穷解的情况,消元失败
上述左边是一个方程组,右边是通过常规方法得到的方程组的解,我们可以看到方程组是有无穷个解的。在几何图像中我们也可以看出,下图是方程组的几何图。
同上原理一样,上边的方程组在几何中的表示如上所示,这是无穷多解的情况,这种情况下也是无法消元的。消元失败。
3. (例子3):当矩阵的主元是零的时候,我们无法消元的情况下,可以通过行变换来进行消元
上边的公式是一个方程组,把他写成矩阵的形式如下所示:
上式第一个主元的位置式0,不能消元,那么我们可以把行进行互换,找出新的主元,这在消元法中叫做行变换,其实就是方程组调换位置。调换之后消元如下:
直接得到我们需要的消元举着,要记得把b也进行调换,直接可以求得消元后的解。
二、矩阵消元
在这一部分,我们就直接给出一个矩阵,来对他进行消元,也就是求出直角三角型状。矩阵如下所示:
消元步骤如下:
上述是通过一步一步的消元得到了最终得消元结果,也就是三角型得结果。那么接下来我们分解步骤。
1.第一次消元其实是在原始矩阵的左边乘了一个消元矩阵,为了消去第二行第一列的数。
2.第二次消元其实是在第一次消元后矩阵的左边乘了一个消元矩阵,为了消去第三行第二列的数。
3.把上面两个步骤合并
是消元之后得到的矩阵,是原始的矩阵,是消元矩阵,两者可以利用乘法的结合率结合到一起。矩阵的结合律证明比较麻烦,我们只要了解就好,这些就交给大数学家们证明吧,我们站在巨人的肩膀上,知道了矩阵也可以用乘法的结合律。所以我们把相乘得到了原始矩阵的消元矩阵,如下所示:
上面就是最终得到的原始矩阵A的消元矩阵。
三、(置换矩阵)补充知识
1.行变换(矩阵左边)
2.列变换(矩阵右边)
四、引申知识(逆矩阵)
上述都是通过消元法找到消元矩阵。都是通过。那么我们思考一下,怎么通过得到呢,按照一般乘除法的思路,就是而已,可这是矩阵,方法不在适用。大家思考一下,我们引出下节要讲的内容:逆矩阵(Inverses)。