色付きノイズが存在する場合の最尤推定


7

測定ノイズ(1)ホワイトガウスノイズ(2)カラーノイズ-ピンク、バイオレットの存在下でシステムの識別をテストしようとしています。パラメータを推定するとき、iidの存在下で推定すると、ゼロは無相関ノイズを意味します。

Q1:色付きノイズが相関しているかどうか知りたいのですが。それらは異なる分布をしていると思いますが、サンプルが相関するかどうかについては、情報を見つけることができませんでした。

Q2:推定では、ノイズは相関のないiidである加法性ホワイトガウスノイズであると想定しています。ノイズがガウスでない場合はどうなりますか?例:バツ=sθ+Coloredose 見積もろうとしているところ θ。パフォーマンス、つまりMSEは、色付きノイズと色なしノイズのレベルによって異なりますか?


1
「着色」ノイズはパワースペクトルではなく、流通や、実際にPDFファイルについてです、条件付きPDFを行い、スペクトルとは何かを持っていますが、ない無条件PDF
ロバート・ブリストウ・ジョンソン

@ robertbristow-johnson:以下の返信から、pdfは常にカラーノイズのガウスと見なされますか?
Ria George

1
常にではないが、通常。単純な反例は、ハイパス三角PDFディザです。y[] 均一なPDF「ホワイト」ノイズから形成 バツ[](良いrand()関数から何が得られるか)
y[]=バツ[]バツ[1]
これはカラーノイズ(低周波数ではスペクトルが少ない)ですが、ガウスpdfではなく、三角pdfです
robert bristow-johnson '18 / 06/18

回答:


5

ノイズプロセスの自己相関関数は、ホワイトノイズの場合のようにデルタ関数ではないため、(異なる時間に取得された)色付きノイズのサンプルは一般 相関確率変数です。したがって、ゼロ平均プロセス(ノイズは通常、その色に関係なく想定される)を仮定すると、2つの信号の共分散は、τ 秒は Rτ どこ Rt=F1Sfプロセスの自己相関関数です(パワースペクトル密度の逆フーリエ変換)。可能であることに注意してくださいRtいくつかの値に対してゼロになるt (例えば Rt=tは有効な自己相関関数です)が、ゼロ以外のすべてに対してゼロにすることはできませんt

限り密度プロセスがガウス分布である場合、任意のサンプルの関数、サンプルプロセスがで濾過されていてもガウス分布であるリニアサンプリングする前にフィルタ。プロセスがある場合でもない ガウス各サンプルはラプラシアンないあらゆる種類のフィルタ処理した後、一般的に、プロセスのサンプルの言うことが同じ缶になり、その後ながら、(それは、私たちは、ラプラシアンとしましょう)。つまり、Gaussianityは線形フィルタリングを生き延びますが、LaPlacismは一般的に生き残りません。


では、サンプルに相関ノイズがある場合、最尤推定はどのように機能しますか?aの未知の平均を推定したい場合を考えてみましょう。Nμ1 確率変数、2つの観測があります バツ そして y。独立した観測の標準的な場合、尤度関数は次のとおりです。

Lμ=12πexp12[バツμ2+yμ2]
_maximum-likelihood estimator for μ 数です μ^ 最大化する Lμ、それは数になることがわかります μ^ その最小化 バツμ2+yμ2。これは二次式ですμ そして、最尤推定値は μ^=バツ+y2。観測値が相関係数と相関している場合ρ、その後
Lμ=12π1ρ2exp121ρ2[バツμ22ρバツμyμ+yμ2]
もう一度見つける必要があります μ^ どこ バツμ22ρバツμyμ+yμ2最小値があります。私たちはまだ二次方程式を持っていますμ しかし今、私たちは次のような用語を得ます バツy係数で。何μ^ あなたがワークアウトするために残されています。

そして、もし 観察 >2?上記すべてが引き続き適用されます。サンプル内の独立して同一に分布したガウスノイズの場合、サンプルの平均 1Σバツ は、の最尤推定値です μ しかし、相関ガウス確率変数の場合、最小化しようとしている二次式は共分散行列の逆に依存し、結果は単純な簡単なではなくデータの非線形関数であるため、非常に厄介な最小化問題が発生します。サンプル平均のような結果を思い出してください。


ノイズがガウスでない場合はどうなりますか?同じ原則が適用されます-尤度関数を設定し、それが最大値に到達する場所を見つけます-しかし、計算はかなり異なります。

弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.