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