MMSE估计(二):从估计到估计器

2024-03-13 18:48
文章标签 估计 mmse

本文主要是介绍MMSE估计(二):从估计到估计器,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

【本文内容摘自"Signals, Systems and Inferences"之"8.2-From Estimate to an Estimator", by Alan V.Oppenheim and George C.Verghese, 2010.】
  上面我们得到对于特定的 X = x X=x X=x,有
y ^ ( x ) = E [ Y ∣ X = x ] . \hat y(x)={\rm E}\left[Y|X=x\right]. y^(x)=E[YX=x].
然而,尽管 X X X是随机变量,但其取值 x x x不是,因此 y ^ ( x ) \hat y(x) y^(x)也不是。
  在我们进入更深入的讨论时,我们有必要在估计随机变量和实现估计的步骤间划分出界限。这就好像对于函数的取值和函数本身进行区分。我们把产生估计的过程或者函数称为估计器(estimator)。
  我们用 y ^ ( X ) \hat y(X) y^(X)来表示MMSE估计器,即
Y ^ = y ^ ( X ) = E [ Y ∣ X ] . \hat Y=\hat y({\bf X})={\rm E}[Y|{\bf X}]. Y^=y^(X)=E[YX].
注意这里 E [ Y ∣ X ] {\rm E}[Y|{\bf X}] E[YX]表示一个随机变量,而非一个数值。
  给定 X {\bf X} X时对Y的MMSE估计器将最小化MSE,这是由于MMSE估计器对于每个特定的值 x x x,都最小化MSE,即
E Y , X ( [ Y − y ^ ( X ) ] 2 ) = E X { E Y , X ( [ Y − y ^ ( X ) ] 2 ∣ X ) } = ∫ − ∞ ∞ E Y , X ( [ Y − y ^ ( X ) ] 2 ∣ X = x ) f X ( x ) d x . {\rm E}_{Y,{\bf X}}\left( [Y-\hat y({\bf X})]^2\right)={\rm E}_{{\rm X}}\left\{{\rm E}_{Y,{\bf X}}\left( [Y-\hat y({\bf X})]^2|{\bf X}\right)\right\}\\ =\int_{-\infty}^{\infty}{\rm E}_{Y,{\bf X}}\left( [Y-\hat y({\bf X})]^2|{\bf X=x}\right)f_{\bf X}({\bf x})d{\bf x}. EY,X([Yy^(X)]2)=EX{EY,X([Yy^(X)]2X)}=EY,X([Yy^(X)]2X=x)fX(x)dx.
由于估计 y ^ ( x ) \hat y({\bf x}) y^(x)对于 X X X的每个值都能够最小化 E Y ∣ X {\rm E}_{Y|{\bf X}} EYX,因此也能够最小化 E X {\rm E}_{\bf X} EX,因为 f X ( x ) f_{\bf X}({\bf x}) fX(x)非负。

正交性
  MMSE估计器的另外一个重要特性时残留误差 Y − y ^ ( X ) Y-\hat y({\rm X}) Yy^(X)与测量的随机变量的任意函数 h ( X ) h({\bf X}) h(X)正交,即
