【矩阵论】4.矩阵运算——广义逆——定义性质与特殊矩阵的广义逆

news2025/7/17 20:13:17

在这里插入图片描述

4.3 广义逆

4.3.1 定义

若 m × n 矩 阵 A = A m × n 与 矩 阵 X = X n × m 满 足 四 个 条 件 ① A X A = A , ② X A X = X , ③ ( A X ) H = A X , ④ ( X A ) H = X A 则 X 为 A 的 加 号 逆 ( 广 义 逆 ) , 记 为 X = A + \begin{aligned} &若 m\times n 矩阵A=A_{m\times n} 与矩阵 X=X_{n\times m} 满足四个条件\\ &①AXA=A,\quad ② XAX=X,\quad ③(AX)^H=AX,\quad ④(XA)^H=XA\\ &则X为A的加号逆(广义逆),记为X=A^+ \end{aligned} m×nA=Am×nX=Xn×mAXA=A,XAX=X(AX)H=AX,(XA)H=XAXA(广)X=A+

a. 性质

  • A与X互为加号逆

  • A + A^{+} A+ 具有唯一性
    证 明 : ( 反 证 法 ) 假 设 X , Y 都 满 足 四 个 条 件 , 有 A X A = A , X A X = X , ( A X ) H = A X , ( A Y ) H = A Y , A Y A = A , Y A Y = Y , ( A Y ) H = A Y , ( Y A ) H = Y A 有 X = X A X = A = A Y A ( 替 换 右 侧 X ) X A Y A X = X ( A Y ) H ( A X ) H = X [ ( A X ) ( A Y ) ] H = A X A = A X ( A Y ) H = X A Y = A Y A = A ( 消 去 左 侧 X ) X A Y A Y = ( X A ) H ( Y A ) H Y = [ ( Y A ) ( X A ) ] H A = ( Y A ) H Y = Y 故 可 得 A + 具 有 唯 一 性 \begin{aligned} &证明:(反证法)\\ &假设X,Y都满足四个条件,有AXA=A,XAX=X,(AX)^H=AX,(AY)^H=AY,\\ &AYA=A,YAY=Y,(AY)^H=AY,(YA)^H=YA\\ &有X=XAX\xlongequal{A=AYA(替换右侧X)}XAYAX=X(AY)^H(AX)^H=X[(AX)(AY)]^H\\ &\xlongequal{AXA=A}X(AY)^H=XAY\xlongequal{AYA=A(消去左侧X)}XAYAY=(XA)^H(YA)^HY=\\ &[(YA)(XA)]^HA=(YA)^HY=Y\\ &故可得A^+具有唯一性 \end{aligned} ()XY,AXA=A,XAX=X,(AX)H=AX,(AY)H=AY,AYA=A,YAY=Y,(AY)H=AY,(YA)H=YAX=XAXA=AYA(X) XAYAX=X(AY)H(AX)H=X[(AX)(AY)]HAXA=A X(AY)H=XAYAYA=A(X) XAYAY=(XA)H(YA)HY=[(YA)(XA)]HA=(YA)HY=YA+

  • A + A^+ A+ 的H穿脱公式: ( A H ) + = ( A + ) H (A^H)^+=(A^+)^H (AH)+=(A+)H

  • 乘积的加号逆无穿脱公式: ( A B ) + ≠ B + A + (AB)^+\neq B^+A^+ (AB)+=B+A+

