矩阵求逆引理:
(Z+UWV)−1=Z−1−Z−1U(W−1+VZ−1U)−1VZ−1(Z+UWV)^{-1} = Z^{-1} - Z^{-1}U(W^{-1}+V Z^{-1}U)^{-1}V Z^{-1} (Z+UWV)−1=Z−1−Z−1U(W−1+VZ−1U)−1VZ−1


你可能知道这个定理是用来加速计算的,但你知道怎么用这个定理吗?什么时候可以用?

下面来解释一下这个定理的应用场景:

Z∈Rn×nZ\in R^{n\times n}Z∈Rn×n, 我们需要对ZZZ求逆,计算Z−1Z^{-1}Z−1. 计算复杂度为 O(n3)O(n^3)O(n3).

假设我们计算好了 Z−1Z^{-1}Z−1,然后发现 ZZZ 受到了扰动变成了 Z^=Z+E\hat{Z}= Z+EZ^=Z+E,需要做什么修正呢?

如果直接计算 Z^−1\hat{Z}^{-1}Z^−1,再花费O(n3)O(n^3)O(n3)的计算量,未免有点浪费,因为之前的计算结果Z−1Z^{-1}Z−1完全没有用到,有什么办法用到已有的计算结果Z−1Z^{-1}Z−1来加速求解Z^−1\hat{Z}^{-1}Z^−1呢?

这就用到了逆矩阵引理!

假设扰动是低秩矩阵,即可分解成
E=UWVE=UWV E=UWV其中 U∈Rn×m,W∈Rm×m,V∈Rm×n,m<<nU\in R^{n\times m}, W \in R^{m \times m}, V\in R^{m \times n}, m<< nU∈Rn×m,W∈Rm×m,V∈Rm×n,m<<n,且 WWW 可逆。例如奇异值分解。

接下来利用矩阵求逆引理,就可快速求解 Z^−1\hat{Z}^{-1}Z^−1 了,看看等式右端的计算复杂度,只需要计算 mmm 阶矩阵的逆和几个矩阵乘法即可,矩阵乘法是可以并行加速的,因此计算时间大大减少!

矩阵求逆引理(Matrix Inversion Lemma)的意义相关推荐

  1. 矩阵求逆引理(Matrix inversion lemma)推导

      矩阵求逆引理(Matrix inversion lemma):   现有矩阵AAA可以写为如下分块矩阵形式: A=[A11A21A12A22](m+n)×(m+n)A=[A11A12A21A22] ...

  2. 矩阵求逆引理(matrix inversion lemma)

    论文:Li J, Stoica P. An adaptive filtering approach to spectral estimation and SAR imaging[J]. IEEE Tr ...

  3. 矩阵行列式引理 Matrix Determinant Lemma

    在看论文时有用到矩阵行列式引理,搜了一下感觉很有趣,记录在这里. wiki:https://en.wikipedia.org/wiki/Matrix_determinant_lemma 数学描述:AA ...

  4. 再探矩阵求逆引理 : Woodbury恒等式的证明

    在之前的许多次接触中感到了矩阵求逆引理的强大可谓在通信中无处不在.然而稍显繁琐的表达式让我总是没能完整记忆,每次都要对着矩阵论的书才得以使用.痛定思痛,恰好今天看了维基上对其的介绍,觉得太过优雅,就以 ...

  5. 逆矩阵(inverse matrix)的概念及其意义

    逆矩阵(inverse matrix)的概念及其意义 2015年09月17日 00:09:10 阅读数:21838 标签: 逆矩阵为何需要逆矩阵逆矩阵应用逆矩阵实例逆矩阵与倒数 更多 版权声明:本文为 ...

  6. Monotonic Matrix(lemma定理)

    lemma定理: 有n对点,设表示到的路径方案数 则行列式表示到的严格不相交路径的方案数 #include<bits/stdc++.h> #define ll long long usin ...

  7. 20220312 矩阵求逆引理

    如果对任一 n×nn \times nn×n 阶非奇异矩阵 A\boldsymbol{A}A 与任意两个 n×mn \times mn×m 阶矩阵 B\mathbf{B}B 和 C\mathbf{C} ...

  8. Inversion Lemma

    来源:概率机器人

  9. 最小二乘估计 Least Squares estimation

    本文主要讲标准最小二乘方法及其常见的变形:加权最小二乘和总体最小二乘算法,关注不同方法之间的逻辑. 一.最小二乘估计(Least Squares estimation,LS) 最小二乘估计方法是一种不 ...

最新文章

  1. vue 函数 路由跳转_vue中通过路由跳转的三种方式
  2. 笔记-项目成本管理-估算成本的步骤
  3. F​P​G​A​工​作​原​理
  4. 互联网 DBA 需要做那些事(转)
  5. 牛客挑战赛53G-同源数组(Easy Version)【NTT】
  6. aix系统java堆_浅谈AIX环境下的Java性能调优
  7. html5 页面加载缓慢,html5体验优化页面加载的14条建议
  8. class属性表种类集合,以及字段详解
  9. 首届 “女生科技体验节” 大爆料!
  10. CentOS 7.0 关闭防火墙
  11. 在OpenCV中使用YOLO v3进行物体检测
  12. 《移动应用开发技术——Android》课程报告-个人记账系统
  13. python文件操作
  14. 第一讲 数学方法论引论
  15. 华为版计算机软件,Huawei华为手机PC客户端软件
  16. 三角函数中的和差化积公式编辑方法
  17. 为何硅谷第一性感女人也没能拯救雅虎?
  18. Excel中替换单元格中的换行符
  19. 使用Dir函数遍历文件和目录
  20. Java系列技术之Mybatis3-钟洪发-专题视频课程

热门文章

  1. UE4中的玩家类UPlayer、ULocalPlayer 和 UNetConnection
  2. 2.机器学习基础(一)
  3. “添翼杯”人工智能创新应用大赛垃圾分类之模型部分
  4. 重头系统的学习,不会咱就学!2014.6.18
  5. python樱花_python画樱花树代码 具体代码介绍
  6. npm run build 打包报错primordials is not defined的解决方法
  7. win 32学习笔记(三) 消息队列
  8. 微信整合CRM系统的好处
  9. 超新星与暗能量的发现--今年诺贝尔物理奖工作的介绍(转)
  10. 对于程序员来说养成规范的打字真的很重要