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.mzph.cn/news/588960.shtml

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

相关文章

老子的《道德经》透露,不努力反而更成功

人类生而自由&#xff0c;但到处都是枷锁。 永远不要怀疑经过慎思且足够投入的一小群人能否改变这个世界。事实上&#xff0c;只有他们才办得到。 优美灵魂的两个发展方向&#xff1a;崇拜道德的天才&#xff0c;对别人实行道德的判断。 一、道 《道德经》开始的名字是《老子…

关键字:try-catch关键字

在 Java 中&#xff0c;try-catch关键字用于异常处理。它们允许编写代码来捕获和处理异常&#xff0c;以确保程序能够在出现问题时合理地处理它们而不会崩溃。 以下是try-catch关键字的基本语法&#xff1a; 在try块中编写可能会抛出异常的代码。如果在try块中的任何代码抛出…

全连接层:神经网络的桥梁

全连接层&#xff1a;神经网络的桥梁 大家好&#xff0c;我是免费搭建查券返利机器人赚佣金就用微赚淘客系统3.0的小编&#xff0c;也是冬天不穿秋裤&#xff0c;天冷也要风度的程序猿&#xff01;今天&#xff0c;让我们一同探讨深度学习中至关重要的一环——“全连接层”&am…

Linux上使用Certbot生成免费SSL证书

Linux上使用Certbot生成免费SSL证书 补充&#xff1a;certbot在大多数linux上是自带的&#xff0c;若没有可执行以下命令下载 Ubuntu&#xff1a; sudo apt-get update sudo apt-get install --only-upgrade certbotCentos sudo yum update sudo yum install certbot一、暂时停…

JavaScript setTimeout和setInterval的用法与区别详解

目录 I. 总述 II. setTimeout()函数 III. setInterval()函数 IV. 新年倒计时案例 Javascript的setTimeOut和setInterval函数应用非常广泛&#xff0c;它们都用来处理延时和定时任务&#xff0c;下面这篇文章主要给大家介绍了关于JavaScript setTimeout和setInterval的用法与…

读取某股票的日线数据

只需修改对应股票的ts_code&#xff0c;start_date&#xff0c;end_date即可获取对应股票的全部数据。 import tushare as ts import pandas as pdpd.set_option(expand_frame_repr, False) # 当列太多时不换行 pd.set_option(display.max_rows, 5000) # 最多显示数据的行数…

结构型设计模式—装饰器模式

装饰器模式&#xff1a;不改变原对象的基础上&#xff0c;灵活地给对象添加额外职责。装饰器相比于生成子类更加灵活。即将实现同一接口的父类当做参数传入包装器对象&#xff0c;动态创建出新的对象。 给对象添加新行为最简单直观的办法就是扩展本体对象&#xff0c;通过继承…

编程笔记 html5cssjs 018 HTML颜色

编程笔记 html5&css&js 018 HTML颜色 一、HTML 颜色二、HTML中设置颜色值注意 颜色是视觉中重要因素&#xff0c;尤其是处理人机界面中&#xff0c;更是要处理颜色设置和搭配。在网页中&#xff0c;提供了设置颜色的一些方案&#xff0c;需要我们认真学习和掌握。 一、…

HTML5和JS实现新年礼花效果

HTML5和JS实现新年礼花效果 2023兔年再见&#xff0c;2024龙年来临了&#xff01; 祝愿读者朋友们在2024年里&#xff0c;身体健康&#xff0c;心灵愉悦&#xff0c;梦想成真。 下面是用HTML5和JS实现新年礼花效果&#xff1a; 源码如下&#xff1a; <!DOCTYPE html>…

MySQL数据库学习一

1 什么是数据库的事务&#xff1f; 1.1 事务的典型场景 在项目里面&#xff0c;什么地方会开启事务&#xff0c;或者配置了事务&#xff1f;无论是在方法上加注解&#xff0c;还 是配置切面。 <tx:advice id"txAdvice" transaction-manager"transactionMa…

个人简历范本(精选5篇)

