MIT18.06线性代数课程笔记8:求解Ax=b、矩阵的秩以及矩阵的逆

2023-12-15 23:18

本文主要是介绍MIT18.06线性代数课程笔记8:求解Ax=b、矩阵的秩以及矩阵的逆,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

课程简介

18.06是Gilbert Strang教授在MIT开的线性代数公开课,课程视频以及相关资料请见https://ocw.mit.edu/courses/mathematics/18-06-linear-algebra-spring-2010/index.htm。

课程笔记

1. 求解Ax=b

这里涉及两个问题:1. 是否有解 2. 如何求解

1.1. 是否有解

有两种方法:1. bC(A)xRn,Ax=b , where Am×n 2. 对增广矩阵消元, b 所在列不是pivot column。

对于第一种方法,如MIT18.06线性代数课程笔记6:vector space,subspace,column space,null space所述,Ax A 的列向量的线性组合,所以AxC(A),而且 vC(A),x,Ax=c 。由此可证明第一个方法。但是这只是理论上的结论,实际应用中不如消元法方便。

如MIT18.06线性代数课程笔记7:使用消元法求解Null space中所述,对矩阵 A 消元,A=LU,可以得到pivot column和free column,且有任意free column都可以表示为所有pivot column的线性组合(更严格的结论是任意free column都可以表示为之前列的线性组合)。证明也很容易,因为if c=AxEc=EAx ,since E is invertible。所以对增广矩阵[A|b]做消元,即 E[A|b]=[U|b] ,若 ui=0 ,即 U 的第i行为0, bi0 ,则 Ax=b 无解,因为 b [U|b] 的pivot column。所以 bC(U)

1.2. 如何求解

如果有解的情况下如何求解的问题在MIT18.06线性代数课程笔记6:vector space,subspace,column space,null space有所阐述。具体地,就是先求解一个特解 xp xp 的获取比较随意,常用方法是对增广矩阵做消元之后令free column取0,这样就可以得到唯一解。然后因为 xRn,Ax=bA(xxp)=0xxpN(A) ,固所有解集为 {x:x=xp+xn,xnN(A)} 。至于 N(A) 的求解方法详见MIT18.06线性代数课程笔记7:使用消元法求解Null space,简单的说就是令free column取值为 nr 个线性无关的基向量,然后求解 Ax=0 ,得到 nr N(A) 的基向量。

2. 矩阵的秩

矩阵的秩现在的定义是消元后主元的个数,也是pivot column的个数,使用 r 表示。

Am×n,容易证明 rm , rn

并且通过 r n,m的关系,可以判断矩阵是否有解,以及解集的维数。

2.1. r==m 必定有解

笔者想到两种证明方法,一是 U 中不存在零行,所以任意b都可以;二是有 m 个pivot column,即有m个线性无关基向量,其线性组合遍布整个 Rm

2.2. r==n 有解必然唯一解

行满秩的情况下,没有free column,所以 N(A)={0⃗ } ,解集为 {xp}

但是不一定有解,因为 n m维向量构成的子空间维度为 n ,不遍布Rm。对于大多数的 b 都无解。

2.3. r==m==n有解且唯一解

行列均满秩的情况下,由上诉两个结论可知有解且唯一解。

2.4. 秩与可逆的关系

目前只讨论左逆元,即 A1A=I ,则有只有 r==n 时有逆元。即列满秩是 A 存在左逆元的充要条件。

  1. 列满秩 A 有左逆元

    如上所述,因为列满秩,所以E,EA=[I0]。即 E1,E2,[E1E2]A=[I0]E1,E1A=I E1 即为 A1 ,所以 A 存在左逆元。

  2. A有左逆元 列满秩

    因为 r(A1A)=nr(A)nr(A)=n

这篇关于MIT18.06线性代数课程笔记8:求解Ax=b、矩阵的秩以及矩阵的逆的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/498256

相关文章

C/C++中OpenCV 矩阵运算的实现

