统计信号估计 (二) 最小方差无偏估计量(MUV)的线性模型

MVU表示最小方差无偏估计,下面将推导一个线性模型的MVUE,即最小方差无偏估计量。

线性模型可以表示为:
x = H θ + w \bf{x = H\theta + w} x=Hθ+w
其中 θ \theta θ 为待估参数,均为向量。当然,这里的MVU估计量也可以被称为BLUE,即最佳线性无偏估计量。
我们假设 w ∼ N ( 0 , C ) \bf{w} \sim N(0,C) wN(0,C) ,要根据克拉美罗界限求得最小方差无偏估计量,依据公式:
∂ ln ⁡ p ( x ; θ ) ∂ θ = I ( θ ) ( g ( x ) − θ ) \frac{\partial \ln p(\mathbf{x} ; \theta)}{\partial \theta}=I(\theta)(g(x)-\theta) θlnp(x;θ)=I(θ)(g(x)θ)
可以得到克拉美罗界限 I ( θ ) I(\theta) I(θ),以及最小方差无偏估计量 θ ^ = g ( x ) \hat{\theta} = g(x) θ^=g(x),所以我们的目的就是求 ∂ ln ⁡ p ( x ; θ ) ∂ θ \frac{\partial \ln p(\mathbf{x} ; \theta)}{\partial \theta} θlnp(x;θ)
p ( x ; θ ) = 1 ( 2 π ) N 2 d e t ( C ) 1 2 e x p ( − 1 2 ( x − H θ ) T C − 1 ( x − H θ ) ) p(\bf{x;\theta}) = \frac{1}{ (2\pi)^{\frac{N}{2}}det(C)^{\frac{1}{2}}}exp(-\frac{1}{2}(x-H\theta)^{T}C^{-1}(x-H\theta)) p(x;θ)=(2π)2Ndet(C)211exp(21(xHθ)TC1(xHθ))
∂ ln ⁡ p ( x ; θ ) ∂ θ = − 1 2 ∂ ∂ θ [ x T C − 1 x + ( H θ ) T C − 1 ( H θ ) − ( H θ ) T C − 1 x − x T C − 1 ( H θ ) ] = − 1 2 [ 2 H T C − 1 H θ − 2 H T C − 1 x ] = − H T C − 1 H θ + H T C − 1 x = ( H T C − 1 H ) ( ( H T C − 1 H ) − 1 H T C − 1 x − θ ) \begin{aligned} \frac{\partial \ln p(\mathbf{x} ; \theta)}{\partial \theta} &= -\frac{1}{2}\frac{\partial}{\partial \theta}[x^TC^{-1}x +(H \theta)^TC^{-1}(H \theta) - (H \theta)^TC^{-1}x-x^TC^{-1}(H \theta)] \\ &= - \frac{1}{2}[2 H^TC^{-1}H\theta -2H^TC^{-1}x] \\ &=-H^T C^{-1}H\theta + H^TC^{-1}x \\ &=(H^TC^{-1}H)((H^TC^{-1}H)^{-1}H^TC^{-1}x - \theta) \end{aligned} θlnp(x;θ)=21θ[xTC1x+(Hθ)TC1(Hθ)(Hθ)TC1xxTC1(Hθ)]=21[2HTC1Hθ2HTC1x]=HTC1Hθ+HTC1x=(HTC1H)((HTC1H)1HTC1xθ)

得到 θ ^ = g ( x ) = ( H T C − 1 H ) − 1 H T C − 1 x \hat{\theta} = g(x) = (H^TC^{-1}H)^{-1}H^TC^{-1}x θ^=g(x)=(HTC1H)1HTC1x
v a r ( θ ) = I ( θ ) − 1 = ( H T C − 1 H ) − 1 var(\theta) = I(\theta)^{-1} = (H^TC^{-1}H)^{-1} var(θ)=I(θ)1=(HTC1H)1
这个结论与最小二乘法相同。

  • 2
    点赞
  • 4
    收藏
    觉得还不错? 一键收藏
  • 打赏
    打赏
  • 0
    评论

“相关推荐”对你有帮助么?

  • 非常没帮助
  • 没帮助
  • 一般
  • 有帮助
  • 非常有帮助
提交
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包

打赏作者

肆拾伍

你的鼓励将是我创作的最大动力

¥1 ¥2 ¥4 ¥6 ¥10 ¥20
扫码支付:¥1
获取中
扫码支付

您的余额不足,请更换扫码支付或充值

打赏作者

实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值