Spacecraft trajectory optimization with GA in Matlab (on/off constant thrust)
19 views (last 30 days)
I would like to kindly ask for support or any advice on how to implement my problem in Matlab, perhaps using the (Global) Optimization Toolbox, and whether it is even possible.
My goal is to find a time-history of control (σ,) during a fuel-optimal spacecraft rendezous with constant low-thrust.
The control variables are defined as:
- σ - total thrust acceleration
- - thrust acceleration projections
And the state vector is:
The objective is to minimise:
- The state equations in state-space representation (CW equations):
- Control variables constraint: and σ can be equal to either 0 or (on/off)
- Initial conditions: given
- Terminal constraints: given
- Final time inequality constraint:
As I understand, this problem can be categorised as a dynamical optimization problem, that involves integer programming. Could it be solved in Matlab, perhaps using the Genetic Algorithms which I believe allow for integer programming?
Are there any available examples on how to implement a spacecraft (or not necessarily spacecraft) trajectory optimization problem in Matlab, using GA? I have been looking for examples for a very long time, but I could not find any. In fact, I could not find any examples even without the on/off thursting constraints, and I would be grateful if anyone could direct me to other spacecraft trajectory optimization implementations in Matlab, perhaps using the fmincon function.
Thank you very much.
Alan Weiss on 21 Aug 2022
You might be interested in this example: Discretized Optimal Trajectory, Problem-Based. The problem formulation is different than yours, so it is probably not directly applicable, but you might be able to make it work for you. One thing to note: I recently found out that this sort of optimal trajectory problem works better when you lower the optimality tolerance, as described here: https://www.mathworks.com/matlabcentral/answers/1774135-possible-bug-with-coneprog.
MATLAB mathematical toolbox documentation