% Eksempel -- vektet linear regresjon

clear all;
m=4;        % antall punkter

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

figure(15)

plot(x, y,'or')
hold on;

d = [3; 3; 2; 1]; W=diag(d) 
% W er vektmatrisen


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

WX=W*X
Wy=W*y

beta = WX'* WX \ WX'* Wy    
% beta gir vektet minste kvadraters losning


t=linspace(-1, 4); 

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

plot(t,vf,'b')  % plotter vektet m.k.-linje 

hold on;

f= 0.6 + 1.1*t;
plot(t,f,'g') % plotter ogsaa den vanlig m.k.-linjen

