load('adaptFilter19taps.mat', 'rAdaptEq');
%% convergence (|error| < 1e-3) at ~ 6000th symbol
symPerFrame = 500;
for i = 1 : symPerFrame : 6000
  subset = rAdaptEq(i : i + symPerFrame - 1);
  h = scatterplot(subset);
  title({'Adaptive filter (CMA) with 19 taps', ...
         strcat(['Symbols ', num2str(i), ' to ', ...
                 num2str(i + symPerFrame - 1)])}, 'Interpreter', 'Latex');
  xlabel('In-Phase', 'Interpreter', 'Latex');
  ylabel('Quadrature', 'Interpreter', 'Latex');
  axis([-1 1 -1 1]);
  set(gcf, ...
      'Units', 'centimeters', ...
      'Position', [0 0 15 15], ...
      'PaperPositionMode', 'auto');
  formatFigure;
  set(gca, ...
      'Units', 'normalized', ...
      'YTick', -1:0.5:1, ...
      'XTick', -1:0.5:1);
  %%formatFigure;

  print(strcat('adaptEqAni_', num2str(ceil(i / symPerFrame))), '-depsc2');
  close(h);
end