tmp603.m

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

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

%signal x and y
%========
x=sin(2*pi*10*t);  
y=0.5*sin(2*pi*10*t+90*pi/180);
[c,lags]=xcorr(x,y,Lag,'unbiased');

subplot(211)
plot(t,x,'r')
hold on
plot(t,y,'b')
xlabel('t');ylabel('x(t) y(t)');
title('Original signal ');
grid
hold off

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