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

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

03.06.2013 Views

Sistema de navegación autónomo: Navegación inercial. Errores en navegación inercial. Modelos de Error Breve recordatorio de estadística Procesos estocásticos. Ruido blanco. Propagación. Medidas del error. Variables aleatorias continuas multidimensionales Sea una variable aleatoria X ∈ Rn continua multidimensional. Cada componente de X sigue una distribución unidimensional. Como en el caso unidimensional, se define una función de distribución conjunta, que se calcula mediante la función de densidad f (x). Igualmente E[g(X )] = Rn g(y)f (y)dy. Los dos casos importantes son: Media: m(X ) = E[X ] = Rn yf (y)dy. Matriz de covarianzas: Cov(X ) = E[(X − m(X ))(X − m(X )) T ] = Σ. Es una matriz simétrica y definida positiva. Los valores de la diagonal representan la varianza de cada componente de X , mientras que los valores fuera de la diagonal la correlación entre dos componentes de X . Se tiene Σ = E[(X X T ] − m(X )m(X ) T . Por ejemplo, para n = 3 y escribiendo X = [X , Y , Z]: 2 6 Σ = 4 σ 2 x E[(X − mx )(Y − my )] E[(X − mx )(Z − mz )] E[(X − mx )(Y − my )] σ 2 y E[(Y − my )(Z − mz )] E[(X − mx )(Z − mz )] E[(Y − my )(Z − mz )] σ 2 z Sistema de navegación autónomo: Navegación inercial. Errores en navegación inercial. Modelos de Error Distribución normal multivariante I 3 7 5 Breve recordatorio de estadística Procesos estocásticos. Ruido blanco. Propagación. Medidas del error. Se escribe X ∼ Nn(m, Σ) y su función de densidad es 1 f (x) = Det(Σ)(2π) n/2 Exp − 1 2 (x − m)T Σ−1 (x − m) . Los intervalos de confianza son ahora regiones de Rn , definidos por P(X ∈ Ω) = PΩ. La forma de estas regiones de confianza es de elipsoides, descritos por la ecuación (x − m) T Σ−1 (x − m) = d 2 , donde d depende de PΩ. Cuanto mayores sean los valores de los autovalores de Σ, mayor será el elipsoide. Las direcciones de los ejes del elipsoide vendrán dados por los autovectores de Σ. 41 / 49 42 / 49 Sistema de navegación autónomo: Navegación inercial. Errores en navegación inercial. Modelos de Error Breve recordatorio de estadística Procesos estocásticos. Ruido blanco. Propagación. Medidas del error. Distribución normal multivariante II Si por ejemplo describimos el error en posición en ejes cuerpo, δr b = [δx δy δz] T , como una normal multivariante con n = 3, de media cero (centrada en el avión) y con matriz de covarianzas: 2 6 Σ = 4 σ 2 x 0 0 0 σ 2 y 0 0 0 σ 2 z Entonces σx representa la magnitud del error ATE (along-track error), σy del error XTE (cross-track error) y σz del error VE (vertical error) y podemos asimilar el movimiento del avión al movimiento del elipsoide, que representa una región de incertidumbre donde se puede encontrar el avión con gran probabilidad. Se verifica que si X ∼ Nn(m x, Σx) e Y ∼ Nn(m y , Σy ) y son independientes, entonces si Z = X + Y resulta Z ∼ Nn(m x + m y , Σx + Σy ). Igualmente AX + b donde A y b son no-aleatorios verifica que AX + b ∼ Nn(Am x + b, AΣxA T ). 43 / 49 Sistema de navegación autónomo: Navegación inercial. Errores en navegación inercial. Modelos de Error Procesos estocásticos. 3 7 5 Breve recordatorio de estadística Procesos estocásticos. Ruido blanco. Propagación. Medidas del error. Un proceso estocástico o variable estocástica no es sino una variable aleatoria X (t) que cambia con el tiempo. Los errores de navegación serán este tipo de variables. Por tanto la media y la covarianza también varían con el tiempo: m(t), Σ(t). Para un proceso, se define la autocorrelación como R(t, τ) = E[X (t)X (τ) T ]. La autocorrelación permite conocer hasta que punto la historia pasada de X influye en su valor actual. Proceso gaussiano: Un proceso gaussiano verifica X (t) ∼ Nn(m(t), Σ(t)), es decir, se distribuye como una normal multivariante cuya media y covarianza varían con el tiempo. 44 / 49

Sistema de navegación autónomo: Navegación inercial. Errores en navegación inercial. Modelos de Error Ruido blanco. Breve recordatorio de estadística Procesos estocásticos. Ruido blanco. Propagación. Medidas del error. Ruido blanco: Se define como ruido blanco un proceso ν(t) que verifica: E[ν(t)] = 0. E[ν(t)ν(t) T ] = σ 2 Id. R(t, τ) = E[ν(t)ν(τ) T ] = δ(t − τ)σ 2 Id, donde δ(x) vale 1 si x = 0 y 0 en cualquier otro caso. La última condición quiere decir que el valor del ruido blanco en un instante es independiente de su valor en cualquier instante anterior. Ruido blanco gaussiano: Es un proceso que cumple las condiciones anteriores, y además es gaussiano. Un buen modelo para las fuentes de error (errores de medida, errores gravitatorios) es δɛ(tk) = b + Dν, donde ν es ruido blanco gaussiano. El valor de b dará la media del error (sesgo, llamado bias en inglés). Sistema de navegación autónomo: Navegación inercial. Errores en navegación inercial. Modelos de Error Propagación del error. Breve recordatorio de estadística Procesos estocásticos. Ruido blanco. Propagación. Medidas del error. 45 / 49 Si en las ecuaciones de propagación del error δx(tk+1) = A(k)δx(tk) + δɛ(tk) sustituimos δɛ(tk) = b + Dν, obtenemos el siguiente modelo de propagación del error: δx(tk+1) = A(k)δx(tk) + b + Dν. Observación: típicamente b también está sometido a un error variable, de forma que b(tk+1) = b(tk) + Dbν b. Para simplificar ignoramos esta variación. Se realizan las siguientes hipótesis: ν es ruido blanco gaussiano con varianza σ2 ν. Inicialmente, δx(t0) ∼ Nn(m0 , Σ0). Si se conocieran perfectamente, entonces Σ0 = 0. Además se tiene la hipótesis de que δx(t0) y ν son independientes. Bajo estas condiciones, se tiene que δx(tk) es un proceso gaussiano, es decir, δx(tk) ∼ Nn(mk, Σk), donde la media y la covarianza verifican la siguiente evolución: Propagación de la media: mk+1 = Amk + b. Propagación de la covarianza: Σk+1 = AΣkAT + σ2 νDDT . 46 / 49 Sistema de navegación autónomo: Navegación inercial. Errores en navegación inercial. Modelos de Error Propagación del error: ejemplo sencillo Breve recordatorio de estadística Procesos estocásticos. Ruido blanco. Propagación. Medidas del error. Supongamos que tuviéramos una ecuación de propagación del error en una dimensión (por ejemplo la posición en el eje x) dada simplemente por: δxk+1 = δxk + ν, donde: La variable temporal k representa minutos, es decir, x6 es el error en posición pasados 6 minutos. ν es ruido blanco gaussiano de varianza σ 2 ν. Inicialmente, δx(t0) = 0. Además δx(tk) y ν son independientes. Entonces aplicando las fórmulas anteriores, δx(tk) ∼ N(mk, σk), donde la media y la varianza verifican: Propagación de la media: mk+1 = amk. Como m0 = 0, se tendrá mk = 0 para todo k. Propagación de la varianza: σ2 k+1 = σ2 k + σ2 nu. Como σ2 0 = 0, se tiene que σ2 k = kσ2 nu. Por tanto la varianza verifica σk = √ kσν. Si por ejemplo x son metros y el ruido blanco tiene σν = 0,1 metros, entonces aunque inicialmente la posición se conoce sin error, pasada una hora σ60 = √ 60 · 0,1 = 0,77, es decir un intervalo 2-σ sería δx ∈ [−1,55, 1,55]. 47 / 49 Sistema de navegación autónomo: Navegación inercial. Errores en navegación inercial. Modelos de Error Medida del error en 2-D Breve recordatorio de estadística Procesos estocásticos. Ruido blanco. Propagación. Medidas del error. Para el caso 2-D (por ejemplo posición sobre un mapa) y si el error está distribuido como X ∼ N2(0, Σ), las regiones de confianza serían elipses: Dado Σ podemos escribir Σ = Pdiag{σ1, σ2}P T donde P es una matriz con autovectores y σi los autovalores. Los autovectores dan la dirección de los ejes de la elipse y los autovalores son proporcionales a su magnitud (cuyo valor exacto dependedel grado de confianza del intervalo). 48 / 49

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 />

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

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

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

Variables aleatorias continuas multidimensionales<br />

Sea una variable aleatoria X ∈ Rn continua multidimensional.<br />

Cada componente <strong>de</strong> X sigue una distribución unidimensional.<br />

Como en el caso unidimensional, se <strong>de</strong>fine una función <strong>de</strong><br />

distribución conjunta, que se calcula mediante la función <strong>de</strong><br />

<strong>de</strong>nsidad f (x).<br />

Igualmente E[g(X )] = <br />

Rn g(y)f (y)dy. Los dos casos<br />

importantes son:<br />

Media: m(X ) = E[X ] = <br />

Rn yf (y)dy.<br />

Matriz <strong>de</strong> covarianzas:<br />

Cov(X ) = E[(X − m(X ))(X − m(X )) T ] = Σ. Es una matriz<br />

simétrica y <strong>de</strong>finida positiva. Los valores <strong>de</strong> la diagonal<br />

representan la varianza <strong>de</strong> cada componente <strong>de</strong> X , mientras<br />

que los valores fuera <strong>de</strong> la diagonal la correlación entre dos<br />

componentes <strong>de</strong> X . Se tiene Σ = E[(X X T ] − m(X )m(X ) T .<br />

Por ejemplo, para n = 3 y escribiendo X = [X , Y , Z]:<br />

2<br />

6<br />

Σ = 4<br />

σ 2 x E[(X − mx )(Y − my )] E[(X − mx )(Z − mz )]<br />

E[(X − mx )(Y − my )] σ 2 y E[(Y − my )(Z − mz )]<br />

E[(X − mx )(Z − mz )] E[(Y − my )(Z − mz )] σ 2 z<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 />

Distribución normal multivariante I<br />

3<br />

7<br />

5<br />

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

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

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

Se escribe X ∼ Nn(m, Σ) y su función <strong>de</strong> <strong>de</strong>nsidad es<br />

1<br />

f (x) =<br />

Det(Σ)(2π) n/2 Exp − 1<br />

2 (x − m)T Σ−1 (x − m) .<br />

Los intervalos <strong>de</strong> confianza son ahora regiones <strong>de</strong> Rn ,<br />

<strong>de</strong>finidos por P(X ∈ Ω) = PΩ.<br />

La forma <strong>de</strong> estas regiones <strong>de</strong> confianza es <strong>de</strong> elipsoi<strong>de</strong>s,<br />

<strong>de</strong>scritos por la ecuación (x − m) T Σ−1 (x − m) = d 2 , don<strong>de</strong> d<br />

<strong>de</strong>pen<strong>de</strong> <strong>de</strong> PΩ. Cuanto mayores sean los valores <strong>de</strong> los<br />

autovalores <strong>de</strong> Σ, mayor será el elipsoi<strong>de</strong>. Las direcciones <strong>de</strong><br />

los ejes <strong>de</strong>l elipsoi<strong>de</strong> vendrán dados por los autovectores <strong>de</strong> Σ.<br />

41 / 49<br />

42 / 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 />

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

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

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

Distribución normal multivariante II<br />

Si por ejemplo <strong>de</strong>scribimos el error en posición en ejes cuerpo,<br />

δr b = [δx δy δz] T , como una normal multivariante con n = 3,<br />

<strong>de</strong> media cero (centrada en el avión) y con matriz <strong>de</strong><br />

covarianzas:<br />

2<br />

6<br />

Σ = 4<br />

σ 2 x 0 0<br />

0 σ 2 y 0<br />

0 0 σ 2 z<br />

Entonces σx representa la magnitud <strong>de</strong>l error ATE<br />

(along-track error), σy <strong>de</strong>l error XTE (cross-track error) y σz<br />

<strong>de</strong>l error VE (vertical error) y po<strong>de</strong>mos asimilar el movimiento<br />

<strong>de</strong>l avión al movimiento <strong>de</strong>l elipsoi<strong>de</strong>, que representa una<br />

región <strong>de</strong> incertidumbre don<strong>de</strong> se pue<strong>de</strong> encontrar el avión con<br />

gran probabilidad.<br />

Se verifica que si X ∼ Nn(m x, Σx) e Y ∼ Nn(m y , Σy ) y son<br />

in<strong>de</strong>pendientes, entonces si Z = X + Y resulta<br />

Z ∼ Nn(m x + m y , Σx + Σy ).<br />

Igualmente AX + b don<strong>de</strong> A y b son no-aleatorios verifica que<br />

AX + b ∼ Nn(Am x + b, AΣxA T ).<br />

43 / 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 />

Procesos estocásticos.<br />

3<br />

7<br />

5<br />

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

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

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

Un proceso estocástico o variable estocástica no es sino una<br />

variable aleatoria X (t) que cambia con el tiempo. Los errores<br />

<strong>de</strong> navegación serán este tipo <strong>de</strong> variables.<br />

Por tanto la media y la covarianza también varían con el<br />

tiempo: m(t), Σ(t).<br />

Para un proceso, se <strong>de</strong>fine la autocorrelación como<br />

R(t, τ) = E[X (t)X (τ) T ]. La autocorrelación permite conocer<br />

hasta que punto la historia pasada <strong>de</strong> X influye en su valor<br />

actual.<br />

Proceso gaussiano: Un proceso gaussiano verifica<br />

X (t) ∼ Nn(m(t), Σ(t)), es <strong>de</strong>cir, se distribuye como una<br />

normal multivariante cuya media y covarianza varían con el<br />

tiempo.<br />

44 / 49

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

Saved successfully!

Ooh no, something went wrong!