A MINIMOS CUADRADOS NILTON LEON

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 36

TOPOGRAFÍA II

COMPENSACIÓN DE POLIGONALES POR


EL MÉTODO DE MÍNIMOS CUADRADOS

ESTUDIANTE: LEÓN VALVERDE NILTON


MÉTODO DE MÍNIMOS CUADRADOS

 Es un proceso en el cual buscamos es obtener una


función y=f(x) que se ajuste lo mejor posible a los
valores experimentales obtenidos en una
determinada labor.

 Estas funciones pueden ser: rectas, polinomios,


funciones potenciales, etc.
Ejemplo:

X 2 4.2 5.8 8.3 9.6


Y 1 2 3 4 5

Veamos:
 Supongamos que la función que mejor se ajusta es
una función lineal Y +aX+b=0
 Siendo e1, e2,e3,e4 ,… ,en los errores
 Para estas mediciones debemos hacer que
F=e12+ e22+e32+ e4 2+… +en2=mínimo…(Método de
mínimos cuadrados).
 Por otra parte tenemos para cada punto tenemos:
 Y1+aX1+b=e1
 Y2+aX2+b=e2
 Y3+aX3+b=e3
 …..
 ….
 ….
 Yn+aXn+b=en
𝜕𝐹 𝜕𝐹
= =0
𝜕a 𝜕b

 Derivada respecto a “a”


 X1(Y1+ax1+b)+ X2(Y2+ax2+b)+
X3(Y3+ax3+b)+….+Xn(Yn+axn+b)=0…..(1)

 Derivada respecto a “b”


 (Y1+aX1+b)+ (Y2+aX2+b)+(Y3+aX3+b)+….+
(Yn+aXn+b)=0……..(2)

 Entonces tenemos dos ecuaciones con dos incógnitas donde


las incógnitas son “a” y “b”
𝑥𝑥 𝑥 𝑎 −𝑦𝑥
 ( )( )=( −𝑦 )
𝑥 𝑛 𝑏
Siendo:

 XX= X12+ X22+ X32+…+ Xn2=216.33


 X= X1+ X2+ X3+…+ Xn=29.9
 YX= X1Y1+ X2Y2+ X3Y3+ ….+ XnYn=109
 Y= Y1+ Y2+ Y3+…+ Yn=15

 POR LO TANTO PARA NUESTRO PROBLEMA TENEMOS:

216.33 29.9 𝑎 −109


 ( )( )=( ) a=-0.514 b=0.075
29.9 5 𝑏 −15
 Entonces la ecuación lineal que mejor se ajusta es:
Y-0.514X+0.075=0
GENERALIZACIÓN

 **Siguiendo con el problema si suponemos una


función cuadrática: Y+ax2+bx+c=0
 Haciendo el proceso del método de mínimos
cuadrados encontraríamos:
𝑥𝑥𝑥𝑥 𝑥𝑥𝑥 𝑥𝑥 𝑎 −𝑦𝑥𝑥
 ( 𝑥𝑥𝑥 𝑥𝑥 𝑥 ) (𝑏 ) =( −𝑦𝑥 )
𝑥𝑥 𝑥 𝑛 𝑐 −𝑦
 Y así podemos generalizar las ecuación matricial
para un polinomio de grado n.
PROBLEMA EN EL CUAL NOS DAN LAS
COORDENADAS Y LA ALTITUD DE VARIOS
PUNTOS

 En este caso nosotros buscaremos la superficie que


más se ajusta a estos puntos.
 Haciendo Z+aX+bY+c=0
 Entonces al seguir el proceso del método de mínimos
cuadrados obtendremos:
𝑦𝑦 𝑥𝑦 𝑦 𝑎 −𝑦𝑧
 (𝑥𝑦 𝑥𝑥 𝑥 ) (𝑏 ) =(−𝑥𝑧)
𝑦 𝑥 𝑛 𝑐 −𝑧
 De esta ecuación obtendremos los parámetros a,b y
