steplen = 5;
figure
hold on
voronoi(real(constellation), imag(constellation))
xlim([-1.5 1.5])
ylim([-1.5 1.5])
axis square

thing_to_plot = MF_output_shifted.*exp(1i*-phzError);
thing_to_plot = thing_to_plot./( sum(abs(samples))./length(samples)); 
thing_to_plot = thing_to_plot.*(1./avg_preamble_energy);

for j = (peak_idx+1)*sps:steplen:length(thing_to_plot)
    
    plot(real(thing_to_plot(j)), imag(thing_to_plot(j)), 'ob')
    pause(0.1)
end