Learning Stabilizing Policies via an Unstable Subspace Representation
Leonardo F. Toso, Lintao Ye, James Anderson
TL;DR
This work addresses learning to stabilize unknown high-dimensional LTI systems by a two-phase, model-free approach that first learns the left unstable subspace and then performs discounted LQR on the learned low-dimensional subspace. By focusing on the unstable modes, the method reduces the effective control dimension and sample complexity, achieving a non-asymptotic guarantee that scales as $\tilde{O}(\ell^2 d_U)$ (plus an $\mathcal{O}(d_X)$ term) rather than $\tilde{O}(d_X^2 d_U)$ for the full-space problem. Key components include finite-sample guarantees for estimating the left unstable representation from adjoint data, a zeroth-order gradient estimator for the low-dimensional problem, and an explicit discount-annealing scheme that leads to a stabilizing controller $K=\theta\widehat{\Phi}^T$ with high probability. The proposed approach accommodates non-diagonalizable dynamics via Jordan form and yields substantial practical gains in instability-driven dimensionality reduction, with demonstrated improvements in sample efficiency in numerical experiments.
Abstract
We study the problem of learning to stabilize (LTS) a linear time-invariant (LTI) system. Policy gradient (PG) methods for control assume access to an initial stabilizing policy. However, designing such a policy for an unknown system is one of the most fundamental problems in control, and it may be as hard as learning the optimal policy itself. Existing work on the LTS problem requires large data as it scales quadratically with the ambient dimension. We propose a two-phase approach that first learns the left unstable subspace of the system and then solves a series of discounted linear quadratic regulator (LQR) problems on the learned unstable subspace, targeting to stabilize only the system's unstable dynamics and reduce the effective dimension of the control space. We provide non-asymptotic guarantees for both phases and demonstrate that operating on the unstable subspace reduces sample complexity. In particular, when the number of unstable modes is much smaller than the state dimension, our analysis reveals that LTS on the unstable subspace substantially speeds up the stabilization process. Numerical experiments are provided to support this sample complexity reduction achieved by our approach.
