Assessing the Quality of a Set of Basis Functions for Inverse Optimal Control via Projection onto Global Minimizers
Filip Bečanović, Jared Miller, Vincent Bonnet, Kosta Jovanović, Samer Mohammed
TL;DR
This work tackles inverse optimal control by evaluating how well a convex-basis set $\mathcal{F}$ expresses a given observation $y$ through projection onto the set of global minimizers $\mathcal{G}$ (unconstrained) or $\mathcal{G}^c$ (constrained). It develops a geometric understanding of these optima sets for convex and quadratic bases, and provides practical upper bounds via bilevel gradient methods and lower bounds via LMIs, enabling validation or rejection of candidate bases. The authors demonstrate the method with toy and small-scale examples, showing how LMI-based bounds can robustly detect incompatible bases and offer quantitative distance measures, while also providing a path to extensions such as piecewise costs, local minimizers, and polynomial optimization. The framework thus offers a principled, computationally tractable approach to assessing IOC bases with potential applicability to biomechanics, robotics, and related domains.
Abstract
Inverse optimization (Inverse optimal control) is the task of imputing a cost function such that given test points (trajectories) are (nearly) optimal with respect to the discovered cost. Prior methods in inverse optimization assume that the true cost is a convex combination of a set of convex basis functions and that this basis is consistent with the test points. However, the consistency assumption is not always justified, as in many applications the principles by which the data is generated are not well understood. This work proposes using the distance between a test point and the set of global optima generated by the convex combinations of the convex basis functions as a measurement for the expressive quality of the basis with respect to the test point. A large minimal distance invalidates the set of basis functions. The concept of a set of global optima is introduced and its properties are explored in unconstrained and constrained settings. Upper and lower bounds for the minimum distance in the convex quadratic setting are implemented by bi-level gradient descent and an enriched linear matrix inequality respectively. Extensions to this framework include max-representable basis functions, nonconvex basis functions (local minima), and applying polynomial optimization techniques.
