Sherman-Morrison-Woodbury formula 证明

文章目录

  • 1. 公式
  • 2. 证明

1. 公式

M = I − u v T ⇒ M − 1 = I + u v T 1 − v T u (1) M=I-uv^T\Rightarrow M^{-1}=I+\frac{uv^T}{1-v^Tu}\tag{1} M=IuvTM1=I+1vTuuvT(1)

2. 证明

  • 定义矩阵E表示如下:
    E = [ I u v T 1 ] , D = 1 − v T u (2) E=\begin{bmatrix}I&&u\\\\v^T&&1\end{bmatrix},D=1-v^Tu\tag{2} E= IvTu1 ,D=1vTu(2)
    我们想求 E − 1 E^{-1} E1,可以通过增广矩阵,进行行变换得到,
  • 第一种方法是:将第一行乘以 v T v^T vT后加到第二行中.
    [ I 0 − v T 1 ] E = [ I u 0 D ] (3) \begin{bmatrix}I&&0\\\\-v^T&&1\end{bmatrix}E=\begin{bmatrix}I&&u\\\\0&&D\end{bmatrix}\tag{3} IvT01 E= I0uD (3)
    E − 1 = [ I u 0 D ] − 1 [ I 0 − v T 1 ] = [ I + u D − 1 v T − u D − 1 − D − 1 v T D − 1 ] (4) E^{-1}=\begin{bmatrix}I&&u\\\\0&&D\end{bmatrix}^{-1}\begin{bmatrix}I&&0\\\\-v^T&&1\end{bmatrix}=\begin{bmatrix}I+uD^{-1}v^T&&-uD^{-1}\\\\-D^{-1}v^T&&D^{-1}\end{bmatrix}\tag{4} E1= I0uD 1 IvT01 = I+uD1vTD1vTuD1D1 (4)
  • 第二种方法是:将第二行乘以 u u u后加到第一行中.
    [ I − u 0 1 ] E = [ I − u v T 0 v T 1 ] (5) \begin{bmatrix}I&&-u\\\\0&&1\end{bmatrix}E=\begin{bmatrix}I-uv^T&&0\\\\v^T&&1\end{bmatrix}\tag{5} I0u1 E= IuvTvT01 (5)
    E − 1 = [ I − u v T 0 v T 1 ] − 1 [ I − u 0 1 ] = [ M − 1 − M − 1 u − v T M − 1 1 + v T M − 1 u ] (6) E^{-1}=\begin{bmatrix}I-uv^T&&0\\\\v^T&&1\end{bmatrix}^{-1}\begin{bmatrix}I&&-u\\\\0&&1\end{bmatrix}=\begin{bmatrix}M^{-1}&&-M^{-1}u\\\\-v^TM^{-1}&&1+v^TM^{-1}u\end{bmatrix}\tag{6} E1= IuvTvT01 1 I0u1 = M1vTM1M1u1+vTM1u (6)
  • 由公式4,6 可得,两个 E − 1 E^{-1} E1相等, M = I − u v T M=I-uv^T M=IuvT
    [ I + u D − 1 v T − u D − 1 − D − 1 v T D − 1 ] = [ M − 1 − M − 1 u − v T M − 1 1 + v T M − 1 u ] (7) \begin{bmatrix}I+uD^{-1}v^T&&-uD^{-1}\\\\-D^{-1}v^T&&D^{-1}\end{bmatrix}=\begin{bmatrix}M^{-1}&&-M^{-1}u\\\\-v^TM^{-1}&&1+v^TM^{-1}u\end{bmatrix}\tag{7} I+uD1vTD1vTuD1D1 = M1vTM1M1u1+vTM1u (7)
  • 由第一个行,第一列可得如下
    M − 1 = I + u D − 1 v T = I + u v T 1 − v T u (8) M^{-1}=I+uD^{-1}v^T=I+\frac{uv^T}{1-v^Tu}\tag{8} M1=I+uD1vT=I+1vTuuvT(8)
  • 结论:
    M = I − u v T ⇒ M − 1 = I + u v T 1 − v T u (9) M=I-uv^T\Rightarrow M^{-1}=I+\frac{uv^T}{1-v^Tu}\tag{9} M=IuvTM1=I+1vTuuvT(9)

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/news/833859.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

js手动实现unshift

js 手动实现数组的unshift unshift是什么? unshift() 方法可向数组的开头添加一个或更多元素,并返回新的长度。 注意: 该方法将改变数组的数目。 语法: array.unshift(item1,item2, ..., itemX)代码实现 首先,在…

【2】STM32·FreeRTOS·任务创建和删除

目录 一、任务创建和删除的API函数 1.1、动态创建任务函数 1.2、静态创建任务函数 1.3、任务删除函数 二、任务创建和删除(动态方法) 三、任务创建和删除(静态方法) 一、任务创建和删除的API函数 任务的创建和删除本质就是…

Flutter笔记:手动配置VSCode中Dart代码自动格式化

Flutter笔记 手动配置VSCode中Dart代码自动格式化 - 文章信息 - Author: 李俊才 (jcLee95) Visit me at CSDN: https://jclee95.blog.csdn.netMy WebSite:http://thispage.tech/Email: 291148484163.com. Shenzhen ChinaAddress of this article:https://blog.csd…

数据结构学习——线性表、顺序表

1.线性表 线性表 ( linear list ) 是n个具有相同特性的数据元素的有限序列。 线性表是一种在实际中广泛使 用的数据结构,常见的线性表:顺序表、链表、栈、队列、字符串… 线性表在逻辑上是线性结构,也就说是连续的一…

Jpa自定义查询结果封装到实体