HR浏览一份简历也就25秒左右&#xff0c;如果你连「好简历」都没有&#xff0c;怎么能找到好工作呢&#xff1f; 如果你不懂得如何在简历上展示自己&#xff0c;或者觉得怎么改简历都不出彩&#xff0c;那请你一定仔细读完。 个人求职简历第 1 篇 男 22 本科 AI简历 市场营…

007、控制流

先看下本篇学习内容&#xff1a; 通过条件来执行 或 重复执行某些代码 是大部分编程语言的基础组成部分。在Rust中用来控制程序执行流的结构主要就是 if表达式 与 循环表达式。 1. if表达式 if表达式允许我们根据条件执行不同的代码分支。我们提供一个条件&#xff0c;并且做出…

【NTN 卫星通信】Oneweb星座以及Oneweb与Starlink比较

1 什么是OneWeb OneWeb于2012年以WorldVu的名义成立&#xff0c;于2020年开始构建其星座。然而&#xff0c;对于这家英国公司来说&#xff0c;这是一个艰难的旅程&#xff0c;OneWeb于2020年3月宣布破产&#xff0c;并认为covid-19大流行是一个主要因素。OneWeb星座当时仅完成…

一次降低进程IO延迟的性能优化实践——基于block层bfq调度器

如果有个进程正频繁的读写文件&#xff0c;此时你vim查看一个新文件&#xff0c;将会出现明显卡顿。即便你vim查看的文件只有几十M&#xff0c;也可能会出现卡顿。相对的&#xff0c;线上经常遇到IO敏感进程偶发IO超时问题。这些进程一次读写的文件数据量很少&#xff0c;正常几…

Redis的缓存过期淘汰策略

&#x1f44f;作者简介&#xff1a;大家好&#xff0c;我是爱吃芝士的土豆倪&#xff0c;24届校招生Java选手&#xff0c;很高兴认识大家&#x1f4d5;系列专栏&#xff1a;Spring源码、JUC源码、Kafka原理、分布式技术原理、数据库技术&#x1f525;如果感觉博主的文章还不错的…

RocketMQ(Linux版本5.1.4)

1、停止之前的运行服务 [roottssvr1-c1 rocketmq-all-4.7.0-bin-release]# sh bin/mqshutdown namesrv No mqnamesrv running. [roottssvr1-c1 rocketmq-all-4.7.0-bin-release]# [roottssvr1-c1 rocketmq-all-4.7.0-bin-release]# [roottssvr1-c1 rocketmq-all-4.7.0-bin-r…

Debezium发布历史39

原文地址&#xff1a; https://debezium.io/blog/2018/09/19/debezium-0-8-3-final-released/ 欢迎关注留言&#xff0c;我是收集整理小能手&#xff0c;工具翻译&#xff0c;仅供参考&#xff0c;笔芯笔芯. Debezium 0.8.3.Final 发布 2018 年 9 月 19 日 作者&#xff1a; …

LOJ #6278 数列分块2题解 2024年第一篇题解

Part #0 . 前言 \text{Part \#0 . 前言} Part #0 . 前言 数列分块1 分块是一种优雅的暴力。 Part #1 . 数列分块入门2 \text{Part \#1 . 数列分块入门2} Part #1 . 数列分块入门2 传送门 观察题目&#xff0c;我们可以发现题目是一个区间查询&#xff0c;区间修改。 首先&a…

数据结构式新年贺词

冒泡排序&#xff0c;选择排序&#xff0c;插入排序&#xff0c;快速排序&#xff0c;堆排序&#xff0c;归并排序&#xff0c;希尔排序&#xff0c;桶排序&#xff0c;基数排序新年帮您排忧解难。 有向图&#xff0c;无向图&#xff0c;有环图&#xff0c;无环图&#xff0c;…

汽车保养软件app开发步骤

“增强您的动力&#xff0c;为您的旅程加油——每一刻都讲述着关爱的故事。构建汽车维护软件app&#xff0c;为您的车辆提供数字化的维修站&#xff0c;从而开启长寿之路。智能驾驶、互联驾驶、自信驾驶。” 疯狂地搜索旧收据并猜测上次换油时间的日子已经一去不复返了。如果您…