Лабораторные / Lab3

x=[0; 2; 4; 6; 10];
y=[5; 10.5; 12.6; 15.9; 18.5];
plot(x,y,['O','r']);
grid on;
n=5;

p1=polyfit(x, y, 1);
f1=p1(1)*x+ p1(2);
plot(x,y,'ob',x,f1,'k');
hold on;
Y1=polyval(p1,x);

p2=polyfit(x, y, 2);
f2=p2(1)*x.^2+p2(2)*x+p2(3);
plot(x, f2);
hold on
Y2=polyval(p2,x);

p3=polyfit(x, y, 3);
f3=p3(1)*x.^3+p3(2)*x.^2+p3(3)*x+p3(4);
plot(x, f3, 'r');
hold on
Y3=polyval(p3,x);

p4=polyfit(x, y, 4);
f4=p4(1)*x.^4+p4(2)*x.^3+p4(3)*x.^2+p4(4)*x+p4(5);
plot(x, f4, 'gr');
hold on
Y4=polyval(p4,x);

pogr1=Y1-y;
pogr2=Y2-y;
pogr3=Y3-y;
pogr4=Y4-y;

s1=sqrt(1/(1+n)*sum((Y1-y).^2));
s2=sqrt(1/(1+n)*sum((Y2-y).^2));
s3=sqrt(1/(1+n)*sum((Y3-y).^2));
s4=sqrt(1/(1+n)*sum((Y4-y).^2));