【矩阵论】4. 矩阵运算——广义逆——减号逆

news2024/11/26 14:55:11

在这里插入图片描述

4.5 减号逆

A = A m × n A=A_{m\times n} A=Am×n X = X n × m X=X_{n\times m} X=Xn×m ,有 A X A = A AXA=A AXA=A ,则称 X = X n × m X=X_{n\times m} X=Xn×m 为A的减号逆(一号逆),记为 X = A − = A ( 1 ) X=A^{-}=A^{(1)} X=A=A(1)

全体 A − A^{-} A 的集合记为 A { 1 } = { X ∣ A X A = A } A^{\{1\}}=\{X\mid AXA=A\} A{1}={XAXA=A}

  • A − ∈ A { 1 } A^{-}\in A^{\{1\}} AA{1}

4.5.1 性质

自反性: A A − A = A AA^{-}A=A AAA=A

幂等性: ( A − A ) 2 = A − A (A^{-}A)^2=A^{-}A (AA)2=AA ,且 ( A A − ) 2 = A A − (AA^{-})^2=AA^{-} (AA)2=AA ,其中A是方阵
在这里插入图片描述

A − A^{-} A 不唯一,可以看做 A − 1 A^{-1} A1 的推广

  • A = ( 1 0 ) A=\left(\begin{matrix}1\\0\end{matrix}\right) A=(10) ,可取 X=(1 0) 或 Y=(1 1) 作为A的减号逆
  • A − A^{-} A 唯一的阵:方阵 A n × n A_{n\times n} An×n可逆,则必有唯一 A − 1 = A + = A − A^{-1}=A^{+}=A^- A1=A+=A

若A为列满秩(高阵),则 A − A = I A^-A=I AA=I;若A为行满秩(低阵),则 A A − = I AA^-=I AA=I

在这里插入图片描述

4.5.2 计算

a. 求解 A X A = A AXA=A AXA=A

求解减号逆 A − A{-} A 即求解 A X A = A AXA=A AXA=A 的全体通解

在这里插入图片描述在这里插入图片描述

由 矩 阵 方 程 A X B = D 的 特 解 X 0 = A + = ( 1 , 0 , 0 ) , 故 通 解 为 A − = X = X 0 + Y − A + A Y A A + = ( 1 , 0 , 0 ) + ( a , b , c ) − ( a , 0 , 0 ) = ( 1 , b , c ) \begin{aligned} &由矩阵方程AXB=D的特解X_0=A^+=\left(1,0,0\right),故通解为A^-=X=X_0+Y-A^+AYAA^+\\ &=\left(1,0,0\right)+\left(\begin{matrix}a,b,c\end{matrix}\right)-\left(\begin{matrix}a,0,0\end{matrix}\right)=\left(\begin{matrix}1,b,c\end{matrix}\right) \end{aligned} AXB=DX0=A+=(1,0,0),A=X=X0+YA+AYAA+=(100)+(abc)(a00)=(1bc)

也可见 A − A^- A 不唯一

  • 对于高阶阵 A − = A + + ( Y − A + A Y A A + ) A^-=A^++(Y-A^+AYAA^+) A=A++(YA+AYAA+) 的计算比较复杂

b. 标准对角形

A = ( I r 0 0 0 ) m × n A=\left(\begin{matrix}I_r&0\\0&0\end{matrix}\right)_{m\times n} A=(Ir000)m×n ,则全体 A − = ( I r B C D ) n × m A^-=\left(\begin{matrix}I_r&B\\C&D\end{matrix}\right)_{n\times m} A=(IrCBD)n×m ,BCD为任一小块

在这里插入图片描述


在这里插入图片描述

SP

P A Q = ( I r 0 0 0 ) m × n PAQ=\left(\begin{matrix}I_r&0\\0&0\end{matrix}\right)_{m\times n} PAQ=(Ir000)m×n ,全体 A − = Q ( I r 0 0 0 ) n × m P A^{-}=Q\left(\begin{matrix}I_r&0\\0&0\end{matrix}\right)_{n\times m}P A=Q(Ir000)n×mP ,BCD为任一小块

c. 初等行,列变换(一般方法)

