麻省理工学院 - MIT - 线性代数学习笔记

学习视频地址

文章目录

  • 1.01方程组的几何解释
  • 2.02矩阵消元
  • 3.03乘法和逆矩阵
    • 乘法
  • 4.04矩阵A的LU分解
  • 5.05转置,置换,向量空间
    • 置换
    • 转置
    • 向量空间
  • 6.06列空间和零空间
  • 7.07求解Ax=0:主变量,特解


1.01方程组的几何解释

对于二元方程组:row picture

在这里插入图片描述

方程转化为矩阵 * 向量 = 向量的形式,并且二维三维可以用图像几何表示出来(解析几何)

用column picture进行表示(线性组合):
x [ 2 − 1 ] + y [ − 1 2 ] = [ 0 3 ] x\begin{bmatrix}2 \\ -1 \\ \end{bmatrix}+y\begin{bmatrix}-1 \\ 2 \\ \end{bmatrix}= \begin{bmatrix}0 \\ 3 \\ \end{bmatrix} \\ x[21]+y[12]=[03]

在这里插入图片描述

x=1,y=2,上述是用向量的角度进行求解。

任意两个不共线的非零向量能够表示平面中的任意向量

对于三元方程组:row picture

在这里插入图片描述

column picture(向量表示):

在这里插入图片描述

x=0,y=0,z=1
A x = b Ax=b Ax=b
是否对于任意b都有上式可解 = 列的线性组合是否能覆盖整个三维(n维)空间

非奇异矩阵/可逆矩阵可达到上面的要求(可逆矩阵被称为非奇异矩阵,因为可逆矩阵的行列式不为0,在数学中如果一个矩阵的行列式为0,则称这个矩阵为奇异矩阵。)

是3个向量都处于同一平面,即维度减一,所以涉及第三维度的向量b无法构造。

对于n个向量也是如果最后消解只有n-1维那也无解(一种平面只能覆盖n维空间的某n-1维平面,最后的求解也只能在这n-1维平面展开)。

方程组的矩阵形式:
A x = b ( A x 看做 A 各列的线性组合 ) Ax=b(Ax看做A各列的线性组合) Ax=b(Ax看做A各列的线性组合)
在这里插入图片描述

2.02矩阵消元

x + 2 y + z = 2 3 x + 8 y + z = 12 4 y + z = 2 A = [ 1 2 1 3 8 1 0 4 1 ] x+2y+z=2\\ 3x+8y+z=12\\ 4y+z=2\\ A=\begin{bmatrix}1 & 2 & 1 \\3 & 8 & 1 \\0 & 4 &1 \end{bmatrix} x+2y+z=23x+8y+z=124y+z=2A= 130284111

消元

在这里插入图片描述

增广矩阵:Ab

回代(也就是加入b进行求解)

在这里插入图片描述

方程从上往下消元,从下往上回代求解

消元矩阵:

xx

在这里插入图片描述

由上式可知,我们把左边的矩阵(初等矩阵)当成3个1*3的行向量,右边的结果的第一行=左边第一行的第一个数 * 中间第一行 + 左边第一行的第二个数 * 中间第二行 + 左边第一行的第三个数 *中间第三行,所以结果为:1,0,0,以此类推第二行为-3,1,0(-3 * 中间第一行+1 * 中间第二行),第三行为0,0,1

在这里插入图片描述

E 21 表示初等矩阵第二行第一个位置发生了变换:变为 0 E_{21}表示初等矩阵第二行第一个位置发生了变换:变为0 E21表示初等矩阵第二行第一个位置发生了变换:变为0
在这里插入图片描述

E 32 ( E 21 A ) = U 相当于 ( E 32 E 21 ) A = U 结合律 E_{32}(E_{21}A)=U \\ 相当于(E_{32}E_{21})A=U\\结合律 E32(E21A)=U相当于(E32E21)A=U结合律
置换矩阵P

20231210160047641.png&pos_id=img-Q3Rr62RZ-1702196289492)

行变换左乘,列变换右乘

在这里插入图片描述

矩阵中基本没有交换律

逆矩阵:从U变回A

