%相關運算函數
%輸入：序列向量，當序列個數為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