信号パラメータの推定

I. 推定される評価基準

a が一般化された定常ランダム信号 x(n) の特徴量であり、\ハットa の推定量を表すものとします。推定された偏差は、次のように定義される推定値の真の値への近さを反映することができます。

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

直感的には、B が小さいほど、\ハット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 aE(\hat a)]^2\rbrace=E[\hat a^2]-[E(\hat a )]^2

推定された平均二乗誤差は、推定値の特性を包括的に反映できます。推定値は次のように定義されます。

E[\チルダ a^2]=E[(\hat aa)^2]

平均二乗誤差が次の条件を満たす場合、それは次のように一貫した推定値と呼ばれます。

  1. サンプル数N\から\infty
  2. \lim_{N\to\infty}E[\tilde a^2]=0

2. コンセンサス推定

一貫した推定値が与えられた場合、バイアスと分散の両方がゼロになる傾向があることを示します。

与えられた平均二乗誤差は次のように単純化されます。

条件によって取得されます:

E[\チルダa^2]=0

したがって、次のものを取得できます。

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

最終的に利用可能なバイアスと分散は両方とも 0 です

1 つのアルゴリズムを表す推定値を使用すると\ハット、他のアルゴリズムの推定値は次のように表されます。

\hat a_1、\hat a_2、\cdots、\hat a_k、\cdots

次の不等式が一定である場合:

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

この推定値は有効推定値と呼ばれます。

3. 推定平均値

観測数を 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_iと がx_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 点で区切られた Nm 個のデータ サンプルのペアがあるため、次のように取得できます。

 信号データに相関がある場合、推定値の分散は共分散に関係しており、一貫した推定値ではありませんが、もちろん N の値を変更すると推定分散を改善できます。

4. 推定分散

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. 自己相関関数の推定

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}(男)y は のr_{xx}(男)一貫した推定値になります。

偏自己相関関数の推定式はフーリエ変換を求めます。

FFT を使用して線形畳み込みを計算するには、次のように x(n) を 2N-1 点のシーケンスに拡張できます。

 l=n+m とすると、次のようになります。

 上式は|X_{2N}(e^{j\omega})|^2有限信号のエネルギースペクトルを表し、Nで割ったものがパワースペクトルを表します。

おすすめ

転載: blog.csdn.net/forest_LL/article/details/124802708