02 MIT线性代数-矩阵消元 Elimination with matrices

一, 消元法 Method of Elimination

消元法是计算机软件求解线形方程组所用的最常见的方法。任何情况下,只要是矩阵A 可逆,均可以通过消元法求得Ax =b的解

eg:

我们将矩阵左上角的1称之为"主元一"(the first pivot),第一步要通过消元将第一列中除了主元之外的数字均变化为0。操作方法就是用之后的每一行减去第一行的适当倍数,此例中第二行应减去第一行的3倍。之后应对第三行做类似操作,本例中三行第一列数字已经为0,故不用进行操作\

处在第二行第二列的主元二为2,因此用第三行减去第二行的两倍进行消元,得到第三个主元为5。

矩阵A 为可逆矩阵,消元结束后得到上三角阵U (Uppertriangular matrix),其左侧下半部分的元素均为0,而主元1,2,5分列在U的对角线上。

Pivot cannot be 0 主元不能为0,如果恰好消元至某行,0出现在了主元的位置上,应当通过与下方一行进行"行交换"使得非零数字出现在主元位置上。如果0出现在了主元位置上,并且下方没有对等位置为非0数字的行,则消元终止,并证明矩阵A为不可逆矩阵,且线性方程组没有唯一解。

二, 回代 Back-Substitution

做方程的高斯消元时,需要对等式右侧的b 做同样的乘法和加减法。手工计算时比较有效率的方法是应用"增广矩阵"(augmented matrix),将b 插入矩阵A 之后形成最后一列,在消元过程中带着b 一起操作。(Matlab是算完系数矩阵再处理b的)

此时我们将原方程Ax =b 转化为了新的方程Ux =c

从最后一行得到z=-2,依次回代可以得到y=1 和 x=2

三, 消元矩阵 Elimination Matrices

矩阵运算的核心内容就是对"行"或者"列"进行独立操作。

如前一节课"列图像"部分所言,系数矩阵乘以未知数向量,相当于对系数矩阵的列向量进行线性组合

matrix times a column is a column

++combination of the columns of the matrix++


a row times matrix is a row


列向量的线性组合非常容易就接受了,左乘行向量这个总觉得很别扭,偏偏行向量在下面介绍消元矩阵时比较重要。后来想想觉得"列"操作就像是把向量开进矩阵,而"行操作"这个就像把向量倒车进入矩阵(如图中箭头所示)

Step 1 substract 3***×***row1 from row2

通过左乘矩阵E21 来实现原矩阵A 的第二行减去第一行的3倍这一过程。E21 的第二行使矩阵A 的行向量进行前述的线性组合,而其它两行为了保持与原矩阵相同,采用同阶单位阵I (identity matrix 和×1无区别) 的行向量。左乘的这个矩阵为"初等矩阵"(Elementary Matrix),因此记做E 。因为所乘行向量的倍数-3出现在E 矩阵的第二行第一列,因此将之标注为21。完成操作后矩阵变为E21A

Step 2 substract 2***×***row2 from row3

矩阵消元的第二步是完成矩阵E21A 的第三行减去第二行的2倍,通过左乘矩阵E32来实现这一过程

3x3矩阵的消元本来应该分三步完成,最终得到E32 (E31 (E21A ))。本例中E31 =I ,所以结果变为E32 (E21 A )=U ,因为矩阵运算符合结合律associative law括号移动,也可写作(E32 E21 )A =U 。可以记作EA =U

四, 置换矩阵 Permutation

exchange rows 1 and 2 (row operation)

乘置换矩阵可以完成原矩阵的行变换

exchange columns 1 and 2 (column operation)

乘置换矩阵则为列变换

三阶矩阵

构造P 矩阵是通过对I矩阵进行"行交换"来实现的

五, 逆矩阵 Inverse

这里主要讨论消元矩阵的逆矩阵。消元矩阵之逆矩阵的实施效果就是抵消原矩阵的消元操作。消元矩阵实现了对原矩阵A的操作,使第二行行向量3,8,1减掉了第一行1,2,1的3倍变为0,2,-2,则逆向操作就应该是把现在的第二行行向量0,2,-2加上第一行1,2,1的3倍,从而变回原来的第二行3,8,1

相关推荐
邓码农0074 小时前
批量发布 + 数据分析横评|聚媒通、融媒宝、新榜小豆芽、矩阵通优缺点全面拆解
矩阵·数据挖掘·数据分析
棒球1号位1 天前
全球赋能:中国棒球联盟矩阵的国际化人才培养路径·棒球1号位
线性代数·矩阵
梅梅绵绵冰1 天前
算法题-矩阵
数据结构·算法·矩阵
郝学胜-神的一滴2 天前
中级OpenGL教程 022:探秘三维世界的血脉传承——物体父子关系与矩阵递归奥义
c++·线性代数·unity·矩阵·游戏引擎·unreal engine·opengl
luoganttcc4 天前
指数分布的完整推导
线性代数·机器学习·概率论
爱刷碗的苏泓舒4 天前
平方根信息滤波:矩阵推导及 GNSS 参数估计应用
线性代数·算法·矩阵·gnss·参数估计·测量平差·平方根信息滤波
智灵鸟科技4 天前
封闭网络怎么安全地接外部能力:三条通道穷举、威胁封堵矩阵与残余风险
网络·安全·矩阵
zandy10114 天前
企业级BI平台选型指南:评估框架与核心能力矩阵
大数据·人工智能·矩阵
AI科技星4 天前
光速螺旋时空曲率挠率拓扑统一场论——四大力全域闭环求导、精算验证与第一性原理完备证明
线性代数·算法·决策树·机器学习·常温超导·ai科技星
棒球1号位4 天前
中国棒球联盟矩阵建设白皮书:解码本土化创新与国际接轨路径·棒球1号位
大数据·人工智能·矩阵