On Koopman Resolvents and Frequency Response of Nonlinear Systems

Ce papier propose une nouvelle formulation de la réponse fréquentielle des systèmes non linéaires dans le cadre de l'opérateur de Koopman, permettant de définir des diagrammes de Bode via la transformée de Laplace de la sortie et en s'appuyant sur la théorie des résolvantes.

L'essentiel

Imaginez que vous essayez de comprendre comment réagit une voiture de course très complexe (un système non linéaire) lorsque vous appuyez sur l'accélérateur de manière rythmée. Pour les voitures simples (les systèmes linéaires), les ingénieurs ont une boîte à outils magique appelée « réponse fréquentielle ». C'est comme un manuel qui vous dit exactement comment la voiture va vibrer ou accélérer en fonction de la fréquence à laquelle vous tapez sur l'accélérateur.

Mais pour les voitures complexes, ce manuel n'existe pas vraiment. Les mathématiques deviennent un labyrinthe effrayant.

C'est là que cette paper (article) intervient. Elle propose une nouvelle façon de voir les choses, en utilisant une idée mathématique appelée l'opérateur de Koopman. Voici l'explication simple, avec des analogies pour tout le monde.

1. Le Problème : Le Chaos vs. La Ligne Droite

Dans le monde réel, la plupart des systèmes sont « non linéaires ». Cela signifie que si vous doublez l'effort, le résultat n'est pas forcément doublé. C'est comme essayer de prédire la météo : un petit changement peut tout faire basculer.
Les ingénieurs adorent les systèmes « linéaires » (droits, prévisibles) car ils peuvent les analyser avec des graphiques simples appelés diagrammes de Bode. Ces graphiques montrent deux choses :

• Le Gain : À quel point le système amplifie le signal (le volume).
• La Phase : Le décalage dans le temps (le retard).

Pour les systèmes complexes, on ne savait pas vraiment comment tracer ces graphiques de manière rigoureuse.

2. La Solution Magique : Le « Projecteur » de Koopman

Les auteurs (Susuki, Katayama, Mauroy, Mezić) disent : « Et si on ne regardait pas le système lui-même, mais une projection de celui-ci ? »

Imaginez que vous avez une marionnette complexe dont les fils s'emmêlent (le système non linéaire). C'est dur à analyser.
L'opérateur de Koopman, c'est comme un projecteur de cinéma très puissant. Il projette le mouvement de cette marionnette sur un grand écran blanc. Sur l'écran, au lieu de voir des fils emmêlés, vous voyez une danse de lignes droites et de cercles parfaits (un système linéaire, mais dans un monde de dimension infinie).

Même si la marionnette fait des mouvements chaotiques, son ombre sur l'écran suit des règles mathématiques simples et prévisibles.

3. La Nouvelle Méthode : La « Résonance »

Le papier propose d'utiliser ce projecteur pour créer un nouveau type de « réponse fréquentielle » pour les systèmes complexes.

Voici comment ils font, étape par étape :

• L'Input (L'entrée) : Ils envoient un signal rythmé (comme un battement de cœur ou une note de musique) dans le système.
• Le Skew-Product (Le couple) : Ils imaginent que le système et le signal rythmé forment une seule équipe. C'est comme si on attachait le battement de cœur directement à la marionnette pour étudier leur danse commune.
• Le Résolvant (Le détecteur) : Ils utilisent un outil mathématique appelé « résolvant de Koopman ». Imaginez-le comme un tuning fork (diapason) géant. Quand vous le faites vibrer à une fréquence précise, il ne résonne que si le système contient cette fréquence cachée.
• Le Résultat (Les Modes) : Si le système résonne, le « diapason » révèle un nombre magique (une valeur complexe). Ce nombre contient exactement le Gain (l'amplitude) et la Phase (le retard) que nous cherchions.

4. Pourquoi c'est génial ?

Grâce à cette méthode, on peut enfin dessiner les fameux diagrammes de Bode pour des systèmes complexes (comme des réacteurs chimiques, des réseaux électriques ou des systèmes biologiques).

• Avant : On disait « C'est trop compliqué, on ne peut pas prédire la réponse à une fréquence précise. »
• Maintenant : On dit « Regardez, si vous tapez à cette fréquence, le système va réagir avec ce gain et ce retard, exactement comme le montre notre graphique. »

5. Les Conditions (Quand ça marche ?)

Les auteurs précisent que cette magie ne fonctionne pas dans tous les cas. Il faut que le système soit « stable » (qu'il ne s'effondre pas) et qu'il ait un comportement périodique (qu'il finisse par se caler sur un rythme).
Ils donnent trois scénarios où cela fonctionne parfaitement :

