% exa060104_xcorr_xcov.m, for example 6.1.4  
% to test xcorr.m and xcov.m ;

clear;
N=256;
f=.1;
a1=5;
a2=3;
Mlag=N/4;

% 產生兩個正弦加白噪聲的數據；
x=a1*sin(2*pi*f*(0:N-1))+2*randn(1,N);
y=a2*sin(2*pi*f*(0:N-1))+randn(1,N);
subplot(2,2,1);
plot(x(1:N/2));
grid on;
title('原始信號1')
subplot(2,2,2);
plot(y(1:N/2));
grid on;
title('原始信號2')

% 求這兩個數據向量的互相關函數；
corrxy=xcorr(x,y,Mlag,'unbiased');
subplot(2,2,3);
plot((-Mlag:1:Mlag),corrxy);
grid on;
title('互相關函數') 

% 求這兩個數據向量的互協(xié)方差函數；
covxy=xcov(x,y,Mlag,'biased');
subplot(2,2,4);
plot((-Mlag:1:Mlag),covxy); 
title('互協(xié)方差函數')