置信传播算法(Belief Propagation)简介_app-bp概率域置信传播-程序员宅基地

技术标签: 算法  5G  python  通信系统  人工智能  

基础知识

条件概率(Conditional Probability)

在这里插入图片描述
在这里插入图片描述
相互独立时,p(A | B) = p(A)

贝叶斯规则

贝叶斯网络(Bayesian Network)定了一个独立的结构:一个节点的概率仅依赖于它的父节点。贝叶斯网络适用于稀疏模型,即大部分节点之间不存在任何直接的依赖关系。
在这里插入图片描述
在这里插入图片描述

联合概率(Joint Probability)

表示所有节点共同发生的概率,将所有条件概率相乘:
在这里插入图片描述
我们最终的目标是计算准确的边缘概率(Marginal Probability),比如计算Hangover的概率,边缘概率为各种状态下所有其他节点对本节点影响的概率的和。

边缘概率(Marginal Probability)

即某个事件发生的概率,而与其它事件无关。边缘概率是这样计算的:在联合概率中,把最终结果中不需要的那些事件合并成其事件的全概率而消失(在两个离散随机变量的条件下,对于其中任一行或任一列求和,得到的概率就是边缘概率)。在本例中,针对不同的Hangover进行求和,得到的就是Hangover的边缘概率:
在这里插入图片描述
如果贝叶斯网络比较小,我们可以很简单的做边缘求和运算,但是如果问题规模较大,整个运算复杂度和数据将会以指数级增长。而利用BP算法去计算这样的网络问题,可以使得运算复杂度只和节点数线性相关。在这种意义上,BP算法在大型贝叶斯网络推断问题中扮演着越来越重要的作用。

马尔科夫随机场(Markov Random Field,MRF)

在这里插入图片描述

在概率图模型中,每个结点表示一个随机变量,结点之间的边表示这些随机变量之间的概率关系。在概率图模型中,所有随机变量的联合概率分布可以表示成若干随机变量子集的乘积。典型的概率图模型包括贝叶斯网和马尔科夫网。贝叶斯网是有向图模型, 用于表示随机变量之间的因果关系,而马尔科夫网是无向图模型, 用于表示随机变量的概率分布和概率推理,或者说是随机变量之间的软约束关系。

BP算法的基础就是建立于MRF上,MRF是一种条件概率模型,它可以被认为是马尔科夫链的一种推广,其对于场内所有节点的相关性都能很有效的描述。

假设我们观察到yi的一些信息,需要利用这些已知信息去推断关于隐含的场景xi的另外一些信息。每个顶点i都有一个状态值xi和一个观测值yi,每个状态值和观测值之间的似然函数为Фi(xi,yi),反映了i处的 xi和 yi存在统计依赖性,表示节点i的联合相容度,相邻邻居节点之间的势能量为Ψij(xi,xj),Ψij(xi,xj)也称为相邻节点之间的不连续代价,反映了节点变量 xi 和 xj 之间的相容性,体现了随机场自身具备的约束条件。

置信度传播(Belief Propagation,BP)

置信度传播算法利用结点与结点之间相互传递信息而更新当前整个MRF的标记状态,是基于MRF的一种近似计算。该算法是一种迭代的方法,可以解决概率图模型概率推断问题,而且所有信息的传播可以并行实现。经过多次迭代后,所有结点的信度不再发生变化,就称此时每一个结点的标记即为最优标记,MRF也达到了收敛状态。对于无环环路的MRF,BP算法可以收敛到其最优解。

BP算法的两个关键过程:(1)通过加权乘积计算所有的局部消息;(2)节点之间概率消息在随机场中的传递。
在这里插入图片描述
有了消息更新规则以及置信度计算公式,就可以先任意初始化每个 b i ( x i ) b_i(x_i) bi(xi),然后迭代的求解 m i j m_{ij} mij b i ( x i ) b_i(x_i) bi(xi)直至收敛,也就是说 m i j m_{ij} mij不再发生变化。也就是说首先对一些初始节点的消息赋初值,然后多次迭代消息传播和置信度更新直到它们稳定,最后就能从置信度中获取相应的概率。

