%例程5－5  計算隨機信號的自相關函數
% e.g.5-5.m for example5-5;
% to test function of xcorr;
 
clear all;
N=10^3;
n=[0:N-1];
Fs=500;
t=n/Fs;
Mlag=100;
% Generate signal plus noise
x1=sin(2*pi*10*t)+0.7*randn(1,length(t));
figure(1);  plot(t,x1);
xlabel('n'); ylabel('x1(n)');
title('Original signal x1');
grid on;
% Computing the unbiased atuocorretaion of x1
[r1,lags]=xcorr(x1,Mlag,'unbiased');
figure(2); plot(lags/Fs,r1);
xlabel('n'); ylabel('Rx1(n)');
title('Autocorrelation of x1');
grid on;