% autocorr.m

%

% ＜　Autocorrelation Function ＞

%

Fs = 200; % サンプリング周波数 (Hz)

T = 1; % 観測時間（1sec）

dt = 1/Fs; % 時間分解能（サンプリング間隔）

N = T/dt; % サンプル数

t = (0:dt:T-dt); % time vector

%

% ----------------------signal---------------------------

rn = 0.7*randn(size(t)); % 正規乱数

x = sin(2*pi*4*t) + rn; % 信号（時系列データ）

%

figure(1)

subplot(311); plot(x); % 観測信号のプロット

axis([0,2*N,-4,4]); legend('signal');

% ----------------------autocorrelation------------------

s = (-N+1:1:N-1); % time-lag vector

w1 = xcorr(x,x,'unbiased');

subplot(312); plot(s,w1); % 自己共分散関数（不偏）

legend('autocorrelation (unbiased)');

w2 = xcorr(x,x,'biased');

subplot(313); plot(s,w2); % 自己共分散関数（バイアスあり）

legend('autocorrelation (biased)');

% -------------------------------------------------------

% 2002.4.28 by K.Tsukada