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

Sintesi Tecnica: Analisi di Stabilità di Tipo Hill delle Soluzioni Periodiche di Equazioni Differenziali di Ordine Frazionario

Enunciato del Problema
Il lavoro affronta l'analisi di stabilità delle soluzioni periodiche nelle equazioni differenziali non lineari di ordine frazionario (FODE). Un limite fondamentale delle FODE di tipo Caputo classico è che non ammettono soluzioni esattamente periodiche; le funzioni periodiche sotto derivate di tipo Caputo sono solo asintoticamente periodiche. Di conseguenza, gli strumenti standard di analisi di stabilità, come la teoria di Floquet e il metodo della matrice di Hill, non possono essere applicati direttamente alla dinamica delle perturbazioni linearizzate delle soluzioni periodiche nei sistemi di tipo Caputo. Sebbene le FODE di tipo Liouville-Weyl (definite con un limite inferiore infinito, $-\infty$) ammettano soluzioni esattamente periodiche, una teoria di Floquet rigorosa per tali sistemi non è stata ancora stabilita. I precedenti tentativi di estendere la teoria di Floquet ai sistemi di ordine frazionario si sono basati su approssimazioni (ad esempio, assumendo identità moltiplicative per le funzioni di Mittag-Leffler) o non sono riusciti a caratterizzare in modo completo i confini di stabilità e i tipi di biforcazione.

Metodologia
Gli autori propongono un quadro basato su FODE di tipo Liouville-Weyl per studiare le soluzioni periodiche. La metodologia centrale comprende:

1. Formulazione del Sistema: Definizione di una classe di equazioni differenziali non lineari di tipo Liouville-Weyl con un limite inferiore infinito ($-\infty$). Il problema ai valori iniziali è formulato con una funzione di inizializzazione $x_0(t)$ definita su $(-\infty, t_0]$, appartenente a uno specifico spazio funzionale $X(-\infty, t_0]$ che garantisce l'esistenza del termine forzante derivato dalla storia del sistema.
2. Linearizzazione: La stabilità di una soluzione periodica $x_p(t)$ è analizzata linearizzando il sistema attorno ad essa, risultando in un'equazione differenziale di ordine frazionario lineare tempo-periodica (LTP-FODE) con una matrice Jacobiana periodica $J(t)$.
3. Proprietà della Derivata di Liouville-Weyl: Sfruttando la proprietà per cui la derivata di Liouville-Weyl di una funzione esponenziale $e^{st}$ è $s^\alpha e^{st}$ (per $\text{Re}(s) \ge 0$), gli autori tentano di estendere l'ansatz di Floquet $y(t) = e^{\lambda t} p(t)$ ai sistemi di ordine frazionario.
4. Matrice di Hill Frazionaria: Sostituendo l'ansatz di Floquet nell'LTP-FODE e utilizzando sviluppi in serie di Fourier per i coefficienti periodici, gli autori derivano un problema agli autovalori non lineare infinito-dimensionale. Ciò porta alla definizione della "matrice di Hill frazionaria", $H^{(\alpha)}_\infty(\lambda)$, dove i blocchi diagonali coinvolgono termini della forma $(\lambda + ik\omega)^\alpha I - J_k$.
5. Analisi dei Sottosistemi Lineari Tempo-Invarianti (LTI): Prima di affrontare il caso periodico, gli autori analizzano le LTP-FODE con memoria infinita. Confrontano le soluzioni del sistema di Liouville-Weyl con quelle del sistema Caputo standard (limite inferiore zero) per comprendere il comportamento delle soluzioni decrescenti.

Contributi e Risultati Chiave

• Quadro per le Soluzioni Periodiche: Il lavoro stabilisce una formulazione rigorosa del problema ai valori iniziali per le FODE di tipo Liouville-Weyl che ammette soluzioni esattamente periodiche, superando il limite delle derivate di tipo Caputo.
• Estensione della Teoria di Floquet: Gli autori derivano un approccio basato sulla matrice di Hill frazionaria per le LTP-FODE. Dimostrano che se esiste una soluzione della forma $y(t) = e^{\lambda t} p(t)$ con $\text{Re}(\lambda) \ge 0$, essa deve soddisfare il problema agli autovalori non lineare $H^{(\alpha)}_\infty(\lambda)p_\infty = 0$.
• Limitazione sulle Soluzioni Decrescenti: Una scoperta centrale è che il metodo della matrice di Hill frazionaria proposto (e l'ansatz esponenziale in generale) può valutare solo le soluzioni esponenzialmente crescenti (instabilità). Non riesce a catturare le soluzioni algebricamente decrescenti.
  • Nel caso LTI, gli autori mostrano che, mentre le soluzioni divergenti corrispondono ad autovalori con parte reale positiva, le soluzioni convergenti (che decadono algebricamente, $O(t^{-\alpha})$) non corrispondono ad autovalori con parte reale negativa nell'equazione caratteristica derivata dall'ansatz esponenziale.
  • Questa limitazione è riconducibile al fatto che le soluzioni alle FODE di tipo Liouville-Weyl con parametri stabili tipicamente esibiscono un decadimento algebrico piuttosto che esponenziale. L'ansatz esponenziale è invalido per questi modi decrescenti poiché non appartengono allo spazio funzionale richiesto affinché la derivata di Liouville-Weyl sia definita nel modo standard.
• Validazione Numerica: Gli autori dimostrano il metodo numericamente su sistemi scalari e 2D (incluso un sistema di tipo Mathieu frazionario). Mostrano che la matrice di Hill frazionaria troncata identifica con successo gli autovalori con parte reale positiva corrispondenti a traiettorie instabili, confermando l'utilità del metodo per rilevare l'instabilità.
• Proprietà della Matrice di Hill Frazionaria: Il lavoro dimostra che gli autovalori della matrice di Hill frazionaria troncata compaiono in coppie coniugate complesse e sono localizzati all'interno di specifiche regioni di tipo Gershgorin nel piano complesso.

Significato e Affermazioni
Il lavoro afferma che il merito principale di questo studio è la consapevolezza che le soluzioni di Floquet esponenzialmente decrescenti non esistono per le LTP-FODE, mentre le soluzioni esponenzialmente crescenti possono essere osservate e valutate tramite un approccio di tipo Hill.

Gli autori dichiarano che il loro lavoro fornisce una via nuova per accedere alle soluzioni di equazioni differenziali frazionarie LTP dove la teoria di Floquet classica fallisce. Sottolineano che, sebbene il metodo sia limitato al rilevamento dell'instabilità (soluzioni divergenti), offre un'alternativa rigorosa all'analisi delle curve di transizione, che identifica solo i confini di stabilità senza caratterizzare il tipo di biforcazione o la natura dei modi instabili.

Il lavoro conclude con modestia, notando che l'assenza di autovalori con parte reale non negativa nel problema della matrice di Hill frazionaria è un forte indicatore di stabilità, ma richiede ulteriori indagini per essere considerato conclusivo. Inoltre, gli autori identificano la costruzione di un operatore di monodromia per le LTP-FODE come una sfida aperta, suggerendo che lo spettro di tale operatore sarebbe necessario per caratterizzare completamente la stabilità, inclusi i modi algebricamente decrescenti che il metodo di Hill attuale non coglie.