I am trying to perform an orbit propagation using the built-in PartialDerivativesEquations class to compute the State Transition Matrix wrt some initial conditions. My intention is to propagate a given covariance ellipsoid and see how it evolves after a certain period of time. This is fairly straigh-forward when working in Cartesian parameters, but would it be possible to do using Keplerian parameters? My goal is to check if the results using Keplerian would be in any way better (given certain statistical figures of merit) than the results using Cartesian, due to a (possible) better linearization of the dynamics.
I have seen an old (circa Aug. 2017) post regarding this subject, but no clear solution was given. For now, what I have tried is to set the OrbitType to “Keplerian” for the NumericalPropagator, without much success (after closer inspection, the “computeDerivatives” method always uses Cartesian anyway, just like the link above suggests).
Thank you in advance. Best regards,