转载自 http://blog.163.com/to_be_myself/blog/static/176060227201101762159227/
Matlab中用于计算自相关函数的指令是xcorr.比如矩阵A=[1 2 3];
xcorr(A)=3.0000 8.0000 14.0000 8.0000 3.0000
自相关函数是信号间隔的函数,间隔有正负间隔,所以n个长度的信号,有2n-1个自相关函数值,分别描述的是不同信号间隔的相似程度。
比如,上面的矩阵,最后得到5个结果,其中第三个是自己和自己相乘,最后相加的结果,值最大1*1+2*2+3*3=14。而第二个和第四个分别是间隔正负1的结果也就是1*2+2*3=8,2*1+3*2=8。第1个和第五个分别是间隔正负2,也就是1*3=3,3*1=3。
xcorr求出的结果仅仅是x(n)*x(n+m)并对其求和,并没有除以前面的N或者是N-|K|。不用这个函数也可以求
for k=0:1:p
t5=0;
for n=0:1:N-k-1
t5=t5+conj(x(n+1))*x(n+1+k);
end
Rxx(k+1)=t5/N;
end
也可以实现,其中N为序列长度