A MINIMOS CUADRADOS NILTON LEON
A MINIMOS CUADRADOS NILTON LEON
A MINIMOS CUADRADOS NILTON LEON
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
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.
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
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:
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)
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)
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:
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:
𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈 𝜕𝑈
= = = = = = = =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: