% ENM 4.6
% Denna kod skrevs på övningen 2017-09-21 i SF1545

clear

% Datapunkter
t = (0:2:10)';
y = [1 1.6 1.4 0.6 0.2 0.8]';

% MKV
A = [ones(size(t)), sin(2*pi*t/12), cos(2*pi*t/12)];
c = A\y

% Plotta
clf % rensa figuren
plot(t,y,'o') % datapunkter
hold on
tfin = linspace(0,10);
yfin = c(1) + c(2)*sin(2*pi*tfin/12) + c(3)*cos(2*pi*tfin/12);
plot(tfin,yfin)

% Räkna ut residualvektorn själva om ni vill!
