150

1. AJUSTAMENTO DE UM POSICIONAMENTO POR ESTAÇÃO LIVRE

Coordenadas dos pontos de referência

𝑋1 = 842,281𝑚 𝑋2 = 1337,544𝑚 𝑋3 = 1505,177𝑚

𝑌1 = 925,523𝑚 𝑌2 = 996,249𝑚 𝑌3 = 882,957𝑚
 151

Observações realizadas

Distâncias observadas
𝐷01 = 244,512𝑚
𝐷02 = 321,570𝑚
𝐷03 = 443,501𝑚

Ângulos observados
𝐴102 = 123°31′01"
𝐴203 = 24°39′27"

1.1. Vetor das observações (Lb)

𝐷01 244,512
𝐷02 321,570
𝐿𝑏 = 𝐷03 = 443,501 Ângulos em radianos
𝐷𝑀𝑆(𝐴102 ) 2,15578
[𝐷𝑀𝑆(𝐴203 )] [0,43035]
 152

1.2. Cálculo das coordenadas aproximadas da estação (0)

𝐷12 = √(𝑋2 − 𝑋1)2 + (𝑌2 − 𝑌1)2 = 500,288

𝑋2 − 𝑋1
𝐴𝑍12 = 𝑎𝑡𝑎𝑛 | | = 81°52′22"
𝑌2 − 𝑌1
2 2 2
𝐷02 − 𝐷01 − 𝐷12
𝐴210 = 𝑎𝑐𝑜𝑠 ( ) = 32°21′23"
−2. 𝐷12 . 𝐷01
𝐴𝑍10 = 𝐴𝑍12 + 𝐴210 = 114°13′46"
𝑋0 = 𝑋1 + 𝐷01 . sin(𝐴𝑍10 ) = 1065,254
𝑌0 = 𝑌1 + 𝐷01 . cos(𝐴𝑍10 ) = 825,178

1065,254
𝑋0 = [ ]
825,178
 153

1.3. Cálculo das observações em função dos parâmetros (L0)