我们要取消这次消元,我们就希望E’EA=A,所以让E’E等于单位矩阵即可,E’是逆矩阵

变换后的矩阵EA,如何变回A?以行变换为例(左侧乘),就是找一个矩阵X使得XEA为A,此时若XE为单位阵,则XEA为A

在这里插入图片描述

E − 1 E = I E^{-1}E=I E1E=I

3.03乘法和逆矩阵

乘法

方法一:行与列之间的点积

在这里插入图片描述

C 34 = ( r o w 3 o f A ) ⋅ ( c o l u m n 4 o f B ) C_{34}=(row3ofA) ·(column4ofB) C34=(row3ofA)(column4ofB)
在这里插入图片描述

A是m * n,B是n * P,C=AB大小为m * P

方法二:把B看做是P个单独的列向量,像上一节课一样对每个列向量单独乘积最后再合并:

在这里插入图片描述

方法三:把A看做m个单独的行向量,向上节课一样对每个行向量单独乘积最后再合并:

在这里插入图片描述

方法四:A中的列 * B中的行

在这里插入图片描述

在这里插入图片描述

矩阵分块:把每个块都看成单独的元素运用方法一就得到下图

在这里插入图片描述

逆矩阵(非奇异矩阵) * 原矩阵 = I(单位矩阵)

在这里插入图片描述

奇异矩阵(非逆)可以找到一个向量x 使得Ax=0,也就是A中的向量组共面

Ax=0 有非零解,说明A中有的列没有贡献,所以不可逆

在这里插入图片描述

反证法:逆乘原得到单位矩阵,单位矩阵乘任一矩阵得原矩阵,所以得到X=0,但是X≠0,所以不成立

如何求逆矩阵:对于下列矩阵本质上就求两个二元一次方程组

在这里插入图片描述

同时求解两个方程:增广矩阵形式

高斯-若尔当消元法:

在这里插入图片描述

原理可理解为,假定有B矩阵,BA=I,那么右乘A逆得BI=A逆,如此可以得到逆矩阵

相当于对AI同时乘以A逆,得到I和A逆
A I = = > A ∗ A ′ I A ′ = = > I A ′ AI==>A*A^{\prime}I{\,}A^{\prime}==>IA^{\prime} AI==>AAIA==>IA

在这里插入图片描述

4.04矩阵A的LU分解

两矩阵相乘,且它们的逆均已知,那么AB的逆是什么 B − 1 A − 1 B^{-1} A^{-1} B1A1

A转置乘以A逆的转置等于单位E,定理A转置乘以A转置的逆等于单位E,所以A逆转置等于A转置的逆

如果还是不明白,因为我们已经知道A * A逆=A逆 * A=E,所以下面那个式子还可以写成,A转置的逆*A的转置=E,逆矩阵又是唯一的,所以A逆的转置=A转置的逆

这个由“行”组成的系数矩阵的每一行应该和之前A-1代表的系数矩阵的“列”相对应,因此就是(A-1)T*AT

将原本的A*A-1=I中左侧的A中的所有列看作用于线性组合的列向量,A-1中的每一列代表每次线性组合的系数组成
A A − 1 = I ( A − 1 ) T A T = I ( A − 1 ) T = ( A T ) − 1 AA^{-1}=I\\ (A^{-1})^TA^T=I\\ (A^{-1})^T=(A^{T})^{-1} AA1=I(A1)TAT=I(A1)T=(AT)1

在这里插入图片描述

把-4倍的第一行和1倍的第二行相加,放到第二行

在这里插入图片描述

E 21 = [ 1 0 − 4 1 ] E_{21}=\begin{bmatrix}1 & 0 \\-4 & 1 \end{bmatrix} E21=[1401]
在这里插入图片描述

L = E 21 − 1 L=E_{21}^{-1} L=E211
外链图片转存失败,源站可能有防盗链机制,建议将图片保存下来直接上传

E 32 E 31 E 21 A = U ( n o r o w e x c h a n g e ) A = E 21 − 1 E 31 − 1 E 32 − 1 U E_{32}E_{31}E_{21}A=U(no{\,\,} row{\,}exchange)\\ A=E_{21}^{-1}E_{31}^{-1}E_{32}^{-1}U E32E31E21A=U(norowexchange)A=E211E311E321U
dd