c.
***COMPENSACIÓN DE ÁNGULOS DE
IGUAL PRECISIÓN

 Ejemplo:
En un triángulo teóricamente la suma de sus ángulos
es 180 grados, pero ello no se da en la práctica ya que
hay errores humanos, instrumentales o por refracción.

 a1+a2+a3-180=error de cierre
 Siendo V1,V2,V3 los errores parciales de la medición
de los ángulos del triangulo.

f1=V1+V2+V3-E=0. Ecuación condicional.

 Por condición del método de mínimos cuadrados


F= V12+V22+V32=mínimo
Ecuación de LaGrange

 Usando la ecuación de LaGrange U=F-λ1f1


 Para que F sea mínimo se debe cumplir lo
𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈
siguiente: = = = =0
𝜕V1 𝜕V2 𝜕V3 𝜕λ1
 -------2V1- λ1=0………..(1)
 -------2V1- λ1=0………….(2)
 -------2V1- λ1=0………….(3)
 -------V1+V2+V3-E=0………………(4)
 Resolviendo las ecuaciones:
 V1= λ1/2
V2= λ1/2 V3= λ1/2 entonces
V1=V2=V3=E/3
***COMPENSACIÓN DE ÁNGULOS DE
DIFERENTE PRECISIÓN

 Sea V1, V2 ,V3 los errores accidentales de los ángulos


internos de un triángulo y P1 ,P2 y P3 los pesos de
cada error que se comete al medir los ángulos.

 f1=V1+V2+V3-E=0 ecuación condicional

 Por condición del método de mínimos cuadrados


F=P1 V12+P2V22+P3V32=mínimo
Ecuación de LaGrange

 Usando la ecuación de LaGrange U=F-λ1f1


 Para que F sea mínimo se debe cumplir lo
𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈
siguiente: = = = =0
𝜕V1 𝜕V2 𝜕V3 𝜕λ1
 -------2P1V1- λ1=0………..(1)
 -------2P2V2- λ1=0………….(2)
 -------2P3V3- λ1=0………….(3)
 -------V1+V2+V3-E=0………………(4)
 Resolviendo las ecuaciones:
 V1= λ1/2 P1 V2= λ1/2 P2 V3= λ1/2 P3
reemplazando en:
 V1+V2+V3-E=0………………(4)
 (λ1/2)( 1/ P1+1/ P2+1/P3)=E
 (λ1/2)K0=E
 Por lo tanto .
 V1= E/ K0P1 V2= E / K0P2 V3= E/ K0P3
 De esto concluimos que en el error angular es
inversamente proporcional al peso.
***RED DE NIVELACIÓN

 Ejemplo:
Se muestra una red de nivelación constituida por dos
circuitos cerrados, se tiene como datos el desnivel
entre cada banco de nivel. Se pide realizar los ajustes
respectivos sabiendo que las longitudes son
iguales.
 Por dato del problema las longitudes son iguales por
lo tanto el peso es la misma en todos los tramos.
 Siendo V1,V2,V3,V4 y V5 los errores de los tramos
AB,BC,CA,DB y CD respectivamente.
 Ecuaciones condicionales:
 Circuito I(ABCA):(+12,179+V1)+(-7.324+V2)+ (-4.870 +V3)=0

f1 =V1+ V2+ V3-0.015=0


 Circuito II(BCDB): (-7.324+V2)+(+15.172+V5)+ (-7.843+V4)=0

f2 =V2+ V4+ V5+0.005=0


Multiplicadores de LaGrange

 Debemos encontrar V1, V2, V3,V4 y V5 tal que


F=V12+V22+ V32+V4 2+V52=mínimo

 U=F-λ1f1-λ2f2 entonces:
𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈
 = = = = = = =0
