矩阵分析与应用-18-Moore-Penrose逆矩阵02

最新推荐文章于 2024-02-01 10:30:20 发布

Xiao__fly

最新推荐文章于 2024-02-01 10:30:20 发布

阅读量1k

点赞数

文章标签：矩阵线性代数 numpy

版权声明：本文为博主原创文章，遵循 CC 4.0 BY-SA 版权协议，转载请附上原文出处链接和本声明。

本文链接：https://blog.csdn.net/Xiao__fly/article/details/126280172

版权

Moore-Penrose 逆矩阵的计算

方程求解法

第一步:分别求解矩阵方程得到 $X^H,Y$ .

$AA^HX^H=A$

$A^HAY=A^H$

第二步:计算广义逆矩阵 $A^\dagger=XAY$

若矩阵A 为Hermitian矩阵，上述方法可以简化，因为以上两个矩阵方程等价为一个矩阵方程

$A^2X^H=A$

若 $A^H=A$ ,Moore-Penrose逆矩阵可以计算为

$A^\dagger=XAX^H$

算法1

步骤1计算矩阵 $B=AA^H$ 。

步骤2求解矩阵方程 $B^2X^H=B$ 得到矩阵 $X^H$ 。

步骤3计算B的Moore-Penrose逆矩阵 $B^\dagger=(AA^H)^\dagger=XBX^H$ 。

步骤4计算矩阵A的Moore-Penrose逆矩阵 $A^\dagger=A^H(AA^H)^\dagger=A^HB^\dagger$ 。

算法2

步骤1计算矩阵 $B=AA^H$ 。

步骤2求解矩阵方程 $B^2X^H=B$ 得到矩阵 $X^H$ 。

步骤3计算的Moore-Penrose逆矩阵 $B^\dagger=(A^HA)^\dagger=XBX^H$ 。

步骤4计算矩阵A的Moore-Penrose逆矩阵 $A^\dagger=(A^HA)^\dagger A^H=B^\dagger A^H$ 。

KL分解法

若A=KL是矩阵Amxn的满秩分解，则

$G=L^H(K^HAL^H)^{-1}K^H$

递推法

算法3

初始值 $A_1^\dagger=a_1^\dagger=(a_1^Ha_1)^{-1}a_1^H$ 。

递推令k= 2,3,…- ,n，进行以下计算:

$d_k=A^\dagger_{k-1}a_k$

$b_k=\begin{cases} (1+d_k^Hd_k)^{-1}d_k^HA_{k-1}^\dagger & d_k^Hd_k\neq -1 \\ (a_k-A_{k-1}d_k)^\dagger &d_k^Hd_k= -1 \end{cases}$

$A_k^\dagger=\begin{bmatrix} A^\dagger_{k-1}d_kb_k\\ b_k \end{bmatrix}$

迹方法

步骤1计算矩阵 $B=AA^H$ 。

步骤2令 $C_1=I$ 。

步骤3计算

$C_{i+1}=\frac{1}{i}tr(C_iB)I-C_iB,i=1,2,...,r-1$

步骤4计算

$A^\dagger=\frac{r}{tr(C_iB)}C_iA^T$

注： $C_{i+1}B=O,tr(C_iB)\neq 0$ .

关注

0
点赞
踩
2

收藏

觉得还不错? 一键收藏
0
评论
复制链接

分享到 QQ

分享到新浪微博

扫一扫

评论

被折叠的条评论为什么被折叠?

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

1.余额是钱包充值的虚拟货币，按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载，可以购买VIP、付费专栏及课程。