MIT线性代数笔记-第34讲-左右逆,伪逆

目录

  • 34.左右逆,伪逆
      • 左右逆
      • 伪逆
    • 打赏

34.左右逆,伪逆

左右逆

之前讲到的逆都是针对可逆方阵而言的,对于长方矩阵,实际上也有广义的逆,那就是左逆和右逆

  1. 左逆

    当矩阵列满秩,即 r = n r = n r=n时,该矩阵有左逆(虽然各列线性无关,但是 r < m r < m r<m,列向量无法组成一组基,所以没有右逆),设该矩阵为 A A A A A A列满秩,所以 A T A A^T A ATA是可逆矩阵,有 ( A T A ) − 1 A T A = I (A^T A)^{-1} A^T A = I (ATA)1ATA=I,所以左逆 A l e f t − 1 = ( A T A ) − 1 A T A^{-1}_{left} = (A^T A)^{-1} A^T Aleft1=(ATA)1AT,且得到的单位矩阵是 n n n阶的

  2. 右逆

    当矩阵 A A A行满秩时, A A T A A^T AAT是可逆矩阵,有 A A T ( A A T ) − 1 = I A A^T (A A^T)^{-1} = I AAT(AAT)1=I,所以有右逆 A r i g h t − 1 = A T ( A A T ) − 1 A^{-1}_{right} = A^T (A A^T)^{-1} Aright1=AT(AAT)1,且得到的单位矩阵是 m m m阶的

  3. 当列满秩时,有 A l e f t − 1 A = I A^{-1}_{left} A = I Aleft1A=I,但是如果把左逆放到右边,就可以得到: A A l e f t − 1 = A ( A T A ) − 1 A T A A^{-1}_{left} = A (A^T A)^{-1} A^T AAleft1=A(ATA)1AT,即 A A A列空间的投影矩阵;同理,当行满秩时, A r i g h t − 1 A A^{-1}_{right} A Aright1A A A A行空间的投影矩阵


