自己相関関数は確率過程を完全に記述していますか?


31

確率過程はその自己相関関数によって完全に記述されていますか?

そうでない場合、どの追加プロパティが必要になりますか?

回答:


44

確率過程の完全な説明とはどういう意味ですか?数学的には、確率過程はランダム変数のコレクションで、インデックスセット各瞬間に1つずつ、通常は実数行全体または正の実数行であり、完全な説明は、各整数および瞬間について、(結合)分布を知っていることを意味しますの確率変数、、{X(t):tT}t TTn1nt1,t2,,tnTnX(t1)X(t2),X(tn)X T T X T 1X T 2、T 1 t 2 X t 1X t。これは、莫大な情報量:我々はCDF知っている必要があります、各時刻のための、の(二次元)共同CDFとの時刻のすべての選択肢についておよび、、、およびなどの(3次元)CDF など。X(t)tX(t1)X(t2)t1t2X(t1)X(t2)X(t3)

したがって、当然のことながら、人々はより単純な説明とより制限的なモデルを探してきました。プロセスが時間原点の変更に対して不変である場合、1つの単純化が発生します。これが意味することは

  • プロセス内のすべてのランダム変数は同一のCDFを持っています:すべてのに対してです。FX(t1)(x)=FX(t2)(x)t1,t2
  • 任意の二つの確率変数は、同じ関節CDFた時間のいくつかの指定された量によって分離された任意の他の確率変数の対により分離同じ時間の量。たとえば、ランダム変数およびは、ランダム変数およびと同様に秒で区切られているため、X(t1)X(t1+τ)τX(t2)X(t2+τ)FX(t1),X(t1+τ)(x,y)=FX(t2),X(t2+τ)(x,y)
  • 任意の3の確率変数、、離間及び離れて同じ関節CDFを有する、、どのようにも離間と離れ、X(t1)X(t1+τ1)X(t1+τ1+τ2)τ1τ2X(t2)X(t2+τ1)X(t2+τ1+τ2)τ1τ2
  • すべての多次元CDF について同様です。たとえば、多次元のケースの詳細については、Peter K.の回答を参照してください。

事実上、ランダムプロセスの確率的記述は、時間軸上の原点と呼ぶものに依存しません。すべての時刻を一定量からだけシフトします。は、確率変数と同じ確率的記述を提供します。この特性は、厳密な意味での定常性と呼ばれ、この特性を享受するランダムプロセスは、厳密に定常的なランダムプロセス、またはより簡単には、定常的なランダムプロセスと呼ばれます。 t1,t2,,tnτt1+τ,t2+τ,,tn+τ

厳密な定常性自体は、特定の形式のCDFを必要としないことに注意してください。たとえば、すべての変数がガウス分布であるとは言いません。

形容詞、より緩やかな形の定常性を定義できることを厳密に示唆しています。場合の関節CDF -order 同じである関節-order CDFすべての選択肢のための及び、ランダムプロセスがあると言われています、順番に静止ととも呼ばれる -order定常ランダムプロセス。次の定常ランダムプロセスも、各正の次定常であることに注意してください。 NthX(t1),X(t2),,X(tN)NthX(t1+τ),X(t2+τ),,X(tN+τ)t1,t2,,tNτNNthNthnn<N N 番目の N 番目の N-N F XX= LIM Y F X YXYN。(これは、引数のが近づくにつれて、次の結合CDFが次のCDF の制限であるためです:一般化 )。厳密に定常的なランダムプロセスは、すべての次数に対して定常的なランダムプロセスです。nthNthNnFX(x)=limyFX,Y(x,y)N

ランダムプロセスが(少なくとも)次数に定常的である場合、すべてのは同じ分布を持つため、平均が存在すると仮定すると、はすべて同じです。。同様に、 はすべてので同じであり、プロセスのパワーと呼ばれます。すべての物理プロセスには有限の力があるため、と仮定するのが一般的 です。この場合、特に古い工学文献では、このプロセスは2次プロセスと呼ばれます。名前の選択は、2次との混乱を招くため、残念です 1X(t)E[X(t)]=μtE[(X(t))2]tE [ X T 2 ] < E [ X T 2 ] T E [ X t 2 ]E[(X(t))2]<定常性stats.SEの私のこの回答を参照)、したがって、ここでは、すべてのに対してが有限であるプロセスを呼び出します(は定数です)有限パワープロセスとして、この混乱を避けます。しかし、もう一度注意してくださいE[(X(t))2]tE[(X(t))2]

