abc Matlab - elektronski priručnik
III dio Matematika u Matlabu - 3.2. Polinomi

3.2.7. Fitovanje podataka u polinom

Ako imamo neke podatke u formi niza (x,y), Matlab može izvršiti tzv. fitovanje najmanjih kvadrata polinoma bilo kog reda koji se izabere za te podatke. U sljedećem primjeru neka podaci budu sinusna funkcija između 0 i π u koju će biti fitovan polinom reda 4. Nakon toga ćemo nacrtati te dve funkcije u istom okviru da vidimo da li je fitovanje bilo dobro.
clear; close all;

x=linspace(0,pi,50);

% napravi sinusnu funkciju sa 1% slucajne greske

f=sin(x)+.01*rand(1,length(x));

% fitovanje po podacima

p=polyfit(x,f,4);

% evaluacija fitovanja

g=polyval(p,x);

% nacrtaj fitovanje i podatke zajedno

plot(x,f,'r*',x,g,'b-')

Evaluacija polinoma    <    Index    >    Interpolacija i ekstrapolacija