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

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

Variables <strong>de</strong> error. Error en actitud<br />

Mo<strong>de</strong>lo <strong>de</strong> propagación linealizado<br />

El canal vertical<br />

Por tanto podremos escribir, como antes,<br />

δ ˙φ = Cφpδp + Cφv δv n + Cφφδφ − δωb b/i .<br />

Es una ecuación lineal en los errores, don<strong>de</strong> las matrices están<br />

<strong>de</strong>finidas en función <strong>de</strong> la estimación <strong>de</strong>l INS, y con un<br />

términos forzante: el error en los giróscopos δωb b/i .<br />

Si ponemos todos los errores juntos, llegamos a:<br />

d<br />

dt<br />

⎡<br />

⎣<br />

δp<br />

δv n<br />

δφ<br />

⎤<br />

⎡<br />

⎦ = ⎣<br />

Cpp Cpv 0<br />

Cvp Cvv Cvφ<br />

Cφp Cφv Cφφ<br />

⎤ ⎡<br />

⎦ ⎣<br />

δp<br />

δv n<br />

δφ<br />

⎡<br />

⎤<br />

⎡<br />

⎦+ ⎣<br />

Caδab NG<br />

−δωb b/i<br />

⎤<br />

0<br />

n<br />

+ δG ⎦<br />

δp<br />

A<strong>de</strong>más estarán las condiciones iniciales: ⎣ δv n ⎦ (t = 0).<br />

δφ<br />

Éste es el mo<strong>de</strong>lo <strong>de</strong> propagación <strong>de</strong>l error <strong>de</strong>l INS. Puesto<br />

que el término forzante es <strong>de</strong>sconocido (y se mo<strong>de</strong>la mediante<br />

la estadística) es una ecuación diferencial estocástica.<br />

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

⎤<br />

Variables <strong>de</strong> error. Error en actitud<br />

Mo<strong>de</strong>lo <strong>de</strong> propagación linealizado<br />

El canal vertical<br />

Ecuación <strong>de</strong>l error en el canal vertical I<br />

Si trabajamos sólo con el error en h y VD, y <strong>de</strong>spreciamos<br />

todos los términos excepto el gravitatorio, llegamos a la<br />

siguiente ecuación:<br />

δ ˙ h = −δVD<br />

δ ˙VD −2<br />

µe<br />

δh.<br />

(Re + ˆh) 3<br />

Por otro lado po<strong>de</strong>mos aproximar en el <strong>de</strong>nominador<br />

Re + ˆ h Re. Teniendo en cuenta que la aceleración <strong>de</strong> la<br />

gravedad al nivel <strong>de</strong>l mar g0 = µe<br />

R 2 e<br />

δ ˙h = −δVD<br />

δ ˙ VD − 2g0<br />

δh.<br />

Re<br />

, tendríamos las ecuaciones:<br />

Escribiéndolo como una única ecuación para δh: δ ¨ h = 2g0<br />

Re δh.<br />

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

Variables <strong>de</strong> error. Error en actitud<br />

Mo<strong>de</strong>lo <strong>de</strong> propagación linealizado<br />

El canal vertical<br />

Ecuación <strong>de</strong>l error en el canal vertical II<br />

La solución <strong>de</strong> la ecuación diferencial es:<br />

q<br />

2g0<br />

Re t + C2e −<br />

q<br />

2g0<br />

δh = C1e<br />

Re t , don<strong>de</strong> las constantes son función<br />

<strong>de</strong> las condiciones iniciales <strong>de</strong> altura y velocidad vertical.<br />

Éstas ecuaciones son inestables! El primer término crece hasta<br />

el infinito.<br />

Físicamente, lo que suce<strong>de</strong> es lo siguiente: si hay un error <strong>de</strong><br />

altitud, p.ej. el INS piensa que el avión está más alto <strong>de</strong> lo<br />

que realmente está, el mo<strong>de</strong>lo <strong>de</strong> gravedad predice que la<br />

gravedad es menor <strong>de</strong> lo que es, con lo que el INS predice que<br />

el avión se eleva, es <strong>de</strong>cir, el error inicial se amplifica!<br />

Éste resultado se mantiene si no se <strong>de</strong>sprecian los términos<br />

que no se han consi<strong>de</strong>rado. Por tanto el canal vertical <strong>de</strong>l INS<br />

es inestable y no se pue<strong>de</strong> usar por sí sólo; empleando otras<br />

medidas (p.ej. barométricas) es posible compensar el canal<br />

vertical y obtener una medida fiable <strong>de</strong> la altura.<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 />

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

Hemos visto que las ecuaciones <strong>de</strong> propagación <strong>de</strong>l error son<br />

<strong>de</strong>l tipo δ ˙x = A(ˆx)δx + δɛ, don<strong>de</strong> δx son las variables <strong>de</strong><br />

navegación (posición, velocidad, actitud) y los δɛ las fuentes<br />

<strong>de</strong> error. Estas fuentes son:<br />

Errores en el mo<strong>de</strong>lo <strong>de</strong> gravedad δg n .<br />

Errores en los sensores inerciales δa b NG , δωb b/i .<br />

Aparte está el error en las condiciones iniciales δx(t0).<br />

Si discretizamos estas ecuaciones en el tiempo, podríamos<br />

escribir un mo<strong>de</strong>lo algo más sencillo:<br />

δx(tk+1) = A(tk)δx(tk) + δɛ(tk).<br />

¿Cómo se mo<strong>de</strong>lan los errores? ¿Cómo se interpretan las<br />

ecuaciones que contienen errores?<br />

Para respon<strong>de</strong>r a estas preguntas es necesario recordar<br />

algunos conceptos estadísticos.<br />

35 / 49<br />

36 / 49

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

Saved successfully!

Ooh no, something went wrong!