% Eksempel -- Lineaer regresjon

clear all; format short;

x=(0:0.2:4)'; m=21; % antall punkter

y1=ones(m,1)+2*x;  

y=y1+ 0.5*(0.5*ones(m,1)-rand(m,1)); 
% antar x og y er koordinatene til maalepunkt

figure(3)
plot(x, y,'or')
         % plotter disse punktene

X=[ones(m,1) x];     % X er designmatrisen

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

t=linspace(-1, 5); 

f= beta(1)+ beta(2)*t;  %gir minste kvadraters linje

figure(4)

plot(t,f,'b')  % plotter m.k.-linjen 

hold on

plot(x, y,'or') % og maalepunktene



    