Saso@1: function rl = reallocate(mag, ests, N_rs, fs)
smusevic@2:   frqs = imag(ests(end,:))/2/pi;
Saso@1:   %frqs_rs = [0:N_rs-1]*fs/N_rs;
Saso@1:   bins_rs = round(frqs * N_rs/fs)-1;
smusevic@2:   bins_rs = bins_rs(bins_rs > 0 & bins_rs <= N_rs);
Saso@1:   rl = zeros(N_rs,1);
Saso@1:   % can this be done without for loop?
Saso@1:   for k=1:length(bins_rs)
Saso@1:     rl(bins_rs(k)) = rl(bins_rs(k)) + mag(k);
Saso@1:   end
Saso@1: end