===== Cvičení 6 =====
[[https://cw.felk.cvut.cz/lib/exe/fetch.php/courses/a4b33opt/cviceni/ls.pdf?id=courses%3Aa4b33opt%3Acviceni%3Astart&cache=cache|Zadání]]

<code matlab>
% 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;

%
</code>


===== druhý příklad =====

<code matlab>
clear
close all

%nacteme
Data = load('teplota.txt','-ascii');
days = Data(:,1);
temp = Data(:,2);
COUNT = size(days, 1);
PI = 3.14159265;

%matice podle PDFka
A = [ones(COUNT, 1) days sin(2*PI*days/365) cos(2*PI*days/365) ];
%brutalne slozity vypocet
x = A\temp

hold on
plot(days, temp, 'ro');
viz = 1:1:1200;  %pouzijeme misto dnu, pro jemnejsi kresbu grafu
plot(viz, x(1,1) + x(2,1)*viz + x(3,1)*sin(2*PI*viz/365) + x(4,1)*cos(2*PI*viz/365), 'b-'); %vykresleni interpolovane funkce
hold off


T = A*x - temp;  
T'*T  %kvadrat kriterialni funkce

</code>

~~DISCUSSION~~