Zadání

% uloha 1:
Data = load('mzdy.txt','-ascii');
Rok = Data(:,1);
Mzdy = Data(:,2);


hold on;
plot(Rok, Mzdy, 'ob');
A = [ ones(length(Rok), 1) Rok ];

% v uloze dve je predchozi radek v tomto tvaru:
% A = [ ones(length(Den), 1) Den sin((2*pi*Den)/365) cos((2*pi*Den)/365) ];
% jinak je postup obdobny..

b = Mzdy;
x = A \ b;
plot([min(Rok) max(Rok)], [x(1) + x(2)*min(Rok) x(1) + x(2) * max(Rok)], 'g');

%odhad druheho kvartalu 2009
odhad = x(1) + x(2) * 2009.25;
plot(2009.25, odhad, 'or');

%kriterialni funkce F()
F = x'*(A')*A*x - 2*b'*A*x + b'*b;

%