Non-Convex Global Optimization as an Optimal Stabilization Problem: Dynamical Properties
Yuyang Huang, Dante Kalise, Hicham Kouhkouh
TL;DR
The paper casts global non-convex optimization as a deterministic optimal stabilization problem, introducing three connected control formulations (evolutive discounted, stationary discounted, evolutive non-discounted) and studying their trajectories via occupation measures. It proves that (quasi-)optimal trajectories globally and practically converge to the set of global minimizers $\mathfrak{M}$: for any tolerance $\eta$, suitable $\lambda$ and horizon $t$ exist such that trajectories stay within $\eta$ of $\mathfrak{M}$ after time $\tau$, with the non-minimizer time vanishing as $\lambda\to0$ or $t\to\infty$. The results hold across all three problem variants, and are complemented by explicit bounds on the time spent outside neighborhoods of $\mathfrak{M}$ via occupation measures, as well as Lyapunov-type stability guarantees. The approach eschews ergodic Hamilton–Jacobi–Bellman limits, offering a direct dynamical mechanism for global optimization grounded in deterministic control and measure-theoretic analysis. This provides a principled interpretation of optimization dynamics and potential guidance for designing robust global-search algorithms grounded in optimal-control theory.
Abstract
We study global optimization of non-convex functions through optimal control theory. Our main result establishes that (quasi-)optimal trajectories of a discounted control problem converge globally and practically asymptotically to the set of global minimizers. Specifically, for any tolerance $η> 0$, there exist parameters $λ$ (discount rate) and $t$ (time horizon) such that trajectories remain within an $η$-neighborhood of the global minimizers after some finite time $τ$. This convergence is achieved directly, without solving ergodic Hamilton-Jacobi-Bellman equations. We prove parallel results for three problem formulations: evolutive discounted, stationary discounted, and evolutive non-discounted cases. The analysis relies on occupation measures to quantify the fraction of time trajectories spend away from the minimizer set, establishing both reachability and stability properties.
