矩阵分析与应用(11)

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

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

        前面学习了一致方程的最小范数解和非一致方程的最小二乘解,现在要寻找一个非一致方程最小二乘解中的具有最小范数的解,这样的解称为非一致方程的最小范数最小二乘解,也称半范数最小二乘解。

1. 定义

        对于非一致方程 A_{m\times n}x_{n\times 1}=y_{m\times 1} ,当矩阵 G 满足条件

\left \| Gy \right \|_n \leqslant \left \| \hat{x} \right \|_n \quad \forall \hat{x}\in \left \{ \hat{x}:\left \| A \hat{x}-y \right \|_m \leqslant \left \| Az-y \right \|_m\forall y\in R^m ,z \in R^n\right \}

时,称矩阵 G 为 A 的最小范数最小二乘广义逆矩阵。式中 \left \| \cdot \right \|_n 和 \left \| \cdot \right \|_m 分别是在 R^n 和 R^m 空间的范数(半范数); \left \{ \cdot \right \} 表示 \hat{x} 是非一致方程 Ax=y 的最小二乘解,而 \left \| Gy \right \|_n\leqslant \left \| \hat{x} \right \|_n 表示 Gy 是在所有的最小二乘解中具有最小范数的解。

2. 定理

        矩阵 G 使得 Gy 是非一致方程 Ax=y 的最小范数最小二乘解,当且仅当 G 满足条件 

AGA=A,\quad (AG)^\#=AG,\quad GAG=G,\quad (GA)^\#=GA

式中, A^\# 表示 A 的伴随矩阵。

         由于伴随矩阵具有性质 B^\#=B^H ,所以式中 (AG)^\#=AG 即 (AG)^H=AG ;(GA)^\#=GA 即 (GA)^H=GA 。因此,上述定理也可以等价描述为:

        矩阵 G 使得 Gy 是非一致方程 Ax=y 的最小范数最小二乘解,当且仅当 G 是 A 的 Moore-Penrose 逆矩阵。

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

1. 左伪逆矩阵的阶数递推

        考虑 n\times m 矩阵 F_m ,并设 F^\dagger =(F_M^HF_M)^{-1} F_m^H 是 F_m 的左伪逆矩阵。

        定理:令

F_m=\left [ F_{m-1}f_m \right ]

式中, f_m 是矩阵 F_m 的第 m 列,且 rank(F_m)=m ,则计算 F_m^\dagger 的递推公式为:

F_m^\dagger =\begin{bmatrix} F_{m-1}^\dagger -F_{m-1}^\dagger f_me_m^H \Delta_m^{-1} \\ e_m^H \Delta_m^{-1} \end{bmatrix}

式中

e_m=\left [ I_n-F_{m-1}F_{m-1}^\dagger \right ]f_m

\Delta_m^{-1}=\left [ f_m^He_m \right ] ^{-1}

且初始值为

F_1^\dagger=f_1^H/(f_1^Hf_1)

2. 右伪逆矩阵的阶数递推

        考虑矩阵 F_m \in C^{n\times m},n< m 。

        定理:记 F_m=\left [ F_{m-1},f_m \right ] ,则右伪逆矩阵 F_m^\dagger =F_m^H(F_mF_m^H)^{-1} 具有以下递推公式:

F_m^\dagger =\begin{bmatrix} F_{m-1}^\dagger -\Delta_mF_{m-1}^\dagger f_me_m \\ \Delta_me_m^H \end{bmatrix}

式中

 e_m^H=f_m^H(I_n-F_{m-1}F_{m-1}^\dagger),\Delta_m=e_m^Hf_m 

递推的初始值为

 F_1^\dagger =f_1^H/(f_1^hf_1) 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

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

余额充值