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

matlab,功率谱估计

%两个周期信号的互相关函数
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