⚡ electrical engineering

A piecewise ellipsoidal reachable set estimation method for continuous bimodal piecewise affine systems

Cet article propose une méthode simplifiée d'estimation de l'ensemble atteignable ellipsoïdal par morceaux pour les systèmes affines par morceaux bimodaux continus en utilisant des fonctions de Lyapunov quadratiques par morceaux et en dérivant des caractérisations par matrices linéaires qui exploitent la structure dynamique spécifique du système.

Auteurs originaux : Thuan Le Quang, Nam Phan Thanh, Simone Baldi

Publié 2026-07-09

📖 6 min de lecture🧠 Analyse approfondie

CC BY 4.0

Auteurs originaux : Thuan Le Quang, Nam Phan Thanh, Simone Baldi

Article original sous licence CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ Ceci est une explication générée par l'IA de l'article ci-dessous. Elle n'a pas été rédigée ni approuvée par les auteurs. Pour une précision technique, consultez l'article original. Lire la clause de non-responsabilité complète

Imaginez que vous essayez de prédire exactement jusqu'où une balle rebondissante pourrait voler si vous la lanciez dans une pièce remplie de murs invisibles et élastiques. Cette balle ne se contente pas de rebondir sur les murs ; elle est aussi poussée par un vent changeant et imprévisible. Dans le monde de l'ingénierie, cette « balle » est une machine ou un robot, les « murs » sont des règles qui changent selon l'endroit où se trouve la machine, et le « vent » est un bruit ou des perturbations aléatoires.

La grande question que les ingénieurs posent est la suivante : « Quelle est la plus grande zone dans laquelle cette machine pourrait éventuellement errer ? » Cette zone est appelée l'ensemble atteignable (reachable set). Si vous pouvez dessiner une bulle invisible et serrée autour de cette zone, vous savez que la machine ne sera jamais capable de s'échapper de cette bulle, peu importe la violence du vent. Cela est crucial pour la sécurité : si votre robot est censé rester dans une zone sûre, vous devez savoir exactement quelle taille doit avoir cette zone.

Le problème des anciennes cartes

Pendant longtemps, essayer de dessiner cette « bulle de sécurité » pour des machines qui passent entre deux comportements différents (appelées systèmes affines par morceaux bimodaux) revenait à cartographier un labyrinthe avec une masse. Les méthodes existantes étaient souvent trop complexes, comme essayer de résoudre un puzzle géant en 3D alors qu'un simple croquis en 2D suffirait. Certaines méthodes étaient si lourdes qu'elles prenaient des heures à calculer, et d'autres rendaient la bulle de sécurité si grande et floue qu'elle n'était plus très utile.

Les auteurs de cet article, Le Quang Thuan, Phan Thanh Nam et Simone Baldi, ont décidé de construire un outil meilleur et plus léger. Ils ne voulaient pas seulement un gros nuage flou ; ils voulaient une bulle serrée et de forme personnalisée.

Le nouvel outil : Un costume élastique en deux pièces

La principale découverte de l'équipe est une nouvelle façon de calculer cette bulnes de sécurité en utilisant ce qu'ils appellent une fonction de Lyapunov quadratique par morceaux.

Pour comprendre cela, imaginez que la bulle de sécurité n'est pas une sphère unique et rigide. Au lieu de cela, c'est comme un costume élastique en deux pièces.

Sur le côté gauche de la pièce, le costume est fait d'un type de caoutchouc extensible (défini par une forme mathématique).
Sur le côté droit, le costume est fait d'un deuxième type de caoutchouc légèrement différent (défini par une seconde forme).
La magie opère précisément à la ligne où les deux côtés se rejoignent. Les auteurs ont trouvé comment recoudre ces deux pièces si parfaitement que le costume est lisse et continu — sans écarts, sans déchirures, sans bosses bizarres.

Ils ont prouvé mathématiquement que si vous recousez ces deux pièces correctement, le « costume » résultant restera toujours serré autour des trajectoires possibles de la machine. Ils ont utilisé un ensemble de règles appelées Inégalités Matricielles Linéaires (LMI) pour trouver la forme parfaite pour chaque pièce du costume. Voyez ces règles comme une recette qui dit à un ordinateur exactement comment étirer le caoutchouc pour qu'il s'adapte parfaitement aux mouvements de la machine.

Ce qu'ils ont exclu

L'article est très clair sur ce que leur méthode n'est pas.

