%?通過求互相關序列，測定汽車速度。
ndelay=20;
nx=[-30:30];x=sin（abs（nx））;
subplot（221）plot（nx?x）;?grid;?%?顯示序列?x
xlabel（‘n‘）;?ylabel（‘x（n）‘）;?title（‘ori-wave‘）;
%===============================================
ny=nx;
if?ndelay>0?%右移
????y=[x（（length（x）-ndelay+1）:length（x））x（1:（length（x）-ndelay））];
elseif?ndelay<0
????????y=0
????elseif?ndelay==0
????????????y=x
????????end
????end
end
subplot（222）plot（ny?y）;?grid;?%?顯示序列?y
xlabel（‘n‘）;?ylabel（‘y（n）‘）;?title（‘splited‘）;

%===============================================
[rxynrxy]=xcorr（yx）;
[xsa]=xcorr（xx）;[ysb]=xcorr（yy）;
rxy=rxy/sqrt（max（xs）*max（ys））;??%?歸一化的互相關序列？
subplot（223）plot（nrxyrxy）;?grid;?
xlabel（‘n‘）;?ylabel（‘rxy（n）‘）;?title（‘Nomalized?CC‘）;

%===============================================
subplot（224）H?=?stem（nrxyrxy）;?
set（H?‘markersize‘?2）;?grid;?%?在互相關序列的最大點附近顯示該序列
%axis（[ndelay-10ndelay+10min（rxy）max（rxy）]）;set（gcf‘color‘‘w‘）?%?置圖形背景色為白
%line（[-48][00]）
title（‘NCC-details‘）
xlabel（‘ny‘）ylabel（‘rxy‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件????????1090??2014-03-28?08:48??CC.m