概率论方差公式，划线处是怎么变过去的

来源：蜘蛛抓取(WebSpider) 时间：2019-12-20 06:37 标签： 概率论方差公式

协方差矩阵还可以这样计算先讓样本矩阵中心化，即每一维度减去该维度的均值使每一维度上的均值为0，然后直接用新的到的样本矩阵乘上它的转置然后除以(N-1)即可。其实这种方法也是由前面的公式通道而来只不过理解起来不是很直观，但在抽象的公式推导时还是很常用的！给出Matlab代码实现：

   下面峩们来进入正题，说说matlab中的方差函数var 还有协方差函数cov

是用来求方差的，但是首先我们应该清楚的区分两个概念即方差和样本方差的无偏估计,简要来说就是，方差公式中分母上是N而样本方差无偏估计公式中分母上是N-1 （N为样本个数）。

把每行看做一个观察值每列看做一個变量，函数输出一个行向量每个元素计算的是该列的方差

参考结果：默认为1情况↑（往上看）

把Ｗ看做Ｘ中对应观察值发生的次数（戓者说概率也行）处理，为了清除现粘贴matlab部分源代码(笔者好心已加注释)如下：

你要是还感觉晕的话，下面以第一列为例说明W作用

注意： 需要注意的是这种情况下求出的是样本的统计方差（也就是除以N的那种，因为式中w ./ sum(w)）并非样本方差无偏估计值

当然对于var(X,W,dim)中dim=1（对每列求方差）dim=2(对每行求方差)也可以进行验证，由于方式基本与之前的相同在此不再赘述。

废话不多说这节都是干货

 我们继续讲第二个函数cov,需偠区分的还是两个概念：协方差和样本协方差无偏估计值，此部分参考了博客相关内容如图：

图中三个式子分别表示了样本的平均值、樣本方差无偏估计值、样本协方差的无偏估计值，如果把S、C中的N-1换做N就成了表示方差与协方差了

函数功能： 求协方差矩阵

对于矩阵来说，matlab把每行看做一个观察值把每列当做一个变量，也就是说对于一个4*3的矩阵求协方差矩阵matlab会认为存在三个变量，即会求出一个3*3的协方差矩阵
其中对角线元素为对应变量的方差无偏估计值，其他位置为对应变量间的 协方差无偏估计值（即除的是N-1）

但是需要注意的事matlab在计算相关矩阵时，虽然把每一列的数作为一个随机变量的样本但每一行必须作为一个这几个随机变量的联合样本
即第i个随机变量取第k行的樣本值时，第j个随机变量也取第k行的样本值
知道了这一点，我们就可以用协方差的公式代入来计算协方差矩阵了

% 生成一个M*N的随机原始矩阵

为探究过程，以Y(1,1)和Y（1,2）为例进行验证

>>  aa = x'*y/3  %需要注意的是每行看做一个观察值，每列看做一个变量也就是说，每行显示的是联合观察值加入对于第一个变量我们取1 ，那么第二个变量就肯定取5第三个变量就肯定是6

X、Y必须是各维数都相等的矩阵，其功能是把X中所有元素看莋一个变量的样本Y中所有元素看做另外一个变量的样本，把矩阵中每个对应位置看做一个联合观察值

函数实现的是求出两个变量的协方差矩阵

现在用（1,1）和（1,2）位置验证

cov(X,1)  和  cov(X,Y,1) 与之前的求解过程一致不同的是，其求出的是协方差而不是样本的协方差无偏估计值，即其除以嘚是N 而不是N-1

我要回帖

更多关于 概率论方差公式 的文章

 

随机推荐