1. Les systèmes déjà simples (linéaires).
2. Les systèmes qui finissent toujours par se stabiliser (comme une balançoire qui finit par s'arrêter ou osciller régulièrement).
3. Les systèmes « ergodiques » (qui explorent tout leur espace de manière équitable, comme un gaz dans une boîte).

En résumé

Cette paper est comme un traducteur universel. Elle prend le langage chaotique des systèmes complexes et le traduit en un langage simple et linéaire que les ingénieurs comprennent parfaitement. Elle permet de voir les « couleurs » (fréquences) cachées dans le chaos et de prédire comment un système complexe va réagir à la musique du monde réel.

C'est une avancée majeure pour pouvoir concevoir et contrôler des systèmes complexes avec la même précision que l'on contrôle une simple machine à café.

Résumé technique

Voici un résumé technique détaillé de l'article « On Koopman Resolvents and Frequency Response of Nonlinear Systems » en français.

Titre : Sur les Résolvantes de Koopman et la Réponse en Fréquence des Systèmes Non Linéaires

Auteurs : Yoshihiko Susuki, Natsuki Katayama, Alexandre Mauroy et Igor Mezić.

---

1. Problématique

La réponse en fréquence est un concept fondamental en ingénierie des commandes pour l'analyse et la synthèse des systèmes linéaires invariants dans le temps (LTI). Elle est définie comme une fonction complexe de la fréquence d'excitation, permettant la construction de diagrammes de Bode (gain et phase).

Cependant, pour les systèmes non linéaires, la définition d'une réponse en fréquence rigoureuse reste un défi. Les approches existantes (méthode de l'équilibre harmonique, fonction descriptive, séries de Volterra) sont principalement basées sur des formulations temporelles ou d'espace d'état, et ne s'inscrivent pas naturellement dans une théorie fréquentielle unifiée comparable à celle des systèmes LTI. L'objectif de cet article est de combler ce vide en proposant une formulation de la réponse en fréquence pour les systèmes non linéaires, guidée par la théorie rigoureuse des opérateurs de Koopman.

2. Méthodologie

Les auteurs utilisent le cadre des opérateurs de Koopman, qui permet d'analyser des dynamiques non linéaires via des opérateurs linéaires (mais de dimension infinie) agissant sur un espace d'observables.

La méthodologie repose sur les étapes clés suivantes :

1. Forme Skew-Produit (Produit Croisé) :
   Pour traiter un système non linéaire soumis à une entrée périodique $u(t) = u_0 e^{i\omega t}$, les auteurs transforment le système non autonome en un système autonome équivalent de dimension supérieure. Ce système augmenté (forme skew-product) inclut la dynamique de l'entrée comme une variable d'état supplémentaire :
   $$\dot{x} = F(x, u), \quad \dot{u} = i\omega u$$
   Cela permet d'appliquer la théorie des opérateurs de Koopman à un système autonome.

2. Utilisation de la Résolvante de Koopman :
   La réponse du système est analysée via la transformée de Laplace de la sortie $y(t)$. Dans le cadre de Koopman, cette transformée est représentée par l'action de la résolvante de Koopman $R(s; \mathcal{L}_{forced}) = (sI - \mathcal{L}_{forced})^{-1}$ sur l'observable $g$, où $\mathcal{L}_{forced}$ est le générateur infinitésimal du système skew-product.

