clear;
close all;

taille_ecran = get(0,'ScreenSize');
L = taille_ecran(3);
H = taille_ecran(4);

% Fenetre d'affichage :
figure('Name','Points situes au voisinage d''un cercle', ...
	   'Position',[0.4*L,0.05*H,0.6*L,0.7*H]);
axis equal;
hold on;
set(gca,'FontSize',20);
hx = xlabel('$x$','FontSize',30);
set(hx,'Interpreter','Latex');
hy = ylabel('$y$','FontSize',30);
set(hy,'Interpreter','Latex');

% Limites de l'affichage des donnees centrees en (0,0) :
taille = 20;
limites_affichage = [-taille taille -taille taille];

% Creation du cercle et des donnees bruitees :
n = 50;
sigma = 0.5;
[x_cercle,y_cercle,theta_cercle,x_donnees_bruitees,y_donnees_bruitees,theta_donnees_bruitees] ...
		= creation_cercle_et_donnees_bruitees(taille,n,sigma);

% Affichage du cercle :
plot(x_cercle([1:end 1]),y_cercle([1:end 1]),'r','LineWidth',3);

% Donnees aberrantes :
proportion = 0.25;
n_donnees_aberrantes = floor(proportion*n);
x_donnees_bruitees = [x_donnees_bruitees taille*(2*rand(1,n_donnees_aberrantes)-1)];
y_donnees_bruitees = [y_donnees_bruitees taille*(2*rand(1,n_donnees_aberrantes)-1)];

% Affichage des donnees bruitees :
plot(x_donnees_bruitees,y_donnees_bruitees,'k+','MarkerSize',10,'LineWidth',2);
axis(limites_affichage);
lg = legend(' Cercle', ...
		' Donnees (bruitees + aberrantes)', ...
		'Location','Best');
grid on;