From Orbit Matching to Docking: A Hierarchical Reinforcement-Learning Approach to Multiphase Rendezvous Trajectory Design
Authors: Mariko Storey-Matsutani, Anatoly Zlotnik, Lee Holguin, Zigfried Hampel-Arias, Richard Linares
Organizations: Department of Aeronautics and Astronautics, Massachusetts Institute of Technology, Cambridge, MA 02139, USA · Los Alamos National Laboratory, Los Alamos, NM, USA
As the complexity of space missions increases, autonomous spacecraft guidance must accommodate transitions between mission phases with distinct dynamics, objectives, and constraints. Rendezvous and proximity operations (RPO) are particularly challenging because they span multiple operating regimes, from orbital-scale matching and far-range approach to close-range maneuvering and terminal docking. Existing reinforcement-learning approaches generally address only a subset of this sequence or rely on manually defined switching conditions between controllers. This study presents a hierarchical reinforcement-learning framework for autonomous, three-degree-of-freedom RPO that spans the complete trajectory from orbit matching to docking. Five specialized policies are trained for orbit matching, far-range approach, mid-range approach, close-range approach, and docking. A learned high-level manager policy then selects among these phase-specific controllers based on the evolving spacecraft state. The hierarchical approach is evaluated over the full end-to-end rendezvous sequence and compared with a fixed sequential-invocation heuristic to assess trajectory performance and transition behavior. Training and inference metrics are also analyzed for the individual phase policies and the complete hierarchical system.
Figures & tables
Notation/Variable
Definition
r
Position
v
Velocity
F
Thrust force (control)
m
Mass
t
Time
[B,C,T]
[Body, Chaser, Target]
TABLE I : Notation and Variable Definitions
Phase
Initial Condition
Terminal Condition
Model Parameters
Orbit Match
Phase: [-90.0, 90.0] ∘
ρtol=10,000 m
Δvmax=10,000 m/s
SMA: [-20k, 20k] m
vtol=50 m/s
Δvk,max=1200 m/s
Ecc.: [0.0, 0.01]
ρmax=25,000 km
Inc.: [-45.0, 45.0] ∘
kmax= 180
RAAN: [-45.0, 45.0] ∘
Far Range
∥rrel,0L∥ : [3000, 10000] m
ρtol=1000 m
Δvmax=180 m/s
TABLE II : Rendezvous Approach Phases
Reward
Definition
Position Penalty
ϕ={rrel/100rrel<rproxotherwise
Velocity Penalty
Constant: ϕ=−∣∣vrel∣∣
Distance-Scaled: ϕ=−max(rrel,rmin)∣∣vrel∣∣
Progress
ϕ=sϵprev−ϵcurr , where
ϵ∈ {altitude difference ∣htgt−h∣ , (s=8M);
eccentricity difference ∣etgt−e∣ , (s=0.005);
TABLE III : Reward Component Definitions
Single-Phase Models
Parameter
Orbit Match
Far Range
Mid Range
Close Range
Dock
Multi Phase
Reward Weights
Position Penalty
-
-
-
-1.0
-1.0
-
Velocity Penalty
Distance-Scaled
-
1.0
1.0
-
-
1.0
Constant
-
-
-
1.0
1.0
-
TABLE IV : Reward Weights Across All Controllers
Parameter
Orbit Match
Far Range
Mid Range
Close Range
Docking
Multi Phase
Network Size
[128, 128]
[128, 128]
Activation
tanh
tanh
Learning Rate
0.0003
0.0003
Batch size
256
256
Discount
0.998
0.99
0.99
0.99
0.99
0.99
GAE λ
0.95
0.95
TABLE V : RL Model Parameters
Phase
Success
Min. Range
Final Range
Final Vel.
Δv
Rate
(m)
(m)
(m/s)
(m/s)
Orbit Match
97.0%
5380.10
5383.52
2.005
5660.64
Far Range
100.0%
249.96
269.97
0.478
26.83
Mid Range
100.0%
53.52
64.48
0.296
8.99
Close Range
100.0%
9.80
9.95
0.023
2.66
Docking
100.0%
0.98
0.98
0.009
0.25
TABLE VI : Single-phase policy Evaluation ( N=1000 episodes each, deterministic policy, held-out seed, all except success rate are median values).
TABLE VII : Comparison of learned manager vs. heuristic baseline performance metrics with N=1000 (Format: Median (Mean ± Std))
Fig. 3 : Comparison of learned manager vs. heuristic baseline for relative position and velocity through each trajectory. Only successful trajectories are displayed.
Phase
Learned Manager
Heuristic
Orbit Match
97.141% ± 1.052%
97.491% ± 0.803%
( 25.8836±7.4155 h)
( 25.9421±6.3255 h)
Far Range
1.161% ± 0.380%
0.970% ± 0.303%
( 0.2868±0.0275 h)
( 0.2456±0.0448 h)
Mid Range
0.484% ± 0.154%
0.375% ± 0.121%
( 0.1194±0.0109 h)
( 0.0945±0.0174 h)
TABLE VIII : Time breakdown of learned manager and heuristic baseline by phase (percentages and hours reported as mean ± std across episodes)
Timesteps
Training
Throughput
Model
Time (h)
(fps)
Orbit Match
Stage 1
491,520
0.06
2,296
Stage 2
430,848
0.05
2,256
Stage 3
2,179,072
0.29
2,072
Total
3,101,440
0.40
2,129
TABLE IX : Training metrics per model
Params
Checkpoint
Latency
Control Step
Model
(KB)
(ms)
Margin
Docking
35,847
462
0.095±0.032
1.05×104
Close Range
35,847
464
0.117±0.072
8.53×103
Mid Range
35,847
464
0.093±0.037
1.08×104
Far Range
35,847
464
0.096±0.053
3.75×104
Orbit Matching
36,871
476
0.110±0.053
9.35×106
TABLE X : Model size and inference latency (N = 1,000)
Metrics
Learned
Heuristic
Success rate
0.994
0.970
Mean episode wall-clock (s)
0.961±0.152
0.917±0.221
Median episode wall-clock (s)
0.932
0.905
Expected wall-clock per success (s)
0.967
0.945
Computational Breakdown
Manager Inference (%)
0.21
-
TABLE XI : Computational breakdown for the learned and heuristic hierarchical managers (N = 1,000 paired episodes)
Spacecraft rendezvous and proximity operations cover the navigation, guidance, and control needed to approach, inspect, dock with, or capture spacecraft safely.
Autonomous rendezvous and proximity operations (RPO) in adversarial orbital environments require guidance architectures balancing target pursuit, safety preservation, and real-time adaptability under dynamically evolving interaction conditions. Although learning-based approaches show promise, their application to safety-critical orbital robotics remains limited by concerns regarding interpretability, robustness, and constraint awareness. This work presents an adaptive Model Predictive Control (MPC) framework for autonomous spacecraft RPO in multi-agent adversarial scenarios. The proposed architecture combines a constrained receding-horizon MPC formulation with a data-driven supervisory tuning layer that adjusts controller parameters from offline closed-loop evaluation and online interaction geometry. Relative motion follows Clohessy-Wiltshire (CW) dynamics, enabling computationally efficient finite-horizon prediction and real-time quadratic optimization. The MPC formulation incorporates actuator limits, predictive keep-out-zone constraints, slack-variable feasibility handling, and optional Control Barrier Function (CBF) safety filtering. Rather than generating thrust commands directly, the adaptive layer modifies interpretable MPC parameters, including tracking weights, safety penalties, minimum-separation objectives, and keep-out-zone objectives. The framework was evaluated in the official Kerbal Space Program Differential Game (KSPDG) Capture-the-Satellite environment through Monte Carlo simulations. Results demonstrate improved closed-loop robustness, adaptive maneuvering behavior, and rendezvous performance compared with fixed-parameter MPC while preserving safety-aware operation and real-time feasibility, providing a modular, interpretable foundation for adaptive spacecraft RPO.
Luca Sportelli, Tyler Barr, Cagri Kilic +1
Department of Aerospace Engineering, Embry–Riddle Aeronautical University, Daytona Beach, FL 32114
Emerging mission classes such as on-orbit servicing, satellite inspection, and active debris removal require trajectory design methods that are adaptable to a variety of mission scenarios. We present a diffusion-based trajectory generation approach for rendezvous and proximity operations (RPO) that enables flexible configuration of mission constraints. First, individual energy-based diffusion models are trained to satisfy distinct constraints such as approach cone and sensor line-of-sight from a set of optimized trajectories. Then, at inference time, the learned energy models can be composed with one another, or with an analytically defined energy field, to enforce specific constraint combinations. We validate this framework with the composition of a learned approach cone model and a learned sensor line-of-sight model, as well as a learned approach cone model and synthetic obstacle avoidance model, both of which yield constraint satisfaction rates that are within 1 percentage point of the single-constraint models or higher. These results indicate that our compositional diffusion framework can provide a modular approach to RPO trajectory design and enable reconfiguration for new constraint combinations without requiring model retraining.
Mariko A. Storey-Matsutani, Richard Linares
Graduate Student, Department of Aeronautics and Astronautics, Massachusetts Inst. of Technology, Cambridge MA · Associate Professor, Department of Aeronautics and Astronautics, Massachusetts Inst. of Technology, Cambridge MA
Autonomous vertical landing of reusable launch vehicles is a critical yet formidable challenge in aerospace engineering, primarily due to highly nonlinear dynamics, stringent precision require- ments, and inherent environmental uncertainties. While reinforcement learning (RL) has emerged as a promising alternative to traditional model-based control, existing approaches frequently rely on flat training paradigms and highly simplified initial conditions, drastically limiting their real-world applicability. This paper proposes a comprehensive curriculum-driven deep reinforcement learning framework designed to master both far-range navigation and precision landing. By employing a three-stage curriculum augmented with a critic warm-up procedure and a rehearsal mechanism to prevent catastrophic forgetting, alongside a high-dimensional observation vector utilizing a full 3 X 3 rotation-matrix, the agent is trained to navigate from an initial horizontal distance of 1,000 to 1,500 units vastly farther than prior benchmarks before executing a soft touchdown. Crucially, we intro- duce a redesigned reward function that reduces unnecessary computational overhead and addresses limitations of earlier formulations, including poorly calibrated reward components and potentially conflicting optimization objectives. The redesigned function bounds the per-step reward magnitude by (O(1)), while its progress term accumulates the total distance closed over the trajectory. The available 100,000-episode results report a 99.98% success rate for the proposed policy and 68.31% for the baseline under the evaluation conditions. In the evaluation, the proposed policy retained 69.1% of its initial fuel capacity and achieved a mean final distance from the pad centre of 7.07 units, compared with 13.10 units for the baseline. Landing visualizations are available in Video 1 and Video 2.
Atin Chowdhury, Piyush Joshi
Department of Computer Science and Engineering, Indian Institute of Information Technology, Sri City, Andhra Pradesh, India