abc Matlab - elektronski priručnik
III dio Matematika u Matlabu - 3.3. Interpolacija i ekstrapolacija

3.3.3. Interpolacija sa polyfit i polyval

Moguće je upotrijebiti Matlabove polinomske naredbe da bi se izgradio interpolacioni polinom. Ovdje je primjer kako se one koriste da bi se našlo fitovanje polinomom 5-og reda za grubo predstavljanje sinusne funkcije.
clear; close all;

% napravi grubi skup podataka sa dx suvise velikim za dobru tacnost

dx=pi/5;
x=0:dx:2*pi;

y=sin(x);

% napravi fit polinomom 5-og reda za ove podatke

p=polyfit(x,y,5);

% napravi finu x-resetku

xi=0:dx/20:2*pi;

% izracunaj fiting polinom na finoj resetki

yi=polyval(p,xi);

% i prikazi fit, podatke i egzaktnu sinusnu funkciju

plot(x,y,'b*',xi,yi,'r-',xi,sin(xi),'c-')
legend('Podaci','Fit','Egzaktna sinusna funkcija')

% prikazi razliku izmedju polinomskog fita i egzaktne sinusne funkcije

figure
plot(xi,yi-sin(xi),'b-')
title('Error in fit')

Kvadratna interpolacija i ekstrapolacija    <    Index    >    Matlab-ove naredbe Interp1 i Interp2