1次定常プロセスは有限パワープロセスである必要はありません

次数定常的なランダムプロセスを考えます。現在、との共同分布はとの共同分布関数と同じであるため、および値はのみに依存します。これらの期待は有限パワープロセスに対して有限であり、その値はプロセスの自己相関関数と呼ばれます:は関数、時間ランダム変数と分離、および依存しない2X(t1)X(t1+τ)X(t2)X(t2+τ)E[X(t1)X(t1+τ)]=E[X(t2)X(t2+τ)]τRX(τ)=E[X(t)X(t+τ)]τX(t)X(t+τ)tE [ X T X T + τ ] = E [ X T + τ X t ] = E [ X t +まったく。注意また、その したがって自己相関関数は引数の偶数関数です。

E[X(t)X(t+τ)]=E[X(t+τ)X(t)]=E[X(t+τ)X(t+ττ)]=RX(τ),

有限パワーの2次定常ランダムプロセスには、次の特性があります。

  1. その平均は定数ですE[X(t)]
  2. その自己相関関数の関数で、ランダム変数の時間間隔および、としにまったく依存しません。RX(τ)=E[X(t)X(t+τ)]τX(t)X(t+τ)t

定常性の仮定により、ランダムプロセスの記述はある程度単純化されますが、実験データからモデルを構築することに関心のあるエンジニアや統計学者にとって、特に1つのサンプルパスのセグメントしかない場合、それらのすべてのCDFを推定するのは簡単な作業ではありません(または実現)測定を行うことができる。エンジニアがワークベンチ(またはソフトウェアライブラリ内のMATLAB / Python / Octave / C ++のプログラム)に必要な機器を既に持っているため、比較的簡単に実行できる2つの測定は、DC値 ofおよび自己相関関数x(t)1 1T0Tx(t)dtx(t)Rx(τ)=1T0Tx(t)x(t+τ)dt(またはそのフーリエ変換、パワースペクトル)。これらの測定値を、有限パワープロセスの平均および自己相関関数の推定値として使用すると、次に説明する非常に有用なモデルが得られます。x(t)


有限パワーランダムプロセスは、一定の平均値とその自己相関関数持つ場合、広義定常(WSS)プロセス(幸いにも同じ初期化WSSも持つ定常ランダムプロセス)と呼ばれます。は、時間差 (または)のみに依存します。RX(t1,t2)=E[X(t1)X(t2)]t1t2t2t1

定義は、プロセスを構成するランダム変数のCDFについて何も述べていないことに注意してください。それは完全にランダム変数の1次および2次モーメントに対する制約です。もちろん、有限パワーの2次定常(または次定常()または厳密に定常の)ランダムプロセス WSSプロセスですが、その逆は真である必要はありません。NthN>2

WSSプロセスは、任意の順序で静止している必要はありません。

例えば、ランダムなプロセスを考える ここで 4つの等しい可能性が高い値をとると。(怖がらないでください:このランダムプロセスの4つの可能なサンプルパスは、QPSK信号の4つの信号波形です)。各は離散確率変数であり、一般に、4つの同等の可能性のある値および、一般的におよび{X(t):X(t)=cos(t+Θ),<t<}Θ0,π/2,π3π/2X(t)cos(t),cos(t+π/2)=sin(t),cos(t+π)=cos(t)cos(t+3π/2)=sin(t)X(t)X(s)分布が異なるため、プロセスは一次定常でさえありません。一方、 間 、ごとに 要するに、プロセスの平均はゼロであり、その自己相関関数は時間差のみに依存するため、プロセス広義の定常です。しかし、それは一次定常ではないため、高次に対しても定常ではありません。

E[X(t)]=14cos(t)+14(sin(t))+14(cos(t))+14sin(t)=0
t
E[X(t)X(s)]=14[cos(t)cos(s)+(cos(t))(cos(s))+sin(t)sin(s)+(sin(t))(sin(s))]=12[cos(t)cos(s)+sin(t)sin(s)]=12cos(ts).
ts

2次定常(または厳密に定常)ランダムプロセスであるWSSプロセスで、ランダム変数の分布の特定の形式についてはほとんど言えません。要するに、

WSSプロセスは必ずしも(任意の順序で)静止しているわけではなく、WSSプロセスの平均および自己相関関数は、プロセスの完全な統計的記述を提供するのに十分ではありません

最後に、確率過程がガウス過程であると想定されると仮定します(合理的な信頼度でこれを「証明」するのは簡単な作業ではありません)。つまり、各に対して、はガウス確率変数であり、すべての正の整数および時刻、、に対して、 ランダム変数、、はガウス確率変数です。結合ガウス密度関数は完全にtX(t)n2nt1t2,tnNX(t1)X(t2),X(tn)ランダム変数の平均、分散、共分散によって決定され、この場合、平均関数知っている広義に必要な定数である必要はない-定常性)および自己相関関数すべての(広義の定常性に必要なのみに依存する必要はありません)プロセスの統計を完全に決定するには十分です。μX(t)=E[X(t)]RX(t1,t2)=E[X(t1)X(t2)]t1,t2t1t2

