有偏和无偏的自协方差序列估计是什么意思?

信息处理 频谱
2022-02-19 16:24:03

在非参数谱估计方法中,有偏和无偏自协方差序列估计是什么意思,为什么使用最常见的有偏估计?我有 Peter Stoica 的第二个问题(信号光谱分析)的答案,但没有得到他想说的。

2个回答

在光谱估计中,人们估计“光谱”和其他值。并且有兴趣知道估计是否x^收敛到“真实”但未知的值 x. 它是未知的,例如,因为一个人只能访问有限的数量N要计算的样本数量,实现的数量有限,因为不能保证遍历性等。所以我们想知道是否 x^Nx作为N. 收敛可以用不同的方式来思考,最常见的是均方意义:

limNE(|x^Nx|2)=0
在哪里E是期望吗?在实践中证明这可能非常复杂。因此,当满足更简单的条件时,人们通常会感到满意。并且有两个必要条件:

limNE(x^N)=x

limNE(|x^N|2)=0

第一个是渐近无偏性,第二个是消失的方差。它们共同定义了一致的估计量。因此,有偏差的自相关估计不会收敛到真正的自相关。派生的估计量(如周期图)也不太可能收敛。

只是为了提供不同的视角,这里有一个基本的解释,试图在没有任何理论的情况下给出一些直觉。

假设我们有N零均值随机过程的样本x(n),n=0...N1我们想估计R(d)=E(x(n)x(nd)). 一种明显的方法是平均x(n)x(nd)对于我们拥有的所有数据:

R^(d)=1Ndn=dN1x(n)x(nd)

这是一个无偏估计,因为随着您使用越来越多的数据(更大N)。但更大滞后的估计d将不如对较小滞后的估计准确,因为我们没有可用的数据来作为对大滞后的估计的基础。

另一种方法是将估计标准化1/N每次:

R^(d)=1Nn=dN1x(n)x(nd)

这是一个有偏差的估计,因为即使您使用大量数据,它也不会收敛到真实值。较大滞后的估计d是基于Nd样本,但我们通过1/N不管滞后,所以估计有点太小了。但是相对于更准确的估计(具有较小滞后的那些),不太准确的估计(具有较大滞后的估计)会有所减弱,这会有所帮助。它通常也更容易计算:事实上,周期图估计器(例如使用 DFT 的估计器)可以被视为使用有偏自相关。