12.07.2015 Views

Teoría de errores - ETSI Topografia (UPM)

Teoría de errores - ETSI Topografia (UPM)

Teoría de errores - ETSI Topografia (UPM)

SHOW MORE
SHOW LESS
  • No tags were found...

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Según las causas que los provocan Errores graves o groseros Son <strong>de</strong>bidos a imperfeccionesen los sentidos <strong>de</strong>l observadoro a distracciones o ten<strong>de</strong>nciasnerviosas en el momento <strong>de</strong> laobservación Errores instrumentales Debidos al equipo utilizado Errores teóricos Debidos a <strong>de</strong>ficienciasen la teoría en que se basanlas observacionesEsta clasificación no se prestaconvenientemente a un estudiomatemático <strong>de</strong> los <strong>errores</strong> Según los efectos que producen Errores constantes Producen el mismo efectosobre todas las observaciones<strong>de</strong> una serie <strong>de</strong> medidas Ej.: error en la orientación<strong>de</strong> un instrumento Errores sistemáticos Son aquellos que siempre que serealice una observación <strong>de</strong> unamisma magnitud en idénticascondiciones se presentan con elmismo valor y en el mismo sentido Errores aleatorios Influyen <strong>de</strong> manera irregulary varían en magnitud y signo<strong>de</strong> una observación a otra


-Tipos:-<strong>errores</strong> <strong>de</strong> escritura-<strong>errores</strong> <strong>de</strong> lectura-interpretaciones erróneas <strong>de</strong> las lecturas-Etc.-Interpretación estadística: Las observaciones con<strong>errores</strong> graves no pertenecen a la muestraaleatoria( y por tanto no siguen distribuciones <strong>de</strong>tipo Normal)


-Detección y eliminación <strong>de</strong> estos <strong>errores</strong> antes <strong>de</strong>l procesado, pueslos Métodos <strong>de</strong> Ajuste basado en los Mínimos Cuadrados no soncapaces <strong>de</strong> <strong>de</strong>tectarlos.-Algunos métodos:-realizar múltiples lecturas-comprobar su consistencia-utilizar técnicas distintas-incrementar el número <strong>de</strong> observaciones-En general se utilizan tests estadísticos sobre los datos para<strong>de</strong>tectarlos( por ejemplo el Test <strong>de</strong> Baarda, muy utilizado en Topografíay Fotogrametría) Veamos un ejemplo <strong>de</strong> una posible metodología para<strong>de</strong>tectarlos(ejercicio propuesto)


!"-Su comportamiento viene <strong>de</strong>terminado por algún sistema“<strong>de</strong>terminístico” cuya expresión funcional es conocida.-Métodos a utilizar:-calibrado <strong>de</strong> instrumentos-procedimientos <strong>de</strong> observación-procesado <strong>de</strong> los datos-Interpretación estadística: Afectan a todas las observaciones<strong>de</strong>l mismo modo.Parámetro <strong>de</strong> localización <strong>de</strong> la distribución: µ (media)-Otra causa <strong>de</strong> su existencia: Elección incorrecta <strong>de</strong>l mo<strong>de</strong>lofuncional Inconsistencia con las observaciones


Valor y signo permanecen constantes: Error sistemático constante- Ejemplo: ¿Errores sistemáticos en las Fotocoor<strong>de</strong>nadas en unfotograma aéreo?-Origen <strong>de</strong> los <strong>errores</strong> s. en Topografía, Fotogrametría, Geo<strong>de</strong>sia:-Causas físicas-Factores instrumentales-Limitaciones <strong>de</strong> la observación humana-Causas físicas: Temperatura, humedad y cambios <strong>de</strong> presión observaciones angulares y <strong>de</strong> distancia-Factores instrumentales: imperfecciones <strong>de</strong> construcción o<strong>de</strong>sajustes en su utilización


-Tratamiento <strong>de</strong> los <strong>errores</strong> sistemáticos en el mo<strong>de</strong>lo funcional:Aplicación <strong>de</strong> sus efectos(si se conocen) sobre las observaciones apriori <strong>de</strong>l ajusteSe consi<strong>de</strong>ran variables <strong>de</strong>sconocidas durante el proceso <strong>de</strong> ajusteSe consi<strong>de</strong>ran como observaciones, con un valor a priori.Ejemplos:1)Medida clásica <strong>de</strong> distancias:-diferencia <strong>de</strong> temperatura entre el momento <strong>de</strong> la calibración y el <strong>de</strong>observación-tensión aplicada en la cinta-extremos a diferentes alturas-corrección por catenaria


