google = matlab aproximace metoda nejmensich ctvercu
http://uprt.vscht.cz/majerova/matlab/lekce9.html
>> u=[1 1.5 2.1 2.5 3 3.1 3.2 3.5]; % namerena data, nezav. prom.
>> v=[7.8 8.15 8.3 8.25 8.1 8.3 8.35 8.2]; % namerena data, zav. prom.
>> p=polyfit(u,v,2) % koeficienty polynomu 2. stupne pro 'u' a 'v'
p =
-0.1684 0.8977 7.1150
>> v_aprox=polyval(p,u); % hodnoty polynomu v 'u'
>> S=sum((v_aprox-v).^2) % soucet ctvercu odchylek
S =
0.0345
>> plot(u,v,'r+') % graf, puvodni data jako cervene krizky
>> hold on % prikreslime dalsi
>> plot(u,v_aprox,'k.-') % graf polynomu jako cerne body spojene carou
>> axis([0.8 3.7 7.7 8.4]) % uprava os
>> title('Data a jejich aproximace parabolou')% nazev grafu
>> legend('data','aproximacni polynom',4) % zobrazime legendu (4 = vpravo dole)
https://files.klaska.net/sites/files.klaska.net/files/manual_files/cvut/matlab/z%20hodin/matlab08.pdf
% připravíme data: fční hodnoty daných fcí v bodech xX = [f1(x), f2(x), f3(x), f4(x)];koef = X\y; % metoda nejmenších čtvercůa=koef(1),b=koef(2),c=koef(3),d=koef(4)
y_aprox = X*koef; % hledanéřešení
---
x = [0 0.2 0.4 0.5];%data
y = [2 3 4 5]; %data
x1 = 0:.1:.5;
spl1 = interp1(x,y,x1, 'spline');
spl2 = spline(x,y,x1);
plot(x,y,'ko', x1,spl1,'g', x1,spl2,'r')
V matlabu nedelam, ale minimalne by to mohl byt dobry zacatek programu a po par upravach by to asi mohlo fungovat.