矩阵论记号约定​zhuanlan.zhihu.com

Schur补与矩阵打洞

考虑正态随机向量

,求
条件分布的一个类似于线性回归的常用技巧是——寻找常值矩阵
适合
;对于正态变量而言独立等价于不相关,因而只需
为零,即
;此时
,从而立得
。形式化上述过程如下:
  • 设分块矩阵

    的子矩阵
    可逆,将

    Schur补记为

  • [打洞技巧[1]]
  • 易得推论[Guttman]
    ;考虑
    可得

    Frobenius不等式

另证:注意到

是良定义的

同态。

  • 对方阵成立

    ;考虑
    可得
    ,进而有

    Sylvester行列式定理

  • 进一步假设
    是Hermite矩阵,特征值为

    惯性指数记为

    ,其中
  • [Haynsworth惯性加法公式] 注意

    可逆,惯性指数在合同变换下保持不变,即得
  • 特别地,
    正定

    当且仅当

    正定。

Sherman–Morrison–Woodbury求逆公式

可逆(
),
,考虑低秩扰动
后逆矩阵的更新

[2],有Woodbury公式

,这里假设了
可逆;当
时,
称为

Sherman–Morrison公式

证明:

打洞,分别化为
,求逆可得
  • [] 我们还得到了分块矩阵求逆公式

另证:可以利用巧妙的恒等式

记忆花招[3]

直接计算可得

  • [习题] 设

    均为正定矩阵, 则
    (半)正定

    当且仅当

    (半)正定.

分块半正定矩阵与行列式(不)等式

  • [Fischer不等式] 设

    是半正定矩阵,则

证明:

正定,则
正定,由Brunn-Minkowski不等式的推论可得
,而
。若
半正定,在
中令
即证。
  • [华罗庚] 设

    ,则
    ;若
    小于等于(严格小于)
    ,则称

    压缩(严格压缩)算子,对于严格压缩算子有

  • 由Brunn-Minkowski不等式的推论可得
    ;对于压缩算子有

证明:

,则
。应用前述分块矩阵求逆公式,对于非奇异情形有
,进而
得到
,即
;一般情形通过连续性论证易得。对于压缩算子,考虑
,同理可证。

参考

  1. ^孟道骥, & 王立云. (2006). 打洞技巧. 高等数学研究, 9(4): 15-19. http://qikan.cqvip.com/Qikan/Article/Detail?id=22377003
  2. ^W. W. Hager (1989). Updating the inverse of a matrix. SIAM review, 31(2): 221-239. https://doi.org/10.1137/1031049
  3. ^In a ring, 1−ab is invertible ⟹ 1−ba is invertible (from scratch) https://math.stackexchange.com/q/2341289/469791

