I am integrating some equations that represent an astrodynamics problem. They are defined by a vector field such as:
where x are the states (e.g. position, velocity, mass...). The dynamics that I'm integrating are fed by a control u, that is computed based on the states and time: . This control u is computed inside the vector field, which we could explicitely state as . I am integrating with an ODE solver (e.g. ode45()), and thus the solver computes the control at every time step.
Now, what I want is to avoid computing the control at every time step, but rather only at certain times, with a given frequency (e.g. let's say I want to compute the control at 0.1 Hz, every 10s). How would you recommend me to approach this?
I hope I was clear enough. Thank you for any help.