协方差的定义

对于一般的分布,直接代入E(X)之类的就可以计算出来了,但真给你一个具体数值的分布,要计算协方差矩阵,根据这个公式来计算,还真不容易反应过来。这里用一个例子说明协方差矩阵是怎么计算出来的吧。

记住,X、Y是一个列向量,它表示了每种情况下每个样本可能出现的数。比如给定

则X表示x轴可能出现的数,Y表示y轴可能出现的。注意这里是关键,给定了4个样本,每个样本都是二维的,所以只可能有X和Y两种维度。所以

用中文来描述,就是:

协方差(i,j)=(第i列的所有元素-第i列的均值)*(第j列的所有元素-第j列的均值)

这里只有X,Y两列,所以得到的协方差矩阵是2x2的矩阵,下面分别求出每一个元素:

所以,按照定义,给定的4个二维样本的协方差矩阵为:

用matlab计算这个例子

z=[1,2;3,6;4,2;5,2]

cov(z)

ans =

2.9167   -0.3333

-0.3333    4.0000

可以看出,matlab计算协方差过程中还将元素统一缩小了3倍。所以,协方差的matlab计算公式为:

协方差(i,j)=(第i列所有元素-i列均值)*(第j列所有元素-j列均值)/(样本数-1

下面在给出一个4维3样本的实例,注意4维样本与符号X,Y就没有关系了,X,Y表示两维的,4维就直接套用计算公式,不用X,Y那么具有迷惑性的表达了。

(3)与matlab计算验证

Z=[1 2 3 4;3 4 1 2;2 3 1 4]

cov(Z)

ans =

1.0000    1.0000   -1.0000   -1.0000

1.0000    1.0000   -1.0000   -1.0000

-1.0000   -1.0000    1.3333    0.6667

-1.0000   -1.0000    0.6667    1.3333

可知该计算方法是正确的。我们还可以看出,协方差矩阵都是方阵,它的维度与样本维度有关(相等)。参考2中还给出了计算协方差矩阵的源代码,非常简洁易懂,在此感谢一下!

covariance matrix相关推荐

  1. R语言使用psych包进行探索性因子分析EFA、使用cov2cor函数将原始数据的协方差矩阵将其转换为相关性矩阵( covariance matrix into correlation matrix)

    R语言使用psych包进行探索性因子分析EFA.使用cov2cor函数将原始数据的协方差矩阵将其转换为相关性矩阵( covariance matrix transform into correlati ...

  2. 【Matlab】错误使用 classify (line 233) The pooled covariance matrix of TRAINING must be positive definite.

    在 Matlab 用 Classify 函数做判别分析时,有时会碰到下面的问题: 错误使用 classify (line 233) The pooled covariance matrix of TR ...

  3. 加权协方差矩阵(weighted covariance matrix)

    国内完全没一个有用的,这里给出了加权协方差矩阵计算函数.用的时候可以将权重先归一化. def weighted_cov(values, weights):"""Compu ...

  4. 散布矩阵(Scatter Matrix)及其与协方差矩阵(The Covariance Matrix)的关系

    在多元统计和概率论中,散点矩阵是一种统计量,用来估计协方差矩阵,例如多元正态分布. In multivariate statistics and probability theory, the sca ...

  5. 实验笔记之——covariance matrix pooling

    本博文为协方差池化的实验笔记.具体理论请参考博客<学习笔记之--covariance pooling> 实验python train.py -opt options/train/train ...

  6. covariance matrix r语言_时间序列分析|ARIMAX模型分步骤详解和R中实践

    这是关于时间序列的第N篇文章,本文将介绍ARIMAX模型,简单来说就是在ARIMA的基础上增加一个外生变量.ARIMAX和ARIMA相比在理论上没有太多新的内容,所以本文直接介绍在R里怎么一步一步跑A ...

  7. 协方差矩阵(covariance matrix)

    Xn×d⇒(XTX)d×d X_{n\times d}\Rightarrow \left (X^TX\right )_{d\times d} (1)协方差矩阵:半正定(semi-positive de ...

  8. A geometric interpretation of the covariance matrix(reproduced)

    转载自:http://www.visiondummy.com/2014/04/geometric-interpretation-covariance-matrix/

  9. Why does the Eigen decomposition of the covariance matrix of a point cloud give its orientation?

    数学分析 PCA 协方差矩阵的几何角度解释 我见过的使用PCA估计点云法向量的最好的解释

最新文章

  1. jvm运行时数据区是干啥的?CPU切换线程会不会突然忘记程序执行到哪一步了
  2. improve php,解析提高PHP执行效率
  3. 【Uva 10934】Dropping water balloons
  4. 使用sshfs挂载linux远程服务器目录到windows
  5. mybatis动态SQL配置文件
  6. 广播中等待较久出现anr问题
  7. 博文视点读书节第八日丨大咖书单升级放送,今日直播助你成为Offer收割机~晒单赢福袋正在派送中!
  8. Python在线编程网站
  9. jsp项目实战--项目介绍
  10. 又一个PS2汉化入门
  11. 视频教程-思科CCNP路由实验专题--EIGRP篇【适用于思科CCNA、CCNP】-网络技术
  12. win10系统蓝牙耳机的连接与连接后无法听到声音
  13. FZU11685 之 跑跑卡丁车
  14. 查询显卡状态 Tesla M40
  15. 南宁西乡塘区的计算机中专学校,南宁市西乡塘区中专学校
  16. 如何解决学习进度缓慢?
  17. 一些句子【未完待续】
  18. 关于PN532 读取二代证UUID____记录自己亲自实验成功!!!
  19. unity跑酷怎么添加金币_【Unity3D实战】零基础一步一步教你制作跑酷类游戏(填坑完整版)...
  20. 读取excel中数据时,数字格式发生改变

热门文章

  1. python开源的etl工具_一个开源的ETL工具Airflow
  2. windows键盘快捷键
  3. 代码静态分析工具-splint的学习与使用[转]
  4. 马克思主义基本原理概论(二)复习资料
  5. 推荐一个 Android 图片压缩框架
  6. 各种未解决配置烦心事:
  7. Windows 10右下角系统时间如何显示长时间?年月日时分秒?星期几?
  8. LVM磁盘格式分区创建挂载并刷新存储
  9. SEO策略全面分析与总结
  10. python验证哥德巴赫猜想之一