√(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2

√(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2 244,512
√(𝑋3 − 𝑋0)2 + (𝑌3 − 𝑌0)2 321,570
𝐿0 = = 443,701
𝑋2 − 𝑋0 𝑋1 − 𝑋0
𝑎𝑡𝑎𝑛 | | − (360° − 𝑎𝑡𝑎𝑛 | |) + 360° 2,15776
𝑌2 − 𝑌0 𝑌1 − 𝑌0 [0,43035]
𝑋3 − 𝑋0 𝑋2 − 𝑋0
𝑎𝑡𝑎𝑛 | | − 𝑎𝑡𝑎𝑛 | |
[ 𝑌3 − 𝑌0 𝑌2 − 𝑌0 ]

1.4. Cálculo do vetor diferença entre observações (L)

−0,000
0,000
𝐿 = 𝐿𝑜 − 𝐿𝑏 = 0,200
0,00199
[−0,00000]
 154

1.5. Cálculo da Matriz (A)

𝜕𝐷01
𝐴1,1 = √(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2
𝜕𝑋0
𝜕𝐷01
𝐴1,2 = √(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2
𝜕𝑌0
𝜕𝐷02
𝐴2,1 = √(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2
𝜕𝑋0
𝜕𝐷02
𝐴2,2 = √(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2
𝜕𝑌0
𝜕𝐷03
𝐴3,1 = √(𝑋3 − 𝑋0)2 + (𝑌3 − 𝑌0)2
𝜕𝑋0
𝜕𝐷03
𝐴3,2 = √(𝑋3 − 𝑋0)2 + (𝑌3 − 𝑌0)2
𝜕𝑌0
𝜕𝐴102 𝑋2 − 𝑋0 𝑋1 − 𝑋0
𝐴4,1 = 𝑎𝑡𝑎𝑛 | | − 𝑎𝑡𝑎𝑛 | |
𝜕𝑋0 𝑌2 − 𝑌0 𝑌1 − 𝑌0
𝜕𝐴102 𝑋2 − 𝑋0 𝑋1 − 𝑋0
𝐴4,2 = 𝑎𝑡𝑎𝑛 | | − 𝑎𝑡𝑎𝑛 | |
𝜕𝑌0 𝑌2 − 𝑌0 𝑌1 − 𝑌0
𝜕𝐴203 𝑋3 − 𝑋0 𝑋2 − 𝑋0
𝐴5,1 = 𝑎𝑡𝑎𝑛 | | − 𝑎𝑡𝑎𝑛 | |
𝜕𝑋0 𝑌3 − 𝑌0 𝑌2 − 𝑌0
𝜕𝐴203 𝑋3 − 𝑋0 𝑋2 − 𝑋0
𝐴5,2 = 𝑎𝑡𝑎𝑛 | | − 𝑎𝑡𝑎𝑛 | |
𝜕𝑌0 𝑌3 − 𝑌0 𝑌2 − 𝑌0
 155

0,9119 −0,4104
−0,8468 −0,5320
𝐴 = −0,9915 −0,1302
0,0000 0,0064
[ 0,0014 −0,0004]

1.6. Atribuição da matriz Peso (P)

Neste caso como não foi dado as incertezas das observações vamos considerar
peso unitário
𝑃𝑖,𝑖 = 𝜎02 . Σ𝐿𝑏 −1 = 1

1 0 0 0 0
0 1 0 0 0
𝑃= 0 0 1 0 0
0 0 0 1 0
[0 0 0 0 1]

1.7. Cálculo da correção nas coordenadas aproximadas (CX)

𝐶𝑋 = −(𝐴𝑇 . 𝐴)−1 . 𝐴𝑇 . 𝐿
0,077
𝐶𝑋 = | |
0,022
 156

1.8. Coordenadas corrigidas ou ajustadas (Xa)

𝑋𝑎 = 𝑋0 + 𝐶𝑋
1065,331
𝑋𝑎 = | |
825,200

Tendo em vista que houve correção devemos continuar aproximando as

coordenadas aproximadas das coordenadas ajustadas até não haver mais
correção, ou seja, 𝐶𝑋 ≅ 0

1ª ITERAÇÃO

1.9. Coordenada aproximada da estação do primeiro ajuste

𝑋𝑎1
𝑋0 = | |
𝑋𝑎2
 157

1.10. Cálculo do vetor L0

√(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2
√(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2 244,573
√(𝑋3 − 𝑋0)2 + (𝑌3 − 𝑌0)2 321,493
𝐿0 = = 443,622
𝑋2 − 𝑋0 𝑋1 − 𝑋0
𝑎𝑡𝑎𝑛 | | − (2𝜋 − 𝑎𝑡𝑎𝑛 | |) + 2𝜋 2,15791
𝑌2 − 𝑌0 𝑌1 − 𝑌0 [0,43045]
𝑋3 − 𝑋0 𝑋2 − 𝑋0
𝑎𝑡𝑎𝑛 | | − 𝑎𝑡𝑎𝑛 | |
[ 𝑌3 − 𝑌0 𝑌2 − 𝑌0 ]

1.11. Cálculo do vetor L

0,061
−0,077
𝐿 = 𝐿𝑜 − 𝐿𝑏 = 0,121
0,00213
[0,00010]

1.12. Cálculo da correção CX

Cálculo da correção nas coordenadas aproximadas (CX)
𝐶𝑋 = −(𝐴𝑇 . 𝐴)−1 . 𝐴𝑇 . 𝐿
 158

−0,000
𝐶𝑋 = | |
0,000
1.13. Coordenadas corrigidas ou ajustadas (Xa)
𝑋𝑎 = 𝑋0 + 𝐶𝑋
1065,331
𝑋𝑎 = | |
825,200

1.14. Cálculo do vetor dos resíduos

0,061 0,061
−0,077 −0,077
𝑉 = 𝐴. 𝐶𝑋 + 𝐿 = 0,121 = 0,121
0,00213 0°07′19,1"
[0,00010] [0°00′19,6"]

𝑉̅ = 0,022 𝜎𝑉 = 0,074

1.15. Cálculo da variância da unidade de peso a posteriori

𝑇
𝑉 . 𝑃. 𝑉
𝜎̂𝑜2 = = 0,0080991
𝑛−𝑢
Obs: a variância a posteriori deve ser próximo de 1
 159

1.16. Avaliação das variâncias em função do teste do qui-quadrado

→ muito baixo! 𝑉está adequado, então o problema está na matriz peso (está
subestimada), ou seja, está menor do que deveria ser.
𝑉 𝑇 . 𝑃. 𝑉
𝑐𝑎𝑙𝑐 =
2
= 0,0243
𝜎02
23;0,025 = 0,216 23;0,975 = 9,348

0,216 ≤ 0,0243 ≤ 9,348

Porém 23;0,025 > 2𝑐𝑎𝑙𝑐 ∴ 𝐻0 é rejeitada

Sendo assim deve-se mudar a matriz peso.

1.17. Matriz Peso

Para passar no teste do qui-quadrado basta usar o 𝜎𝑜𝑏𝑠 = 𝑉

 160

𝑉1−2 0 0 0 0
0 𝑉2−2 0 0 0
𝑃= 0 0 𝑉3−2 0 0
0 0 0 (𝑉4 . 𝐷02)−2 0
[ 0 0 0 0 (𝑉5 . 𝐷03)−2 ]

270,35 0 0 0 0
0 170,44 0 0 0
𝑃= 0 0 67,90 0 0
0 0 0 2,13 0
[ 0 0 0 0 560,92]

1.18. Cálculo da correção CX

Cálculo da correção nas coordenadas aproximadas (CX)
𝐶𝑋 = −(𝐴𝑇 . 𝑃. 𝐴)−1 . 𝐴𝑇 . 𝑃. 𝐿
−0,043
𝐶𝑋 = | |
0,002

1.19. Coordenadas corrigidas ou ajustadas (Xa)

𝑋𝑎 = 𝑋0 + 𝐶𝑋
 161

1065,287
𝑋𝑎 = | |
825,202

1ª ITERAÇÃO

1.20. Coordenada aproximada da estação do primeiro ajuste

𝑋𝑎1
𝑋0 = | |
𝑋𝑎2

1.21. Cálculo do vetor L0

√(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2
√(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2 244,533
√(𝑋3 − 𝑋0)2 + (𝑌3 − 𝑌0)2 321,529
𝐿0 = = 443,665
𝑋2 − 𝑋0 𝑋1 − 𝑋0
𝑎𝑡𝑎𝑛 | | − (2𝜋 − 𝑎𝑡𝑎𝑛 | |) + 2𝜋 2,15792
𝑌2 − 𝑌0 𝑌1 − 𝑌0 [0,43039]
𝑋3 − 𝑋0 𝑋2 − 𝑋0
𝑎𝑡𝑎𝑛 | | − 𝑎𝑡𝑎𝑛 | |
[ 𝑌3 − 𝑌0 𝑌2 − 𝑌0 ]
 162

1.22. Cálculo do vetor L

0,021
−0,041
𝐿 = 𝐿𝑜 − 𝐿𝑏 = 0,164
0,00214
[0,00004]

1.23. Cálculo da correção CX

Cálculo da correção nas coordenadas aproximadas (CX)
𝐶𝑋 = −(𝐴𝑇 . 𝐴)−1 . 𝐴𝑇 . 𝐿
−0,000
𝐶𝑋 = | |
0,000

1.24. Coordenadas corrigidas ou ajustadas (Xa)

𝑋𝑎 = 𝑋0 + 𝐶𝑋
1065,287
𝑋𝑎 = | |
825,202
 163

1.25. Cálculo do vetor dos resíduos

0,021 0,021
−0,041 −0,041
𝑉 = 𝐴. 𝐶𝑋 + 𝐿 = 0,164 = 0,164
0,00214 0°07′21,7"
[0,00004] [0°00′07,3"]

𝑉̅ = 0,029 𝜎𝑉 = 0,079

Cálculo da variância a posteriori

2
𝑉 𝑇 𝑃𝑉 𝑉 𝑇 𝑃𝑉
𝜎̂0 = = = 0,742
𝑛−𝑢 3

Agora a 𝜎̂02 ≅ 𝜎02

2𝑐𝑎𝑙𝑐 = 𝑉 𝑇 𝑃𝑉 = 2,227
0,216 ≤ 2,227 ≤ 9,348 -> 𝐻0 não é rejeitada
 164

1.26. MVC de 𝑿𝒂

𝑋𝑎 = 𝜎̂02 (𝐴𝑇 𝑃𝐴)−1

0,0018 0,0003
𝑋𝑎 = [ ]
0,0003 0,0079

Sigma com 95% de nível de confiança

𝜎95 = 1,96

𝜎𝑋0 = 1,96. √𝑋𝑎1,1 = ±0,072

𝜎𝑌0 = 1,96. √𝑋𝑎2,2 = ±0,150

𝑋0 1065,287 𝑚 0,072 𝑚
[ ]=[ ]±[ ]
𝑌0 825,202 𝑚 0,150 𝑚

1.27. Vetor das observações ajustadas (𝑳𝒂)

𝐿𝑎 = 𝐿𝑏 + 𝑉
 165

244,573 244,573
321,493 321,493
𝐿𝑎 = 443,622 = 443,622
2,15791 123°38′20"
[0,43045] [ 24°39′47" ]

1.28. MVC da 𝑳𝒂

𝐿𝑎 = 𝐴. Σ𝑋𝑎. 𝐴𝑇

𝜎𝐿𝑎𝑖 = 1,96. √𝐿𝑎𝑖,𝑖 𝑖 = 1…6

 166

0,134
0,126
𝜎𝐿𝑎 = 0,091
0,00137
[0,00016]

244,573 𝑚 0,134 𝑚
321,493 𝑚 0,126 𝑚
𝐿𝑎 = 443,622 𝑚 ± 0,091 𝑚
123°38′20" 0°04′43"
[ 24°39′47" ] [0°00′34"]
 167

2. AJUSTAMENTO DE UM POSICIONAMENTO POR ESTAÇÃO LIVRE COM

INCERTEZAS DAS MEDIDAS

Coordenadas dos pontos de referência

𝑋1 = 842,281𝑚 𝑋2 = 1337,544𝑚 𝑋3 = 1505,177𝑚

𝑌1 = 925,523𝑚 𝑌2 = 996,249𝑚 𝑌3 = 882,957𝑚
 168

Observações realizadas
O equipamento utilizado foi uma estação total com uma precisão angular de ±9"
e linear de ±(5𝑚𝑚 + 5𝑝𝑝𝑚).

Distâncias observadas
𝐷01 = 244,512𝑚
𝐷02 = 321,570𝑚
𝐷03 = 443,501𝑚

Cálculo das incertezas das distâncias

𝜎𝐷01 = 0,005 + 𝐷01 × 5 × 10−6 = ± 0,006𝑚
𝜎𝐷02 = 0,005 + 𝐷02 × 5 × 10−6 = ± 0,007𝑚
𝜎𝐷03 = 0,005 + 𝐷03 × 5 × 10−6 = ± 0,007𝑚

Ângulos observados
𝐴102 = 123°31′01"
𝐴203 = 24°39′27"
 169

9"

A102

Cálculo das incertezas dos ângulos

𝜎𝐴102 = 𝐷02 × 9. 𝑠𝑒𝑛(1") = 0,014

𝜎𝐴203 = 𝐷03 × 9. 𝑠𝑒𝑛(1") = 0,019
 170

2.1. Vetor das observações (Lb)

𝐷01 244,512
𝐷02 321,570
𝐿𝑏 = 𝐷03 = 443,501 Ângulos em radianos
𝐴102 2,15578
[𝐴203 ] [0,43035]

2.2. Cálculo das coordenadas aproximadas da estação (0)

𝐷12 = √(𝑋2 − 𝑋1)2 + (𝑌2 − 𝑌1)2 = 500,288

𝑋2 − 𝑋1
𝐴𝑍12 = 𝑎𝑡𝑎𝑛 | | = 81°52′22"
𝑌2 − 𝑌1
2 2 2
𝐷02 − 𝐷01 − 𝐷12
𝐴210 = 𝑎𝑐𝑜𝑠 ( ) = 32°21′23"
−2. 𝐷12 . 𝐷01
𝐴𝑍10 = 𝐴𝑍12 + 𝐴210 = 114°13′46"
𝑋0 = 𝑋1 + 𝐷01 . sin(𝐴𝑍10 ) = 1065,254
𝑌0 = 𝑌1 + 𝐷01 . cos(𝐴𝑍10 ) = 825,178
 171

1065,254
𝑋0 = [ ]
825,178

2.3. Cálculo das observações em função dos parâmetros (L0)

√(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2

√(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2 244,512
√(𝑋3 − 𝑋0)2 + (𝑌3 − 𝑌0)2 321,570
𝐿0 = = 443,701
𝑋2 − 𝑋0 𝑋1 − 𝑋0
𝑎𝑡𝑎𝑛 | | − (2𝜋 − 𝑎𝑡𝑎𝑛 | |) + 2𝜋 2,15776
𝑌2 − 𝑌0 𝑌1 − 𝑌0 [0,43035]
𝑋3 − 𝑋0 𝑋2 − 𝑋0
𝑎𝑡𝑎𝑛 | | − 𝑎𝑡𝑎𝑛 | |
[ 𝑌3 − 𝑌0 𝑌2 − 𝑌0 ]

2.4. Cálculo do vetor diferença entre observações (L)

−0,000
0,000
𝐿 = 𝐿𝑜 − 𝐿𝑏 = 0,200
0,00199
[−0,00000]
 172

2.5. Cálculo da Matriz (A)

𝜕𝐷01 𝑋0−𝑋1
𝐴1,1 = √(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2 =
𝜕𝑋0 𝐷01
𝜕𝐷01 𝑌0−𝑌1
𝐴1,2 = √(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2 =
𝜕𝑌0 𝐷01
𝜕𝐷02 𝑋0−𝑋2
𝐴2,1 = √(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2 =
𝜕𝑋0 𝐷02
𝜕𝐷02 𝑌0−𝑌2
𝐴2,2 = √(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2 =
𝜕𝑌0 𝐷02
𝜕𝐷03 𝑋0−𝑋3
𝐴3,1 = √(𝑋3 − 𝑋0)2 + (𝑌3 − 𝑌0)2 =
𝜕𝑋0 𝐷03
𝜕𝐷03 𝑌0−𝑌3
𝐴3,2 = √(𝑋3 − 𝑋0)2 + (𝑌3 − 𝑌0)2 =
𝜕𝑌0 𝐷03
𝜕𝐴102 𝑋2−𝑋0 𝑋1−𝑋0 𝑌2−𝑌0 𝑌1−𝑌0
𝐴4,1 = 𝑎𝑡𝑎𝑛 | 𝑌2−𝑌0 | − 𝑎𝑡𝑎𝑛 | 𝑌1−𝑌0 | = − ( 2 )+( 2 )
𝜕𝑋0 𝐷02 𝐷01
𝜕𝐴102 𝑋2−𝑋0 𝑋1−𝑋0 𝑋2−𝑋0 𝑋1−𝑋0
𝐴4,2 = 𝑎𝑡𝑎𝑛 | 𝑌2−𝑌0 | − 𝑎𝑡𝑎𝑛 | 𝑌1−𝑌0 | = ( 2 )−( 2 )
𝜕𝑌0 𝐷02 𝐷01
𝜕𝐴203 𝑋3−𝑋0 𝑋2−𝑋0 𝑌3−𝑌0 𝑌2−𝑌0
𝐴5,1 = 𝑎𝑡𝑎𝑛 | | − 𝑎𝑡𝑎𝑛 | | = −( 2 )+( 2 )
𝜕𝑋0 𝑌3−𝑌0 𝑌2−𝑌0 𝐷03 𝐷02
𝜕𝐴203 𝑋3−𝑋0 𝑋2−𝑋0 𝑋3−𝑋0 𝑋2−𝑋0
𝐴5,2 = 𝑎𝑡𝑎𝑛 | 𝑌3−𝑌0 | − 𝑎𝑡𝑎𝑛 | 𝑌2−𝑌0 | = ( 2 )−( 2 )
𝜕𝑌0 𝐷03 𝐷02
 173

0,9119 −0,4104
−0,8468 −0,5320
𝐴 = −0,9915 −0,1302
0,0000 0,0064
[ 0,0014 −0,0004]

2.6. Atribuição da matriz Peso (P)

Neste caso como temos as incertezas das observações vamos considera-las na
atribuição da matriz peso
𝑃 = 𝜎02 . Σ𝐿𝑏 −1
𝜎02 = 1
−2
𝜎𝐷01 0 0 0 0
−2
0 𝜎𝐷02 0 0 0
Σ𝐿𝑏 −1 = 0 0 −2
𝜎𝐷03 0 0
−2
0 0 0 𝜎𝐴102 0
[ 0 0 0 0 𝜎𝐴−2
203 ]
 174

2.7. Cálculo da correção nas coordenadas aproximadas (CX)

𝐶𝑋 = −(𝐴𝑇 . 𝑃. 𝐴)−1 . 𝐴𝑇 . 𝑃. 𝐿
0,066
𝐶𝑋 = | |
0,026
2.8. Coordenadas corrigidas ou ajustadas (Xa)
𝑋𝑎 = 𝑋0 + 𝐶𝑋
1065,320
𝑋𝑎 = | |
825,204

Tendo em vista que houve correção devemos continuar aproximando as

coordenadas aproximadas das coordenadas ajustadas até não haver mais
correção, ou seja, 𝐶𝑋 ≅ 0
 175

2.8.1. 1ª Iteração

2.9. Coordenada aproximada da estação do primeiro ajuste

𝑋𝑎1
𝑋0 = | |
𝑋𝑎2

2.10. Cálculo do vetor L0

√(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2

√(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2 244,561
321,500
2
√(𝑋3 − 𝑋0) + (𝑌3 − 𝑌0) 2
𝐿0 = = 443,633
𝑋2−𝑋0 𝑋1−𝑋0 2,15793
𝑎𝑡𝑎𝑛 | 𝑌2−𝑌0 | − (360° − 𝑎𝑡𝑎𝑛 | 𝑌1−𝑌0 |) + 360°
[0,43043]
𝑋3−𝑋0 𝑋2−𝑋0
[ 𝑎𝑡𝑎𝑛 | 𝑌3−𝑌0 | − 𝑎𝑡𝑎𝑛 | 𝑌2−𝑌0 | ]
 176

2.11. Cálculo do vetor L

0,049
−0,070
𝐿 = 𝐿𝑜 − 𝐿𝑏 = 0,132
0,00216
[0,00008]

2.12. Cálculo da correção CX

Cálculo da correção nas coordenadas aproximadas (CX)

𝐶𝑋 = −(𝐴𝑇 . 𝑃. 𝐴)−1 . 𝐴𝑇 . 𝑃. 𝐿
−0,000
𝐶𝑋 = | |
0,000
2.13. Coordenadas corrigidas ou ajustadas (Xa)

𝑋𝑎 = 𝑋0 + 𝐶𝑋
1065,320
𝑋𝑎 = | |
825,204
 177

2.14. Cálculo do vetor dos resíduos

0,049
−0,070
𝑉 = 𝐴. 𝐶𝑋 + 𝐿 = 0,132
0°07′25"
[0°00′16"]

𝑉̅ = 0,023 𝜎𝑉 = 0,074

2.15. Cálculo da variância da unidade de peso a posteriori

2
𝑉 𝑇 . 𝑃. 𝑉
𝜎̂𝑜 = = 168,635
𝑛−𝑢
Obs: a variância a posteriori deve ser próximo de 1

2.16. Avaliação das variâncias em função do teste do qui-quadrado

→ muito alto! 𝑉 está adequado, então o problema está na matriz peso (está
superestimada), ou seja, está maior do que deveria ser.
𝑇
𝑉 . 𝑃. 𝑉
2𝑐𝑎𝑙𝑐 = 2 = 505,905
𝜎0
23;0,025 = 0,216 23;0,975 = 9,348
 178

0,216 ≤ 505,905 ≤ 9,348

Porém 2𝑐𝑎𝑙𝑐 > 23;0,975 ∴ 𝐻0 é rejeitada

2.17. Alteração da matriz peso

Para melhorar a matriz peso vamos considerar 𝜎𝑜𝑏𝑠 = 𝑉

𝑉1−2 0 0 0 0
0 𝑉2−2 0 0 0
𝑃= 0 0 𝑉3−2 0 0
0 0 0 (𝑉4 . 𝐷02)−2 0
[ 0 0 0 0 (𝑉5 . 𝐷02)−2 ]
 179

2.18. Cálculo da correção nas coordenadas aproximadas (CX)

𝐶𝑋 = −(𝐴𝑇 . 𝑃. 𝐴)−1 . 𝐴𝑇 . 𝑃. 𝐿
0,043
𝐶𝑋 = | |
0,005
2.19. Coordenadas corrigidas ou ajustadas (Xa)
𝑋𝑎 = 𝑋0 + 𝐶𝑋
1065,297
𝑋𝑎 = | |
825,183

Tendo em vista que houve correção devemos continuar aproximando as

coordenadas aproximadas das coordenadas ajustadas até não haver mais
correção, ou seja, 𝐶𝑋 ≅ 0
 180

2.20. 1ª Iteração

2.21. Coordenada aproximada da estação do primeiro ajuste

𝑋𝑎1
𝑋0 = | |
𝑋𝑎2

2.22. Cálculo do vetor L0

√(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2

√(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2 244,549
321,531
2
√(𝑋3 − 𝑋0) + (𝑌3 − 𝑌0) 2
𝐿0 = = 443,658
𝑋2−𝑋0 𝑋1−𝑋0 2,15780
𝑎𝑡𝑎𝑛 | 𝑌2−𝑌0 | − (360° − 𝑎𝑡𝑎𝑛 | 𝑌1−𝑌0 |) + 360°
[0,43041]
𝑋3−𝑋0 𝑋2−𝑋0
[ 𝑎𝑡𝑎𝑛 | 𝑌3−𝑌0 | − 𝑎𝑡𝑎𝑛 | 𝑌2−𝑌0 | ]
 181

2.23. Cálculo do vetor L

0,037
−0,039
𝐿 = 𝐿𝑜 − 𝐿𝑏 = 0,157
0,00202
[0,00006]

2.24. Cálculo da correção CX

Cálculo da correção nas coordenadas aproximadas (CX)
𝐶𝑋 = −(𝐴𝑇 . 𝑃. 𝐴)−1 . 𝐴𝑇 . 𝑃. 𝐿
−0,020
𝐶𝑋 = | |
0,017
2.25. Coordenadas corrigidas ou ajustadas (Xa)
𝑋𝑎 = 𝑋0 + 𝐶𝑋
1065,277
𝑋𝑎 = | |
825,199

2.26. 2ª Iteração
 182

2.27. Coordenada aproximada da estação do ajuste anterior

𝑋𝑎1
𝑋0 = | |
𝑋𝑎2

2.28. Cálculo do vetor L0

√(𝑋1 − 𝑋0)2 + (𝑌1 − 𝑌0)2
√(𝑋2 − 𝑋0)2 + (𝑌2 − 𝑌0)2 244,524
321,539
2
√(𝑋3 − 𝑋0) + (𝑌3 − 𝑌0) 2
𝐿0 = = 443,676
𝑋2−𝑋0 𝑋1−𝑋0 2,15790
𝑎𝑡𝑎𝑛 | 𝑌2−𝑌0 | − (360° − 𝑎𝑡𝑎𝑛 | 𝑌1−𝑌0 |) + 360°
[0,43038]
𝑋3−𝑋0 𝑋2−𝑋0
[ 𝑎𝑡𝑎𝑛 | 𝑌3−𝑌0 | − 𝑎𝑡𝑎𝑛 | 𝑌2−𝑌0 | ]

2.29. Cálculo do vetor L

0,012
−0,031
𝐿 = 𝐿𝑜 − 𝐿𝑏 = 0,175
0,00212
[0,00002]
 183

2.30. Cálculo da correção CX

Cálculo da correção nas coordenadas aproximadas (CX)
𝐶𝑋 = −(𝐴𝑇 . 𝑃. 𝐴)−1 . 𝐴𝑇 . 𝑃. 𝐿
−0,000
𝐶𝑋 = | |
0,000
2.31. Coordenadas corrigidas ou ajustadas (Xa)
𝑋𝑎 = 𝑋0 + 𝐶𝑋
1065,277
𝑋𝑎 = | |
825,199

2.32. Cálculo do vetor dos resíduos

0,012
−0,031
𝑉 = 𝐴. 𝐶𝑋 + 𝐿 = 0,175
0°07′18"
[0°00′05"]

𝑉̅ = 0,032 𝜎𝑉 = 0,082

2.33. Cálculo da variância a posteriori

 184

𝑉 𝑇 𝑃𝑉 𝑉 𝑇 𝑃𝑉
𝜎̂02 = = = 0,672
𝑛−𝑢 3

Agora a 𝜎̂02 ≅ 𝜎02

2𝑐𝑎𝑙𝑐 = 𝑉 𝑇 𝑃𝑉 = 2,017
0,216 ≤ 0,672 ≤ 9,348 -> 𝐻0 é aceita

2.34. MVC de 𝑿𝒂

𝑋𝑎 = 𝜎̂02 (𝐴𝑇 𝑃𝐴)−1

0,0013 0,0005
𝑋𝑎 = [ ]
0,0005 0,0054

Sigma com 95% de nível de confiança

𝜎95 = 1,96

𝜎𝑋0 = 1,96. √𝑋𝑎1,1 = ±0,070

 185

𝜎𝑌0 = 1,96. √𝑋𝑎2,2 = ±0,144

𝑋0 1065,277 𝑚 0,070 𝑚
[ ]=[ ]±[ ]
𝑌0 825,199 𝑚 0,144 𝑚

2.35. Vetor das observações ajustadas (𝑳𝒂)

𝐿𝑎 = 𝐿𝑏 + 𝑉

244,524 244,524
321,539 321,539
𝐿𝑎 = 443,676 = 443,676
2,15790 123°38′19"
[0,43038] [ 24°39′32" ]

2.36. MVC da 𝑳𝒂

𝐿𝑎 = 𝐴. Σ𝑋𝑎. 𝐴𝑇
 186

𝜎𝐿𝑎𝑖 = 1,96. √𝐿𝑎𝑖,𝑖 𝑖 = 1…6

0,077
0,106
𝜎𝐿𝑎 = 0,076
0,00092
[0,00010]
 187

244,524 0,077 𝑚
321,539 0,106 𝑚
𝐿𝑎 = 443,676 ± 0,076 𝑚
123°38′19" 0°03′10"
[ 24°39′32" ] [0°00′21"]