6bijiao correlation.m

来自「信号自相关进行比较」· M 代码 · 共 42 行

%互相关函数
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;