Computos para El Ajuste
Computos para El Ajuste
Computos para El Ajuste
DEPARTAMENTO DE GEODESIA
LEVANTAMIENTOS DE CAMPO
INDICE
1. INTRODUCCION A LAS MATRICES
Introduccin. Definicin de una matriz. Tamao o dimensiones de una matriz.
Tipos de matrices. Igualdad de matrices. Suma o resta matricial. Multiplicacin de
una matriz por un escalar. Multiplicacin de una matriz por otra. Tarea
autodidctica.
2. SOLUCION DE ECUACIONES POR METODOS MATRICIALES
Introduccin. La inversa de una matriz. Inversa de una matriz de 2x2. Inversas por
adjuntos. Inversas por transformaciones de fi1as. Ejemplos de problemas. Tarea
autodidctica.
3. ELIMINACION GAUSIANA; TEOREMA DE TAYLOR
Eliminacin gausiana: Introduccin; Explicacin del procedimiento; Ejemplo
numrico. Linearizacin de ecuaciones no lineales de la serie de Taylor:
Introduccin; Aproximaciones de la serie de. Taylor; Ejemplo numrico.
Linearizacin de ecuaciones de la observacin de distancias. Linearizacin de
ecuaciones de la observacin de ngulos. Tarea autodidctica.
4. TEORIA GENERAL DE ERRORES
Clasificacin de errores. Definicin de trminos. Teora de la probabilidad.
Representacin grafica de distribuciones de los errores randomizados/al azar.
Ejemplo numrico del histograma. Polgono de frecuencia. Propiedades de la
curva de distribucin normal. Mediciones de precisin. El error probable de 50 por
ciento. El error probable de 90 por ciento. Ejemplo numrico de clculos de
precisin. Trazado de una curva de distribucin normal. Tarea autodidctica.
5. PROPAGACION DE ERRORES FORTUITOS
Derivacin de la ecuacin bsica. Funciones especficas frecuentemente
encontradas. Ejemplos numricos. Tarea autodidctica.
6. PESOS DE LAS OBSERVACIONES
Generalidades. Media ponderada (peso promedio). Relacin entre el peso y el
error estndar. Pesos en la nivelacin y medicin de ngulos. Ejemplo numrico.
Errores estndar de las observaciones ponderadas. Ejemplos numricos. Tarea
autodidctica.
Numero
de
ecuaciones
Introduccin. Encajado de una lnea recta a los datos. Encajado de una parbola a
los datos. Encajado de un polinomio a los datos. Tarea autodidctica.
APENDICES
LA CURVA DE DISTRIBUCION NORMAL DEL ERROR
DE
MINIMOS
CUADRADOS
TOR
CAPITULO 1
INTRODUCCION A LAS MATRICES
1.1 Introduccin
El uso del lgebra matricial brinda por lo menos dos importantes ventajas al
matemtico, ingeniero o agrimensor modernos: (1) proporciona los medios para
reducir sistemas complicados de ecuaciones a expresiones sencillas que pueden
ms fcilmente visualizarse y manejarse, y (2) proporciona un mtodo sistemtico
de tratamiento matemtico que esta muy bien adaptado para las actuales
computadoras digitales electrnicas de alta velocidad. Hoy da se tienen ms
y ms problemas en agrimensura, geodesia y fotogrametra que requieren la
solucin de .grandes sistemas de ecuaciones.
En este texto estudiaremos en
especial los cmputos de ajustes que involucran observaciones redundantes y
satisfaccin de ciertas condiciones.
Tales observaciones y condicin
frecuentemente resultan en enormes sistemas de ecuaciones.
Cuando esas ecuaciones se resuelven conforme al mtodo de mnimos la solucin
rinde los mejores clculos para las observaciones ajustadas y los parmetros
desconocidos. La anotacin matricial es la ms adecuada para los cmputos de
mnimos cuadrados y, por ende se emplea a fin de facilitar el manejo implicado
en el anlisis y solucin de esos sistemas ecuaciones.
1.2 Definicin de una Matriz
Una matriz es un juego de nmeros o smbolos arreglado en una ordenacin
cuadrada o rectangular de filas "m" y columnas "n". El aarreglo es tal que se
pueden efectuar ciertas operaciones matemticas definidas de una manera
sistemtica y eficaz.
Como ejemplo de una representacin matricial, considrese el siguiente sistema
de tres ecuaciones "lineales involucrando a tres incgnitas:
a11 x1 + a12 x 2 + a13 x 3 = c1
a 21 x1 + a 22 x 2 + a 23 x 3 = c 2
a 31 x1 + a 32 x 2 + a 33 x 3 = c3
a11 a 12 a 13 x 1
c1
a 21 a 22 a 23 x 2 = c 2
a 31 a 32 a 33 x 3
c3
Esto a su vez puede representarse en la notacin matricial compacta como:
A X = C . . . . . . . . . . . . (Lc)
Se usan letras maysculas (algunas veces subrayadas) como (A, X,C) para indicar
una matriz o una ordenacin de nmeros o smbolos.
Hablamos de la posicin de un elemento de matriz con subscritos ij, jk, etc., en
donde la primera letra representa la fila y la segunda letra representa la columna.
1.4 Tipo de matrices (se puede usar una variedad de smbolos para designar
las matrices abajo ilustradas).
1.4.1 Matriz columnar (el numero de filas puede ser cualquier entero, mas
el numero de columnas es uno (1)).
3
[A] = (A) = 2
3X1
3X1
5
3X1
3X1
1.4.2 Matriz por filas (El numero de columnas puede ser cualquier entero de filas
es uno (1))
A <A> A A -6 4 2
3X1
1X3
1X 3
1X 3
1X3
A A
2X 3
2X 3
3 2 6
=
5 4 1
2
A A
3X 3
4 2 -5
-7 3 4
6
-1
9
3
3X 3
6
3
5
3
A
3X3
7 0 0
0 -3 0
0
0
6
3
3X3
[I]
4X4
1 0 0 0
0 1 0 0
0 0 1 0
0 0 0 1
4
Restarse trmino por trmino. Si dos matrices son de igual dimensin. Se dice que
son conformes para la suma o resta.
A +B = C
2X 3
2X3
2X3
8 8 5
7 3 -1
1 5 6
2 -5 6 + -4 -2 8 = -2 -7 9
Suponiendo que las dos matrices A y B sean conformes para la suma y la resta,
entonces lo siguiente sera verdad:
a) A + B = B + A
(Ley conmutativa)
b) A + (B+C) = (A+B) + C
(Ley asociativa)
KA +AK = C
4X 2
4X2
4X2
3
2
-1
6
12 -4
8 24
16 28
20 12
[A]
B = C
B = C
B = C
4X 2
1X 3
3X 3
2X3
4X3
3X1
1X1
3X1
3X1
2X 3
6X 2
4X2
6X3
2X 3
a 11 a 12 a 13
a 21 a 22 a 23
B = C
3X2
2X2
b 11 b 12
b 21 b 22
b 31 b 32
c11 c12
c c
21 22
4 4
1 2 3
6
2
4 2 7
5 3
2X 3
c11 c12
c c
21 22
2X 2
31 21
63 57
2X 2
3X 2
2X 2
A = A
2X2
2X2
1 0
0 1
5 6
5 6
7 8 = 7 8
1.8.3 Suponiendo que las matices A,B y C sean conformes para la multiplicacin,
y en el orden indicado, luego lo que sigue es verdadero:
a) A(B +C) = AB + AC (Primera ley distributiva)
b) (A + B) C = AC + BC (Segunda ley distributiva)
c) A(BC) = (AB) C
(ley asociativa)
Tambin se declaran las siguientes advertencias:
d) AB no es generalmente igual a BA y BA quiz no sea conforme.
e) Si AB = 0, ni A ni B necesariamente = 0
f) Si AB = AC, B ni es necesariamente = c
Permtase que:
-1 -1
B=
2 2
2 1
A=
1 2
1 1
C=
1 1
Ejemplo de (a):
A
1
2
AB
0 0
3 3
2
1
Ejemplo de (b):
Ejemplo de (C):
(A + B)
1 0
3 4
AB
3 3
3 3
(B + C)
0
0
3 3
AC
3 3
3 3
C
1 1
1 1
AC
-2 -2
+
4 4
BC
3 3
6 6
3 3
6 6
1 1
7 7
1 1
7 7
2
1
1
2
AB
-2 -2
4 4
0 0
3 3
Ejemplo de (d):
AB
0 0
3 3
1 1
1 1
0 0
6 6
0 0
6 6
BA
-3 -3
6 6
Ejemplo de (e):
Permita que
2 2
A=
1 1
-1 -2
B=
1 2
Luego AB = 0, pero ni A ni B = 0
Ejemplo de (f):
2 2
Permita que A =
;
1 1
1 2
B=
1 3
AB
4 10
-2 -5
y
=
4 10
C=
-2 -5
AC pero B C
4 10
-2 -5
nxp
sigue:
(Ntese que son conformes para la multiplicacin)
mxn
mxn
mxn
k p
k
=
C IK =
ij
j 1
X B JK ; i=1 k=1
Paso 2.
Aumentar k sucesivamente por 1 colocando el resultado en C IK continuar aumento
de k y repetir paso 2 hasta k=p.
Paso 3.
Aumentar i de una a dos y repetir paso 1 y 2 en su totalidad. Al terminar con i=2,
aumentar I a 3 y repetir paso 1 y 2 en su totalidad. Continuar este proceso hasta
i=m. Esto completa la multiplicacin matricial.
En FORTRAN, todo lo anterior se puede expresar sencillamente como:
Do100, I 1, m
Do100, k 1, p
C(I , K ) 0
Do100, J 1, N
100
C ( I , K ) C ( I , K ) A( I , J ) * B ( J , K )
Una nota final sobre multiplicacin: es muy importante especificar el orden como
han de multiplicarse las matrices.
Tarea autodidctica No. 1
1) Supngase que el siguiente sistema de ecuaciones lineal ha de ser
representado en forma matricial como AX=B. formlense las tres matrices.
2 X 1 X 2 5X 3 X 4 5
X 1 X 2 3 X 3 4 X 4 1
3 X 1 6 X 2 2 X 3 X 4 8
2 X 1 2 X 2 2 X 3 3 X 4 2
3 -4
1 2 1
X 1 5
d)
4 0 2
-2
2
A.
CAPITULO 2
SOLUCION DE ECUACIONES POR METODO MATRICIAL
2.1 Introduccin.
Como antes se menciono, una de las ventajas ofrecidas por el algebra matricial es
que se adapta extremadamente bien para su uso en computadoras digitales. Los
sistemas de ecuaciones lineales simultneas de cualquier tamao hasta los ms
grandes pueden programarse para la solucin general de computadoras usando
solo unos cuantos pasos sistemticos. La sencillez de la programacin de
multiplicaciones matriciales, por ejemplo se presento en la seccin 1.8.4.
Para resolver el sistema de ecuaciones usando mtodos matriciales, primero es
necesario definir y computar la inversa de una matriz.
2.2 La inversa de una matriz
2.2.1 Si una matriz cuadrada no es singular (no singular), entonces posee una
matriz inversa. (Ntese que si un sistema de m ecuaciones lineales simultaneas
conteniendo m ecuaciones y m incgnitas se expresa como: CX=B, C la matriz de
coeficientes es una matriz cuadrada de mxm dimensiones).
Considere este referido sistema de ecuaciones lineales:
C X
Definimos
= B . . . . . . . . . . . . . . . (2a)
I
X
C
X
=
X
=
C
1
, luego
B . . . . . . . . . . . . . . . . . (2b)
Por definicin: A
Permtase que A
a b
=
c
d
I . . . . . . . . . . . . . . . (2c)
2
w x
y z
a b
w x 1 0
c d X y z 0 1
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (2d)
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (2e)
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (2f)
. . . . . . . . . . . . . . . . . . . . . . . . . . . . . . (2g)
e I
A y es igual a ad bc
1 aw
cw
; De (2F) y =
b
b
d
1 aw
d
cw
Luego
=
Reduciendo: w =
=
A
b
da bc
b
De (2d) y =
1 cw
ax
; De (2g) z =
b
b
b
1 cx
b
ax
Luego
=
Reduciendo: x =
=
A
d
da bc
b
De (2e) z =
1 by
dy
; De (2f) w =
a
c
c
1 by
c
dy
Luego
=
Reduciendo: y =
=
A
a
ad cb
c
De (2d) w =
bz
1 dz
; De (2g) x =
a
c
a
bz
1 dz
a
Luego
=
Reduciendo: z =
=
A
a
c
ad bc
De (2e) x =
-b
A
d
A
-c
A
1
A
d -b
-c a
c. Ejemplo numrico:
Si:
2 3
4 1 ,
1 1 -3
10 -4 1
A =
1 1 -3 2 3 1 0
10 -4 1 4 1 0 1
Comprobar:
Adjunto de A =
Adjunto de A
Determinante de A A
a 11 a 12 a 13
a 21 a 22 a 23
a 31 a 32 a 33
4 3 2
3 4 1
2 3 4
(1)
Co-factor de a 31 =
(3x1 2x4) = -5
Matriz de cofactores =
13 -10 1
-6 12 6
-5 2 7
Adjunto de A
13 -6 -5
-10 12 2
1 -6 7
13 -6 -5
1
-10 12 2
=
24
1 -6 7
13 -1 -5
24 4 24
-5 1 1
12 2 12
-1 7
1
24 4 24
3X 3
-1
A = I
3X3
3X3
13 -1 -5
24 4 24
-5 1 1
12 2 12
-1 7
1
24 4 24
4 3 2
3 4 1
2 3 4
1 0 0
0 1 0
0 0 1
A =
4 3 2
3 4 1
2 3 4
1 0 0
0 1 0
0 0 1
1
1
o , y colquese en la fila 1
a11
4
a uno (1).
, respectivamente. Esto
a cero.
, respectivamente. Esto
a cero.
3
4
7
4
3
2
1
2
-1
2
son:
4 0 0
-3 1 0
4
-1
0 1
2
, por
1
4
o , y colquese en la fila 2 de
a 22
7
2. Multipliquese la fila 2 de
y colquese en la fila 1 de
3. Multipliquese la fila 2 de
y colquese en la fila 3 de
12
32
1 0 7
0 1 -2
0 0 24
7
4
7
-3
7
1
7
son:
-3
0
7
4
0
7
-6
1
7
Ntese que como resultado de esta segunda serie de tres pasos, la segunda columna de
ha sido hecha para concordar con la segunda de una matriz idntica de 3x3.
1
7
o
, y colquese en la fila 3 de
a 33
24
, por
, respectivamente.
5
2. Multipliquese la fila 3 de
y colquese en la fila 1 de
3. Multipliquese la fila 3 de
13
a 23 o
-2
, rstese de la fila 2 de
7
1 0 0
= I = 0 1 0
0 0 1
A e I por
I , y colquese en la fila 2 de A e I .
-1
23
24
-5
12
1
24
-1
4
1
2
-2
4
son:
-5
24
1
12
7
24
-1
-1
304
302
305
306
303
307
DO 307 K= 1,N
DO 302 K= 1,N
IF (1-K) 304, 302,304
QXX (K,J) = QXX (K,J) / QXX (K,K)
CONTINUE
QXX (K,K) = 1, / QXX (K,K)
DO 307 I = 1,N
IF (I-K) 305, 307, 305
DO 303 J = 1,N
IF (J-K) 306, 303, 306
QXX (I,J) = QXX (I,J) - QXX (I,K) * QXX (K,J)
CONTINUE
QXX (I,K) = - QXX (I,K) * QXX (K,K)
CONTINUE
Fig. 2.1
Valores medidos:
X
X X
X X
125.127
259.60
395.85
1 0 0
1 1 0
1
1
1
3
X1
X 2
X 3
X =
-1
L =
125.27
259.60
395.85
3
1
2
3
1 0 0 1 0 0 1 0 0 1 0 0
1 1 0 0 1 0 0 1 0 -1 1 0
1 1 1 0 0 1 0 0 1 0 -1 1
X =
-1
L =
1 0 0
-1 1 0
0 -1 1
X
X
X
125.27
259.60
395.85
= 125.27
= 134.33
= 136.25
125.27
134.33
136.25
3 -1 -1
-1 3 -1
-1 -1 3
CAPITULO 3
ELIMINACION GAUSIANA; TEOREMA DE TAYLOR
a 21 x1 + a 22 x 2 + ........+a 2n x n +k 2 = 0
1
2
a n1 x1 + a n2 x 2 + ........+a nn x n +k n = 0
a11 x1 + a 12 x 2 + ........+a 1n x n +k 1 = 0
21
11
y se resta de la ecuacin
2 ,
1
que elimina
de
2 .
1
Llmenos a la
ecuacin resultante
y se resta de
3
1
2 . La ecuacin
2
, multiplquese por
41
y se resta de
31
X
X
n 1
yX
n2
n 1
y X n.
n2
n 1
A. solucin delantera
1
2
3
1
1
11
y3 +
=0
2
2
3
48
5y 2 y 3 +
=0
2
2
y1 + 2y 2 +
Sistema de
(n-1)=2
1 X
1
1
=
2
2 1
1
*
X 1 2
-7y + 4y 32 = 0 3 1 X 2 3
1
y 0.3y +4.3=0 2 X
= 2
5
1.9y +1.9=0 3 2 X 7 3
*
Sistema
de n-2=1
B. solucin trasera
y 3 = 1.0
3
de la ecua. 3
de la ecua. 3
de la ecua.
-7y 2 4(1)+32 = 0
2y1 3( 4) 5(1) 21=0
L = F (x,y) .. (3a)
Por el teorema de Taylor la ecuacin puede ser representada as:
En
(3b),
son
la
x , y dx F " x , y dx
x , y = F x , y 1!
2!
F x , y dx
F ' x , y dy F " x , y dy
+ .......... +
F'
L=F
n!
+ .......... +
Fy n
x , y dy
0
1!
2!
............................(3b)
n!
ecuacin
y
y
x x 0 dx
y y dy
. (3c)
x , y = F x , y Lx
0
L
dx dy...............(3d )
y 0
0
Segunda iteracin
F = 2.25 + 2.75 - 2(2.75) + dx + [1 - 4(2.75)] dy = -4
G = (2.25) + (2.75) + 2(2.25) dx + 2(2.75) dy = 8
De las dos soluciones anteriores, dx = -0.25 y dy = -0.64; de lo cual:
x x 0 dx 2.25 0.25 2.00
y y dy 2.75 0.64 2.11
0
Tercera iteracin
F = 2 + 2.11- 2(2.11) + dx + [1 - 4(2.11)] dy = -4
G = (2) + (2.11) + 2(2) dx + 2(2.11) dy = 8
De las dos soluciones anteriores, dx = 0.0 y dy = -0.11; de lo cual:
x x 0 dx 2.00 0.00 2.00
y y dy 2.11 0.11 2.11
0
Cuarta iteracin
F = 2 + 2 - 2(2) + dx + [1 - 4(2)] dy = -4
G = (2) + (2) + 2(2) dx + 2(2) dy = 8
De las dos ecuaciones anteriores, las correcciones dx y dy son cero al centsimo
ms cercano. Por tanto, la serie ha convergido segn la evidencia por las
correcciones cero y los valores X = 2.00 y Y = 2.00 son las incgnitas deseadas.
3.2.4 Linearizacion de ecuaciones de observacin de distancias.
Los dos tipos de observaciones implcitas en levantamientos de control horizontal
son las mediciones angulares obtenidas usando teodolitos, y la medicin de
distancias obtenidas por la medicin con cinta o por medios electrnicos. Al ajustar
las posiciones horizontales por el mtodo de la ecuacin de observacin de
mnimos cuadrados, las ecuaciones se escriben relacionando las cantidades
observadas y sus errores inherentes fortuitos con las coordenadas X y Y mas
probables de los puntos involucrados. Este procedimiento a menudo se conoce
como el mtodo de la variacin de coordenadas. Haciendo referencia a la Fig.
3.1, se puede escribir la siguiente ecuacin de observacin de distancias para
cualquier lnea I-J:
L V
IJ
IJ
X X Y Y
2
....................(3e)
L V
IJ
IJ
X X Y Y
2
F X ,X ,Y ,Y F X ,X ,Y ,Y F
x
I
I0
J0
I0
J0
dxi
F
F
F
dyi
dy j ...................................(3 g )
dx j
y
x j
y
0
i 0
j 0
xi x i dxi
0
yi y
i0
F
En la ecuacin (3g),
es la derivada parcial de F con respecto a
dyi
x i 0
X ,X ,Y ,Y
dX , dX , dY ,dY
I0
J0
I0
J0
son
yi y dyi
y j y dy j
i0
j0
j0
y y
j0
i0
( dyi ) X j0 X
(
dxi
)
ij0
ij0
y y
j
i0
( dyj ).............................(3h)
0
ij0
Donde:
K Lij L
ij
ij0
i0
( dxi )
ij0
Xi
Y
2
j0
Yi
X X Y Y
2
F 1
x i 2
X X Y Y
2
xi
1/ 2
X X
j
-X j X i
X j Xi
Y j Y i
X( 1)
1/ 2
xi
ij0
3.2.5 Linearizacin de ecuaciones de observacin de ngulos con referencia
a la Fig. 3.2 puede escribirse la siguiente ecuacin de observacin de
ngulos:
-1
+
=
tan
jik V jik
y y
X X
j
j
jik
tan -1
y y
X X
k
................(3i)
, el ngulo observado y
jik
el error residual
jik
Donde U
+V U
jik
X ,Y ,X ,Y ,X ,X ............................(3 j)
I
tan
-1
y y
X X
j
j
tan -1
y y
X X
k
U
U X I,Y I,X J,Y J,X ,X k U X I ,Y I ,X ,Y J ,X ,Y k +
dX i
k
k0
0
0
0
0
J0
x i 0
U
U
U
U
U
dyi
dyk ................(3K )
dy j
dX j
dX k
y
x j
y
x
k 0
0
i 0
j 0
k 0
Yi
X i0
Yj
Xj
K jik V jik 0 2 0 2 dX i 0 2
dYi
2
(IJ 0 ) (IK 0 )
(IJ 0 ) (IK 0 )
Y j Yi
X j Xi
Yk Yi
Xi Xk
0 2 0 dX j 0 2 0 dY j 0 2 0 dX k 0 2 0 dYk
(IJ 0 )
(IJ 0 )
(IK 0 )
(IK 0 )
Donde
..(3m)1
X j0 X i0
X k0 X i0
o tan1
ijk
IJ 0 ( X j0 X i0 ) 2 (Y j0 Yi0 ) 2
IK 0 ( X k0 X i0 ) 2 (Yk0 Yi0 ) 2
En la ecuacin (3m), V , K estn medidos en radin. Puesto que es ms comn trabajar en el sistema sexagesimal en
jik
jik
este pas y, ya que las magnitudes de los residuales de ngulos generalmente estn en la gama de segundo, pueden
convertirse las unidades de la ecuacin (3m) a segundos multiplicando el lado derecho de la ecuacin por p (Rho), el
nmero de segundos por radin , que es 206.265
Y j Yi Yk Yi
Xi X j
Yi Y j
X j Xi
Yk Yi
Xi Xk
X i X k0
0
0
0
0
0
K jik V jik 0 20 0
dxi 0
0
dyi 0
dx j 0 2 0 dy j 0
dxk 0
dyk (3m ) *
2
2
2
2
2
2
IK 0
IK 0
IJ 0
IK 0
IK 0
IJ 0
IJ 0
IJ 0
Donde
X j0 X i0
X k0 X i0
0 tan 1
jik
IJ 0 ( X j0 X i0 ) 2 Y j0 Yi0
IK 0 ( X k0 X i0 ) 2 Yk0 Yi0
X 2 Y 2 2 XY Z 4
X Y Z3 4
2 X Y Z 3 23
CAPITULO 4
TEORIA GENERAL DE ERRORES
Clasificacin de errores
Durante el proceso de medicin de cualquiera cantidad fsica, ciertos factores tales
como limitaciones humanas, imperfecciones e inestabilidades instrumentales
harn inexactos los valores medidos.
De hecho, puede declararse
incondicionalmente hablando que cada valor medido contiene errores.
Antes de proceder con un anlisis de errores, ser til clasificar los tipos de
errores que puedan ocurrir en las siguientes categoras:
Equivocacin: Un disparate por parte del observador como resultado de su
descuido o confusin. En general, las equivocaciones no se clasifican como
errores. Slo pueden eliminarse mediante una cuidadosa verificacin para aislar
el disparate y luego corregirlo por medio de nuevas mediciones. Ejemplos de
equivocaciones: 1) aadiendo un pie en una distancia medida cuando se usa una
extensin de cinta o al dejar de registrar una longitud de cinta; 2) equivocaciones
al leer escalas graduadas; y 3) disparates al registrar, escribir 27.55 por 25.75, etc.
Error sistemtico: Error en un valor medido que sigue alguna ley fsica o que
puede expresarse como una funcin matemtica. Si se miden las condiciones
causantes del error, puede calcularse una correccin y se elimina el error
sistemtico. Como ejemplos: 1) la temperatura no es estndar al medir con cinta,
2) un error ndice vernier al medir un ngulo vertical, 3) uso de una mira de nivel
que no tiene la longitud estndar.
...............(4 a )
n 1
n 1 .. (4b)
Donde S = desviacin estndar
S = variacin
n 1 = los grados de libertad
V = la suma de los cuadrados de los residuales.
Teora de la probabilidad
El ajuste de las cantidades medidas que contienen es uno de los problemas
principales en las ciencias exactas. Ciertamente, esto es de importancia global
para los agrimensores, geodesias y fotogrametristas. En los restantes artculos de
este texto, asumiremos que todas las equivocaciones y errores sistemticos han
sido eliminados de los valores medidos, y que solo quedan errores fortuitos. Tales
errores fortuitos siguen las leyes de la probabilidad y, por ende, as sern tratados.
La probabilidad se define como la razn / relacin del numero de veces que
debiera ocurrir un suceso o evento comparado al nmero total de posibilidades.
Por ejemplo, considrese el tiro de un dado, donde hay un sexto de probabilidad
que salga un dos, lo que sencillamente significa que hay seis posibilidades en total
y slo una de estas es un dos. En general, si ocurriera un suceso en forma m, y
dejara de ocurrir en forma n, entonces la probabilidad de ocurrencia sera m /
(m+n) y la probabilidad de fracaso sera
n / (m+n).
La probabilidad se percibe como una fraccin que se extiende en valor del cero al
uno, el cero indicando imposibilidad y el uno indicando certeza. Puesto que un
evento debe ocurrir o dejar de suceder, la suma de las probabilidades de cualquier
evento es uno. Por tanto, si 1 / 6 es la probabilidad de sacar un dos con un tiro de
un dado, entonces (1-1 / 6) 5 / 6 es la probabilidad de que no salga un 2.
En terminologa de probabilidad, un evento compuesto es la ocurrencia simultnea
de dos o ms eventos independientes. La ocurrencia compuesta de errores es el
caso ms frecuente confrontado en la prctica de levantamientos; por ejemplo, la
ocurrencia de errores compuestos en un levantamiento de poligonales que impide
su cierre.
La probabilidad de la ocurrencia simultnea de dos eventos
independientes es el producto de las probabilidades individuales de esos eventos
independientes.
Para ilustrar esta condicin, considrese el simple ejemplo de dos cajas que
contienen una combinacin de bolas rojas y blancas. Si la caja A contiene 4
bolas-una roja y 3 blancas- y si la caja B contiene 5 bolas-2 rojas y 3 blancascul sera la probabilidad de sacar 2 bolas rojas si una bola se saca de cada
caja? El nmero total de posibilidades de pares es 4*5 20, por haber sacado
digamos la bola roja de caja A,
Cualquiera de las 5 bolas en la caja B podra haberse sacado para completar el
par. Habiendo sacado una cierta bola blanca de la caja A, entonces cualquiera de
las bolas en la caja B podra completar ese par, etc. El numero total de maneras
como podran obtenerse 2 bolas rojas es 2; la bola roja de la caja A slo puede
equiparse con cualquiera de las bolas rojas de la caja B. Luego, la probabilidad de
obtener simultneamente dos bolas rojas es: x 2/5 2/20, que es el producto de
las probabilidades individuales de sacar una bola roja de las cajas respectivas. La
probabilidad de sacar simultneamente dos bolas blancas es x 3/5 9/20, y la
probabilidad de obtener una roja y una blanca es entonces 1 (2/20 + 9/20)
9/20. Esta ultima probabilidad de 9/20 puede verificarse considerando que si se
saca la bola roja de la caja A, cualquiera de las tres bolas blancas de la caja B
Eventos independientes.
eventos, tales como:
P P1 * P2 *...........* Pn *.................(4c)
(1)
1
(2)
+1
-1
(3)
1
1
(4)
2
2
0
-2
3
1
-1
-3
1
2
1
1
3
3
1
4
2
0
-2
-4
5
3
1
-1
-3
-5
1
4
6
4
1
1
5
10
10
5
1
16
32
1/2
1/4
1/8
3/8
3/8
1/8
1/16
3/8
1/16
1/32
5/32
10/32
10/32
5/32
1/32
Los resultados de la tabla 4.1 pueden trazarse en una grafica de barras segn la
Fig. 4.1, de (a) a (e). El tamao del error resultante se traza como la abscisa
versus sus probabilidades en las ordenadas que se trazan como los centros de
barras del ancho unitario.
La ecuacin de la curva de distribucin normal del error es: (se da una deduccin
de esta ecuacin en el apndice A).
h e h x
2 2
valores medidos.
Se obtienen histogramas de formas varias para distintas personas y para las
variaciones en el equipo usado. Por ejemplo, el histograma de muchas mediciones
de distancias por geodmetro, probablemente ser una curva muy angosta,
inclinada en los lados con un alto valor en el centro. Un histograma de la misma
distancia medida por una barra subtensa y trazada a la misma escala
probablemente sera ms ancha, con los lados no muy inclinados y el valor central
no tan grande.
ngulo
ngulo
ngulo
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
532656.7
54.8
51.1
50.9
49.8
49.6
48.9
48.1
47.6
46.2
45.7
45.0
44.9
44.3
43.9
43.2
16.62
14.72
11.02
10.82
9.72
9.52
8.82
8.02
7.52
6.12
5.62
4.92
4.82
4.22
3.82
3.12
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
532642.8
42.4
42.1
41.6
41.1
40.8
40.3
40.1
39.8
39.8
39.3
39.3
39.0
38.8
38.6
38.5
2.72
2.32
2.02
1.52
1.02
0.72
0.22
0.02
-0.28
-0.28
-0.78
-0.78
-1.08
-1.28
-1.48
-1.58
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
532637.7
37.3
36.9
36.4
36.0
35.9
35.6
34.7
34.2
33.9
33.4
32.9
31.6
30.4
29.6
27.8
24.7
19.8
-2.38
-2.78
-3.18
-3.68
-4.08
-4.18
-4.48
-5.38
-5.88
-6.18
-6.68
-7.18
-8.48
-9.68
-10.48
-12.28
-15.38
-20.28
Tabla 4.2
De los datos de la tabla 4.2, el VMP (la media) es:
VMP
ngulos 53 26 40.08
"
n 1
7."377
Los residuales se calculan restando el VMP de cada valor medido. Los residuales
se tabulan en la columna v de la tabla 4.2.
Luego se clasifican los residuales en grupos o clases seleccionando algn
intervalo de clase compatible con los datos. Es importante seleccionar un
intervalo de clase que sea del tamao adecuado para exhibir correctamente la
distribucin. Para los datos de la tabla 4.2, se seleccionan 3 segundos como i- el
intervalo de clase-. El primero se centra en cero pero para que se cuente el
nmero de residuales en el intervalo de 3 segundos entre -15 y +15. Igualmente
se sigue el mismo procedimiento para los residuales de los intervalos entre -1.5 y
+4.5, +1.5 y +4.5,-4.5 y -7.5, +4.5 y +7.5, etc. Finalmente, el histograma con el
tamao de residual en la abscisa VS el nmero de residuales y las frecuencias de
ocurrencia como ordenada, se traza segn la figura 4.3.
h h 2 x 2 dx 1
h e
.. (4e)
Si diferenciamos la ecuacin (4d) obtenemos:
2 2
y
h
2h 2 x[( )e h x ]
x
h
2h 2 .
curva de distribucin normal ocurren en x =
h
0
Medidas de precisin
La precisin es la medida de la consistencia o grado de refinamiento de un grupo
de mediciones. Esto difiere de exactitud que se define como la proximidad al valor
verdadero de cualquiera medicin o grupo de mediciones. En agrimensura, los
trminos mas comnmente usados para expresar precisin son variacin y error
estndar ( desviacin estndar).
Variacin y error estndar
En el artculo 4.2 el trmino variacin se defini como el promedio de los
cuadrados de todos los errores. El error estndar se defini como la raz cuadrada
de la variacin. En cmputos de ajuste, es til relacionar la variacin y el error
estndar con la curva de distribucin normal.
Puede tambin definirse la variacin de una variable fortuita como su momento de
segunda rea, que puede expresarse como:
x da
da (4h)
Esta expresin produce la media de los cuadrados de las distancias de todas las
partes elementales de una figura desde una lnea central, que en el caso de la
curva de distribucin normal, djese que da = ydx, luego:
h h2 x2 2
2
x
ydx
e x dx
2
h2 x2 2
he
x dx
ydx
.. (4i)
h h2 x2 2
e
x dx
(4j)
De la ecuacin (4d) para y, y reduciendo:
d 2 y 4h 5 x 2 h2 x 2 2h 3 h 2 x 2
e
dx 2
d 2 y 4h5
dx 2
Pero
xe
2 h2 x2
2h 3
e h x dx
2 2
.. (4k)
d y
dy
dx 0
2
dx
dx
2
d2y
dx 2 dx
Y reordenando:
h
x 2 e h x dx
2 2
1
2h
e h x dx
2 2
. (4l)
h h 2 x 2 dx 1
h e
2 2
Y= 2
.. (4n)
Ntese que ya hemos establecido que el punto de inflexin de la curva de
1
2
distribucin normal ocurre en 2h ; por tanto, el valor de x en el punto de
inflexin es la raz cuadrada de la variacin o del error estndar.
Para cualquier grupo de mediciones normalmente distribuidas, el rea entre los
limites de a y a se da por la expresin
a
a
a
h h2 x2
a ydx 20 ydx 20 e dx
(4p)
Usemos la ecuacin anterior para determinar el rea dentro de los limites de .
1
dt
2h
x 0
t hx
et
dt
2
h
2h
1
2
e t dt
2
t 0
rea =
(4q)
la ecuacin (4q) igual a un rea de 0.50 y encontrar que valor del lmite superior
1
2 , luego t
x50
Por lo tanto para el error probable del 50%, simbolizado como E50 , 0.477= 2 y
x50 E50
=
= 2(0.477) 0.6745 (4s)
El error probable de 90 por ciento
Dicho error E90 , error tal que debiera caer tericamente dentro del 90% de los
errores de un grupo de observaciones, es bastante popular con los agrimensores
como un medio para expresar la precisin.
Extendiendo el razonamiento usado para el error probable del 50%, se escribe una
ecuacin semejante a (4r) para el error de 90%:
t
2
2
0.90
e t dt
0
0.90=
(4t)
Nuevamente de la tabla de Rainsford, el lmite superior de t que corresponde a un
rea de 0.90 es 1.163. Por tanto, el error probable del 90% es:
x90 E 2 (1.163) 1.6449
90
(4u)
Como previamente se indic, la desviacin estndar (S) y el error estndar ( )
son intercambiables en las anteriores ecuaciones.
212.23
15
VMP= n
Para calcular S, debemos calcular los residuales y sus cuadrados. Ello se hace
con facilidad en forma tabular como sigue:
Dist.
v
212.22 -.01 -.0001
212.25 .02 .0004
212.23 0
0
212.15 -.08 .0064
212.23 0
0
212.11 -.12 .0144
212.29 .06 .0036
212.34 .11 .0121
Dist.
22
24
19
25
27
20
25
3183.44
v
de de la
v
-.01 .0001
.01 .0001
.04 .0016
.02 .0004
.04 .0016
-.03 .0009
.02 .0004
.0421
Al extraer el valor
resuelve como sigue:
.0421
0.055
n 1
14
Pie
Por la ecuacin (4s)
E50 0.6745S .6745(0.055) 0.037
S
Pie
50(3)
7.377 2
v2
e 2S
v2
2
Ke 2 S
2
K
v2
2 s2
Donde:
50(3)
K
8.1
7.377 2
v
0
5
10
15
20
25
v2
0
25
100
225
400
625
v2
2S 2
0
.23
.92
2.07
3.68
5.75
v2
2S2
1.0
1.26
2.5
7.9
39.6
314
8.1
6.4
3.2
1.0
0.20
0.03
Los resultados de este cmputo se trazan de la Fig. 4.6. Ntese 2que adems de
trazar los valores y (ordenadas) vs. los valores v (abscisas), es til dibujar la curva
para ubicar tambin los puntos de inflexin que ocurren en una abscisa igual a la
desviacin estndar, o en el valor de 7.377.
value
1.677
1.676
n
16
17
value
1.663
1.665
n
31
32
value
1.664
1.690
n
46
47
value
1.664
1.685
3
4
5
6
7
8
9
10
11
12
13
14
15
1.657
1.667
1.673
1.671
1.673
1.670
1.675
1.664
1.664
1.668
1.671
1.664
1.651
18
19
20
21
22
23
24
25
26
27
28
29
30
1.670
1.671
1.651
1.665
1.677
1.662
1.660
1.667
1.660
1.660
1.667
1.667
1.652
33
34
35
36
37
38
39
40
41
42
43
44
45
1.649
1.671
1.675
1.653
1.654
1.665
1.683
1.668
1.658
1.657
1.690
1.666
1.671
48
49
50
51
52
53
54
55
56
57
58
59
60
1.667
1.655
1.679
1.682
1.663
1.662
1.672
1.667
1.667
1.663
1.670
1.667
1.669
Calclese lo siguiente:
El VMP del area medida.
La desviacin estandar
Trazar el histograma (de residuales) para los datos anteriores usando un intervalo
de clase de 0.003. Superponer el polgono de frecuencia en el histograma.
Calcular y trazar la CDN.
Calcular E50 .
Calcular E90.
CAPITULO 5
PROPAGACION DE ERRORES FORTUITOS
y1
zT ( y1ii x1ii ) ( y2ii x2ii ) y1ii y2ii x1ii x2ii z ii x1ii x2ii
zT ( y1iii x1iii ) ( y2iii x2iii ) y1iii y2iii x1iii x2iii z iii x1iii x2iii
, etc.
(5a)
zT z x x
Por nuestra definicin de variacin:
iii
iii
1
iii
2
2
2
n y n x
Luego para el caso bajo consideracin:
x x
2
i
1
i i
ii
Para una muestra grande, la suma de los productos 2 x1 x2 , 2 x1 , etc. Se convierte en
cero, puesto que los errores fortuitos pueden ser tanto positivos como negativos.
Por tanto,
z2
2
1
n
Finalmente:
2
2
y12 y22
z y2 y2
. (5b)
En general, si Z (un valor computado= es una funcin de y1, y2 yp. (Valores
medidos), entonces z , el error estndar de la funcin, se da por la siguiente
1
expresin:
dz
dz
dz
z y1
y2 ...
yp
dy1
dy2
dy p
PRUEBA
Supngase que Z= f(y1, y2, yp) en donde y1,y2, , yp son cantidades dadas
y i , y ii , y iii ,..., y1n
independientemente observadas. Tambin deje que 1 1 1
sean n
i
ii
iii
n
i
ii
x , x , x ,..., x1
y , y , y iii ,..., y2n
mediciones de y1 conteniendo errores 1 1 1
; permita que 2 2 2
x2i , x2ii , x2iii ,..., x2n
sean n mediciones de y2 conteniendo errores
y permita que
i
ii
iii
n
i
ii
y p , y p , y p ,..., y p
x , x , xiii ,..., x np
sean n mediciones de yp conteniendo errores p p p
.
Entonces el valor verdadero de Z, designado como zT , para cada dete4rminacion
independiente es:
Etc.
Tratando las y como aproximaciones iniciales y las x como correcciones, estas
ecuaciones pueden sustituirse en una serie de Taylor y linealizarlas desistiendo de
las potencias de orden dos o ms altas, como:
dz i
dz i dz i
zT z i
x p
x1
x2 ...
dy
dy
dy
1
2
p
dz ii
dz ii dz ii
zT z ii
x p
x1
x2 ...
dy
dy
dy
1
2
p
dz iii
dz iii dz iii
zT z iii
x p
x1
x2 ...
dy
dy
dy
1
2
p
Etc.
Etc.
Reconociendo las ecuaciones (5c) como equivalentes a (5a) de la prueba anterior,
entonces siguiendo la misma manipulacin general resulta:
2
dz
dz
dz
z
y1
y2 ...
yp
dy
dy1
dy2
p
(5d)*
2
funcin. Ntese que la ecuacin (5b) es implcita en la ecuacin (5d) como sigue:
Permtase que A = B + C
Luego por la ecuacin (5d):
dA
dA
1y
1
dB
dC
Sustituyendo:
1 B
1 C
B1
B2
B , B ,..., B
1
.. (5e)*
... Bn
, ,..., Bn
Si cada uno de los valores B1 B2
es igual a, digamos, B , entonces la
ecuacin (5e) se convierte en:
A n B .. (5f)*
*el error estndar y
ecuaciones.
z1 z2 ... zn
n
1
1
z1 z2 ... zn
n
n
n 2
2
n
n (5g)*
Ejemplos numricos:
Midamos las mediciones de un tanque (vase la Fig.), y obtenemos:
dv dv
dv
S v Sl
Sw Sh
dl dw dh
Sv
wh 0.05
2
lh 0.03 lw 0.02
2
cu.ft.
Ejemplo 2:
Supngase que el ngulo vertical de A al punto B se mide desde el punto A
como 300, S = +/- 1`; y que la distancia inclinada D, de A a B se mide como
100.00 pies, S D =+/- 0.05 ft (pie).
dH
dH
cos ;
D sin
dD
d
Sustituyendo en la ecuacin (5d):
H D cos ;
dH 2 dH 2
SH
SD
S
dD
d
D sin 1
S H cos 0.05
[57.296 60 ]2
2
pie
En el ejemplo anterior, ntese que la mayor fuente contributiva de error (el nmero
mas grande bajo el radical) es 0.0018, y que ese es el error asociado con la
medicin de la distancia. Por tanto, si el error resultante de +/- 0.045 pies es
demasiado grande para tolerarse, es adoptar un mtodo mas preciso a fin de
medir la distancia.
seccin
AB
BC
CD
DE
2
El volumen de un cono se da por V=1/12 D h. La altura medida del cono es de
10 pulgadas. Sh =+/- 0.20 pulg. El dimetro medido es 6 pulg. S D =0.10 pulg. La
desviacin estndar del volumen es:
CAPITULO 6
PESOS DE LAS OBSERVACIONES
Generalidades
El peso de una observacin es vala relativa de aquel valor observado segn se
compara a cualquier otro valor. En otras palabras, los pesos son clculos o
expresiones de la relativa confiabilidad de las observaciones.
Una alta precisin, indicada por un pequeo error estndar, implica una buena
observacin y, de esto, un peso alto. A la inversa, una baja precisin, indicada por
un error estndar grande, implica una observacin deficiente y, por tanto, un peso
pequeo. Es obvio que los pesos son inversamente proporcionales a alguna
potencia del error estndar. De ah que el tamao de la correccin tambin tendr
una relacin inversa con el peso.
Media ponderada (peso promedio)
Supngase que tomamos M observaciones independientes (z1, z2, z3,, zm) de
una cantidad z, y que cada una de estas observaciones tiene un error estndar s.
Luego la media de las observaciones tiene un error estndar s. Luego la media de
las observaciones es:
m
Z
z i
i 1 m . (6a)
Si estas m observaciones se separan ahora en dos juegos, uno de observaciones
ma y el otro con mb, de tal manera que ma +mb = m, entonces los promedios de
esos juegos son:
ma
Z
z a i
i 1 ma
z b
. (6b)
Zi
i m 1
b
. (6b)
m
zi
ma
zi
m
Pero de (6b) y (6c):
i 1
i m 1
i 1
i ma 1
ma mb
zi
ma
za ma zi
zb mb
i 1
y
Luego:
z m zb mb
z a a
ma mb
zi
i ma 1
za
y zb , podemos escribir la ecuacin general:
pa za pb zb
z
pa pb
pz
p (6d)
za
en la ecuacin (6b)
dza 2
dza 2 dza 2
...
dz1
dz2
dzma
. (6e)
Sustituyendo las derivadas parciales en (6e) tenemos:
za 2
1
1
1
1
1 2
z 2 ... 2 ma 2
m
m
m
m
m
a
a
a
a
a
Finalmente:
1
1
za 2 2
zb 2 2
ma Similarmente:
mb . (6f)
2
a
.
.
.
.
S n2 ln / d S 2
Por la ecuacin (6g), los pesos son inversamente proporcionales a las variaciones
y ya que son relativos, los pesos de la lnea de nivel son:
p1
1
d
1
2
2
S1 l1S
l1
p2
1
d
1
2
2
S 2 l2 S
l2
pn
1
d
1
2
2
S n ln S
ln
.
.
.
.
sn2 ( n / d ) s 2
Por la ecuacin (6g), los pesos son inversamente proporcionales a las variaciones
1
s1
d
1
2
1s
1
d
1
2
1s
2
d
1
2
ns
n
1
s2
pn
1
sn
Resumiendo, hemos probado que los pesos de las lneas de nivel son
inversamente proporcionales a sus longitudes, y puesto que la longitud es
proporcional al nmero de puestas en estacin los pesos tambin se hallan
inversamente proporcionales al nmero de puestas en estacin.
6.4.2 Pesos en la medicin de ngulos
Digamos que los tres ngulos 1, 2 y 3 en un tringulo plano se miden n1, n2 y
n3 veces con el mismo instrumento. Cules son los pesos relativos de los
ngulos?
Solucin:
Permita que S sea la desviacin estndar de una sola medicin de ngulo. Las
medias de cada uno de los tres ngulos son:
n1
; 2
n2
; 3
n3
1 2 2
1
1
) s ; s 2 ( ) s 2 ; s2 3 ( ) s 2
n1
n2
n3
Una vez ms, puesto que los pesos son inversamente proporcionales a las
variaciones y relativos, los pesos de los 3 ngulos son:
p1
p2
p3
n1
n1
s2
n2
n2
s2
n3
1
( s 1 )
1
( s 2 )
1
( s 3 )
s2
n3
BMX
Circuito de nivel
Solucin:
Del
Artculo 6.4.1, los pesos de las lneas de nivel son inversamente
proporcionales a sus longitudes. Por tanto, los pesos de las lneas 1, 2 y 3 son
1/2, 1/3 y 1/4, respectivamente. Puesto que los pesos son relativos,
arbitrariamente multiplicamos estas fracciones por 12 para obtener pesos de 6, 4 y
3, respectivamente.
(Ntese que si se hubieran desatendido los pesos, el promedio simple sera
+21.240).
6.6
6.6.1
Por definicin, se dice que una observacin tiene un peso p, donde p es un entero,
cuando su precisin iguala la precisin de la media de observaciones p de peso
unitario. Djese que 0 sea el error estndar de una observacin de peso unitario.
Si y1, y2, yn son observaciones que tienen errores estndar 1, 2, n,
y pesos p1, p2, ...... pn, entonces conforme a la definicin anterior, de la ecuacin
(5g) tenemos:
0
p1
; 2
0
p2
;......, n
(6h)
p n1
Ahora bien, en el caso cuando las observaciones no son de igual peso, debemos
modificar la ecuacin anterior para calcular el error estndar del peso unitario, a
saber:
p1 x1 p 2 x 2 ..... p n x n
n
2
px
n
(6i)
Y una vez modificada para la desviacin estndar del peso unitario tenemos:
p1 v1 p 2 v 2 ..... p n x n
n 1
2
s0
6.6.2
pv
(6j)
n 1
La relacin entre el error estndar del peso unitario y el error estndar del peso p
se dio en la ecuacin (6h). Combinando esa ecuacin con (6i), pueden obtenerse
ecuaciones para errores estndar del peso p en trminos de S0 como sigue:
px 2 1
px 2
0
1
; 1
*
p
n
np1
p1
1
0
p2
; 2
px
*
p
2
px
*
p
n
px
(6k)
np 2
0
p1
; n
px
n
np n
s1
pv
p1 ( n 1)
;s 2
pv
p 2 ( n 1)
....; s n
pv
p n ( n 1)
(6l)
px
pn
(6m)
sM
pv
p(n 1)
6.6.3
Ejemplos numricos
(6n)
Ejemplo 1:
Suponiendo que una distancia se midi como 625.79 pies por barra subtensa y se
le dio un peso de 1,
Se midi otra vez como 625.71 pies usando una
cinta y se le asign un peso de 2; la tercera medicin fue + 625.69 pies por
EDME (equipo electrnico para medir distancias) y recibi un peso de 4. Calcular
el VMP de la longitud (media ponderada) y determinar desviacin estndar de la
media ponderada.
a) Por la ecuacin (6d) la media ponderada es:
M
pv
p(n 1)
2
sM
Donde:
0.0080
0.024
7 ( 2)
Ejemplo 2:
En la nivelacin desde BMA hasta BMB, se tomaron cuatro diferentes rutas de
distintas longitudes. Se obtuvieron los datos siguientes:
Ruta
1
2
3
4
Longitud(mi.)
1
2
3
4
Elevacin
25,25
25,41
25,38
25,3
a. Hallar la media
(el VMP para la elevacin ):
Por la ecuacin (6d):
P
12
6
4
3
=25
Ruta
1
2
3
4
P
12
6
4
3
v
-0,065
0,095
0,065
-0,015
v2
0,00422
0,00902
0,00422
0,00022
pv2
0,0506
0,0541
0,0168
0,0007
0,1222
ponderada
0.1222
0.201
3
sM
0.1222
0.040
( 25)3
0.1222
0.058
12 (3)
s6
0.1222
0.082
6(3)
s4
0.1222
0.101
4(3)
s3
0.1222
0.117
3(3)
Valor
A
B
C
14 25' 20"
58 16' 00"
107 19' 10"
2
3
6
BM
Meas. Elev.
11
12
440.98 ft.
464.25 ft. La elevacin fija
482.23 ft, de BM 13 es 482.10
pies
13
CAPITULO 7
PRINCIPIOS DE LOS MNIMOS CUADRADOS
7.1 El principio fundamental de los mnimos cuadrados no ponderados.
En el Artculo 4.2.4, se declar la condicin fundamental de los mnimos
cuadrados; i. e., que el VMP era aquella cantidad, que haca de la suma de los
cuadrados de los residuales un mnimo. Vamos a probar esta condicin
considerando un caso especfico. Supngase que uno obtiene z1, z2 ., zn, como
n valores independientemente medidos de la misma cantidad. Vamos a designar
M como el VMP y luego por definicin:
z1 M V1
z 2 M V2
(7a)
z n M Vn
he h
2 2
Ke h
h2v 2
donde :
K
p1 y1v Ke h
p2 y2 v Ke h
v1
v
2
v2
v1 n 2
pn yn v Ke h
P ( Ke h v1 v)( Ke h v2 v)....( Ke h vn v)
2
Reduciendo
P K (v) n e h
( v1 2 v 2 2 ...v n 2 )
(7b)
(7c)
( z1 M ) 2 ( z2 M ) 2 ..... ( zn M ) 2
(7d)
d ( v 2 )
dM
0 2( z1 M )( 1) 2( z2 M )( 1) .... 2( z n M )( 1)
Reduciendo:
z1 M z2 M .... zn M 0
nM z1 z2 .... zn
(7e)
M ( z1 z2 .... zn ) / n
La ecuacin (7e) indica que el VMP para una cantidad que ha sido
independientemente observada varias veces, es sencillamente la media aritmtica.
En el proceso de esta prueba, se desarroll el principio fundamental del ajuste de
mnimos cuadrados: que el VMP hace de la suma de los cuadrados de los
residuales un mnimo.
p 2 v 2 2 .... p n v n 2 )
(7f)
Sustituyendo en la (7f), los valores para los residuales dados en (7a), tenemos:
p1 ( z M ) 2 p2 ( z2 M ) 2 .... pn ( zn M ) 2 Minimo
pz pM ;
Y de la cual:
pz
p
(7g)
Las ecuaciones que relacionan las cantidades medidas con sus errores
observacionales residuales y con parmetros independientes y desconocidos, se
llaman ecuaciones de observacin. Una tal ecuacin se escribe para cada
observacin. Para una solucin singular de las incgnitas, el nmero de
ecuaciones debe ser igual al nmero de incgnitas. En el caso normal hay ms
ecuaciones que incgnitas, lo que permite la determinacin de los VMP de las
incgnitas, basndose en los principios de mnimos cuadrados.
7.3.1 Ejemplo elemental del ajuste de la ecuacin de observacin.
Como ejemplo elemental del ajuste de mnimos cuadrados por el mtodo de la
ecuacin de observacin, considrense las tres siguientes ecuaciones de
observacin:
1) x y 3.0
2 ) 2 x y 1 .5
(7h)
3) x y 0.2
Estas tres ecuaciones relacionan las dos incgnitas x y y, con las cantidades
observadas. Los valores para x y y podran obtenerse de cualquiera de dos de
esas ecuaciones, para que la ecuacin restante fuera redundante. Ntese que
si las dos primeras ecuaciones de (7h) se resuelven, x=1.5 y y=1.5. Sin embargo,
si se resuelven las dos ltimas ecuaciones, x=1.3 y y=1.1 Adems, si se
resuelven la primera y ltima ecuaciones, x. = 1,6 y y = 1,4. Estas tres ecuaciones
han sido escritas como resultado de mediciones y se hace aparente que de esas
discrepancias observadas, las mediciones contienen errores. Tales ecuaciones
pueden, reescribirse como sigue:
4) x y 3.0 v1
5) 2 x y 1.5 v2
(7i)
6) x y 0.2 v3
Las ecuaciones (7i) relacionan las incgnitas con las ecuacones y con los errores
en las observaciones. Obviamente es posible seleccionar valores de v1, v2, y v3
que producirn los mismos valores y y, no importa el par de ecuaciones resuelto.
A fin de obtener estabilidades por todas las ecuaciones, podramos dejar que v1 =
0, v2 = 0, y v3= -0.2. Luego x = 1.5 y y= 1.5, no importa el par de ecuaciones
resuelto. Esto parece ser una solucin bastante buena.
Se pueden tener otros valores para las v, que harn suma de sus cuadrados un
valor ms pequeo. De hecho, la solucin que la suma de los cuadrados de
las v un mnimo, produce los valores x y probables y se llama la solucin de
mnimos cuadrados.
Para llegar a la solucin de mnimos cuadrados, se elevan cuadrado las
ecuaciones para los residuales y estas expresiones se suman juntas para dar una
funcin f(x,y) = Ev2, como sigue:
v1 ( x y 3.0) 2
2
v2 ( 2 x y 1.5) 2
2
(7j)
v3 ( x y 0.2) 2
2
( x y 3.0) 2 (2 x y 1.5) 2
(7k)
Luego se toman las derivadas parciales de esta funcin con respecto a cada
incgnita, y se establecen igual a cero. Entonces resultar una ecuacin para cada
incgnita, como:
v2
1)6 x 2 y 6.2 0
2 ) 2 x 3 y 1 .3 0
(7l)
Solucin aproximada
v1=0
v2=0
v3=0,2
v1 = -0,44
v2 = 0,086
v3 = -0,128
v1=0
v2=0
v3=0,4
v2 = 0,4
v1 = 0,002
v2 = 0,007
v3 = 0,016
v2 = 0,25
Una comparacin de la suma de los cuadrados de los residuales les indica que la
solucin de mnimos cuadrados produce el total ms pequeo; por ende, esto
resultara en un mejor ajuste; de hecho, es el mejor ajuste posible.
7.4 Formulacin sistemtica de las ecuaciones normales
En grandes sistemas de ecuaciones de observacin, es muy til emplear
procedimientos sistemticos para formular las ecuaciones nrmale! A fin de
desarrollar tales procedimientos, consideremos el sistema siguiente de ecuaciones
lineales de observacin de igual peso:
a1 A b1 B c1C ..... n1 N L1 v1
a2 A b2 B c2C ..... n2 N L2 v2
(7m)
am A bm B cmC ..... nm N Lm vm
v2 ( a2 A b2 B c2C ..... n2 N L2 )
2
(7n)
v1 ( am A bm B cmC ..... nm N Lm )
2
am A bm B cmC ...nm N Lm
(7p)
Las ecuaciones normales se forman como siempre tomando las
parciales:
derivadas
v2
v2
(7q)
Multiplicando a travs de cada una de las ecuaciones anteriores por 1/2, se
elimina de la constante 2. Las ecuaciones (7q ) son las normales. Vamos a
reordenar el orden de sus trminos como sigue:
a n a n ..... a n N a L a L ..... a L
0 b a b a ..... b a A b b ..... b B b c b c ..... b c C ...
b n b n ..... b n N b L b L ..... b L
0 c a c a ..... c a A c b c b ..... c b B c c ..... c C ....
c n c n ..... c n N c L c L ..... c L
0 a1 a 2 ..... a m A a1b1 a 2 b2 ..... a m bm B a1c1 a 2 c 2 ..... a m c m C ...
2
1 1
1 1
1 1
1 1
2 2
1 1
1 1
1 1
2 2
m m
m m
n
1
n 2 ..... n m N n1 L1 n 2 L2 ..... n m Lm
2
(7r)
Pueden ahora formularse ecuaciones generalizadas expresando las ecuaciones
normales (7r) como:
aa A ab B ac C ...........an N aL
ba A bb B bc C ...........bn N bL
ca A cb B cc C ...........cn N cL
(7s)
na A nb B nc C ...........nn N nL
En este sistema de ecuaciones (7s), las a, b, c, .... n, son los coeficientes de las
incgnitas A, B, C, .... N; L es el trmino constante y el smbolo ( significa la
suma de los productos).
Asimismo, puede probarse matemticamente que las ecuaciones normales se
pueden formar sistemticamente de ecuaciones de observacin ponderadas de la
manera siguiente:
paa A pab B pac C ........... pan N paL
pba A pbb B pbc C ........... pbn N pbL
pca A pcb B pcc C ........... pcn N pcL
(7t)
En las ecuaciones (7t) los trminos son los mismos segn descritos arriba para las
ecuaciones (7s), con la adicin de las p que son los pesos relativos de las
observaciones.
7.5
(7u)
x y 0.2 v3
Equn
1
2
3
a
1
2
1
b
1
-1
-1
L
3
1,5
0,2
aa
1
4
1
[aa]
6
ab
1
-2
-1
[ab]
-2
bb
1
1
1
[bb]
3
aL
3
3
0,2
[aL]
6,2
bL
3
-1,5
-0,2
[bL]
1,3
Sustituyendo los trminos apropiados para [aa], [ab], [bb], [aL] y [bL] de la Tabla
7.1 en las ecuaciones (7s), se obtienen las siguientes ecuaciones normales:
6 x 2 y 6 .2
2 x 3 y 1 .3
(7v)
Ntese que las ecuaciones (7v) son exactamente iguales a las (7l).
Se advierte basndose en estudios previos, que el nmero de ecuaciones
normales en un ajuste de mnimos cuadrados es igual al nmero de incgnitas que
han de ser ajustadas. Tratando con grandes sistemas de ajuste conteniendo
muchas incgnitas, resultara en un gran nmero de ecuaciones normales, por lo
que la solucin simultnea de las ecuaciones plantea un problema bastante
formidable. Tales ecuaciones pueden resolverse por mtodos matriciales segn el
Captulo 2, o por eliminacin gaussiana segn el Captulo 3.
Tarea autodidctica No. 7
1) Considerando que las siguientes observaciones son de igual peso, calcular
los VMP de A y B en las ecuaciones de observacin a continuacin, mtodo de
mnimos cuadrados no ponderados: (Usar el mtodo tabular para formar las
ecuaciones normales).
a ) A 2 B 16.0 v
b) 2 A 3B 12.5 v2
c ) 2 A B 0.5 v3
CAPITULO 8
AJUSTE DE REDES DE NIVEL
8.1 Introduccin
Ahora se considerar la aplicacin del mtodo de mnimos cuadrados por
ecuaciones de observacin al problema del ajuste del circuito de nivel. En el
Artculo 7.1 , se desarrollo la condicin de mnimos cuadrados para calcular el
VMP de un grupo de observaciones directas de la misma cantidad. En el ajuste de
circuitos de nivel , se trata con un grupo de observaciones directas que no es todo
de la misma cantidad, sino ms bien cada observacin representa la medicin de
la diferencia en elevacin de una lnea separada del circuito de nivel. No obstante,
tales observaciones se tratan de la misma manera y las correcciones buscadas
son aqullas que producen la ms alta probabilidad de su ocurrencia simultnea
en el circuito de nivel.
8.2 Ejemplo considerando pesos iguales
Considrese ahora un ejemplo especfico de un circuito de nivel, cuyo croquis y
datos del levantamiento se muestran en la Fig. 8.1. -En dicho ejemplo, se asumir
que cada observacin es de igual peso. Ms tarde se considerar el concepto de
las observaciones de peso relativo. En la figura las flechas indican la direccin de
la nivelacin. Por lo tanto para la lnea No. 1, la nivelacin procedi desde BMX
hasta A y la diferencia observada en elevacin fue de 5,10 pies.
Fig. 8.1
Lnea No
Diferencia
observada
en elevacin
1
2
3
4
5
6
7
5,1
2,34
-1,25
-6,13
-0,68
-3
1,7
(8a)
B A 0.68 v5
B BMY 3.00 v6
C B 1.70 v7
(8b)
v 42 (106 .13 C ) 2
v 52 ( B A 0.68) 2
v 62 ( B 104 .50 ) 2
v 72 (C B 1.70 ) 2
Al
resolver simultneamente las tres ecuaciones normales reducidas
anteriores, se producen las siguientes elevaciones ajustadas de las
Incgnitas:
A = 105.14
B = 104.48
C = 106.19
8.3
(8c)
Tabla 8. 1
Lnea
Longitud
Peso
relativo
1
2
4
3
3
4
4
5
6
7
3
2
2
2
4
6
6
6
p1 A p1 ( BMX 5.10 v1 )
p 2 BMY p 2 ( A 2.34 v 2 )
p 3 C p 3 ( BMY 14.25 v3 )
p 4 BMX p 4 (C 6.13 v 4 )
p 5 B p 5 ( A 0.68 v5 )
p 6 B p 6 ( BMY 3.00 v 6 )
p 7 B p 7 ( B 1.70 v 7 )
p 3 v3 6(C 106.25) 2
2
p 4 v 4 4(106.13 C ) 2
2
(8d)
p 5 v5 6( B A 0.68) 2
2
p 6 v 6 6( B 104.50) 2
2
p 7 v 7 6(C B 1.70) 2
2
(8e)
F
2 * 3( A 105.10) 2 * 4(105.16 A) 2 * 6( B A 0.68) 0
A
F
2 * 6( B A 0.68) 2 * 6( B 104.50) 2 * 6(C B 1.70) 0
B
F
2 * 6(C 106.25) 2 * 4(106.13 C ) 2 * 6(C B 1.70) 0
C
13A - 6B = 740.02
-A + 3B - C = 102.12
-3B + 8C = 536.11
= 105.10 + v1
= -105.16 + v2
C = 106.25 + v3
4)
-C
5) -A + B
6)
B
7)
-B + C
= -106.13 + v4
= -0.68 + v5
= 104.50 + v6
= 1.70 + v7
(8f)
Tarea autodidctica No 8
1)
En el croquis del circuito de nivel, calcular as elevaciones ms probables
para BM1 y BM2. Todas las lneas son de igual longitud y, por tanto, se aplica el
procedimiento de mnimos cuadrados de igual peso. Los valores observados se
dan en la tabla adjunta. Usar una solucin tabular para formular las ecuaciones
normales.
A=100.00
A=101.00
(2)
(1)
BM1
(3)
BM1
(4)
(5)
A=102.00
A=103.00
Lnea
Longitud
Difer.
elevacin
Lnea
Difer.
elevacin
1
2
3
4
5
1
0,08
0,92
-1
0,08
1
2
3
4
5
6
4
4
2
2
1
3
+1.05
-0.95
+2.10
-1.95
+0.10
+0.05
A=10
(5)
x3
(4)
(6)
x2
(1
(3)
A=101.00
CAPITULO 9
METODOS MATRICIALES PARA EL AJUSTE DE MINIMOS CUADRADOS
Am n X 1 n L1 n V1
(9a)
donde:
a1
a
2
A=
a m
L1
L
2
L3
L =
.
.
Lm
b1
c1
b2
c2
.
.
.
bm
cm
n1
n2
n m
A
B
C
X =
.
.
N
V1
V
2
V3
V =
.
.
V m
AT A X = AT L
(9b)
9.2
. .
'
Estas ecuaciones de observacin en forma matricial son:
7
Donde:
A3 3 X 1 7 L1 7 V1
0
1 0
1 0
0
0 0
1
A = 0 0
1
1 1 0
0
0 1
0 1 1
A
X = B ;
C
105 .10
105 .16
106 .25
L = 106 .13 ;
0.68
104 .50
1.70
v1
v
2
v3
V = v 4
v5
v 6
v
7
0 0
1
1 1 0 0 1 0 0
3 1 0
AT A = 0 0 0 0 1 1 1 * 0 0
1 = 1 3 1
1 1 0
0 0 1 1 0 0 1
0 1 3
0
0 1
0 1 1
Ntese que estos elementos matriciales son los coeficientes de las ecuaciones
normales. Vase el Artculo 8.2.
8 3 1
1
3 9 3
(AT A)-1 =
21
1 3 8
105 .10
105 .16
106 .25
1 1 0 0 1 0 0
210 .94
104 .50
1.70
Ntese que stos son los trminos constantes de las ecuaciones normales.
D.
X=
(AT
A) -1
AT
8 3 1
210 .94
105 .14
1
21
214 .08
1 3 8
106 .19
E.
0
1 0
105 .14
1 0
105 .14
0
0 0
106 .19
1
105 .14
A X = 0 0
1 * 104 .48 = 106 .19
1 1 0
0.66
106 .19
0
0 1
104 .48
0 1 1
1.71
105 .14
105 .10
0.04
105 .14
105 .16
0.02
106 .19
106 .25
0.06
104 .48
104 .50
0.02
1.71
1.70
0.01
(9b)
Pm ( a m A bm B .... n m N ) Pm ( Lm v m )
Usando los principios del lgebra matricial, resolvemos las X desconocidas corno
sigue:
(9g)
PAX = PL + PV
3
0
P = 0
0
0
0
0 0 0 0 0 0
4 0 0 0 0 0
0 6 0 0 0 0
0 0 4 0 0 0 ;
0 0 0 6 0 0
0 0 0 0 6 0
0 0 0 0 0 6
0
0
1
1 0
0
0
0
1
A = 0
0 1
1 1
0
1
0
0
0 1 1
v1
v
2
v3
V = v 4
v5
v 6
v
7
105 .10
105 .16
106 .25
L = 106 .13 ;
0.68
104 .50
1.70
3
0
0
1
1
0
0
1
0
0
1 1 1 * 0
ATP= 0 0 0 0
0
0 0 1 1 0 0 1
0
0
0 0 0 0 0 0
4 0 0 0 0 0
0 6 0 0 0 0
0 0 4 0 0 0 =
0 0 0 6 0 0
0 0 0 0 6 0
0 0 0 0 0 6
3 4 0 0 6 0 0
0 0 0 0
6 6 6
0 0 6 4 0 0 6
0
0
1
1 0
0
0
0
1
3 4 0 0 6 0 0
13 6 0
T
6 6 6 * 0
A PA = 0 0 0 0
0 1 = 6 18 6
1 1
0 0 6 4 0 0 6
0 6 16
0
1
0
0
0 1 1
Ntese que estos son los coeficientes de las incgnitas en las ecuaciones
normales ponderadas. Vase el articulo8.3.
(ATPA)-1
0.0933 0.0355
= 0.0355 0.0770
0.0133 0.0289
0.0133
0.0289
0.0733
105 .10
105 .16
106 .25
3 4 0 0 6 0 0
740 .02
T
6 6 6 * 106 .13 = 612 .72
A PL = 0 0 0 0
0.68
0 0 6 4 0 0 6
1072 .22
104 .50
1.70
Ntese que estos son los trminos constantes de las ecuaciones normales.
0.0933 0.0355
X = (ATPA)-1 ATPL = 0.0355 0.0770
0.0133 0.0289
0.0133
740 .02
105 .15
1072 .22
106 .20
0.0733
Resumiendo las elevaciones ajustadas para BMA, BMB y BMC son 105.15,
p 2 a 21 X 2 a 22 X 2 ..... a 2 n X n p 2 2 v 2
(9h)
p m a m1 X 1 a m 2 X 2 ..... a mn X n p m m v m
(9i)
donde:
p1
P=
p2
.
.
.
.
;A=
p m
1
2
.
L = . ;
.
.
m
a11
a
21
a
m1
a12
a 22
.
.
.
.
am2
a1n
a2n
;X=
a mn
X1
X
2
.
.
.
.
X
n
v1
v
2
.
V = .
.
.
v
m
(9j)
Pero de (8p)
V = AX L
(9k)
V ( AX L)
A
X
X
(9m)
(9n)
(9p)
Multiplicando (9p)
0 = ATPTAX - ATPTL
Pero p = pT puesto que es una matriz diagonal
Luego ATPAX = ATPL
(9q)
(9l)
X AT PA
AT PL (9r)
Si se supone que los pesos son iguales, la matriz P es una matriz identidad y la
solucin para X es:
X AT A
AT L
.(9s)
Ntese que las ecuaciones (9r) y (9s) son idnticas a las (9c) y (9g),
respectivamente.
9.6
Programa de Computadora
0.000000
0.000000
0.000000
0.000000
-1.000000
1.000000
1.000000
0.000000
0.000000
1.000000
1.000000
0.000000
0.000000
-1.000000
106.250000
106.130000
.680000
104.500000
-1.700000
MATRIZ DE PESO
1.000
1.000
1.000
1.000
1.000
MATRIZ DE ECUACIONES NORMALES
3.000000
-1.000000
0.000000
-1.000000
3.000000
-1.000000
0.000000
-1.000000
3.000000
MATRIZ DE COVARIACION
.380952
.142857
.047619
.142857
.428571
.142857
.047619
.142857
.380952
INCOGNITAS X
105.140900
104.482819
108.187592
1.000
1.000
RESIDUALES V
.040909
-.019089
-.062408
.057602
-.021919
-.017181
-.004773
DESVIACION ESTANDAR DEL PESO UNITARIO
.050123
DESVIACION ESTANDAR DE VALORES AJUSTADOS
.030938
-.032813
.030938
EJEMPLO DE RED DE NIVEL PONDERADA, COMPUTOS DE AJUSTE DE
WOLF
MATRIZ A
1.000000
-1.000000
0.000000
0.000000
-1.000000
0.000000
0.000000
0.000000
0.000000
0.000000
0.000000
1.000000
1.000000
-1.000000
0.000000
0.000000
1.000000
-1.000000
0.000000
0.000000
1.000000
MATRIZ L
105.100000
-105.160000
106.250000
-106.130000
-.680000
104.500000
1.700000
MATRIZ DE PESO
3.000
4.000
6.000
4.000
6.000
MATRIZ DE ECUACIONES NORMALES
13.000000
-6.000000
0.000000
-6.000000
18.000000
-6.000000
0.000000
-6.000000
18.000000
MATRIZ DE COVARIACION
.093333
.035556
.013333
.035556
.077037
.028889
.013333
.028889
.073333
INCOGNITAS X
105.150375
104.489182
106.197189
6.000
RESIDUALES V
.050385
.009628
-.052811
-.067184
.018806
-.010818
.008008
DESVIACION ESTANDAR DEL PESO UNITARIO
.107220
DESVIACION ESTANDAR DE VALORES AJUSTADOS
.032756
.029759
.029035
Tarea Autodidctica N 9
1. Resolver el problema 1 del juego de problemas N 8 usando mtodos
matriciales
2. Resolver problema 2 del juego de problemas N 8 usando mtodos
matriciales.
CAPITULO 10
PRECISIONES DE LAS CANTIDADES INDIRECTAMENTE DETERMINADAS
10.1 Desviacin estndar del peso unitario
En el artculo 6.6.1, la ecuacin (6j) expreso la desviacin estndar del peso
unitario como:
S0
pv
n 1
pv
mn
..(10a)
X AT PA
..(10b)
AT PL
Por lo tanto, los valores X son -al parecer- los errores en los valores ajustados de
mnimos cuadrados. Considrese ahora algn pequeo cambio incremental L en
los valores medidos L requeridos para cambiar X a X , de manera tal que la
ecuacin (10b) se convierte en:
X X A T PA
A T P L L
Reduciendo
X X A T PA
A
1
PL A T PA
A T PL
Pero de (10b),
X AT PA
X A T PA
A T P L
AT PL
as que
.(10c)
X A T PA
A T PV
PA
AT P B
Luego:
X BV
.(10d)
De la cual:
X X BVV T B T
T
.(10e)
X X
X 2
1
X 2 X 1
X 3 X 1
X X
n 1
X 1 X 2
X 2
X 1 X 3
X 2 X 3
X 3 X 2
X 3
X n X 2
X n X 3
X 1 X n
X 2 X n
X 3 X n
(10f)
X n
B
nm
v v
1 1
v 2 v1
v m v1
v1v 2
v2 v2
vm v2
..........v1v m
.........v 2 v m
.........v m v m
m m
BT
(10g)
mn
1
X X T
x1 2
,
a
x x
2 1 ,
a
x n x1
x x
1
x
a
x x
.....
1
2
2
.........
x
..........................
a
2
n
n n
(10h)
S 2
1
S 2 S1
S S
n 1
S1 S 2
S2
SnS2
........S1 S n
......S 2 S n
2
.......S n
n n
.(10i)
a
a
a
2
v v
v2
v 2 v m
2 1
....
a
a
a
B T (10j)
B
2
v m v1
v m v 2 ...... v m
a
a
a
En la matriz anterior, las v son residuales en cantidades medidas, y debido a que
estos residuales concuerdan con las distribuciones gausianas, probablemente
sern tanto mas como menos. Siendo ese el caso, para un gran valor de a , los
elementos no diagonales van a cero. Reconociendo los trminos diagonales como
variaciones de las cantidades observadas, reescribimos la matriz (10j) como sigue:
S 2
1
2
S2
2
S3
B
zeros
zeros
BT
2
Sn
.(10k)
S0
Pi
.(10l)
S0 B
P 1
1
1
P2
.................0
1
..............Pm
B T S 0 BP 1 B T
2
.(10m)
S 0 BP 1 B T S 0 AT PA
2
AT PP 1 P T A AT PA
1 T
A PA
T
1 T
A T PA
(10n)
S 0 AT PA
A
1
PA AT PA
S 0 AT PA
2
(10p)
S X S 0 AT PA
2
S0 Q
2
..(10q)
.(10r)
2
2 1
T
S 0 A A S 0 3 9 8
2l
1 9 8
Utilizando la ecuacin (10a), la estimada desviacin estndar del peso unitario es:
pv
S0
mn
pv
pv
Hance, S 0
0.0101
.0101
0.030 ft
73
Ahora bien, por la (10r) las estimadas desviaciones estndar de las elevaciones
desconocidas de las marcas de cota fija (BM), A, B y C son:
S A S 0 Q AA .05
8
.031 ft
21
S B S 0 Q BB .05
9
.033 ft
21
S C S 0 QCC .05
8
.031 ft
21
pv
mn
.0459
0.107 ft
73
2
2
S 0 A T PA S 0 .0332 .0770 .0289
..(10s)
SX
dx
dx
dx
Sl1
Sl 2 .............
Sl n
dl1
dl 2
dl n
.(10t)
Puesto que los coeficientes de las variaciones individuales en la (10t) son las
derivadas parciales con respecto a cada observacin, estos son sencillamente los
coeficientes lineales de las observaciones. Por tanto, la (10t) puede escribirse en
S X CS 2 C T ..(10u)
2
C C11C12 ............C1n ;
Donde:
Y:
SX
S 2
1
2
S2
zeros
CT
C11
C
12
C1 n
zeros
2
Sn
S AE
d AE
d AE
d AE
d AE
S d1
S d 2
S d 3
S d 4
dd 1
dd 2
dd 3
dd 4
S AE
S AE 0.13 ft
AAE CS 2 C T
2
2
S AE 1
S AE 0.0168;
2
.06 2
.08
from
which
.06 2
1
1
1
1
.08 2
S AE 0,13 ft
q
aa
q ba
q ca
q
na
q ab
q ac
q bb
q bc
q cb
q cc
q nb
q nc
q ab , etc.
....................n
................q an
................q bn
................q cn
..............q nn
ecuacin:
r a ,c
cov a , c
Sa Sc
.(10v)
ra ,c
S 0 q ac
S 0 q aa S 0
q cc
q ac
(10w)
q aa q cc
S a ra ,c S 0 q aa
Sc
S 0 q cc
q ac
q cc
..(10x)
q vx 0.105
0.41
q xx 0.258
bv , z
q vz 0.117
0.22
q xz 0.529
La Fig. 10.1 (b) bsicamente es la misma red de nivel que la de la Fig. 10.1 (a),
excepto que algunas de las lneas se han eliminado. La matriz Q obtenida como
resultado de este ajuste de mnimos cuadrados tambin esta ubicada a la
izquierda de la figura. Ntese por ejemplo que la incgnita V no esta directamente
conectada con ninguna otra incgnita, y por tanto no puede correlacionarse con
ninguna otra incgnita. Esto puede ilustrarse calculando su coeficiente de
regresin con respecto a las otras incgnitas, a saber:
bv , x
q vx 0.000
0.000,
q xx 0.318
bv , z
q vz 0.000
0.000
q zz 0.590
and
Otra vez la Fig. 10.1 (c) es la misma red bsica de nivel con aun mas lneas
eliminadas. La matriz Q obtenida de su ajuste de mnimos cuadrados se da a la
izquierda. Advirtase lo siguiente en esa figura:
1. La figura ya no es simtrica y hay mas lneas que van hacia la incgnita W
de las que van a la incgnita Y. Por eso W deber determinarse a una
mayor precisin que Y. Esto es verdadero ya que el elemento q ww es mas
pequeo (0.440) que el elemento q yy (0.593).
2. La incgnita Z esta directamente conectada con la incgnita X, y cualquier
cambio en elevacin que ocurriera a X afectara directamente la elevacin
de Z. En otras palabras, Z depende totalmente de X. Esto deber producir
un coeficiente de regresin de 1.0 y as sucede segn el calculo siguiente:
b z , x
V
.270
.123
Q
.105
.123
.117
.123
.105
.123
.366
.135
.116
.135
.258
.135
.116
.135
.366
.205
.176
.205
q zx 0.372
1.000
q xx 0.372
Z V
.117
W
.205
.176 X
.205
Y
.529
V
.500
.000
Q
.000
.000
.000
V
.322
.152
Q
.135
.067
.135
.000
.000
.000
.484
.181
.151
.181
.318
.181
.151
.181
.484
.227
.272
.152
.145
.067
.440
.169
.084
.169
.372
.186
.084
.186
.593
.169
.372
.186
.272
Fig. 10.1
Z V
.000
W
.272
.227 X
.272
Y
.590
.135
W
.169
.372 X
.186
Y
1.372
Z
Ejemplos de Correlaciones
saber:
X 1 l1 l 3
X 2 l1 l 2
.(10y)
X n a n1l1 a n 2 l 2 ............ a nm l m
S 2
l1
S l22
S
=
mm
S l2m
S X2 CS 2 C T
.(10z)
dX 1
dl
1
dX 2
dl1
Donde C =
dX n
dl
1
dX 1
dl 2
...........
dX 2
dl 2
..........
dX n
dl 2
..........
dX 1
dl m
dX 2
dl m
dX n
dl m
n m
f x
.......
1 n
X1
X 2
X n
n1
FX
(10aa)
f x A B
f x es:
S A B FCS 2 C T F T FS X F T
2
13
De (10q): S X S 0 Q
2
2 1
3 9 3
S0 2Q .050
21
1 3 8
Entonces
8 3 1
1
.05
3 9 3 1
1 0
21
1 3 8 0
S A B 2 1
S A B 2
.05
21
11
.0275
21
21
.036ft
X a X b 2 Ya Yb 2
.(10ab)
Tarea Autodidctica N 10
100.00
200.00
300.00
99.92
200.04
299.96
CAPITULO 11
AJUSTE DE LA TRILATERACION
11.1 Introduccin
La tcnica de la extensin del control de trilateracin se ha hecho sumamente
popular en la ultima dcada. La invencin y perfeccionamiento de los equipos
electrnicos para medir distancias (EDME), han sido responsables por esa enorme
popularidad. El EDME ha hecho posible la rpida y extremadamente precisa
medicin de lneas. La longitud o extensin de la lnea tiene poco efecto sobre el
tiempo y dificultad de una medicin, ya que lneas de hasta 50 millas o mas se
miden fcilmente con el equipo terrestre EDME que funciona transmitiendo
energa electromagntica. Los cursos sobre agua o terreno abrupto/irregular se
miden con comodidad. Es pues importante que se dedique bastante tiempo al
tema del ajuste de trilateracin de mnimos cuadrados en un texto moderno sobre
ajustes de levantamientos. En esta seccin en particular, vamos a considerar el
ajuste por el mtodo de la ecuacin de observacin, con cmputos efectuados en
un sistema de coordenadas rectangulares planas tal como los Sistemas de
Coordenadas Planas Estatales.
11.2 La ecuacin de observacin de distancia bsica
En el Articulo 3.2.4, se presento la siguiente ecuacin de observacin de la
distancia bsica:
Lij V Lij X j X i Y j Yi
Lau AU 0 V L
au
Lbu BU 0 V L
bu
Lcu CU 0 V L
cu
..(11a)
Xu Xa
Yu0 Ya
0
dX u
dYu
AU 0
AU 0
Xu Xb
Yu0 Yb
0
dX u
dYu
BU 0
BU 0
Xu Xc
Yu0 Yc
0
dX u
dYu
CU 0
CU 0
.(11c)
donde :
AU 0
BU 0
CU 0
X
X
X
u0
Xa
u0
Xb
u0
Xc
Y
Y
Y
2
u0
u0
u0
Y
Y
Ya
Lau ,
Lbu
X u0
Donde:
A es la matriz o coeficientes de las incgnitas.
X es la matriz o las correcciones desconocidas dX u
dYu .
X AT A
AT K
X b 2432 .55
Ya 4527 .15
Yb 2047 .25
X c 2865 .22
Yc 27 .15
Solucin
A. PRIMERA ITERACION
1. Calcular las aproximadas iniciales para X u0
Yu0
a. Calcular el acimut AB
X Xa
Az AB 180 tan 1 b
Yb Ya
2432 .55 865 .40
Az AB 180 tan 1
180 tan 1 .631940
b. Calcular longitud AB
AB 2 X b X a 2 Yb Ya 2
c. Calcular acimut AU 0
De la ley de cosenos: c 2 a 2 b 2 2ab cos C
cos UAB
UAB 50 6'50"
Az AU 0 147 42 '34"50 6'50" 97 35'44"
d. Calcular X u0
Yu0
2933.58 ft
2. Calcular AU 0 ,
AU 0
X u0
y
y
BU 0 ,
CU 0
Por tanto:
AU 0 6049 .00
BU 0 4736 .83
5446 .29 ft
Donde:
a12
a 22
a 32
4736.83
3727.59 27.15 0.680
5446.29
dYu
X AT A
AT K
0.991 0.933
T
A A
0.132 0.354
2.391
A A
0.699
T
A A
T
0.991
0.735
0.933
0.680
0.735
0.132
0.354
0.680
0.699
0.603
1 0.603
0.95 0.699
0.991 0.933
A K
0.132 0.354
T
1 0.603
X
0.95 0.699
0.699
2.391
0.00 0.136
0.680 0.20
SEGUNDA ITERACION
1. Calcular AU 0 ,
BU 0 ,
CU 0
BU 0
CU 0
6048 .968
4736 .915
5446 .451
3. La solucin matricial
.991 .933 .735 .032
.077 .011
A K
.005
1 .603
X
0.95 .699
N a de la estacin de control A
dY a en las ecuaciones de
dYb
dYc
dX c
dYd
dX d
AB
Yb Ya
AB
AC
Yc Ya
AC
Xc Xa
AC
AD
Yd Ya
AD
Xd Xa
AD
BC
BD
CD
Yb Yc
BC
Yb Yd
BD
Con excepcin de X a
Yc Yb
BC
Xc Xb
BC
Yc Yd
CD
Xc Xd
CD
Yd Yb
BD
Yd Yc
CD
Xd Xb
BD
Xd Xc
CD
AC AC 0
AD AD 0
K
BC BC 0
BD BD
0
CD CD 0
61
dYb
dYc
X dX c
dYd
dX
d
Vab
Vac
Vad
V
Vbc
V
bd
Vcd
51
61
IJ
dE j ,
dN j es:
.(11d)
.(11f)
..(11g)
E b E a tan N b0 N a
K Lab V Lab c
dN b (11h)
AB 0
dN b
dN c
dE c
dN d
dEd
Ec Ea
AC
Nd Na
AD
Ed Ea
AD
Nd Nb
BD
E d Eb
BD
Nd Nc
CD
Ed Ec
CD
E b E a tan N b N a
AB
AC
Nc Na
AC
AD
Nc Nb
BC
Ec Eb
BC
Nc Nd
CD
Ec Ed
CD
E b E d tan N b N d
BC
BC
E b E d tan N b N d
BD
BD
CD
AD ( AD )o
K
BC ( BC )o
BD ( BD )o
CD (CD )o
VLab
VLac
VLad
V
VLbc
VLbd
VLcd
dNb
dNc
X dEc
dNd
dEd
x
0
x
x
x
0
0
0
0
0
0
0
0
x
0
x
0
0
0
0
0
x
0
x
0
x
0
0
x
0
0
0
0
x
0
0
x
0
0
0
0
0
0
x
0
0
x
x
0
0
0
x
0
0
x
x
DE
EF
0
0
0
0
x
0
x
0
x
x
x
x
0
x
0
x
En este problema, los puntos Bucky y Badger son estaciones de control y cuyas
coordenadas se mantienen fijas; por tanto la direccin de la lnea Bucky-Badger
tambin se mantiene fija. Las 5 distancias observadas son:
LINEA
LONGITUD OBSERVADA (PIES)
Badger - Wisconsin
5.870,30
Badger - Campus
7.297,59
Wisconsin - Campus
3.616,43
Wisconsin - Bucky
5.742,88
Campus - Bucky
5.123,76
Las coordenadas de control de la estacin Badger son X = 10.000 y Y= 10.000, y
las de Bucky son X= 11.820 y Y = 6.881.222.
El Apndice C ofrece una detallada descripcin de la manera como se la da
entrada a los datos en la computadora para este problema. En las pginas
siguientes se lista la salida de computadora. En la pgina 168, se lista la
informacin de entrada. Esto es muy aconsejable en trabajos con computadoras
para asegurarse de recibir los correctos datos de entrada por la computadora. Las
pginas 169, 170 y 171 listan las matrices generadas por la computadora en cada
una de las tres iteraciones que se requirieron para lograr la convergencia. Ntese
en particular en esas pginas lo siguiente:
LONGITUDES Y ANGULOS
APROXIMACIONES INICIALES.
LONGITUD
PRELIMINARES
PARA
ANGULO
5870 80 0 0.000
3616 82 0 0.000
CALCULAR
ESTACION
BADGER
WIS
CAMPUS
BUCKY
Y
10000
15780.621
1698.231
11820.756
10000
11019.319
7580.3
6881.222
DATOS DE ENTRADA
LINEA
BADGER - WIS
BADGER - CAMPUS
WIS - CAMPUS
WIS - BUCKY
LONGITUD
5870.302
5871.302
5872.302
5873.302
LA MATRIZ A
0.98481
0.00000
0.30902
0.69139
0.00000
0.00000
0.17365
0.94363
0.00000
0.030902
0.95106
0.00000
0.72248
0.99065
0.00000
0.00000
0.33100
0.95106
0.00000
0.13640
LA MATRIZ L
0.30200
12.71572
0.43400
15.22937
1.61452
LA MATRIZ X
-0.09985
-4.96800
25.00019
24.26058
LA MATRIZ A
0.98405
0.00000
0.30874
0.68884
0.00000
0.00000
0.17790
0.94459
0.00000
0.030874
0.95115
0.00000
0.72491
0.98999
0.00000
0.00000
0.32825
0.95115
0.00000
0.14117
LA MATRIZ L
0.5635
0.2848
0.00129
0.03378
0.06457
LA MATRIZ X
-0.04630
-0.04970
-0.01519
-0.02065
LA MATRIZ A
0.98405
0.00000
0.30874
0.68884
0.00000
0.00000
0.17790
0.94459
0.00000
0.030874
0.95115
0.00000
0.72491
0.98999
0.00000
0.00000
0.32825
0.95115
0.00000
0.14117
LA MATRIZ L
0.0809
0.01169
0.00544
0.00912
0.00946
LA MATRIZ X
-0.00000
0.00000
0.00000
0.00000
LA MATRIZ Q
1.19891
0.09946
1.16024
1.40248
0.09946
1.16024
0.58294
0.19326
0.19326
2.63389
0.45964
2.72503
1.40248
0.45964
2.72503
3.96217
SX
10000
15776.675
16893.213
11820.756
Y
0
0.022
0.015
0
SY
10000
11044.304
7604.54
5881.222
0
0.022
0.015
0
LINEA
LONGITUD RESIDUAL (PIES)
BADGER - WIS
5870.31
0.008
BADGER - CAMPUS
7297.576
-0.012
WIS - CAMPUS
3616.439
0.005
WIS - BUCKY
5742.869
-0.009
CAMPUS - BUCKY
5123.769
0.009
PRELIMINARS
LONGITUD
16353
10993
34955
14469
22939
ANGULO
26
30
11
30
120
0
-170
-15
43
0
0
0
0
0
0
DATOS DE ENTRADA
Y
833766.69
841063.359
847831.323
860925.742
857855.7
869584.246
830040.4
886672.11
887981.23
454925.95
469560.813
478223.418
445813.724
459953.272
479667.199
467467.97
473334.4
453589.28
PARA
CALCULAR
LINEA
PAUL - BILL
RAM - TOM
PAUL - TOM
PAUL - JHON
PAUL - RAM
LUU - PAUL
TOM - JHON
BILL - JHON
TOM - TIM
PAT - TOM
TOM - BILL
BILL - TOM
LOU - BILL
BILL - PAT
TOM - JIM
PAUL- JIM
BILL - JIM
LONGITUD
30319.17
34954.73
31376.2
45986.22
10992.63
16353.41
39446.64
17445.22
29712.49
37437.93
35162.5
35162.49
43007.65
41305.83
14168.7
19896.23
22939.17
CODIGO
0
0
0
2
0
1
2
2
2
1
0
0
1
2
0
0
0
SX
833766.69
840831.926
848382.302
859667.382
857555.236
869923.518
830040.4
886672.11
Y
0
0.175
0.403
0.171
0.312
0.155
0
0
SY
45925.95
469674.295
477663.646
444580.718
458894.734
478214.199
467467.97
473334.4
0
0.231
0.282
0.198
0.261
0.248
0
0
LINEA
PAUL - BILL
RAM - TOM
PAUL - TOM
PAUL - JHON
PAUL - RAM
LUU - PAUL
TOM - JHON
BILL - JHON
TOM - TIM
PAT - TOM
TOM - BILL
BILL - TOM
LOU - BILL
BILL - PAT
TOM - JIM
PAUL- JIM
BILL - JIM
PUNTO
A
B
C
Y
0
1000
0
0
0
1000
CAPITULO 12
AJUSTE DE LA TRIANGULACION
12.1 Introduccin
La triangulacin es un mtodo de extensin del control horizontal en el cual las
mediciones de ngulos son las observaciones bsicas. Luego se computan las
posiciones de los puntos extensivamente espaciados basados en los ngulos
medidos, y slo a un nmero mnimo de distancias medidas se les llaman Lneas
Base. Se toman todas las precauciones posibles en las mediciones de las lneas
base para asegurarse de su ptima precisin. Antes de advenimiento del EDME,
con frecuencia la triangulacin era el mtodo ms conveniente, econmico y
preferido para extender el control horizontal bsico.
Fue el mtodo usado por el Servicio Geodsico y de Costas de EE.UU. para
extender mucha de la red nacional. Hoy da, la triangulacin contina siendo
popular con muchos agrimensores.
El ajuste de triangulacin por mnimos cuadrados puede efectuarse por
ecuaciones de condicin, ecuaciones de observacin de direccin o ecuaciones
de observacin de ngulos. En esta seccin, consideraremos la triangulacin por
el mtodo de la ecuacin de observacin de ngulos mediante mnimos
cuadrados. Los cmputos se harn en un sistema de coordenadas rectangulares
planas tal como los Sistemas de Coordenadas Planas Estatales. Los tipos
especficos de triangulacin que se tratarn son interseccin, reseccin y la
tradicional triangulacin cuadrilateral.
12.1 La ecuacin bsica de observacin de ngulos.
Xk Xi
Xj Xi
Kjik Vjik Yj
dXi
dYi
2
2
2
( Iko )
( Iko ) 2
( IJo )
( IJo )
o
Yio Yjo
Xj o Xi o
YK o Yio
Xi o Xk o
dXj
dYj
dXk
dYk .......(12b)
2
2
2
2
( IJo )
( IJo )
( IK o )
( IK o )
En las ecuaciones (12) y 12b), los smbolos se definen como en el Artculo 3.2.5.
La ecuacin (12b) se deriv para un ngulo en el primer cuadrante segn se
muestra en la Fig. 3.2. En la fig. 12.1 se muestran las doce posibles
combinaciones de cuadrantes para ngulos inferiores a 180.
La ecuacin (12b) se aplica exactamente a los ngulos de la Fig. 12.1 de (g) a (1),
excepto que Kjik se calcula como sigue:
V
41
3=984117
4=591701
Xb= 2432.55
Xc= 2865.22
Ya= 4527.15
Yb= 2047.25
Yc= 27.15
Solucin
1.
ABSeno 2
2933.58Seno10130'47"
6049 ft.
Seno(180 1 2)
Seno28 27'23"
Xb Xa
1 2432.55 865.40
AzimutAB tan 1
tan
147 42'34"
Yb Ya
2047.25 4527.15
AzimutAUo 147 42'34"50 06'50" 97 35'44"
Xuo Xa AUoSeno( AzimutAUo ) 865.40 6049.00 Seno97 35'44" 6861.35
Yuo Ya AUoCos ( AzimutAUo ) 4527.15 6049.00Cos97 35'44" 3727.59
BUo (6861.35 2432.55) 2 (3727.59 27.15) 2 4736.83 ft.
CUo (6861.35 2865.22) 2 (3727.59 27.15) 2 5446.29 ft.
2.
Yuo Yb
Xb Xuo
1 Yuo Yb
1 Ya Yb
(UBo ) 2 (dXu ) (UBo ) 2 (dYu ) 2 180 tan Xuo Xb tan Xa Xb v 2
Yb Yuo
Xuo Xb
1 Yuo Yb
1 Yc Yb
(BUo ) 2 (dXu ) (BUo ) 2 (dYu ) 3 tan Xuo Xb tan Xc Xb v 3
Yuo Yc
Xc Xuo
1 Yuo Yc
1 Yb Yc
(UCo ) 2 (dXu ) (UCo ) 2 (dYu ) 4 180 tan Xuo Xc tan Xb Xc v 4
(6049 .00 ) 2
A p
( 2047 .25 3727 .59 )
( 4736 .83) 2
(5446 .29 ) 2
(6049 .00 ) 2
(6049 .00 ) 2
15 .447
(5446 .29 ) 2
33 .800
40 .713
40 .713
27 .788
50 05'50"
101 30'47"
K
98 41'17"
5917'01"
3727.59 4527.15
1 2047.25 4527.15
tan 1
tan
6861.35 865.40
2432.55 865.40
0.00"
4527.15 2047.25
1 3727.59 2047.25
180 tan 1
tan
865.40 2432.55
2865.22 2432.55 0.00"
0.69"
3727.59 2047.25
1 27.15 2047.25
tan 1
tan
6861.35 2432.55
2865.22 2432.55 20.23"
2047.25 27.15
1 3727.59 27.15
180 tan 1
ta
n
2432.55 2865.22
6861.35 2865.22
1820 .5
5229 .7
0.00190
X ( A T A) 1 ( A T K )
0.00066
0.00190
( A T A) 1
0.00066
0.00066
0.00042
0.00066 509 .9
dXu
.
dYu
0.00042
534
.
1
donde :
25 .732
33 .800
0.00"
6.5"
5.1"
40 .713 0.62
0.00"
.
5.8"
40 .713
0.11 0.069"
27 .788
20 .23"
7.3"
V T V 6.5
6 .5
5 .1
155 .2 Sec . 2
7.3.
5 . 8
7 .3
So
V TV
mn
5 .1
5 .8
155 .2
8.8"
42
509 .9
AT K
534 .1
dXu 0.62 ft
dYu 0.11 ft
observaciones mltiples.
12.7 Ejemplo del problema
Para la figura a continuacin, se dispone de las siguientes observaciones y datos
conocidos:
ngulos observados (suponer pesos iguales)
Coordenadas fijas:
Xa= 4270.33
Ya= 8448.90
Xd= 7610.58
Yd= 4568.75
Se us el programa de computadora del Apndice D para resolver el problema
anterior. Se emplearon pesos iguales. En el programa, los ngulos son entrada
en orden del 1 al 8. La matriz X en la solucin es:
dYb
dXb
X
dYc
dXc
4270 .330
4568 .750
7610 .580
(Ya )
( Xa )
(Yd )
( Xd )
5599 .549
(Yb )
( Xb )
( Xc )
0
0
0
0
0 0 0 0 0 0 0
1 0 0 0 0 0 0
0 1 0 0 0 0 0
0 0 1 0 0 0 0
0 0 0 1 0 0 0
0 0 0 0 1 0 0
0 0 0 0 0 1 0
0 0 0 0 0 0 1
20.108 16.201
7.745
6.602
0.000
0.000
0.285
22.830
2.722
16.485
2.722 16.485
12.255
12.255
22.830
0.000
4.824
10.575
0.000
7.430
9.682
9.682
0.285
0.000
3.086
9.967
0.000
12.769
MATRIZ L (VECTOR)
4.650
4675
4.262
4.291
7.805
4.230
7.904
8.593
0.007
0.008
INCOGNITAS VECTOR X
0.646
VECTOR DE RESIDUOS V
2.152
ITERACIN N 2
MATRIZ A
5.010
4.162
1.534
24.229
3.879
0.000
0.000
20.111 16.200
7.744
6.601
0.000
0.000
0.284
22.832
2.722
16.484
2.722 16.484
12.254
12.254
22.832
0.000
4.824
10.579
0.000
7.430
9.681
9.681
0.284
0.000
3.086
9.966
0.000
12.767
MATRIZ L (VECTOR)
2.152
6.509 1.534
0.007
0.008
INCOGNITAS VECTOR X
0.000
VECTOR DE RESIDUOS V
2.152
5.010
4.162
1.534
5599 .356
(Yb )
( Xb )
( Xc )
SYB= 0.500
SXB= 0.227
SYC= 0.048
SXC= 0.580
ANGULOS AJUSTADOS
4235'31.15"
8735'15.61"
7954'37.94"
1828'20.92"
2129'25.70"
391'29.98"
3120'52.31"
3934'26.37"
Las
Y
10,000
10,000
10,000
10,000
10,000
7,500
2,500
0
CAPITULO 13
AJUSTE DE POLIGONALES
13.1 Introduccin
De los muchos mtodos para el ajuste de poligonales, la caracterstica que
distingue el ajuste de poligonales por mnimos cuadrados de los mtodos
aproximados, es que las observaciones de distancia y direccin se ajustan
simultneamente en el mtodo de mnimos cuadrados que con cualquier otro
mtodo.
Adems, el mtodo de mnimos cuadrados permite la asignacin de pesos
relativos a las observaciones basndose en sus confiabilidades relativas
esperadas.
13.2 Las ecuaciones de observacin
En este captulo se presenta el ajuste de poligonales por mnimos cuadrados
escribiendo ecuaciones de observacin. En este procedimiento, se escribe una
ecuacin de observacin para cada distancia y ngulo medido ( o direccin). La
ecuacin de observacin bsica para la distancia, una ecuacin no lineal, se
formul y lineariz en el cpitulo 11 como la ecuacin (11b).
En esta presentacin se usarn ecuaciones de observacin de ngulos, aunque
como una alternativa, es posible formular ecuaciones de observacin de direccin.
Se ha seleccionado el concepto angular porque es algo mas simple de usar, y
tambin porque el nmero de incgnitas, en la solucin es menor. La ecuacin
bsica de observacin de ngulos se formul en el Cpitulo 3 como la ecuacin
(3m). Se repiti en el Captulo 12 como la (12b). En el Captulo 12 se demostr
su uso en el ajuste de la triangulacin.
El hecho que las observaciones de ngulos y distancias con unidades diferentes
se combinan en un ajuste, se resuelve mediante el uso de pesos apropiados.
En este procedimiento se asigna un peso relativo a cada observacin de acuerdo
con la inversa de su variacin y priori (el cuadrado de la desviacin estndar).
Como antes se indic, la condicin fundamental que se hace cumplir en el ajuste
por mnimos cuadrados es la minimizacin de los trminos pv. De ah que si se
ponderan las observaciones conforme a la inversa de sus variaciones a priori, y si
se usan las mismas unidades para la desviacin residual y estndar, todos los
trminos pv quedarn sin dimensin alguna y, por tanto, compatibles para el
ajuste simultneo. Este procedimiento se demuestra en el problema de la Seccin
13.4.
Observaciones
AB = 200.00
BC = 100.00
1 240 00 '
2 150 00 '
3 240 01'
Fig. 13.2.
A. Calcular las coordenadas aproximadas iniciales para la estacin B
X bo 100 200 Sen60 1173.20
ybo 100 200Cos 60 1100.00
d xb
X
d yb
K Lab
K Lbc
K K1
K
2
K3
200.00 200.00 ft
0.00
100.00 99.81 ft
0.19
* Ntese
1100 1000
1173.20 1000
200
200
1100 1186.50
1173.20 1223.00
100
100
1173.20 1000
1000 1100
2
2
200
200
2
2
2
2
200
100
200
100
1223.00 1173.20
1100 1186.50
2
2
100
100
0.866
0.498
A 515.7
2299.8
1784.2
0.500
0.865
893.1
1920.3
1027.2
D. Formular la matriz P.
Para resolver unidades dismiles, se usan pesos de acuerdo con las inversas de
las variaciones. Puesto que antes del ajuste se desconocen las variaciones, stas
tienen que estimarse a priori; y son importantes porque influyen en el ajuste. Los
PLij
Para longitud:
P jik
Para ngulo:
S Lij
jik
AB 0.05 ft
BC 0.08 ft
1 30 sec.
2 30 sec.
3 30 sec.
Basndose en esas desviaciones estndar, la matriz P es:
1
0.052
1
0.082
1
30
1
30
1
30
400
156
0.001
0.001
0.001
X AT PA AT PL
1
dX b 0.11 ft
dYb 0.01 ft
Los residuales sin como sigue:
Vab 0.10 ft
Vbc 0.12 ft
V 1 49.0 sec
V 2 17.4 sec
V 3 604 sec
dX b
dYb ).
MATRIZ A
0.866000
0.498000
515.700000
2299.800000
1784.200000
MATRIZ L
0.500000
0.865000
893.100000
1920.300000
1027.200000
0.000000
0.190000
3.000000
249.000000
192.000000
MATRIZ DE PESO
400.000
MATRIZ DE COVARIACIN
0.000553 0.000622
0.000622 0.000872
INCGNITAS X
-0.111379
-0.012762
RESIDUALES Y
-0.102835
-0.123494
49.040077
-17.358521
-6.387482
DESVIACIN ESTNDAR DEL PESO UNITARIO
1.765992
ngulo
Valor
observado
Desviacin
estndar(seg)
Longitud
Valor
observado
XAB
633430"
AB
2731.25
Desviacin
estndar
(pies)
0.29
4512'15"
10
BC
3005.81
0.32
ABC
943920"
10
CA
4222.01
0.44
BCA
400905"
10
BAC
Acimut AX = 345173
Coordenadas de A:
X = 10,000.00
Y = 10,000.00
CAPITULO 14
AJUSTE POR CONDICIONES
14.1 Introduccin
Hay una clase de problemas en los cmputos de ajustes en la cual ciertas
condiciones tienen que satisfacerse exactamente en el ajusta de las
observaciones. Las ecuaciones que expresan esas condiciones se llaman
ecuaciones de condicin. En el ajuste de mnimos cuadrados por el mtodo de
ecuacin de condicin, se halla la solucin que produce las cantidades ajustadas
ms probables y al mismo tiempo satisface exactamente todas las condiciones.
Claro est que la solucin que da las cantidades ajustadas ms probables, es la
que minimiza la suma de los pesos multiplicada por los residuales elevados al
cuadrado.
El clsico ejemplo de este tipo de condicin es el ajuste de los tres ngulos
observados en un tringulo plano. La condicin que se debe satisfacer es que la
suma de los tres ngulos ajustados debe ser exactamente igual a 180 grados. El
ajuste por mnimos cuadrados es el que produce simultneamente los valores
ajustados que tienen la ms alta probabilidad.
En redes de triangulacin, hay muchas condiciones geomtricas que deben
cumplirse y, por ende, frecuentemente se ajustan por el mtodo de condiciones.
Tambin pueden ajustarse las poligonales por dicho mtodo. Por ejemplo, el
mtodo Crandall es un ajuste donde despus del ajuste del cierre angular de
manera tal que cada ngulo recibe una correccin igual, todo
el error restante se supone que este contenido en las observaciones de distancia.
Se escriben ecuaciones que expresan las dos condiciones de cierre de latitud y de
partida.
Una de las dificultades confrontadas en el ajuste por condiciones es la de
reconocer y formular todas las ecuaciones de condicin. Debe recalcarse que
todas las condiciones deben considerarse para obtener el ajuste correcto. Para
ciertos problemas de ajuste, se puede reducir la cantidad de clculos si el ajuste
se efecta por ecuaciones de condicin en lugar de observaciones. Sin embargo,
para otros problemas, lo opuesto es cierto. En esta era de computadoras
electrnicas, la cantidad de clculos parece proporcionar muy poca Justificacin
para el ajuste mediante condiciones. En pro del mtodo de la ecuacin de
observacin se encuentran la fcil, sistemtica formulacin de las ecuaciones, y el
conveniente y sistemtico anlisis del error proporcionado por la matriz de
covariacin.
b11
v1 b12
v2
b1m
b21
v1 b22
v2
b2 m
m
m
vm c1
v m c2
. . . .(14a)
bq1
v1 bq 2
v2
bqm
vm cq
y sus
. . . . . . . . . . . . . . . . (14c)
Luego BV = W
La expresin matricial para la solucin de los residuales que tienen la probabilidad
ms alta (solucin por manimos cuadrados) se da por:
Vm1 BT mq Bqm BT mq Wq1 * . . . . . . . . . . . (14d)
1
. . . . . . . . . (14g)
Estacin
A
B
C
ngulo
observado
681040
6634 22"
451449
BB
T
T
31
V B
13
T
31
BB
T
1 / 3
1
1
1/ 3
W 1/ 3
1/ 3
1/ 3
1/ 3 1/ 3
1/ 3
3 ''
9 '' 3 ''
3 ''
respectivas observaciones:
A ' A V A 6810 ' 40 '' 03 '' 6810 ' 43 ''
B ' B VB 66 34 ' 22 '' 03 '' 66 34 ' 25 ''
C ' C VC 4514 ' 49 '' 03 '' 4514 '52 ''
-------------------180 00 00
Lnea No.
Dif. elev.
observada
1
2
3
4
5
6
7
5.10
2.34
-1.25
-6.13
-0.68
-3.00
1.70
Fig 14.2
En el ajuste de redes de nivel, el nmero de condiciones es igual al numero de
observaciones menos el numero de Incgnitas. Por tanto, en este problema hay 73, 4 ecuaciones de condicin que deben formularse. Tales condiciones son:
V1
V3
V4 7.50
V2
V5
V2 7.50
V3
V6 0
V6
V7 0
-1
-1
0
0
7.50 1
7.50 0
W
0.00 0
0.00 0
-1
-1
0.06
0.12
W
0.02
0.05
5.10
2.34
1.25
6.13
0.68
3.00
1.70
BB T
BB
T
BB
T
2
0
0.619
0.048
0.238
0.095
0.619
0.381
0.048
0.048
0.238
0.143
0.095
0
1
0.048
0.238
0.619
0.095
0.095
0.476
0.238
0.190
0.048
0.238
0.047
0.238
0.381
0.095
0.619
0.095
0.095
0.476
0.143
0.286
0.238
0.190
0.095
0.228
0.190
0.476
0.095
0.095
0.238
0.238
0.190
0.286
0.467
0.041
0.019
0.06
0.062
0.12
0.058
V B T BB T W B T BB T
0.02
0.022
0.05
0.017
0.005
0.0101
0.050 ft
7-3
m-n
P 0
0
0
0
0
0
0
0
0
6
1.0000000
0.0000000
0.0000000
0.0000000
0.0000000
0.0000000
1.0000000
0.0000000
1.0000000
-1.0000000
0.0000000
1.0000000
0.0000000
0.0000000
-1.0000000
1.0000000
-0.1200000
-0.0200000
1.0000000
1.0000000
1.0000000
1.0000000
MATRIZ W
0.0600000
0.0500000
MATRIZ P
1.0000000
1.0000000
1.0000000
RESIDUALES
0.0409524
0.0190476
-0.0623810
-0.0576190
0.0219048
-0.0219048
0.0047619
EJEMPLO DE RED DE NIVEL PONDERADA, CMPUTOS DE AJUSTE DE
WOLF
MATRIZ B
1.0000000
0.0000000
0.0000000
0.0000000
1.0000000
0.0000000
0.0000000
0.0000000
0.0000000
0.0000000
0.0000000
1.0000000
1.0000000
0.0000000
0.0000000
0.0000000
0.0000000
0.0000000
1.0000000
0.0000000
1.0000000
-1.0000000
0.0000000
1.0000000
0.0000000
0.0000000
-1.0000000
1.0000000
-0.1200000
-0.0200000
4.0000000
6.0000000
6.0000000
6.0000000
MATRIZ W
0.0600000
0.0500000
MATRIZ P
3.0000000
4.0000000
6.0000000
RESIDUALES
0.0504000
0.0096000
-0.0528000
-0.0672000
0.0188000
-0.0108000
0.0080000
14.6 Ajuste de poligonales por el Mtodo de ecuaciones de condicin
(mtodo Crandall)
Algunas veces se le llama errneamente al mtodo Crandall como el mtodo de
mnimos cuadrados. En realidad, una vez que se distribuye el error angular de
cierre en porciones iguales a todos los ngulos, el mtodo Crandall mantiene fijos
esos ngulos y coloca a todas las correcciones restantes en las mediciones
lineales mediante un procedimiento de mnimos cuadrados ponderados. Sin
embargo, no es un verdadero ajuste por mnimos cuadrados, ya que las
mediciones angulares y lineales no se ajustan simultneamente.
El mtodo CrandalI es vlido para ajustar poligonales en donde se espera que las
mediciones lineales contengan errores fortuitos ms grandes que las mediciones
angulares. Como ejemplo de una situacin ideal para el mtodo Crandall es el
ajuste de una poligonal de estada.
La poligonal de la Fig. 14.3 se ajustara usando el mtodo Crandall. Dicha poligonal
contiene adrede un cierre absurdamente grande, para que el procedimiento sea
mas claro. El proceso computacional del mtodo Crandall se compone de la
formulacin de las dos ecuaciones que expresen las condiciones de que, las
sumas algebraicas de las latitudes y partidas en una poligonal cerrada debern ser
igual a cero.
Primero se calculan las marcaciones ajustadas, seguidas del clculo de latitudes y
partidas, segn la tabla:
Curso
AB
BC
CD
DE
EF
FA
Extensin
748.25
532.18
659.74
440.33
345.90
445.02
Marcacin
ajustada
Due North
N 44 48' E
S 40 14' E
S 24 14' W
S 89 52' W
S 49 47' W
Lat.
748.25
448.58
-503.66
-401.53
- 0.81
-287.34
N= 3.49
Part.
0.00
445.46
426.13
-180.73
-345.90
-339.82
E= 5.14
Las dos condiciones que han de cumplirse, exclusivo de los pesos, se expresan
en las siguientes ecuaciones:
Sumas algebraicas de latitudes = 0
Sumas algebraicas de partidas = 0. . . . . . . . . . . . . (14g)
AB Cos (Bng AB) + BC Cos (Bng BC) + CD Cos (Bng CD) + DE Cos (Bng DE) +
EF Cos (Bng EF) + FA Cos (Bng FA) = 0
. . . . . . . . . . (14h)
AB Sin (Bng AB) + BC Sin (Bng BC) + CD Sin (Bng CD) + DE Sin (Bng DE)
EF Sin (Bng EF) + FA Sin (Bng FA) = 0
ngulos interiores
Punto Observado Ajustado
A
4945
4947
B
13510
13512
C
8500
8502
D
11530
11532
E
11420
11422
F
22003
22005
Fig. 14.3
. . (14i)
Donde
1.000
B
0.000
7.10
0.763
0.912
0.002
0.705
0.646
0.410
1.000
0.646
0.764
W1 0 3.49 3.49
W2 0 5.14 5.14
En el mtodo Crandall se supone que los pesos sean inversamente proporcionales
a las longitudes. Por ende, la matriz inversa de peso para esta poligonal es:
748.25
0
0
0
0
0
632.18
639.74
440.33
345.90
445.02
MATRIZ B
1.0000000
-0.9120000
0.0000000
-0.4100000
0.7100000
-0.0020000
0.7050000
-1.0000000
-0.7630000
-0.6460000
0.6460000
-0.7640000
MATRIZ W
-3.4900000
-5.1400000
MATRIZ P
1.3400000
2.2700000
RESIDUALES
-0.7615914
-2.1829762
-1.0643644
1.0800547
1.2843576
1.5526733
1.5300000
2.8900000
1.5200000
2.2500000
Fig. 14.4
Inherente en estas tres ecuaciones de ngulos es la ecuacin para el cuarto
AD
Sin 2
Sin 5
and BC
CD
Sin 8
Sin 2
Sustituyendo;
BC = AD Sin 2 Sin 8
Sin 3 Sin 5
Calcular tambin BC de los tringulos DAB y ABC:
AB
AD
Sin 7
Sin 4
and BC
AB
Sin 1
Sin 6
Sustituyendo:
BC
AD Sin 1 Sin 7
Sin 4 Sin 6
Sin 1 Sin 7
Sin 4 Sin 6
AD
Sin 2 Sin 8
Sin 3 Sin 5
Considrese el problema del cuadriltero en la Fig. 12.4 que fue resuelto antes por
la computadora usando el mtodo de la ecuacin de observacin. A continuacin,
la formacin de la ecuacin de condicin lateral (14k) para ese cuadriltero, en
forma tabular:
Angle
Value
Log Sin
1
3
5
7
42 35
290
79 54421
21 29239
31 20458
Diff
l
9.8304381
9.9932329
9.5638824
9.7161751
2.29
.38
5.34
3.46
Angle
Value
Log Sin
Diff
l*
2
4
6
8
87 35106
18 28224
39 01354
39 34279
9.9996145
9.5008618
9.7991197
9.8041938
.08
6.30
2.60
2.54
39.1037285
39.1037898
0
BL
0
2.29
0
0
0
1
1
1
1
0
0
0
1
1
1
1
1
2
3
4
5
6
7
8
La matriz W es:
W
W
W
W
1
2
3
4
180 (1 3 4 6) 08"9
180 (3 5 6 8) 09"3
180 (2 5 7 8)
11"8
61.3
0 ( 61.3)
0. 0 0 0 0 0 0 0
0. 0 0 0 0 0 0 0
0. 0 0 0 0 0 0 0
0. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
-0. 0 0 0 0 0 0 0
5. 3 4 0 0 0 0 0
-2. 5 4 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
0. 0 0 0 0 0 0 0
0. 0 0 0 0 0 0 0
0. 3 8 0 0 0 0 0
-2. 6 0 0 0 0 0 0
MATRIZ W
-8. 9 0 0 0 0 0 0
61. 3 0 0 0 0 0 0
-9. 3 0 0 0 0 0 0
11. 8 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
MATRIZ P
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
1. 0 0 0 0 0 0 0
RESIDUALES
2. 0 9 8 2 5 2 6
5. 0 3 4 6 3 9 0
-4. 1 8 1 8 0 5 9
-1. 4 1 5 6 2 4 2
1. 7 5 3 0 3 4 2
-5. 4 0 0 8 2 2 5
6. 4 8 2 7 3 2 6
-1. 4 7 0 4 0 5 8
14.8 Deduccin de la ecuacin para resolver ecuaciones de condicin por
mnimos cuadrados
14.8.1 Solucin por los multiplicadores de Lagrange / lagrangianos
Suponiendo que tenemos un sistema de ecuaciones de condicin tal como el de
(14a). Si desarrollamos estas ecuaciones tendremos:
.............. b2 m
b1m vm
m
m
b2 m vm
.
........(14 )
.
.
cq bq1 1 bq1v1 bq 2
bq 2 v2 ............... bqm
bqm vm
.. ........(14 )
w1 c1 (b11 1 b12
w2 c2 (b21 1 b22
............... b1m
.............. b2 m
.............. bqm
.
.
.
wq cq (bq1 1 b22
pv 2 min imun.................................................................(14n)
Ya que las (14m) todas son iguales a cero, podemos multiplicar cada ecuacin por
una constante (-2) y, por un "multiplicador lagrangiano", o correlacionar (K) y
sumarlos todos a la ecuacin (14n) sin cambiar el valor de
.............................................(14 p)
.
.
Dividiendo a travs de las (14p) por 2, resolviendo para las v, y lueqo sustituyendo
estos valores v en las (l4m), obtenemos:
................
p1
p1
p1
p2
p2
.............
b12 k q bq 2
p2
...........................
b1m k1b1m
pm
b1m k q bqm
b1m k 2b2 m
..............
w1 0
pm
pm
.
.
.
.
..............
p1
p1
p1
p2
p2
b22 k q bq 2
p2
...............................
..............
b2 m k q bqm
pm
w2 0
.
.
.
.
b2 m k1b1m b2 m k1b2 m
pm
pm
bq1k1b11
p1
bq1k 2b21
bq 2 k q bq 2
p1
p2
...............
bq1k q bq1
p1
bq 2 k q bq1
................................
..............
bqm k q bqm
pm
p2
bqm k1b1m
pm
bq 2 k 2b22
p2
.............
bqm k 2b2 m
p2
wq 0
Las ecuaciones (14q) son q en nmero e implican como incgnitas slo a los
correlativos o valores K, K1, K2.., Kq. Estas ecuaciones se resuelven para los
correlativos, con lo cual esos valores correlativos luego se sustituyen de nuevo en
las ecuaciones (l4p) y se calculan los residuales. Teniendo los residuales, que
son valores de mnimos cuadrados, podemos aplicarlos a las observaciones
originales para calcular las observaciones corregidas. De stas corregidas, uno
puede calcular los valores ajustados para los parmetros desconocido requerido.
14.8.2 Solucin por mtodos matriciales
El procedimiento bsico anterior que acabamos de presentar, puede realizarse
usando mtodos matriciales. Las ventajas obviamente son la conveniencia de
expresar las ecuaciones y el procedimiento sistemtico computacional que se
adaptan idealmente a cmputos mediante computadoras electrnicas.
Comenzando otra vez con el sistema de ecuaciones de condicin (14),
representamos este sistema en forma matricial como:
B (L + V ) = C . . . . . . . . . . . . . . . . . . . ( 14r)
De la cual despus del desarrollo tenemos:
BL + BV = C
Y reordenando resulta:
BV = C - BL
Let BV = W = C BL.
Entonces BV - W = O. . . . . . . . . . . . . . . . . . . . . . (14 S)
' V T PV 2 K T ( BV W ) . . .. .. . . . . . . . . . . . . . . . . . . . . . (14v)
La expresin (14v) es equivalente a la ecuacin algebraica (14o). Para hacer
cumplir la minimizacin en la solucin por mnimos cuadrados, tomamos las
'
derivadas parciales de con respecto a cada residual y se establece el resultado
igual a cero. Esto se expresa en forma matricial como:
'
2V T P 2 K T B 0
V
. .. . . . . . . . . . . . . . . . . . . . . (14w)
V P 1 BT K . . . . . . . . . . . . . . .(14x)
Sustituyendo (14x) en (14s) tenemos:
B P 1 B T K W 0 .. . . . . . . (14y)
La ecuacin (14y) es equivalente a la expresin algbrica (14q).
Resolviendo (14y) para K, los correlativos, resulta:
K ( B P 1 B T ) 1 W
CAPITULO 15
AJUSTE DE POLIGONALES POR CONDICIONES
15.1 Introduccin
En cualquiera poligonal cerrada completamente levantada, se pueden cumplir en
los ajustes las tres siguientes condiciones: 1) cierre de partida o cierre en
coordenadas X; 2) cierre de latitud o cierre coordenadas Y; y 3) cierre de
acimut/acimutal. Estas condiciones existen sin importar cuantos lados estn
presentes en la poligonal. La manera de formar esas ecuaciones de condicin se
da en los artculos siguientes.
15.2 Las condiciones de partida Y latitud
La Fig. 15.1 ilustra una lnea IJ, y con referencia a la figura se pueden escribir las
ecuaciones siguientes:
X D sin . . . . . . . . . . . . . . . . . . . . . . (15a)
Y D cos . . . . . . . . . . . . . . . . . . . . . . . . (15b)
Fig. 15.1
En las 15a) y (15b) X y Y son la partida y latitud respectivamente del
curso, D es la longitud o largo de la lnea, y es el acimut medido del curso.
Tomando la derivada de la ecuacin (15a), un pequeo error de partida d X
d X Yd
X
dD . . . . . . . . . . . . . . . . . (15d)
D
d Y Xd
dD / D . . . . . . . . (15f)
vs
d X Yva
X
vs . . . . . . . . . . . . . . (15g)
D
d Y Xva
X
vs . . . . . . . . . . . . . (15h)
D
va2
, el Y que
Fig 15.2
( X 2 l X 1 )vs1
D1,2
( X 4 l X 3l )vs3
D3,4
( X 3l X 1l )vs2
D2,3
( X 5 l X 4 l ) v s4
D4,5
. . . . . . . . (15i)
Simplificando la (15i):
n 1
n 1
( X j X i )vsi
i 1
i 1
Dij
( X n X n ) (Yn Yi )vai
l
Fig 15.2
. . . . .. . . .
(15j)
( X 2 l X 1 )vs1
D1,2
( X 4 l X 3l )vs3
D3,4
( X 3l X 1l )vs2
D2,3
( X 5 l X 4 l ) v s4
D4,5
. . . . (15i)
Simplificando la (15i):
n 1
n 1
( X j X i )vsi
i 1
i 1
Dij
( X n X n l ) (Yn Yi )vai
. . . . . . . . . . (15j)
( Yn - Yn ) =
n 1
i 1
(Y j Y1 )
n 1
( X n X 1 ) v a1
i 1
1)
aij
vs1...................
(15k)
F I (1 va1 ) ( 2 va 2 ) ( 3 va 3 ) ( 4 va 4 ) ( 5 va 5 )............(15l )
En la ecuacin (15l ) i A j 180 , I es el acimut de control inicial y F es el
acimut de control final o de cierre. La (15l ) puede simplificarse como sigue:
n
i 1
i 1
F I i vai ................................................................................(15 m)
Los valores w en las (15 p ) y (15m) son ngulos de desviacin, y todas las
unidades estn en radianes.
15.4 Ejemplo del problema del ajuste de poligonales
Para ilustrar los principios discutidos en los artculos precedentes. La muy sencilla
poligonal presentada en el Art. 13.4 la cual se ajust por el mtodo de la ecuacin
de observacin, ser ahora ajustada por el mtodo de la ecuacin de condicin.
Como abajo se indica, las soluciones tabulares son muy tiles para formar las
ecuaciones de condicin.
15.4.1 Condicin acimutal
F I 1 va1 va 2 va 3 ................................................................(15 n)
i 1
Estacin
E
A
B
C
D
01
0.000291rad.
3438
ngulos
observados (A)
ngulos de
desviacin (W)
24000
15000
24001
6000
3000
6001
i 9001'
Acimuts
preliminares
000
6000
3000
9001
( X 'b X a )
( X 'c X 'b )
v s1
vs 2 ................(15o)
Dab
Dbc
(Y 'b Ya )
(Y 'c Y 'b )
v s1
vs 2 ..................(15 p )
Dab
Dbc
Es de utilidad una solucin tabular para calcular los coeficientes y cierres para las
(15o) y (15 p ) , a saber:
Estaci
n
Acimut
prelimina
r
Distancia
medida
Partida
Latitud
6000
200.00
173.205
100.00
3000
100.00
50.000
86.603
A
B
C
Coordena
das
X
1000.000
Coordena
da
Y
1000.000
1173.205
1100.000
1223.205
1186.603
0.866
Dab
200 .00
X 'c X 'b
50 .000
0.500
Dbc
100 .000
0.500
Dab
200.000
Y 'c Y ' b
86.603
0.866
Dbc
100.000
w 3 Yc Y 'c 1186.500 1186.603 0.103ft .
B=
v a1
va 2
va 3
V=
v s1
v
s2
vs 3
0.000291
W = 0.205
0.103
A
B
C
Ntese
Residual
-49
-17
+6
Valor ajustado
2395911
1495943
2400106
200.00
100.00
-0.11
-0.12
199.89
99.88
Distancia
Coordenadas Coordenadas
Partida Latitud
ajustada
ajustadas X
ajustadas Y
1000.00
1000.00
595911 199.89 173.09 99.99
1173.09
1099.99
295854
99.88
49.91 86.51
1186.50
1223.50
que las coordenadas ajustadas para el punto B concuerdan exactamente
Estacin
Acimut
ajustado
Valor medido
2400000
1500000
2400100
i 1
i 1
Fig. 15-4
MATRIZ B
1.0000000
0.0000000
186.60300
0.8660000
-223.2050
0.5000000
1.0000000
0.0000000
86.603000
0.5000000
-50.00000
0.8660000
1.0000000
0.0000000
0.0000000
0.0000000
0.0000000
0.0000000
MATRIZ W
[-0.0002910
-0.2050000
-0.1030000]
MATRIZ P
7300000.00
400.000000
47300000.00
156.0000000
47300000.00
0.000000000
RESIDUALES
0.0002350
0.0000834 Radianes
0.0000275
0.1069571
Pies
0.1225815
Observaciones
Angulo
PQA
QAB
ABC
BCR
CRS
Valor
8517
26711
9434
17930
30803
X r 11 .813,99
Yq 10 .000
Yr 10 .458,84
Desviaciones estndar:
ngulos medidos = 60"
Distancias medidas = 0.20 pies
Curso
QA
AB
BC
CR
Longitud
453.85
621.20
598.35
671.91
CAPITULO 16
TRANSFORMACION CONFORME DE COORDENADAS BIDIMENSIONALES
16.1 Introduccin
Se usa tal transformacin para transformar las coordenadas de puntos de un
sistema de coordenadas planas a otro sistema paralelo de coordenadas planas.
La caracterstica de una transformacin conforme es que mantiene una forma
cierta/verdadera entre las posiciones relativas de los puntos transformados. Una
tpica situacin en agrimensura donde se utiliza dicha transformacin, es la de
colocar dos levantamientos separados (efectuados en distintos sistemas de
coordenadas) en uno de los sistemas.
Las transformaciones conformes de coordenadas bidimensionales involucran un
proceso de tres pasos: (1) hacer el cambio utilizando una constante
(desmultiplicar), para crear iguales tamaos de los dos sistemas de coordenadas;
(2) rotacin, para hacer paralelos los ejes de referencia de los dos sistemas; y (3)
traslacin, a fin de crear un origen comn para los dos sistemas de coordenadas.
Para efectuar dicha transformacin, es necesario tener un mnimo de 2 puntos
comnmente llamados puntos de control, cuyas coordenadas son conocidas en
ambos sistemas. No obstante, es deseable tener ms de dos puntos de control.
16.2 desarrollo de las ecuaciones
Las Figs. 16.1 y 16.2 ilustran a dos sistemas de coordenadas independientes con
tres puntos de control: A, B y C, se requieren transformar los puntos adicionales 1,
2, 3,.., n del sistema x-y al sistema X-Y. las ecuaciones necesarias se
deducen como sigue:
Paso 1 Desmultiplicacin.
Para hacer las longitudes de lneas segn definidas por las coordenadas x-y,
iguales a sus longitudes segn definidas por las coordenadas X-Y, es necesario
multiplicar las coordenadas x-y por un factor de escala s. las coordenadas
desmultiplicadas x y y pueden, por tanto, definirse como:
Paso
2
x = sx
Rotacin
y = sy. (16a)
Fig. 16.3
Paso 3 Traslacin
Para llegar finalmente a las coordenadas X y Y de cualquier punto, es necesario
trasladar desde el origen del sistema X-Y hasta el origen del sistema X-Y.
Haciendo referencia a la Fig. 16.3, se indica que esta traslacin para cualquier
punto (vg., punto 4), se logra sumando factores de traslacin como sigue:
X 4 X 'TX
Y4 Y 'TY .................................................................................................(16 c )
.. (16f)
A=
a
X = b
c
d
XA
Y
A
XB
L=
YB
X
C
YC
VXA
VYA
VXB
V=
VYB
V
XC
VYC
E
1,049,422.40
1,049,413.95
1,049,244.95
N
51,089.20
49,659.30
49,884.95
X
121.622
141.228
175.802
174.148
513.520
754.444
972.788
1068.118
Y
-128.066
187.718
135.728
-120.262
-192.130
-67.706
120.994
161.624
PUNTOS DE CONTROL
PUNTO
A
B
C
ABCISA
1049422.400
1049413.950
1049244.950
ORDENADA
51089.200
49659.300
49884.950
COORDENADAS MEDIDAS
PUNTO
X
A
121.622
B
141.224
C
175.802
1
174.144
2
513.520
3
754.444
4
972.788
5
1064.118
Y
-124.066
187.714
135.724
-120.262
-192.130
-67.706
120.994
161.624
ABCISA
1049422.404
1049414.051
1049244.845
1049187.361
1047637.713
1046582.113
1045644.713
1045224.845
ORDENADA
51089.171
49659.223
49885.056
51040.629
51278.829
50656.241
49749.336
49541.807
ECUACIONES DE TRANSFORMACION
X TIERRA A1 A2 ( XMEDIDA ) A3(YMEDIDA )
X
3000.00
1000.00
3000.00
Y
2000.00
3000.00
4000.00
x
1767.77
707.11
1060.66
1414.21
353.52
1414.21
353.52
y
1060.66
707.11
1767.77
1414.21
1060.66
-353.52
-353.52
2) Para el problema 16.1, cules son los valores de: (a) el ngulo de rotacin,
(b) el factor de escala, (c) la traslacin en X, y (d) la traslacin en Y?
CAPITULO 17
ELIPSES DE ERRORES
17.1 Introduccin
Como se demostr en el capitulo 10, despus de realizar los ajustes por mnimos
cuadrados para obtener las posiciones ms probables de los puntos en un sistema
de coordenadas de referencia X-Y, las desviaciones posicionales estndar SX i y
SYi para el punto i en las direcciones X y Y pueden calcularse fcilmente de los
elementos de la matriz de covariacin. Esas desviaciones estndar proporcionan
estimaciones de errores en las direcciones de los ejes de referencia que tienen un
68 por ciento de probabilidad. En forma grfica, representan las medias
dimensiones de un rectngulo de error estndar alrededor de cada punto.
sen 2tQXY
2
2
2
2
2
QYY QYY
Q
Q XX
QUU
Q XX
Finalmente,
QUU
D. Con una manipulacin similar, todos los elementos de la matriz QZZ pueden
representarse como: ecuacin (17h)
QYY Q XX QYY Q XX
cos 2t Q XY sen 2t
2
2
Q Q
YY
sen 2t Q XY cos 2t
XX
2
QXX QYY
sen 2t Q XY cos 2t
2
2
Reduciendo:
(QYY Q XX ) sen 2t 2Q XY cos 2t
sen 2t
2Q XY
tan 2t.........................................................................................(17 i )
cos 2t Q yy Q XX
Cuadrante
1
2
3
4
Permita que K
QYY QXX
cos 2t
QYY QXX
.........................................................................(17 j )
K
2QXY
Y de la ecuacin ( 17i ) : sen 2t
.....................................................(17k )
K
Luego: cos 2t
QYY QXX 2
2Q
) ( XY ) 2
K
k
Resolviendo para K:
K (QYY QXX ) 2 ( 2QXY ) 2
1 / 2 ........................................................(17l )
QYY QXX
Q QXX QYY QXX
Q
( YY
)(
) QXY (2)( XY )......................(17m)
2
2
K
K
Reduciendo:
QUU
QUU
Finalmente:
QYY QXX
1 (Q QXX ) 2
YY
2(QXY ) 2
2
K
2
QYY Q XX
2
1 K2
K 2
QUU
1
(QYY QXX k )..........................................................................................(17n)
2
tan 2t.........................................................................................(17 o)
cos 2t (QYY QXX )
G. De previas consideraciones:
QYY QXX )
cos 2t
Q QXX
cos 2t ( YY
)..........................................................................................................(17 p )
K
K (
2QXY
....................................................................................................................(17q )
K
(QYY Q XX ) 2 4(Q XY ) 2
K2
K (QYY Q XX ) 2 4(Q XY ) 2
1/ 2
QVV
QYY QXX
Q QXX
Q QXX
Q (2)QXY
( YY
)( YY
) XY
2
2
K
K
Reduciendo:
QVV
(Q XY ) 2
QYY QXX 1 (QYY QXX ) 2
2
K
2
QVV
QYY QXX 1 K 2
2
K 2
Finalmente:
QVV
1
(QYY QXX K ).............................................................................................(17r )
2
La ecuacin (17r) produce un valor numrico de QVV que cuando se multiplica por
s 2 da la variacin a lo largo del eje V. La raz cuadrada de esta variacin es el
semieje menor de la elipse.
(1)
2Q XY
sen 2t
QYY Q XX
cos 2t
(2)
QUU
1
(QYY QXX K )
2
(3)
QVV
1
(QYY QXX K )
2
Donde:
K (QYY Q XX ) 2 4(Q XY ) 2
1 / 2
La ecuacin (1) da el ngulo t, que el eje U hace con el eje Y. La ecuacin (2) da
el valor numrico para QUU que al multiplicarse por s 2 , da la variacin a lo largo del
eje U. La raz cuadrada de la variacin es el semieje mayor para la elipse de error
estndar.
La ecuacin (3) da el valor numrico para QVV que al multiplicarse por s 2 da la
variacin a lo largo del eje V. la raz cuadrada de esta variacin da el semieje
menor de la elipse de error estndar.
Ntese que no es menester buscar las funciones trigonomtricas, excepto para 2t
del valor numrico obtenido en la ecuacin (1).
Q XX
X=
dX 2
dX
3
dY2
dY3
1)
tan 2t
2( 1.66 )
2.333;2t 293 40' ; t 146 50'
2.634 1.199
2)
1/ 2
K 2.059 11 .022
1 / 2
13 .08 3.617
3)
SU S 0 QUU 0.020
SV S 0 QVV 0.020
4)
5)
S X 0.020 1.199 0.022 ft .
6)
SY 0.020 29364 0.033 ft .
2)
2(0.460 )
0.270;2t 15 08' ; t 7 34'
3.962 0.583
1 / 2
12 .26 3.502
3)
SU 0.020
SV 0.020
4)
5)
S X 0.020 0.583 0.015 ft .
6)
SY 0.020 3.962 0.040 ft .
Fig. 17.3
Elipses de error estndar para el cuadriltero trilaterado del Art. 11.9
17.3.2 Calcular la elipse de error estndar para el punto B del problema de la
poligonal del Art. 13.4
A. Para el ajuste, S0 1.76 pies , y las matrices X y Q son:
dX
X=
dY
Q XX Q XY 0.000553 0.000622
QB
=
0.000622 0.000872
Q XY QYY
1)
2(0.000622 )
0.000872 0.00553
2t 75 .6
tan 2t
t 37 .8
2)
1/ 2
3)
SU 1.76
SV 1.76
4)
5)
S X 1.76 0.000553 0.41 ft .
6)
SY 1.76 0.000872 0.052 ft .
Fig. 17.4
Elipse de error estndar para el ejemplo de poligonal Art 13.4
Estadstica F
200
19.0
9.55
6.94
5.79
4.10
3.68
3.49
3.39
3.19
1. Las estaciones Sand y Birch en ambas figuras tienen las precisiones ms altas.
Claro est que esto se esperaba, debido a su proximidad a la estacin de
control Bug y por la densidad de observaciones hechas a esas estaciones.
2. El gran tamao de las EE en las estaciones White, Baever, Bunker y Bee de la
Fig. 17.6 indica que estas estaciones tienen precisiones ms bajas.
3. Las estaciones White y Schutt (Fig. 17.5) tienen precisiones relativamente
altas de este a oeste y precisiones relativamente bajas de norte-sur. El examen
de la red indica que esto era de esperar.
4. Las estaciones Beaver y Bunker (Fig. 17.5) y White (Fig. 17.6) tienen una
precisin relativamente baja este-oeste. Una vez ms el examen de de las
redes indican que esto se esperaba.
5. La estacin Bunker (Fig. 17.69) tiene una precisin uniforme total en todas las
direcciones segn se ilustra por la casi forma circular de sus EE.
6. Las EE generalmente ms pequeas en la Fig. 17.5 indican que el
levantamiento por trilateracin espera producir precisiones superiores al
levantamiento por triangulacin de la Fig. 17.6.
Estos ejemplos sirven para ilustrar algunos de los valores de las EE de estacin.
Tales observaciones se hicieron fcil y rpidamente mediante una comparacin de
las EE. Esa misma informacin hubiera sido difcil, sino imposible, de averiguar
slo con las desviaciones estndar en las direcciones X y Y.
Variando el levantamiento, sera posible el encontrar finalmente un diseo que
proporcionar ptimos resultados en trminos de llenar los requisitos de una
precisin uniformemente aceptable con la mayor economa posible.
Fig. 17.5
Fig. 17.6
Tarea autodidctica No. 17
CAPITULO 18
ENCAJADO/ACOMODAMIENTO DE CURVAS POR MNIMOS
CUADRADOS
18.1 Introduccin:
Con frecuencia en trabajos de ingeniera, es deseable o necesario
encajar/acomodar una lnea recta o curva en un juego de datos. Los datos normalmente se deducen de mediciones o de un experimento. El problema esta en la
determinacin de la ecuacin de la lnea o curva que mejor encaje los datos. La
decisin de usar una lnea recta, parbola u otra curva de orden ms alto,
generalmente puede tomarse despus de trazar los datos y estudiar su forma.
Puede evaluarse la fineza del encaje de la lnea o curva seleccionada basndose
en los residuales despus de haber computado la ecuacin.
(18a)
Fig. 18.1
Mxa + B = ya + vya
Mxb + B = yb + vyb
Mxc + B = yc + vyc
Mxd + B = yd + vyd
(18b)
xd
1
1
1
ya
M y b
X L
B y c
yd
vy a
vy
b
V
vy c
vy d
(18 d)
(18e)
x a2
2
xb
A x c2
2
xd
x 2
e
xa
xb
xc
xd
xe
1
1
1
1
A
X B
C
ya
y
b
L yc
yd
y e
va
v
b
V vc
v d
ve
Distancia radial, r
(mm)
0.000
0.000
20.072
0.004
40.855
0.18
63.155
0.047
88.034
0.062
116.995
0.035
152.472
-0.064
k 3 1.97388 10 11
k 2 1.94801 10 7
k 4 0.43934 10 15
Usando estas k en la (18g), pueden calcularse fcilmente las distorsiones radiales
del lente para cualquier valor de r. Para comprobar la validez de esas k, se
calcularon las distorsiones usndolas a incrementos de 10 mm de r, y los
resultados se trazaron segn la Fig. 18.4 (curva quebrada). Superpuesto en la
misma figura se halla el trazo de los datos reales dados para el problema (lnea
Longitud (pies)
B
464.8
503.1
820.0
917.3
329.8
287.4
345.9
Estacin
18017'
C
D
E
F
G
H
ngulo a la derecha
179051'
17928'
18033'
17910'
17959'
18102'
AB
2) Calcular una parbola con e1 mejor encajado que se encaje a los siguientes
datos obtenidos en un levantamiento de una curva vertical existente. La curva
empieza en la estacin 10+00 y termina en la estacin 18+00.
Estacin
10+00
11+00
12+00
13+00
14+00
Elevacin
51.2
49.5
48.2
47.3
46.8
Estacin
15+00
16+00
17+00
18+00
Elevacin
46.9
47.3
48.3
49.6
APNDICE A
LA CURVA DE DISTRIBUCIN NORMAL DEL ERROR
que esto sea la ecuacin para la curva de probabilidad, vamos a determinar ahora
la forma de f(x). Las probabilidades de ocurrencia de errores dentro de las gamas
de (x1 y x1 +dx), (x2 y x2 + dx), etc., son f(x1) dx1 , f(x2) dx2,...., f(xn) dxn. El rea
total bajo la curva de probabilidad representa la probabilidad total o simplemente el
entero uno. Luego para un nmero finito de posibles errores,
f(x1) dx1 + f(x2)dX2 + .... + f(xn)dXn = 1
(Aa)
f ( x)dx
f ( x ) dx
ydx
(Ab)
Supongamos que hemos medido una cantidad M que es igual a alguna funcin de
n parmetros desconocidos z1 , z2,...., Zn para que M = f(z1, z2.zn). Deje
tambin que x1. x2,...., xn sean los errores de observaciones m1 , M2, ...., Mm y
deje que f(x1) dx1 , f(x2) dx2, .f(xm) dxm sean las probabilidades de los errores
cayendo dentro funciones gamas de (x1 y x1 + dx1) , (x2 y x2 + dx2), etc. Por la
ecuacin (4c) la probabilidad P de 1a ocurrencia simultnea de todos estos errores
es igual al producto de las probabilidades individuales, por ende:
P f ( x1 )dx1 f ( x 2 )dx2 ............... f ( x m )dxm
.......
0
p Z 1
f ( x1 ) dx1 dz 1 f ( x 2 ) dx 2 dz 1
f ( x m ) dx m dz 1
1 P
1 df ( x1 ) dx1
1 df ( x 2 ) dx 2
1 df ( x m ) dx m
.......
0
p Z 2
f ( x1 ) dx1 dz 2 f ( x 2 ) dx 2 dz 2
f ( x m ) dx m dz 2
.
.
.
1 P
1 df ( x1 ) dx1
1 df ( x 2 ) dx 2
1 df ( x m ) dx m
.......
0
p Z n
f ( x1 ) dx1 dz n
f ( x 2 ) dx 2 dz n
f ( x m ) dx m dz n
df ( x )
dx
(Ae)
f `( x m ) dx m
f `( x1 ) dx1 f `( x 2 ) dx 2
..........
0
f ( x1 ) dz1
f ( x 2 ) dz1
f ( x m ) dz1
f `( x m ) dx m
f `( x1 ) dx1 f `( x 2 ) dx 2
..........
0
f ( x1 ) dz 2
f ( x 2 ) dz 2
f ( x m ) dz 2
.
.
.
f `( x m ) dx m
f `( x1 ) dx1 f `( x 2 ) dx 2
..........
0
f ( x1 ) dz n
f ( x 2 ) dz n
f ( x m ) dz n
(Af)
(Ag)
Luego para este caso especial, sustituyendo (Ag) y (Ah) en las ecuaciones (Af), se
reducen a una sola ecuacin:
f `( z * M 1 )
f ( z * M 1 )
f `( z * M 2 )
.........
f ( z * M 2 )
(Ai)
f `( z * M m )
f ( z * M m )
Y la ecuacin (Ai) para este caso especial es tambin general para cualquier valor
de m y para cualesquier valores observados m1, M2,...., Mm. Por tanto, digamos
que los valores de M son:
M 2 M 3 ....... M m M 1 M n
donde N se escoge por conveniencia como
N
M1 M 2
m
La media aritmtica es el valor ms probable para este caso de una sola cantidad
que ha sido observada varias veces; por tanto, f el valor ms probable en este
caso es:
z* ( M 1 M 2 ......... M m ) / m
z* M 1 ( m 1)( M 1 mN ) / m
z* M 1 mN N
z* M 1 N ( m 1)
z * M 1 N ( m 1) N (1 m)
(Aj)
N
Recuerde que
M1 M 2
m
de lo cual M1=mN + M2
Sustituyendo en (5j)
z * ( mN m 2 ) N (1 m )
z * M 2 N
M1 M 3 M1 M 4
.......etc .
m
m
z * M 3 N
z * M 4 N
etc.
0
f N (1 m)
f (N )
(Ak)
reordenando lo anterior:
f `N (1 m)
f `( N )
Nf (1 m)(1 m) f ( N ) N
xf ( x )
una constante = k
(A2)
df ( x)
dx
De lo cual
df ( x )
k dx
f ( x)
(Am)
Integrando (Am)
log e f ( x )
1
c1
2kx 2
f ( x ) e c1 e kx
/2
luego:
f ( x) Ce kx
/2
(An)
En (An), ya que f(x) disminuye mientras que x aumenta, luego K/2 debe ser
negativo. Cambiando arbitrariamente esta constante 1/2K a h2, resulta:
f ( x) Ce h
2 2
(Ap)
Ce
h2 x2
dx
C / h e t dt 1
2
, del cual:
C
1
h
h
C
(Aq)
e h
2 2
(Ar)
Rainsford.
APNDICE B
Programa general para el ajuste de minimos cuadrados mediante
ecuaciones de observacion
B.1 Introduccin al programa
Este programa de computadora dirige la solucin del algoritmo de mnimos
cuadrados dado en la ecuacin (9g). El algoritmo no ponderado de la (9c) tambin
puede resolverse si los pesos de unidad o uno (1) se ingresan (entrada). El
programa acomodara tantas como 30 observaciones involucrando 30 incgnitas.
Para usar este programa, es necesario calcular manualmente e ingresar en la
computadora las matrices A, L y P. Tal programa es especialmente til para
ajustar las redes de nivel como las de los ejemplos en los Arts. 9.2 y 9.4y
tambin para una variedad de otros problemas de ajuste por mnimos cuadrados.
Adems de calcular los valores ms probables para las incgnitas, el programa
dirige la computadora a resolver para los residuales, la desviacin estndar del
peso unitario y las desviaciones estndar estimadas de las cantidades ajustadas.
Fig. B.1
B.4 Salida
El Art. 9.5 brinda ejemplos resultantes de la salida de computadora de este
programa. Esta salida es para los ejemplos de la red de nivel "no ponderada'' y
"ponderada" en los Arts. 9.2 y 9.4, respectivamente. Los datos en la salida se
identifican en las hojas/copias impresas que son auto explicativas. Se da otro
ejemplo en el Art. 13.4 que es para el ajuste de una poligonal por a1 mtodo de la
ecuacin de observacin.
B.5 Listado de1 programa
En las pginas siguientes se ofrece el listado del programa. La ltima pgina del
mismo se compone de los datos de entrada para los problemas en los Art. 9.2, 9.4
y 13.4. Dichos listados de entrada deben estudiarse con cuidado junto con los Art.
B.2 y B.3 para asegurarse de ingresar correctamente tos datos en la computadora
al usar este programa
PROGRAMA GENERAL PARA EL AJUSTE DE MNIMOS CUADRADOS POR
LAS ECUACIONES DE OBSERVACIN NO PONDERADAS O PONDERADAS.
Programado por Pal R. Wolf. PH.D.
LENGUAJE DEL PROGRAMA FORTRAN
NOMBRES VARIABLES DE ENTRADA DEL PROGRAMA
NO = NUMERO DE SOLUCIONES ESTA PASADA
TITLE = CUALESQUIER NOMBRES O NMEROS IDENTIFICACIN DEL
TRABAJO
M AND N = EL NUMERO DE ECUACIONES E INCGNITAS
RESPECTIVAMENTE
A(I,J) = LA MATRIZ DE COEFICIENTE
EL(I,1) = LA MATRIZ CONSTANTE
OLL(I,J)' = LA MATRIZ DE PESO (LOS PESOS SE INGRESAN COMO 1 SI
LA SOLUCIN ES IGUALMENTE PONDERADA O NO PONDERADA)
DIMENSION A(30,309, EL(30,1), OLL(30,30). AT(30,30), AQ(30,30),QX1X(30,30),
AQL(30,l), X(30,l), V(30,l), VAR(30),title(1 6)
LEER DATOS DE ENTRADA
READ 500,NO
DO 510 II= l,NO
PRINT 509
READ 710,TITLE
PRINT 710,TITLE
READ 500, M ,N
PRINT 500, M, N
DO 502 I = 1,M
VAR(I)=SQRT (QXX(1,1)*S1GMA**2)
PRINT 90
PRINT 508,(VAR(I),I=l,N)
FORMAT(l6X,9H A MATRIZ)
FORMAT(11X, 9H L MATRIX)
FORMAT(10X, 14H WEIGHT MATRIX)
FORMAT(10X, 18H COVARIANCE MATRIX)
FORMAT(10X, 11H UNKNOWNS X)
FORMAT(10X, 12H RESIDUALS V)
FORMAT(10X, 5FI5, 3)
FORMAT(10X, 34H STANDARD DEVIATION OF UNIT WEIGHT)
FORMAT(11X, STANDARD DEVIATION OF ADJUSTED VALUES)
FORMAT(10X, 2I5)
FORMAT(8F10, 3)
FORMAT(10X, NORMAL EQUATIONS MATRIX)
FORMAT(10X, 8F14, 6)
FORMAT(1H1)
FORMAT(10X, 8F10, 3)
FORMAT(10X, 16A5)
CONTINUE
PRINT 509
STOP
END
EJEMPLO RED DE NIVEL NO PONDERADA, COMPUTOS DE AJUSTE
DE WOLF
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
1.
7
0.
0.
0.
0.
-1.
1.
1.
1.0
5,
1.6
3
0.
0.
1.
1.
0.
0.
-1.
10
6,
25
10
6,
13
1.
1.
1.
0.
68
10
4,
50
1.
70
1.0
1.
5,1
0
1.
1.
1.
.
EJEMPLO RED DE NIVEL PONDERADA, COMPUTOS DE AJUSTE DE
WOLF
7
3
1
0.
0.
.
1
0.
0.
.
1
0.
1.
.
1
0.
-1.
.
1
1.
0.
.
1
1.
0.
.
1
-1.
1.
.
1
10
10
1.7
.
1.0
6,
10
0.6
4,
0
5,
25
6,
8
50
1.6
13
1
6.
.
6.
6.
6.
6.
6.
EJEMPLO DE POLIGONACION, COMPUTOS DE AJUSTE DE WOLF
.866
10.498
515.7
5
0.500
-0.865
893.1
2299.8
2784.2
0.0
1.
10.0
1920.3
1027.2
0.19
-1.
156.0
-3.0
249.0
192.0
0.001
0.001
APNDICE C
PROGRAMA PARA EL AJUSTE POR MNIMOS CUADRADOS DE
REDES DE
TRILATERACION
C.1 Introduccin al programa
Este programa est diseado para ajustar levantamientos de trilateracin por el
mtodo de la ecuacin de observacin de manimos cuadrados. Podr acomodar
cualquiera configuracin de la red, sin importar su complejidad. Segn la
determinacin de medidas, puede manejar una red conteniendo hasta 20 lneas
medidas e incluyendo hasta 10 estaciones nuevas. Cualquier numero de
estaciones (control) fijas pueden existir en la red. El programa puede fcilmente
ampliarse para acomodar mayores cantidades de lneas medidas y estaciones
desconocidas sencillamente modificando las declaraciones/sentencias de
dimensin.
Conociendo los correctos datos de entrada segn los Arts. C.3 y C.4, este
programa dirige la computadora a resolver las aproximaciones iniciales para las
coordenadas de estacin desconocidas. Luego procede a formular las matrices A
y L de 1a ecuacin (9g). Sin embargo, la matriz de peso debe leerse en la
computadora.
Siguiendo a la formulacin de las matrices, e1 programa resuelve la ecuacin (9g)
por mnimos cuadrados de forma iterativa hasta lograr una convergencia
satisfactoria. Habiendo alcanzado la convergencia, computa los valores ms
probables para las coordenadas ajustadas de los puntos desconocidos, los
residuales en las longitudes medidas, la desviacin estndar del peso unitario, las
desviaciones estndar de las coordenadas ajustadas y las longitudes ajustadas de
las lneas medidas. Estos datos los lista la computadora en su salida.
C.2 Nombres y descripciones de las variables de entrada
En este programa se han asignado los siguientes nombres a las variables de
entrada:
NON - nmero de redes procesadas en esta pasada de computadora.
TITLE - identificacin del trabajo.
LL - numero de estaciones nuevas en la red.
MM - nmero de lneas medidas en la red.
NC - nmero de puntos de control en la red.
NAM1 - primeros 3 caracteres de nombres de estacin.
NAM2 - siguientes 3 caracteres de nombres de estacin; por tanto, los nombres de
estacin pueden tener hasta 6 caracteres. Pueden ser alfa o numricos. Deben
justificarse a la derecha en las columnas 2 al 7.
XX(1) - coordenada X del punto de control inicial No. 1.
YY(1) - coordenada Y del punto de control inicial No. 1.
SI -. longitud medida de la lnea 1.
AZD - grados de acimut de la lnea 1.
AZM - minutos de acimut de la lnea 1.
AZS - segundos de acimut de la lnea 1.
NP - nmeros de puntos de estaciones de control del dibujo de entrada.
XP coordenadas X de estaciones de control salvo la estacin No. 1.
YP - coordenadas Y de estaciones de control salvo la estacin No. 1.
S - longitudes medidas de las lneas en sucesin de la estacin No. 2 a No. 3, No.
3 a No. 4, No. 4 a No. 5, etc., hasta haber incluido a toda estacin nueva.
DEFD - esta variable, y las siguientes dos variables, se relacionan con las
porciones de grados, minutos y segundos, respectivamente, de los ngulos de
desviacin. Se ingresan en sucesin para las lneas de la estacin No. 2 a No. 3,
No. 3 a No. 4, No. 4 a No. 5, etc., hasta incluir a cada estacin nueva. Ntese que
un ngulo de desviacin derecho viene precedido por un signo ms (+), y un
ngulo de desviacin izquierdo esta precedido por un signo menos (-).
DEFM - la porcin de minutos de los ngulos de desviacin arriba descritos.
OEFS - La porcin de segundos de los ngulos de desviacin arriba descritos.
DIST - las longitudes medidas de las lneas.
IEND - el nmero del dibujo que corresponde a los extremos I de las lneas
medidas.
JEND - el nmero del dibujo que corresponde a los extremos J de las lneas
medidas.
ICOOE - nmeros de cdigo indicando a la computadora si est fijo el extremo I, el
extremo J o ningn extremo de cada lnea de entrada. Un ICODE de uno indica
que est fijo el extremo I, un dos indica que est fijo el extremo J y un cero indica
que ningn extremo de la lnea est fijo.
P - pesos de las lneas medidas.
C.3 Preparacin del lote de datos para la entrad
El croquis siguiente indica la manera de preparar las tarjetas para la lectura de los
datos de entrada en la computadora:
Fig. C.1
C.4 Discusin del procedimiento de entrada
C.4.1 Dibujo de la red
podran escalarse del dibujo con este propsito. Sin embargo, se prefiere el uso de
longitudes medidas porque esto producir valores ms exactos para las
aproximaciones iniciales de las coordenadas de estaciones desconocidas, lo que
ayudara a la computadora a lograr la convergencia con menos iteraciones y, por
tanto, ahorrando tiempo y dinero.
A fin de obtener una mejor comprensin sobre el procedimiento de numeracin de
estacin y e1 mtodo para dar entrada a los datos y calcular asi las coordenadas
iniciales, se resumen los ejemplos de las Figs. 11.5 y 11.6. En la Fig. 11.5, se
ingresaron tas coordenadas de la estacin 1 (Badger) junto con la longitud
aproximada (5870 pies) y el acimut , (80 00' 00") de la lnea 1-2 (Badgerwisconsin). Esto permite calcular las coordenadas de la estacin 2 (Wisconsin).
Luego se ingresan la longitud aproximada (3616 pies) de la Lnea 2-3 (WisconsinCampus) y el ngulo de desviacin (82" 00" 00") de la lnea 2-3. Esto permite
computar las coordenadas para la estacin 3 (Campus).
En la Fig. 11.6, se ingresaron las coordenadas de la estacin de control Lou junto
con la longitud y acimut de la lnea Lou-Paul para permitir calcular las
coordenadas de la estacin Pal. Luego se ingresaron la longitud de 1a lnea
Paul-Ram y el ngulo de desviacin "a" (ver Fig.) para permitir computar las
coordenadas para la estacin Ram. Despus se ingresaron la longitud de la lnea
Ram-Tom y e1 ngulo de desviacin "b" a fin de obtener coordenadas para Tom.
Luego se ingresaron la longitud de 1a lnea Tom-Jim y el ngulo de desviacin "c"
a fin de determinar las coordenadas para la estacin Jim. (Ntese que el ngulo
de desviacin "c" est a la izquierda o al sinistrorso, por lo que se le da un signo
menos (-.) en 1a entrada). Finalmente, se leyeron en la computadora la longitud
de la lnea Jim-Bitl y el ngulo de desviacin "d" a fin de obtener las coordenadas
para la estacin Bill.
C.4.3 Procedimiento de entrada
Dicho procedimiento, incluyendo los nombres de variables, secuencia de las
tarjetas de entrada, datos requeridos de tarjetas y tos necesarios formatos de
lectura se discuten en los Arts. C.2 y C.3. Ser de utilidad una explicacin
adicional.
La primera tarjeta le dice a la computadora cuantas redes se vienen resolviendo
en cualquiera pasada en particular. Este nmero, con el nombre variable NON, se
justifica a la derecha en la columna 5 de la tarjeta.
La segunda tarjeta contiene la identificacin del proyecto. Este nombre variable es
TITLE, y cualesquier datos alfa o numricos pueden perforarse dondequiera en
esta tarjeta desde la columna 2 hasta la 80.
La tercera tarjeta contiene las variables LL, MM y NC en tal orden. Ellas son el
nmero total de estaciones nuevas en la red, el nmero total de lneas medidas en
la red y el nmero total de puntos de control en la red, respectivamente. Esos
nmeros se perforan justificados a la derecha en 1as columnas 5, 10 y 15,
respectivamente.
El cuarto juego de tarjetas da los nombres de las estaciones en orden consecutivo:
un nombre de estacin por tarjeta. Se justifican a la derecha en la columna 7.
Pueden usarse caracteres alfa o numricos de hasta seis caracteres para
identificar cualquiera estacin.
La tarjeta siguiente da las coordenadas X y Y de la estacin No. 1, la longitud de la
lnea 1-2 y las porciones de los grados, minutos y segundos del acimut de la lnea
1-2. Estos nombres de variables son XX, YY, SI, AZD, AZM y AZS,
respectivamente. La variable XX se perfora en alguna parte dentro de las
columnas 1 al 15, YY en las columnas 16 al 30, SI en las columnas 31 al 40, AZD
de la 41 a la 50, AZM de la 51 a 1a 60 y AZS de la 61 a la 70.
El juego siguiente de tarjetas da los restantes nmeros de estacin y sus
coordenadas X y Y. Tales nombres variables son NP , XP y YP , respectivamente.
Se perfora una tarjeta para cada estacin de control (excepto la No. 1, que ya ha
sido ingresada en la computadora). En cada tarjeta el nmero de estacin se
justifica a la derecha en la columna 5, su coordenada X se perfora dentro de las
columnas 6 a 20 y su coordenada Y en las columnas 21 a 35.
El siguiente juego de tarjetas da las longitudes y ngulos de desviacin
necesitados para calcular las coordenadas iniciales de todas las estaciones
desconocidas excepto la No. 2. Los datos perforados en estas tarjetas constan de
longitudes de lnea, las porciones de grados, minutos y segundos del ngulo de
desviacin. Los nombres de estas variables son S, DEFD, DEFM y DEFS, y se
perforan dentro de las columnas 1 a 10, 11 a 20, 21 a 30 y 31 a 40,
respectivamente. Hay LL-1 de estas tarjetas, o una menos que el nmero total de
estaciones nuevas en la red.
El siguiente juego de tarjetas da a la computadora las longitudes medidas, el
nmero del extremo I y del extremo J de cada lnea, y un nmero de cdigo. Se
perfora una tarjeta para cada longitud medida. Los nombres de variables son
DIST, IEND, JEND e ICODE, respectivamente. Los extremos I y J de las lneas
corresponden a los nmeros asignados de estacin del dibujo. El cdigo le dice a
la computadora si la lnea tiene su extremo I fijo, en cuyo caso el cdigo es uno; si
est fijo el extremo J, en cuyo caso el cdigo es dos; o si ningn extremo est fijo,
en cuyo caso el cdigo es cero. La distancia se perfora en las columnas 1 a 10. El
y los nmeros ICODE se justifican a la derecha en las columnas 15, 20 y 25,
respectivamente.
E1 juego final de tarjetas contiene los pesos de las lneas medidas. El nombre de
la variable es P. Cada tarjeta puede contener hasta 8 pesos perforados en campos
de 10.
C.5 -Listado del programa
En las paginas siguientes se da el listado del programa. La pgina final consta de
los datos de entrada para los problemas de los Arts. 11.9 y 110.De esos datos de
entrada, se produjo la salida dada en aquellos artculos. Seria til comparar estos
datos de entrada con las Figs. 11.5 y 11.5.
AJUSTE DE TRILATERACION-COORDENADAS PLANAS MTODO
ECUACIN DE OBSERVACIN-MNIMOS CUADRADOS
POR PAUL R WOLF
IMPLICIT DOUBLE PRECISIN (A=M,0=Z)
DIMENSIN A(30,20),P(30),Q(20,20),AL(20),FL.(30),XX(15),YY(15)
DIMENSIN X(20),DIST(20},IEND(20),JEND(20),ICODE(20),V(30)
DIMENSIN TITLE(l6),NAM1(30),NAM2(30),EX(20),EY(20)
SIN(X)=DSIN(X)
C0S(X)=DCOS(X)
ABS(X)=DABS(X)
SQRT(X)=DSQRT(X)
LEER 10, NON
NON ES EL NUMERO DE REDES PROCESADO EN ESTA PASADA.
DO 695 IM=1,NON
KRUN=1
READ 11,TITLE
READ 10,LL,MM,NC
TITLE ES IDENTIFICACIN DEL TRABAJO , ETC.
LL ES NUMERO DE ESTACIONES NUEVAS
MM ES NUMERO DE DISTANCIAS MEDIDAS
NC ES NUMERO DE PUNTOS DE CONTROL EN LA RED
PRINT 999
PRINT 11,TITLE
PRINT 12,LL
PRINT 13,M
PRINT 16,NC
PI=4..*ATAN(1.)
RAD=180./PI
M=MM
N=2*LL
PRINT 670
PRINT 998
LEER NOMBRES ESTACIN
NTOT=LL+NC
DO 36 I=L, NTOT
READ 87,NAM1(I),NAM2(I)
CONTINU
CALCULAR APROXIMACIONES INICIALES PARA ESTACIONES .
READ 15,XX(1),YY(1),S1,AZD,AZM,AZS
XX (1) ES COORDENADA X DEL PUNTO DE CONTROL INICIAL
YY(1) ES COORDENADA Y DEL PUNTO DE CONTROL INICIAL SI ES
LONGITUD MEDIDA DE LINEA 1 AZD, AZM Y AZS SON GRADOS MINUTOS Y
SEGUNDOS DEL ACIMUT DE LA LINEA 1
IAZD=AZD
IAZM=AZM
PRINT 997,S1,IAZD,IAZM,AZS
AZ=(AZD+AZM/60.*AZS/3600.)/RAD
XX(2)=XX(1)+S1*SIN(AZ)
YY(2)=YY(1)+S1*C0S(AZ)
LEER COORDENADAS DE PUNTOS DE CONTROL ADEMAS DE NUMERO 1
IF(NC=1)925,925,924
CONTINU
DO 860 I=2,NC
READ 46,NP,XP,YP
XX(NP)=XP
YY(NP)=YP
CONTINU
IF(LL=1)738,738,739
CONTINU
DO 75 I=2,LL
II=I+1
READ 20,S,DEFD,DEFM,DEFS
DEFD, DEFM Y DEFS SON GRADOS, MIM., SEG DE NGULOS DE
DESVIACIN COMENZANDO CON EST. 2, + SI DESVIACIN DERECHA,
- SI DESVIACIN IZQUIERDA S ES LONGITUD MEDIDA COMENZANDO
DESDE EST. 2 HASTA 3
IDEFD=DEFD
IDEFM=DEFM
PRINT 997,S,IDEFD,IDEFM,DEFS
DEF=(DEFD+DEFM/60,*DEFS/3600.)/RAD
AZ=AZ+DEF
IF(AZ=2.*PI)72,71,71
71 AZ=AZ=2.*PI
GO TO 74
72 IF(AZ)73,74,74
73 AZ=AZ+2,*PI
70 CONTINU
XX(II)=XX(I)+S*SIN(AZ)
YY(II)=YY(I)+S*COS(A2)
75 CONTINU
738 CONTINU
PRINT 670
PRINT 996
DO 910 I=1,NTOT
910 PRINT 995,NAMI(I),NAM2(I),XX(I),YY(I)
FORMULAR LAS MATRICES DE COEFICIENTES Y CONSTANTES
DO 90 I=1,M
EL(I)=0,
DO 90 J=1,M
90 A(I,J)=0.
CONTINUE
PRINT 994
DIST (I) SON LONGITUDES CURSOS MEDIDOS EN SUCESIN IEND ES
NUMERO EST. EN EXTREMO I DEL CURSO JEND ES NUMERO EST. EN
EXTREMO J DEL CURSO ICODE IDENTIFICA TIPO DE LINEA
1 SI EXTREMO I DE LINEA ESTA FIJO
2 SI EXTREMO J DE LINEA ESTA FIJO
EN BLANCO SI NINGN EXTREMO DE LINEA ESTA FIJO
DO 99 I=1,MM
READ 40,DIST(I),IENDM.,JEND(I),ICODE(I)
IE=IEND(I)
JE=JEND(I)
99 PRINT 993,NAM(IE),NAM2(IE),NAM1(JE),NAM2(JE),DIST(I),ICODE(I)
70 CONTINU
P(I) SON PESOS DE OBSERVACIONES, EN SUCESIN, LONGITUDES
PRIMERO
READ 20,(P(I).I=L,M)
509 CONTINU
DO 100 I=1,MM
IJK=IEND(I)
IJL=JEND(I)
DO=(XX(IJK)=XX(IJL))**2+(YY(IJK)-YY(IJL))**2
DO=SQRT(DD)
KEND=I0END(I)=1
II=LL+KEND
LEND=JEND(I)=1
JJ=LL*LEND
IF(ICODE(I)=1)50,60,50
50 A(I,KEND)=(XX(IJK)=XX(IJL))/DD
A(I,II)=(YY(IJK)=YY(IJL))/DD
IF(ICOOE(I)-1)60.60,100
60 A(I,LEND)=-XX(IJK)-XX(IJL))/DD
A(I,JJ)=YY(LJK)-.YY(IJL))/DD
EL(I)=((DIST (I) DD)
PRINT 999
PRINT 660
DO 102 I=1,M
102 PRINT 45,(A(I,J)',J=1,N)
PRINT 661
PRINT 45.(EL(I),I=1,M)
LA RUTINA POR MNIMOS CUADRADOS
DO 1 I1,N
DO 1 J=1,I
Q(I,J)=0
DO 3 K=1,M
3 1Q(I,J)=Q(I,J)+A(K,I)*P(K)*A(K,J)
Q(J,I)=(Q(I.J)
1 CONTINU
.
DO 307 K=1,N
DO. 302 J=1,N
IF_(J=K)304,302 304
Q(K,J)=Q(K,J)/Q(K,K)
302 CONTINU
Q(K,K)=1./Q(K,K)
DO 307 I=1,N
IF(I=K) 305,307,305
DO 303 J=1,N
IF(J=K) 306,303,306
306 Q(I,J)=Q(I,J)-Q(I,K)*Q(K,J)
303 CONTINU
Q(I,K)=Q(I,K)*Q(K,J)
307 CONTINU
DO 2 J=1,N
AL(J)=O.
DO 2 K=1,M
2 AL(J)=AL(J)*A(K,J)*P(K)*EL(K)
DO 900 I=L,N
X(1)=0.
DO 900 J=1,N
900 X(I)=X(I)+Q(I,J)*AL(J)
PRINT 662
PFTINT 45,(X(I),I=1,N)
DO 51S I=1,LL
J=I+1
K=I+LL
XX(J)=XX(J)+X(I)
515 YY(J)=YY(J)+X(K)
PROBAR PARA LA NO CONVERGENCIA
KRUN=KRU+1
IF(KRUN=5)702,702,701
702 CONTINU
PROBAR PARA LA CONVERGENCIA AL DESEADO ESTNDAR
EXACTITUDNCOR=2*LL
I=1
507 CONTINU
IF(ABS(X(-I )) =.001) 505, 505,500
505 CONTINU
I=L+1
IF(I=NCOR)507,507,506 506
CONTINU
GO TO 508
506 CONTINU
IMPRIMIR LA MATRIZ DE COVARIACION
PRTNT 693
DO 692 I=1,N
692 PRINT 45, ((Q(I,J),J1L*N)
COMPUTAR RESIDUALES
DO 400 I=1,M
V(I)=O
DO 400 K=L,N
400 V(I)=V(I)+A(I,K)*X(K)
DO 401 I=L,M
V(I)=V(I)-EL(I)
COMPUTAR DESVIACIN ESTNDAR DEL PESO UNITARIO
SIG=O.
DO 404 I=I,M
002 SIG=SIS+V(I)**2*P(L)
DM=M
DN=N
SIGMA=SQRT(SIG/DOM-DN))
COMPUTAR DESVIACIN ESTNDAR DE COORDENADAS AJUSTADAS
EX(L)=0.
EY(L)=0.
ILL=L*LL
DO 410 I=ILL,NTOT
EX(I)=0.
EY(I)=0.
410 CONTINU
DO 415 I=2,ILL
IMN=I-1
IM2=IM1+LL
EX(I)=SIGMA*SQRT( (Q(IM1, IM1))
EY(I)=SIGMA*SQRT(Q(IM2*IM2))
415 CONTINU
COMPUTAR LONGITUDES DE LINEA AJUSTADAS
DO 403 I=I,MM
J=IEND(I)
J=IENO(I)
403 DIST(I)=SQRTC(XX(J)=.XX(K))**2(YY(J)=YY(X))**23
IMPRIMIR COORDENADAS AJUSTADAS Y SUS DESVIACIONES ESTNDAR
PRINT 999
PRINT 405
DO 413 I=1,NTOT
413 PRINT 995,NAMI(I),NAM2(L),XX(I),EX(I),YY(I),EY(I)
IMPRIMIR DESVIACIN ESTNDAR DEL PESO UNITARIO
PRINT 406,SIGMA
IMPRIMIR LONGITUDES DE LINEA AJUSTADAS Y SUS RESIDUALES
PRTNT 407
DO 444 I=1,MM
IE=IEND(I) ,
JE=JEND(I)
444 PRINT 408,NAMF (IE),NAM2(IE),NAM1(JE),NAM2(JE)DIS(I),V(I)
PRINT 999
701 CONTINU
695 CONTINU
PRINT 999
FORMATOS
10 FORMAT(16I5)
11 FOMAT(16A5)
12 FORMAT_(//,10X,20HNUMBER OF NEW STATIONS =,13)
13 FORMAT(//,10X,30MNUMBER OF DISTANCES MEASUREO A,I=13)
15 FORMAT(2F10.3)
16 FORMAT(//,10X,37MNUMBES OF CONTROL POINTS IN THE NET =,I3)
20 FORMAT(8FL0.3)
40 FORMAT(F10.3.3I5)
45 FORMAT(1X,8F12,5)
46 FORMAT(I5.2F15.3)
87 FOBMAT(1X,2A3)
405 FORMAT(//,11X,44HADJUSTED COORDINATES ANO
STANDARDEVIATIONS,/,9
IX,7HSTATION,8X,IHX.17X,2HSX,11X,IHY,15X,2HSY,//) '
406 FORMAT(/.20X,36H STANDARD DEVIATION OF UNIT WEIGHT =F6,3,3H
FT)
407 FORMATO(//,30X,30HADJUSTED LENGTHS AND
RESIDUALS,/,18X,4HL1NE,20X,6 IHLENGTH,11X',12HRESIDUAL(FT),//)
408 FORMAT(11X,2A3,3H . ,2A3,2X,2F20.3)
660 FORMAT(//,25X,12HTHE A MATRIX,//)
661 FORMATO(//,25X,12HTHE L MATRIX,//)
662 FORMAT(//,25X,L2HTME X MATRIX,//)
670 FORMAT(//)
693 FORMAT(//,25X,12HTHE Q MATRIX,//)
993 FORMAT(11X,2A3=3H ,2A3,F28.3,115)
994 FORMAT(1HL,40X,10HIMPUT
.DATA,//,IIX,4HLINE,23X,6HLENGTH,13X,AHCCOE I//)
995 FORMAT(9X,2A3,4F15.3)
996 FORMAT(7X, 44MFIXED COORDINATES AND INITIAL
APPROXIMATIONS,//,9X, 17HSTATION,AX,LHX,15X,1HY,//)
997 FORMAT(20X,F17.3, 1X,2I5,F83.3)
998 FORMAT(10X,'PRFLIMINARY LENGTHS AND ANGLES TO CALCLATE' 1'
INITIAL APPROXTMATIONS'-//30X, LENGTH,20X,'ANGLE'//)
999 FORMAT(1H1)
STOP
END
2
EJEMPLO CUADRILATERO
TRILATERADO
5
2
2
BADGER
WIS
CAMPUS
BUCKY
10000,000
10000,000
5780,000 80. 0. 0.
4
11820,756
6881,222
3616,000 82.
0.
0.
5870,302 1
2
1
7297,588
3616,434
5742,878
5123,760
1.
5
LOU
PAUL
RAM
TAM
JIM
BILL
PAT
JOHN
TIM
10993.
34955.
144.69.
22939.
30319,17
34954,73
31376,2
45986,22
10992,63
16,353,41
39446,64
17445,22
29712,49
37437,93
35162,5
35162,49
1
2
2
3
1.
3
3
4
4
1.
1
0
2
2
1.
1.
EJEMPLO DE
TRILATERACIN
17
4
833766,690
7830040,4
8886672,11
9887981,21
11.
120.
-170
43.
454925,95
467467,97
473334,4
453589,28
30.
00.
-15
00.
2
3
2
2
2
1
4
6
4
7
4
6
16353. 26.
00.
00.
0
00.
6
4
4
8
3
2
8
8
9
4
6
4
0
0
0
2
0
1
2
2
2
1
0
0
30.
00.
43007,65
41305,83
14468,7
19896,23
22939,17
1.
1.
1.
1.
1
6
4
2
6
1.
1.
6
7
5
5
5
1.
1.
1
2
0
0
0
1.
1.
1.
1.
1.
1.
1.
1.
APENDICE D
PROGRAMA PARA EL AJUSTE POR MINIMOS CUADRADOS
DE UN CUADRILATERO TRIANGULADO
D.1
Introduccin al programa
A fin de que este programa sea general para todos los cuadrilteros, el ajuste se
realizar en un sistema giratorio de eje provisional X-Y. El eje Y se coloca de
manera tal que corte a travs de los ngulos 4 y 8 segn muestra la Fig. 12.4.
Esto hace que tales ngulos concuerden con los casos (g) y (h) respectivamente,
de la Fig. 12.1. Esto es conveniente para permitir a la computadora que forme la
matriz L sin importar la orientacin del cuadriltero. Despus de terminar el ajuste,
las coordenadas ajustadas se vuelven a girar a travs del ngulo (360-a) al
sistema del eje N-E. La formula de rotacin para la rotacin inicial a travs de a, y
la rotacin final a travs de (360-a) son las ecuaciones (16e). Se obtiene mejor el
valor del ngulo trazando el cuadriltero a escala, dibujando el eje provisional Y, y
desmultiplicando su valor con un transportador. Slo se necesita darlo al grado
ms cercano.
D.3 Nombres y descripciones de las variables de entrada para este
programa
NO el nombre de ajustes al cuadriltero que se viene efectuando en esta pasada
TITLE el nombre o descripcin del problema de ajuste que se est resolviendo
M el nmero de ngulos observados (usualmente 8)
N el nmero de incgnitas (siempre 4)
YA la coordenada Y de la estacin de control A
XA la coordenada X de la estacin de control A
YD la coordenada Y de la estacin de control D
XD la coordenada X de la estacin de control D
ALPHA ngulo de rotacin
DEG la porcin de grados de los ngulos observados
AMIN la porcin de los minutos de los ngulos observados
SEC la porcin de los segundos de los ngulos observados
QLL la matriz de peso. Estos son los trminos P de la ecuacin (9g). Slo se
leen los elementos diagonales de la matriz de peso; por lo tanto, no se
asume correlacin entre ngulos observados. Se asignan valores cero a
todos los elementos no diagonales en el programa.
D.4
El siguiente croquis indica la manera de preparar las tarjetas para la lectura de los
datos de entrada en la computadora:
Fig. D.1
D.5
C
C
C
C
205
133
READ 110,(QLL(J,J),J=1,M)
BASE=SQRT((YA-YD)**2+(XA-XD)**2)
DIR=PI/2.-ATAN((YD-YA)/(XD-XA))
AB=(BASE*SIN(ANGLE(7)))/SIN(ANGLE(4))
BC=(AB*SIN(ANGLE(1)))/SIN(ANGLE(6))
YB=AB*COS(DIR-ANGLE(1)-ANGLE(2))+YA
XB=AB*SIN(DIR-ANGLE(1)-ANGLE(2))+XA
YC=BC*COS(DIR-ANGLE (1)-ANGLE (2)+PI-ANGLE(3)ANGLE(4))+YB
XC=BC*SIN(DlR-ANGLE(1)-ANGLE(2)+Pl-ANGLE(3)-ANGLE(4))+XB
PRINT 3
PRINT 9
PRINT 71,YB,XB,YC,XC
PRINT 3
PRINT 13
DO 20 I=1,M
20
PRINT 11,(QLL(I,J),J=1,M)
CALL MATINV(QLL,M,30)
3
CONTINUE
PRINT 4
PRINT 15,TITLE
PRINT 3
PRINT 14,II
C
C
GENERAR LAS MATRICES
C
EL(1,1)=3ANGLE(1)-ATAN((YB-YA)/(XR-XA))+ATAN((YC-YA)/(XC-XA))
EL(2,1)=ANGLE(2)-ATAN((YC-YA)/(XC-XA))-ATAN((YA-YD)/(XD-XA))
EL(3,1)=ANGLE(3)-ATAN((YB-YD)/(XD-XB) )+ATAN((YB-YC)/(XC-XB))
EL(4,1)=ANGLE(4)-ATAN((XD-XB)/(YB-YD))-ATAN((XB-XA)/(YB-YA))
EL(5,1)=ANGLE'(5)-ATAN((XC-XA)/(YC-YA))+ATAN((XC-XD)/(YC-YD))
EL(6.1)=ANGLE(6)-ATAN((YB-YC)/(XC-XB))-ATAN((YC-YA)/(XC-XA))
EL(7,1)=ANGLE(7)-ATAN((YB=YD)/(XD-XB))+ATAN((YA-YD)/(XD-XA))
EL(8,1)=ANGLE(8)-TAN((XD-XB)/(YB-YD))-ATAN((XC-XD)/(YC-YD))
DO 77 I=1,8
EL(I,1)=EL(I,1)*(648000./PI)
AC=SQRT((YC-YA)**2+(XC-XA)**2)
BD=SQRT((YD YB)**2+(XD-XB)**2)
CD=SQRT((YD YC)**2+(XD-XC)**2)
AB=SQRT((XB-XA)**2+(YB-YA)**2)
CB=SQRT((XC-XB)**2+(YC-YB)**2)
DO 500 I=1,M
DO 500 J=1,N
500
A(I,J)=0.
PRINT 3
PRINT 7
A(1,1)= (XB-XA)/AB**2
A(1,2)= -1.*(YB-YA)/AB**2
A(1,3)= -1.*(XC-XA)/AC**2
A(1,4)= (YC-YA)/AC**2
A(2,3)= (XC-XA)/AC**2
A(2,4)= -l.*(YC-YA)/AC**2
A(3,1)= (XD-XB)/BD**2-(XC-XB)/BC**2
A(3,2)= (YB-YD)/YD)**2-(YB-YC)/BC**2
A(3,3)= (XB-XC)/BC**2*(-1.)
A(3,4)= (YB-YC)/BC**2
A(4,1)= (XB-XD)/BD**2-(XB-XA)/AB**2
A(4,2)= -l.*(YB-YD)/BD**2+(YB-YA)/AB**2
A(5,3)= -l.*(XC-XA)/AC**2+(XC-XD/CD**2
A(5,4)= (YC-YA)/AC**2-(YC-YD)/CD**2
A(6,1)= (XC-XB)/BC**2
A(6,2)= (YB-YC)/BC**2
A(6,3)= -l.*(XC-XB)/BC**2+(XC-XA)/AC**2
A(6,4)= -l.*(YB-YC)/BC**2-(YC-YA)/AC**2
A(7,1)= (XD-XB)/BD**2
A(7,2)= (YB-YD)/BD**2
A(8,1)= -1.*(XD-XB)/BD**2
A(8,2)= -1.*(YB-YD)/BD**2
l.*(XC-XD)/CD**2
A(8,4)= (YC-YD)/CD**2
DO 78 I=1,M
DO 78 J=1,N
78
A(I,J)=A(I,J)*(648000./PI)
DO 10 I=1,M
10
PRINT 41,(A(I,J),J=l,N)
PRINT 3
PRINT 12
PRINT 11,(EL(I,J),I=1,M)
C
C
LLAMAR LA SUBRUTINA DE MINIMOS CUADRADOS
C
CALL MALIME(A,X,EL,QLL,V,QXX,M,N)
PRINT 3
PRINT 21
DO 30 I=1,N
PRINT 41,(QXX(I,J),J=1,N)
PRINT 31
A(8,3)= -
PRINT 11,(X(I,1),I=1,N)
PRINT 3
PRINT 32
PRINT 11,(V(I,1),I=1,M)
YB= YB*X(1,1)
XB= XB*X(2,1)
YC= YC+X(3,1)
XC= XC+X(4,1)
II=II+1
C
C
C
I=1
CONTINUE
If(X(I,1).GT.0.0001)GO TO 34
I=I+1
IF(N=I)36,201,20l
CONTINUE
ANGLE(1)= ATAN((YB-A)/(XB-XA))-ATAN((YC-YA)/(XC-XA))
ANGLE(2)= ATAN((YC-YA)/(XC-XA))+ATAN((YA-YD)/(XD-XA))
ANGLE(3)= ATAN((YB-YD)/(XD-XB))-ATAN((YB-YC)/(XC-XB))
ANGLE(4)= ATAN((XD-XB)/(YB-YD))+ATAN((XB-XA)/(YB-YA))
ANGLE(5)= ATAN((XC-XA)/(YC-YA))-ATAN((XC-XD)/(YC-YD))
ANGLE(6)= ATAN((YB-YC)/(XC-XB)+ATAN((YC-YA)/(XC-XA))
ANGLE(7)= ATAN((YB-YD)/(XD-XB))-ATAN((YA-YD)/(XD-XA))
ANGLE(8)= ATAN((XD-XB)/(YB-YD))+ATAN((XC-XD)/(YC-YD))
THETA=2.*PI-ALPHA
TEMP1=YB
TEMP2=XB
TEMP3=YC
TEMP4=XC
YB=TEMP1*COS(THETA)+TEMP2*SIN(THETA)
XB=TEMP2*COS(THETA)-TEMP1*SIN(THETA)
YC=TEMP3*COSfTHETA)+TEMP4*SIN(THETA)
XC=TEMP4*COS(THETA)-TEMP3*SIN(THETA)
C
C
IMPRIMIR LOS VALORES AJUSTADOS
C
PRINT 4
PRINT 15,TITLE
PRINT 3
PRINT 108
PRINT 71,YB,XB,YC,XC
201
901
576
902
C
C
C
1
2
3
4
6
7
8
9
11
12
13
14
15
CALL MATRAN(V,VT,M,1,30,1,1,30)
CALL MATMUL(VT,QLL,VTP,1,M,M,1,30,30,30,1,30)
CALL MATMUL(VTP,V,VTPV,1,M,1,1,30,30,1,1,1)
AM= M
AN=N
S=SQRT(VTPV(1,1)/(AM-AN))
PRINT 3
PRINT 35,S
PRINT 3
PRINT 900
DO 901 I=1,N
SD=SQRT(QXX(1,1)*S**2)
PRINT 111,SD
PPINT 3
PRINT 843
DO 576 I=1,M
DEG(I)= ANGLE(I)*(180./PI)
IDEG(I)= DEG(I)
BDEG(I)= IDEG(I)
AMIN(I)= (DEG(I)-BDEG(I))*60.
IMIN(l)= AMIN(l)
BMIN(I)= IMIN(l)
SEC(I)= (AMIN(I)-BMIN(I))*60.
PRINT 845,IDEG(I),IMIN(I),SEC(I)
KK=KK+1
IF(KK-NO)5,5,902
CONTINUE
FORMATOS
FORMAT(10X,'ROTATION ANGLE =',F5.1,' DEGREES')
FORMAT(10X,'OBSERVED ANGLES')
FORMAT(/)
FORMAT(1H1)
FORMAT(2I10)
FORMAT(10X,'A MATRIX')
FORMAT(10X,'CONTROL POINT CO-ORDINATES')
FORMAT(10X,'APPROX. CO-ORDINATES OF UNKNOWNS')
FORMAT(8F8.3)
FORMAT(10X, 'L MATRIX (VECTOR)')
FORMAT(10X,'VARIANCE - COVARIANCE MATRIX QLL')
FORMAT(10X,ITERATION No. ,I4//)
FORMAT(15A6)
21
31
32
35
41
71
108
110
111
171
843
845
900
C
C
C
C
C
C
C
4
5
2
6
7
8
3
1
END
SUBROUTINE MATINV(A,N,M)
DIMENSION A(M,M)
DO 1 K=1, N
DO 2 J=1, N
IF(J=K) 4,5,4
A(K,J)=A(K,J)/A(K,K)
CONTINUE
CONTINUE
A (K, K)= 1./(K,K)
DO 1 I=1,N
IF(I=K)6,1,6
DO 3 J=1, N
IF (J = K) 7, 8, 7
A (I, J) =A (I, J) = A (I, K) * A (K, J)
CONTINUE
CONTINUE
A(I,K)=-A(I,K)*A(K,K)
CONTINUE
RETURN
END
SUBROUTINE MATMUL(A,B,C,L,M,N,LA,MA,MB,NB,LC,NC)
DIMENSION A (LA,MA),B(MB,NB),C(LC,NC)
DO 1 I=1,L
DO 1 J=1,N
C(I,J)=0.
DO 1 K=1,M
C(I=J)= C(I,J)+A(I,K)*B(K,J)
RETURN
END
TRANSPONER MATRIZ B= A-TRANSPONER
SUBROUTINE MATRAN(A,B,M,N,MA,NA,MB,NB)
DIMENSION A(MA,NA), B(MB,NB)
DO 1 I=1,M
DO 1 J=1,N
B(J,I)= A(I,J)
RETURN
END
1
EJEMPLO DE TYRIANGULACIN, CALCULOS DE AJUSTE DE
WOLF
8
4
8448,9
0.
42.
42.1
1.
1.
4270.33
35.
18.
35.4
1.
29.
28.
31.
1.
4568.75
87.
22.4
20.
1.
35.
21.
45.8
1.
7610.58
10.6
29.
39.
1.
79.
23.9
34.
1.
54.
39.
27.9
1.
APENDICE E
PROGRAMA PARA EL AJUSTE DE MINIMOS CUADRADOS POR
ECUACIONES DE CONDICION
E.1
Introduccin al programa
El croquis siguiente indica la manera de preparar las tarjetas para la lectura de los
Fig. E.1
E.4
Salida
Las siguientes pginas contienen listados del programa y de los datos de entrad
para los problemas de Arts. 14.5, 14.6, 14.7 y 15.4.4.
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
FORMATO (15)
2DA TARJETA-TITLE CULESQUIER CARACTERES ALFA O NUMERICOS
N COLS 2 A 30
3RA TARJETA-MY N, EL NMERO DE ECUACUINOES DE COND. Y
OBSERVACIONES RESPECTIVAMENTE. FORMAT0 (215)
MATRIZ W-UN (1) VALOR PARA CADA ECUACION DE CONDICIN
MATRIZ B-COEFICIENTES DE OBSERVACIONWES CON MXIMO DE
CUATRO VALORES POR TARJET. LEEER TODO LOS VALORES PARA
CADA ECUACIN, LUEGO COMENZAR CON UNA TARJETA NUEVA.
FORMATO (4F20.3)
MATRIZ P-PESOS DE OBSERVACIONES CON MXIMO DE CUATRO
VALORES POR TARJETA. LEER SOLO ELEMENTOS DIAGONALES Y
EN SECUENCIA CORRECTA. FORMATO (4F20.3)
IMPLICIT DOUBLE PRECISION (A=H,O-Z)
DIMENSION W(20,1),B(20,50),C(20,20),D(50,20),V(50,1)
DIMENSION WW(20),P(50,50), TITLE(16),PT(50,20)
C
C
C
155
15
16
145
120
18
C
C
C
C
C
44
C
C
C
C
130
20
C
C
C
11
C
C
C
4
2
PRINT 111,(B(I,J)),I=1,M)
DO 18 I=1,N
DO 18 J=1,N
P(I,J)=0.D0
READ 110,(P(J,J)J=1,N)
PRINT 17
PRINT 190
PRINT 111,(P(J,J),J=1,N)
RESOLVER EL ALGORITMO DE MINIMOS CUADRADOS
P INVERSA = P
DO 44 I=1,N
P(I,J)=1.D0/P(I,I)
TRANSPONER B A BT
DO 130 I=1,M
DO 130 J=1,N
MAULTIPLICAR P*BT=PT
DO 20 I=1,N
DO 20 K=1,M
PT(I,K)=0.D0
DO 20 J=1,N
PT(I,K)=PT(I,K)+P(I,J)*BT(J,K)
MULTIPLICADOR B*PT = C
DO 11 I=1,M
DO 11 J=1,M
C(I,J)=0.D0
DO 11 K=1,N
C(I,J)=C(I,J)+B(I,K)+PT(K,J)
INVERTIR C
DO 1 K=1,M
DO 2 J=1,M
IF(J=K)4,2,4
C(K,J)=C(K,J)/C(K,K)
CONTINUE
5
6
3
1
C
C
C
21
C
C
C
31
140
105
110
111
160
170
180
190
150
C(K,K)=1.D0/C(K,K)
DO 1 I=1,M
IF(I=K)5,1,5
DO 3 J=1,M
IF(J=K)6,3,6
C(I,J)=C(I,J)-C(I,K)*C(K,J)
CONTINUE
C(I,K=-C(I,K))*C(K,K)
CONTINUE
MULTIPLICADOR PT * C = D
DO 21 I=1,N
DO 21 J=1,M
D(I,J)=0.D0
DO 21 K=1,M
D(I,J)=D(I,J)+PT(I,K)*C(K,J)
MULTIPLICADAR D*W=V
DO 31 I=1,N
DO 31 J=1,1
V(I,J)=0.D0
DO 31 K=1,M
V(I,K)=V(I,J)+D(I,K)*W(K,J)
PRINT 17
PRINT 180
DO 140 I=1,N
PRINT 111,(V(I,J),J=1,1)
KK=KK+1
IF(KK-N0)115,155,150
FORMAT(215)
FORMAT(4F20.3)
FORMAT(3F20.7)
FORMAT(10X,9H B MATRIX)
FORMAT(10X,9H N MATRIX)
FORMAT(10X,10H RESIDUALS)
FORMAT(10X,9H P MATRIX)
CONTINUE
PRINT 16
STOP
END
EJEMPLO RED DE NIVEL NO PONDERADA, COMPUTOS DE
WOLF
4
.06
1.
-.12
1.
-1
1.
1.
1.
1.
1.
1.
1.
.06
1.
-1.
3.
6.
-8.9
1.
-.02
.05
1.
1.
-1.
1.
1.
1.
1.
.0.5
1.
4.
61.3
1.
1.
1.
1.
2.29
1.34
1.
1.
-3.49
1.
1.
1.
1.
-.08
.38
-6.30
-2.60
3.46
-2.54
1.
1.
1.
1.
1.
1.
METODO CRANDALL PARA EJEMPLO POLIGONAL, COMPUTOS
AJUSTE DE WOLF
2
6
-5.14
1.
-.002
.710
-.763
-.646
.705
.646
1.
-.764
1.34
1.53
1.52
2.89
2.25
POLIGONAL DE LA FIG. 13.2, COMPUTOS AJUSTE
DE WOLF
3 5
-0.000291
-0.205
-0.130
1.00
1.0
1.0
0.0
186.603
86.603
0.0
0.5
-223.205
-50.00
0.0
0.866
4730000.0
4730000.0
4730000.0
156.00
-.912
-.410
2.27
0.0
0.866
0.5
400.0
APENDICE F
PROGRAMA PARA LA TRANSFORMACIN CONFORME DE COORDENADAS
BIDIMENSIONALES
F.1
Introduccin al programa
F.2
El croquis siguiente indica la manera de preparar las tarjetas para la lectura de los
datos de entrada en la computadora:
Fig. F.1
F.4
Salida
Las pginas siguientes contienen 1os listados del programa de computadora y los
datos de entrada para e1 problema del Art. 16.4.
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
C
77
C
C
C
4
C
C
C
PRINT 77
FORMAT(1H1)
READ 32,TITLE
PRINT 32,TITLE
READ 1,NP,NCON
READ 2,(P(I),Q(I),XCON(I),YCON(I),I=1,NCON)
PRINT 20
PRINT 6,(P(I),Q(I),XCON(I),YCON(I),I=1,NCON)
READ 7,(P(I),Q(I),XMEAS(I),YMEAS(I),I=1,NP)
PRINT 21
PRINT 8,(P(I),Q(I),XMEAS(I),YMEAS(I),I=1,NP)
FORMAR LAS MATRICES
DO 3 I=1,NCON
A(I,1)=1.
A(I,2)=XMEAS(I)
A(I,3)= YMEAS(I)
A(I,4)=0.
EL(I,1)=XCON(I)
J=NCON+1
K=2*NCON
DO 4 I=J,K
II=I=NCON
A(I,1)=0.
A(I,2)=YMEAS(II)
A(I,3)=XMEAS(II)
A(I,4)=1.
EL(I,1)=YCON(II)
RESOLVER EL ALGORITMO DE MINIMOS CUADRADOS
5
C
C
C
1
2
6
7
8
20
21
22
30
31
32
C
C
C
C
C
304
302
305
306
303
307
C
C
10
C
C
C
100
DO 307 K=1,N
DO 302 J=1,N
IF(J=K) 304,302,304
QXX(K,F)=QXX(K,J)/QXX(K,K)
CONTINUE
QXX(K,K)=1./QXX(K,K)
DO 307 I=1,N
IF(I=K) 305,307,305
DO 303 J=1,N
IF(J=K) 306,303,306
QXX(I,J)=QXX(I,J)-QXX(I,K)*QXX(K,J)
CONTINUE
QXX(I,K)=-QXX(I,K*QXX(K,K))
CONTINUE
RETURN
END
SUBROUTINE MATR (A,AT,M,N,NMAX,NMAX)
A=MATRIZ ORIGINAL, AT=MATRIZ TRANSPUESTA
M=NMERO DE FILAS EN A, N = NMERO DE COLUMNAS EN A
DOUBLE PRECISION A(MMAX,NMAX),AT(NMAX,MMAX)
DO 10 I=1,M
DO 10 J=1,N
AT(J,I)=A(I,J)
CONTINUE
RETURN
END
SUBROUTINE MATM (A,B,C,M,N,L,MMAX,NMAX,LMAX)
A= PRIMERA MATRIZ DE DIMENMSIONES (M,N)
B= PRIMERA MATRIZ DE DIMENMSIONES (N,L)
C= PRODUCTO DE A VECES B
DOUBLE PRECISION A(MMAX,NMAX),B(NMAX,LMAX),C(MMAX,LMAX)
DO 100 I=1,M
DO 100 H=1,L
C(I,K)=0
DO 100 J=1,N
C(I,K)=C(I,K)+A(I,J)*B(J,K)
CONTINUE
RETURN
END
REFERENCIAS
Barry, B.A., Engineering Measurements, John Wiley & sons. New York, 1964.
Benjamin, J.R. and C.A. Cornell, Probability, Statistics and Decision for Engineers,
McGraw-Hill Book Co., Inc., New York, 1970.
Bjerhammar. A., Theory of Errors and Generalized Matrix Inverses, Elsevier
Scientific Publishing Co., New York, 1973.
Brinker, R.C. and P.R. Half, Elementary Surveying , 6th Ed., Harper and Row
Publishing Co., Inc., New York, 1977.
Cooper, M.A.R. Fundamentals of Survey Measurement and Analysis, Crosby,
Lockwood, Staples, London, 1973.
Crandall, K.C. and R.W. Sea bloom, Engineering Fundamentals in Measurements.
Probability, Statistics and Dimensions. McGraw-Hill Book Co., Inc., New
York, 1970
Deming, W.E. Statistical Adjustment of Data. Dover Publications, Inc., New York,
1943.
Hamilton, W.C., Statistics in Physical Science, the Ronald Press Co., New York,
1964.
Hirvonen, R.A., Adjustment by Least Squares Geodesy and Photogrammetry,
Frederick Unger Publishing Co., New York, 1965.
Maksoudina, Y.L., Probability and Statistics with Applications, International
Textbook Co., Scranton, 1969.
Richardus, P., project Surveying, North Holland Publishing Co., Amsterdam John
Wiley & Sons, Inc.., New York, 1966.
Snedecor, G.W. and W.C. Cochran, Statistical Methods, Iowa State Univ. Press,
Ames, Iowa, 1973.
Wolf, P.R., Elements of Photogrammetry, McGraw-Hill Book Co., Inc., New York,
1974.
Veress, S.A., Adjustment by Least Squares, American Congress on Surveying and
Mapping, Washington, D.C., 1974.
Aguilar, A.M, Principles of Survey Error- Analysis and Adjustment," Journal of the
the surveying and Mapping Division, ASCE, Vol. 99, No. SU1 Sept., 1973.
Baker,J.R., Adjusting Radial Surveys," Journal of the Surveying and Mapping
Division, ASCE, Vol. 98, No. SU2, Nov., 1972.
Fubara, D.M.J., "Three-Dimensional Adjustment of Terrestrial Geodetic Networks,
Canadian Surveyor. Vol. 26, No. 4, Sept., 1972.
Gale, L.A., Theory of Adjustments by Least Squares," Canadian Surveyor, Vol. IX,
No. 1, Mar., 1965.
Konecny, G., Classical Concepts of Least Squares Adjustment," Canadian
Surveyor, Vol. IX, No. 1, Mar., 1965.
Kouba, G Principal Property and Least Squares Adjustment Design," Canadian
Surveyor, Vol. 26, No. 2, June, 1972.
Madkour, M.F., Non-linear Least Squares Algorithm," journal of the Surveying and
Mapping Division, ASCE, Vol. 98. No. SU2, Nov., 1972.
Navasky, M., Adjusting Observed Directions by Angles Using Equitable Wights,
Surveying and Mapping, Vol. XXXIV, No. 4, Dec., 1974.
Pinch, M.C and A.E. Peterson, "A Methods for Adjusting Survey Networks," in
sections, Canadian Surveyor, Vol, 28, No. 1, Mar., 1974.
E. Schmid, "A Generalized Least Squares Solution for hybrid Measuring Systems,"
Canadian Surveyor, Vol. IX, No. 1, Mar., 1965