php随机数字不重复使等式成立_Schur补与矩阵打洞,SMW求逆公式,分块矩阵与行列式(不)等式...相关推荐

  1. php随机数字不重复使等式成立_当随机数遇上量子

    导读 小墨读大学的时候,坊间流传着一段顺口溜,用来描述大学几门特别难学的课,其中有一句:"-随机过程随机过,量子力学量力学...",可见这两门课消耗了多少学子的脑细胞.那么,当&q ...

  2. php随机数字不重复使等式成立_php随机数不重复生成的五种方式

    在前一篇文章中我介绍了php随机数的生成方法,我们都知道php随机数不管是在开发中,还是在移动应用中,都是不可少的一部分,平时日常开发中经常会跟php随机数打交道,那么今天就给大家介绍php随机数不重 ...

  3. php随机数字不重复,php生成N个不重复的随机数

    php生成N个不重复的随机数 生成N个不重复的随机数,如何在php中实现呢?本文分享的这例php代码,可以实现随机数的'生成,生成多个不重复的随机数,有兴趣的朋友参考下. php生成N个不重复的随机数 ...

  4. 34.写一个排序算法1-100随机数字进行排序

    编程题目: 34.写一个排序算法 1-100随机数字 进行排序,要求效率. 示例代码: package program.calculation.exercise34;/*** 34.写一个排序算法 1 ...

  5. 用1、2、3、4、5、6、7、8、9这9个数字,填入□ 中使等式□□×□□□ = □□□□ 成立,每个数字恰好只用一次。

    用1.2.3.4.5.6.7.8.9这9个数字,填入□ 中使等式□□×□□□ = □□□□ 成立,每个数字恰好只用一次. 以下是所有的7 个答案供参考: 12*483=5796 18*297=5346 ...

  6. 无重复数字的随机数字数组

    这几天看MFC,头晕眼花,也觉得没有什么可以写的. 今天学习之余,看到一个面试题"1000个范围0-2000的数字,进行排序打印出来".就想着实现一下. 第一步先要建立1000个范 ...

  7. c语言----移动n根火柴使等式成立

    一.项目背景 近日在观看某短视频软件时,发现某些用户直播移动火柴游戏,笔者有时会很快找出方法,有时却毫无头绪,如5+7=9,怀疑无解,故以此编程. 二.编程环境 笔者用的是 Visual Studio ...

  8. 移动一根火柴使等式成立js版本(递归)

    修改成递归版本 思路: 1.设定规则数组,比如:1加一根火柴只可以变成7. 2.设定方法数组,比如:一个数增加了一根火柴,其他的数必然减少一根火柴. 3.增加Array方法,由元素名和方法,得到规则对 ...

  9. Mysql 生成不重复的随机数字

    在网上查找Mysql 生成不重复的随机数字 ,竟然没找到合适的例子. 其实思路很简单,利用MySQL现有的函数,然后进行加工处理,达到预期的结果.可以用到的MySQL函数为rand() ,以及 rou ...

最新文章

  1. python制作动图-用Python制作迷宫GIF
  2. 虚拟机中web服务器的搭建,虚拟机下搭建web服务器
  3. 【Linux】一步一步学Linux——kill命令(124)
  4. 电脑桌面图标异常 计算机只显示C盘,Windows7系统桌面图标异常如何解决?
  5. 1-2docker-基本的使用
  6. python混沌时间序列分析_用Python进行时间序列分析
  7. Linux系统学习: 进程、重定向和管道指令:xargs 指令的作用是
  8. 阶乘的传统流程图 c语言,C语言算法与流程图.ppt
  9. ipv6一致性测试(IPv6 Specification)内核修改
  10. ros + gazebo未报错却加载不出来机器人模型 解决办法
  11. 网页中人民币符号的正确输入姿势
  12. T1292:宠物小精灵之收服
  13. P35-前端基础-新闻列表
  14. 搭建多节点Linux环境
  15. 算法竞赛入门经典(刘汝佳)——代码笔记
  16. python写的百度贴吧邮箱采集(带界面)
  17. 怎么下载linux sio镜像,在linux下如何安裝ftdi_sio USB轉串口驅動
  18. 28推论坛邀请注册攻略
  19. 目标检测网络之三叉戟TridentNet
  20. flash高科技php网站源码下载,Flash小游戏PHP源码

热门文章

  1. HTML5 代码要怎样凭“魅力”吸引搜索引擎的注意?
  2. JavaScript 败北,TypeScript 大势所趋?
  3. 华为双 11 发 20 亿奖金!?
  4. 重磅!Google推出了Python最牛X的编辑器....
  5. 亲测!这款耳机性价比堪比 AirPods
  6. 干不过支付宝的微信,还有翻盘的可能吗?| 畅言
  7. iPhone 大降价;谷歌再爆丑闻;京东云金山云回应合并传闻 | 极客头条
  8. 前端开发者如何利用 CSS 实现酷炫的变色方案?
  9. 蚂蚁金服双 11 大促全面揭秘:百万支付、容器化和平台智能化
  10. 中兴正式解禁;小米市值与索尼相当;CMU教授炮轰谷歌 | 极客头条