Лабораторная работа 6

Лабораторная работа №6 «Интерполирование алгебраическими многочленами»

>> x=[0 0.5 1 1.5 2]

>> y=[1.4346 0.3421 0.0674 0.0601 0.0586]

>> xzv=0.83

>> P1=polyfit(x(2:3),y(2:3),1)

P1 = -0.5494 0.6168

>> P2=polyfit(x(2:4),y(2:4),2)

P2 =0.5348 -1.3516 0.8842

>> P3=polyfit(x(1:4),y(1:4),3)

P3 =-0.7339 2.7364 -3.3697 1.4346

>> z1=polyval(P1,xzv)

z1 =0.1608

>> z2=polyval(P2,xzv)

z2 =0.1308

>> z3=polyval(P3,xzv)

z3 = 0.1032

>> xi1=x(2):0.05:x(3);

>> xi2=x(2):0.05:x(4);

>> xi3=x(1):0.05:x(4);

>> y1=polyval(P1,xi1);

>> y2=polyval(P2,xi2);

>> y3=polyval(P3,xi3);

>> plot(x,y,'*k',xi1,y1,xi2,y2,xi3,y3);grid

>> P4=polyfit(x(1:5),y(1:5),4)

P4 =

0.1925 -1.3115 3.2659 -3.5141 1.4346

>> z4=polyval(P4,xzv)

z4 =0.1092

>> eps1=abs(z1-z2)

eps1 =

0.0300

>> eps2=abs(z2-z3)

eps2 =

0.0276

>> eps3=abs(z3-z4)

eps3 =

0.0060