% EX 4.6

% Tidvattenstånd i Nordsjön ungefär 12-h-cykler
% h(t) = h0 + a1*sin(2*pi*t/12) + a2*cos(2*pi*t*12)

% Mätningar
t = (0:2:10)';
h = [1.0 1.6 1.4 0.6 0.2 0.8]';

% Minstakvadratanpassa, hitta c = [h0 a1 a2]' som bäst passar A*c = h
A = [ones(length(t), 1) sin(2*pi*t/12) cos(2*pi*t/12)];
c = A\h;

tt = 0:0.1:12;
hh = c(1) + c(2)*sin(2*pi*tt/12) + c(3)*cos(2*pi*tt/12);

figure(1)
clf
plot(tt, hh, 'r')
hold on
plot(t, h, 'k*')
grid on