3.6 矩阵自反广义逆

醒过来摸鱼

已于 2023-05-19 14:44:05 修改

阅读量614

点赞数

分类专栏：线性代数【更新中】文章标签：矩阵线性代数算法

于 2023-02-07 11:12:24 首次发布

本文链接：https://blog.csdn.net/m0_66201040/article/details/128911882

版权

线性代数【更新中】专栏收录该内容

71 篇文章 38 订阅

订阅专栏

定义

自反广义逆reflexive generalized inverse是广义逆加上一个自反条件。广义逆 $A^-$ 不唯一，所以在这么多广义逆里可以加上各种条件缩小广义逆的范围，如果 $G$ 是 $A$ 的广义逆，那么再加上下面这个条件就成了自反广义逆：
$G A G = G$
自反广义逆的全体叫做 $A^{\{1,2\}}$ ,这个上标的1和2是下面两个条件:
$AGA=A\\ ~GAG=G$
从这里也不知道 $G$ 是 $A$ 的广义逆，但是A不一定是 $G$ 的广义逆，举个例子：
$A=\begin{pmatrix}0 & 0 & 1\\ 1 & 2 & -2\\ 2 & 4 & -2\\ \end{pmatrix}\\ G=\begin{pmatrix}2 & 5 & -2\\ 0 & -2 & 1\\ 1 & 0 & 0\\ \end{pmatrix}\\ AGA=\begin{pmatrix}0 & 0 & 1\\ 1 & 2 & -2\\ 2 & 4 & -2\\ \end{pmatrix}\\ GAG=\begin{pmatrix}-2 & 1 & 0\\ 2 & 0 & 0\\ 1 & 0 & 0\\ \end{pmatrix}\\$

计算方法

还记得计算广义逆的步骤吗？就是这一步：
$Q\begin{pmatrix} E_{r} & \ast\\ \ast & \ast \end{pmatrix}_{n \times m}P$
这一步中间的矩阵决定了是否为自反广义逆，如果中间这个矩阵和原矩阵的秩一样的话，那么得到的矩阵就是自反广义逆，也就是下面的公式：
$Q\begin{pmatrix} E_{r} & \ast\\ \ast & \bold 0 \end{pmatrix}_{n \times m}P$
如果右下角不为0，那么矩阵的秩就改变了，得到的就不是自反广义逆。因为我写广义逆的算法中把右下角全部变成零了，得到的就是自反广义逆，所以这篇文章我就不贴代码了。