[ 1 0 0 2 1 0 0 5 1 ] = L ( l e f t o f U ) E A = U A = L U \begin{bmatrix}1 &0 & 0\\2 & 1&0\\0&5&1 \end{bmatrix}=L(left{\,}of{\,}U)\\ EA=U A=LU 120015001 =L(leftofU)EA=UA=LU
原因:然而L是恢复为original矩阵的逆过程,也就是从下向上处理的,同样使用较上的行来处理较下的行,然而自下而上不会更改较上的行。

原因:E是从上向下处理的,使用较上的行来处理较下的行,然而较上的行不一定是original的,有可能是经过处理修改的。

  1. L中的矩阵相乘的顺序很好,从左往右从上往下

  2. 初等矩阵都是单位阵来的,对角线都是1

在这里插入图片描述

类似于时间复杂度如果n=100,对于n * n的矩阵A需要进行多少次操作(每个元素一次乘法+一次减法进行小消元)变成下三角矩阵?

首先把矩阵留下第一行第一列的主元,第n行第一列为0,需要大约100^2次运算

同理,把矩阵留下第二行第二列的主元,第n行第二列为0,需要大约99^2次运算…

总操作数是:大约是 1 / 3 ∗ n 3 1/3*n^3 1/3n3的总操作次数
a l l o p r a t i o n s = n 2 + ( n − 1 ) 2 + . . . 1 2 = n ( n + 1 )( 2 n + 1 ) / 6 all{\,}oprations=n^2+(n-1)^2+...1^2=n(n+1)(2n+1)/6 alloprations=n2+(n1)2+...12=nn+1)(2n+1/6

在这里插入图片描述

对于增广矩阵右侧的b,需要n^2次操作次数

n * n的矩阵有 n ! n! n!种置换矩阵
P − 1 = P T P^{-1}=P^{T} P1=PT

在这里插入图片描述

5.05转置,置换,向量空间

置换

permutations P:execute row exchanges

A=LU的分解其实就是L是A行变换的过程U是行变换的结果

如果要进行 行变换A=LU 就要变成PA=LU

置换矩阵是行重新排列了的单位矩阵

n * n的矩阵有n!种置换矩阵,就是各行重新排列后所有可能的数目

P都是可逆的,并且满足以下式子
P − 1 = P T P T P = I P^{-1}=P^{T}\\ P^{T}P=I P1=PTPTP=I
在这里插入图片描述

转置

R^TR 得出来的矩阵都是对称的

在这里插入图片描述

( R T R ) T = R T R T T = R T R (R^TR)^T=R^{T}R^{TT}=R^TR (RTR)T=RTRTT=RTR

向量空间

在这里插入图片描述

向量空间对数乘和加法两种运算(线性组合)是封闭的

子空间内必须包含0向量

二维空间的子空间是过原点的一条直线

在这里插入图片描述

在这里插入图片描述

这一讲的核心思想是,通过某些向量构成一个向量组成的空间如果这些空间属于R3,那么它们构成的空间也在R3内,关键是,对其进行线性组合后仍然在子空间内

子空间:过原点、点(1,2,4)、点(3,3,1)的平面

如何构造矩阵列空间:取出每一列,线性组合,所有的线性组合就能构成列空间

6.06列空间和零空间

子空间S和T的交集是一个子空间

交集必含有零向量,线与线的交集是原点,线与面的交集要么是原点要么是线本身,面与面的交集要么是一条过原点的线要么是面本身

S和T的交集取向量v和w,则v,w既属于S也属于T,而S,T中的v,w是满足“加和数乘”规则的

在这里插入图片描述

在这里插入图片描述

在这里插入图片描述

三个列向量的线性组合无法充满整个四维空间

那么理解:四维空间里的三个向量最多只能组成四维空间里的三维空间,而要让B有解,就必须充满四维空间才行,因为B是四维空间里任意一个向量

问题:什么b让方程组有解:Ax=b

只有b是各列的线性组合时,Ax=b才有解

在这里插入图片描述

对于A,去掉某一列得到同样的列空间:去掉列三,因为列一和列二的线性组合能得到列三,也就是说列三没有贡献,列一列二为主列。(同理也可以去掉列一)

A的零空间是指是的Ax=0成立的向量空间

在这里插入图片描述

在这里插入图片描述

C ( A ) i n R 4 , x = [ x 1 x 2 x 3 ] i n R 3 C(A) in R^4, x=\begin{bmatrix}x_1 \\ x_2\\x_3 \end{bmatrix}in R^3 C(A)inR4x= x1x2x3 inR3

在这里插入图片描述

Ax=0包括 c [ 1 1 − 1 ] c\begin{bmatrix}1 \\ 1 \\-1 \end{bmatrix} c 111 的一条直线

%5Cimage-20231215145831155.png&pos_id=img-dMEAMKkN-1702624096791)