A = A m × n A=A_{m\times n} A=Am×n ,令 ( A I m I n 0 ) → 列 变 换 行 变 换 ( ( I r 0 0 0 ) m × n P Q 0 ) \left(\begin{array}{c:c}A&I_m\\\hdashline I_n&0\end{array}\right)\xrightarrow[列变换]{行变换}\left(\begin{array}{c:c}\left(\begin{matrix}I_r&0\\0&0\end{matrix}\right)_{m\times n}&P\\\hdashline Q&0\end{array}\right) (AInIm0) (Ir000)m×nQP0 ,则有 A − = Q ( I r B C D ) n × m P A^-=Q\left(\begin{matrix}I_r&B\\C&D\end{matrix}\right)_{n\times m}P A=Q(IrCBD)n×mP

eg

在这里插入图片描述

4.5.3 矩阵方程求解

前置知识:正规方程求解

a. 特解

在这里插入图片描述

b. 解空间

N ( A ) N(A) N(A) A Y = 0 AY=0 AY=0 的通解为 Y = ( I n − A − A ) y Y=(I_n-A^-A)y Y=(InAA)y ∀ y ∈ C n \forall y\in C^n yCn

N ( A ) = { Y ∣ A Y = 0 } N(A)=\{Y\vert AY=0\} N(A)={YAY=0} X = ( I n − A − A ) y , y = ( y 1 ⋮ y n ) ∈ C n X=(I_n-A^-A)y,y=\left(\begin{matrix}y_1\\\vdots\\y_n\end{matrix}\right)\in C^n X=(InAA)y,y=y1ynCn

设y的值域为 R,则 N ( A ) = R ( I n − A − A ) N(A)=R(I_n-A^-A) N(A)=R(InAA) ,维数 d i m N ( A ) = n − r ( A − A ) dim N(A)=n-r(A^-A) dimN(A)=nr(AA)

c. 通解

在这里插入图片描述

在这里插入图片描述

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

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

相关文章

【自然语言处理概述】百度百科数据爬取

【自然语言处理概述】百度百科数据爬取 作者简介:在校大学生一枚,华为云享专家,阿里云专家博主,腾云先锋(TDP)成员,云曦智划项目总负责人,全国高等学校计算机教学与产业实践资源建设…

氮化镓(GaN)功率半导体之预测

前言 氮化镓(GaN)是一种非常坚硬且在机械方面非常稳定的宽带隙半导体材料。由于具有更高的击穿强度、更快的开关速度,更高的热导率和更低的导通电阻,氮化镓基功率器件明显比硅基器件更优越。 氮化镓晶体可以在各种衬底上生长&…

充分利用自动化测试的 10 个最佳实践

