% Vypise odchylky namerenych a predpovezenych hodnot uhlu pro
% mikrostrukturu X
addpath('./skripty', './funkce');
clear all;

orientace2
twins

AB=[4,6,2,1,1];
CD=[2,6,2,1,1];
% AB=[1,6,2,1,2];
% CD=[3,6,2,1,2];

UB=U{AB(2)};
UD=U{CD(2)};
QAB=Q{AB(1),AB(2),AB(3),AB(4),AB(5)};
QCD=Q{CD(1),CD(2),CD(3),CD(4),CD(5)};
nAB=normc(inv(QAB*UB)'*n{AB(1),AB(2),AB(3)});
nCD=normc(inv(QCD*UD)'*n{CD(1),CD(2),CD(3)});
mAB=m{AB(1),AB(2),AB(3),AB(4),AB(5)};
mCD=m{CD(1),CD(2),CD(3),CD(4),CD(5)};

t=0;

disp '-----------------------'
disp '-----------------------'
disp '(mAB, nCD)|A'; c=chyba(AC, mAB, nCD, 2); disp(c); t=t+abs(c);
disp '(mAB, nCD)|C'; c=chyba(AC, mAB, nCD, 2); disp(c); t=t+abs(c);
disp '(mAB, mCD)|A'; c=chyba(AC, mAB, mCD, 44); disp(c); t=t+abs(c);
disp '(mAB, mCD)|C'; c=chyba(AC, mAB, mCD, 43.4); disp(c); t=t+abs(c);
disp '(nAB, nCD)|A'; c=chyba(AC, nAB, nCD, 49); disp(c); t=t+abs(c);
disp '(nAB, nCD)|C'; c=chyba(AC, nAB, nCD, 49); disp(c); t=t+abs(c);
disp '(nAB, mCD)|A'; c=chyba(AC, nAB, mCD, 3); disp(c); t=t+abs(c);
disp '(nAB, mCD)|C'; c=chyba(AC, nAB, mCD, 4); disp(c); t=t+abs(c);
disp '(mAB, nAB)|A'; c=chyba(AC, nAB, mAB, 47); disp(c); t=t+abs(c);
disp '(mAB, nAB)|C'; c=chyba(AC, nAB, mAB, 46.9); disp(c); t=t+abs(c);
disp '(mAB, nAB)|D'; c=chyba(BD, nAB, mAB, 1.3); disp(c); t=t+abs(c);
disp '(mCD, nCD)|A'; c=chyba(AC, nCD, mCD, 46); disp(c); t=t+abs(c);
disp '(mCD, nCD)|B'; c=chyba(BD, nCD, mCD, 0); disp(c); t=t+abs(c);
disp '(mCD, nCD)|C'; c=chyba(AC, nCD, mCD, 45); disp(c); t=t+abs(c);
disp '(mAB, BD)|A'; c=chyba(AC, mAB, BD, 67); disp(c); t=t+abs(c);
disp '(mAB, BD)|C'; c=chyba(AC, mAB, BD, 67); disp(c); t=t+abs(c);
disp '(mAB, AC)|D'; c=chyba(BD, mAB, AC, 48.2); disp(c); t=t+abs(c);
disp '(nAB, BD)|A'; c=chyba(AC, nAB, BD, 66.2); disp(c); t=t+abs(c);
disp '(nAB, BD)|C'; c=chyba(AC, nAB, BD, 65.5); disp(c); t=t+abs(c);
disp '(nAB, AC)|D'; c=chyba(BD, nAB, AC, 47.2); disp(c); t=t+abs(c);
disp '(mCD, BD)|A'; c=chyba(AC, mCD, BD, 69.4); disp(c); t=t+abs(c);
disp '(mCD, AC)|B'; c=chyba(BD, mCD, AC, 46.5); disp(c); t=t+abs(c);
disp '(mCD, BD)|C'; c=chyba(AC, mCD, BD, 68.9); disp(c); t=t+abs(c);
disp '(nCD, BD)|A'; c=chyba(AC, nCD, BD, 65); disp(c); t=t+abs(c);
disp '(nCD, AC)|B'; c=chyba(BD, nCD, AC, 46.6); disp(c); t=t+abs(c);
disp '(nCD, BD)|C'; c=chyba(AC, nCD, BD, 65); disp(c); t=t+abs(c);
disp '-----------------------'
disp 'Prumerna odchylka: '; disp(t/26);
disp '-----------------------'

clear all;