工具类 import cn.hutool.core.convert.Convert; import cn.hutool.core.text.CharSequenceUtil;import java.lang.reflect.Field; import java.lang.reflect.Modifier; import java.util.ArrayList; import java.util.List; import java.util.Map;/*** 查询结果集转换工具类 …

在WPS表格(Excel)中,每10行增加一个特定的值

注:如下为WPS表格操作演示 例如1-15的数值是1,16-30就变为2,31-45就变为3,类推! 1、在B1单元格输入一个起始值,B2单元格输入公式IF(MOD(ROW(),15)0,B11,B1) 然后鼠标放到B2单元格右下角小点处&…

利用生成式AI重新构想ITSM的未来

对注入 AI 的生成式 ITSM 的需求,在 2023 年 Gartner AI 炒作周期中,生成式 AI 达到预期值达到顶峰后,三分之二的企业已经将生成式 AI 集成到其流程中。 你问为什么这种追求?在预定义算法的驱动下,IT 服务交付和管理中…

trivy使用方法

trivy使用方法 1、将镜像tar上传至服务器。 2、在tar包目录下,运行 docker load -i XXX.tar 3、docker images 查看镜像是否成功上传 4、trivy image 仓库名 --timeout 12h(每天首次扫描都会更新 无法跳过 耐心等待 后续扫描可不加timeout参数&#xff…

[HBCPC2023] Sakura(笛卡尔树)

Given A 1 , A 2 , ⋯ , A n A_1,A_2,⋯,A_n A1​,A2​,⋯,An​, please count the number of valid pairs of ( l , r l,r l,r) where l ≤ r l≤r l≤r and A l A r m a x i l r A i A_lA_rmax_{il}^rA_i Al​Ar​maxilr​Ai​. Input format: The first line contai…

C++学习第二十七课:STL中的位标志(Bitset)使用指南

C学习第二十七课:STL中的位标志(Bitset)使用指南 在C标准模板库(STL)中,std::bitset是一个固定大小的位集合,它提供了一种紧凑且方便的方式来存储和操作二进制位。本课将详细介绍std::bitset的…

代码随想录算法训练营第三十八天|动态规划理论基础,509. 斐波那契数,70. 爬楼梯,746. 使用最小花费爬楼梯

目录 动态规划理论基础509. 斐波那契数思路代码 70. 爬楼梯思路代码 746. 使用最小花费爬楼梯思路代码 动态规划理论基础 文档讲解:代码随想录 视频讲解:从此再也不怕动态规划了,动态规划解题方法论大曝光 !| 理论基础 |力扣刷题总…

React 学习-2

1.React State(状态) 每当 Clock 组件第一次加载到 DOM 中的时候,我们都想生成定时器,这在 React 中被称为挂载。 同样,每当 Clock 生成的这个 DOM 被移除的时候,我们也会想要清除定时器,这在 React 中被称为卸载。 …

对NI系统和PLC系统的应用比较

以下是对这两种系统的基本比较: 1. 设计和功能性 NI系统: 通常基于LabVIEW等软件平台,提供强大的数据采集、信号处理和图形界面开发能力。高度模块化和可扩展,支持各种传感器和信号类型。适合进行复杂的数据分析和高级控制算法的…

第七届机电、机器人与自动化国际会议(ICMRA 2024)即将召开!

第七届机电、机器人与自动化国际会议(ICMRA 2024)将于2024年9月20日-22日在中国武汉举行。ICMRA 2024为各国专家学者提供一个学术交流的平台,讨论机电、机器人和自动化领域的最新研究成果和未来的研究方向,旨在能够建立起国家间&a…

Python 基础知识:入门指南

Python 是一种简单易学、功能强大的编程语言,适用于各种用途,从简单的脚本编写到大型应用程序开发。如果你是初学者,以下是一份 Python 基础知识的入门指南,帮助你开始学习这门语言。 1. 安装 Python 首先,你需要在你…

Ansible剧本playbook之--------Templates 模块、roles角色详细解读

目录 一、Templates 模块 1.1准备模板文件并设置引用的变量 1.2修改主机清单文件,使用主机变量定义一个变量名相同,而值不同的变量 1.3编写 playbook 1.4ansible主机远程查看修改参数 1.5验证 二、tags 模块 always应用 三、Roles 模块 3.1ro…

在国企分公司做信息宣传新闻投稿的经验分享

作为一名国企分公司的信息宣传工作者,我亲历了从传统投稿方式到数字化转型的全过程,这段经历既充满了挑战,也收获了成长。回首最初的日子,那些用邮箱投稿的时光,至今仍让我感慨万千。 初尝辛酸,邮箱投稿的艰难岁月 刚接手信息宣传工作时,我满腔热情,却很快被现实的冷水浇了个透…

包管理器——apt篇

先给出两个官方文档。 1. PackageManagement - Debian Wiki 2.sourcelist 的编写规范 第 6 章 维护和更新:APT 工具 - 6.1. 写入sources.list文件 - 《Debian 8 管理员手册(Debian Jessie 从入门到精通)》 - 书栈网 BookStack 在 Debian 和…

RabbitMQ的五种模式

一、简单模式 简单模式(Simple):一个生产者,一个消费者 package com.qiangesoft.rabbitmq.mode.simple;import lombok.extern.slf4j.Slf4j; import org.springframework.amqp.rabbit.annotation.Queue; import org.springframe…

事业单位向媒体投稿发文章上级领导交给了我投稿方法

作为一名事业单位的普通职员,负责信息宣传工作,我见证了从传统投稿方式到智能化转型的全过程,这段旅程既是一次挑战,也是一次宝贵的成长。回想起初涉此领域的日子,那些通过邮箱投稿的时光,至今仍然历历在目,其中的酸甜苦辣,构成了我职业生涯中一段难忘的经历。 邮箱投稿:费时费…