你试过吃带壳的坚果吗?如果是,我想这不是一次非常成功的经历。 虽然大家都知道坚果是非常健康和有营养的。矛盾的是,许多公司在实施测试自动化时没有考虑细微差别。每个人都知道它对身体有好处(就像吃坚果对我们健康有益一样&…

线上环境不要使用console.log,会导致页面卡顿,内存泄漏的原因

在我的这篇文章中分享了一个经验,就是vite打包之后会自动屏蔽所有的console.log打印的日志。 那么,在线上环境中不要使用 console.log 打印日志,已经成为共识,你知道是为什么么? 原因当然在于这个罪魁祸首 console 方…

深拷贝、浅拷贝的方法

目录 浅拷贝 深拷贝 与浅拷贝概念: 深浅拷贝出现的前提: 应用类型的数据(对象和数组) 深拷贝就是把目标对象里面的数据一个一个都复制下来,创建出一个一模一样的,存放地 址不一样,互不影响。…

【白嫖】如何底价续费服务器

目录 背景 问题 缓解方案 背景 现在各大云服务商的学生价服务器都已经关闭了,华为云、阿里云、百度云,以前都有学生价服务器,一年只要99。现在我找半天都没找到入口,而原价的一年得500块起步。。。 但是!&#xff0…

MySQL事务详解

目录 引例 什么是事务 一个完整事务所具有的四大属性 为什么会出现事务 事务常见操作方式 事务隔离级别 为什么要存在隔离级别 一致性 引例 如下图,是一个火车售票系统,当客户端A发现还有一张票时,将票卖掉,还没执行更新…

C++实现UDP可靠传输(二)

声明:禁止以任何形式转载本文章。本文章仅供个人学习记录与交流探讨,文章中提供的思路只是一种解决方案,代码也并非完整代码,如有需要,请自行设计协议并完成编程任务。 食用本文章之前,推荐阅读&#xff…

RKMEDIA--VENC/VDEC使用

前面两篇已经介绍了VI和VO的使用,本章节来介绍rkmedia且也是瑞芯微平台重点部分:编解码。 目录 一、简介 二、编解码能力 三、编码 编码初始化: 在初始化VENC时需要注意几点: venc编码帧率控制: 编码添加osd位图…

软件设计师考试重点1 计算机组成与体系结构

软件设计师考试重点1 计算机组成与体系结构一、 数据的表示1. 进制转换2. 原码/反码/补码/移码3. 数值表示范围4. 浮点数的运算二、运算器与控制器1. 计算机结构2. 计算机五大组成部分:3. CPU组成部分:4. 控制器组成部分:5. 运算器组成部分&a…

如何构建一个自动化油田注水站监控系统?

一、应用背景 目前,在大部分的钻井平台中,维护油田注水站的工作状态主要通过人工方式进行,这种方式不仅作业效率低且对工人的经验有着较高要求。此外,油田注水站的工作环境恶劣,为了能够有效地掌握各个设备的工作状态…

【测试沉思录】20. 如何做好测试需求分析?

作者:刘亚茹 编辑:毕小烦 我们都知道测试用例是软件测试中保障质量的必要手段,而测试需求作为用例编写的主要依据却往往被很多人忽视。到底什么是测试需求?又如何做好测试需求分析呢?本文带你了解一下。 1. 测试需求到…

Java项目如何导出数据为 PDF 文件?

文章目录Java项目如何导出数据为 PDF 文件?一、代码结构如下二、代码说明1、添加依赖 pom.xml3、添加字体4、PDF 导出工具类三、效果图结语Java项目如何导出数据为 PDF 文件? 一个小需求,需要将页面上的数据导出为PDF,正常情况下…

11月业务安全月报 | 台湾2300万人信息泄露;黑客两分钟即可破解安卓锁屏;乌克兰“IT军团”入侵俄罗斯中央银行

导语:随着数字化的深入普及,业务愈加开放互联。企业的关键数据、用户信息、基础设施、运营过程等均处于边界模糊且日益开放的环境中,涉及利益流和高附加值的业务面临多样的安全隐患,随时可能遭遇损失,进而影响企业运营…

Android 8.0网络DNS

1 Linux DNS规范 Linux上并没有一个单独的方法可以完成DNS查询工作;没有一个有这样的明确接口的核心系统调用system call。不过,glibc (nss)的getaddrinfo (3), gethostbyname (3)等相关API (RFC3493)提供了DNS查询功能。 1)不支持nscd&#…

关于使用STM32CubeMx配置串口出现的一些问题

一、使用CubeMX配置好串口中断的工程,打开工程添加相关代码后串口没有数据输出或者看不见数据但串口助手的RX会一直增加的问题。 参考:这里 大概意思就是时钟的原因,需要把stm32f4xx_hal_conf.h文件中的时钟改一下,将25MHz改为8…

Ansys Zemax | 使用 OpticStudio 进行闪光激光雷达系统建模(中)

在消费类电子产品领域,工程师可利用激光雷达实现众多功能,如面部识别和3D映射等。尽管激光雷达系统的应用非常广泛而且截然不同,但是 “闪光激光雷达” 解决方案通常都适用于在使用固态光学元件的目标场景中生成可检测的点阵列。凭借具有针对…

我要涨知识——TypeScript 经典高频面试题(二)

又是一个年底来了,好大一批人可能又准备跑路了,翻了翻掘金和 CSDN 发现好多大佬都有大厂 Offer ,看着看着我心动了! 话不多说,赶紧开干,给自己整了一个前端面试小助手——微信小程序内搜索 “WEB学习学习加…

jar启动指定JDK/JRE 安装路径教程

前言 因为疫情在家办公的缘故,有个老项目,需要改个接口,然后需要前端联调,因为外网服务器没有多余的空间了,想着把jar给前端让前端开发人员,在自己的本机启动后端服务,进行接口联调,…

streptavidin-PEG-6-FAM 链霉亲和素-聚乙二醇-6-羧甲基荧光素

产品名称:链霉亲和素-聚乙二醇-6-羧甲基荧光素 英文名称:streptavidin-PEG-6-FAM 纯度:95% 存储条件:-20C,避光,避湿 外观:固体或粘性液体,取决于分子量 PEG分子量可选:350、550、75…