𝜕V1 𝜕V2 𝜕V3 𝜕v4 𝜕v5 𝜕λ1 𝜕λ2
 2 V1- λ1=0…………(1)
 2 V2- λ1- λ2=0……..(2)
 2 V3- λ1=0……..(3)
 2 V4- λ2=0……..(4)
 2 V5- λ2=0……..(5)
 Entonces solucionando las ecuaciones:
 V1=λ1/2 V2=( λ1+λ2)/2 V3=λ1/2 V4=λ2/2 V5= λ2/2
 Sustituyendo V1, V2, V3,V4 y V5 en la ecuaciones
condicionales.
 En f1: (3λ1+ λ2)/2 =0.015
 En f2: (λ1+ 3λ2)/2 =-0.005
3 1 λ1 0.03
 ( )( )=( )
1 3 λ2 −0.01
 λ1=0.0125 λ2=-0.0075 …. por lo tanto: V1=0.006
V2=0.003 V3=0.006 V4=-0.004 V5= -0.004
UTILIZANDO DISTANCIAS DIFERENTES
DE LA MEDICIÓN:

 Se puede deducir que para una mayor longitud se


tendrá un menor peso.
 Entonces: P1=1/L1 , P2=1/L2, P3=1/L3 son los pesos
de cada tramo del triángulo.

 Por lo tanto en general tenemos

 F=(1/L1 )V12 +(1/L2 )V22 +(1/L3 )V32 +……+(1/Ln )Vn2 =MÍNIMO


COMPENSACIÓN DEL CUADRILÁTERO

1
2
8
3

4
5
7
6
CONDICIÓN DE ÁNGULOS
 SUMA DE ÁNGULOS OPUESTOS POR EL VÉRTICE
1+2=6+5+e1……………..(1)
3+4=7+8+e2……………(2)

 En los triángulos:
ABC: 1+2+3+4=180+e3…..(3)
ADC:5+6+7+8=180+e4 …(4)
ADB:1+2+7+8=180+e5…(5)
BCD:3+4+5+6=180+e6..(6)

 SUMA TOTAL DE ÁNGULOS


1+2+3+4+5+6+7+8=360+e7…..(8)
 SUMA DE ANGULOS OPUESTOS POR EL VÉRTICE
f1:v1+v2-v5-v6-e1=0……………..(1)
f2:v3+v4-v7-v8-e2=0……………(2)

 En los triángulos:
ABC: f3:v1+v2+v3+v4-e3=0…..(3)
ADC: f4:v5+v6+v7+v8-e4 =0…(4)
ADB: f5:v1+v2+v7+v8-e5=0…(5)
BCD: f6:v3+v4+v5+v6-e6=0..(6)

 SUMA TOTAL DE ÁNGULOS


f7:v1+v2+v3+v4+v5+v6+v7+v8-e7=0…..(7)
POR CONDICION DEL LADO.

Un cuadrilátero puede quedar determinado


completamente por exactamente cinco de los seis lados
(sus cuatro lados y sus dos diagonales). Por lo que el
último lado resultaría redundante, que puede ser
utilizado como comprobación. Esta comprobación es
expresada mediante una ecuación de condición.
sen 1 sen 3 sen 5 sen 7 = sen 2 sen 4 sen 6 sen 8
sen(1 + 𝑣1) sen(3 + 𝑣3) (sen 5 + 𝑣5) sen(7 + 𝑣7)
= sen(2 + 𝑣2) sen(4 + 𝑣4) sen(6 + 𝑣6) sen( 8 + 𝑣8)

log sen 1 + 𝑣1 + 𝑙𝑜𝑔 sen(3 + 𝑣3) + 𝑙𝑜𝑔 (sen 5 + 𝑣5) + 𝑙𝑜𝑔 sen(7 + 𝑣7)
− 𝑙𝑜𝑔 sen 2 + 𝑣2 − 𝑙𝑜𝑔 sen 4 + 𝑣4 − 𝑙𝑜𝑔 sen 6 + 𝑣6
− 𝑙𝑜𝑔 sen( 8 + 𝑣8) = 0 … … … … … … … … … … … … … … … … … … … . (8)
 Como los errores son pequeños entonces de
