function r=corr(seq1,seq2)
if nargin==1    %輸入參數的個數
    seq2=seq1;
end
N=length(seq1);
for k=-N+1:-1
    seq2_shift=[seq2(k+N+1:N) seq2(1:k+N)];
    r(N+k)=seq1*seq2_shift';
end
for k=0:N-1
    seq2_shift=[seq2(k+1:N) seq2(1:k)];
    r(N+k)=seq1*seq2_shift';
end