%
% File per la verifica dei risultati
% e discussione
%

% Funzione teorica

u = linspace(0,1,100);

y = sin(2*pi*u).*cos(4*pi*u).^2;

figure, plot(u,y,'r')
xlabel('Ingresso indipendente')
ylabel('y*')
title('Funzione incognita da approssimare')

% Generazione dati

[Psims,Is]=sort(Psim);
Tsims=Tsim(Is);

figure, plot(Psims,Tsims,'b')
xlabel('Ingresso indipendente')
ylabel('y_n')
title('Campioni della funzione incognita con rumore')

figure, plot(Psims,Tsims,'b--',u,y,'r-')
xlabel('Ingresso indipendente')
ylabel('y_n')
title('Confronto funzione incognita e dati con rumore')

[Pvals,Iv]=sort(Pval);
Tvals=Tval(Iv);


figure, plot(Pvals,Tvals,'b--',u,y,'r-')
xlabel('Ingresso indipendente')
ylabel('y_n')
title('Confronto funzione approssimata e funzione teorica')