一 广义逆矩阵A-
1) 定义:设A为n×m矩阵,秩R(A)= r<min(m,n), 满足如下方程的G
AGA = A (4-2-1)
定义为A的广义逆,G为m×n矩阵,并记为A-,一般不唯一,称为A-型广义逆。
仅当A为m阶非奇异方阵,凯利逆A-必然满足定义式(
2) A-型广义逆有如下性质:
(1) (其中之一,即
)
(2)
(3)
(4)
(4)若矩阵P正定,则
(6) G为ATA的广义逆,则GT也是ATA的广义逆。
3)广义逆A- 的计算
A-的计算有许多种方法,这里仅介绍一种常用的简便方法。
当的秩R(A) = r <min(n,m), 可得矩阵
A分块写成
其中R(A11)= r, A11为非奇异方阵。则其广义逆为
(4-2-2)
在证明此A-为A的广义逆之前,先证明如下等式:
(4-2-3)
根据A的奇异性质,R(A)=R(A11)= r,故有
(A
此式说明 (A,则(
可见(
例4-1. 设有矩阵
R(A) = 2
取
, 则
== A
二 广义逆 A+
1.定义
如果对A-作某些限制,就可得到一种唯一的广义逆,称为伪逆,并用A+表示。A+定义伪满足下列四个方程:
(4-2-4)
的广义逆。伪逆A+也称为Moore-Penrose广义逆。
A+唯一,证明如下。
设G1和G2为两个A+,按定义(
由此,由
广义逆A+也是一个A-,是一个同时满足(
2.广义逆A+的计算
在一般情况下,,在测量计算A+常用如下方法:
(1) 当A为对角阵时,则有
,
(4-2-4)
例4-2
设,则
(2) (4-2-6)
证:
(
例4-3
,
(3)当N为对称方阵,则由(
N+=N(NN)-N(NN)-N (4-2-7)
(4)设N=ATA,则
A+=N+AT (4-2-8)
证: A
N+AT AN+AT=
(AN+AT)T= AN+AT
(N+ATA)T=( N+N)T=( N(NN)-N(NN)-NN)T
= (N(NN)-N)T= N(NN)-N=N+ATA
三 相容线形方程组的广义逆解
广义逆理论的提出与解线形方程组理论有关,测量平差的基本计算归根到底是解线形方程组问题。所以下面介绍用广义逆解线形方程组的几本方法。
设由相容线形方程组
(4-2-9)
其解存在,设为,则有
按A-定义,有
或 (4-2-10)
此式称为方程组(
(4-2-11)
为(
其次方程组AX=O 的一般解为
(4-2-12)
I为单位阵,M为任意向量。事实上
因此 非其次方程组的一般解等于它的任一特解与对应其次方程组一般解之和,故(
(4-2-13)
这是相容线形方程组解的一般公式。
当A为m阶非奇异方阵时,则A- = A-1,即A-等于凯利逆,此时
(4-2-14)