3. Définition de la Réponse en Fréquence :
   Les auteurs définissent la réponse en fréquence non pas comme un simple rapport entrée/sortie, mais comme les modes de Koopman associés aux pôles de la résolvante.

   • Cas Fondamental et Harmonique : Pour une fréquence de sortie $n\omega$, la réponse $H_n(\omega; g)$ est liée au résidu de la résolvante au pôle $s = in\omega$.
   • Cas Subharmonique : Pour une fréquence $\omega/n$, la réponse $H_{1/n}(\omega; g)$ est liée au pôle $s = i\omega/n$.

4. Conditions d'Existence :
   L'article établit des conditions suffisantes pour l'existence de ces réponses (pôles simples et isolés de la résolvante) pour trois classes de dynamiques :

   • Dynamiques LTI (cas de référence).
   • Dynamiques globalement stables (existence d'un cycle limite stable).
   • Dynamiques ergodiques (sur un attracteur compact).

3. Contributions Clés

• Nouvelle Formulation Théorique : Proposition d'une définition de la réponse en fréquence pour les systèmes non linéaires basée sur la théorie spectrale des opérateurs de Koopman. Cette approche généralise la méthode classique des systèmes LTI.
• Lien avec les Modes de Koopman : Démonstration que la réponse en fréquence d'un système non linéaire correspond directement aux modes de Koopman (les coefficients d'amplitude associés aux fonctions propres) du système skew-product.
• Visualisation Fréquentielle : La réponse étant une fonction complexe de la fréquence, elle permet de tracer des diagrammes de Bode (gain et phase) pour des systèmes non linéaires, offrant une intuition physique similaire à celle utilisée en automatique linéaire.
• Estimation Numérique : Puisque la réponse est un mode de Koopman, elle peut être estimée numériquement à l'aide de méthodes établies comme la Décomposition en Modes Dynamiques (DMD).

4. Résultats Principaux

• Théorèmes 1 et 2 : Ils établissent des formules explicites pour calculer les réponses en fréquence $H_n$ et $H_{1/n}$ en utilisant le résidu de la transformée de Laplace de la sortie autour des pôles $in\omega$ et $i\omega/n$.
  $$H_n(\omega; g) = u_0^{-n} \lim_{s \to in\omega} (s - in\omega) [R(s; \mathcal{L}_{forced})g](x_0, u_0)$$
• Exemples Résolus :
  • Cas Linéaire 1D : La méthode redonne exactement la fonction de transfert classique $b/(i\omega - a)$.
  • Cas Non Linéaire 2D : Pour un système avec un terme quadratique ($x_2^2$), la méthode montre que la sortie $x_1$ ne contient que la composante à $2\omega$(réponse du second ordre), tandis que$x_2$contient la composante à$\omega$. Les diagrammes de Bode générés confirment le comportement de type "retard de 3ème ordre" pour la composante non linéaire.
• Validité des Conditions : L'article prouve que pour les systèmes LTI, les systèmes à cycle limite stable (avec champs analytiques) et les systèmes ergodiques, les conditions nécessaires pour appliquer ces théorèmes (existence de pôles simples isolés) sont satisfaites.

5. Signification et Perspectives

Cet article représente une avancée majeure en reliant l'analyse fréquentielle classique à la théorie moderne des opérateurs de Koopman.

• Unification : Il offre un cadre unifié pour analyser les systèmes linéaires et non linéaires dans le domaine fréquentiel, dépassant les approches purement temporelles ou phénoménologiques.
• Outils de Synthèse : En permettant de visualiser gain et phase pour des systèmes non linéaires, cette méthode ouvre la voie à de nouvelles stratégies de synthèse de contrôleurs non linéaires basées sur des critères fréquentiels.
• Futur : Les auteurs suggèrent d'explorer la continuité et la différentiabilité de ces réponses, ainsi que leur application à l'analyse de stabilité et de passivité des systèmes non linéaires. Ils prévoient également d'étendre la théorie aux spectres continus (cas des mélanges faibles).

En résumé, ce travail transforme la réponse en fréquence d'un concept purement linéaire en un outil puissant et général pour l'ingénierie des systèmes non linéaires, en s'appuyant sur la puissance de la linéarisation infinie offerte par l'opérateur de Koopman.