2)Técnicas electrónicas o electro-ópticas:-cambios en la <strong>de</strong>nsidad <strong>de</strong>l aire variación <strong>de</strong> la frecuencia-centrado imperfecto <strong>de</strong>l instrumento-la señal no sigue exactamente una línea recta en supropagación (factores ambientales) Precisión = + a, b : Parámetros <strong>de</strong>l distanciómetrod la distancia medida3)Observación angulares (horizontales y verticales):-centrado imperfecto <strong>de</strong>l círculo horizontal-graduaciones no uniformes en los círculos-eje horizontal <strong>de</strong>l telescopio ⊥ eje vertical <strong>de</strong>l instrumento-eje <strong>de</strong>l telescopio || eje <strong>de</strong>l nivel <strong>de</strong> burbuja


Métodos <strong>de</strong> reducción <strong>de</strong>l efecto <strong>de</strong> los <strong>errores</strong> sistemáticossobre las observacionesEstablecer la ley <strong>de</strong> aparición <strong>de</strong>l error y corregir lasobservacionesEjemplo:Una vez realizada la calibración <strong>de</strong> un distanciómetro utilizandoun patrón <strong>de</strong> medida <strong>de</strong> distancias, se obtiene una expresión <strong>de</strong>ltipo: Corrección = F(Longitud, temperatura, presión,etc)Establecer un programa <strong>de</strong> observaciones <strong>de</strong> tal forma que los<strong>errores</strong> no actúen <strong>de</strong> forma unilateralEjemplos:-En la observación <strong>de</strong> direcciones acimutales : lecturasdiametralmente opuestas para eliminar el efecto <strong>de</strong> laexcentricidad.


Observaciones geodésicas: Realizadas en horas diurnas ynocturnas durante un periodo <strong>de</strong> tiempo prolongado , con elobjetivo <strong>de</strong> <strong>de</strong>bilitar la influencia <strong>de</strong> la refracción.Re<strong>de</strong>s <strong>de</strong> nivelación : Cálculo <strong>de</strong> los <strong>de</strong>sniveles directos einversos para minimizar el efecto <strong>de</strong> los <strong>errores</strong> sistemáticosproducidos por una <strong>de</strong>terminación no <strong>de</strong>masiado precisa <strong>de</strong>lCoeficiente <strong>de</strong> Refracción.Metodología especial en la elaboración o tratamiento <strong>de</strong> losdatos:Ejemplo histórico : En el siglo XVII y bajo la dirección <strong>de</strong>lastrónomo Cassini comienza la construcción <strong>de</strong>l ObservatorioAstronómico <strong>de</strong> París. Era necesario <strong>de</strong>terminar el meridiano <strong>de</strong>París y por ello la latitud : ϕ PARIS


Observaciones realizadas : Distancias cenitales <strong>de</strong> un conjunto <strong>de</strong>estrellas a su paso por el meridiano(z*), y <strong>de</strong> este modo calcular lalatitud utilizando la fórmula:ϕ = δ* - z*don<strong>de</strong> δ* = <strong>de</strong>clinación <strong>de</strong> la estrella(buscadas en el catálogo <strong>de</strong>estrellas <strong>de</strong>l año)


DATOS Z ϕ 10 0 48 0 50’ 11”8 63 0 13”5 1 0 10”9 74 0 14”3 40 022 0 13”211”9 31 0 12”1 53 0 13”3 Dispersión : ¡¡¡¡¡¡ (10”9-14”3) !!!!!!! ¿Podría ocurrir que existiera un error sistemático en todas las z <strong>de</strong>forma que: error sistemático = f(z) ?Efectivamente, se comprobó que:corrección(z) = -(3”) sen z Una vez aplicadas las correcciones a todas las distancias cenitales: (10”6 , 11”4) dispersión < 1”


# *.,.*./. La teoría <strong>de</strong> <strong>errores</strong> clásica fue <strong>de</strong>sarrollada por Gauss (1777 – 1855) Se refiere única y exclusivamente aa los <strong>errores</strong> acci<strong>de</strong>ntales Las observaciones utilizadas se suponenlibres <strong>de</strong> <strong>errores</strong> constantes y sistemáticos Esta teoría parte <strong>de</strong> la siguiente premisa:Se consi<strong>de</strong>ra una serie <strong>de</strong> gran número <strong>de</strong> observaciones,efectuadas en idénticas condiciones y con el mismo esmeroY está basada en los siguientes postulados:1 Los distintos <strong>errores</strong> posibles se presentan con tanta mayorfrecuencia cuanto menor sea su valor absoluto2 Los <strong>errores</strong> <strong>de</strong>l mismo valor absoluto pero <strong>de</strong> distinto signose presentan con igual frecuenciani=1εi=0, n→ ∞#$%$&%$'() '()* +, -


