03.06.2013 Views

Geodesia. Cartografía. Sistemas de referencia. Tiempos.

Geodesia. Cartografía. Sistemas de referencia. Tiempos.

Geodesia. Cartografía. Sistemas de referencia. Tiempos.

SHOW MORE
SHOW LESS

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

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

Sistema <strong>de</strong> navegación autónomo: Navegación inercial.<br />

Errores en navegación inercial.<br />

Mo<strong>de</strong>los <strong>de</strong> Error<br />

Ruido blanco.<br />

Breve recordatorio <strong>de</strong> estadística<br />

Procesos estocásticos. Ruido blanco. Propagación.<br />

Medidas <strong>de</strong>l error.<br />

Ruido blanco: Se <strong>de</strong>fine como ruido blanco un proceso ν(t)<br />

que verifica:<br />

E[ν(t)] = 0.<br />

E[ν(t)ν(t) T ] = σ 2 Id.<br />

R(t, τ) = E[ν(t)ν(τ) T ] = δ(t − τ)σ 2 Id, don<strong>de</strong> δ(x) vale 1 si<br />

x = 0 y 0 en cualquier otro caso.<br />

La última condición quiere <strong>de</strong>cir que el valor <strong>de</strong>l ruido blanco<br />

en un instante es in<strong>de</strong>pendiente <strong>de</strong> su valor en cualquier<br />

instante anterior.<br />

Ruido blanco gaussiano: Es un proceso que cumple las<br />

condiciones anteriores, y a<strong>de</strong>más es gaussiano.<br />

Un buen mo<strong>de</strong>lo para las fuentes <strong>de</strong> error (errores <strong>de</strong> medida,<br />

errores gravitatorios) es δɛ(tk) = b + Dν, don<strong>de</strong> ν es ruido<br />

blanco gaussiano. El valor <strong>de</strong> b dará la media <strong>de</strong>l error (sesgo,<br />

llamado bias en inglés).<br />

Sistema <strong>de</strong> navegación autónomo: Navegación inercial.<br />

Errores en navegación inercial.<br />

Mo<strong>de</strong>los <strong>de</strong> Error<br />

Propagación <strong>de</strong>l error.<br />

Breve recordatorio <strong>de</strong> estadística<br />

Procesos estocásticos. Ruido blanco. Propagación.<br />

Medidas <strong>de</strong>l error.<br />

45 / 49<br />

Si en las ecuaciones <strong>de</strong> propagación <strong>de</strong>l error<br />

δx(tk+1) = A(k)δx(tk) + δɛ(tk) sustituimos δɛ(tk) = b + Dν,<br />

obtenemos el siguiente mo<strong>de</strong>lo <strong>de</strong> propagación <strong>de</strong>l error:<br />

δx(tk+1) = A(k)δx(tk) + b + Dν.<br />

Observación: típicamente b también está sometido a un error<br />

variable, <strong>de</strong> forma que b(tk+1) = b(tk) + Dbν b. Para<br />

simplificar ignoramos esta variación.<br />

Se realizan las siguientes hipótesis:<br />

ν es ruido blanco gaussiano con varianza σ2 ν.<br />

Inicialmente, δx(t0) ∼ Nn(m0 , Σ0). Si se conocieran<br />

perfectamente, entonces Σ0 = 0.<br />

A<strong>de</strong>más se tiene la hipótesis <strong>de</strong> que δx(t0) y ν son<br />

in<strong>de</strong>pendientes.<br />

Bajo estas condiciones, se tiene que δx(tk) es un proceso<br />

gaussiano, es <strong>de</strong>cir, δx(tk) ∼ Nn(mk, Σk), don<strong>de</strong> la media y la<br />

covarianza verifican la siguiente evolución:<br />

Propagación <strong>de</strong> la media: mk+1 = Amk + b.<br />

Propagación <strong>de</strong> la covarianza: Σk+1 = AΣkAT + σ2 νDDT . 46 / 49<br />

Sistema <strong>de</strong> navegación autónomo: Navegación inercial.<br />

Errores en navegación inercial.<br />

Mo<strong>de</strong>los <strong>de</strong> Error<br />

Propagación <strong>de</strong>l error: ejemplo sencillo<br />

Breve recordatorio <strong>de</strong> estadística<br />

Procesos estocásticos. Ruido blanco. Propagación.<br />

Medidas <strong>de</strong>l error.<br />

Supongamos que tuviéramos una ecuación <strong>de</strong> propagación <strong>de</strong>l<br />

error en una dimensión (por ejemplo la posición en el eje x)<br />

dada simplemente por: δxk+1 = δxk + ν, don<strong>de</strong>:<br />

La variable temporal k representa minutos, es <strong>de</strong>cir, x6 es el<br />

error en posición pasados 6 minutos.<br />

ν es ruido blanco gaussiano <strong>de</strong> varianza σ 2 ν.<br />

Inicialmente, δx(t0) = 0.<br />

A<strong>de</strong>más δx(tk) y ν son in<strong>de</strong>pendientes.<br />

Entonces aplicando las fórmulas anteriores,<br />

δx(tk) ∼ N(mk, σk), don<strong>de</strong> la media y la varianza verifican:<br />

Propagación <strong>de</strong> la media: mk+1 = amk. Como m0 = 0, se<br />

tendrá mk = 0 para todo k.<br />

Propagación <strong>de</strong> la varianza: σ2 k+1 = σ2 k + σ2 nu. Como σ2 0 = 0, se<br />

tiene que σ2 k = kσ2 nu. Por tanto la varianza verifica σk = √ kσν.<br />

Si por ejemplo x son metros y el ruido blanco tiene σν = 0,1<br />

metros, entonces aunque inicialmente la posición se conoce<br />

sin error, pasada una hora σ60 = √ 60 · 0,1 = 0,77, es <strong>de</strong>cir un<br />

intervalo 2-σ sería δx ∈ [−1,55, 1,55]. 47 / 49<br />

Sistema <strong>de</strong> navegación autónomo: Navegación inercial.<br />

Errores en navegación inercial.<br />

Mo<strong>de</strong>los <strong>de</strong> Error<br />

Medida <strong>de</strong>l error en 2-D<br />

Breve recordatorio <strong>de</strong> estadística<br />

Procesos estocásticos. Ruido blanco. Propagación.<br />

Medidas <strong>de</strong>l error.<br />

Para el caso 2-D (por ejemplo posición sobre un mapa) y si el<br />

error está distribuido como X ∼ N2(0, Σ), las regiones <strong>de</strong><br />

confianza serían elipses:<br />

Dado Σ po<strong>de</strong>mos escribir Σ = Pdiag{σ1, σ2}P T don<strong>de</strong> P es<br />

una matriz con autovectores y σi los autovalores. Los<br />

autovectores dan la dirección <strong>de</strong> los ejes <strong>de</strong> la elipse y los<br />

autovalores son proporcionales a su magnitud (cuyo valor<br />

exacto <strong>de</strong>pen<strong>de</strong>rá <strong>de</strong>l grado <strong>de</strong> confianza <strong>de</strong>l intervalo).<br />

48 / 49

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

Saved successfully!

Ooh no, something went wrong!