considera que:

 log sen 1 + 𝑣1 = 𝑙𝑜𝑔𝑠𝑒𝑛1 + 𝑉1δ1


 log sen 1 + 𝑣2 = 𝑙𝑜𝑔𝑠𝑒𝑛2 + 𝑉2δ2
 log sen 1 + 𝑣3 = 𝑙𝑜𝑔𝑠𝑒𝑛3 + 𝑉3δ3
 …..
 …..
 log sen 8 + 𝑣8 = 𝑙𝑜𝑔𝑠𝑒𝑛8 + 𝑉8δ8
EJEMPLO:

 log sen 1 + 𝑣1 = 𝑙𝑜𝑔𝑠𝑒𝑛1 + 𝑉1 δ1


Supongamos que 1=60grados y v1=5´´

δ1= log sen 60 + 1/3600 − log sen 60 =


-0.06246815269-(-0.0624693683)=0.000001215

log sen 60 + 5/3600 = 𝑙𝑜𝑔𝑠𝑒𝑛60 + 5δ1


-0.06246329036=-0.0624693683+5x0.000001215
-0.06246329036= -0.0624632933
 REEMPLAZANDO:

 𝑙𝑜𝑔𝑠𝑒𝑛1 + 𝑉1δ1+𝑙𝑜𝑔𝑠𝑒𝑛3 + 𝑉3δ3+𝑙𝑜𝑔𝑠𝑒𝑛5 + 𝑉5δ5


+𝑙𝑜𝑔𝑠𝑒𝑛7 + 𝑉7δ7-𝑙𝑜𝑔𝑠𝑒𝑛2 − 𝑉2δ2-𝑙𝑜𝑔𝑠𝑒𝑛4 − 𝑉4δ4-
𝑙𝑜𝑔𝑠𝑒𝑛6 + 𝑉6δ6-𝑙𝑜𝑔𝑠𝑒𝑛 − 𝑉8δ8

f8: log(sen1)+log(sen3)+log(sen5)+log(sen7)-log(sen2)
-log(sen4)-log(sen6)-log(sen8)+ 𝑉1 δ1+𝑉3 δ3+𝑉5δ5+𝑉7
δ7−𝑉2δ2-𝑉4δ4−𝑉6 δ6−𝑉8δ8=0………(8)
 Entonces debemos encontrar los errores
V1,V2,V3,V4,V5,V6,V7 yV8 haciendo que:

 F=V12+ V22+ V32+ V42+V52 +V62 +V72 +V82=mínimo


 POR LA ECUACIÓN DE LaGrange
U=F-λ1f1-λ2f2-λ3f3-λ4f4-λ5f5-λ6f6-λ7f7-λ8f8

𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈
= = = = = = = =0
𝜕v1 𝜕v2 𝜕v3 𝜕v4 𝜕V5 𝜕v6 𝜕v7 𝜕v8
𝜕𝑈

𝜕v1
=2v1-λ1(1)-λ2(0)-λ3(1)-λ4(0)-λ5(1)-λ6(0)-λ7(1)- λ8(δ1)=0
𝜕𝑈
 =2v2-λ1(1)- λ2(0)-λ3(1)-λ4(0)-λ5(1)-λ6(0)-λ7(1)- λ8(-δ2)=0
𝜕v2
𝜕𝑈
 =2v3-λ1(0)-λ2(1)-λ3(1)-λ4(0)-λ5(0)-λ6(1)-λ7(1)- λ8(δ3)=0
𝜕v3
𝜕𝑈
 =2v4-λ1(0)-λ2(1)-λ3(1)-λ4(1)-λ5(0)-λ6(1)-λ7(1)- λ8(-δ4)=0
𝜕v4
𝜕𝑈
 =2v5-λ1(-1)-λ2(0)-λ3(0)-λ4(1)-λ5(0)-λ6(1)-λ7(1)- λ8(δ5)=0
𝜕v5
𝜕𝑈
 =2v6-λ1(-1)-λ2(0)-λ3(0)-λ4(1)-λ5(0)-λ6(1)-λ7(1)-λ8(-δ6)=0
