Sherman-Morrison-Woodburg 定理

接触到这个定理的缘由是碰到矩阵秩一校正的问题,下面给出定理的内容。
定理
A n×n 非奇异矩阵, U n×m 矩阵, C m×m 非奇异矩阵, V m×n 矩阵,如果 C1+UA1V 非奇异,那么 A+UCV 也是非奇异的,且

(A+UCV)1=A1A1U(C1+UA1V)1VA1

1.证明:
U+UCVA1U=UC(C1+VA1U)=(A+UCV)A1U

(A+UCV)1UC=A1U(C1+VA1U)1

现在:
A1=(A+UCV)1(A+UCV)A1=(A+UCV)1(I+UCVA1)=(A+UCV)1+(A+UCV)1UCVA1=(A+UCV)1+A1U(C1+VA1U)1VA1

这个定理有一个有用的推论以及应用,下面我们稍微介绍一下。
推论
(Sherman-Morrison 定理)设 ARn×n 是非奇异矩阵, u,vRn 是任意向量,若 1+vTA1u0 ,那么 A 的秩一校正 A+uvT 非奇异,且其逆矩阵为

(A+uvT)1=A1A1uvTA11+vTA1u

在拟牛顿法的应用
定理经常在无约束最优化问题的拟牛顿法中,用来求解一个核心的方程的根。具体的说,给定一个向量值函数 f:RnRn ,利用牛顿法求解 f(x)=0 的一个近似解,是通过如下的迭代:

xk+1=xkJ(xk)1f(xk)

其中 J 表示 f 的雅可比矩阵,当然有时候 J 的性质可能并不是很好,我们通常对它进行一些校正,设 Ak J(xk) 的近似,在第(k+1)步迭代中,我们令

Ak+1=Akukvk

uk,vk 的选取是使得 Ak+1 近似 J(xk+1) Ak 近似 J(xk) 得更好。通常, uk,vk 是这样选取的,解如下方程:
Ak+1(xk+1xk)=fk+1fk

其中 xk+1=xkA1kfk,fk=f(xk) ,于是我们令
vTk=xk+1xk,uk=AkvTkgkvkvTk

其中 gk=fk+1fk .这个过程源于下面的事实:
使
(AkPk)(xk+1xk)=fk+1fk

成立的且使 ||Pk||2 最小的是一个秩一矩阵。
为了实现拟牛顿法,我们需要求得 A1K ,通过前面的定理,我们可以递归地求出
A1k+1=A1k+(vTkA1kgk)vkA1kvkA1kgk

Powell 发现利用上面的迭代公式,可以减少一定的计算误差。

你可能感兴趣的:(Sherman-Morrison-Woodburg 定理)