% vypocet zapisniku
clear; format compact; format long g;

% zadani ------------------------------------------------------------------

sm = load('meas_dir11.txt'); % seznam mereni
poloh = 2;
skupin = 3;

% vypocet -----------------------------------------------------------------
n = size(sm,1) / poloh / skupin; % pocet smeru
for sk = 1:skupin
    bR = 2*(sk-1)*n+1;   % zacatek I rady v sm
    eR = 2*(sk-1)*n+n;   % konec I rady v sm
% prumer z poloh
    osnova(:,sk) = (  sm( bR:eR, 3 ) + sm( eR+n:-1:bR+n, 3 )-200  )./2;
% redukce
    os_red(:,sk) = osnova(:,sk)-osnova(1,sk);
end

% prumer ze skupin
mean = mean (os_red, 2);

% opravy prumeru
err = repmat(mean,1,skupin)-os_red;

% smerodatna odchylka radkoveho prumeru
smo = std(osnova')/sqrt(skupin);

% output table ------------------------------------------------------------
out_matrix( :, [1:2,3:2:2*skupin+1,4:2:2*skupin+2,2*skupin+3,10] ) = [ sm(1:n,1:2), os_red, err, mean, smo' ]

fprintf( ' st.| cil|     1.sk |   chyba |     2.sk |   chyba |     3.sk |   chyba |   prumer |    std\n')
fprintf( '% 3.0f |% 3.0f |% 9.4f |% 8.4f |% 9.4f |% 8.4f |% 9.4f |% 8.4f |% 9.4f |% 7.4f\n', out_matrix' )
