このページの内容は最新ではありません。最新版の英語を参照するには、ここをクリックします。
サンプルの自己相関の信頼区間
この例では、ホワイト ノイズ過程の自己相関列の信頼区間を作成する方法を示します。長さ サンプルのホワイト ノイズ過程の実現を作成します。サンプルの自己相関をラグ 20 まで計算します。ホワイト ノイズ過程に対し 95% 近似の信頼区間でサンプルの自己相関をプロットします。
ホワイト ノイズ ランダム ベクトルを作成します。再現性のある結果を得るために、乱数発生器を既定の状態に設定します。正規化された自己相関サンプルを 20 ラグまで求めます。
rng default L = 1000; x = randn(L,1); [xc,lags] = xcorr(x,20,'coeff');
標準偏差が である 正規分布について、95% 信頼限界の下限および上限を求めます。95% 信頼区間に対する棄却限界値は で、信頼区間は以下となります。
vcrit = sqrt(2)*erfinv(0.95)
vcrit = 1.9600
lconf = -vcrit/sqrt(L); upconf = vcrit/sqrt(L);
95% 信頼区間でサンプルの自己相関をプロットします。
stem(lags,xc,'filled') hold on plot(lags,[lconf;upconf]*ones(size(lags)),'r') hold off ylim([lconf-0.03 1.05]) title('Sample Autocorrelation with 95% Confidence Intervals')
上の Figure から、ホワイト ノイズ過程に対して予想したとおり、自己相関の値が 95% 信頼区間の外側にあるのはラグ 0 だけであることがわかります。この結果に基づき、このデータがホワイト ノイズ過程の実現であることがわかります。