clear
clc
format long G
%výpoèet smìrodatné odchylky laserového dálkomìru
A = load ('3dalk_20cm.txt');
B = load ('3dalk_30cm.txt');
C = load ('3dalk_40cm.txt');
D = load ('3dalk_50cm.txt');
E = load ('3dalk_60cm.txt');
F = load ('3dalk_70cm.txt');
G = load ('3dalk_80cm.txt');
H = load ('3dalk_90cm.txt');
I = load ('3dalk_100cm.txt');
J = load ('3dalk_110cm.txt');
K = load ('3dalk_120cm.txt');

%% vzdálenost 20cm
n = length(A);   
%prùmìr
prum1 = sum(A(:,1))/n;
prum2 = sum(A(:,2))/n;
prum3 = sum(A(:,3))/n;
%výpoèet oprav
for i = 1:n
      o1(i) = prum1 - A(i,1);
      o2(i) = prum2 - A(i,2);
      o3(i) = prum3 - A(i,3);
end
%smìrodatná odchylka
sm_od1 = sqrt(sum(o1.^2)/(n-1));
sm_od2 = sqrt(sum(o2.^2)/(n-1));
sm_od3 = sqrt(sum(o3.^2)/(n-1));
%graf namìøených hodnot
figure(1)
plot(A(:,1),'r')
title('Test laserových dálkomìrù (20cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(A) 190 215])
hold on
plot(A(:,2),'b')
plot(A(:,3),'g')
hold off

%% vzdálenost 30cm
n1 = length(B);   
%prùmìr
prum4 = sum(B(:,1))/n1;
prum5 = sum(B(:,2))/n1;
prum6 = sum(B(:,3))/n1;
%výpoèet oprBv
for i = 1:n1
      o4(i) = prum4 - B(i,1);
      o5(i) = prum5 - B(i,2);
      o6(i) = prum6 - B(i,3);
end
%smìrodatná odchylka
sm_od4 = sqrt(sum(o4.^2)/(n1-1));
sm_od5 = sqrt(sum(o5.^2)/(n1-1));
sm_od6 = sqrt(sum(o6.^2)/(n1-1));
%graf namìøených hodnot
figure(2)
plot(B(:,1),'r')
title('Test laserových dálkomìrù (30cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(B) 190 215])
hold on
plot(B(:,2),'b')
plot(B(:,3),'g')
hold off
%% vzdálenost 40cm
n2 = length(C);   
%prùmìr
prum7 = sum(C(:,1))/n2;
prum8 = sum(C(:,2))/n2;
prum9 = sum(C(:,3))/n2;
%výpoèet oprCv
for i = 1:n2
      o7(i) = prum7 - C(i,1);
      o8(i) = prum8 - C(i,2);
      o9(i) = prum9 - C(i,3);
end
%smìrodatná odchylka
sm_od7 = sqrt(sum(o7.^2)/(n2-1));
sm_od8 = sqrt(sum(o8.^2)/(n2-1));
sm_od9 = sqrt(sum(o9.^2)/(n2-1));
%graf namìøených hodnot
figure(3)
plot(C(:,1),'r')
title('Test laserových dálkomìrù (40cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(C) 190 218])
hold on
plot(C(:,2),'b')
plot(C(:,3),'g')
hold off
%% vzdálenost 50Dm
n3 = length(D);   
%prùmìr
prum10 = sum(D(:,1))/n3;
prum11 = sum(D(:,2))/n3;
prum12 = sum(D(:,3))/n3;
%výpoèet oprDv
for i = 1:n3
      o10(i) = prum10 - D(i,1);
      o11(i) = prum11 - D(i,2);
      o12(i) = prum12 - D(i,3);
end
%smìrodatná odchylka
sm_od10 = sqrt(sum(o10.^2)/(n3-1));
sm_od11 = sqrt(sum(o11.^2)/(n3-1));
sm_od12 = sqrt(sum(o12.^2)/(n3-1));
%graf namìøených hodnot
figure(4)
plot(D(:,1),'r')
title('Test laserových dálkomìrù (50cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(D) 1120 2111])
hold on
plot(D(:,2),'b')
plot(D(:,3),'g')
hold off
%% vzdálenost 60cm
n4 = length(E);   
%prùmìr
prum13 = sum(E(:,1))/n4;
prum14 = sum(E(:,2))/n4;
prum15 = sum(E(:,3))/n4;
%výpoèet oprav
for i = 1:n4
      o13(i) = prum13 - E(i,1);
      o14(i) = prum14 - E(i,2);
      o15(i) = prum15 - E(i,3);
end
%smìrodatná odchylka
sm_od13 = sqrt(sum(o13.^2)/(n4-1));
sm_od14 = sqrt(sum(o14.^2)/(n4-1));
sm_od15 = sqrt(sum(o15.^2)/(n4-1));
%graf namìøených hodnot
figure(5)
plot(E(:,1),'r')
title('Test laserových dálkomìrù (60cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(E) 1420 2141])
hold on
plot(E(:,2),'b')
plot(E(:,3),'g')
hold off
%% vzdálenost 70cm
n5 = length(F);   
%prùmìr
prum16 = sum(F(:,1))/n5;
prum17 = sum(F(:,2))/n5;
prum18 = sum(F(:,3))/n5;
%výpoèet oprav
for i = 1:n5
      o16(i) = prum16 - F(i,1);
      o17(i) = prum17 - F(i,2);
      o18(i) = prum18 - F(i,3);
end
%smìrodatná odchylka
sm_od16 = sqrt(sum(o16.^2)/(n5-1));
sm_od17 = sqrt(sum(o17.^2)/(n5-1));
sm_od18 = sqrt(sum(o18.^2)/(n5-1));
%graf namìøených hodnot
figure(6)
plot(F(:,1),'r')
title('Test laserových dálkomìrù (70cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(F) 1720 2171])
hold on
plot(F(:,2),'b')
plot(F(:,3),'g')
hold off
%% vzdálenost 80cm
n6 = length(G);   
%prùmìr
prum19 = sum(G(:,1))/n6;
prum20 = sum(G(:,2))/n6;
prum21 = sum(G(:,3))/n6;
%výpoèet oprav
for i = 1:n6
      o19(i) = prum19 - G(i,1);
      o20(i) = prum20 - G(i,2);
      o21(i) = prum21 - G(i,3);
end
%smìrodatná odchylka
sm_od19 = sqrt(sum(o19.^2)/(n6-1));
sm_od20 = sqrt(sum(o20.^2)/(n6-1));
sm_od21 = sqrt(sum(o21.^2)/(n6-1));
%graf namìøených hodnot
figure(7)
plot(G(:,1),'r')
title('Test laserových dálkomìrù (80cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(G) 2020 2201])
hold on
plot(G(:,2),'b')
plot(G(:,3),'g')
hold off
%% vzdálenost 90cm
n7 = length(H);   
%prùmìr
prum22 = sum(H(:,1))/n7;
prum23 = sum(H(:,2))/n7;
prum24 = sum(H(:,3))/n7;
%výpoèet oprav
for i = 1:n7
      o22(i) = prum22 - H(i,1);
      o23(i) = prum23 - H(i,2);
      o24(i) = prum24 - H(i,3);
end
%smìrodatná odchylka
sm_od22 = sqrt(sum(o22.^2)/(n7-1));
sm_od23 = sqrt(sum(o23.^2)/(n7-1));
sm_od24 = sqrt(sum(o24.^2)/(n7-1));
%graf namìøených hodnot
figure(8)
plot(H(:,1),'r')
title('Test laserových dálkomìrù (90cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(H) 2323 2231])
hold on
plot(H(:,2),'b')
plot(H(:,3),'g')
hold off
%% vzdálenost 100cm
n8 = length(I);   
%prùmìr
prum25 = sum(I(:,1))/n8;
prum26 = sum(I(:,2))/n8;
prum27 = sum(I(:,3))/n8;
%výpoèet oprav
for i = 1:n8
      o25(i) = prum25 - I(i,1);
      o26(i) = prum26 - I(i,2);
      o27(i) = prum27 - I(i,3);
end
%smìrodatná odchylka
sm_od25 = sqrt(sum(o25.^2)/(n8-1));
sm_od26 = sqrt(sum(o26.^2)/(n8-1));
sm_od27 = sqrt(sum(o27.^2)/(n8-1));
%graf namìøených hodnot
figure(9)
plot(I(:,1),'r')
title('Test laserových dálkomìrù (100cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(I) 2626 2531])
hold on
plot(I(:,2),'b')
plot(I(:,3),'g')
hold off
%% vzdálenost 110cm
n9 = length(J);   
%prùmìr
prum28 = sum(J(:,1))/n9;
prum29 = sum(J(:,2))/n9;
prum30 = sum(J(:,3))/n9;
%výpoèet oprav
for i = 1:n9
      o28(i) = prum28 - J(i,1);
      o29(i) = prum29 - J(i,2);
      o30(i) = prum30 - J(i,3);
end
%smìrodatná odchylka
sm_od28 = sqrt(sum(o28.^2)/(n9-1));
sm_od29 = sqrt(sum(o29.^2)/(n9-1));
sm_od30 = sqrt(sum(o30.^2)/(n9-1));
%graf namìøených hodnot
figure(10)
plot(J(:,1),'r')
title('Test laserových dálkomìrù (110cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axJs([0 length(J) 2929 2831])
hold on
plot(J(:,2),'b')
plot(J(:,3),'g')
hold off
%% vzdálenost 120cm
n10 = length(K);   
%prùmìr
prum31 = sum(K(:,1))/n10;
prum32 = sum(K(:,2))/n10;
prum33 = sum(K(:,3))/n10;
%výpoèet oprav
for i = 1:n10
      o31(i) = prum31 - K(i,1);
      o32(i) = prum32 - K(i,2);
      o33(i) = prum33 - K(i,3);
end
%smìrodatná odchylka
sm_od31 = sqrt(sum(o31.^2)/(n10-1));
sm_od32 = sqrt(sum(o32.^2)/(n10-1));
sm_od33 = sqrt(sum(o33.^2)/(n10-1));
%graf namìøených hodnot
figure(11)
plot(K(:,1),'r')
title('Test laserových dálkomìrù (120cm)')
xlabel('Poèet mìøení')
ylabel('Mìøená hodnota [mm]')
% axis([0 length(K) 3232 3131])
hold on
plot(K(:,2),'b')
plot(K(:,3),'g')
hold off
%prùmìrné mìøené vzdálenosti
prum = [prum1;prum2;prum3;prum4;prum5;prum6;prum7;prum8;prum9;prum10;prum11;prum12;prum13;prum14;prum15;prum16;...
      prum17;prum18;prum19;prum20;prum21;prum22;prum23;prum24;prum25;prum26;prum27;prum28;prum29;prum30;prum31;prum32;prum33];
%smìrodatné odchylky
sm_od = [sm_od1;sm_od2;sm_od3;sm_od4;sm_od5;sm_od6;sm_od7;sm_od8;sm_od9;sm_od10;sm_od11;sm_od12;sm_od13;sm_od14;sm_od15;sm_od16;...
      sm_od17;sm_od18;sm_od19;sm_od20;sm_od21;sm_od22;sm_od23;sm_od24;sm_od25;sm_od26;sm_od27;sm_od28;sm_od29;sm_od30;sm_od31;sm_od32;sm_od33];
%výpis výsledkù
vysledky = [prum sm_od]
