% esempio di approssimazione mediante retta dei minimi quadrati
close all
clear all
clc

% dati lungo x
x = [1,2,3,4,5];

% dati lungo y
y = [1,2.5,2.5,4.5,4.5];

% costruzione della interpolante polinomiale
p = polyfit(x,y,1);


% visualizzazione dei dati
xpoly = linspace(1,5,1000); % punti su cui valutare il polinomio
ypoly = polyval(p,xpoly);   % valutazione del polinomio su un insieme di punti

% disegno polinomio e i sui punti dove ho interpolato
plot(xpoly,ypoly,'b-');
hold on;
plot(x,y,'o');
hold off;