两正弦信号的自相关函数.m

来自「matlab,功率谱估计」· M 代码 · 共 24 行

%两个周期信号的互相关函数
clf;
N=1000;
n=0:N-1;
Fs=500;
t=n/Fs;
Lag=200;
x=sin(2*pi*10*t);
y=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,'r')
xlabel('t');ylabel('Rxy(t)')
title('correlation')
grid