摘自Steven M. Kay,《Fundamentals of Statistical Signal Processing: Estimation Theory》。

6.4 找到BLUE

  BLUE估计,其实就是无偏的具有最小方差线性估计。因此如果MVU估计是线性的,则BLUE就是MVU,否则就不是。并且BLUE有可能不存在。其中根据无偏约束,有
E(θ^)=∑n=1N−1anE(x[n])=θ.(6.2)\tag{6.2} {\rm E}(\hat \theta)=\sum_{n=1}^{N-1}a_n{\rm E}(x[n])=\theta. E(θ^)=n=1∑N−1​an​E(x[n])=θ.(6.2)θ^\hat \thetaθ^的方差为
var(θ^)=E[(∑n=1N−1anx[n]−E(∑n=1N−1anx[n]))2].{\rm var}({\hat \theta})={\rm E} \left[ \left(\sum_{n=1}^{N-1}a_nx[n] -{\rm E}\left(\sum_{n=1}^{N-1}a_nx[n]\right) \right)^2 \right]. var(θ^)=E⎣⎡​(n=1∑N−1​an​x[n]−E(n=1∑N−1​an​x[n]))2⎦⎤​.令a=[a0a1…aN−1]T{\bf a}=[a_0\ a_1\ \ldots \ a_{N-1}]^{\rm T}a=[a0​ a1​ … aN−1​]T,则θ^=aTx\hat \theta={\bf a}^{\rm T}{\bf x}θ^=aTx,因此可以得到
var(θ^)=E[(aTx−aTE(x))2]=E[(aT(x−E(x)))2]=E[aT(x−E(x))(x−E(x))Ta]=aTCa.(6.3)\tag{6.3} \begin{aligned} {\rm var}({\hat \theta})&={\rm E} \left[ \left({\bf a}^{\rm T}{\bf x}-{\bf a}^{\rm T}{{\rm E}(\bf x)} \right)^2 \right]\\ &={\rm E} \left[ \left({\bf a}^{\rm T}({\bf x}-{{\rm E}(\bf x))} \right)^2 \right]\\ &={\rm E} \left[ {\bf a}^{\rm T}({\bf x}-{{\rm E}(\bf x))} ({\bf x}-{\rm E}(\bf x))^{\rm T}{\bf a} \right]\\ &={\bf a}^{\rm T}{\bf C}{\bf a}. \end{aligned} var(θ^)​=E[(aTx−aTE(x))2]=E[(aT(x−E(x)))2]=E[aT(x−E(x))(x−E(x))Ta]=aTCa.​(6.3)
  由于满足无偏估计,因此E(x[n]){\rm E}(x[n])E(x[n])对于θ\thetaθ一定是线性的,即满足
E(x[n])=s[n]θ,(6.4)\tag{6.4} {\rm E}(x[n])=s[n]\theta, E(x[n])=s[n]θ,(6.4)这里的s[n]s[n]s[n]是已知的。否则就不可能满足无偏约束。例如,如果E(x[n])=cos⁡θE(x[n])=\cos\thetaE(x[n])=cosθ,则无偏估计应该是∑n=0N−1cos⁡anθ=θ\sum_{n=0}^{N-1}\cos a_n\theta=\theta∑n=0N−1​cosan​θ=θ。显然,不存在能够使得无偏约束成立的系数ana_nan​。
  注意如果我们将x[n]x[n]x[n]表示为
x[n]=E(x[n])+[x[n]−E(x[n])],x[n]=E(x[n])+[x[n]-E(x[n])], x[n]=E(x[n])+[x[n]−E(x[n])],如果将x[n]−E(x[n])x[n]-E(x[n])x[n]−E(x[n])看作噪声w[n]w[n]w[n],则有
x[n]=θs[n]+w[n].x[n]=\theta s[n]+w[n]. x[n]=θs[n]+w[n].(6.4)意味着BLUE适用于噪声中对已知信号进行幅度估计。为了让其应用更广泛,我们可以将其进行非线性变换(如前文所述的功率估计的例子。)
  基于(6.4),下面我们来总结估计问题。为了找到BLUE,我们需要满足无偏估计约束条件下,最小化方差
