信号参数估计

一. 估计的评价准则

假设a是一个广义平稳随机信号x(n)的一个特征量,\hat a代表a的一个估计量。估计的偏差可以反应估计量与真值的接近程度,定义如下:

B=E[a-\hat a]=a-E[\hat a]

直观上,B越小,\hat a对a的估计就越好。理论上当样本数N趋于无穷大时,会形成逐渐无偏估计,如下:

\lim_{N\to \infty}B=a-\lim_{N\to\infty}E[\hat a]=0

估计的方差可以表示各次估计值相对估计均值的分散程度。估计的方差定义如下:

var[\hat a]=\sigma^2_{\hat a}=E\lbrace [\hat a-E(\hat a)]^2\rbrace=E[\hat a^2]-[E(\hat a)]^2

估计的均方误差可以综合反映估计的特性,定义如下:

E[\tilde a^2]=E[(\hat a-a)^2]

如果均方误差满足如下条件,则称其为一致估计,如下:

  1. 样本数N\to \infty
  2. \lim_{N\to\infty}E[\tilde a^2]=0

二. 一致估计

给定一致估计,证明偏差和方差都趋于零。

给定均方误差化简如下:

由条件可得:

E[\tilde a^2]=0

所以可得:

B^2+\sigma_{\hat a}^2=0

最终可得偏差和方差均为0

利用\hat a代表某算法的估计值,其他算法的估计值表示如下:

\hat a_1,\hat a_2,\cdots,\hat a_k,\cdots

如果以下不等式恒成立:

E(\hat a-a)^2\leq E(\hat a_k-a)^2

则称该估计为有效估计。

三. 估计均值

利用N代表观察次数,则平稳信号序列x(n)的观察样本如下:

x_0,x_1,\cdots,x_{N-1}

由此可计算均值的估计:

\hat m_x=\frac{1}{N}\sum_{i=0}^{N-1}x_i

3.1 偏差

首先可得:

E[\hat m_x]=E(\frac{1}{N}\sum_{i=0}^{N-1}x_i)=\frac{1}{N}\sum_{i=0}^{N-1}E[x_i]=m_x

由此可计算偏差:

B=m_x-E[\hat m_x]=0

所以此方式是一种无偏估计。

3.2 方差

根据定义可得:

E(\hat m_x^2)=E[(\frac{1}{N}\sum_{i=0}^{N-1}x_i)(\frac{1}{N}\sum_{j=0}^{N-1}x_j)]=\frac{1}{N^2}\sum_{j=0}^{N-1}\sum_{i=0}^{N-1}E[x_ix_j]

(1)当x_ix_j互不相关时,有:

E[x_ix_j]=E[x_i]E[x_j]=m_x^2

代入原式子进行化简可得:

所以该估计的方差为:

\sigma_{m_x}^2=E[\hat m_x^2]-m_x^2=\frac{1}{N}E[x_i^2]-\frac{1}{N}m_x^2=\frac{1}{N}\sigma_x^2

可得如下极限:

\lim_{N\to\infty}\sigma_{m_x}^2=0

所以该估计为无偏一致估计。 

(2)当x_ix_j相关时,有:

\sigma_{m_x}^2=E\lbrace [\hat m_x-E(\hat m_x)]^2\rbrace

进一步化简可得:
 

当i和j相差m时,可得:

E[(x-m_x)(x_j-m_x)]=cov(m)

由于在N个数据中,相距m点的数据样本有N-m对,所以可得:

 当信号数据存在相关性时,估计值的方差与协方差相关,不是一致估计,当然改变N值可以改善估计方差。

四. 估计方差

4.1 均值已知

当信号均值m_x已知,方差估计可计算得:

\hat \sigma_x^2=\frac{1}{N}\sum_{i=0}^{N-1}(x_i-m_x)^2

证明此式子为无偏一致估计。

解:

(1)首先验证偏差:

(2)接着验证一致性: 

所以,估计的方差计算为:

 4.2 均值未知

当估计的均值未知时,m_x用估值\hat m_x代替,方差可估计如下:

\hat \sigma_x^2=\frac{1}{N}\sum_{i=0}^{N-1}(x_i-\hat m_x)^2

   (1) 证明该偏差为有偏估计

(2)修改原式子,形成无偏估计

解:

(1)

很明显此为有偏估计。

(2)无偏估计的形式如下:

\hat \sigma_x^{'2}=\frac{1}{N-1}\sum_{i=0}^{N-1}(x_i-\hat m_x^2)^2

以下证明此式子为无偏估计:

显然可得:

\hat \sigma_x^{'2}=\frac{N}{N-1}\hat \sigma_x^2

对上式子两边求均值,可得:

E(\hat \sigma_x^{'2})=\frac{N}{N-1}E(\hat \sigma_x^2)=\sigma_x^2

所以B=0,此为无偏估计。 

五. 估计自相关函数

5.1 无偏自相关函数估计

估计公式为:

\hat r_{xx}(m)=\frac{1}{N-|m|}\sum_{n=0}^{N-|m|-1}x(n)x(n+m)

首先可计算:

E[\hat r_{xx}(m)]=\frac{1}{N-|m|}\sum_{n=0}^{N-|m|-1}E[x(n)x(n+m)]=r_{xx}(m)

由此可计算偏差为:

B=r_{xx}(m)-E[\hat r_{xx}(m)]=0

所以此估计为无偏估计。

估计方差的计算较复杂,可以近似可得:

当N满足如下时,方差趋于0:

N>>m,\quad N\to\infty

5.2 有偏自相关函数估计

估计公式如下:

\hat r_{xx}(m)=\frac{1}{N}\sum_{n=0}^{N-|m|-1}x(n)x(n+m)

首先可计算:

所以估计的偏差为:

B=r_{xx}(m)-E[\hat r_{xx}(m)]=\frac{|m|}{N}r_{xx}(m)

接着可计算:

若x(n)是零均值的实高斯信号,估计的方差为:

 显然可得如下极限:

 所以对于固定的m,\hat r_{xx}(m)r_{xx}(m)的一致估计。

有偏自相关函数估计式求傅氏变换:

为了利用FFT计算线性卷积,可以将x(n)扩展到2N-1点的序列,如下:

 令l=n+m,可得:

 上式子中|X_{2N}(e^{j\omega})|^2代表有限长信号的能量谱,除以N后代表功率谱。

猜你喜欢

转载自blog.csdn.net/forest_LL/article/details/124802708