Přetransformováno

     Dneska od devíti přišla na řadu Geodézie v pořadí Čtvrtá. Byla jsem docela dost nervní. Hrnec s gulášem usazený v mé hlavě mi včera téměř začal přetékat. Něco jsem se snažila nadupat do hlavy, ale ve čtyři jsem už jenom koukala doblba. Nechala jsem snahy přivést se do hrobu a šla si raději Zatrůnit... Game of thrones s02e09 ... už jenom jeden díl!!! áááááá nééé

     Po půl devátý jsem dorazila do devátého patra. Nikdo tam nebyl. Jenom jeden človíček ze třeťáku, který nešel na geodézii 4 ale na katastr. Za chvíli se přišoural Skořepa a otevřel mi učebnu B-977 - s počítači. Zalezla jsem si k jenomu. Po pěti nebo deseti minutách dorazili další dva spolutrpitelé. Podle školního ptáka - KOSu měl být tento termín plný - předpokládám 8 míst jako bývalo. Domnívali jsme se, že bude zkouška na dlouho, zvlášť ústní, když každého zkoušel alespoň hodinu. Byli jsme smířeni s tím, že se domu  půjde až tak v šest, nejdřív. Najednou bylo devět a nás pořád trojice. Nikdo už nedorazil.

     Písemka byla na Helmertovu transformaci - vypočítat opravy vx a vy, aposteriorní jednotkovou chybu a míru identity, poklesy sum čtverců - podle toho vybrat bod, který by se měl vyloučit z transformace a nakonec nějakou divnou chybu, která se spočítala jako odmocnina z dvojnásobku aposteriorní mínus nějaká divná chyba prvního bodu z rajonu... neslyšela jsem teda o ní nikdy dřív, ale vzoreček byl uveden na papíře se zadáním...

     Měli jsme na to čas hoďku a čtvrt - do půl jedenáctý. Zanedlouho jsme dostali takovej další testík, dost podivnej, trochu jsem ho nepochopila chvílema... ale něco vyřešila a dokonce i dobře XD ... Hned na to to s náma probral, co jsme ten test odevzdali a posléze plynule přešel na ústní. Ptal se na Helmertovku obecně, u ní si zas vykecával játra trochu, až došel k obyčejný podobnostní transformaci a shodnostní a zakončil to tim, že se mě zeptal na afinní. Pak to už neprotahoval a dal nám všem B... asi Blbá nálada XD

skriptík co jsem vyplodila: (ta dvojka - lomeno dvěma tam prej nemá být)

clear
clc
format long g
n=6;

M=[
% psi    s
34.4325 109.750 %1
38.0851 116.838 %2
40.1090 113.723 %3
39.2505 111.869 %4
42.5186 107.407 %5
39.8077 101.748 %6
];

s=[
%y   x   
 12  9.6 %1   
 12  0   %2
 7.2 0   %3
 7.2 2.4 %4
 0   2.4 %5
 0   9.6 %6
];

%S=zeros(n,2);
%for i=1:n
%   S(i,1)=M(i,2)*sin(M(i,1)*pi/200);
%   S(i,2)=M(i,2)*cos(M(i,1)*pi/200);
%end
%S;

S =[
    %Y        X
   56.5078   94.0847
   65.8018   96.5466
   67.0021   91.8893
   64.6849   91.2718
   66.5197   84.3291
   59.5571   82.4961
];

Yt=sum(S(:,1))/n;
Xt=sum(S(:,2))/n;
yt=sum(s(:,1))/n;
xt=sum(s(:,2))/n;

Yr=S(:,1)-Yt;
Xr=S(:,2)-Xt;
yr=s(:,1)-yt;
xr=s(:,2)-xt;


l1=sum(xr.*Xr+yr.*Yr)/sum(xr.^2+yr.^2);
l2=sum(xr.*Yr-yr.*Xr)/sum(xr.^2+yr.^2);

tx=Xt-(l1*xt-l2*yt);
ty=Yt-(l1*yt+l2*xt);

Xpo=l1.*s(:,2)-l2.*s(:,1)+tx;
Ypo=l1.*s(:,1)+l2.*s(:,2)+ty;

vx=S(:,2)-Xpo
vy=S(:,1)-Ypo

sv=sqrt(sum(vx.^2+vy.^2)/n);

so=sqrt(sum(vx.^2+vy.^2)/(2*n-4));

delta=zeros(6,1);
for i=1:6
  delta(i,1)=(vx(i,1)^2+vy(i,1)^2)/( (n-1)/n - (xr(i,1)^2+yr(i,1)^2)/(sum(xr.^2+yr.^2)) ); 
end
delta

%delta=(vx.^2+vy.^2)/( (n-1)/n - (xr.^2+yr.^2)/(sum(xr.^2+yr.^2)) )

%rajon
sd=0.002;
spsi=0.0010*(pi/200);

sigX=sqrt( (cos(M(1,1)*(pi/200)))^2 * sd^2 + (M(1,2)*sin(M(1,1)*(pi/200)))^2 * spsi^2  )
sigY=sqrt( (sin(M(1,1)*(pi/200)))^2 * sd^2 + (M(1,2)*cos(M(1,1)*(pi/200)))^2 * spsi^2  )

sxy=sqrt((sigX^2+sigY^2)/2)

roz=sqrt(2*so^2-sxy^2)