Woodbury矩阵恒等式证明及推论

Woodbury矩阵恒等式又称矩阵求逆引理,由该引理可以证明push-through矩阵恒等式、Sherman-Morrison 定理等,下面给出Woodbury矩阵恒等式的具体形式。

     {(A + BCD)^{ - 1}} = {A^{ - 1}} - {A^{ - 1}}B{({C^{ - 1}} + D{A^{ - 1}}B)^{ - 1}}D{A^{ - 1}}

其中,{(A + BCD)^{ - 1}}AC{C^{ - 1}} + D{A^{ - 1}}B可逆,为了方便,我们设Am \times mBm \times nCn \times nDn \times m

下面给出证明:

首先考虑{A + BCD},我们将右乘{A^{ - 1}},有:

                    (A + BCD){A^{ - 1}} = I + BCD{A^{ - 1}}

为了和右式的{C^{ - 1}} + D{A^{ - 1}}B建立联系,我们右乘B,有:

                 (I + BCD{A^{ - 1}})B = B + BCD{A^{ - 1}}B

由于C可逆,所以:

              (A + BCD){A^{ - 1}}B = BC({C^{ - 1}} + D{A^{ - 1}}B)

由于{C^{ - 1}} + D{A^{ - 1}}B可逆,所以:

            BC = (A + BCD){A^{ - 1}}B{({C^{ - 1}} + D{A^{ - 1}}B)^{ - 1}}

为了与Woodbury恒等式的左边相联系,同时发现右边的一项右乘了DA^{-1},自然的想法就是,我们将BC配成A+BCD,有:

    A+BCD = A+(A + BCD){A^{ - 1}}B{({C^{ - 1}} + D{A^{ - 1}}B)^{ - 1}}D

再右乘A^{-1},就出现了等式右边的那一项:

(A + BCD){A^{ - 1}} = I + (A + BCD){A^{ - 1}}B{({C^{ - 1}} + D{A^{ - 1}}B)^{ - 1}}D{A^{ - 1}}

现在就很显然了,为了匹配右边的那一项,我们肯定要左乘(A+BCD)^{-1},于是有:

      {A^{ - 1}} = {(A + BCD)^{ - 1}} + {A^{ - 1}}B{({C^{ - 1}} + D{A^{ - 1}}B)^{ - 1}}D{A^{ - 1}}

移项即可,有:

       {(A + BCD)^{ - 1}}= {A^{ - 1}} - {A^{ - 1}}B{({C^{ - 1}} + D{A^{ - 1}}B)^{ - 1}}D{A^{ - 1}}

证明完毕。

一个简单的记忆和推导方式

只要清楚(A+BCD)^{-1}的逆是A^{-1}+X的形式,通过解关于X的方程即可推出Woodbury矩阵恒等式。下面给出某博主的求解思路:

矩阵求逆引理(matrix inversion lemma)冬瓜班小朋友的博客-CSDN博客矩阵求逆引理

特殊情形

● 当AC是单位阵时,Woodbury矩阵恒等式可以变成:

               {(I + BD)^{ - 1}} = I - B{(I + DB)^{ - 1}}D

这个等式可以让我们联想到:
         {(I + P)^{ - 1}} = I - P{(I + P)^{ - 1}} = I - {(I + P)^{ - 1}}P

以及push-through等式:

                   {(I + BD)^{ - 1}}B = B{(I + DB)^{ - 1}}

关于这个不等式的证明只需左乘I+BD,提取公因式,化简即可。当然Woodbury矩阵恒等式也可以通过这两个等式配凑得到。

● 当BD是单位阵时,Woodbury矩阵恒等式可以变成:

          {(A + C)^{ - 1}} = {A^{ - 1}} - {A^{ - 1}}{({C^{ - 1}} + {A^{ - 1}})^{ - 1}}{A^{ - 1}}

化简可以得到:

               {(A + C)^{ - 1}} = {C^{ - 1}}{({C^{ - 1}} + {A^{ - 1}})^{ - 1}}{A^{ - 1}}

证明思路是将右边第二项的A^{-1}凑成-C^{-1}+C^{-1}+A^{-1},两边展开化简,即可得到。

推论

● Sherman-Morrison 定理(秩1校正定理)
A \in {R^n}可逆,向量u,v \in {R^n},若1 + {v^T}{A^{ - 1}}u \ne 0,则秩一校正矩阵为A + u{v^T}可逆,其逆矩阵为:
                {(A + u{v^T})^{ - 1}} = {A^{ - 1}} - \frac{{{A^{ - 1}}u{v^T}{A^{ - 1}}}}{{1 + {v^T}{A^{ - 1}}u}}

