Matlab Hogyan - Hidrodinamikai Rendszerek Tanszék

8. fejezet
Görbeillesztés
Az alap Matlab csomagban csak a polinomiális görbeillesztés található meg előre kész eljárás formájában. A
p = polyfit(x,y,n)
parancs eredménye p vektor, mely az n-edfokú polinom együtthatóit tartalmazza csökkenű kitevő szerint.
Ha egyszer ezeket az együtthatókat kiszámítottuk, tetszőleges xi pontokban
yi = polyval(p,xi)
paranccsal számíthatjuk ki a függvény értékét.
A polyfit pelda.m program az f(x) = x 2 sin(5x) függvényhez illeszt poinomot a (0, π/2) intervallumon,
majd kiszámítja a kapott polinom eltérését az osztáspontokban a norm parancs segítségével.
polyfit pelda.m
function polyfit_pelda
X=0:0.01:pi/2; Y=sin(5*X).*X.^2;
for i=1:1:10
p=polyfit(X,Y,i);
yy=polyval(p,X);
hiba=norm(Y-yy,2);
fprintf(’\nfokszam: %2d, hiba=%5.3e’,i,hiba);
end
>> polyfit_pelda
fokszam: 1, hiba=9.571e+00
fokszam: 2, hiba=6.386e+00
fokszam: 3, hiba=2.374e+00
fokszam: 4, hiba=2.341e+00
fokszam: 5, hiba=8.731e-01
fokszam: 6, hiba=2.577e-01
fokszam: 7, hiba=1.375e-01
fokszam: 8, hiba=1.321e-02
fokszam: 9, hiba=1.030e-02
fokszam: 10, hiba=4.593e-04
29