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块中的任何代码抛出…

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

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

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星座当时仅完成…

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…

数据结构式新年贺词

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

汽车保养软件app开发步骤

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

JavaScript 基础通关

快速熟悉 JavaScript 的基础语法&#xff0c;比较高级的比如事件放在后面的笔记中。 JavaScript 1. JavaScript 介绍 1.1 JavaScript 基本介绍 JavaScript 是一门运行在客户端&#xff08;浏览器&#xff09;的编程语言&#xff0c;实现人机交互的效果。实现网页特效、表单验…

Python入门学习篇(十四)——模块文件操作

1 模块 1.1 理解 包: python中带有__init__.py文件的文件夹 模块: 文件名(不包含.py后缀),如python官方的time.py中time就是模块1.2 示例代码 import datetime# 调用datetime模块中的datetime类的now()方法 t datetime.datetime.now() # 格式化输出日期和时间 strftime(&qu…

Vue2中使用echarts,并从后端获取数据同步

一、安装echarts npm install echarts -S 二、导入echarts 在script中导入&#xff0c;比如&#xff1a; import * as echarts from "echarts"; 三、查找要用的示例 比如柱状图 四、初始化并挂载 <template><div id"total-orders-chart" s…

开发Python网络爬虫应用,爬取链家新房楼盘信息保存到mongodb中,并分析相关数据

这里写自定义目录标题 爬取代码分析数据问题 爬取代码 import requests import time from lxml import html from pymongo import MongoClient import randomBASEURL https://cq.fang.lianjia.com/loupan/# 获取某市区域的所有链接 def get_areas(url):print(获取区县列表)# …

git解决冲突场景

文章目录 git解决冲突场景 git解决冲突场景 假设我们在公司开发了一个功能修改了一个文件 我们现在模拟修改文件之后提交一个版本到本地&#xff0c;但是不上传到远程仓库 假设我们现在回到家开发代码&#xff0c;需要拉去最新的代码 提示已经更新。根本没有最新的代码改动&am…

Codeforces Round 917 (Div. 2) A~F

A.Least Product(思维) 题意&#xff1a; 给出一个数组 a 1 , a 2 , . . . , a n a_1, a_2, ..., a_n a1​,a2​,...,an​&#xff0c;你可以进行若干次以下操作&#xff1a; 选择数组中的一个元素 a i a_i ai​&#xff0c;将这个数字修改为 0 ∼ a i 0 \sim a_i 0∼ai​之…

PHP与Angular详细对比 帮助你选择合适的项目技术

开发可有效扩展并提供诺克斯堡级安全性的Web应用程序和网站是每个开发人员的梦想。而使用这样的产品是每个用户的愿望。因此&#xff0c;为您的项目选择最合适和可靠的技术非常关键。 虽然PHP和Angular是完全不同的技术——PHP与JavaScript是一个更恰当的比较——但它们都广泛…

ACES 增强版不丹水稻作物地图(2016-2022 年)

ACES 增强版不丹水稻作物地图&#xff08;2016-2022 年&#xff09; 用于改善粮食安全决策的 2016-2022 年年度作物类型稻米地图仍然是不丹的一项挑战。这些地图是与不丹农业部和 SERVIR 合作开发的。通过专注于发展不丹的科学、技术、工程和数学 (STEM)&#xff0c;我们共同开…