# *.,.*./. Sin embargo, en la primera exposición <strong>de</strong> su teoría(en la obra “Theoria Motus Corporum Celestium” )no aparecía el segundo postulado como se ha enunciadoAparecía con el siguiente enunciado:2’ Obtenidos por observación n valores distintos <strong>de</strong> una magnitud,siempre que todos hayan sido obtenidos en las mismas circunstanciasy con el mismo esmero, el valor más probable <strong>de</strong> la magnitu<strong>de</strong>s la media aritmética <strong>de</strong> los valores observados Este postulado es consecuencia <strong>de</strong>l segundo, puesto quesi x 1 , . . . x n son los valores observadoscada observación vendrá afectada <strong>de</strong> un cierto error i<strong>de</strong> forma que si llamamos x al verda<strong>de</strong>ro valor resulta: i = x – x i , i = 1, . . . n#$%$&%$'() '()* +, 0


# *.La <strong>de</strong>finición dada para el error verda<strong>de</strong>ro presenta una dificultadinsuperable para su cálculo:No se conoce el valor verda<strong>de</strong>ro <strong>de</strong> la magnitudSe introduce, por ello, el residuo que es la diferencia entre el valormás probable(lo que llamaremos valor ajustado) y el valorobservado, es <strong>de</strong>cirv i = x – x in·ii = 1 i = 1x − xini = 1, . . . , nEn el caso <strong>de</strong> medidas repetidas <strong>de</strong> una misma magnitud este valormás probable coinci<strong>de</strong> con la media aritmética.Entonces:Y sumando los residuos resultani = 1vi=n( x − xi) = n·x − x = 0 i = 1nnvi= 0#$%$&%$'() '()* +, 1


# *.Se <strong>de</strong>muestra que partiendo <strong>de</strong> una muestra aleatoria simplesuficientemente gran<strong>de</strong>x 1 , x 2 , . . . x nY consi<strong>de</strong>rando los residuosLa función f toma la formaf( v )=σque correspon<strong>de</strong> a una distribuciónNormal(0,σ)2v1 −2σ· e2π2v i = x – x iSu gráfica es <strong>de</strong> la formaSiendo más puntiaguda o achatadasegún el valor <strong>de</strong> σ#$%$&%$'() '()* +, &


2 .(.Primer Estudio <strong>de</strong> los <strong>errores</strong> aleatorios en medida directa(sinconsi<strong>de</strong>rar la Distribución Normal) Sean X 1 ,X 2 ,...,X n las observaciones y e 1 ,e 2 ,...,e n los <strong>errores</strong>aleatorios presentes en dichas observacionesVeamos diferentes formas <strong>de</strong> aproximar estos <strong>errores</strong>:Error probable: Or<strong>de</strong>nados los <strong>errores</strong> <strong>de</strong> menor a mayormagnitud(sin signo), ocupa la posición central|e1| < |e2|


-Error Medio Aritmético:Media aritmética <strong>de</strong> todos los <strong>errores</strong>(sinsigno)| ei|i=1e.m.a. = ≈nnError Máximo: Establece una tolerancia para el valor <strong>de</strong> los<strong>errores</strong>e i ≤ e maxni=1| v |n∀ i=1...nie max ≈ v maxError Medio Cuadrático: e c2eie c = 1 ≈nnn1v2in −1(Se pue<strong>de</strong> <strong>de</strong>mostrar: Corrección <strong>de</strong> Bessel)


• Error Medio Cuadrático <strong>de</strong> la Media Aritmética:Sean X 1...X nlas observaciones y su mediaxm 2 2 2c=1+2+ +ne = 1 e + e + ... + enm 2 2 2c 1 2 nPero no conocemos el valor exacto <strong>de</strong> los <strong>errores</strong>, pero po<strong>de</strong>mosacotarlos:xe i ≤ e cEntonces:e 1 e e ... ene = 1 e + e + ... + enm 2 2 2c c c cmecn1 2e 1 2cn nec1 n= =e¿cuántas observacionesc= será recomendable hacer?n


Segundo Estudio <strong>de</strong> los <strong>errores</strong> aleatorios en medidadirecta(consi<strong>de</strong>rando la Distribución Normal) Ahora estudiaremos estos conceptos <strong>de</strong>s<strong>de</strong> la perspectiva <strong>de</strong> la D.Normal1 f(x) =σ 2πSe <strong>de</strong>nomina Módulo <strong>de</strong> precisión a la expresión:e2x−2 σ2h= y por tanto quedará:σ12f(x)=2 2 h e−h xπ


Error medio aritmético: e a∞ e a ≈ E[|X|] = ≈ =−∞| x | f(x)dx1h πσ2π e a = 0.797σError probable: ep. Es aquél que cumple lo siguiente: P[|X| < e p ] = 0.5Demostremos que:e p≈ 0.68σ


Error Medio Cuadrático: e c∞e c2 ≈ E[X 2 21] = 2 x f(x)dx ≈ 2 24h01 1ec≈h 212σ 2Error Máximo: e maxemax ≈ 4e p ≈ 4(0.68σ) ≈ 2.5σemax ≈ 2.5e c= e c ≈ σ


Las expresiones más utilizadas son finalmente:e p≈ 0.68σ ≈e a ≈ 0.797σ ≈2345e ce c

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!