伪逆

  1. 伪逆

    可以发现,矩阵在左或右无逆是分别由零空间和左零空间不只有 0 ⃗ \vec{0} 0 引起的,因为如果左乘了零空间的非零向量得到 0 ⃗ \vec{0} 0 ,乘上任何矩阵都无法得到原来的向量,右乘同理

    但是如果不考虑这两个零空间呢,比如只让矩阵乘上它行空间中的向量,很明显这会得到它列空间中的向量,接下来思考这个过程是否可逆

    因为行空间和列空间的维数一致,所以二者中的向量可能存在一定的对应关系,也就是映射,考虑行空间中的两个不同向量 x ⃗ , y ⃗ \vec{x} , \vec{y} x ,y ,让它们分别对应列空间中的 A x ⃗ , A y ⃗ A \vec{x} , A \vec{y} Ax ,Ay

    证明 A x ⃗ ≠ A y ⃗ A \vec{x} \ne A \vec{y} Ax =Ay

    ​    若 A x ⃗ = A y ⃗ A \vec{x} = A \vec{y} Ax =Ay ,即 A ( x ⃗ − y ⃗ ) = 0 ⃗ A (\vec{x} - \vec{y}) = \vec{0} A(x y )=0 ,则 x ⃗ − y ⃗ \vec{x} - \vec{y} x y 属于 A A A的零空间,又 x ⃗ , y ⃗ \vec{x} , \vec{y} x ,y 都属于 A A A的行空间,所以着 x ⃗ − y ⃗ \vec{x} - \vec{y} x y 同时属于零空间和行空间,而行空间与零空间互为正交补,所以 x ⃗ − y ⃗ \vec{x} - \vec{y} x y 只能为 0 ⃗ \vec{0} 0 ,即 x ⃗ = y ⃗ \vec{x} = \vec{y} x =y ,与条件矛盾,假设不成立,因而 A x ⃗ ≠ A y ⃗ A \vec{x} \ne A \vec{y} Ax =Ay

    由此可以说明行空间和列空间中向量的一一对应关系,而由列空间中的向量得到行空间中对应向量所用的矩阵即为原矩阵的伪逆,记作 A + A^+ A+,即 x ⃗ = A + ( A x ⃗ ) \vec{x} = A^+ (A \vec{x}) x =A+(Ax )

    考虑 A + A A^+ A A+A对任意向量的作用,对于任意一个 n n n维向量 x ⃗ \vec{x} x ,它一定由零空间和行空间中的向量线性组合而成,零空间中的向量乘 A A A得到 0 ⃗ \vec{0} 0 ,而行空间中的向量乘 A + A A^+ A A+A得到其本身,所以 A + A A^+ A A+A可视为一个投影矩阵,用于把任意 n n n维向量投影至 A A A的行空间;同理, A A + A A^+ AA+也可视为一个投影矩阵,用于把任意 m m m维向量投影至 A A A的列空间

  2. 求伪逆

    使用 S V D SVD SVD求伪逆,先将一个任意矩阵 A A A分解为 U Σ V T U \Sigma V^T UΣVT

    Σ = [ σ 1 ⋯ 0 0 ⋯ 0 ⋮ ⋱ ⋮ ⋮ ⋯ ⋮ 0 ⋯ σ r 0 ⋯ 0 0 ⋯ 0 0 ⋯ 0 ⋮ ⋯ ⋮ ⋮ ⋱ ⋮ 0 ⋯ 0 0 ⋯ 0 ] \Sigma = \begin{bmatrix} \sigma_1 & \cdots & 0 & 0 & \cdots & 0 \\ \vdots & \ddots & \vdots & \vdots & \cdots & \vdots \\ 0 & \cdots & \sigma_r & 0 & \cdots & 0 \\ 0 & \cdots & 0 & 0 & \cdots & 0 \\ \vdots & \cdots & \vdots & \vdots & \ddots & \vdots \\ 0 & \cdots & 0 & 0 & \cdots & 0 \end{bmatrix} Σ= σ10000σr0000000000 ,很容易求得 Σ \Sigma Σ的伪逆为 [ 1 σ 1 ⋯ 0 0 ⋯ 0 ⋮ ⋱ ⋮ ⋮ ⋯ ⋮ 0 ⋯ 1 σ r 0 ⋯ 0 0 ⋯ 0 0 ⋯ 0 ⋮ ⋯ ⋮ ⋮ ⋱ ⋮ 0 ⋯ 0 0 ⋯ 0 ] \begin{bmatrix} \dfrac{1}{\sigma_1} & \cdots & 0 & 0 & \cdots & 0 \\ \vdots & \ddots & \vdots & \vdots & \cdots & \vdots \\ 0 & \cdots & \dfrac{1}{\sigma_r} & 0 & \cdots & 0 \\ 0 & \cdots & 0 & 0 & \cdots & 0 \\ \vdots & \cdots & \vdots & \vdots & \ddots & \vdots \\ 0 & \cdots & 0 & 0 & \cdots & 0 \end{bmatrix} σ110000σr10000000000

    因为 U , V T U , V^T U,VT均为正交矩阵,是可逆的,且伪逆在一定条件下本质上对应一个逆操作,所以 A + = V Σ + U T A^+ = V \Sigma^+ U^T A+=VΣ+UT


打赏

制作不易,若有帮助,欢迎打赏!
赞赏码

支付宝付款码

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

如若内容造成侵权/违法违规/事实不符,请联系我们进行投诉反馈qq邮箱809451989@qq.com,一经查实,立即删除!

相关文章

x-cmd-pkg | deno - 更快更强的 JS 和 TS 运行时

目录 简介首次用户技术特点进一步阅读 简介 Deno 是一个基于 V8 引擎和 Rust 语言构建的 JavaScript 和 TypeScript 运行时环境&#xff0c;于 2018 年由 Ryan Dahl 在演讲中宣布&#xff0c;并在 2020 年正式发布 1.0&#xff0c;目标是提供一个高效且安全的脚本环境。 安全…

