La formación de satélites enfrenta desafíos significativos para mantener sus configuraciones deseadas debido a diversas perturbaciones orbitales, particularmente en entornos de órbita baja terrestre. Este artículo presenta un enfoque novedoso para generar trayectorias de referencia óptimas en combustible para formaciones de satélites en pista, incorporando tanto la oblaticidad de la Tierra (perturbación) como la no linealidad inherente del problema de dos cuerpos. El problema de control óptimo indirecto resultante se resuelve utilizando Redes Neuronales de Pontryagin (PoNNs). El método propuesto transforma el problema convencional de valor en la frontera de dos puntos en un problema de programación matemática, lo que permite el cálculo eficiente de trayectorias óptimas. La efectividad de nuestro enfoque se valida a través de extensas simulaciones numéricas en diferentes inclinaciones del satélite principal (0-90 grados) y distancias de separación cruzada (1-400 km), demostrando reducciones significativas en el consumo anual de combustible en comparación con enfoques convencionales. La viabilidad de estas trayectorias óptimas se verifica mediante simulaciones en bucle cerrado utilizando un controlador PD, confirmando su aplicabilidad práctica en escenarios de misión realistas. Esta investigación contribuye a mejorar la sostenibilidad a largo plazo de las misiones de vuelo en formación de satélites al optimizar la eficiencia del combustible mientras se mantienen formaciones precisas.