Shortest-path recovery from signature with an optimal control approach
Marco Rauscher, Alessandro Scagliotti, Felipe Pagginelli Patricio
TL;DR
This work treats the problem of reconstructing the shortest path consistent with a given signature as a sub-Riemannian geodesic problem on the signature manifold $G^N(\mathbb{R}^d)$. The authors embed the problem into optimal control with dynamics $\dot{\xi}=\sum_{i=1}^d a_i U_i(\xi)$ and connect path-length to the control energy, establishing a $\Gamma$-convergence bridge from a soft-endpoint formulation to the hard-endpoint geodesic problem. They derive a Pontryagin Maximum Principle framework for signatures and implement a Sakawa-type iterative scheme, including variants with final-time optimization and with additional bracket vector fields, to numerically compute length-minimizing signatures and recover the associated $H^1$ paths. Numerical experiments validate the approach, showing effective recovery of shortest-signature-preserving paths and revealing the practical sparsity of higher-order bracket directions. The results provide a scalable, geometry-driven method for inverse signature problems with potential applications in stochastic-process path inference and signature-based learning, without restrictive assumptions on signature order or path dimension.
Abstract
In this paper, we consider the signature-to-path reconstruction problem from the control theoretic perspective. Namely, we design an optimal control problem whose solution leads to the minimal-length path that generates a given signature. In order to do that, we minimize a cost functional consisting of two competing terms, i.e., a weighted final-time cost combined with the $L^2$-norm squared of the controls. Moreover, we can show that, by taking the limit to infinity of the parameter that tunes the final-time cost, the problem $Γ$-converges to the problem of finding a sub-Riemannian geodesic connecting two signatures. Finally, we provide an alternative reformulation of the latter problem, which is particularly suitable for the numerical implementation.