如果b不为零,Ax=b中x的解不构成子空间(0不是解),它是不过原点的一条直线。这条直线是 [c+1, c, -c],过(1,0,0)

在这里插入图片描述

7.07求解Ax=0:主变量,特解

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

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

相关文章

Scratch自制:《袁坤》游戏攻略

大家好!我也是很久没有动静了,这次我又来了,并且还带来了一个用Scratch制作的游戏,大家还记得我很久以前用C制作的《袁坤》吗?(详见C自制游戏《袁坤》1.2版本发布!-CSDN博客)这次它又…

种类并查集

最近玩的太嗨了,都忘了自己还有三篇博客还在拖更,也是今天一更到底好吧,边更新边写题,让看官老爷有更多的样题去联系 引入—— 在学这个之前,我相信各位应该已经接触过了并查集了吧,嗯?什么&a…

Kubernetes 1.24 版弃用 Dockershim 后如何迁移到 containerd 和 CRI-O

在本系列的上一篇文章中,我们讨论了什么是 CRI 和 OCI,Docker、containerd、CRI-O 之间的区别以及它们的架构等。最近,我们得知 Docker 即将从 kubernetes 中弃用!(查看 kubernetes 官方的这篇文章)那么让我…

VSCODE 下 openocd Jlink 的配置笔记

