Processing math: 100%

3. 雑音を含む周期性信号の場合

f[i] は定常過程と仮定しましたので(標本)自己相関関数 R[n] の値は時刻 i によらずに遅延時刻 n だけで決まります。
そして R[0]R[0] から R[N1] の中で常に最大で、全く同じ信号間の自己相関を表します。
また R[n] の値がある程度大きい値である場合は f[i]f[i+n] の間に強い相関がみとめられます。

以上から次の結論が導かれます。

「ある n に対して R[n]R[2n]R[3n]、・・・の値が R[0] に近ければ、信号 f[i] は周期 n周期性信号である可能性が高い」

ちなみに通常 f[i] には雑音が含まれるので各 R[in]R[0] の値が完全に一致することは稀なのですが、もしこれらの値が一致した場合は信号は周期 n の周期性信号で、かつ雑音が全く含まれていないことを意味します。
例えば前ページで雑音を含んでいない周期 5 のサイン波の例では R[0]=R[5]=1/2 となっていました。

そこでここでは例として周期 5 のディジタル信号

f[i]=3+1sin(2π(51)i)+2cos(2π(53)i)

に雑音として平均 0、分散 1 の正規乱数を各時刻で付加すると自己相関関数がどう変化するかみてみましょう。
なお信号長は N=10000 とします。

まず、雑音を加える前の信号は以下の図1です。
確かに周期 5 の信号となっている事が分かります。

図1: 雑音を加える前の信号 (最初の 20 点だけ表示)

 

そして雑音を加えて実際に得られた信号が図2です。

図2: 図1の信号に正規乱数 N(0,1) を加えた信号 (最初の 20 点だけ表示)

 

この信号の自己相関関数のグラフを図3に示します。

図3: 図2の信号の自己相関関数 (最初の 20 点だけ表示)

 

この様に R[5]R[10]R[15] の値が R[0] に近く、他の自己相関関数値は相対的に小さかったので、雑音を加える前の信号が周期 5 の周期性信号であることが分かります。