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

技术摘要：分数阶微分方程周期解的希尔型稳定性分析

问题陈述
本文探讨了非线性分数阶微分方程（FODEs）中周期解的稳定性分析。经典 Caputo 型分数阶微分方程的一个根本局限性在于它们不存在严格的周期解；在 Caputo 导数意义下，周期函数仅是渐近周期的。因此，标准的稳定性分析工具（如 Floquet 理论和希尔矩阵方法）无法直接应用于 Caputo 系统中周期解的线性化扰动动力学。虽然 Liouville-Weyl 型分数阶微分方程（定义为具有无限下限 $-\infty$）确实存在严格的周期解，但针对这些系统的严格 Floquet 理论尚未建立。此前将 Floquet 理论推广至分数阶系统的尝试，要么依赖于近似（例如，假设 Mittag-Leffler 函数具有乘法恒等式），要么未能全面刻画稳定性边界和分岔类型。

方法论
作者提出了一种基于 Liouville-Weyl 型分数阶微分方程的框架来研究周期解。核心方法论包括：

1. 系统构建：定义一类具有无限下限（$-\infty$）的非线性 Liouville-Weyl 微分方程。初值问题通过定义在 $(-\infty, t_0]$ 上的初始化函数 $x_0(t)$ 来构建，该函数属于特定的函数空间 $X(-\infty, t_0]$，以确保由系统历史推导出的强迫项存在。
2. 线性化：通过围绕周期解 $x_p(t)$ 对系统进行线性化来分析其稳定性，从而得到一个具有周期雅可比矩阵 $J(t)$ 的线性时变周期分数阶微分方程（LTP-FODE）。
3. Liouville-Weyl 导数性质：利用指数函数 $e^{st}$ 的 Liouville-Weyl 导数为 $s^\alpha e^{st}$（对于 $\text{Re}(s) \ge 0$）这一性质，作者尝试将 Floquet 假设 $y(t) = e^{\lambda t} p(t)$ 推广至分数阶系统。
4. 分数阶希尔矩阵：通过将 Floquet 假设代入 LTP-FODE，并利用周期系数的傅里叶级数展开，作者推导出了一个无限维非线性特征值问题。这导致了“分数阶希尔矩阵”$H^{(\alpha)}_\infty(\lambda)$ 的定义，其中对角块包含形如 $(\lambda + ik\omega)^\alpha I - J_k$ 的项。
5. 线性时不变（LTI）子系统分析：在解决周期情况之前，作者分析了具有无限记忆的 LTI-FODEs。他们比较了 Liouville-Weyl 系统与标准 Caputo 系统（零下限）的解，以理解衰减解的行为。

主要贡献与结果

• 周期解框架：本文建立了 Liouville-Weyl 分数阶微分方程的严格初值问题表述，该表述承认严格周期解的存在，从而克服了 Caputo 导数的局限性。
• Floquet 理论的推广：作者推导出了针对 LTP-FODEs 的分数阶希尔矩阵方法。他们证明，如果存在形式为 $y(t) = e^{\lambda t} p(t)$ 且 $\text{Re}(\lambda) \ge 0$ 的解，则该解必须满足非线性特征值问题 $H^{(\alpha)}_\infty(\lambda)p_\infty = 0$。
• 对衰减解的局限性：一个核心发现是，所提出的分数阶希尔方法（以及一般的指数假设）仅能评估指数增长的解（不稳定性）。它无法捕捉代数衰减的解。
  • 在 LTI 情况下，作者表明，虽然发散解对应于具有正实部的特征值，但收敛解（表现为代数衰减，$O(t^{-\alpha})$）并不对应于由指数假设导出的特征方程中具有负实部的特征值。
  • 这一局限性归因于这样一个事实：具有稳定参数的 Liouville-Weyl 分数阶微分方程的解通常表现出代数衰减而非指数衰减。指数假设对于这些衰减模态是无效的，因为它们不属于 Liouville-Weyl 导数以标准方式定义所需的函数空间。
• 数值验证：作者在标量和二维系统（包括分数阶 Mathieu 型系统）上数值演示了该方法。他们表明，截断的分数阶希尔矩阵成功识别了对应于不稳定轨迹的正实部特征值，证实了该方法在检测不稳定性方面的实用性。
• 分数阶希尔矩阵的性质：本文证明了截断分数阶希尔矩阵的特征值以复共轭对的形式出现，并位于复平面中特定的 Gershgorin 型区域内。

意义与主张
本文声称，这项工作的主要价值在于洞察到LTP-FODEs 中不存在指数衰减的 Floquet 解，而指数增长的解可以通过希尔型方法进行观察和评估。

作者指出，他们的工作提供了一条通往 LTP 分数阶微分方程解的新途径，在这些方程中经典的 Floquet 理论失效。他们强调，虽然该方法仅限于检测不稳定性（发散解），但它为过渡曲线分析提供了一种严格的替代方案，后者仅能识别稳定性边界，而无法刻画分岔类型或不稳定模态的性质。

本文结论较为谦逊，指出分数阶希尔问题中不存在非负实部特征值是稳定性的强烈指标，但需要进一步研究才能被视为确凿结论。此外，作者指出为 LTP-FODEs 构建单值算子（monodromy operator）是一个未解决的挑战，表明此类算子的谱对于全面刻画稳定性（包括当前希尔方法所遗漏的代数衰减模态）是必要的。