• 矩阵分析与应用


    广义逆矩阵的定义与性质

    考虑m\times n阶的一般矩阵A,其秩k小于或等于min\left \{ m,n \right \}

    我们的问题是:是否存在某种合适意义下的逆矩阵,使得线性方程组Ax=y的解可以用这种逆矩阵表示?

    定义1:有线性方程A_{m\times n}x_{n\times 1}=y_{m\times 1},若矩阵A的行之间存在的线性关系也存在于向量y的对应元素中,则该方程组称为一致方程(即至少存在一个解能够严格满足该方程组)。

    如:

    x_{1}+x_{2}=4

    3x_{1}+3x_{2}=12

    矩阵A的第二行是第一行的3倍,向量y的第二个元素也是第一个元素的3倍。

    定义2:仅当线性方程组为一致方程时。这一线性方程组方可求解。

    定义3:当且仅当增广矩阵\left [ A,y \right ]的秩等于矩阵A_{m \times n}的值,即rank([A,y])=rank(A)时,线性方程A_{m\times n}x_{n\times 1}=y_{m\times 1}是一致方程。

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

    定义5:当且仅当AGA=A时,一致方程Ax=y对于y\neq 0有解x=Gy

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

    m*n矩阵A的广义逆矩阵G用符号A^{-}表示,即G=A^{-}

    广义逆矩阵有着以下重要性质。

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

    2.A^{-}存在\Leftrightarrow H= A^{-}A为幂等矩阵(即H^{2}=H)且rank(H)=rank(A)

    3.A^{-}存在\Leftrightarrow F=AA^{-}为幂等矩阵(即F^{2}=F)且rank(F)=rank(A)

    定义:m*n的矩阵A的广义逆矩阵是一个满足AA^{-}A=A的n*m的矩阵A^{-}

    定义:m*n的矩阵A的广义逆矩阵是满足下列条件之一的n*m矩阵A^{-}

    1.A^{-}A为幂等矩阵,且rank(A^{-}A)=rank(A)

    2.AA^{-}为幂等矩阵,且rank(AA^{-})=rank(A)

    逆矩阵,左逆矩阵,右逆矩阵都可以视为广义逆矩阵的一个特例。

    1.逆矩阵A^{-1}满足AA^{-1}A=A

    2.左逆矩阵L满足ALA=A,因为LA=I

    3.右逆矩阵R满足ARA=A,因为AR=I

    广义逆矩阵的计算

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

    命题:一个秩为r的m*n矩阵A可以分解为

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

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

    证明:存在m*m的非奇异矩阵P和n*n非奇异矩阵Q使得:


    PAQ=\begin{bmatrix} I_{r} &O_{r\times (n-r)} \\ O_{(m-r)\times r }& O_{(m-r)\times (n-r )} \end{bmatrix}

    或等价为:

    A=P^{-1}\begin{bmatrix} I_{r} &O_{r\times (n-r)} \\ O_{(m-r)\times r }& O_{(m-r)\times (n-r )} \end{bmatrix}Q^{-1}

    将逆矩阵P^{-1}Q^{-1}分块为:

    P^{-1}=[K_{m\times r},W_{m\times (n-r)}]

    Q^{-1}=\begin{bmatrix} L_{r\times n}\\Z_{(n-r)\times n} \end{bmatrix}


    于是,有:

    A=[K,W]\begin{bmatrix} I_{r} &O \\ O& O \end{bmatrix}\begin{bmatrix} L\\Z \end{bmatrix}=[K,O]\begin{bmatrix} L\\Z \end{bmatrix}=K_{m\times r}L_{r\times n}

    由于P是非奇异的,因此P^{-1}也是非奇异的,即P^{-1}的列是线性无关的,特别地,矩阵K的r列线性无关,固有rank(K)=r,即K具有满列秩。类似的,L具有满行秩。

    矩阵的满秩分解算法:

    步骤1:利用行初等变换将矩阵A化为阶梯型。

    A\rightarrow \begin{bmatrix} G_{r\times n}\\ O_{(m-r)\times n} \end{bmatrix}

    步骤2 对单位矩阵执行逆行初等变换,得到逆矩阵P^{-1}

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

    步骤4 书写满秩分解结果A=FG

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

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

    是A的一个广义逆矩阵。

    证明如下:

    AA^{-}A=AG^{T}(F^{T}AG^{T})^{-1}F^{T}A\\=AG^{T}((F^{T}F)(GG^{T}))^{-1}F^{T}A\\=AG^{T}(GG^{T})^{-1}(F^{T}F)^{-1}F^{T}A\\=FGG^{T}(GG^{T})^{-1}(F^{T}F)^{-1}F^{T}FG\\=FG\\=A

    广义逆矩阵的计算总结如下:

    1.计算矩阵的满秩分解A=FG

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

  • 相关阅读:
    一文了解超融合信创云基础设施
    【面试题精讲】如果一个类没有声明构造方法,该程序能正确执行吗?
    编译器一日一练(DIY系列之中间代码生成)
    vue3+element-plus国际化
    鸿鹄工程项目管理系统em Spring Cloud+Spring Boot+前后端分离构建工程项目管理系统
    Oracle/PLSQL: NumToDSInterval Function
    分组统计查询
    开源组件安全风险及应对
    Kubernetes体系结构及其主要组件
    linux luasocket 使用
  • 原文地址:https://blog.csdn.net/weixin_46007132/article/details/126010824