置信度替换为概率:
在这里插入图片描述
b i ( x i ) b_i(x_i) bi(xi)为节点 i i i的联合概率分布,其中 m j i ( x i ) m_{ji}(x_{i}) mji(xi)代表隐含节点 j j j传递给隐含节点i的消息,表明了隐含节点 i i i对隐含节点 j j j当前状态的影响。 Ф i ( x i , y i ) Ф_i(x_i,y_i) Фi(xi,yi) 表示节点 i i i的局部证据,表示节点 i i i的联合相容度。节点i的置信度 b i ( x i ) b_i(x_i) bi(xi) i i i的邻域向 i i i传递的所有消息的乘积成正比,同时也正比于 Ф i ( x i , y i ) Ф_i(x_i,y_i) Фi(xi,yi) 1 / z i 1/z_i 1/zi为归一化常数,可使置信度的和为1, N ( i ) N(i) N(i) 为节点i的MRF一阶邻域。

消息传播的信息:
在这里插入图片描述
包含所有其他传入节点i的消息乘积, N ( j ) / i N(j) / i N(j)/i表示节点j的MRF一阶邻域中排除掉目标节点i的邻域。

有了消息更新规则以及置信度计算公式,就可以先任意初始化每个 b i ( x i ) b_i(x_i) bi(xi),然后迭代的求解 m i j m_{ij} mij b i ( x i ) b_i(x_i) bi(xi)直至收敛,也就是说 m i j m_{ij} mij不再发生变化。也就是说首先对一些初始节点的消息赋初值,然后多次迭代消息传播和置信度更新直到它们稳定,最后就能从置信度中获取相应的概率。

置信度传播算法中迭代运算步骤可以表示如下:

(1)随机选择相邻的隐含节点 x i x_i xi, x j x_j xj

(2)从 x i x_i xi x j x_j xj发送消息 m i j m_{ij} mij
在这里插入图片描述
在这里插入图片描述
(3)更新节点 x j x_j xj的置信度
在这里插入图片描述
在这里插入图片描述

(4)跳至步骤(1),直到算法收敛

在以此为规则的计算中,从无环图的边缘节点开始传播,然后如果一个节点所有相邻节点的消息都已经计算出来,则计算该节点的消息。易得整个无环图仅仅只需计算一遍就可以得到所有隐含节点的边缘概率分布。可以看出,BP 算法相对于一般的算法,时间复杂度上是大幅下降的。

在这里插入图片描述

版权声明:本文为博主原创文章,遵循 CC 4.0 BY-SA 版权协议,转载请附上原文出处链接和本声明。
本文链接:https://blog.csdn.net/daijingxin/article/details/115800989

智能推荐

c# 调用c++ lib静态库_c#调用lib-程序员宅基地

文章浏览阅读2w次,点赞7次,收藏51次。四个步骤1.创建C++ Win32项目动态库dll 2.在Win32项目动态库中添加 外部依赖项 lib头文件和lib库3.导出C接口4.c#调用c++动态库开始你的表演...①创建一个空白的解决方案,在解决方案中添加 Visual C++ , Win32 项目空白解决方案的创建:添加Visual C++ , Win32 项目这......_c#调用lib

deepin/ubuntu安装苹方字体-程序员宅基地

文章浏览阅读4.6k次。苹方字体是苹果系统上的黑体,挺好看的。注重颜值的网站都会使用,例如知乎:font-family: -apple-system, BlinkMacSystemFont, Helvetica Neue, PingFang SC, Microsoft YaHei, Source Han Sans SC, Noto Sans CJK SC, W..._ubuntu pingfang

html表单常见操作汇总_html表单的处理程序有那些-程序员宅基地

文章浏览阅读159次。表单表单概述表单标签表单域按钮控件demo表单标签表单标签基本语法结构<form action="处理数据程序的url地址“ method=”get|post“ name="表单名称”></form><!--action,当提交表单时,向何处发送表单中的数据,地址可以是相对地址也可以是绝对地址--><!--method将表单中的数据传送给服务器处理,get方式直接显示在url地址中,数据可以被缓存,且长度有限制;而post方式数据隐藏传输,_html表单的处理程序有那些

PHP设置谷歌验证器(Google Authenticator)实现操作二步验证_php otp 验证器-程序员宅基地

