% Tund 11, ülesanne 1
x=[30,33,34,35,39,44,45];
y=[1.11,1.21,1.23,1.25,1.3,1.4,1.42];
kaalud=[1,1,1,1.5,1.5,2,2];
plot(x,y,'o')
hold on
% 1) lähendame andmeid lineaarfunktsiooniga
% a) vähimruutude mõttes kaalusid arvestades
% Phi(x)=c1*x+c0, c1,c0=?
% AC=B
A=[sum(kaalud.*x.^2) sum(kaalud.*x.^1)
    sum(kaalud.*x.^1) sum(kaalud.*x.^0)];
B=[sum(kaalud.*y.*x.^1)
    sum(kaalud.*y.*x.^0)];
C=inv(A)*B % C=[c1,c0]
% Phi(x)=0.0187*x+0.5810
y1kaalud=polyval(C,x);
plot(x,y1kaalud,'m-')
hold on
% b) lähendame andmeid lineaarsplainiga
y1splain=interp1(x,y,x,'linear');
plot(x,y1splain,'b-')
hold off
% 2) tööstusmasina täpsus, kui t=50 h
% x=50
polyval(C,50)
interp1(x,y,50,'linear') % ei saa leida, kuna 50 jääb ette antud piirkonnast välja