矩阵分析与应用(10)

Butterfffly

已于 2022-07-08 21:36:31 修改

阅读量398

点赞数

分类专栏：矩阵分析与应用文章标签：线性代数

于 2022-07-07 19:02:23 首次发布

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

本文链接：https://blog.csdn.net/qq_40206924/article/details/125663612

版权

矩阵分析与应用专栏收录该内容

22 篇文章 9 订阅

订阅专栏

本文详细介绍了Moore-Penrose逆矩阵的四种计算方法：方程求解法、KL分解法、递推法和迹方法。对于不同情况，如矩阵的行数与列数的关系，可以选择合适的算法进行计算。此外，还提供了具体步骤和简化公式，适用于Hermitian矩阵的特殊情况。

摘要由CSDN通过智能技术生成

学习来源：《矩阵分析与应用》张贤达清华大学出版社

Moore-Penrose 逆矩阵的计算

1. 方程求解法

第一步：分别求解矩阵方程

$AA^HX^H=A$

$A^HAY=A^H$

得到 $X^H$ 和 $Y$ 。

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

若矩阵 $A$ 为 Hermitian 矩阵，则上述方程可以简化为

$A^2X^H=A$

$A$ 的广义逆矩阵 $A^\dagger =XAX^H$ 。由此，可以得到两种算法：

算法一：

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$ 。

算法二：

1）计算矩阵 $B=A^HA$ 。

2）解矩阵方程 $B^2X^H=B$ 得到矩阵 $X^H$ 。

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

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

若矩阵 $A_{m\times n}$ 的列数大于行数，则矩阵乘积 $AA^H$ 的维数比 $A^HA$ 的维数小，所以选择算法一；反之，选择算法二。

2. KL分解法

在矩阵满秩分解中提到了 Moore-Penrose 逆矩阵的一种方法。若 $A=KL$ 是矩阵 $A_{m\times n}$ 的满秩分解，则

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

满足 Moore-Penrose 逆矩阵的四个条件，所以 $G$ 是 $A_{m\times n}$ 的 Moore-Penrose 逆矩阵。

3. 递推法

对矩阵 $A_{m\times n}$ 的前 $k$ 列进行分块 $A_k=\left [ A_{k-1},a_k \right ]$ ，其中 $a_k$ 是矩阵 $A$ 的第 $k$ 列。因此，分块矩阵 $A_k$ 的 Moore-Penrose 逆矩阵 $A_k^\dagger$ 可以由 $A_{k-1}^\dagger$ 递推计算。

算法：

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

递推：令 $k=2,3,\cdots ,n$ ，有

$d_k=A_{k-1}^\dagger a_k$

$b_k=\left\{\begin{matrix} (1+d_k^Hd_k)^{-1}d_k^HA_{k-1}^\dagger , & \quad d_k^Hd_k\neq -1\\ (a_k-A_{k-1d_k}), &\quad d_k^Hd_k=-1 \end{matrix}\right.$

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

4. 迹方法

已知矩阵 $A_{m\times n}$ 的秩为 $r$ ，求 Moore-Penrose 逆矩阵的迹方法为：

1）计算 $B=AA^T$ 。

2）令 $C_1=I$ 。

3）计算

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

4）计算

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

其中 $C_{i+1}B=O$ ，且 $tr(C_iB)\neq 0$ 。

关注

0
点赞
踩
0

收藏

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

分享到 QQ

分享到新浪微博

扫一扫

专栏目录

评论

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

到【灌水乐园】发言

查看更多评论

添加红包

成就一亿技术人!

hope_wisdom

发出的红包

实付元

使用余额支付

点击重新获取

扫码支付

钱包余额 0

抵扣说明：

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