𝜕v6
𝜕𝑈
 =2v7-λ1(0)-λ2(-1)-λ3(0)-λ4(1)-λ5(1)-λ6(0)-λ7(1)- λ8(δ7)=0
𝜕v7
𝜕𝑈
 =2v8-λ1(0)-λ2(-1)-λ3(0)-λ4(1)-λ5(1)-λ6(0)-λ7(1)- λ8(-δ8)=0
𝜕v8
 2v1=λ1+λ3+λ5+λ7+ λ8(δ1)
 2v2=λ1+λ3+λ5+λ7-λ8(δ2)
 2v3=λ2+λ3+λ6+λ7+λ8(δ3)
 2v4= λ2+λ3+λ6+λ7-λ8(δ4)
 2v5=-λ1+λ4+λ6+λ7+λ8(δ5)
 2v6=-λ1+λ4+λ6+λ7-λ8(δ6)
 2v7=- λ2+λ4+λ5+λ7+λ8(δ7)
 2v8=- λ2+λ4+λ5+λ7-λ8(δ8)
 En (1): v1+v2-v7-v8-e1=0
2λ1+0λ2+λ3- λ4+λ5-λ6+0λ7+λ8(δ1-δ2-δ5+δ6)/2-e1=0
 En (2): v3+v4-v7-v8-e2=0
0λ1+2λ2+λ3-λ4-λ5+λ6+0λ7+λ8(δ3-δ4–δ7+δ8)/2-e2=0
 En (3): v1+v2+v3+v4-e3=0
λ1+λ2+2λ3+0λ4+λ5+λ6+2λ7+λ8(δ1-δ2+δ3-δ4)/2-e3=0
 En (4): v1+v2+v5+v6-e1=0
-λ1-λ2+0λ3+2λ4+λ5+λ6+2λ7+λ8(δ5-δ6+δ7-δ8)/2-e4=0
 En (5): v1+v2+v7+v8-e5=0
λ1-λ2+λ3+λ4+2λ5+0λ6+2λ7+λ8(δ1-δ2+δ7-δ8)/2-
e5=0
 En (6): v3+v4+v5+v6-e6=0
-λ1+λ2+λ3+λ4+0λ5+2λ6+2λ7+λ8(δ3-δ4+δ5-δ6)/2-
e6=0
 En (7): v1+v2+v3+v4+ v5+v6+v7+v8-e7=0
0λ1+0λ2+2λ3+2λ4+2λ5+2λ6+4λ7+λ8(δ1-δ2+δ3-δ4
+δ5 –δ6+δ7-δ8)/2-e7=0
 En (8):
log 𝑠𝑒𝑛 𝑖𝑚𝑝𝑎𝑟𝑒𝑠 − log 𝑠𝑒𝑛 𝑝𝑎𝑟𝑒𝑠 +
λ1(δ1-δ2-δ5+δ6)/2+λ2(δ3-δ4-δ7+δ8)/2+
λ3(δ1-δ2+δ3-δ4)/2+λ4(δ5-δ6+δ7-δ8)/2+
λ5(δ1-δ2+δ7-δ8)/2+λ6(δ3-δ4+δ5-δ6)/2+
λ7(δ1-δ2+δ3-δ4+δ5-δ6+δ7-δ8)/2
λ8(δ12+δ22+δ32+δ42+δ52+δ62+δ72+δ82)/2
 Y vemos que tenemos 8 ecuaciones con 8
incógnitas(λ1, λ2,λ3,λ4,λ5,λ6,λ7 y λ8).
 La ecuación matricial es A x λ=e
A:
 λ: e:

 Al resolver la ecuación matricial obtenemos los valores de λ1,


λ2,λ3,λ4,λ5,λ6,λ7 y λ8 por ende encontramos los errores
V1,V2,V3,V4,V5,V6,V7 yV8 ya que estos dependen de los lamda.
GRACIAS POR SU
ATENCIÓN!!!

También podría gustarte

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy