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

news2024/11/20 10:30:56

一, 消元法 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) 

乘置换矩阵则为列变换

三阶矩阵

构造矩阵是通过对矩阵进行“行交换”来实现的

五, 逆矩阵 Inverse

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

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.coloradmin.cn/o/1039164.html

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈,一经查实,立即删除!

相关文章

【Xilinx】基于MPSoC的OpenAMP实现(一)

【Xilinx】基于MPSoC的OpenAMP实现(一) 一、开发环境1、开发思路2、下载官方bsp包 二、编译Linux1、配置petalinux环境变量2、创建工程3、进入目录4、设置缓存目录(重点:可离线编译,加快编译速度)5、配置u-…

浏览器截图扩展增加快捷键

Tabshot – 下载 🦊 Firefox 扩展(zh-CN) 最近一个用户找到我,想要这个浏览器扩展有一个快捷键截图功能。 我找了一下,发现火狐扩展的确支持快捷键 研究源码 about:support 配置文件夹,打开文件夹。 附…

MATLAB APP纯小白入门 两数相加

万事开头难,最怕第一次。使用matlab APP 实现两数求和,如下图所示,c a b,输入数字后,按 “” 就计算。 步骤 拖拽三个 Edit Field(Numeric) 过来,并且双击名字分别改为 a,b,c。注意修改名字后右边会有点变…

第6讲:v-for使用

目录 1.循环遍历 2.v-for遍历整形变量(99乘法表) 3.v-for遍历普通数组 4.v-for遍历数组对象 1.循环遍历 v-for指令基于一个数组渲染一个列表,它和JavaScript的遍历语法相似: v-for”item in list” list 是一个数组, i…

文件审计及文件完整性监控

什么是文件审核 对文件服务器中发生的所有事件的检查称为文件审核。这包括监视文件访问,其中包含谁访问了什么文件、何时以及从何处访问的详细信息;对访问最多和修改的文件的分析;成功和失败的文件访问尝试;等等。文件服务器审核过程的主要目标是跟踪在配置的服务器…

复习Day03:数组part03:76 . 最小覆盖子串、438. 找到z字符串z中所有字母异位词

之前的blog链接:https://blog.csdn.net/weixin_43303286/article/details/131700482?spm1001.2014.3001.5501 我用的方法是在leetcode再过一遍例题,明显会的就复制粘贴,之前没写出来就重写,然后从拓展题目中找题目来写。辅以Lab…

PY32F003F18之窗口看门狗

一、PY32F003F18窗口看门狗特点: 即使窗口看门狗被禁止,窗口看门狗的"递减计数器"也会继续递减计数。 二、窗口看门狗复位的条件: 1、将"控制寄存器WWDG_CR"中的WDGA1,激活"窗口看门狗计数器等于0x3F"时,则产…

次时代摸鱼骚操作:人在办公室轻松观看家里电脑上的4k电影(移动端公网访问本地群辉存储视频文件)

如何使用iPhone15在办公室观看家里电脑上的4k电影? 文章目录 如何使用iPhone15在办公室观看家里电脑上的4k电影?1.使用环境要求:2.下载群晖videostation:3.公网访问本地群晖videostation中的电影:4.公网条件下使用电脑…

Unity3D 简易音频管理器

依赖于Addressable 依赖于单例模板&#xff1a;传送门 using System.Collections.Generic; using System.Security.Cryptography; using System; using UnityEngine; using UnityEngine.AddressableAssets;namespace EasyAVG {public class AudioManager : MonoSingleton<…

2023第十二届中国智能产业高峰论坛之文档大模型的探索与思考

文章目录 前言合合信息多模态大模型与文档图像智能理解文档图像分析识别与理解的技术难题文档图像分析与预处理文档解析与识别版面分析与还原文档信息抽取与理解AI安全知识化&存储检索和管理 文档图像的分析识别与理解和大模型的关系文档图像大模型的进展LayoutLMUDOPDonut…

科技云报道:勒索金额再创新高,企业应如何防范?

科技云报道原创。 今年上半年的数据表明&#xff0c;勒索软件活动和赎金金额有望创下历史新高。 Check Point Research在《2023 年年中安全报告》中指出&#xff0c;今年上半年&#xff0c;随着新的勒索软件团伙不断涌现&#xff0c;勒索软件攻击态势持续升级。 区块链分析公…

【Proteus仿真】【STM32单片机】基于单片机的智能晾衣架控制系统

文章目录 一、功能简介二、软件设计三、实验现象联系作者 一、功能简介 系统运行后&#xff0c;LCD1604显示传感器检测的温湿度、光线强度和风速&#xff0c;工作模式&#xff0c;以及相应阈值&#xff0c;系统工作状态等&#xff1b;系统默认为自动模式&#xff0c; 可通过K4…

企业如何寻找适合的CRM系统软件

在数字化转型步伐不断加快的今天&#xff0c;CRM软件成为企业必不可少的一项工具。如果您并不清楚如何寻找合适的CRM系统软件&#xff1f;不妨通过下列几点来寻觅。 1.根据CRM系统行业和类别 例如企业重视营销功能&#xff0c;搜寻营销型CRM软件&#xff1b; 企业希望梳理好…

ideogram.ai 不同风格的效果图

https://ideogram.ai/ 提示词&#xff1a; French bulldog with sunglasses, playing skateboarding, speed up, happiness, front viewPhoto 相片 正常照片 Poster 海报 偏绘画&#xff0c;清晰的勾线 3D Render 3D 渲染 胶质感&#xff0c;像 3D 模型 Typography …

王珊教授:坚定信念走自主可控之路

在人大金仓2023新品发布会上&#xff0c;我国数据库泰斗、CCF最高科学技术奖获得者、原人民大学信息学院院长、原中国计算机学会&#xff08;CCF&#xff09;副理事长、人大金仓首任董事长王珊教授发表了致辞。 去年&#xff0c;王珊教授同样见证了人大金仓KESV9产品的发布。一…

氨基酸代谢:从基础到应用,揭示其在健康与疾病的角色

氨基酸(Amino Acid, 简称AA ) 是蛋白质合成的基石。它们是正常细胞生长、分化和功能所必需的细胞结构元素和能量来源。 蛋白质是长链氨基酸。人的身体有数千种不同的蛋白质&#xff0c;每种蛋白质都有重要的作用。每种蛋白质都有自己的氨基酸序列。该序列使蛋白质呈现不同的形状…

freemarker自定义模板

模板编程器指南 <dependency><groupId>org.freemarker</groupId><artifactId>freemarker</artifactId><version>2.3.31</version> </dependency>freemarker官网参考&#xff1a; https://freemarker.apache.org/docs/pgui_qu…

C# OpenCvSharp Yolov8 Cls 图像分类

效果 项目 代码 using OpenCvSharp; using OpenCvSharp.Dnn; using System; using System.Collections.Generic; using System.ComponentModel; using System.Data; using System.Drawing; using System.Linq; using System.Text; using System.Windows.Forms;namespace OpenC…

【OpenPLC学习】RK3568上运行OpenPLC

1 下载Runtime源码 git clone https://github.com/thiagoralves/OpenPLC_v3.git2 安装程序 ./install.sh linux3 在RK3568上运行Runtime sh start_openplc.sh4 在windows网页端登录 账号&#xff1a;openplc 密码&#xff1a;openplc 6 下载OpenPLC Editor https://git…

​旅行季《乡村振兴战略下传统村落文化旅游设计》许少辉八一著作想象和世界一样宽广

​旅行季《乡村振兴战略下传统村落文化旅游设计》许少辉八一著作想象和世界一样宽广