Hill-Type Stability Analysis of Periodic Solutions of Fractional-Order Differential Equations

Résumé technique : Analyse de stabilité de type Hill des solutions périodiques d'équations différentielles d'ordre fractionnaire

Énoncé du problème
L'article traite de l'analyse de stabilité des solutions périodiques dans les équations différentielles non linéaires d'ordre fractionnaire (FODE). Une limitation fondamentale des FODE de type Caputo classiques est qu'elles n'admettent pas de solutions strictement périodiques ; les fonctions périodiques sous les dérivées de Caputo ne sont qu'asymptotiquement périodiques. Par conséquent, les outils standards d'analyse de stabilité, tels que la théorie de Floquet et la méthode de la matrice de Hill, ne peuvent pas être appliqués directement à la dynamique des perturbations linéarisées des solutions périodiques dans les systèmes de Caputo. Bien que les FODE de type Liouville-Weyl (définies avec une borne inférieure infinie, $-\infty$) admettent des solutions strictement périodiques, une théorie de Floquet rigoureuse pour ces systèmes n'a pas encore été établie. Les tentatives précédentes d'étendre la théorie de Floquet aux systèmes d'ordre fractionnaire se sont appuyées sur des approximations (par exemple, en supposant des identités multiplicatives pour les fonctions de Mittag-Leffler) ou ont échoué à caractériser de manière exhaustive les frontières de stabilité et les types de bifurcation.

Méthodologie
Les auteurs proposent un cadre basé sur les FODE de type Liouville-Weyl pour étudier les solutions périodiques. La méthodologie centrale comprend :

1. Formulation du système : Définition d'une classe d'équations différentielles non linéaires de type Liouville-Weyl avec une borne inférieure infinie ($-\infty$). Le problème aux valeurs initiales est formulé avec une fonction d'initialisation $x_0(t)$ définie sur $(-\infty, t_0]$, appartenant à un espace fonctionnel spécifique $X(-\infty, t_0]$ qui garantit l'existence du terme de forçage dérivé de l'historique du système.
2. Linéarisation : La stabilité d'une solution périodique $x_p(t)$ est analysée en linéarisant le système autour de celle-ci, aboutissant à une équation différentielle linéaire d'ordre fractionnaire à coefficients périodiques (LTP-FODE) avec une matrice jacobienne périodique $J(t)$.
3. Propriétés de la dérivée de Liouville-Weyl : En utilisant la propriété selon laquelle la dérivée de Liouville-Weyl d'une fonction exponentielle $e^{st}$ est $s^\alpha e^{st}$ (pour $\text{Re}(s) \ge 0$), les auteurs tentent d'étendre l'ansatz de Floquet $y(t) = e^{\lambda t} p(t)$ aux systèmes d'ordre fractionnaire.
4. Matrice de Hill fractionnaire : En substituant l'ansatz de Floquet dans la LTP-FODE et en utilisant des développements en série de Fourier pour les coefficients périodiques, les auteurs dérivent un problème de valeurs propres non linéaire de dimension infinie. Cela conduit à la définition de la « matrice de Hill fractionnaire », $H^{(\alpha)}_\infty(\lambda)$, où les blocs diagonaux impliquent des termes de la forme $(\lambda + ik\omega)^\alpha I - J_k$.
5. Analyse des sous-systèmes linéaires invariants dans le temps (LTI) : Avant d'aborder le cas périodique, les auteurs analysent les LTP-FODE avec mémoire infinie. Ils comparent les solutions du système de Liouville-Weyl avec celles du système de Caputo standard (borne inférieure nulle) pour comprendre le comportement des solutions décroissantes.

Contributions et résultats clés

• Cadre pour les solutions périodiques : L'article établit une formulation rigoureuse du problème aux valeurs initiales pour les FODE de type Liouville-Weyl qui admet des solutions strictement périodiques, surmontant ainsi la limitation des dérivées de Caputo.
• Extension de la théorie de Floquet : Les auteurs dérivent une approche par matrice de Hill fractionnaire pour les LTP-FODE. Ils prouvent que si une solution existe sous la forme $y(t) = e^{\lambda t} p(t)$ avec $\text{Re}(\lambda) \ge 0$, elle doit satisfaire le problème de valeurs propres non linéaire $H^{(\alpha)}_\infty(\lambda)p_\infty = 0$.
• Limitation concernant les solutions décroissantes : Une découverte centrale est que la méthode de Hill fractionnaire proposée (et l'ansatz exponentiel en général) ne peut évaluer que les solutions à croissance exponentielle (instabilité). Elle échoue à capturer les solutions à décroissance algébrique.
  • Dans le cas LTI, les auteurs montrent que, tandis que les solutions divergentes correspondent à des valeurs propres à partie réelle positive, les solutions convergentes (qui décroissent algébriquement, $O(t^{-\alpha})$) ne correspondent pas à des valeurs propres à partie réelle négative dans l'équation caractéristique dérivée de l'ansatz exponentiel.
  • Cette limitation est attribuée au fait que les solutions des FODE de type Liouville-Weyl avec des paramètres stables présentent généralement une décroissance algébrique plutôt qu'exponentielle. L'ansatz exponentiel est invalide pour ces modes décroissants car ils n'appartiennent pas à l'espace fonctionnel requis pour que la dérivée de Liouville-Weyl soit définie de la manière standard.
• Validation numérique : Les auteurs démontrent la méthode numériquement sur des systèmes scalaires et 2D (incluant un système de type Mathieu fractionnaire). Ils montrent que la matrice de Hill fractionnaire tronquée identifie avec succès les valeurs propres à partie réelle positive correspondant aux trajectoires instables, confirmant l'utilité de la méthode pour détecter l'instabilité.
• Propriétés de la matrice de Hill fractionnaire : L'article prouve que les valeurs propres de la matrice de Hill fractionnaire tronquée apparaissent par paires conjuguées complexes et sont situées dans des régions spécifiques de type Gershgorin dans le plan complexe.

Importance et affirmations
L'article affirme que le mérite principal de ce travail réside dans la constatation que les solutions de Floquet à décroissance exponentielle n'existent pas pour les LTP-FODE, tandis que les solutions à croissance exponentielle peuvent être observées et évaluées via une approche de type Hill.

Les auteurs déclarent que leur travail offre une voie nouvelle pour accéder aux solutions d'équations différentielles fractionnaires LTP où la théorie de Floquet classique échoue. Ils soulignent que, bien que la méthode soit limitée à la détection de l'instabilité (solutions divergentes), elle offre une alternative rigoureuse à l'analyse des courbes de transition, qui n'identifie que les frontières de stabilité sans caractériser le type de bifurcation ni la nature des modes instables.

L'article conclut avec modestie, notant que l'absence de valeurs propres à partie réelle non négative dans le problème de Hill fractionnaire est un indicateur fort de stabilité, mais nécessite des investigations supplémentaires pour être considéré comme concluant. De plus, les auteurs identifient la construction d'un opérateur de monodromie pour les LTP-FODE comme un défi ouvert, suggérant que le spectre d'un tel opérateur serait nécessaire pour caractériser pleinement la stabilité, y compris les modes à décroissance algébrique que la méthode de Hill actuelle ne capture pas.