UNIVERSIDAD DE LAS FUERZAS ARMADAS. “E.S.P.
E”
DEPARTAMENTO DE LAS CIENCIAS DE LA TIERRA Y DE LA CONSTRUCCIÓN
Asignatura:
CÁLCULO DE COMPENSACIÓN
NRC:
2294
Docente:
Ing. Alexander Robayo
Estudiante:
Katty Asimbaya
Tema:
Ajuste por el método paramétrico
Fecha de entrega del informe: 13/11/2019
Sangolquí-Ecuador
CÁLCULO DE COMPENSACIÓN
Tema: Ajuste por el método paramétrico
I. Objetivo general:
Calcular el valor verdadero de las coordenadas de un punto desconocido a través de
observaciones medidas y su respectivo ajuste a través del método paramétrico.
II. Objetivos específicos:
Determinar el verdadero valor de un punto P y de las mediciones verdaderas
obtenidas de las distancias hacia ese punto, mediante iteraciones.
Analizar las iteraciones para determinar mejores resultados y variabilidad
presentada en cada caso.
Determinar las precisiones para las matrices V ( residual), Xa ( parámetros
ajustados), La ( observaciones de los valores ajustados), mediante la matriz
varianza y covarianza (∑ 𝐿𝑎, ∑ 𝑋𝑎, ∑ 𝑉) para el ejercicio 2.
III. Marco teórico.
Método Paramétrico de ajuste
Características:
Utiliza incógnitas para su desarrollo en base a un modelo matemático.
Es necesario conocer un valor estimado.
Si el modelo es lineal no es posible realizar iteraciones, para los demás grados, sí.
El número de ecuaciones es igual al número de observaciones.
La metodología de resolución impartida en clases pretende encontrar un valor cercano al
valor verdadero mediante ecuaciones de condición las cuales garantizan un resultado más
real si estás son mayores que las incógnitas, esto hará posible mejorar la precisión y se
podrá tener un control fiable de las mediciones en general.
Página 2 de 17
CÁLCULO DE COMPENSACIÓN
Procedimiento:
Paso 1: Se obtiene las ecuaciones del modelo matemático.
Paso 2: Se establecen los valores estimados de Xo.
Paso 3: Calculo de la matriz Lo en función del modelo matemático establecido
anteriormente.
Paso 4: Calculo de la matriz L.
Paso 5: Calculo de la matriz de derivadas parciales en base a las ecuaciones del modelo
matemático.
Paso 6: Calculo de la matriz de pesos:
Paso 7: Calculo de la matriz de corrección de los parámetros X.
Paso 8: Calculo de la matriz de residuos V.
Paso 9: Calculo de la matriz de observaciones ajustadas La.
Paso 10: Calculo de la matriz de parámetros ajustados Xa.
Paso11: Calculo del Indicativo de la precisión (Varianza posteriori).
Paso 12: Realizar iteraciones en caso de que el modelo no sea lineal.
IV. Cálculos:
Ejercicio Nº 1:
Se conocen las coordenadas de los vértices P1, P2, P3, P4, se conoce además el ángulo
del vértice P1P2, además en la libreta de campo se tiene lb1, lb2, lb3, lb4, lb5. Calcular
las coordenadas ajustadas (XPa, YPa).
Página 3 de 17
CÁLCULO DE COMPENSACIÓN
Tabla Nº 1. Coordenadas de los vértices.
Coordenadas conocidas
Punto X(m) Y(m)
P1 842,281 925,523
P2 1337,544 996,249
P3 1831,727 723,962
P4 840,408 658,345
Fuente: Katty Asimbaya (2019)
Tabla Nº 2 Mediciones obtenidas de las distancias al punto P.
Distancias observadas
Lbi Distancias Precisión
Lb1 244,512 m 0,012 m
Lb2 321,57 m 0,016 m
Lb3 773,154 m 0,038 m
Lb4 279,992 m 0,014 m
Lb5 2,1578 2 ’’
Fuente: Katty Asimbaya (2019)
Matriz Lb:
244,512
321,570
Lb= 773,154
279,992
2,1578
Paso 1: Obtención del modelo matemático y determinación del número de ecuaciones e
incógnitas del ejercicio.
Tabla Nº 3. Fórmulas del modelo matemático.
Lo Fórmulas
La1 √(XP1 − XP0 )2 + (YP1 − YP0 )2
La2 √(XP2 − XP0 )2 + (YP2 − YP0 )2
La3 √(XP3 − XP0 )2 + (YP3 − YP0 )2
La4 √(XP4 − XP0 )2 + (YP4 − YP0 )2
XP2 − XP0 XP1 − XP0
La5 tan−1 ( ) − tan−1 ( )
YP2 − YP0 YP1 − YP0
Fuente: Katty Asimbaya (2019)
Página 4 de 17
CÁLCULO DE COMPENSACIÓN
Número de ecuaciones: Cinco (La1 , La2 , La3 , La4 , La5 )
Número de incógnitas: Dos (XP0 , YP0)
Paso 2: Obtención de los valores estimados de los parámetros 𝑿𝒐.
Valor estimado
XPo 1070
Xo
YPo 900
Paso 3: Cálculo de la matriz de observaciones ajustadas ( Lo).
Se reemplazan las coordenadas de la tabla Nº 1 en las ecuaciones del modelo
matemático, para calcular La1 , La2 , La3 , La4 , La5 .
𝐿0 = 𝑓(𝑋0 )
La1
La2
𝐿0 = La3
La4
La5
229,144859
284,330199
Lo= 781,80394
333,331105
2,684642
Paso 4: Cálculo de la matriz de derivadas parciales (A) en base a las ecuaciones del
modelo matemático de la tabla Nº 3.
𝝏𝑭
𝐀= |
𝝏𝑿𝒂 𝑿𝒐
Página 5 de 17
CÁLCULO DE COMPENSACIÓN
Tabla Nº 4. Derivadas parciales del modelo matemático.
𝜕𝐹 𝜕𝐹
𝜕𝑋𝑃𝑎 𝜕𝑌𝑃𝑎
𝜕𝐹1 𝑋𝑃1 − 𝑋𝑃𝑎 𝜕𝐹1 𝑌𝑃1 − 𝑌𝑃𝑎
= −( ) = −( )
𝜕𝑋𝑃𝑎 √(𝑋𝑃1 − 𝑋𝑃𝑎 )2 + (𝑌𝑃1 − 𝑌𝑃𝑎 )2 𝜕𝑌𝑃𝑎 √(𝑋𝑃1 − 𝑋𝑃𝑎 )2 + (𝑌𝑃1 − 𝑌𝑃𝑎 )2
𝜕𝐹2 𝑋𝑃2 − 𝑋𝑃𝑎 𝜕𝐹2 𝑌𝑃2 − 𝑌𝑃𝑎
= −( ) = −( )
𝜕𝑋𝑃𝑎 √(𝑋𝑃2 − 𝑋𝑃𝑎 )2 + (𝑌𝑃2 − 𝑌𝑃𝑎 )2 𝜕𝑌𝑃𝑎 √(𝑋𝑃2 − 𝑋𝑃𝑎 )2 + (𝑌𝑃2 − 𝑌𝑃𝑎 )2
𝜕𝐹3 𝑋𝑃3 − 𝑋𝑃𝑎 𝜕𝐹3 𝑌𝑃3 − 𝑌𝑃𝑎
= −( ) = −( )
𝜕𝑋𝑃𝑎 √(𝑋𝑃3 − 𝑋𝑃𝑎 )2 + (𝑌𝑃3 − 𝑌𝑃𝑎 )2 𝜕𝑌𝑃𝑎 √(𝑋𝑃3 − 𝑋𝑃𝑎 )2 + (𝑌𝑃3 − 𝑌𝑃𝑎 )2
𝜕𝐹4 𝑋𝑃4 − 𝑋𝑃𝑎 𝜕𝐹4 𝑌𝑃4 − 𝑌𝑃𝑎
= −( ) = −( )
𝜕𝑋𝑃𝑎 √(𝑋𝑃4 − 𝑋𝑃𝑎 )2 + (𝑌𝑃4 − 𝑌𝑃𝑎 )2 𝜕𝑌𝑃𝑎 √(𝑋𝑃4 − 𝑋𝑃𝑎 )2 + (𝑌𝑃4 − 𝑌𝑃𝑎 )2
1 1 (𝑋𝑃2 − 𝑋𝑃𝑎 ) (𝑋𝑃1 − 𝑋𝑃𝑎 )
𝜕𝐹5 (𝑌𝑃2 − 𝑌𝑃𝑎 ) (𝑌𝑃1 − 𝑌𝑃𝑎 ) 𝜕𝐹5 (𝑌𝑃2 − 𝑌𝑃𝑎 )2 (𝑌𝑃1 − 𝑌𝑃𝑎 )2
= −( )+( ) =( )−( )
𝜕𝑋𝑃𝑎 𝑋𝑃 − 𝑋𝑃𝑎 2 𝑋𝑃 − 𝑋𝑃𝑎 2 𝜕𝑌𝑃𝑎 𝑋𝑃2 − 𝑋𝑃𝑎 2 𝑋𝑃1 − 𝑋𝑃𝑎 2
1+( 2 ) 1+( 1 ) 1+( ) 1+( )
𝑌𝑃2 − 𝑌𝑃𝑎 𝑌𝑃1 − 𝑌𝑃𝑎 𝑌𝑃2 − 𝑌𝑃𝑎 𝑌𝑃1 − 𝑌𝑃𝑎
Fuente: Katty Asimbaya (2019)
0,993777 -0,111384
-0,940962 -0,338511
A = -0,974320 0,225169
0,688781 0,724970
-0,000704 0,007646
Paso 5: Calculo de la matriz de pesos P.
𝟏
𝑷=
𝝈𝟐
6944,4444 0 0 0 0
0 3906,2500 0 0 0
P= 0 0 692,5208 0 0
0 0 0 5102,0408 0
0 0 0 0 0,2500
Paso 6: Calculo de la matriz L.
𝐿 = 𝐿0 − 𝐿𝑏
Página 6 de 17
CÁLCULO DE COMPENSACIÓN
-15,367
-37,240
L= 8,650
53,339
0,52683
Paso 7: Cálculo de la matriz N-1.
𝑁 −1 = (𝐴𝑇 . 𝑃. 𝐴)−1
9,2094E-05 -8,1353E-05
N-1 =
-8,1353E-05 0,000379517
Paso 8: Calculo de la matriz U.
𝑈 = 𝐴𝑇 . 𝑃. 𝐿
212434,814
U=
259769,888
Paso 9: Cálculo del vector de corrección (X).
X= -(N-1) .U
1,5690
X=
-81,3049
Paso 10: Cálculo de la matriz Xa, parámetros ajustados.
Xa=Xo+X
1071,5690
Xa=
818,6951
Paso 11: Cálculo de la matriz V, de residuales.
V=Ax+L
-4,752
-11,194
V= -11,186
-4,524
-0,096
Paso 12: Cálculo de la matriz La, observaciones de los valores ajustados.
Página 7 de 17
CÁLCULO DE COMPENSACIÓN
𝑳𝒂 = 𝑳𝒃 + 𝑽
239,7601
310,3764
La= 761,9679
275,4682
2,0619
Paso 13: Cálculo de la varianza aposteriori.
𝑉 𝑇 .𝑉
𝜎̂02 =
𝑛−𝑢
(# de observaciones) n=5
(# de incógnitas) u=2
𝜎̂02 = 837308,4111
3
2
𝜎̂0 =279102,8037
Paso 14: Realizar iteraciones.
Como las ecuaciones planteadas no son lineales es necesario realizar iteraciones, en este
caso se realizaron cuatro iteraciones hasta que los valores de la matriz de parámetros
ajustados Xa y de la matriz de observaciones ajustadas La, fueran las mismas.
Ejercicio Nº 2:
Se conocen las coordenadas de los vértices P1, P2, P3, P4, se conoce además el ángulo
del vértice P1P2, además en la libreta de campo se tiene lb1, lb2, lb3, lb4, lb5. Calcular
las coordenadas ajustadas (XPa, YPa).
Ver tabla Nº 1 y Nº 2.
Datos estimados
Valor estimado
XPo 1065.2 m
Xo
YPo 825.2 m
Página 8 de 17
CÁLCULO DE COMPENSACIÓN
Matriz Lb.
244,512
321,570
Lb= 773,154
279,992
2,1578
Paso 1: Obtención del modelo matemático y determinación del número de ecuaciones e
incógnitas del ejercicio.
Ver tabla Nº 3.
Número de ecuaciones: Cinco (La1 , La2 , La3 , La4 , La5 )
Número de incógnitas: Dos (XP0 , YP0)
Paso 2: Cálculo de la matriz de observaciones ajustadas ( Lo).
Se reemplazan las coordenadas de la tabla Nº 1 en las ecuaciones del modelo
matemático, para calcular La1 , La2 , La3 , La4 , La5 .
𝐿0 = 𝑓(𝑋0 )
La1
La2
𝐿0 = La3
La4
La5
244,453646
321,603816
Lo= 773,183532
279,950057
2,15790508
Paso 3: Cálculo de la matriz de derivadas parciales (A) en base a las ecuaciones del
modelo matemático de la tabla Nº 3.
Ver tabla Nº 4.
Página 9 de 17
CÁLCULO DE COMPENSACIÓN
0,91190704 -0,41039682
-0,84683075 -0,53186247
A = -0,99139075 0,13093657
0,8029718 0,59601703
2,5051E-05 0,00636354
Paso 4: Calculo de la matriz de pesos P.
6944,44444 0 0 0 0
0 3906,25 0 0 0
P= 0 0 692,5207756 0 0
0 0 0 5102,040816 0
0 0 0 0 0,25
Paso 5: Calculo de la matriz L.
𝐿 = 𝐿0 − 𝐿𝑏
-0,05835431
0,03381642
L= 0,02953214
-0,04194322
0,00008956
Paso 6: Cálculo de la matriz N-1.
𝑁 −1 = (𝐴𝑇 . 𝑃. 𝐴)−1
8,34142E-05 -3,07761E-05
N-1 =
-3,07761E-05 0,000255322
Paso 7: Calculo de la matriz U.
𝑈 = 𝐴𝑇 . 𝑃. 𝐿
-673,510416
U=
-28,8154986
Paso 8: Cálculo de la matriz X.
X= -(N-1) .U
0,055293512
X=
-0,01337077
Paso 9: Cálculo de la matriz Xa.
Página 10 de 17
CÁLCULO DE COMPENSACIÓN
Xa=Xo+X
1065,255294
Xa=
825,1866292
Paso 10: Cálculo de la matriz V.
V=Ax+L
-0,0024444
-0,0058964
V= -0,0270361
-0,0055133
0,0000059
Paso 11: Cálculo de la matriz La.
244,50955555
321,56410359
La= 773,12696394
279,98648670
2,15782138
Paso 12: Cálculo de la varianza posteriori.
𝑉 𝑇 .𝑉
𝜎̂02 = 𝑛−𝑢
(# de observaciones) n=5
(# de incógnitas) u=2
𝜎̂02 = 0,83858767
3
𝜎̂02 = 0,279529223
Ejercicio Nº 3:
Hallar la matriz varianza y covarianza (∑ 𝐿𝑎, ∑ 𝑋𝑎, ∑ 𝑉) de acuerdo a los valores de los
resultados en el ejercicio 2.
Matriz varianza y covarianza de La.
∑ 𝐿𝑎 = 𝜎̂2
0 (𝐴 ∗ 𝑁
−1
∗ 𝐴𝑇 )
Página 11 de 17
CÁLCULO DE COMPENSACIÓN
3,5526E-12 -1,16174E-13 -2,76349E-12 -2,08339E-13 -2,21241E-14
-1,16174E-13 2,73751E-12 1,0341E-12 -2,85942E-12 -1,83605E-14
∑ 𝐿𝑎 = -2,76349E-12 1,0341E-12 2,47546E-12 -8,2742E-13 1,06217E-14
-2,08339E-13 -2,85942E-12 -8,2742E-13 3,01739E-12 2,13065E-14
-2,21241E-14 -1,83605E-14 1,06217E-14 2,13065E-14 2,71004E-16
Matriz varianza y covarianza de Xa.
∑ 𝑥𝑎 = 𝜎̂2
0 ∗𝑁
−1
2,18845E-12 -8,07647E-13
∑ 𝑋𝑎 =
-8,07647E-13 6,69979E-12
Matriz varianza y covarianza de V.
∑ 𝑉 = 𝜎̂2
0 (𝐴 ∗ 𝑁
−1
∗ 𝐴𝑇 − 𝑃−1 )
2,25484E-13 1,16174E-13 2,76349E-12 2,08339E-13 2,21241E-14
1,16174E-13 3,97908E-12 1,0341E-12 2,85942E-12 1,83605E-14
∑ 𝑉 = 2,76349E-12 1,0341E-12 3,54103E-11 8,2742E-13 1,06217E-14
2,08339E-13 2,85942E-12 8,2742E-13 2,125E-12 2,13065E-14
2,21241E-14 1,83605E-14 1,06217E-14 2,13065E-14 1,04947E-07
V. Resultados:
En la tabla Nº 5 se puede apreciar como los valores para los parámetros ajustados así
como para las observaciones de los parámetros ajustados, en la tercera y cuarta iteración
ya no cambian.
Tabla Nº 5. Resultado de las iteraciones para el ejercicio Nº 1.
Sin Primera Segunda Tercera Cuarta
iteración iteración iteración iteración iteración
Xpo 1071,569002 1064,1710661 1064,067252 1064,067250 1064,067250
Parámetros
ajustados Xa
Ypo 818,6951463 831,4522763 831,4648379 831,4648357 831,4648357
Página 12 de 17
CÁLCULO DE COMPENSACIÓN
Lb1 239,7601328 240,8594618 240,9067838 240,9067846 240,9067846
Lb2 310,3764428 318,8637831 319,2856837 319,2856930 319,2856930
Observaciones
de los
La
parámetros Lb3 761,9678984 774,9569936 775,1505333 775,1505349 775,1505349
ajustados.
Lb4 275,4682359 282,5243125 282,8319885 282,8319952 282,8319952
Lb5 2,061856659 2,1963518 2,198195568 2,19819545 2,19819545
Fuente: Katty Asimbaya
En la tabla Nº 6 se puede apreciar como los valores para los parámetros ajustados así
como para las observaciones de los parámetros ajustados, en primera y segunda iteración
ya no cambian.
Tabla Nº 6. Resultado de las iteraciones para el ejercicio Nº 2.
Sin Primera Segunda
iteración iteración iteración
Xpo 1065,255294 1065,2552931 1065,255293
Parámetros
ajustados Xa
Ypo 825,1866292 825,1866283 825,1866283
Lb1 244,5095556 244,5095558 244,5095558
Lb2 321,5641036 321,564107 321,564107
Observaciones de los
parámetros La
Lb3 773,1269639 773,1269643 773,1269643
ajustados.
Lb4 279,9864867 279,9864892 279,9864892
Lb5 2,157821383 2,157821352 2,157821352
Fuente: Katty Asimbaya
Página 13 de 17
CÁLCULO DE COMPENSACIÓN
Precisiones para las matrices V, Xa y La, obtenidas de la matriz varianza y covarianza
(∑ 𝐿𝑎, ∑ 𝑋𝑎, ∑ 𝑉) para el ejercicio N º2.
1,97897E-07 ± 4,74851E-07
3,43017E-06 ± 1,99476E-06
V= 3,38126E-07 ± 5,95066E-06
2,51991E-06 ± 1,45774E-06
-3,12214E-08 ± 0,000323955
1065,255293 ± 1,47934E-06
Xa=
825,1866283 ± 2,5884E-06
244,50955575 ± 1,88483E-06
321,56410702 ± 1,65454E-06
La= 773,12696428 ± 1,57336E-06
279,98648922 ± 1,73706E-06
2,15782135 ± 1,64622E-08
VI. Conclusiones:
Los valores estimados en el ejercicio Nº 1 fueron 𝑋𝑝𝑜 (1070 ,900), para su
correspondiente ajuste, se realizó cuatro iteraciones que se pueden observar en
la Tabla N º 5.
Como resultado de la cuarta iteración, los valores de Xa y La, fueron:
Tabla N º 7. Resultados del valor verdadero para las coordenadas del punto Pa,
así como para las distancia entre las de los puntos, obtenidos en la cuarta
iteración. .
Cuarta
iteración
Xa Xpo 1064,067250m
Parámetros
ajustados Ypo 831,4648357m
Lb1 240,9067846m
Lb2 319,2856930m
Observaciones de los
parámetros ajustados.
La Lb3 775,1505349m
Página 14 de 17
CÁLCULO DE COMPENSACIÓN
Lb4 282,8319952m
Lb5 2,19819545mm
Fuente: Katty Asimbaya
Los valores estimados en el ejercicio Nº 2 fueron 𝑋𝑝𝑜 (1065.2 ,825.2), para
su correspondiente ajuste, se realizaron dos iteraciones.
Como resultado de la segunda iteración, los valores de Xa y La, fueron:
Tabla N º 8. Resultados del valor verdadero para las coordenadas del punto Pa,
así como para las distancia entre las de los puntos, obtenidos en la segunda
iteración. .
Segunda
iteración
Xpo 1065,255293 m
Parámetros Xa
ajustados Ypo 825,1866283m
Lb1 244,5095558m
Lb2 321,564107m
Observaciones de La Lb3 773,1269643m
los parámetros
ajustados.
Lb4 279,9864892m
Lb5 2,157821352m
Fuente: Katty Asimbaya
Los valores verdaderos se tomaron en cuenta dependiendo del numero de
iteraciones, para saber cuantas se debe realizar, se ha realizado la diferencia
entre las mismas hasta obtener un valor de 0.
Para el ejercicio 1:
Página 15 de 17
CÁLCULO DE COMPENSACIÓN
Tabla N º 9. Diferencias entre iteraciones hasta encontrar el valor verdadero de
Xo y La, es decir que la diferencia entre los valores dados por las iteraciones,
sea 0.
Diferencias
Xa Xpo= S. I. - P.I. P.I. -S.I. S.I. -T.I. T.I. -C.I.
Ypo= 7,39793602 0,103814448 1,93025E-06 0,000000000
-12,75713 -0,012561547 2,21167E-06 0,000000000
La lb1 -1,09932898 -0,047321959 -8,26664E-07 0,000000000
lb2 -8,48734031 -0,421900555 -9,2804E-06 0,000000000
lb3 -12,9890951 -0,193539779 -1,60734E-06 0,000000000
lb4 -7,05607666 -0,307676015 -6,6593E-06 0,000000000
lb5 -0,13449513 -0,001843785 1,16555E-07 0,000000000
Fuente: Katty Asimbaya (2019)
SI : Sin iteración
PI: Primera iteración
SI: Segunda iteración
TI: Tercera iteración
CI: Cuarta iteración
Tabla N º 10. Diferencias entre iteraciones hasta encontrar el valor verdadero de
Xo y La, es decir que la diferencia entre los valores dados por las iteraciones,
sea 0.
Diferencias
Xa S. I. - P.I. P.I. -S.I.
Xpo= 4,3679E-07 0
Ypo= 9,0357E-07 0
La lb1 -1,979E-07 0
lb2 -3,4302E-06 0
lb3 -3,3813E-07 0
lb4 -2,5199E-06 0
lb5 3,1221E-08 0
Fuente: Katty Asimbaya (2019)
SI : Sin iteración
PI: Primera iteración
SI: Segunda iteración
Página 16 de 17
CÁLCULO DE COMPENSACIÓN
El valor de la varianza posteriori depende de los grados de libertad en función del
número de observaciones e incógnitas, para los dos ejercicios cinco y dos.
Tabla 9. Valor varianza posteriori
𝜎𝑜^2 Ejercicio 1 Ejercicio 2
279102.804 0.279529
Fuente: Katty Asimbaya (2019)
[Link]:
Se debe plantear correctamente el modelo matemático del ejercicio, ya que si
aumentamos el numero de incógnitas de manera incorrecta el procedimiento se
complicaría.
Para realizar un buen ajuste de datos mediante el método paramétrico, se
recomienda tener una abundancia de observaciones, aumentando así, el nivel de
confianza para obtener una precisión correcta.
Página 17 de 17