文章浏览阅读1.2k次。使用说明:开启Google的登陆二步验证(即Google Authenticator服务)后用户登陆时需要输入额外由手机客户端生成的一次性密码。实现Google Authenticator功能需要服务器端和客户端的支持。服务器端负责密钥的生成、验证一次性密码是否正确。客户端记录密钥后生成一次性密码。下载谷歌验证类库文件放到项目合适位置(我这边放在项目Vender下面)https://github.com/PHPGangsta/GoogleAuthenticatorPHP代码示例://引入谷_php otp 验证器

【Python】matplotlib.plot画图横坐标混乱及间隔处理_matplotlib更改横轴间距-程序员宅基地

文章浏览阅读4.3k次,点赞5次,收藏11次。matplotlib.plot画图横坐标混乱及间隔处理_matplotlib更改横轴间距

docker — 容器存储_docker 保存容器-程序员宅基地

文章浏览阅读2.2k次。①Storage driver 处理各镜像层及容器层的处理细节,实现了多层数据的堆叠,为用户 提供了多层数据合并后的统一视图②所有 Storage driver 都使用可堆叠图像层和写时复制(CoW)策略③docker info 命令可查看当系统上的 storage driver主要用于测试目的,不建议用于生成环境。_docker 保存容器

随便推点

网络拓扑结构_网络拓扑csdn-程序员宅基地

文章浏览阅读834次,点赞27次,收藏13次。网络拓扑结构是指计算机网络中各组件(如计算机、服务器、打印机、路由器、交换机等设备)及其连接线路在物理布局或逻辑构型上的排列形式。这种布局不仅描述了设备间的实际物理连接方式,也决定了数据在网络中流动的路径和方式。不同的网络拓扑结构影响着网络的性能、可靠性、可扩展性及管理维护的难易程度。_网络拓扑csdn

JS重写Date函数,兼容IOS系统_date.prototype 将所有 ios-程序员宅基地

文章浏览阅读1.8k次,点赞5次,收藏8次。IOS系统Date的坑要创建一个指定时间的new Date对象时,通常的做法是:new Date("2020-09-21 11:11:00")这行代码在 PC 端和安卓端都是正常的,而在 iOS 端则会提示 Invalid Date 无效日期。在IOS年月日中间的横岗许换成斜杠,也就是new Date("2020/09/21 11:11:00")通常为了兼容IOS的这个坑,需要做一些额外的特殊处理,笔者在开发的时候经常会忘了兼容IOS系统。所以就想试着重写Date函数,一劳永逸,避免每次ne_date.prototype 将所有 ios

如何将EXCEL表导入plsql数据库中-程序员宅基地

文章浏览阅读5.3k次。方法一:用PLSQL Developer工具。 1 在PLSQL Developer的sql window里输入select * from test for update; 2 按F8执行 3 打开锁, 再按一下加号. 鼠标点到第一列的列头,使全列成选中状态,然后粘贴,最后commit提交即可。(前提..._excel导入pl/sql

Git常用命令速查手册-程序员宅基地

文章浏览阅读83次。Git常用命令速查手册1、初始化仓库git init2、将文件添加到仓库git add 文件名 # 将工作区的某个文件添加到暂存区 git add -u # 添加所有被tracked文件中被修改或删除的文件信息到暂存区,不处理untracked的文件git add -A # 添加所有被tracked文件中被修改或删除的文件信息到暂存区,包括untracked的文件...

分享119个ASP.NET源码总有一个是你想要的_千博二手车源码v2023 build 1120-程序员宅基地

文章浏览阅读202次。分享119个ASP.NET源码总有一个是你想要的_千博二手车源码v2023 build 1120

【C++缺省函数】 空类默认产生的6个类成员函数_空类默认产生哪些类成员函数-程序员宅基地

文章浏览阅读1.8k次。版权声明:转载请注明出处 http://blog.csdn.net/irean_lau。目录(?)[+]1、缺省构造函数。2、缺省拷贝构造函数。3、 缺省析构函数。4、缺省赋值运算符。5、缺省取址运算符。6、 缺省取址运算符 const。[cpp] view plain copy_空类默认产生哪些类成员函数

推荐文章

热门文章

相关标签