ccorr.m

关于通信基本原理的仿真

function r = ccorr(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