MÉTODO DE
COMPENSACIÓN POR
MÍNIMOS CUADRADOS
TEMA 3
Método de compensación por
mínimos cuadrados Introducción
Lo aplicaremos cuando realicemos más
observaciones de las necesarias con el objetivo de
aumentar la precisión y seguridad en el resultado.
Es el método de ajuste de observaciones más
adecuado.
Consiste en calcular el valor que satisfaga la
condición de que la suma de los cuadrados de los
residuos sea mínima.
Método de compensación por
mínimos cuadrados Introducción
Clases de sistemas de medición:
Método de compensación por
mínimos cuadrados Medidas directas
Medidas de Medidas de
igual precisión distinta precisión
MEDIA MEDIA
ARITMÉTICA PONDERADA
Método de compensación por
mínimos cuadrados Medidas indirectas
P
Z 4
m
,0
99
N
51 82
g
43g 04
(0,0,0) E
En topografía clásica
obtenemos la distancia
reducida y el desnivel entre
dos puntos a partir de la Dr=Dg*senV
distancia geométrica y el
ángulo cenital entre ellos.
Método de compensación por
mínimos cuadrados Medidas indirectas
Lo normal no será el sencillo ejemplo anterior sino
funciones con más de una incógnita.
Vamos a suponer que tratamos de determinar tres
magnitudes x, y, z a partir de una serie de
observaciones directas.
Las observaciones y las magnitudes a determinar
están relacionadas por una función:
ax by cz k 0
A la que se denomina ecuación de observación.
A, b, c y k son coeficientes conocidos, extraídos directa o indirectamente de cada medición.
Método de compensación por
mínimos cuadrados Medidas indirectas
Para encontrar una solución a cada incógnita
habrá que formar tantas ecuaciones como
datos desconocidos.
a1 x b1 y c1 z k1 0
a2 x b2 y c2 z k2 0
a3 x b3 y c3 z k3 0
Pero así sólo tenemos un resultado de cada incógnita.
No podemos contrastarlo ni valorar su precisión.
Método de compensación por
mínimos cuadrados Medidas indirectas
Para ello establecemos más ecuaciones de
observación independientes.
a1 x b1 y c1 z k1 0
a2 x b2 y c2 z k2 0
...................................
an x bn y cn z kn 0
Ahora el problema está en resolverlo, ya que es sobredeterminado.
Esto es, con más ecuaciones independientes que incógnitas,
por tanto, imposible de encontrar valores de x, y, z que
satisfagan rigurosamente el sistema.
Método de compensación por
mínimos cuadrados Medidas indirectas
Recurriremos a hallar los valores más
probables de las incógnitas.
Método de ajuste
de ecuaciones de observación
Primero damos
a1 x b1 y c1 z k1 r1 Al resolver cada
ecuación, seguramente
un valor aproximado a2 x b2 y c2 z k2 r2 no se cumplirá la
a cada incógnita igualdad, sino que
................................... quedará un residuo
a x b y c z k r distinto de cero.
n n n n n
Método de compensación por
mínimos cuadrados Medidas indirectas
Al ser el método de ajuste elegido el de
mínimos cuadrados, adoptaremos como
valores más probables de x, y, z los que
cumplan: r1 r2 ..... rn Mínimo, o sea:
2 2 2
a1 x b1 y c1z k1 2 a2 x b2 y c2 z k2 2 ... an x bn y cn z kn 2 Mínimo
(Derivamos e igualamos a cero)
x, y, z
2a1 a1 x b1 y c1 z k1 2a2 a2 x b2 y c2 z k2 ... 2an an x bn y cn z kn 0
2b1 a1 x b1 y c1 z k1 2b2 a2 x b2 y c2 z k2 ... 2bn an x bn y cn z kn 0
2b1 a1 x b1 y c1 z k1 2b2 a2 x b2 y c2 z k2 ... 2bn an x bn y cn z kn 0
Método de compensación por
mínimos cuadrados Medidas indirectas
2a1 a1 x b1 y c1 z k1 2a2 a2 x b2 y c2 z k2 ... 2an an x bn y cn z kn 0
2b1 a1 x b1 y c1 z k1 2b2 a2 x b2 y c2 z k2 ... 2bn an x bn y cn z kn 0
2b1 a1 x b1 y c1 z k1 2b2 a2 x b2 y c2 z k2 ... 2bn an x bn y cn z kn 0
El sistema resultante tiene resueltos los
dos problemas iniciales:
Cumple la condición
Es compatible y determinado
de los mínimos cuadrados
Método de compensación por
mínimos cuadrados Medidas
condicionales
En topografía las mediciones se extienden
sobre grandes superficies.
Lo que provoca que las medidas entre
tramos se relacionen entre sí formando
figuras geométricas: g
La suma de los tres ángulos de cada b
f
triángulo son 200g.
Cumplir el teorema de los senos. a c
La suma de los ángulos en a es 400g. e
La suma de los ángulos interiores del
d
polígono es (n-2)200g. ...
Método de compensación por
mínimos cuadrados Medidas condicionales
El objetivo del ajuste es determinar el valor de las correcciones a
aplicar a cada dato de campo para que se cumplan estrictamente
todas las condiciones matemáticas impuestas añadiendo la
condición de mínimos cuadrados.
Para el triángulo ABC: A ˆ Bˆ Cˆ 200 g.00.16
La condición geométrica impone que:
Aˆ x1 Bˆ x2 Cˆ x3 200 g.00.00
Restándolas: x1 x2 x3 0 g.00.16 0
Obteniendo la primera ecuación condicional en función
de los residuos.
Método de compensación por
mínimos cuadrados Medidas
condicionales
Generalizando la expresión tendremos :
a1 x1 a2 x2 a3 x3 ... an xn k1 0 f
Y para cada condición independiente:
f1 a1 x1 a2 x2 a3 x3 ... an xn k1 0
f 2 b1 x1 b2 x2 b3 x3 ... bn xn k2 0 Ecuaciones
condicionales
...............................................................
f n n1 x1 n2 x2 n3 x3 ... nn xn kn 0
Método de compensación por
mínimos cuadrados Medidas condicionales
f1 a1 x1 a2 x2 a3 x3 ... an xn k1 0
f 2 b1 x1 b2 x2 b3 x3 ... bn xn k 2 0 Ecuaciones
............................................................... condicionales
f n n1 x1 n2 x2 n3 x3 ... nn xn k n 0
Es un sistema de ecuaciones lineales e independientes.
El número de ecuaciones es igual al número de redundancias.
Redundancia = datos medidos - datos independientes
ecuaciones datos incógnitas
No existe una única solución. Hay que aplicar adicionalmente el criterio de mínimos cuadrados.
Para ello determinaremos una función que los cumpla.
Método de compensación por
mínimos cuadrados Medidas condicionales
Condición de mm.cc. x12 x22 x32 ... xn2 Mínimo
Se recurre a la Función de Lagrange:
F 21 f1 22 f 2 23 f 3 ... 2n f n
No es más que multiplicar por unos factores a cada
ecuación condicional y añadir su suma a f.
li son coeficientes indeterminados denominados
operadores o multiplicadores de Lagrange, y hay
tantos como ecuaciones de condición.
Método de compensación por
mínimos cuadrados Medidas condicionales
f3 0
Cuando se satisfagan
f1 0 las ecuaciones de fn 0
condición
f2 0
F 21 f1 22 f 2 23 f 3 ... 2n f n F
Al ser equivalentes, cuando una tome valor mínimo
también lo tomará la otra
Fmín mín
Método de compensación por
mínimos cuadrados Medidas condicionales
Haciendo nulas las primeras derivadas parciales de F:
2 x1 21a1 22 a1 ... 2n a1 0
2 x2 21a2 22 a2 ... 2n a2 0 Ecuaciones
..................................................... correlativas
2 xm 21am 22 am ... 2n am 0
Despejando las incógnitas x1, x2,…,xm, quedarán en
función de los multiplicadores de Lagrange.
x1 1a1 2 a1 ... n a1 f1 a1 x1 a2 x2 a3 x3 ... an xn k1 0
f 2 b1 x1 b2 x2 b3 x3 ... bn xn k2 0
x2 1a2 2 a2 ... n a2 Sustituyéndolas
...............................................................
....................................... en las ecuaciones f n n1 x1 n2 x2 n3 x3 ... nn xn k n 0
xm 1am 2 am ... n am condicionales
ECUACIONES NORMALES
Método de compensación por
mínimos cuadrados Medidas condicionales
De igual número de ecuaciones que incógnitas,
pudiendo resolverse para obtener los valores li.
Que sustituimos finalmente en las ecuaciones
correlativas para obtener los valores de las
correcciones x1, x2, …,xm.
Con los que corregimos las medidas efectuadas
en campo y resulta ajustada por mínimos
cuadrados la figura.
Método de compensación por
mínimos cuadrados Medidas condicionales
El complejo cálculo ha
sido imposible hasta la
aparición de los potentes
ordenadores actuales,
teniendo en cuenta la
complejidad de las figuras
que se pueden dar en
Topografía.