DINI22

clc

format short

disp (' **********************************************************')
disp (' *** Ajustamento de uma linha de nivelamento geometrico ***')
disp (' **********************************************************')
disp (' **********************************************************')
disp (' ******************percurso 1 DINI22***********************')
disp (' **********************************************************')

pontos=cellstr(['T1';'P7';'C';'P';'B';'A';'I1';'I2';'I3';'I4';'I5';'I6';'I7';'I8';'I9';'I10';'M3']);

%cota_T1+desnivel_T1P7+residuo_T1P7-cota_P7=0
%cota_P7+desnivel_P7C+residuo_P7C-cota_C=0
%cota_C+desnivel_CP+residuo_CP-cota_P=0
%cota_P+desnivel_PB+residuo_PB-cota_B=0
%cota_B+desnivel_BA+residuo_BA-cota_A=0
%cota_A+desnivel_AI1+residuo_AI1-cota_I1=0
%                   ...
%cota_I10+desnivel_I10M3+residuo_I10M3-cota_M3=0
%
%
%Os pontos começados com I são pontos intermédios

A=zeros(16,15);

for i=1:15
    A(i,i) = -1;
    A(i + 1,i) = 1;
end

W(1,1)= 80.108-1.8486;
W(2,1)=-0.0569;
W(3,1)=-0.2131;
W(4,1)=0.5814;
W(5,1)=-0.1922;
W(6,1)=-0.3585;
W(7,1)=-0.5723;
W(8,1)=-0.2508;
W(9,1)=-0.278;
W(10,1)=-0.5356;
W(11,1)=-0.3216;
W(12,1)=-0.1517;
W(13,1)=-0.3064;
W(14,1)=-0.2759;
W(15,1)=-0.2911;
W(16,1)=0.2109-75.23;
%
%
%
dist(1,1)=53.85;
dist(2,1)=49.69;
dist(3,1)=30.23;
dist(4,1)=37.62;
dist(5,1)=57.45;
dist(6,1)=70.28;
dist(7,1)=54.59;
dist(8,1)=56.10;
dist(9,1)=49.03;
dist(10,1)=60.47;
dist(11,1)=65.17;
dist(12,1)=57.42;
dist(13,1)=59.88;
dist(14,1)=60.22;
dist(15,1)=56.47;
dist(16,1)=9.73;
%
%Dini22 "Standard deviation on 1 km of double levelling" + "foldable bar code staff"
%sigma_Dini22 = 1.3mm
%
%Dini22 "Accuracy of distance measuremen" + "foldable bar code staff"
%sigma_Dini22 = 30mm
%
sigma_dini22=30/1000;

sigma=zeros(16,16);

for i=1:16
    sigma(i,i)=sigma_dini22*sqrt(2*dist(i,1)/1000);
end

sigma2=sigma^2;

sigma02=1;

QL=sigma2/sigma02;

P=inv(QL);
%
%*************AJUSTAMENTO**************************
%
At= A.';

N=At*P*A;

U=At*P*(-W);
%
% parametros ajustados
%
X=inv(N)*U
cotas=X;
%
% matriz cofactor dos parametros
%
Qx=inv(N);
%
% PRECISAO DAS COTAS
%
for i=1:15
    precisao_cotas(i,1)=sqrt(Qx(i,i));
end
%
% RESIDUOS
%
v=A*X+W;
%
%
%
Qv=(inv(P)-A*inv(N)*At);
%
% OBSERVACOES AJUSTADAS
%
La(1,1)=-1.8486+v(1,1);
La(2,1)=-0.0569+v(2,1);
La(3,1)=-0.2131+v(3,1);
La(4,1)=0.5814+v(4,1);
La(5,1)=-0.1922+v(5,1);
La(6,1)=-0.3585+v(6,1);
La(7,1)=-0.5723+v(7,1);
La(8,1)=-0.2508+v(8,1);
La(9,1)=-0.278+v(9,1);
La(10,1)=-0.5356+v(10,1);
La(11,1)=-0.3216+v(11,1);
La(12,1)=-0.1517+v(12,1);
La(13,1)=-0.3064+v(13,1);
La(14,1)=-0.2759+v(14,1);
La(15,1)=-0.2911+v(15,1);
La(16,1)=0.2109+v(16,1);

La
%
%
%
QLa=inv(P)-Qv;
%
% variancia a posteriori por unidade de peso
%
s02=v.'*P*v/(16-15);

% distribuicao do qui-quadrado a 0.05%

qui(1,1)=3.841;
qui(2,1)=5.991;
qui(3,1)=7.815;
qui(4,1)=9.488;
qui(5,1)=11.071;
qui(6,1)=12.592;
qui(7,1)=14.067;
qui(8,1)=15.507;
qui(9,1)=16.919;
qui(10,1)=18.307;
qui(11,1)=19.675;
qui(12,1)=21.026;
qui(13,1)=22.362;
qui(14,1)=23.685;
qui(15,1)=24.996;
qui(16,1)=26.296;
qui(17,1)=27.587;
qui(18,1)=28.869;
qui(19,1)=30.144;
qui(20,1)=31.410;
qui(21,1)=32.671;
qui(22,1)=33.924;
qui(23,1)=35.172; 
qui(24,1)=36.415;
qui(25,1)=37.652; 
qui(26,1)=38.885; 
qui(27,1)=40.113; 
qui(28,1)=41.337; 
qui(29,1)=42.557; 
qui(30,1)=43.773; 
qui(31,1)=44.985;
qui(32,1)=46.194; 
qui(33,1)=47.400; 
qui(34,1)=48.602; 
qui(35,1)=49.802;
qui(36,1)=50.998; 
qui(37,1)=52.192; 
qui(38,1)=53.384; 
qui(39,1)=54.572; 
qui(40,1)=55.758;
qui(41,1)=56.942; 
qui(42,1)=58.124; 
qui(43,1)=59.304; 
qui(44,1)=60.481; 
qui(45,1)=61.656; 

% teste da razao de variancias

y=s02*(4-3)/sigma02;

if (y < qui(4-3,1)) 
% sigma02 e s02 estatisticamente semelhantes    
    'Aceitar o ajustamento.'
else
    'Rejeitar o ajustamento.'
end