码农知识堂 - 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

  • 相关阅读:
    【虹科案例】极高的精度水平——虹科数字化仪在大型强子对撞机机器保护系统中的应用
    VS2017下第一个C#控制台程序helloworld
    Java多线程(5):CAS
    我在滴滴做开源
    解决SQL错误(208):对象名‘string_split‘无效的正确方法,亲测有效!!!
    ssh连win10报错:Permission denied (publickey,keyboard-interactive).
    单链表OJ题(2):反转链表(三指针法)、找中间节点(快慢指针)
    ORACLE 基础
    iOS中的锁
    国外发达国家码农是真混得好么?
  • 原文地址:https://blog.csdn.net/Xiao__fly/article/details/126234564
  • 最新文章
  • 攻防演习之三天拿下官网站群
    数据安全治理学习——前期安全规划和安全管理体系建设
    企业安全 | 企业内一次钓鱼演练准备过程
    内网渗透测试 | Kerberos协议及其部分攻击手法
    0day的产生 | 不懂代码的"代码审计"
    安装scrcpy-client模块av模块异常,环境问题解决方案
    leetcode hot100【LeetCode 279. 完全平方数】java实现
    OpenWrt下安装Mosquitto
    AnatoMask论文汇总
    【AI日记】24.11.01 LangChain、openai api和github copilot
  • 热门文章
  • 十款代码表白小特效 一个比一个浪漫 赶紧收藏起来吧!!!
    奉劝各位学弟学妹们,该打造你的技术影响力了!
    五年了,我在 CSDN 的两个一百万。
    Java俄罗斯方块,老程序员花了一个周末,连接中学年代!
    面试官都震惊,你这网络基础可以啊!
    你真的会用百度吗?我不信 — 那些不为人知的搜索引擎语法
    心情不好的时候,用 Python 画棵樱花树送给自己吧
    通宵一晚做出来的一款类似CS的第一人称射击游戏Demo!原来做游戏也不是很难,连憨憨学妹都学会了!
    13 万字 C 语言从入门到精通保姆级教程2021 年版
    10行代码集2000张美女图,Python爬虫120例,再上征途
Copyright © 2022 侵权请联系2656653265@qq.com    京ICP备2022015340号-1
正则表达式工具 cron表达式工具 密码生成工具

京公网安备 11010502049817号