correlate.m

输入为一个序列或两个序列

%相关运算函数
%输入：序列向量，当序列个数为1时，为自相关运算；当序列个数为2时，为互相关运算
%输出：相关运算结果的向量
function r=correlate(seq1,seq2)
if nargin==1
    seq2=seq1;
end
for k=-length(seq1)+1:-1
    seq2_shift=[seq2(k+length(seq1)+1:length(seq1)) seq2(1:k+length(seq1))];
    r(length(seq1)+k)=seq1*seq2_shift'/length(seq1);
end
for k=0:length(seq1)-1
    seq2_shift=[seq2(k+1:length(seq1)) seq2(1:k)];
    r(length(seq1)+k)=seq1*seq2_shift'/length(seq1);
end