Minimal upper bounds in the löwner order and application to invariant computation for switched systems
Majorants minimaux dans l'ordre de löwner et application au calcul d'invariants de systèmes commutés
Résumé
The computation of invariant sets for dynamical systems is a crucialelement of program verification and control theory, as such setscertify the absence of unwanted behaviours. We consider in particularswitched systems, for which the computation of invariants is alreadydifficult. Recently, several approaches based on optimizationtechniques such as semi-definite programming have been appliedsuccessfully to compute piecewise quadratic invariants of switchedsystems. However, their high computational cost becomes prohibitive onlarge instances.In this thesis, we develop a new class of algorithms to computepiecewise quadratic invariants. These algorithms rely on geometrical and metricproperties of the space of positive semidefinite matrices equippedwith the Löwner order.First, we characterize minimal upper bounds in this order. We show inparticular that the set of minimal upper bounds of two matrices can beidentified to a quotient of an indefinite orthogonal group, providinga ``quantitative'' refinement of a theorem of Kadison. More generally,we characterize minimal upper bounds with respect to a cone ordering,and show that for a wide family of cones, there is a canonicalselection of a minimal upper bound, defined in terms of the generatingfunction of the cones. This extends the construction of the Löwnerellipsoid. In the case of the cone of positive semidefinite matrices,we show that this canonical selection satifies several matrixinequalities, and we estimate its Lispchitz constant with respect toconvenient invariant metrics defined on the interior of the cone(Riemannian metric, Thompson metric).Then, we apply these results to the computation of piecewise-quadraticinvariants. We formulate the latter as a non-linear fixed pointproblem over a product of spaces of positive definite matrices. Thisproblem involves an operator which may be thought of as the tropicalanalogue of the Kraus maps (quantum channels) arising in quantuminformation theory. This leads to a class of fast iterative numericalschemes, avoiding the recourse to linear matrix inequalities, which weshow to converge, under some restrictions. We implemented thisapproach, by developing the tool MEGA (``Minimal ellipsoid geometricanalyser''), and report experimental results, on switched linear andaffine systems, showing an improvement of several orders of magnitudein terms of scalability, on some instances (with e.g., approximationsof the joint spectral radius in dimension 500). We also applied thismethod to the approximation of the value function of switched linearquadratic optimal control problems.
Le calcul d'ensembles invariants est un élément crucial en vérification de programme et en théorie du contrôle, car de tels ensembles certifient l'absence de comportements indésirables. Nous étudions en particulier les systèmes commutés, pour lesquel le calcul d'un ensemble invariant est déjà difficile. Plusieurs approches récentes utilisant des techniques d'optimisation telles la programmation semidéfinie ont été appliquées avec succès au calcul d'invariants quadratiques par morceaux pour des systèmes commutés. En revanche, ces méthodes ne sont pas utilisables en grande dimension car elles nécessitent trop de ressources informatiques.Nous développons dans cette thèse une nouvelle classe d'algorithme pour calculer des invariants quadratiques par morceaux. Ces algorithmes reposent sur les propriétés géométriques et métriques de l'espace des matrices positive semidéfinies équipées de l'ordre de Löwner.Tout d'abord, nous caractérisons l'ensemble des majorants minimaux dans cet ordre. Nous montrons que l'ensemble des majorants minimaux de deux matrices s'identifie au quotient d'un groupe orthogonal indéfini, donnant ainsi un raffinement "quantitatif" d'un théorème de Kadison. Plus généralement, nous caractérisons les majorants minimaux dans un ordre défini par un cône et nous prouvons qu'il existe pour une grande famille de cônes une sélection de majorant minimal canonique, définie à partir des fonctions génératrices de ces cônes. Ceci généralise la définition de l'ellipsoïde de Löwner. dans le cas du cône des matrices positives semidéfinies, nous montrons que cette sélection canonique satisfait plusieurs inégalités matricielles et nous donnons des estimations de sa constante de Lipschitz par rapport à plusieurs métriques convenables définies sur l'intérieur du cône (métrique Riemannienne, métrique de Thompson).Nous appliquons ensuite ces résultats au calcul d'invariants quadratiques par morceaux. Nous formulons ce dernier comme un problème de point fixe non-linéaire sur un produit de cônes de matrices positives semidéfinies. Ce problème fait intervenir un opérateur qui peut s'interpréter comme l'analogue tropical d'une application de Kraus (un canal quantique) qui apparaît en théorie d'information quantique. Nous obtenons ainsi une classe de schémas itératifs rapides, n'utilisant les inégalités linéaires matricielles, dont nous prouvons la convergence sous quelques restrictions. Nous avons implémenté cette approche en développant l'outil MEGA (``Minimal ellipsoid geometricanalyser''). Nos résultats expérimentaux démontrent une amélioration de l'ordre de quelques ordres de grandeur en termes de scalabilité (par exemple, approximation du "rayon spectral joint" en dimension 500). Nous avons aussi appliqué cette méthode à l'approximation de la fonction valeur d'un problème de contrôle optimal commutant entre des modèles linéaires-quadratiques.
Origine : Version validée par le jury (STAR)