title: VSCODE 下 openocd Jlink 的配置笔记 tags: STM32HalCubemax 文章目录 内容VSCODE 下 openocd Jlink 的配置笔记安装完成后修改jlink的配置文件然后修改你的下载器为jlink烧录你的项目绝对会出现下面的问题那么打开下载的第一个软件 (点到这个jlink右键&…

Kafka架构详解之分区Partition

目录 一、简介二、架构三、分区Partition1.分区概念2.Offsets(偏移量)和消息的顺序3.分区如何为Kafka提供扩展能力4.producer写入策略5.consumer消费机制 一、简介 Apache Kafka 是分布式发布 - 订阅消息系统,在 kafka 官网上对 kafka 的定义…

【11】微服务链路追踪SkyWalking

1、skywalking是什么 1.1 链路追踪介绍 对于一个大型的几十个、几百个微服务构成的微服务架构系统,通常会遇到下面一些问题,比如: 如何串联整个调用链路,快速定位问题?如何缕清各个微服务之间的依赖关系?…

【AI学习】LLaMA 系列模型的进化(二)

在前面LLaMA 系列模型的进化(一)中学习了LLama模型的总体进化发展,再来看看其中涉及的一些重要技术。 PreLayerNorm Layer Norm有Pre-LN和Post-LN两种。Layer Normalization(LN)在Transformer架构中的放置位置对模型…

基于PaddleOCR + NLP实现证件识别

基于PaddleOCR NLP实现证件识别 PaddleOCR识别paddleOCR安装安装 anconda虚拟环境(可参考yolov5的安装教程) paddleOCR识别PaddleNLP模型信息抽取paddle打包exe 进行ocr识别 什么是PaddleOCR? PaddleOCR 旨在打造一套丰富、领先、且实用的 OCR 工具库,助…

【HarmonyOS开发】Navigation使用

简介 Navigation是路由容器组件,包括单栏(Stack)、分栏(Split)和自适应(Auto)三种显示模式。适用于模块内和跨模块的路由切换。 在页面跳转时,应该使用页面路由router,在页面内的页面跳转时,建议使用Navigation达到更好的转场动效…

Spring Boot集成Spring Batch快速入门Demo

1.什么是Spring Batch? Spring Batch 是一个轻量级的开源框架,它提供了一种简单的方式来处理大量的数据。它基于Spring框架,提供了一套批处理框架,可以处理各种类型的批处理任务,如ETL、数据导入/导出、报表生成等。S…

Armv8/Armv9架构的学习大纲-学习方法-自学路线-付费学习路线

本文给大家列出了Arm架构的学习大纲、学习方法、自学路线、付费学习路线。有兴趣的可以关注,希望对您有帮助。 如果大家有需要的,欢迎关注我的CSDN课程:https://edu.csdn.net/lecturer/6964 ARM 64位架构介绍 ARM 64位架构介绍 ARM架构概况…

uniapp,vue3上传图片组件封装

首先创建一个 components 文件在里面进行组件的创建 下面是 vip组件的封装 也就是图片上传组件 只是我的命名是随便起的 <template><!--图片 --><view class"up-page"><!--图片--><view class"show-box" v-for"(item,ind…

蓝桥杯Python算法竞赛常用的函数库

博客主页&#xff1a;音符犹如代码系列专栏&#xff1a;Python关注博主&#xff0c;后期持续更新系列文章如果有错误感谢请大家批评指出&#xff0c;及时修改感谢大家点赞&#x1f44d;收藏⭐评论✍ ​ 目录 math collectcions heapq functool itertools 常用的库函数 m…

web服务器1

&#xff08; 1 &#xff09;仅提供用户浏览的单向静态网页 单纯是由服务器单向提供数据给客户端&#xff0c; Server 不需要与 client 端有互动&#xff0c;所以你可以到该网站上去浏 览&#xff0c;但是无法进行数据的上传。 &#xff08; 2 &#xff09;提供用户互动接口的…

深入理解Linux网络(三):TCP对象创建

深入理解Linux网络&#xff08;三&#xff09;&#xff1a;TCP对象创建 TCP对象创建inet_createsock_init_data TCP对象创建 常见的三句TCP编程&#xff1a; int main() {int sk socket(AF_INET, SOCK_STREAM, 0);connect(sk, ...)recv(sk, ...) }简单的两三⾏代码&#xff…

十年前的老电脑能装win10吗_十年前的老电脑用U盘安装win10教程

十年前的老电脑能装win10吗&#xff1f;十年前的老电脑只要满足win10最低要求的配置都可以安装win10。安装win10方法很多&#xff0c;有一键重装方法、U盘安装、硬盘安装等方式&#xff0c;但最靠谱的方式还是U盘安装。十年前的老电脑用U盘安装win10首先要将u盘制作成u盘启动盘…

react自定义校验报错问题修复 ProFormText

1、以下是tsx组件 自定义校验告警导致表单无法提交问题修复 修改如下&#xff1a;

体态识别开发记录

1.多标签输出分类算法 https://aiuai.cn/aifarm1372.html 2.体态识别算法 体态对比 体态判断 2.1 3D建模的方案 stable diffusion 插件地址&#xff1a; ControlNet https://github.com/Mikubill/sd-webui-controlnet 3D Openpose https://github.com/nonnonstop/sd…

VsCode SSH远程连接服务器【内网穿透公网连接】

文章目录 1.前言2.VS code的安装和设置2.1 VS code的下载安装2.2 OpenSSH的启用2.3 为VS code配置ssh2.4 局域网内测试VS code的ssh连接2.5 Cpolar下载安装 3.Cpolar端口设置3.1 Cpolar云端设置3.2 Cpolar本地设置 4.公网访问测试5.结语 1.前言 记得笔者小时候看电视&#xff…

仿源码大师主界面UI的iAPP源文件

仿源码大师首页主界面的布局 首页&#xff0c;分类&#xff0c;需求&#xff0c;我的 就只有这几个界面内容而已 资源静态 没有任何动画和功能 纯UI布局 纯UI布局 他的最新版已经不是这个UI布局 放心使用 以学习参考为目的&#xff0c;如有不妥望告知 原创&#xff0c;纯…