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

Resumo Técnico: Análise de Estabilidade do Tipo Hill de Soluções Periódicas de Equações Diferenciais de Ordem Fracionária

Enunciado do Problema
O artigo aborda a análise de estabilidade de soluções periódicas em equações diferenciais de ordem fracionária (FODEs) não lineares. Uma limitação fundamental das FODEs clássicas do tipo Caputo é que elas não admitem soluções estritamente periódicas; funções periódicas sob derivadas de Caputo são apenas assintoticamente periódicas. Consequentemente, ferramentas padrão de análise de estabilidade, como a teoria de Floquet e o método da matriz de Hill, não podem ser aplicadas diretamente à dinâmica de perturbação linearizada de soluções periódicas em sistemas de Caputo. Embora as FODEs do tipo Liouville-Weyl (definidas com um limite inferior infinito, $-\infty$) admitam soluções estritamente periódicas, uma teoria de Floquet rigorosa para esses sistemas ainda não foi estabelecida. Tentativas anteriores de estender a teoria de Floquet para sistemas de ordem fracionária basearam-se em aproximações (por exemplo, assumindo identidades multiplicativas para funções de Mittag-Leffler) ou falharam em caracterizar de forma abrangente os limites de estabilidade e os tipos de bifurcação.

Metodologia
Os autores propõem um framework baseado em FODEs do tipo Liouville-Weyl para estudar soluções periódicas. A metodologia central envolve:

1. Formulação do Sistema: Definição de uma classe de equações diferenciais não lineares de Liouville-Weyl com um limite inferior infinito ($-\infty$). O problema de valor inicial é formulado com uma função de inicialização $x_0(t)$ definida em $(-\infty, t_0]$, pertencendo a um espaço funcional específico $X(-\infty, t_0]$ que garante a existência do termo de forçamento derivado da história do sistema.
2. Linearização: A estabilidade de uma solução periódica $x_p(t)$ é analisada linearizando o sistema em torno dela, resultando em uma equação diferencial de ordem fracionária linear de período temporal (LTP-FODE) com uma matriz Jacobiana periódica $J(t)$.
3. Propriedades da Derivada de Liouville-Weyl: Utilizando a propriedade de que a derivada de Liouville-Weyl de uma função exponencial $e^{st}$ é $s^\alpha e^{st}$ (para $\text{Re}(s) \ge 0$), os autores tentam estender o ansatz de Floquet $y(t) = e^{\lambda t} p(t)$ para sistemas de ordem fracionária.
4. Matriz de Hill Fracionária: Ao substituir o ansatz de Floquet na LTP-FODE e utilizar expansões em série de Fourier para os coeficientes periódicos, os autores derivam um problema de autovalor não linear de dimensão infinita. Isso leva à definição da "matriz de Hill fracionária", $H^{(\alpha)}_\infty(\lambda)$, onde os blocos diagonais envolvem termos da forma $(\lambda + ik\omega)^\alpha I - J_k$.
5. Análise de Subsistemas Lineares Invariantes no Tempo (LTI): Antes de abordar o caso periódico, os autores analisam LTP-FODEs com memória infinita. Eles comparam as soluções do sistema de Liouville-Weyl com as do sistema de Caputo padrão (limite inferior zero) para compreender o comportamento das soluções decrescentes.

Principais Contribuições e Resultados

• Framework para Soluções Periódicas: O artigo estabelece uma formulação rigorosa de problema de valor inicial para FODEs de Liouville-Weyl que admite soluções estritamente periódicas, superando a limitação das derivadas de Caputo.
• Extensão da Teoria de Floquet: Os autores derivam uma abordagem de matriz de Hill fracionária para LTP-FODEs. Eles provam que, se uma solução existir na forma $y(t) = e^{\lambda t} p(t)$ com $\text{Re}(\lambda) \ge 0$, ela deve satisfazer o problema de autovalor não linear $H^{(\alpha)}_\infty(\lambda)p_\infty = 0$.
• Limitação em Soluções Decrescentes: Uma descoberta central é que o método da matriz de Hill fracionária proposto (e o ansatz exponencial em geral) pode apenas avaliar soluções de crescimento exponencial (instabilidade). Ele falha em capturar soluções de decaimento algébrico.
  • No caso LTI, os autores mostram que, enquanto soluções divergentes correspondem a autovalores com partes reais positivas, soluções convergentes (que decaem algebricamente, $O(t^{-\alpha})$) não correspondem a autovalores com partes reais negativas na equação característica derivada do ansatz exponencial.
  • Essa limitação é rastreada ao fato de que as soluções de FODEs de Liouville-Weyl com parâmetros estáveis tipicamente exibem decaimento algébrico em vez de decaimento exponencial. O ansatz exponencial é inválido para esses modos de decaimento porque eles não pertencem ao espaço funcional necessário para que a derivada de Liouville-Weyl seja definida da maneira padrão.
• Validação Numérica: Os autores demonstram o método numericamente em sistemas escalares e bidimensionais (incluindo um sistema do tipo Mathieu de ordem fracionária). Eles mostram que a matriz de Hill fracionária truncada identifica com sucesso autovalores com partes reais positivas correspondentes a trajetórias instáveis, confirmando a utilidade do método para detectar instabilidade.
• Propriedades da Matriz de Hill Fracionária: O artigo prova que os autovalores da matriz de Hill fracionária truncada vêm em pares de conjugados complexos e estão localizados dentro de regiões específicas semelhantes a Gershgorin no plano complexo.

Significância e Afirmações
O artigo afirma que o mérito principal deste trabalho é a percepção de que soluções de Floquet de decaimento exponencial não existem para LTP-FODEs, enquanto soluções de crescimento exponencial podem ser observadas e avaliadas por meio de uma abordagem do tipo Hill.

Os autores afirmam que seu trabalho fornece um caminho novel para acessar soluções de equações diferenciais fracionárias LTP onde a teoria clássica de Floquet falha. Eles enfatizam que, embora o método seja limitado à detecção de instabilidade (soluções divergentes), oferece uma alternativa rigorosa à análise de curvas de transição, que apenas identifica limites de estabilidade sem caracterizar o tipo de bifurcação ou a natureza dos modos instáveis.

O artigo conclui com modéstia, observando que a ausência de autovalores com partes reais não negativas no problema da matriz de Hill fracionária é um forte indicador de estabilidade, mas requer investigação adicional para ser considerado conclusivo. Além disso, os autores identificam a construção de um operador de monodromia para LTP-FODEs como um desafio em aberto, sugerindo que o espectro de tal operador seria necessário para caracterizar completamente a estabilidade, incluindo os modos de decaimento algébrico que o método de Hill atual perde.