%互相關(guān)函數(shù)
clear;
load a;
load b;
load c;
load a;
load b;
load c;
a=a';
b=b';
a=a';
b=b';
a=a';
b=b';
a=a(1,1:10000);
b=b(1,1:10000);
a=a(1,1:10000);
b=b(1,1:10000);
a=a(1,1:10000);
b=b(1,1:10000);
N=1000;
n=0:N-1;
Fs=500;
Lag=100;
[c,lags]=xcorr(a,Lag,'unbiased');
[d,lags]=xcorr(b,Lag,'unbiased');
[c,lags]=xcorr(a,Lag,'unbiased');
[d,lags]=xcorr(b,Lag,'unbiased');
[c,lags]=xcorr(a,Lag,'unbiased');
[d,lags]=xcorr(b,Lag,'unbiased');
figure(1);
hold on;
plot(lags/Fs,c,'r');
plot(lags/Fs,d,'b');
plot(lags/Fs,c,'r');
plot(lags/Fs,d,'b');
plot(lags/Fs,c,'r');
plot(lags/Fs,d,'b');
xlabel('t');
ylabel('y');
title('compare correlation');
grid;