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

Resumen Técnico: Análisis de Estabilidad de Tipo Hill de Soluciones Periódicas de Ecuaciones Diferenciales de Orden Fraccionario

Planteamiento del Problema
El artículo aborda el análisis de estabilidad de soluciones periódicas en ecuaciones diferenciales de orden fraccionario (FODEs) no lineales. Una limitación fundamental de las FODEs clásicas de tipo Caputo es que no admiten soluciones exactamente periódicas; las funciones periódicas bajo derivadas de Caputo son únicamente asintóticamente periódicas. En consecuencia, las herramientas estándar de análisis de estabilidad, como la teoría de Floquet y el método de la matriz de Hill, no pueden aplicarse directamente a la dinámica de perturbación linealizada de soluciones periódicas en sistemas de Caputo. Aunque las FODEs de tipo Liouville-Weyl (definidas con un límite inferior infinito, $-\infty$) sí admiten soluciones exactamente periódicas, una teoría de Floquet rigurosa para estos sistemas no ha sido establecida. Los intentos previos de extender la teoría de Floquet a sistemas de orden fraccionario han dependido de aproximaciones (por ejemplo, asumiendo identidades multiplicativas para funciones de Mittag-Leffler) o han fallado en caracterizar exhaustivamente los límites de estabilidad y los tipos de bifurcación.

Metodología
Los autores proponen un marco basado en FODEs de tipo Liouville-Weyl para estudiar soluciones periódicas. La metodología central incluye:

1. Formulación del Sistema: Definición de una clase de ecuaciones diferenciales no lineales de tipo Liouville-Weyl con un límite inferior infinito ($-\infty$). El problema de valor inicial se formula con una función de inicialización $x_0(t)$ definida en $(-\infty, t_0]$, perteneciente a un espacio funcional específico $X(-\infty, t_0]$ que garantiza la existencia del término forzante derivado de la historia del sistema.
2. Linealización: La estabilidad de una solución periódica $x_p(t)$ se analiza linealizando el sistema alrededor de ella, resultando en una ecuación diferencial de orden fraccionario lineal con coeficientes periódicos en el tiempo (LTP-FODE) con una matriz jacobiana periódica $J(t)$.
3. Propiedades de la Derivada de Liouville-Weyl: Utilizando la propiedad de que la derivada de Liouville-Weyl de una función exponencial $e^{st}$ es $s^\alpha e^{st}$ (para $\text{Re}(s) \ge 0$), los autores intentan extender el ansatz de Floquet $y(t) = e^{\lambda t} p(t)$ a sistemas de orden fraccionario.
4. Matriz de Hill Fraccionaria: Al sustituir el ansatz de Floquet en la LTP-FODE y utilizar expansiones en series de Fourier para los coeficientes periódicos, los autores derivan un problema de autovalores no lineal de dimensión infinita. Esto conduce a la definición de la "matriz de Hill fraccionaria", $H^{(\alpha)}_\infty(\lambda)$, donde los bloques diagonales involucran términos de la forma $(\lambda + ik\omega)^\alpha I - J_k$.
5. Análisis de Subsistemas Lineales Invariantes en el Tiempo (LTI): Antes de abordar el caso periódico, los autores analizan LTP-FODEs con memoria infinita. Comparan las soluciones del sistema de Liouville-Weyl con las del sistema de Caputo estándar (límite inferior cero) para comprender el comportamiento de las soluciones decrecientes.

Contribuciones y Resultados Clave

• Marco para Soluciones Periódicas: El artículo establece una formulación rigurosa del problema de valor inicial para FODEs de tipo Liouville-Weyl que admite soluciones exactamente periódicas, superando la limitación de las derivadas de Caputo.
• Extensión de la Teoría de Floquet: Los autores derivan un enfoque de matriz de Hill fraccionaria para LTP-FODEs. Demuestran que si existe una solución en la forma $y(t) = e^{\lambda t} p(t)$ con $\text{Re}(\lambda) \ge 0$, esta debe satisfacer el problema de autovalores no lineal $H^{(\alpha)}_\infty(\lambda)p_\infty = 0$.
• Limitación en Soluciones Decrecientes: Un hallazgo central es que el método de Hill fraccionario propuesto (y el ansatz exponencial en general) solo puede evaluar soluciones exponencialmente crecientes (inestabilidad). No logra capturar soluciones algebraicamente decrecientes.
  • En el caso LTI, los autores muestran que, mientras que las soluciones divergentes corresponden a autovalores con partes reales positivas, las soluciones convergentes (que decaen algebraicamente, $O(t^{-\alpha})$) no corresponden a autovalores con partes reales negativas en la ecuación característica derivada del ansatz exponencial.
  • Esta limitación se atribuye al hecho de que las soluciones a FODEs de tipo Liouville-Weyl con parámetros estables típicamente exhiben decaimiento algebraico en lugar de decaimiento exponencial. El ansatz exponencial es inválido para estos modos decrecientes porque no pertenecen al espacio funcional requerido para que la derivada de Liouville-Weyl esté definida de la manera estándar.
• Validación Numérica: Los autores demuestran el método numéricamente en sistemas escalares y bidimensionales (incluyendo un sistema de tipo Mathieu fraccionario). Muestran que la matriz de Hill fraccionaria truncada identifica con éxito los autovalores con partes reales positivas correspondientes a trayectorias inestables, confirmando la utilidad del método para detectar inestabilidad.
• Propiedades de la Matriz de Hill Fraccionaria: El artículo demuestra que los autovalores de la matriz de Hill fraccionaria truncada aparecen en pares conjugados complejos y se ubican dentro de regiones específicas similares a las de Gershgorin en el plano complejo.

Significado y Afirmaciones
El artículo afirma que el mérito principal de este trabajo es la comprensión de que las soluciones de Floquet que decaen exponencialmente no existen para LTP-FODEs, mientras que las soluciones exponencialmente crecientes pueden observarse y evaluarse mediante un enfoque de tipo Hill.

Los autores declaran que su trabajo proporciona una vía novedosa para acceder a soluciones de ecuaciones diferenciales fraccionarias LTP donde la teoría de Floquet clásica falla. Subrayan que, aunque el método se limita a detectar inestabilidad (soluciones divergentes), ofrece una alternativa rigurosa al análisis de curvas de transición, el cual solo identifica límites de estabilidad sin caracterizar el tipo de bifurcación ni la naturaleza de los modos inestables.

El artículo concluye con modestia, señalando que la ausencia de autovalores con partes reales no negativas en el problema de Hill fraccionario es un fuerte indicador de estabilidad, pero requiere mayor investigación para considerarse concluyente. Además, los autores identifican la construcción de un operador de monodromía para LTP-FODEs como un desafío abierto, sugiriendo que el espectro de tal operador sería necesario para caracterizar completamente la estabilidad, incluyendo los modos algebraicamente decrecientes que el método de Hill actual pasa por alto.