function corr = Mcorr(sig1, sig2, length_sig, length_c)
if length_c > length_sig
    corr=0;
    return;
end
for k = 1:length_c
    corr(k) = 0;
    for m = 1:length_sig - k
        corr(k) = corr(k) + sig1(m) * sig2(m+k-1);
    end
    corr(k) = corr(k)/length_sig; %normirovka
end