var(θ^)=aTCa.{\rm var}(\hat \theta)={\bf a}^{\rm T}{\bf C}{\bf a}. var(θ^)=aTCa.根据(6.2)和(6.4),可以得到
∑n=1N−1anE(x[n])=θ\sum_{n=1}^{N-1}a_n{\rm E}(x[n])=\theta n=1∑N−1​an​E(x[n])=θ∑n=1N−1ans[n]θ=θ\sum_{n=1}^{N-1}a_ns[n]\theta=\theta n=1∑N−1​an​s[n]θ=θ∑n=1N−1ans[n]=1\sum_{n=1}^{N-1}a_ns[n]=1 n=1∑N−1​an​s[n]=1aTs=1{\bf a}^{\rm T}{\bf s}=1 aTs=1根据Appendix 6A,可以得到最优解为
aopt=C−1ssTC−1s,{\bf a}_{\rm opt}=\frac{{\bf C}^{-1}{\bf s}}{{\bf s}^{\rm T}{\bf C}^{-1}{\bf s}}, aopt​=sTC−1sC−1s​,因此BLUE为
θ^=sTC−1xsTC−1s,(6.5)\tag{6.5} \hat \theta=\frac{{\bf s}^{\rm T}{\bf C}^{-1}{\bf x}}{{\bf s}^{\rm T}{\bf C}^{-1}{\bf s}}, θ^=sTC−1ssTC−1x​,(6.5)最小方差为
var(θ^)=1sTC−1s.(6.6)\tag{6.6} {\rm var}(\hat \theta)=\frac{1}{{\bf s}^{\rm T}{\bf C}^{-1}{\bf s}}. var(θ^)=sTC−1s1​.(6.6)注意到根据(6.4),由于E(x)=θsE({\bf x})=\theta{\bf s}E(x)=θs,BLUE为无偏的
E(θ^)=sTC−1θssTC−1s=θ.(6.5)\tag{6.5} \begin{aligned} E(\hat \theta)&=\frac{{\bf s}^{\rm T}{\bf C}^{-1}\theta {\bf s}}{{\bf s}^{\rm T}{\bf C}^{-1}{\bf s}}\\ &=\theta. \end{aligned} E(θ^)​=sTC−1ssTC−1θs​=θ.​(6.5)正如前面我们强调的,确定BLUE,只需知道

  1. s\bf ss或者标量均值;
  2. 协方差C\bf CC。
    或者说一、二阶矩,而非所有PDF。

