- Grafico de balistica con diferentes valores de velocidad del viento
- En el grafico de la izquierda se observan los comportamientos que tiene el satelite con la orbita predicha "basica" y la real entregada en el archivo EOF. En el grafico de la derecha se observa que la diferencia final entre ambas orbitas calculadas es de: 597,9 km.
- En el grafico se observa la deriva con un numero de subdivisiones igual a 1, con las distintas soluciones de "odeint" y "eulerint". Se ve que "eulerint"(color azul) demuestra un comportamiento oscilante que se aleja cada vez mas de la Tierra, siendo muy poco efectivo en su desempeño, tanto asi, que la deriva entre ambas soluciones al momento final del tiempo es de 19.689,96 km.
- Los tiempos de realizacion de las soluciones son: Odeint = 0.30618 segundos , Eulerint = 0.99744 segundos
- Al realizar varias iteraciones para analizar el comportamiento de la funcion "Eulerint", se vio que al sobrepasar el Nsub = 1000 se veia un porcentaje de error cercano al buscado. Hasta que se realizó con un Nsub = 1700, el cual arrojó el resultado mas positivo de los experimentados (1000, 1200 y 1500). La iteracion con Nsub = 1700, con Eulerint tardó 2022,356 segundos, los cuales son cerca de 33 min de demora. Al momento que python mostró el grafico, vi que en el codigo habia establecido ciertos limites que no permiten visualizar bien el comportamiento de ambas curvas, eso si, se puede ver que la linea azul esta muy cerca del 0, esto significa que el error es muy pequeño.
- Se puede observar que el grafico de la izquierda no hay grandes cambios visuales con respecto a los obtenidos en la Parte 1, pero con el grafico de la derecha queda demostrado que la incorporacion de J2 en la funcion ayudó a disminuir en un 50% nuestro error aproximadamente con un Delta = 319.89 km, lo cual es bueno, pero no lo suficiente como ocurrió con los graficos presentados en Canvas por el profesor.
- En este grafico tambien queda demostrado que hay mejores en el rendimiento del codigo, ya que la curva roja (Odeint) disminuye su "pendiente" y eso indica que el calculo mediante esta solucion esta mas cerca que sin J2.
- Este codigo se demoro en correr 3.793 segundos
- Se puede observar que el grafico de la izquierda no hay grandes cambios visuales con respecto a los obtenidos en la Parte 1, pero con el grafico de la derecha se ve que la incorporacion de J3 ayudó a disminuir de manera leve el error con solo J2 con un Delta = 317.59 km, lo cual es bueno, ya que se respeta lo presentado en los graficos de Canvas.
- En este grafico, no se nota visualmente una diferencia con el anterior, pero si hay un diferencias en los calculos, viendose reflejado en un disminucion en el punto mas alto de la curva roja, lo que minimiza el error.
- Curiosamente, este codigo se demoro en correr 3.657 segundos, siendo menor al tiempo empleado con solo J2 incluido.