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
|