Gaussianプロセス WSSプロセスである場合、厳密に定常的な Gaussianプロセスでもあります。エンジニアとシグナルプロセッサにとって幸いなことに、多くの物理ノイズプロセスはWSSガウスプロセス(したがって厳密に定常的なプロセス)としてうまくモデル化できるため、自己相関関数の実験的観測はすべてのジョイント分布を容易に提供します。さらに、ガウスプロセスは線形システムを通過するときにガウス特性を保持するため、出力自己相関関数はとして入力自己相関関数に関連しています。

Ry=hh~RX
出力統計も容易に決定できるように、WSSプロセス一般および特にWSSガウスプロセスは、エンジニアリングアプリケーションで非常に重要です。


その意味で、「ホワイトノイズ」についてコメントしてください。定義により、での自己相関 は確率変数の分散です。AWGN(Additive White Gaussian Noise)の分散が無限であることを意味しますか?通常、人々はと書くので、それは間違っていますか?と書く必要がありますか?ありがとう。τ=0n(t) N(0,N0/2)n(t) N(0,δ(0)N0/2)
-Royi

1
@Drazick別の質問をしてください。
ディリップサルワテ

1
これは、定常プロセスの定義における素晴らしいミニコースです。私はそれのようなものを見たことがない-とても系統的かつ明確にレイアウトされています。コミュニティWiki?
-abalter

@Dilip Sarwate失礼しました。例では。すべてのtでE [X(t)] = 0なのはなぜですか?エルゴディシティを想定しましたか?期待値を計算するために、シータの確率密度関数からX(t)の確率密度関数をどのように導きましたか?E [X(t)X(s)] = E [cos(t + theta)* cos(s + theta)]でしょ?この表現を単純化し、書いたものにたどり着くためにどのステップを取りましたか?ありがとう
-VMMF

1
@VMMFありませんNO使用エルゴード性が。は離散確率変数です。なぜなら、は離散確率変数であり、等しい確率で値およびを取るからです。。エルゴ、。 の値をとる(、及びに等しい確率で。したがって、X(t)=cos(t+Θ)Θ±cos(t)±sin(t)14E[X(t)]=0X(t)X(s)cos(t)cos(s)(cos(t))(cos(s))=cos(t)cos(s)sin(t)sin(s)1(sin(t))(sin(s))=sin(t)sin(s) E[XtXs]=114E[X(t)(X(s)]=12(cos(t)cos(s)+sin(t)sin(s))=12cos(ts)。したがって、
ディリップ・サーワテ
弊社のサイトを使用することにより、あなたは弊社のクッキーポリシーおよびプライバシーポリシーを読み、理解したものとみなされます。
Licensed under cc by-sa 3.0 with attribution required.