Representación de satélites

Índice

1.- Representación de orbitas ECI        2

2.- Representación de orbitas ECEF        3

3.- Proyección sobre superficie        4

4.- Bonus: Animación        6

1.- Representación de orbitas ECI

        Aprovechando el código proporcionado por el profesor, realizamos un bucle for a lo largo de todos los SVs del almanaque para representar todas las orbitas.

Para representar la posición exacta del satélite es suficiente con eliminar la parte de variación temporal de la función (‘+900*t’), obteniendo la posición del satélite en el instante en el que se generó el fichero RINEX.

El resultado que obtenemos es el siguiente:

[pic 2]

Como vemos es un modelo ideal, ya que no tiene en cuenta el movimiento de rotación de la tierra.

2.- Representación de orbitas ECEF

Modificamos el fichero ‘satposin.m’ formulando una nueva ecuación para Omega, en la que se tenga en cuenta el movimiento de la tierra.

[pic 3]

Sin modificar nada más volvemos a modificar el script de Matlab y obtenemos la representación de las orbitas considerando el movimiento de la tierra (ECEF)

[pic 4]

Como vemos y sabemos, las orbitas tienen forma sinusoidal respecto al ecuador. Además tenemos los polos sin cubrir por los satélites.

3.- Proyección sobre superficie

        Aprovechando el código de ECEF, simplemente en la función del plot, añadimos una nueva línea en la que reducimos la órbita al radio de la tierra y obtener así la proyección[pic 5]

La proyección de los satélites sobre la superficie es la siguiente:

Primero representamos solo 2 satélites para que se vea claramente, para ello modificamos el valor de ‘n’ con n = 2 en el bucle for.

[pic 6]

Como se aprecia claramente la órbita describe una especie de seno centrado en el ecuador. Y de amplitud que no llega a los polos, ya que las zonas polares no disponen de cobertura.

...