估计理论(6):如何确定BLUE?相关推荐

  1. 测试信号分析-估计理论的基本原理

    估计理论的基本原理 随机信号序的统计量估计 随机信号序列 x(n),有限长度 均值.相关函数.功率谱密度函数的估计式为 m ^ = 1 N ∑ x ( n ) \hat m = \frac{1}{N} ...

  2. 【检测与估计理论(Detection and Estimation Theory)】一、估计简介(Introduction of Estimation)

    [检测与估计理论(Detection and Estimation Theory)]一.估计简介(Introduction of Estimation) 估计理论的应用 估计问题 估计量的性能评估 参 ...

  3. 滤波估计理论(一)——贝叶斯滤波

    滤波估计理论(一)--贝叶斯滤波(Bayesian Filtering) 估计问题的建模 状态空间模型 概率模型 贝叶斯估计方法 批处理贝叶斯估计 预测,滤波还是平滑? 贝叶斯滤波 状态预测 量测更新 ...

  4. NormalEstimation法向量估计理论和代码---PCL源码笔记

    NormalEstimation 确定表面一点法线的问题近似于估计临近点所拟合表面的一个相切面法线的问题,因此转换过来以后就变成了一个最小二乘法平面拟合的问题.因此估计表面法线的解决方案就变成了分析一 ...

  5. 信号检测与估计理论_论文解读 | 利用脑功能连接实现疲劳驾驶检测

    ©PaperWeekly 原创 · 作者|张玮玮 学校|东北大学硕士生 研究方向|脑电情绪识别 论文标题:Driving Fatigue Recognition with Functional Con ...

  6. PCL:点云数据基于法线的边界提取(从最初的法线估计理论推导到最终的边界提取)

    该边界提取采用PCL库里边的方法,基于法线估计来实现的边界检测与提取: 首先从原始点云上计算出法线,再由法线结合数据估计出边界.(这样解释还是特别抽像吧) ------------法线求解:(平面的法 ...

  7. 信号检测与估计理论 pdf_CVPR2020|行人检测与重识算法推荐论文源码大盘点

    本文收集了CVPR 2020 一些行人检测与人员重识别优秀论文,我们知道在视频监控相关领域这些技术方向可以得到很好得广泛应用. 行人检测及人群计数从内容来看主要解决行人与行人.行人与物体间的遮挡透视, ...

  8. 信号检测与估计理论 pdf_目标检测的性能上界讨论

    加入极市专业CV交流群,与6000+来自腾讯,华为,百度,北大,清华,中科院等名企名校视觉开发者互动交流!更有机会与李开复老师等大牛群内互动! 同时提供每月大咖直播分享.真实项目需求对接.干货资讯汇总 ...

  9. 估计理论(3):充分统计量的完备性

    文章目录 1.CRLB [定理3.1:标量参数的CRLB] 2.线性模型 [定理4.1:线性模型的最小方差无偏估计] 3.一般最小方差无偏估计 [定理5.2:Rao-Blackwell-Lehmann ...

  10. 【 Notes 】Best linear unbiased estimator(BLUE) approach for time-of-arrival based localisation

    目录 Abstract Introduction BLUE-based positioning BLUE-LSC algorithm BLUE-LLS algorithm Abstract A com ...

最新文章

  1. Ubuntu安装okular PDF阅读器
  2. 简单记线程安全的集合类
  3. 大牛书单 | 春节特辑
  4. MVC下用C#实现Excel导出
  5. mooc作业怎么上传附件_怎么以最低价格入门hifi音响?
  6. Posix线程编程指南(4) 线程终止
  7. [C#] 如何分析stackoverflow等clr错误
  8. 个人信用报告“数字解读”与妙用
  9. Hive近百个常规函数详解
  10. 基于虚拟日志压缩的数据同步方案
  11. 分布式监控系统Zabbix--使用Grafana进行图形展示
  12. 2016/1/15代码
  13. Des加密后传参被特殊字符(如+)截断
  14. oracle 10g for solaris准备工作和实例的删除
  15. 客户机操作系统已禁用 cpu_强实时工业互联网虚拟化操作系统Intewell
  16. 2023宁波大学计算机考研信息汇总
  17. IDC发布2021年中国云计算10大预测;Docker 桌面为 M1 推出技术预览版
  18. 结合可变形注意力的视觉Transformer
  19. re:从零开始的Qt世界冒险生活
  20. C# 反射类Assembly用法举例

热门文章

  1. 算法设计与分析复习——第四章:贪心算法
  2. mysql 数据备份 crontab
  3. 3.4 常用的两种 layer 层 3.7 字体与文本
  4. Ubuntu下PHP开发环境LAMP搭建
  5. Spring已集成jsp的环境下同时集成Velocity
  6. 14.1宽屏Intel集显调分辨率:无解
  7. python实用例子_Python实用案例 - 随笔分类 - 一入测试深似海 - 博客园
  8. string不能输入空格,如何输入有空格字符串呢
  9. mysql只导出表数据_MySQL 如何只导出 指定的表 的表结构和数据 ( 转 )
  10. mysql 数据增量备份_mysqlmysqldump数据备份和增量备份