clear
% Harjutustund nr. 9, ülesanne 2
% andmed
x=[2014:2021];
y=[22.2,19.5,18.0,17.9,25.2,18.6,18.7,29.3];
% joonis
plot(x,y,'o')
hold on
grid on
% a) lähendame andmeid lineaarfunktsiooniga
k1=polyfit(x,y,1)
yk1=polyval(k1,x)
disp([x;y;yk1])
% esitame joonisel leitud sirge
plot(x,yk1,'r')
% lähendame andmeid lineaarsplainiga S^(1,0) ehk igas osalõigus sirgega
ylineaarsplain=interp1(x,y,x,'linear');
plot(x,ylineaarsplain,'b')
hold on
% b) lähendame andmeid neljanda astme polünoomiga
k4=polyfit(x,y,4);
xvork=[2014:0.1:2021];
yk4=polyval(k4,xvork);
plot(xvork,yk4,'k')
hold off
% c) prognoosime
polyval(k1,2023) % lineaarse sirgega
interp1(x,y,2023,'linear') % lineaarsplainiga
polyval(k4,2023) % neljanda astme polünoomiga