《C/C++中OpenCV矩阵运算的实现》本文主要介绍了C/C++中OpenCV矩阵运算的实现,包括基本算术运算(标量与矩阵)、矩阵乘法、转置、逆矩阵、行列式、迹、范数等操作,感兴趣的可以了解一下... 目录矩阵的创建与初始化创建矩阵访问矩阵元素基本的算术运算 ➕➖✖️➗矩阵与标量运算矩阵与矩阵运算 (逐元

利用Python快速搭建Markdown笔记发布系统

《利用Python快速搭建Markdown笔记发布系统》这篇文章主要为大家详细介绍了使用Python生态的成熟工具,在30分钟内搭建一个支持Markdown渲染、分类标签、全文搜索的私有化知识发布系统... 目录引言:为什么要自建知识博客一、技术选型:极简主义开发栈二、系统架构设计三、核心代码实现(分步解析

numpy求解线性代数相关问题

《numpy求解线性代数相关问题》本文主要介绍了numpy求解线性代数相关问题,文中通过示例代码介绍的非常详细,对大家的学习或者工作具有一定的参考学习价值,需要的朋友们下面随着小编来一起学习学习吧... 在numpy中有numpy.array类型和numpy.mat类型,前者是数组类型,后者是矩阵类型。数组

【学习笔记】 陈强-机器学习-Python-Ch15 人工神经网络(1)sklearn

系列文章目录 监督学习:参数方法 【学习笔记】 陈强-机器学习-Python-Ch4 线性回归 【学习笔记】 陈强-机器学习-Python-Ch5 逻辑回归 【课后题练习】 陈强-机器学习-Python-Ch5 逻辑回归(SAheart.csv) 【学习笔记】 陈强-机器学习-Python-Ch6 多项逻辑回归 【学习笔记 及 课后题练习】 陈强-机器学习-Python-Ch7 判别分析 【学

系统架构师考试学习笔记第三篇——架构设计高级知识(20)通信系统架构设计理论与实践

本章知识考点:         第20课时主要学习通信系统架构设计的理论和工作中的实践。根据新版考试大纲,本课时知识点会涉及案例分析题(25分),而在历年考试中,案例题对该部分内容的考查并不多,虽在综合知识选择题目中经常考查,但分值也不高。本课时内容侧重于对知识点的记忆和理解,按照以往的出题规律,通信系统架构设计基础知识点多来源于教材内的基础网络设备、网络架构和教材外最新时事热点技术。本课时知识

hdu 4565 推倒公式+矩阵快速幂

题意 求下式的值: Sn=⌈ (a+b√)n⌉%m S_n = \lceil\ (a + \sqrt{b}) ^ n \rceil\% m 其中: 0<a,m<215 0< a, m < 2^{15} 0<b,n<231 0 < b, n < 2^{31} (a−1)2<b<a2 (a-1)^2< b < a^2 解析 令: An=(a+b√)n A_n = (a +

线性代数|机器学习-P36在图中找聚类

文章目录 1. 常见图结构2. 谱聚类 感觉后面几节课的内容跨越太大,需要补充太多的知识点,教授讲得内容跨越较大,一般一节课的内容是书本上的一章节内容,所以看视频比较吃力,需要先预习课本内容后才能够很好的理解教授讲解的知识点。 1. 常见图结构 假设我们有如下图结构: Adjacency Matrix:行和列表示的是节点的位置,A[i,j]表示的第 i 个节点和第 j 个

论文阅读笔记: Segment Anything

文章目录 Segment Anything摘要引言任务模型数据引擎数据集负责任的人工智能 Segment Anything Model图像编码器提示编码器mask解码器解决歧义损失和训练 Segment Anything 论文地址: https://arxiv.org/abs/2304.02643 代码地址:https://github.com/facebookresear

数学建模笔记—— 非线性规划

数学建模笔记—— 非线性规划 非线性规划1. 模型原理1.1 非线性规划的标准型1.2 非线性规划求解的Matlab函数 2. 典型例题3. matlab代码求解3.1 例1 一个简单示例3.2 例2 选址问题1. 第一问 线性规划2. 第二问 非线性规划 非线性规划 非线性规划是一种求解目标函数或约束条件中有一个或几个非线性函数的最优化问题的方法。运筹学的一个重要分支。2

hdu 6198 dfs枚举找规律+矩阵乘法

number number number Time Limit: 2000/1000 MS (Java/Others)    Memory Limit: 32768/32768 K (Java/Others) Problem Description We define a sequence  F : ⋅   F0=0,F1=1 ; ⋅   Fn=Fn