clear all
% chiudiamo tutti i grafici
close all
n=20;
aoriginale=[-3,1,0];
x=linspace(0,4*pi,n);
y=funzione(aoriginale,x);
% perturbo artificialmente i dati
y=y+(rand(size(x))-0.5)/5;
plot(x,y,'*')
A=[sin(x'),cos(x'),ones(n,1)];
b=y';
% risolvo il sistema
a=A\b;
xx=linspace(0,4*pi,100);
yy=funzione(a,xx);
hold on
plot(xx,yy)
disp('errore minimi quadrati')
norm(y-funzione(a,x))
disp('errore funzione originale')
norm(y-funzione(aoriginale,x))