Ils n'utilisent pas explicitement les anciennes méthodes trop complexes conçues pour les machines possédant de nombreux modes de commutation différents (systèmes multimodaux). Ils soutiennent que pour les machines qui ne basculent qu'entre deux modes (systèmes bimodaux), ces anciennes méthodes reviennent à utiliser un supercalculateur pour résoudre un Sudoku — c'est trop de travail pour la tâche.
Ils ne supposent pas non plus que le « vent » (les perturbations) finit par s'arrêter de souffler. Certaines méthodes plus anciennes exigeaient que le vent tombe à zéro pour que les calculs fonctionnent. Cette nouvelle méthode fonctionne même si le vent continue de souffler éternellement, tant qu'il reste dans une certaine limite de force.
Ils ne reposent pas sur une forme « commune » unique pour toute la pièce. Dans leurs tests numériques, ils ont montré qu'utiliser une seule forme (une fonction quadratique commune) crée une bulle beaucoup plus grande et lâche qui gaspille de l'espace. Leur costume en deux pièces est bien plus ajusté.

À quel point en sont-ils sûrs ?

Les auteurs n'ont pas seulement deviné ; ils ont prouvé que leur méthode fonctionne grâce à des mathématiques rigoureuses. Ils ont dérivé des conditions spécifiques (les LMI mentionnées plus haut) qui garantissent que la bulle de sécurité tiendra bon.

Pour démontrer l'efficacité de leur approche, ils ont mené des simulations (expériences informatiques) sur deux scénarios différents :

Un exemple mathématique simple : Ils ont comparé leur nouveau costume en deux pièces à l'ancien costume en une seule pièce. Le résultat ? Le nouveau costume était beaucoup plus serré. Il a capturé la nature « orthogonale » (presque à angle droit) du mouvement de la machine, ce que l'ancien costume avait manqué.
Un système de chariot mécanique : Ils ont simulé un système mécanique réel composé de deux chariots et d'un ressort. Ici, la différence était énorme.
- Leur nouvelle méthode a calculé la bulle de sécurité en environ 0,1 seconde.
- Ils l'ont comparée à un outil industriel populaire appelé CORA, qui tentait de faire le même travail en utilisant une technique différente (les zonotopes). CORA a mis plus de 200 secondes (plus de 3 minutes) pour obtenir un résultat, et même ainsi, rendre le résultat plus précis ne faisait que rallonger le temps de calcul sans réellement améliorer la taille de la bulle.

Les auteurs affirment que pour ce système d'ordre 4, leur méthode est environ 3 ordres de grandeur plus rapide que l'outil existant. Ils soulignent que, bien que leur méthode soit une amélioration significative pour les systèmes bimodaux (deux modes), elle n'est pas encore une baguette magique pour les systèmes à trois modes ou plus ; cela relève de la recherche future.

Ce qu'il faut retenir

En résumé, cet article propose une façon plus intelligente, plus rapide et plus précise de dessiner la « bulle de sécurité » autour des machines qui basculent entre deux comportements. En traitant la zone de sécurité comme deux pièces élastiques cousues ensemble plutôt que comme une forme rigide unique, les auteurs ont créé une méthode mathématiquement solide et d'une rapidité informatique fulgurante. C'est un rappel que, parfois, la meilleure façon de résoudre un problème complexe n'est pas de construire une machine plus grande, mais de coudre un costume mieux ajusté.

Résumé technique : Une méthode d'estimation de l'ensemble atteignable ellipsoïdal par morceaux pour les systèmes affines par morceaux bimodaux continus

Énoncé du problème
L'article traite du problème de l'estimation des ensembles atteignables pour les systèmes affines par morceaux (PWA) bimodaux continus soumis à des perturbations exogènes. Bien que les systèmes PWA soient largement utilisés pour modéliser des phénomènes physiques (ex. : relais, hystérésis, saturation) et approximer des dynamiques non linéaires, le calcul de leurs ensembles atteignables est généralement un problème indécidable. Les méthodes existantes pour les systèmes hybrides PWA généraux souffrent souvent d'une complexité de calcul élevée ou sont mal adaptées à la structure dynamique spécifique des systèmes bimodaux. De plus, bien que des techniques de bornage ellipsoïdal existent pour les systèmes linéaires et commutés, elles reposent souvent sur des fonctions de Lyapunov communes ou nécessitent des décompositions matricielles spécifiques qui peuvent conduire à des estimations trop conservatrices ou échouer à gérer efficacement le basculement dépendant de l'état avec des perturbations non nulles.

Méthodologie
Les auteurs proposent une nouvelle approche dans le cadre du bornage ellipsoïdal, utilisant des fonctions de Lyapunov quadratiques par morceaux adaptées à la nature bimodale du système. Le système est défini par deux dynamiques affines commutant selon un hyperplan $c^T x + f = 0$ . Les perturbations sont contraintes à l'intérieur d'une ellipsoïde.

La méthodologie suit ces étapes techniques :

