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

Mecanización en ejes e<br />

La IMU: sensores inerciales<br />

Mecanización en ejes n y en ejes e<br />

Alineamiento inicial<br />

En ocasiones, por motivos <strong>de</strong> integración INS-GPS, conviene<br />

mecanizar las ecuaciones en los ejes e (en los que trabaja el<br />

GPS).<br />

Se llega a las siguientes ecuaciones para velocidad y posición:<br />

Velocidad: d<br />

dt v e <br />

= −2 ωe ×<br />

e/i v e + ae NG + g e =<br />

<br />

−2 ωe ×<br />

e/i v e + (C n e ) T (C b n ) T ab e<br />

NG + g<br />

Posición: d<br />

dt r e = v e .<br />

Habría que escribir C n e en función <strong>de</strong> r e y v e , escribir un<br />

mo<strong>de</strong>lo <strong>de</strong> g e , y escribir la ecuación <strong>de</strong> la actitud, y se llegaría<br />

a un esquema similar al anterior.<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 />

Alineamiento inicial I<br />

La IMU: sensores inerciales<br />

Mecanización en ejes n y en ejes e<br />

Alineamiento inicial<br />

Supongamos que tenemos el avión en reposo en un<br />

aeropuerto, y es necesario inicializar el INS con un “fix”.<br />

¿Cómo se haría?<br />

Evi<strong>de</strong>ntemente, se tiene que φ, λ y h son las <strong>de</strong>l aeropuerto, o<br />

incluso con mayor precisión, las tomadas <strong>de</strong> un sistema GPS.<br />

Puesto el avión está en reposo, v n = 0.<br />

Queda encontrar el valor inicial <strong>de</strong> actitud, es <strong>de</strong>cir,<br />

C b n (t = 0). Para ello se usa la medida obtenida <strong>de</strong> giróscopos<br />

y acelerómetros (en reposo).<br />

De la ecuación fundamental <strong>de</strong> la navegación se tiene:<br />

<br />

0 = − ωn n/e + 2ωn ×<br />

e/i 0 + an NG + g n , luego an NG = −g n y<br />

por tanto ab NG = C b n an NG = −C b n g n .<br />

Por otro lado es claro que ωb b/n = ωb b/i − ωb e/i − ωb n/e y<br />

evi<strong>de</strong>ntemente ωb b/n = 0 y ωb n/e = 0.<br />

Por tanto:ωb b/i = ωb e/i = C b n ωn e/i .<br />

17 / 49<br />

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

Alineamiento inicial II<br />

La IMU: sensores inerciales<br />

Mecanización en ejes n y en ejes e<br />

Alineamiento inicial<br />

Tenemos por tanto dos ecuaciones: a b NG = −C b n g n y<br />

ω b b/i = C b n ω n e/i . Llamando a las medidas x b 1 = ab NG y<br />

x b 2 = ωb b/i<br />

x n 2 = ωn e/i<br />

, y <strong>de</strong>notando los mo<strong>de</strong>los como y n<br />

, se tiene que<br />

x n 1 = C b n (0)y b<br />

1 , x n 2 = C b n (0)y b<br />

2<br />

1 = −g n y<br />

Tendríamos 6 medidas (las componentes <strong>de</strong> dos vectores) para<br />

9 grados <strong>de</strong> libertad (las entradas <strong>de</strong> la matriz).<br />

Es necesario pues “generar” una medida adicional<br />

in<strong>de</strong>pendiente.<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 />

Alineamiento inicial III<br />

La IMU: sensores inerciales<br />

Mecanización en ejes n y en ejes e<br />

Alineamiento inicial<br />

Llamemos x 3 = x 1 × x 2. Obsérvese que este vector se pue<strong>de</strong><br />

escribir como x b ×<br />

1 x b<br />

2 en el sistema <strong>de</strong> <strong>referencia</strong> b, don<strong>de</strong> X<br />

es la matriz antisimétrica que representa el producto vectorial.<br />

Por otro lado se tiene que x b <br />

×<br />

1 = C b<br />

n (0) y n<br />

×<br />

C<br />

1<br />

n b (0). Por<br />

<br />

tanto x b 3 = x b ×<br />

1 x b<br />

2 = C b n (0) y n<br />

×<br />

C<br />

1<br />

n b (0)C b n (0)y n<br />

2 =<br />

C b <br />

n (0) y n<br />

×<br />

y<br />

1<br />

n<br />

2 . Por tanto <strong>de</strong>notando y 3 = y 1 × y , se tiene<br />

2<br />

que x b 3 = C b n (0)y n<br />

3 .<br />

Escribiendo la matriz A como la matriz cuyas columnas son<br />

x b 1 , x b 2 y x b 3 , y la matriz B como la matriz cuyas columnas son<br />

y n,<br />

y n<br />

19 / 49<br />

y y n<br />

1 2 3 , se tiene: A = C b n (0)B y por tanto C b n (0) = AB−1 .<br />

No se han tenido en cuenta los errores <strong>de</strong> medida: C b n (0)<br />

probablemente no saldría ortonormal (habría que emplear un<br />

algoritmo más sofisticado que tuviera en cuenta los errores <strong>de</strong><br />

medida).<br />

20 / 49

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

Saved successfully!

Ooh no, something went wrong!