PostgreSQL16.1(Windows版本)

1、卸载原有的PostgreSQL &#xfeff; &#xfeff; 点击Next即可。 &#xfeff;&#xfeff; 点击OK即可。 卸载完成。 2、安装 &#xff08;1&#xff09; 前两部直接Next&#xff0c;第二部可以换成自己想要安装的路径。 &#xff08;2&#xff09; 直接点击Next。…

记矩阵基础概念

转自up&#xff1a;Naruto_Qcsdn&#xff1a;三维空间几何变换矩阵 先贴个站里分享的基础概念。 learn form 肥猫同学VFX b站&#xff1a;会用transform就会用矩阵 移动 旋转 缩放 1.transofrm ——输出变化矩阵 可以移动transform查看变化去理解 位移 缩放 旋转 由此—…

Flink学习-时间和窗口

在流数据处理应用中&#xff0c;一个很重要、也很常见的操作就是窗口计算。所谓的“窗口”&#xff0c;一 般就是划定的一段时间范围&#xff0c;也就是“时间窗”&#xff1b;对在这范围内的数据进行处理&#xff0c;就是所谓的 窗口计算。所以窗口和时间往往是分不开的。 时…

【js】js解析Token:

一、效果&#xff1a; 二、实现&#xff1a; export function getTokenObject(token) {//通过split()方法将token转为字符串数组,数组中的第二个字符进行解析return token ? JSON.parse(decodeURIComponent(escape(window.atob(token.split(".")[1].replace(/-/g &…

2023年亲身经历总结 | 记录从大学到现在历程

现在是2023年12月31日&#xff0c;回想2023年&#xff0c;是充满挑战的一年&#xff0c;也是成就斐然的一年。回首过去&#xff0c;仿佛一幅画卷展开&#xff0c;每一笔每一刻都镌刻着成长的印记。 以下文章从大学经历说起&#xff0c;到现在发展情况&#xff0c;希望我的经历对…

大数据应用领域:数据驱动一切

大数据出现的时间只有十几年&#xff0c;被人们广泛接受并应用只有几年的时间&#xff0c;但就是这短短几年的时间&#xff0c;大数据呈现出爆炸式增长的态势。在各个领域&#xff0c;大数据的身影几乎无处不在。今天我们通过一些大数据典型的应用场景分析&#xff0c;一起来看…

ALSA学习(5)——设备中的alsa

参考博客&#xff1a; https://blog.csdn.net/DroidPhone/article/details/7165482 &#xff08;一下内容基本是原博主的博客转载&#xff09; 文章目录 一、ASOC的由来二、硬件架构三、软件架构四、数据结构五、内核对ASoC的改进 一、ASOC的由来 ASoC–ALSA System on Chip …

Vue-Setup

一、setup概述 小小提示&#xff1a;vue3中可以写多个根标签。 Person.vue中内容 <template><div class"person"><h2>姓名&#xff1a;{{name}}</h2><h2>年龄&#xff1a;{{age}}</h2><!--定义了一个事件&#xff0c;点击这…

【网络面试(5)】收发数据及断开服务器(四次挥手)

前面了解到服务器和客户端在创建套接字&#xff0c;建立连接后&#xff0c;就可以进入到下一步&#xff0c;双发可以互相发送和接收数据&#xff0c;本篇博客就来学习一下这个过程。  我们印象里&#xff0c;发送数据应该是我们在浏览器输入网址&#xff0c;敲击回车的一瞬间&…

L1-077:大笨钟的心情

有网友问&#xff1a;未来还会有更多大笨钟题吗&#xff1f;笨钟回复说&#xff1a;看心情…… 本题就请你替大笨钟写一个程序&#xff0c;根据心情自动输出回答。 输入格式&#xff1a; 输入在一行中给出 24 个 [0, 100] 区间内的整数&#xff0c;依次代表大笨钟在一天 24 小时…