b. 广义逆特例

  • A = 0 A=0 A=0 ,则 A + = 0 A^{+}=0 A+=0

  • 若方阵 A = A n × n A=A_{n\times n} A=An×n 可逆 ( ∣ A ∣ ≠ 0 \vert A \vert\neq 0 A=0) ,则 A + = A − 1 A^+=A^{-1} A+=A1

  • A = ( a ) A=(a) A=(a) 为1阶阵,即 A = ( 复 数 a ) A=(复数a) A=(a) ,则有 ( a ) + = a + = { a − 1 , a ≠ 0 0 , a = 0 (a)^+=a^+=\left\{\begin{aligned}a^{-1},&a\neq 0\\0,&a=0\end{aligned}\right. (a)+=a+={a1,0,a=0a=0

  • 对角阵 D = ( a 1 a 2 ⋱ a n ) D=\left(\begin{matrix}a_1&&&\\&a_2&&\\&&\ddots\\&&&a_n\end{matrix}\right) D=a1a2an ,则 D + = ( a 1 + a 2 + ⋱ a n + ) D^+=\left(\begin{matrix}a_1^+&&&\\&a_2^+&&\\&&\ddots\\&&&a_n^+\end{matrix}\right) D+=a1+a2+an+

    eg:

    在这里插入图片描述

    D + D = D D + D^+D=DD^+ D+D=DD+ ,且 ( D + ) k = ( D k ) + (D^+)^k=(D^k)^+ (D+)k=(Dk)+ ,(k=1,2,…)

    证明:

    在这里插入图片描述

  • A为Hermite幂等阵,即 A H = A = A 2 A^H=A=A^2 AH=A=A2 ,则 A + = A A^+=A A+=A

  • 若A为列优阵( A H A = I A^HA=I AHA=I) ,则 A + = A H A^+=A^H A+=AH
    证 明 : 由 于 A 是 列 半 U 阵 , 则 A H A = I A H A A H = A H , A A H A = A , ( A H A ) H = I = A H A , ( A A H ) H = A A H \begin{aligned} &证明:由于 A 是列半U阵,则 A^HA=I\\ &A^HAA^H=A^H,AA^HA=A,(A^HA)^H=I=A^HA,(AA^H)^H=AA^H \end{aligned} :AUAHA=IAHAAH=AH,AAHA=A,(AHA)H=I=AHA,(AAH)H=AAH

    A H A^H AH 为列优阵(A为行优阵)( A A H = I AA^H=I AAH=I) ,则 A + = A H A^+=A^H A+=AH

    证 明 : 由 于 A 是 行 半 U 阵 , 则 A A H = I A H A A H = A H , A A H A = A , ( A H A ) H = A H A , ( A A H ) H = I = A A H \begin{aligned} &证明:由于 A 是行半U阵,则 AA^H=I\\ &A^HAA^H=A^H,AA^HA=A,(A^HA)^H=A^HA,(AA^H)^H=I=AA^H \end{aligned} :AUAAH=IAHAAH=AH,AAHA=A,(AHA)H=AHA,(AAH)H=I=AAH

  • 若A为U阵(单位正交阵),则 A + = A H = A − 1 A^+=A^H=A^{-1} A+=AH=A1

    若A为U阵,则 A H A = A A H = I , A H = A − 1 A^HA=AA^H=I,A^H=A^{-1} AHA=AAH=I,AH=A1

    A H A A H = A H , A A H A = A , ( A A H ) H = I = A A H , ( A H A ) H = A H A A^HAA^H=A^H,AA^HA=A,(AA^H)^H=I=AA^H,(A^HA)^H=A^HA AHAAH=AH,AAHA=A,(AAH)H=I=AAH,(AHA)H=AHA ,故 A + = A H A^+=A^H A+=AH

  • 若A为单阵(可相似对角化)且可逆 ,则有谱分解 A = λ 1 G 1 + ⋯ + λ k G k A=\lambda_1G_1+\cdots+\lambda_kG_k A=λ1G1++λkGk A + = λ 1 − 1 G 1 + ⋯ + λ k − 1 G k A^+=\lambda_1^{-1}G_1+\cdots+\lambda_k^{-1}G_k A+=λ11G1++λk1Gk

    若不满足可逆条件,则未必成立

  • 若A为正规阵(正规方阵),则 A + A^+ A+ 也为正规阵

    有正规分解 A = Q D D H A=QDD^H A=QDDH ,则 A + = Q D + Q H A^+=QD^+Q^H A+=QD+QH

    有谱分解 A = λ 1 + G 1 + ⋯ + λ k + G k A=\lambda_1^+G_1+\cdots+\lambda_k^{+}G_k A=λ1+G1++λk+Gk

c. A A + AA^+ AA+ A + A A^+A A+A 相关结论

  • A A + AA^+ AA+ A + A A^+A A+A 都是Hermite阵

  • r ( A + A ) = r ( A A + ) = r ( A ) r(A^+A)=r(AA^+)=r(A) r(A+A)=r(AA+)=r(A)

    由于矩阵的秩越乘越小,则 r ( A A + ) ≤ r ( A ) r(AA^+)\le r(A) r(AA+)r(A) A = A A + A A=AA^+A A=AA+A ,则 r ( A ) ≤ r ( A A + ) r(A)\le r(AA^+) r(A)r(AA+) ,故有 r ( A A + ) ≤ r ( A ) ≤ r ( A A + ) r(AA^+)\le r(A)\le r(AA^+) r(AA+)r(A)r(AA+)

  • 幂等性: ( A A + ) 2 = A A + (AA^+)^2=AA^+ (AA+)2=AA+ ( A + A ) 2 = A + A (A^+A)^2=A^+A (A+A)2=A+A

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

( 1 ) A = ( 1 0 2 1 ) 为 可 逆 阵 , 则 A + = A − 1 = ( 1 0 − 2 1 ) , A = ( a b c d ) = 1 ∣ A ∣ ( d − b − c a ) ( 2 ) ( 3 ) ( 4 ) 中 A 为 对 角 阵 \begin{aligned} &(1)A=\left( \begin{matrix}1&0\\2&1 \end{matrix} \right)为可逆阵,则A^+=A^{-1}=\left( \begin{matrix}1&0\\-2&1 \end{matrix} \right),A=\left( \begin{matrix} a&b\\c&d \end{matrix} \right)=\frac{1}{\vert A\vert}\left( \begin{matrix} d&-b\\-c&a \end{matrix} \right)\\ &(2)(3)(4)中A为对角阵 \end{aligned} (1)A=(1201)A+=A1=(1201),A=(acbd)=A1(dcba)(2)(3)(4)A

  • λ ( A A + ) 与 λ ( A + A ) \lambda(AA^+) 与 \lambda(A^+A) λ(AA+)λ(A+A) 都有r(A) 个1

    A ∈ C m , n , A + ∈ C n , m A\in C^{m,n},A^+\in C^{n,m} ACm,n,A+Cn,m 由于 A A + AA^+ AA+ A + A A^+A A+A 都是 Hermite 阵,且 A + A A^+A A+A 是幂等阵, ( A A + ) 2 = A A + , ( A + A ) 2 = A + A (AA^+)^2=AA^+,(A^+A)^2=A^+A (AA+)2=AA+,(A+A)2=A+A ,由Caley定理 x 2 = x x^2=x x2=x 则其0化式为 x 2 − x = 0 x^2-x=0 x2x=0 A A + AA^+ AA+ 有两个不同的根1和0

    已知 r ( A + A ) = r ( A A + ) = r ( A ) = r r(A^+A)=r(AA^+)=r(A)=r r(A+A)=r(AA+)=r(A)=r

    A A + AA^+ AA+ 有r个正根,m-r 个0根; A + A A^+A A+A 有r个正根,n-r 个0根

  • A + A A^+A A+A A A + AA^+ AA+ 都是半正定阵 ( A A + ≥ 0 , A + A ≥ 0 AA^+ \ge 0,A^+A\ge 0 AA+0,A+A0)

    已知 A + A A^+A A+A 为幂等阵,则 ( A A + ) 2 = ( A A + ) H ( A A + ) = P = A A + P H P ≥ 0 (AA^+)^2=(AA^+)^H(AA^+)\xlongequal{P=AA^+}P^HP\ge 0 (AA+)2=(AA+)H(AA+)P=AA+ PHP0 ,故 A A + AA^+ AA+ 为半正定阵,同理 A A + AA+ AA+ 为半正定阵

    由Caley定理, x 2 = x x^2=x x2=x ,则0化式 x ( x − 1 ) = 0 x(x-1)=0 x(x1)=0 ,故特根 ≥ \ge 0 ,为半正定阵

  • 一般有 A A + ≠ I AA^+\neq I AA+=I , A + A ≠ I A^+A\neq I A+A=I , 且 A A + ≠ A + A AA^+\neq A^+A AA+=A+A

  • 若 A 为正规阵,则 A + A = A A + A^+A=AA^+ A+A=AA+ ,且 ( A + ) k = ( A k ) + , k = 0 , 1 , 2 , ⋯ (A^+)^k=(A^k)^+,k=0,1,2,\cdots (A+)k=(Ak)+,k=0,1,2,
    由 于 A 是 正 规 阵 , 则 A A H = A H A , A 有 正 规 分 解 A = Q D Q H = Q ( λ 1 ⋱ λ n ) Q H , 其 中 Q 为 U 阵 ⇒ A + = ( Q D Q H ) + = ( Q H ) + D + Q + = Q + = Q H Q D + Q H , ⇒ ( A + ) k = ( Q D + Q H ) k = Q ( D + ) k Q H = Q ( D k ) + Q H = ( A k ) + 且 A + A = ( Q D + Q H ) Q D Q H = ( Q D + D Q H ) = Q D D + Q H = ( Q D Q H ) ( Q D + Q H ) = A A + \begin{aligned} &由于A是正规阵,则AA^H=A^HA,A有正规分解A=QDQ^H=Q\left( \begin{matrix} \lambda_1&&\\&\ddots&\\&&\lambda_n \end{matrix} \right)Q^H,\\ &其中Q为U阵 \\ &\Rightarrow A^+=(QDQ^H)^+=(Q^H)^+D^+Q^+\xlongequal{Q^+=Q^H}QD^+Q^H,\\ &\Rightarrow (A^+)^k=(QD^+Q^H)^k=Q(D^+)^kQ^H=Q(D^k)^+Q^H=(A^k)^+\\ &且A^+A=(QD^+Q^H)QDQ^H=(QD^+DQ^H)=QDD^+Q^H=(QDQ^H)(QD^+Q^H)=AA^+ \end{aligned} AAAH=AHA,AA=QDQH=Qλ1λnQH,QUA+=(QDQH)+=(QH)+D+Q+Q+=QH QD+QH,(A+)k=(QD+QH)k=Q(D+)kQH=Q(Dk)+QH=(Ak)+A+A=(QD+QH)QDQH=(QD+DQH)=QDD+QH=(QDQH)(QD+QH)=AA+

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

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

相关文章

以太坊:轻松理解EIP-4844

以太坊:轻松理解EIP-4844 以太坊网络在不断发展,多年来已经提出了许多技术提案。其中一个提案是由Vitalik Buterin提出的EIP-4844,它关注的是分片技术。 什么是EIP-4844,它将如何影响以太坊网络及其用户的未来?以下是…

体系结构28_多处理机(2)

对称式共享存储器体系结构 多处理机Cache一致性 **多个处理器共享一个存储器 **当处理器规模较小时,这种机器十分经济 **支持对共享数据和私有数据的Cache缓存 私有数据供一个单独的处理器使用,而共享数据供多个处理器使用 ** 共享数据进入Cache产生了一…

LeetCode 654.最大二叉树 617合并二叉树 700二叉搜索树中的搜索 98验证二叉搜索树

文章目录654最大二叉树c代码实现python 代码实现617合并二叉树c 代码实现python 代码实现700二叉搜索树中的搜索c代码实现python 代码实现98验证二叉搜索树c 代码实现python 代码实现654最大二叉树 给定一个不重复的整数数组 nums 。 最大二叉树 可以用下面的算法从 nums 递归…

[附源码]java毕业设计在线购物商城

项目运行 环境配置: Jdk1.8 Tomcat7.0 Mysql HBuilderX(Webstorm也行) Eclispe(IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持)。 项目技术: SSM mybatis Maven Vue 等等组成,B/S模式 M…

[附源码]java毕业设计疫情居家隔离服务系统

项目运行 环境配置: Jdk1.8 Tomcat7.0 Mysql HBuilderX(Webstorm也行) Eclispe(IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持)。 项目技术: SSM mybatis Maven Vue 等等组成,B/S模式 M…

阿里云与信通院邀您参与云原生安全用户调研

云原生安全建设已经成为企业云原生平台建设、应用云原生化改造进程中的必备项。在云原生的大势所趋下,云的分布式架构、容器化部署、边界消失等特点也带来了有别于传统硬件安全的模式:镜像漏洞、容器逃逸、租户隔离等安全问题威胁着企业的云原生平台和应…

Unity动态创建Avatar骨骼映射

目录前言1 了解Avatar骨骼映射2 动态创建Avatar骨骼映射2.1 寻找相关APIAvatarBuilder.BuildHumanAvatarDeclarationParametersReturnsDescriptionHumanDescriptionDescriptionProperties2.2 创建Skeleton数据2.3 创建Human映射关系2.4 创建Avatar3 总结前言 为了让单个动画可…

Linux(基于Centos7)(四)

文章目录一、任务目标二、任务资讯三、任务实施3-1.RPM软件包管理3-2.YUM方式安装软件一、任务目标 实施该工单的任务目标如下: 知识目标 1.了解RPM提供的功能。 2.了解YUM相对于RPM所具有的优点。 能力目标 1.能够通过RPM安装及管理软件包。 2.能够通过YUM安装及管…

MCE | ATM 激酶活化变单体后的神奇开挂!

ATM (Ataxia-telangiectasia mutated proteins) 是一种丝氨酸-苏氨酸蛋白激酶,是 DNA 损伤应答 (DDR) 的关键调节因子。ATM 是位名副其实的“指挥官”,在 DNA 双链断裂 (DSB) 中,参与细胞周期检查点维护、DNA 损伤修复和端粒维护等&#xff0…

基于Android的个人健康管理系统

目 录 基于Android的个人健康管理系统 Personal Health Management System Based On Android 1 引言 1 1.1 课题背景 1 1.2 编写目的 1 1.3 关于Android 1 1.4 关于MVC框架 3 2 可行性研究 6 2.1 技术可行性 6 2.2 经济可行性 6 2.3 时间可行性 6 3 需求分析 7 4 总体设计 8 …

栈——算法专项刷题(六)

六、栈 6.1后缀表达式 原题链接 根据 逆波兰表示法,求该后缀表达式的计算结果。 有效的算符包括 、-、*、/ 。每个运算对象可以是整数,也可以是另一个逆波兰表达式。 说明: 整数除法只保留整数部分。给定逆波兰表达式总是有效的。换句话…

【kafka】五、kafka工作流程

kafka工作流程 工作流程 kafka中消息是以topic进行分类的,生产者生产消息,消费者消费消息,都是面向topic的。 topic是逻辑上的概念,而partition是物理上的概念,每一个partition对应一个log文件,该log文件…

基建融资与预算软约束(2015年A股大牛市与“流动性堰塞湖”)-中国视角下的宏观经济

基建融资与预算软约束(2015年A股大牛市与“流动性堰塞湖”) – 潘登同学的宏观经济学笔记 文章目录基建融资与预算软约束(2015年A股大牛市与“流动性堰塞湖”) -- 潘登同学的宏观经济学笔记2015年A股大牛市定向宽松的货币政策导向定向宽松的货币政策原因投资项目的预算软约束政…

python基础之循环嵌套

文章目录一、break和continue二、while的循环嵌套知识点print:例题1打印星星:例题2打印九九乘法表一、break和continue break当某一条件满足时直接跳出当前循环; continue当某一条满足时,不执行后续循环重新开始新一轮循环; i1 …

[附源码]java毕业设计婴幼儿玩具共享租售平台

项目运行 环境配置: Jdk1.8 Tomcat7.0 Mysql HBuilderX(Webstorm也行) Eclispe(IntelliJ IDEA,Eclispe,MyEclispe,Sts都支持)。 项目技术: SSM mybatis Maven Vue 等等组成,B/S模式 M…

python基础之字典

文章目录一、字典1.dictionary(字典)2.和列表的区别:二、实例三、应用场景一、字典 1.dictionary(字典) 是除了列表之外最灵活的数据类型,字典同样可以用来存储多个数据(通常用于存储描述一个物体的相关信息); 2.和…

数据分析软件的使用

一 数据分析概述 1 概念 数据分析是利用数学,统计学理论相结合的科学统计分析方法,对Excel数据,数据库中的数据,收集的大量数据,网页抓取的数据进行分析,从中提前有价值的信息并形成结论进行展示的过程。…

用于调整PID控制器增益的遗传算法的实现(Matlab代码实现)

🍒🍒🍒欢迎关注🌈🌈🌈 📝个人主页:我爱Matlab 👍点赞➕评论➕收藏 养成习惯(一键三连)🌻🌻🌻 🍌希…

第八章 兼容多种模块标准的软件包封装

第八章 如何封装兼容多种JS模块标准的软件包? 为了方便用户使用,一款成熟的类库都会提供多种模块封装形式,比如大家最常用到的 Vue,就提供了cjs、esm、umd 等多种封装模式,并且还会提供对应的压缩版本,方便…

Python异或运算符示例

目录 异或 异或的性质 示例1:值交换 示例2:找出现一次的元素 示例2代码: 异或 英文为exclusive OR,缩写成xor,符号是^ aba^b0假0假0假0假1真1真1真0假1真1真1真0假 异或的性质 1、a ^ a0 任何数字和自己异或结…