E Y , X { [ Y − y ^ ( X ) ] h ( X ) } = 0. ( 1 ) {\rm E}_{Y,X}\{[Y-\hat y ({\bf X)}]h({\bf X})\}=0.\qquad \qquad (1) EY,X{[Yy^(X)]h(X)}=0.(1)
因此
E Y , X { y ^ ( X ) h ( X ) } = E Y , X { Y h ( X ) } . ( 2 ) {\rm E}_{Y,X}\{\hat y ({\bf X)}h({\bf X})\}={\rm E}_{Y,X}\{Yh({\bf X})\}.\qquad \qquad (2) EY,X{y^(X)h(X)}=EY,X{Yh(X)}.(2)
特别地,选择 h ( X ) = 1 h({\bf X})=1 h(X)=1,有
E Y , X { y ^ ( X ) } = E Y { Y } . {\rm E}_{Y,X}\{\hat y ({\bf X)}\}={\rm E}_{Y}\{Y\}. EY,X{y^(X)}=EY{Y}.
因此,估计器被称为无偏的(unbiased):它的期望值等于估计的随机变量的期望值。我们可以用无偏性来理解(1),即MMSE估计器的估计误差与用来构造估计器的随机变量的任意函数都是无关的。
(2)的证明如下
E Y , X { y ^ ( X ) h ( X ) } = E X [ E Y ∣ X ( Y ∣ X ) h ( X ) } = E X [ E Y ∣ X [ Y h ( X ) ∣ X } = E Y , X { Y h ( X ) } {\rm E}_{Y,X}\{\hat y({\bf X})h({\bf X})\}={\rm E}_{X}[ {\rm E}_{Y|X}(Y|{\bf X})h({\bf X})\}\\={\rm E}_{X}[ {\rm E}_{Y|X}[Yh({\bf X})|{\bf X}\}={\rm E}_{Y,X}\{Yh({\bf X})\} EY,X{y^(X)h(X)}=EX[EYX(YX)h(X)}=EX[EYX[Yh(X)X}=EY,X{Yh(X)}
由此可以得到(1)中的正交性。

这篇关于MMSE估计(二):从估计到估计器的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/805845

相关文章

数据集 3DPW-开源户外三维人体建模-姿态估计-人体关键点-人体mesh建模 >> DataBall

3DPW 3DPW-开源户外三维人体建模数据集-姿态估计-人体关键点-人体mesh建模 开源户外三维人体数据集 @inproceedings{vonMarcard2018, title = {Recovering Accurate 3D Human Pose in The Wild Using IMUs and a Moving Camera}, author = {von Marc

数据集 Ubody人体smplx三维建模mesh-姿态估计 >> DataBall

Ubody开源人体三维源数据集-smplx-三维建模-姿态估计 UBody:一个连接全身网格恢复和真实生活场景的上半身数据集,旨在拟合全身网格恢复任务与现实场景之间的差距。 UBody包含来自多人的现实场景的1051k张高质量图像,这些图像拥有2D全身关键点、3D SMPLX模型。 UBody由国际数字经济学院(IDEA)提供。 (UBody was used for mesh r

Kaggle克隆github项目+文件操作+Kaggle常见操作问题解决方案——一文搞定,以openpose姿态估计项目为例

文章目录 前言一、Kaggle克隆仓库1、克隆项目2、查看目录 二、安装依赖三、文件的上传、复制、转移操作1.上传.pth文件到input目录2、将权重文件从input目录转移到工作目录 三、修改工作目录里的文件内容1、修改demo_camera.py内容 四、运行! 前言 想跑一些深度学习的项目,但是电脑没有显卡,遂看向云服务器Kaggle,这里可以每周免费使用30h的GP

【译】PCL官网教程翻译(18):估计一组点的视点特征直方图(VFH)签名 - Estimating VFH signatures for a set of points

英文原文查看 估计一组点的视点特征直方图(VFH)签名 本文描述了视点特征直方图([VFH])描述符,这是一种针对聚类(如对象)识别和6DOF姿态估计问题的点簇表示方法。 下图展示了一个VFH识别和姿态估计的例子。给定一组火车数据(除最左边的点云外,最上面一行、最下面一行),学习一个模型,然后使用一个云(最左边的部分)查询/测试模型。匹配的结果按从最好到最差的顺序从左到右从左下角开始。有关更多

分歧时间估计与被子植物的年代-文献精读43

Ad fontes: divergence-time estimation and the age of angiosperms 回归本源:分歧时间估计与被子植物的年代 摘要 准确的分歧时间对于解释和理解谱系演化的背景至关重要。在过去的几十年里,有关冠被子植物推测的分子年龄(通常估计为晚侏罗世至二叠纪)与化石记录(将被子植物置于早白垩纪)之间的差异,引发了广泛的争论。如果冠被子植物早在二

SLAM ORB-SLAM2(29)PnP估计姿态

SLAM ORB-SLAM2(29)PnP估计姿态 1. PnP问题2. EPnP算法2.1. 计算4对控制点的世界坐标2.2. 计算齐次质心坐标2.3. 计算4对控制点的相机坐标2.3.1. 构造M矩阵2.3.2. 计算 M T M M^TM MTM的0特征值对应的特征向量2.3.3. 计算零空间的秩2.3.4. 计算线性组合的系数 2.4. 选择最小重投影误差 3. 标题

损失函数、成本函数cost 、最大似然估计

一、损失函数 什么是损失函数? 【深度学习】一文读懂机器学习常用损失函数(Loss Function)-腾讯云开发者社区-腾讯云 损失函数(loss function)是用来估量模型的预测值f(x)与真实值Y的不一致程度,它是一个非负实值函数,通常使用L(Y, f(x))来表示,损失函数越小,模型的鲁棒性就越好。损失函数是经验风险函数的核心部分,也是结构风险函数重要组成部分。模型的结构风险函

End-to-End视觉里程计新突破:从运动模糊图像中精确估计相机姿态

更多优质内容,请关注公众号:智驾机器人技术前线 1.论文信息 论文标题:MBRVO: A Blur Robust Visual Odometry Based on Motion Blurred Artifact Prior 作者:Jialu Zhang, Jituo Li*, Jiaqi Li, Yue Sun, Xinqi Liu, Zhi Zheng, and Guodong Lu

【MATLAB源码-第253期】基于matlab的8PSK调制载波+相位+符号定时联合估计仿真,输出星座图等。

操作环境: MATLAB 2022a 1、算法描述 1. 系统背景和目标 8PSK是一种调制方式,其中信号的相位被分成8个不同的状态,每个状态代表3比特的数据。这个过程涉及将比特序列转换为相应的相位,经过调制后传输给接收端。在接收端,由于信号通过信道传播过程中可能受到噪声、频偏、符号定时偏差等因素的影响,接收到的信号相对于原始信号会有一定的失真。因此,接收端需要进行一系列处理来恢复原始数据

熵的相关概念及相互关系(信息熵,条件熵,相对熵,交叉熵,最大似然估计)

熵:系统混乱程度的度量,系统越混乱,熵越大。 信息熵:信息量的大小的度量,用于描述随机变量的不确定度。事件的不确定性越大,则信息量越大,信息熵越大。定义如下: 条件熵:表示在已知随机变量X的条件下随机变量Y的不确定性。定义如下: 另外,,说明描述X和Y所需的信息(H(X,Y) ࿰