% Eksempel -- Lineaer modell harmonisk svingning

m=4;        % antall punkter


x=[0; 1; 2; 3]; y=[-3; 3; 5; -2]; %x og y koordinatene til de gitte punktene

figure(11)

plot(x, y,'or') %plotter de gitte punktene



X=[ones(m,1) sin((0.5)*pi*x) cos((0.5)*pi*x)]     % X er designmatrisen

beta=X'*X \ X'*y    %gir minste kvadraters losning

t=linspace(-2, 6); 

f= beta(1)+ beta(2)*sin((0.5)*pi*t) + beta(3)*cos((0.5)*pi*t);  %


figure(12)

plot(t,f,'b') 

hold on;

plot(x, y,'or')


