码农知识堂 - 1000bd
  •   Python
  •   PHP
  •   JS/TS
  •   JAVA
  •   C/C++
  •   C#
  •   GO
  •   Kotlin
  •   Swift
  • 矩阵分析与应用-16-广义逆矩阵


    目录

    左逆矩阵与右逆矩阵

    广义逆矩阵的定义及性质

    广义逆矩阵的计算

    矩阵的满秩分解算法

    广义逆矩阵的计算

    一致方程的最小范数解

    非一致方程的最小二乘解

    左逆矩阵与右逆矩阵

    满足LA=I,但不满足AL=I的矩阵L称为矩阵A的左逆矩阵(left inverse)。

    满足AR=I,但不满足RA=L的矩阵称为矩阵A的右逆矩阵(right inverse)o

    仅当m \leq n时,矩阵A \in C^{m\times n}T可能有右逆矩阵。

    广义逆矩阵的定义及性质

    线性方程A_{m\times n}x_{n\times 1}=y_{m \times 1}称为一致方程,若矩阵A的行之间存在的线性关系也存在于向量y 的对应元素中。

    令A是一个m x n矩阵,具有任意秩。矩阵A的广义逆矩阵是-个n x m矩阵G,并使得当Ax=y 为一致方程时,x=Gy是线性方程Ax=y的解。

    一致方程Ax=y对y≠0有解x=Gy,当且仅当AGA = A。

    m × n矩阵A的广义逆矩阵是一个满足

    AA^-A=A

    的nx m矩阵A^-。

    m x n矩阵A的广义逆矩阵是满足下列条件之一的n × m 矩阵A^-

     (1)A^-A为幂等矩阵,且 rank(.A^-A) rank(A):

    (2)AA^-是幂等矩阵,且rank(AA^-) =rank(A)。

    性质:

    方程Ax=0的解与矩阵A的任意行正交,并且线性无关。

    A^{-}存在\Leftrightarrow AA^-A=A

    定理1:一线性方程组可以求解,当且仅当这些方程为一致方程。

    定理2:线性方程Ax=y 是一致的,当且仅当增广矩阵|A,y|的秩等于矩阵A_{m\times n}的秩,即rank(|A,y|)=rank(A)。

    广义逆矩阵的计算

    令A_{m\times n}具有秩r。若A=FG,其中,F_{m \times r}的秩为r (满列秩矩阵),且 G_{r\times n}的秩也为r(满行秩矩阵),则称A=FG为矩阵A的满秩分解(full-rank decomposition)。

    一个秩为r的mx n矩阵A可以分解为

    A=K_{m\times r}L_{r\times n}

    式中,K和L分别具有满列秩和满行秩。

    若矩阵A_{m\times n}具有秩r,且其满秩分解为A=FG,其中,F_{m \times r}为满列秩G_{r\times n}为满行秩,则

     A^{-}=G^T(F^TAG^T)^{-1}F^T

     是A的一个广义逆矩阵。 

    矩阵的满秩分解算法

    1.利用行初等变换将矩阵A化为阶梯型

    A\xrightarrow[E_1]{}[]\xrightarrow[E_2]{}...\xrightarrow[E_k]{}\begin{bmatrix} G_{r\times n}\\ O_{(m-r)\times n} \end{bmatrix}

    2.对单位矩阵执行逆行初等变换,得到逆矩阵

    I\xrightarrow[E_k^{-1}]{}[ ]\xrightarrow[E_{k-1}^{-1}]{}....\xrightarrow[E_1^{-1}]{}[]=P^{-1}

    3.利用逆矩阵P^{-1}的前r列构造矩阵F。

    4.书写满秩分解结果A=FG。

    广义逆矩阵的计算

    1.计算矩阵A_{m\times n}的满秩分解A=FG。

    2.求广义逆矩阵A^{-}=G^T(F^TAG^T)^{-1}F^T。

    一致方程的最小范数解

    令n x m 矩阵A^-是m×n矩阵A的任意一个广义逆矩阵,则

    (1)齐次方程Ax=0的一个通解为x=(I-A^-A)z,其中,z是n×1任意向量。

    (2非齐次方程Ax=y为-致方程的充分必要条件是

    AA^-y=y

    (3)非齐次方程Aa = g 的一个通解为

    x=A^-y+(I-A^-A)z

    式中,z为n×1任意向量。

    伴随矩阵具有以下性质

    (1)(A^\#)^\#=A

    (2)(AB)^\#=B^\#A^\#

    (3)<Ax,By>=0,\forall x,y\Leftrightarrow A^\#B=O

    (4)A^\#=A^T(若A为实矩阵)或A^\#=A^H(若A为复矩阵)

    Cy是一致方程Ax=y的最小范数解,当且仅当

    AGA=A,(GA)^\#=GA

    非一致方程的最小二乘解

    令G是某个矩阵,则\hat{x}=Gy是非一致方程Ax=y的最小二乘解,当且仅当

    A^\#AG=A^\#

    或

    AGA=A,(AG)^\#=AG

  • 相关阅读:
    搜广推火线入门
    MySQL优化(1):B+树与索引
    漏洞复现 - - -Docker未授权访问漏洞
    基于Python的电子病历实体识别系统
    功率放大器主要性能指标是什么(功率放大器工作状态的分类)
    【工具】电脑网络连接正常,但是有些页面无法登录,如何解决?
    Matlab----程序流程控制
    谈思生物直播—GENOVIS张洪妍抗体特异性酶切技术助力抗体药物结构表征
    [JavaEE系列] Thread类的基本用法
    我试图给你分享一种自适应的负载均衡。
  • 原文地址:https://blog.csdn.net/Xiao__fly/article/details/126234564
  • 最新文章
  • 沪漂五周年了:我越来越迷茫了
    Agentic Skill Routing 实战:别再把所有 Skill 塞进 AI Agent 上下文
    MySQL-Seconds_behind_master的精度误差
    [MAF预定义ChatClient中间件-03]CachingChatClient——利用缓存省钱省时间
    AI的至暗历史:从万众期待到被政府撤资,AI的两次死亡徘徊
    Agent OS :五种驯服不确定性的范式
    PortSwigger SQL注入LAB11
    数据库即时编译JIT
    [Begin]AI Learn Data Day 0
    深度学习进阶(二十七)现代 LLM 的核心架构设计其二:SwiGLU
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
小工具 小游戏
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1

京公网安备 11010502049817号