Conditions de positivité de l'espace moitié : Les auteurs dérivent d'abord des conditions (Lemme 3.1) pour garantir qu'une fonction quadratique est positive sur un espace moitié spécifique. Ceci est réalisé en caractérisant l'espace moitié à l'aide d'une transformation matricielle spécifique impliquant les paramètres de la surface de commutation.
Contraintes de continuité : Pour assurer la continuité de la fonction de Lyapunov par morceaux à travers l'hyperplan de commutation, les auteurs établissent des conditions (Lemme 3.2) reliant les paramètres des fonctions quadratiques des deux modes.
Homogénéisation et procédure S : Les auteurs emploient le concept d'homogénéisation d'ensembles et une procédure S généralisée (Lemme 3.5) pour gérer simultanément les contraintes imposées par l'ellipsoïde de perturbation et la surface de commutation.
Formulation LMI : En combinant ces lemmes avec une caractérisation de l'ensemble atteignable par une fonction de Lyapunov (Lemme 2.3), ils dérivent des conditions d'inégalités matricielles linéaires (LMI) (Théorème 3.6 et Corollaire 3.7). Ces conditions garantissent que la dérivée de la fonction de Lyapunov par morceaux est définie négativement en dehors d'une région spécifique, encadrant ainsi l'ensemble atteignable dans un ensemble ellipsoïdal par morceaux.

Contributions clés

Caractérisations LMI novatrices : L'article introduit de nouvelles conditions LMI (Lemmes 3.1, 3.2, 3.4 et Théorème 3.5) spécifiquement conçues pour exploiter la dynamique bimodale. Contrairement aux méthodes précédentes qui supposent une matrice de Lyapunov commune ou nécessitent des décompositions spéciales (ex. : $P_i = F_i^T T F_i$ ), cette approche construit l'ellipsoïde par morceaux sans de telles hypothèses restrictives.
Gestion des perturbations non nulles : La méthode accepte des perturbations bornées et non nulles, une caractéristique où une partie de la littérature existante (ex. : [28]) exige que les perturbations tendent vers zéro à l'infini.
Commutation dépendante de l'état : L'approche est formulée pour une commutation dépendante de l'état, se distinguant ainsi des méthodes reposant sur une commutation dépendante du temps.
Efficacité de calcul : La méthode proposée aboutit à un problème d'optimisation convexe qui peut être résolu efficacement, offrant une réduction significative du temps de calcul par rapport aux approches basées sur les zonotopes pour les systèmes d'ordre supérieur.

Résultats
La méthode proposée est validée par deux exemples numériques :

Exemple 4.1 (Système 2D) : Un système bimodal avec des plans de phase « orthogonaux » est analysé. Les résultats démontrent que la fonction de Lyapunov quadratique par morceaux fournit une borne plus étroite (moins conservatrice) sur l'ensemble atteignable par rapport à une fonction de Lyapunov quadratique commune. La méthode par morceaux capture avec succès la géométrie du système que la méthode commune manque.
Exemple 4.2 (Système mécanique 4D) : Un système mécanique avec un ressort unidirectionnel est utilisé pour comparer la méthode proposée à la boîte à outils CORA (qui utilise des zonotopes).
- Précision : L'estimation ellipsoïdale par morceaux est montrée comme étant plus petite (plus étroite) que l'estimation ellipsoïdale commune.
- Performance : La méthode proposée résout le système d'ordre 4 en environ 0,1 seconde. En revanche, la boîte à outils CORA a nécessité plus de 200 secondes pour un zonotope d'ordre 40 et plus de 215 secondes pour un zonotope d'ordre 50. Notamment, l'augmentation de l'ordre du zonotope dans CORA augmentait le temps de calcul sans améliorer significativement l'estimation, tandis que la méthode proposée a atteint une haute précision avec un coût de calcul minimal.

Signification et affirmations
L'article affirme que sa principale importance réside dans l'exploitation de la structure dynamique spécifique des systèmes PWA bimodaux continus pour proposer une approche plus simple et moins conservatrice que les méthodes générales pour systèmes hybrides. Les auteurs soulignent que leurs conditions sont novatrices car elles évitent le besoin de décompositions spéciales des matrices de Lyapunov, qui sont souvent difficiles à dériver de manière optimale dans la littérature existante.

Les auteurs notent modestement que, bien que la méthode soit efficace pour les systèmes PWA bimodaux continus, l'étendre à des systèmes PWA multimodaux généraux (trimodaux ou supérieurs) est non trivial et nécessite des recherches futures, impliquant potentiellement des représentations spéciales comme les hyperplans articulés (hinging hyperplanes). Ce travail positionne la méthode comme un outil spécialisé et efficace pour une classe de systèmes incluant d'importantes applications d'ingénierie comme les systèmes de relais linéaires.

Le problème des anciennes cartes

Le nouvel outil : Un costume élastique en deux pièces

Ce qu'ils ont exclu

À quel point en sont-ils sûrs ?

Ce qu'il faut retenir

Articles similaires