I. ABSTRACTLow-thrust trajectories with variable radial thrust is studied in this paper. The problem is tackled by solving the HamiltonJacobi-Bellman equation via State Dependent Riccati Equation(STDE) technique devised for nonlinear systems. Instead of solving the two-point boundary value problem in which the classical optimal control is stated, this technique allows us to derive closed-loop solutions. The idea of the work consists in factorizing the original nonlinear dynamical system into a quasi-linear state dependent system of ordinary differential equations. The generating function technique is then applied to this new dynamical system, the feedback optimal control is solved. We circumvent in this way the problem of expanding the vector field and truncating higher-order terms because no remainders are lost in the undertaken approach. This technique can be applied to any planet-to-planet transfer; it has been applied here to the Earth-Mars low-thrust transfer.
II. INTRODUCTIONHistorically, the optimal low-thrust transfers have been tackled first with indirect and then with direct methods. The former stem from the Pontryagin's maximum principle that uses the calculus of variations [1,2]; the latter aim at solving the problem via a standard nonlinear programming procedure [3]. Even if it can be demonstrated that both approaches lead to the same result [4], the direct and indirect methods have different advantages and drawbacks, but they require in any case the solution of a complex set of equations: the Euler-Lagrange equations (indirect methods) and the Karush-Kuhn-Tucker equations (direct methods). The guidance designed with these methods is obtained in an open-loop context. In other words, the optimal path, even if minimizing the prescribed performance index, is not able to respond to any perturbation that could alter the state of the spacecraft. Furthermore, if the initial conditions are slightly varied (e.g. the launch date changes), the optimal solution needs to be recomputed again. The outcome of the classical problem is in fact a guidance law. expressed as a function of the time,u = u(t), t ∈ [t0, tf ], being t0 and tf the initial and final time, and u the control vector, respectively. This paper deals with the optimal feedback control problem applied to the low thrust interplanetary trajectory design. With this approach the solutions that minimize the performance index are also functions of the generic initial state x0; the outcome is in fact a guidance law written as u = u(x0, t0, t), t t ∈ [t0, tf ]. This represents a closed-loop solution: given the initial conditions (t0,x0) it is possible to extract the optimal control law that solves the optimal control problem. Moreover, if for any reason the state is perturbed and assumes the new value (t0 , x0 ) = (x0 + δx, t0 + δt), we are able to compute the new optimal solution by simply evaluating so avoiding the solution of another optimal control problem. This property holds by virtue of the closed loop characteristics of the control law that can be vie...