ccorr.m

来自「各种调制仿真程序M语言源代码」· M 代码 · 共 13 行

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