变量消元(Varible Elimination)和概率边缘化(Marginalization)的关系
先定义消元:在解线性方程组的时候,把一个变量带入到另外一个变量中,达到减少变量的就结果。
- 虽然方程数少了,但是单个方程变复杂了,所以其实方程组携带的信息并没有减少。
- 如果把方程组写成矩阵形式,就对应之前线性代数学的高斯消元。
- 消元的一个作用是把方程变成上三角形式,就可以很轻松的计算出方程组的解。能计算出解,那么也能求出对应矩阵的逆。这就是为什么消元很重要的原因。
再定义边缘化:边缘化的普通定义我就不多说了,就是求积分那个。
- 当我们是边缘化一个多元高斯分布的随机变量的时候,如果我们知道这个随机变量的均值和协方差矩阵()。那么边缘化就是直接把中对应于要保留的分量拿出来就行。
- 假设a是要保留的,b是要被边缘化的。那么边缘化后不带b的概率分布就是只保留中和a有关的那一部分。
多元高斯分布的最大后验值(MAP)
- 我们想知道x取什么的时候,这个函数的值最大。很显然就是x=的时候
- 高斯函数还有另外一种表示方法
- 其实就是把exp里面的展开就能得到这个表达:,
- 这里面叫做信息向量,就是传说中的信息矩阵了
- 如果我们知道的表达形式是这种形式,就不能直接通过读均值知道MAP的结果了。
- 这种表示情况下MAP的结果是。也就是我们需要求信息矩阵的逆才能得到想要的东西。矩阵求逆等价于求方程,所以就和消元扯上关系了。
高斯分布条件下消元和边缘化的关系
- 基于上面的分析,边缘化是求中只和某些变量相关block。也就是求中的元素和中要求的block的关系。和是逆的关系,所以中的每一个元素都和中的每一个元素有关。既然是求逆,也就和消元搭上关系了。
- 既然和中的所有元素都有关,所以不是简单的丢弃性息。
和状态求解的关系
- 如果我们之表达相互独立的每次观察的结果的时候,可以直接用第一种表达写出观察两的分布。并且我们还知道我们真正关心的值和观察量之间有一定关系。也就是基于现实所知的信息能写出的一个表达通常是。
- x是一个随即变量,是我们想要求得的分布。A和b是已知的数据。Ax-b的到的随机变量分布也是已知的,是均值为0,协防差是的高斯分布。
- Mahalanobis距离是可以转化为L2距离:->,所以最终高斯函数的表示变为。多个这样的表示连乘求对数,对应为多个相加。
- Bx-d的结果是一个向量,就是向量点乘的直。
- 的形式为多个项的平方加在一起。B有多少行,这个平方和就有多少项。多个加一起,就是把这些平方全部加一起。所以可以组成一个大矩阵。这个矩阵的列数等于x的维度。行数为所有B的行数的和。
- 只有一个因子内部会有交叉项,因子之间不会出现交叉项。这就是为什么分布图里面对角上面会很密集。
- 我们知道当=0的时候,对应概率密度函数最大。所以就是要求解这个方程了。
- 因为B不是方阵,所以不能简单的用B的逆来求解,而要用广义逆。
- 在注意这里不能直接和上面的信息矩阵对应
和矩阵分解的关系
- 在求逆的时候,会使用矩阵分解的方法。
- 而矩阵的顺序对矩阵分解的计算量影响很大,比如slam中把和3d点有关的误差项排列在后面一起可以打打减少计算量。所以怎么排列矩阵就是一个大学问,isam中通过对图的变换,自动的实现最优的排列。并且当有新的变量加入的时候,通过对图的分析可以知道,那些因子是受新变量影响的哪些不是。从而实现部分更新。
变量消元(Varible Elimination)和概率边缘化(Marginalization)的关系相关推荐
- 【人工智能】— 贝叶斯网络、概率图模型、全局语义、因果链、朴素贝叶斯模型、枚举推理、变量消元
[人工智能]- 贝叶斯网络 频率学派 vs. 贝叶斯学派 贝叶斯学派 Probability(概率): 独立性/条件独立性: Probability Theory(概率论): Graphical mo ...
- 贝叶斯网专题4:概率推理中的变量消元方法
文章目录 第一部分:贝叶斯网基础 [1.1 信息论基础](https://blog.csdn.net/deepbodhi/article/details/119823055) [1.2 贝叶斯网基本概 ...
- 代数几何:消元(Elimination)、扩展(Extension)
参考文献:Ideals, Varieties, and Algorithms (4th ed.) [Cox, Little & O'Shea 2015-06-14] 前置文章:仿射簇 和 Gr ...
- Matrix 高斯消元Gaussian elimination 中的complete pivoting和partial pivoting
首先科普下Pivoting的含义 一般翻译为"主元",在对矩阵做某种算法时,首先进行的部分元素.在线性规划的单纯形法中常见. wiki的解释如下: Pivot element (t ...
- ACM数学模板1 高斯消元 Gaussian Elimination
声明: 题目来源: https://www.luogu.com.cn/problem/P3389 https://www.luogu.com.cn/problem/P2447 https://ac.n ...
- 论文笔记——基于因子图消元优化的多传感器融合定位算法
创新点: 为了提高抗干扰性和容错能力,在链式因子图模型中加入滑动窗口用于保留窗口内历史状态信息:同时为了避免高维矩阵运算,引入消元算法将因子图转化为贝叶斯网络,依次边缘化历史状态,实现矩阵降维. 提出 ...
- BZOJ 3270: 博物馆 [概率DP 高斯消元]
http://www.lydsy.com/JudgeOnline/problem.php?id=3270 题意:一张无向图,一开始两人分别在$x$和$y$,每一分钟在点$i$不走的概率为$p[i]$, ...
- BZOJ:4820: [Sdoi2017]硬币游戏BZOJ:1444: [Jsoi2009]有趣的游戏(高斯消元求概率)
1444: [Jsoi2009]有趣的游戏 4820: [Sdoi2017]硬币游戏 这两道题都是关于不断随机生成字符后求出现给定字符串的概率的问题. 第一题数据范围较小,将串建成AC自动机以后,以A ...
- BZOJ 3270: 博物馆 1778: 驱逐猪猡 【概率DP+高斯消元】
题目描述: 中文题面,不多解释.1778传送门 3270 传送门 (博物馆)题目分析: 也许很多人做概率题的时候都有种虚幻感..感觉莫名其妙就得出一个期望.概率,一知半解... 所以我在这里仔细地剖析 ...
- [高斯消元 概率 KMP] BZOJ 4820 [Sdoi2017]硬币游戏
一个直观的想法 是建AC自动机 然后消元 但是这样变量个数是O(nm)O(nm) 然后我就不会做了 概率题都好妙啊 一个精妙的设计是再定义一个状态N 表示当前串不包含任何人的概率 举个例子 来自这里 ...
最新文章
- 设计模式之简单工厂模式(Simply Factory)摘录
- 面试官:Java反射是什么?我回答不上来!
- C# ICSharpCode.SharpZipLib.Zip 的使用
- 简单几步教你去除开机出现扫描硬盘!
- HALCON示例程序train_characters_ocr.hdev使用SVM分类器训练字体
- java开发安装mysql_从零开始搭建Java开发环境第二篇:如何在windows10里安装MySQL
- 第十一章 图形视图、动画、状态机框架
- SQL:select case when(转)
- 弱鸡儿长乐爆肝旅Day8
- Pandas数据帧(DataFrame)
- RHEL5.6环境下yum安装MySQL
- linux gpio-led
- sql top加不加括号的区别_SQL易错点大作战
- CAD制图中该怎么查找文字呢?
- Lambda将搬砖变成艺术
- js 数组转json,json转数组
- STM32-IIC模拟从模式
- 长腔耗散孤子锁模光纤激光器
- 电商产品设计:促销活动设计解析
- python刷火车票_五一要来了,教你用Python动刷新抢12306火车票,附源码
热门文章
- Java学习实战教程之mysql优化
- Chrome 无法正常打开站点
- cmd配网关_win7下用DOS命令修改IP、网关和DNS
- ckfinder php,php,ckfinder_ckfinder mac下报错 linux正常,php,ckfinder - phpStudy
- 自主导航与路径规划无人机研究现状
- CloneZilla再生龙 | 打造自己的深度学习镜像
- 用户帐户控制组策略和注册表项设置
- 万年历黄历星座查询v3.6.7
- 装两个独立的双系统,在笔记本电脑上
- mysql slave_pending_jobs_size_max_3分钟解决MySQL 1864 主从错误