%Samp9_4
clf
N=1000;n=0:N-1;Fs=500;t=n/Fs; %數據個數采樣頻率和時間序列
Lag=200;    %最大延遲單位數
x1=90*sinc(pi*(n-0.1*Fs));  %第一個原始信號，延遲0.1s
y1=50*sinc(pi*(n-0.3*Fs));  %第二個原始信號,延遲0.3s
[c,lags]=xcorr(x1,y1,Lag,'unbiased');   %計算兩個函數的互相關
subplot(2,1,1),plot(t,x1,'r');    %繪第一個信號
hold on;plot(t,y1,'b:');   %在同一幅圖中繪第二個信號
legend('信號x', '信號y');  %繪制圖例
xlabel('時間/s');ylabel('x(t)  y(t)');
title('信號x和y');hold off
subplot(2,1,2),plot(lags/Fs,c,'r'); %繪制互相關信號
xlabel('時間/s');ylabel('Rxy(t)');
title('信號x和y的相關');