Matlab的var函数错了吗
2016-05-17 20:10
351 查看
今天算星座符号的功率,以前分析数据的时候计算方差都直接用matlab的var函数,用的没有发现问题。今天计算的星座点数比较小,怎么算都觉得奇怪。
参考百度知道的公式,
![](http://s3.51cto.com/wyfs02/M00/7D/36/wKioL1bieyrhdMegAAAexgGGcas664.png)
没有错。难道是因为均值不为0?
对于简单的
a=[-3 -1 1 3];均值是0
求var(a)=6.6667。W..H..Y..
help var 看了一下,
方差,按照高数里面的定义,可以分为有偏估计(biased estimator)和无偏估计(unbiased estimator),差别就在于上述公式的分母是除以n还是n-1。
Matlab默认的var是计算无偏估计的。回头一看我的测试例,可不是吗,分母就是照着3除的。解决方法就是给加一个参数1,修改为有偏估计的 。
![](http://s5.51cto.com/wyfs02/M01/7D/38/wKiom1bifFCwTuLZAAAd5D_Gx5s243.png)
关于无偏估计和自由度之间的关系也可以再说道说道。
自由度(degree of freedom, df)是指当以样本的统计量来估计总体的参数时,样本中独立或能自由变化的数据的个数称为该统计量的自由度。
例如,在估计总体的平均数时,样本中的n个数全部加起来,其中任何一个数都和其他数据相独立,从其中抽出任何一个数都不影响其他数据(这也是随机抽样所要求的)。因此一组数据中每一个数据都是独立的,所以自由度就是估计总体参数时独立数据的数目,而平均数是根据n个独立数据来估计的,因此自由度为n。
但是为什么用样本估计总体的方差时,方差的自由度就是(n-1)?假设有n个样本,计算方差s^2的公式如下:
![](http://img.blog.csdn.net/20160518141332185?watermark/2/text/aHR0cDovL2Jsb2cuY3Nkbi5uZXQv/font/5a6L5L2T/fontsize/400/fill/I0JBQkFCMA==/dissolve/70/gravity/Center)
其中样本均值x_eq的计算中包括了n个独立的样本的值;或者说在已知n个独立样本值的情况下x_eq是固定的,已知x_eq的情况下,互相独立的样本值只有n-1个(另一个可以由均值推算出来),所以在这里出现了自由度减少的情况。所以无偏估计中的方差要除以N-1而不是N。
参考博文:
http://blog.163.com/wanglei2146073@126/blog/static/90689607201132294156177/ http://blog.sina.com.cn/s/blog_566762d2010097g2.html
本文出自 “妙木山大汉” 博客,请务必保留此出处/article/4416601.html
参考百度知道的公式,
![](http://s3.51cto.com/wyfs02/M00/7D/36/wKioL1bieyrhdMegAAAexgGGcas664.png)
没有错。难道是因为均值不为0?
对于简单的
a=[-3 -1 1 3];均值是0
求var(a)=6.6667。W..H..Y..
help var 看了一下,
方差,按照高数里面的定义,可以分为有偏估计(biased estimator)和无偏估计(unbiased estimator),差别就在于上述公式的分母是除以n还是n-1。
Matlab默认的var是计算无偏估计的。回头一看我的测试例,可不是吗,分母就是照着3除的。解决方法就是给加一个参数1,修改为有偏估计的 。
![](http://s5.51cto.com/wyfs02/M01/7D/38/wKiom1bifFCwTuLZAAAd5D_Gx5s243.png)
关于无偏估计和自由度之间的关系也可以再说道说道。
自由度(degree of freedom, df)是指当以样本的统计量来估计总体的参数时,样本中独立或能自由变化的数据的个数称为该统计量的自由度。
例如,在估计总体的平均数时,样本中的n个数全部加起来,其中任何一个数都和其他数据相独立,从其中抽出任何一个数都不影响其他数据(这也是随机抽样所要求的)。因此一组数据中每一个数据都是独立的,所以自由度就是估计总体参数时独立数据的数目,而平均数是根据n个独立数据来估计的,因此自由度为n。
但是为什么用样本估计总体的方差时,方差的自由度就是(n-1)?假设有n个样本,计算方差s^2的公式如下:
其中样本均值x_eq的计算中包括了n个独立的样本的值;或者说在已知n个独立样本值的情况下x_eq是固定的,已知x_eq的情况下,互相独立的样本值只有n-1个(另一个可以由均值推算出来),所以在这里出现了自由度减少的情况。所以无偏估计中的方差要除以N-1而不是N。
参考博文:
http://blog.163.com/wanglei2146073@126/blog/static/90689607201132294156177/ http://blog.sina.com.cn/s/blog_566762d2010097g2.html
本文出自 “妙木山大汉” 博客,请务必保留此出处/article/4416601.html
相关文章推荐
- 常见MATLAB小技巧
- YUVtoJPG matlab
- 使用MATLAB工具箱TOOLBOX_calib标定摄像头过程
- 记一次c#调matlab时初始化异常的解决过程
- matlab2014a代码示例 逻辑运算符不等于
- matlab2014a代码示例 一次性关闭所有的figure
- Matlab和外部文件进行数据交换相关的命令
- 【MATLAB 学习笔记】 SimMechanics 流程攻略(1)
- Matlab 内存管理
- HOG matlab练习
- Deep Learning学习 之 CNN代码解析(MATLAB)
- matlab rem与mod 的区别
- STL与Dicom坐标系对应关系(Matlab版)
- paper 69:Haar-like矩形遍历检测窗口演示Matlab源代码[转载]
- Matlab 画单点 不同线型 不同点型
- python 学习中的小技术(与matlab不同的地方)
- MATLAB地图工具箱学习总结(三)地图工具箱的基本知识
- Sobel边缘检测(matlab)
- Matlab图像处理常用函数归类总结
- matlab 字符分割