% exercice 8 du TD 1 du Master IAC 2010
% tronc commun ASO (Apprentissage Statistique et Optimisation)

% a appeler depuis la ligne de commande, apres avoir fait tourner compareInversion
function [err1 err2] =  regresseInversion

indices = 100:100:800;
[t1 t2] = compareInversion(indices);

for degre=2:6
[p,ErrorEst] = polyfit(indices, t1,degre);
approx1 = polyval(p,indices);
[p,ErrorEst] = polyfit(indices, t2,degre);
approx2 = polyval(p,indices);

disp(['Erreurs pour le degre ' num2str(degre)]);
err1 = sumsq(approx1-t1)
err2 = sumsq(approx2-t2)
% on les trace
plot(indices, t1,'+', indices, t2,'ro', indices, approx1,'-',  indices, approx2,'r-');
disp ('appuyez sur une touche');
pause
end