关于这个定理的推广,可见这两篇文章:

秩一校正定理 - 知乎 (zhihu.com)

Sherman-Morrison公式及其应用_山阴少年-CSDN博客

● 拟牛顿法,见这位博主写的文章:

Sherman-Morrison-Woodburg 定理_Liang_Ling的博客-CSDN博客

©著作权归作者所有,转载或内容合作请联系作者
  • 序言:七十年代末,一起剥皮案震惊了整个滨河市,随后出现的几起案子,更是在滨河造成了极大的恐慌,老刑警刘岩,带你破解...
    沈念sama阅读 205,033评论 6 478
  • 序言:滨河连续发生了三起死亡事件,死亡现场离奇诡异,居然都是意外死亡,警方通过查阅死者的电脑和手机,发现死者居然都...
    沈念sama阅读 87,725评论 2 381
  • 文/潘晓璐 我一进店门,熙熙楼的掌柜王于贵愁眉苦脸地迎上来,“玉大人,你说我怎么就摊上这事。” “怎么了?”我有些...
    开封第一讲书人阅读 151,473评论 0 338
  • 文/不坏的土叔 我叫张陵,是天一观的道长。 经常有香客问我,道长,这世上最难降的妖魔是什么? 我笑而不...
    开封第一讲书人阅读 54,846评论 1 277
  • 正文 为了忘掉前任,我火速办了婚礼,结果婚礼上,老公的妹妹穿的比我还像新娘。我一直安慰自己,他们只是感情好,可当我...
    茶点故事阅读 63,848评论 5 368
  • 文/花漫 我一把揭开白布。 她就那样静静地躺着,像睡着了一般。 火红的嫁衣衬着肌肤如雪。 梳的纹丝不乱的头发上,一...
    开封第一讲书人阅读 48,691评论 1 282
  • 那天,我揣着相机与录音,去河边找鬼。 笑死,一个胖子当着我的面吹牛,可吹牛的内容都是我干的。 我是一名探鬼主播,决...
    沈念sama阅读 38,053评论 3 399
  • 文/苍兰香墨 我猛地睁开眼,长吁一口气:“原来是场噩梦啊……” “哼!你这毒妇竟也来了?” 一声冷哼从身侧响起,我...
    开封第一讲书人阅读 36,700评论 0 258
  • 序言:老挝万荣一对情侣失踪,失踪者是张志新(化名)和其女友刘颖,没想到半个月后,有当地人在树林里发现了一具尸体,经...
    沈念sama阅读 42,856评论 1 300
  • 正文 独居荒郊野岭守林人离奇死亡,尸身上长有42处带血的脓包…… 初始之章·张勋 以下内容为张勋视角 年9月15日...
    茶点故事阅读 35,676评论 2 323
  • 正文 我和宋清朗相恋三年,在试婚纱的时候发现自己被绿了。 大学时的朋友给我发了我未婚夫和他白月光在一起吃饭的照片。...
    茶点故事阅读 37,787评论 1 333
  • 序言:一个原本活蹦乱跳的男人离奇死亡,死状恐怖,灵堂内的尸体忽然破棺而出,到底是诈尸还是另有隐情,我是刑警宁泽,带...
    沈念sama阅读 33,430评论 4 321
  • 正文 年R本政府宣布,位于F岛的核电站,受9级特大地震影响,放射性物质发生泄漏。R本人自食恶果不足惜,却给世界环境...
    茶点故事阅读 39,034评论 3 307
  • 文/蒙蒙 一、第九天 我趴在偏房一处隐蔽的房顶上张望。 院中可真热闹,春花似锦、人声如沸。这庄子的主人今日做“春日...
    开封第一讲书人阅读 29,990评论 0 19
  • 文/苍兰香墨 我抬头看了看天上的太阳。三九已至,却和暖如春,着一层夹袄步出监牢的瞬间,已是汗流浃背。 一阵脚步声响...
    开封第一讲书人阅读 31,218评论 1 260
  • 我被黑心中介骗来泰国打工, 没想到刚下飞机就差点儿被人妖公主榨干…… 1. 我叫王不留,地道东北人。 一个月前我还...
    沈念sama阅读 45,174评论 2 352
  • 正文 我出身青楼,却偏偏与公主长得像,于是被迫代替她去往敌国和亲。 传闻我的和亲对象是个残疾皇子,可洞房花烛夜当晚...
    茶点故事阅读 42,526评论 2 343

推荐阅读更多精彩内容