% ENM 4.6
% Kod till övningen 2017-09-28 i SF1514/18/19

clear

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

clf
plot(t, y, '*') % rita datapunkter
hold on

% Minstakvadrat-anpassning
A = [ones(size(t)), sin(2*pi*t/12), cos(2*pi*t/12)]

c = A\y;
h0 = c(1)
a1 = c(2)
a2 = c(3)

% Rita ut kurvan
T = linspace(0, 10);
H = h0 + a1*sin(2*pi*T/12) + a2*cos(2*pi*T/12);
plot(T, H)

% (Man kan även räkna ut residualvektorn om man vill:)
r = y - A*c
