矩阵分析与应用-19-Moore-Penrose逆矩阵03

Xiao__fly

已于 2022-08-17 21:30:58 修改

阅读量373

点赞数

文章标签：矩阵算法 python

于 2022-08-12 13:11:32 首次发布

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

版权

非一致方程的最小范数最小二乘解

对丁非一致方程 $gif.latex?A_%7Bm%20%5Ctimes%20n%7Dx_%7Bn%5Ctimes%201%7D%3Dy_%7Bm%5Ctimes%201%7D$ ，矩阵 $gif.latex?G$ 称为 $gif.latex?A$ 的最小范数最小二乘广义逆矩阵，若 $gif.latex?G$ 满足条件

$gif.latex?%7C%7CGy%7C%7C_n%5Cleq%20%7C%7C%5Chat%7Bx%7D%7C%7C_n%20%2C%5Cforall%20%5Chat%7Bx%7D%5Cin%20%7B%5Chat%7Bx%7D%3A%20%7C%7CA%5Chat%7Bx%7D-y%7C%7C_m%20%5Cleq%20%7C%7CAz-y%7C%7C_m%20%5Cforall%20y%20%5Cin%20R%5Em%2Cz%5Cin%20R%5EN%20%7D$

其中， $gif.latex?%7C%7C%5Ccdot%20%7C%7C_n%2C%7C%7C%5Ccdot%7C%7C_m$ 是 $gif.latex?R%5En$ 和 $gif.latex?R%5Em$ 空间的范数(半范数):花括号{·}表示à是非一致方程 $gif.latex?Ax%3Dy$ 的最小二乘解，而 $gif.latex?%7C%7CGy%7C%7C_n%5Cleq%20%7C%7C%5Chat%7Bx%7D%7C%7C_n$ 。表示 $gif.latex?Gy$ 是在所有的最小二乘解中具有最小范数的那个解。

矩阵使得 $gif.latex?Gy$ 是非一致方程 $gif.latex?Ax%3Dy$ 的最小范数最小二乘解，当且仅当满足条件

$gif.latex?AGA%3DA%2C%28AG%29%5E%5C%23%3DAG%2CGAG%3DG%2C%28GA%29%5E%5C%23%3DGA$

其中， $gif.latex?A%5E%5C%23$ 是 $gif.latex?A$ 的伴随矩阵。

广义逆矩阵的阶数递推计算

左伪逆矩阵的阶数递推

考虑 $gif.latex?n%5Ctimes%20m$ 矩阵 $gif.latex?F_m$ ，并设 $gif.latex?F%5E%5Cdagger%3D%28F%5EH_mF_m%29%5E%7B-1%7DF%5EH_m$ 是 $gif.latex?F_m$ 的左伪逆矩阵。

令

$gif.latex?F_m%3D%5BF_%7Bm-1%7Df_m%5D$

$gif.latex?f_m$ 是矩阵 $gif.latex?F_m$ 的第 $gif.latex?m$ 列，且 $gif.latex?rank%28F_m%29%3Dm$ ，则计算 $gif.latex?F%5E%5Cdagger_m$ 的递推公式由下式给出:

$gif.latex?F%5E%5Cdagger_m%3D%5Cbegin%7Bbmatrix%7D%20F%5E%5Cdagger_%7Bm-1%7D-F%5E%5Cdagger_%7Bm-1%7Df_me%5EH_m%5CDelta%20_m%5E%7B-1%7D%5C%5C%20e%5EH_m%5CDelta_m%5E%7B-1%7D%20%5Cend%7Bbmatrix%7D$

式中

$gif.latex?e_m%3D%5BI_n-F_%7Bm-1%7DF%5E%5Cdagger_%7Bm-1%7D%5Df_m$

$gif.latex?%5CDelta_m%5E%7B-1%7D%3D%5Bf%5EH_me_m%5D%5E%7B-1%7D$

且初始值为

$%28f%5EH_1f_1%29$

右伪逆矩阵的阶数递推

考虑矩阵 $gif.latex?F_m%5Cin%20C%5E%7Bn%5Ctimes%20m%7D$ ， $gif.latex?n%20%3C%20m$

记 $gif.latex?F_m%3D%5BF_%7Bm-1%7D%2Cf_m%5D$ ，则右伪逆矩阵 $gif.latex?F%5E%5Cdagger%3DF%5EH_m%28F%5EH_mF_m%29%5E%7B-1%7D$ 具有以下递推公式:

$gif.latex?F%5E%5Cdagger_m%3D%5Cbegin%7Bbmatrix%7D%20F%5E%5Cdagger_%7Bm-1%7D-%5CDelta%20_mF%5E%5Cdagger_%7Bm-1%7Df_me_m%5C%5C%20%5CDelta_me%5EH_m%20%5Cend%7Bbmatrix%7D$