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

orientace1
twins

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

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)};
FCD=F{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)};
mX=normc(inv(FCD)'*rnull(FCD-QAB*UB));

t=0;

disp '-----------------------'
disp '-----------------------'
disp '(mAB, nCD)|A'; c=chyba(AC, nCD, mAB, 9); disp(c); t=t+abs(c);
disp '(mAB, nCD)|C'; c=chyba(AC, mAB, nCD, 11); disp(c); t=t+abs(c);
disp '(mAB, mX)|A'; c=chyba(AC, mX, mAB, 14); disp(c); t=t+abs(c);
disp '(mAB, mX)|C'; c=chyba(AC, mAB, mX, 17.3); disp(c); t=t+abs(c);
disp '(mAB, mCD)|A'; c=chyba(AC, mAB, mCD, 37); disp(c); t=t+abs(c);
disp '(mAB, mCD)|C'; c=chyba(AC, mAB, mCD, 35.7); disp(c); t=t+abs(c);
disp '(nAB, nCD)|A'; c=chyba(AC, nCD, nAB, 40); disp(c); t=t+abs(c);
disp '(nAB, nCD)|C'; c=chyba(AC, nAB, nCD, 42); disp(c); t=t+abs(c);
disp '(nAB, mCD)|A'; c=chyba(AC, mCD, nAB, 7); disp(c); t=t+abs(c);
disp '(nAB, mCD)|C'; c=chyba(AC, nAB, mCD, 4); disp(c); t=t+abs(c);
disp '(nAB, Mx)|A'; c=chyba(AC, mX, nAB, 45); disp(c); t=t+abs(c);
disp '(nAB, mX)|C'; c=chyba(AC, nAB, mX, 50.2); disp(c); t=t+abs(c);
disp '(mAB, nAB)|A'; c=chyba(AC, mAB, nAB, 30.3); disp(c); t=t+abs(c);
disp '(mAB, nAB)|C'; c=chyba(AC, mAB, nAB, 32); disp(c); t=t+abs(c);
disp '(mAB, nAB)|D'; c=chyba(BD, nAB, mAB, 16.2); 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, 71); disp(c); t=t+abs(c);
disp '(mCD, nCD)|C'; c=chyba(AC, mCD, nCD, 46.3); disp(c); t=t+abs(c);
disp '(mCD, mX)|A'; c=chyba(AC, mX, mCD, 51); disp(c); t=t+abs(c);
disp '(mCD, mX)|B'; c=chyba(BD, mCD, mX, 23); disp(c); t=t+abs(c);
disp '(mCD, mX)|C'; c=chyba(AC, mCD, mX, 55.2); disp(c); t=t+abs(c);
disp '(nCD, mX)|A'; c=chyba(AC, mX, nCD, 5.3); disp(c); t=t+abs(c);
disp '(nCD, mX)|B'; c=chyba(BD, nCD, mX, 87.3); disp(c); t=t+abs(c);
disp '(nCD, mX)|C'; c=chyba(AC, nCD, mX, 7); disp(c); t=t+abs(c);
disp '(mAB, BD)|A'; c=chyba(AC, mAB, BD, 80.8); disp(c); t=t+abs(c);
disp '(mAB, BD)|C'; c=chyba(AC, mAB, BD, 80); disp(c); t=t+abs(c);
disp '(mAB, AC)|D'; c=chyba(BD, mAB, AC, 47); disp(c); t=t+abs(c);
disp '(nAB, BD)|A'; c=chyba(AC, nAB, BD, 49); disp(c); t=t+abs(c);
disp '(nAB, BD)|C'; c=chyba(AC, nAB, BD, 48); disp(c); t=t+abs(c);
disp '(nAB, AC)|D'; c=chyba(BD, nAB, AC, 62); disp(c); t=t+abs(c);
disp '(mCD, BD)|A'; c=chyba(AC, mCD, BD, 43);disp(c); t=t+abs(c);
disp '(mCD, AC)|B'; c=chyba(BD, mCD, AC, 65); disp(c); t=t+abs(c);
disp '(mCD, BD)|C'; c=chyba(AC, mCD, BD, 43.7); disp(c); t=t+abs(c);
disp '(nCD, BD)|A'; c=chyba(AC, nCD, BD, 88.3); disp(c); t=t+abs(c);
disp '(nCD, AC)|B'; c=chyba(BD, nCD, AC, 45); disp(c); t=t+abs(c);
disp '(nCD, BD)|C'; c=chyba(AC, nCD, BD, 88); disp(c); t=t+abs(c);
disp '(mX, BD)|A'; c=chyba(AC, mX, BD, 86); disp(c); t=t+abs(c);
disp '(mX, AC)|B'; c=chyba(BD, mX, AC, 42); disp(c); t=t+abs(c);
disp '(mX, BD)|C'; c=chyba(AC, mX, BD, 82); disp(c); t=t+abs(c);
disp '-----------------------'
disp 'Prumerna odchylka: '; disp(t/39);
disp '-----------------------'

clear all;