function R=length_normalized_cor(x)
    R=zeros(length(x));
    helper=x.*flipdim(flipdim(conj(x),2),1);
    for k=1:length(x)
        R(k)=sum(helper(1:(length(x)-k+1)))/length(x)
    end
endfunction

function R=k_normalized_cor(x)
    R=zeros(length(x));
    helper=x.*flipdim(flipdim(conj(x),2),1);
    for k=1:length(x)
        R(k)=sum(helper(1:(length(x)-k+1)))/(length(x)-k+1)
    end
endfunction