【LLM 】7个基本的NLP模型,为ML应用程序赋能

在上一篇文章中&#xff0c;我们已经解释了什么是NLP及其在现实世界中的应用。在这篇文章中&#xff0c;我们将继续介绍NLP应用程序中使用的一些主要深度学习模型。 BERT 来自变压器的双向编码器表示&#xff08;BERT&#xff09;由Jacob Devlin在2018年的论文《BERT:用于语言…

C/C++ 函数重载

函数多态是C在C语言的基础新增的功能。默认参数能够使用不同数目的参数调用同一个函数&#xff0c;而函数多态(函数重载)让您能够使用多个同名的函数。术语“多态”指的是有多种形式&#xff0c;因此函数多态允许函数可以有多种形式。类似地&#xff0c;术语“函数重载”指的是…

【时钟】分布式时钟HLC|Logical Time|Vector Clock|True Time

目录 简略 详细 附录 1 分布式系统不能使用NTP的原因 简略 分布式系统中不同于单机系统不能使用NTP(网络时间协议&#xff08;Network Time Protocol&#xff09;)来获取时间&#xff0c;所以我们需要一个特别的方式来获取分布式系统中的时间&#xff0c;mvcc也是使用time保证读…

2024最全面且有知识深度的web3开发工具、web3学习项目资源平台

在Web3技术迅速发展的时代&#xff0c;寻找一个综合且深入的Web3开发工具和学习项目资源平台变得至关重要。今天&#xff0c;我将向大家介绍一个非常有价值的网站&#xff0c;它就是https://web3x.world 。 Web3X是一个全面而深入的Web3开发者社区&#xff0c;为开发者们提供了…

最优化方法Python计算:无约束优化应用——神经网络回归模型

人类大脑有数百亿个相互连接的神经元&#xff08;如下图(a)所示&#xff09;&#xff0c;这些神经元通过树突从其他神经元接收信息&#xff0c;在细胞体内综合、并变换信息&#xff0c;通过轴突上的突触向其他神经元传递信息。我们在博文《最优化方法Python计算&#xff1a;无约…

跳跃表原理及实现

一、跳表数据结构 跳表是有序表的一种&#xff0c;其底层是通过链表实现的。链表的特点是插入删除效率高&#xff0c;但是查找节点效率很低&#xff0c;最坏的时间复杂度是O(N)&#xff0c;那么跳表就是解决这一痛点而生的。 为了提高查询效率&#xff0c;我们可以给链表加上索…

打破成本壁垒,免费SSL证书为中小企业保驾护航

HTTPS&#xff0c;这个曾经看似遥远的技术词汇&#xff0c;如今已与我们每个人的网络生活息息相关。而实现HTTPS加密传输的关键一环——SSL证书&#xff0c;正以其独特的安全性能&#xff0c;为网站筑起一道坚实的防护墙。更令人惊喜的是&#xff0c;免费SSL证书服务已经到来&a…

数据结构与算法教程,数据结构C语言版教程!(第二部分、线性表详解:数据结构线性表10分钟入门)三

第二部分、线性表详解&#xff1a;数据结构线性表10分钟入门 线性表&#xff0c;数据结构中最简单的一种存储结构&#xff0c;专门用于存储逻辑关系为"一对一"的数据。 线性表&#xff0c;基于数据在实际物理空间中的存储状态&#xff0c;又可细分为顺序表&#xff…

自动化网络故障修复管理

什么是故障管理 故障管理是网络管理的组成部分&#xff0c;涉及检测、隔离和解决问题。如果实施得当&#xff0c;网络故障管理可以使连接、应用程序和服务保持在最佳水平&#xff0c;提供容错能力并最大限度地减少停机时间。专门为此目的设计的平台或工具称为故障管理系统。 …