关于多元正态分布的条件概率密度

关于多元正态分布的条件概率密度原文来自师兄的博客 http blog csdn net wjj article details 多元正态分布多元正态分布的条件密度多元正态分布多元正态分布的密度函数如下 fx x1 xn 1 2 k 1 2exp 12 x T 1 x f x x 1 x n frac 1 sqrt 2 pi

原文来师兄的博客:http://blog.csdn.net/wjj/article/details/

多元正态分布



fx(x1,...xn)=1(2π)k|Σ|1/2exp(12(xμ)TΣ1(xμ))
(1)


其对应的矩母函数(也有称动差函数)为 exp(μTt+12tTΣt) 。事实上,如果随机向量 [X1,...Xn] 满足上面的动差函数,那么我们就称随机向量 [X1,...Xn] 服从多元高斯分布。具体地证明可以看 这里。

多元正态分布的条件密度

令随机向量 [X1,...Xn] 服从多元高斯分布。我们可以推导 Xn 在给定 X1,...Xn1 的情况下的条件密度分布:



f(xn|x1,...,xn1)=f(x1,...,xn1,xn)f(x1,...,xn1)
(2),

其中


f(x1,...,xn)=(2π)n/2(|Σ|1/2)exp[12ni,j=1yiqijyj]
(3)

其中 Q=Σ1=[qij],yi=xiμi 。同样地,



f(x1,...,xn1)=f(x1,...,xn1,xn)dxn=B(y1,...,yn1)
(4).

现在我们将公式(3)中的求和项进行分解,有:



ni,j=1yiqijyj=n1i,j=1yiqijyj+ynn1j=1qnjyj+ynn1i=1qinyj+qnny2n
(5)

因此,最终地条件分布具有如下的形式:



A(y1,...,yn1)B(y1,...,yn1)exp[(Cy2n+D(y1,...,yn1)yn)]
(6)

其中 C=(1/2)qnn ,因为 Q=Σ1 是对称矩阵,所以 D=n1j=1qnjyj=n1i=1qinyi .(6)式又可以进一步表示称如下的式子:



[ABexp(DD24C)]exp[(yn+D2C)2]1C
(7)

从公式(7)很容看出 xn 的条件密度函数是服从正态分布的。
所以条件分布的方差为: 2Var(Xn|X1,...,Xn1)=1/C ,进一步有: Var(Xn|X1,...,Xn1)=12C=1qnn
均值为:



E(Xn|X1,...,Xn1)=μnD2C=μn1qnnn1j=1qnj(Xjμj)


这就说明了再抽样多元正态分布时,如果已知了其它维度的随机变量值, 剩下的那个维度的随机变量也是服从正态分布
















版权声明:本文内容由互联网用户自发贡献,该文观点仅代表作者本人。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如发现本站有涉嫌侵权/违法违规的内容, 请联系我们举报,一经查实,本站将立刻删除。

发布者:全栈程序员-站长,转载请注明出处:https://javaforall.net/232359.html原文链接:https://javaforall.net

(0)
全栈程序员-站长的头像全栈程序员-站长


相关推荐

  • 电子设备日常使用总结

    电子设备日常使用总结

    2021年5月17日
    122
  • mysql磁盘阵列部署_部署磁盘阵列

    mysql磁盘阵列部署_部署磁盘阵列1、mdadm命令用于管理Linux系统中的软件RAID硬盘阵列,格式为mdadm[模式][选项][成员设备名称]常用命令:mdadm-D/dev/md0#md0为自定义设备名,查看详细信息2、mdadm命令的常用参数和作用-a  检测设备名称-n  指定设备数量-l  指定RAID级别-C  创建(阵列卡)-v  显示过程-f   模拟设备损坏-r   …

    2022年6月12日
    40
  • ZTE E700 自用感受及软件推荐

    ZTE E700 自用感受及软件推荐ZTEE700自用感受及软件推荐声明声明:此帖是转的  前两天看到移动08年新出的承诺话费换手机业务,换了个折扣最高的ZTEE700,首先声明我用的手机并不是很多,所以说错了希望大家见谅。这个机器我非常满意,非常好用!质量没得说,我从手里掉到水泥地上机器翻了3个跟头拿起来一点问题都没有。内置功能很强大,不过扩展功能则比较残废,好在机器本身的功能就能满足90%以上人的需求了,游戏…

    2022年7月11日
    12
  • java线程通信的三种方式「建议收藏」

    java线程通信的三种方式「建议收藏」1、传统的线程通信。在synchronized修饰的同步方法或者修饰的同步代码块中使用Object类提供的wait(),notify()和notifyAll()3个方法进行线程通信。关于这3个方法的解释:wait():导致当前线程等待,直到其他线程调用该同步监视器的notify()方法或notifyAll()方法来唤醒该线程。notify():唤醒在此…

    2022年6月19日
    25
  • 基于Redis6.2.6版本部署Redis Cluster集群

    基于Redis6.2.6版本部署Redis Cluster集群基于 Redis6 2 6 版本部署 RedisCluster 集群 1 Redis6 2 6 简介以及环境规划在 Redis6 x 版本中主要增加了多线程的新特性 多线性对于高并发场景是非常有必要的 Redis6 x 新特性如下 多线程 IO 重新设计了客户端缓存功能 RESP3 协议支持 SSLACL 权限控制提升了 RDB 日志加载速度发布官方的 Redis 集群代理模块 RedisCluster 集群原理可以查看之前发布的文章 环境规划 IP 主机名端口号

    2025年10月3日
    2
  • -Bash: Unzip: Command Not Found解决方法 安装unzip

    -Bash: Unzip: Command Not Found解决方法 安装unzip

    2021年10月18日
    68

发表回复

您的邮箱地址不会被公开。 必填项已用 * 标注

关注全栈程序员社区公众号