tmp604.m

来自「信号分析与处理的一些源代码」· M 代码 · 共 28 行

%MATLAB PROGRAM 6-4
clf
N=1000;
n=0:N-1;
Fs=500;
t=n/Fs;
Lag=200;

%signal x1 y1
%============
x1=90*sinc(pi*(n-0.1*Fs));
y1=50*sinc(pi*(n-0.3*Fs));
[c,lags]=xcorr(x1,y1,Lag,'unbiased');
subplot(211)
plot(t,x1,'r')
hold on
plot(t,y1,'b');
xlabel('t');ylabel('x(t) y(t)');
title('Original signal ');
grid
hold off

subplot(212)
plot(lags/Fs,c,'r');
xlabel('t');ylabel('Rxy(t)');
title('Correlation ');
grid