clear;
N=10000;
n=0:N-1;
Fs=500;
t=n/Fs;
Lag=20000;%最大延遲單位數
x1=sin（pi*2*0.1*（t-0.1））;
x2=sin（pi*2*0.1*（t-2.4））;
%x1=90*sinc（pi*（n-0.1*Fs））;%第一個原始信號，延遲0.1s
%x2=50*sinc（pi*（n-0.3*Fs））;%第二個原始信號，延遲0.3s
[clags]=xcorr（x1x2Lag‘unbiased‘）;%計算兩個函數的互相關


subplot（2