Chaos is not an invariant of fractional-order systems: the clock shapes orbits, attractors and bifurcations

技术摘要：“混沌不是分数阶系统的不变量：时钟塑造轨道、吸引子与分岔”

问题陈述
在常微分方程（整数阶系统）的定性理论中，自庞加莱建立以来，一个基本原则是：通过一个正的、与状态相关的函数（即“时钟”）对向量场进行乘法运算，只会改变轨道经过的速度，而不会改变轨道本身的几何形状。这种重参数化不变性意味着相图、极限环、吸引子和分岔结构是向量场本身的属性，而与时间律无关。

然而，分数阶系统（特别是通过 Caputo 导数定义的系统）具有一个权重分配给过去历史状态的非局部记忆核。本文研究了经典的重参数化不变性是否适用于这些系统。具体而言，研究的问题是：如果一个分数阶系统 $D^\alpha x = f(x)$ 被修改为 $D^\alpha x = g(x)f(x)$（其中 $g(x)$ 是一个正的、与状态相关的函数），其解是否仍然描述相同的几何轨道，仅仅是经过的速度不同？作者指出，对于分数阶 $\alpha \in (0, 1)$，这种不变性失效了。

研究方法
本研究结合了解析展开与数值验证的方法：

1. 解析推导： 作者通过分析 Caputo 微分方程 $D^\alpha x = F(x)$ 的短时行为，利用等效的 Volterra 积分方程进行研究。他们针对参考系统（$F=f$）和时钟系统（$F=gf$），在$t^\alpha$的幂次上进行类似泰勒展开。通过比较原系统解$x(t)$与经变换后的时间$s(t) = \int_0^t g(x(\tau)) d\tau$ 重参数化后的参考解，他们分离出了差异项。
2. 普适函数： 分析得到了两个关于阶数 $\alpha$ 的普适函数，分别记作 $c(\alpha)$ 和 $d(\alpha)$，用于量化与不变性的偏离程度。
3. 数值验证： 理论预测在四个不同的系统中得到了测试：
   • 一个非线性振子（用于验证系数 $c(\alpha)$ 和 $d(\alpha)$）。
   • 分数阶范德波尔（van der Pol）振子（用于检查极限环的形变）。
   • 分数阶 Rössler 系统（用于研究倍周期级联和分岔阈值）。
   • 分数阶 Chen 系统（用于测试混沌机制和周期窗口的稳定性）。
     使用了两种数值方案：分数阶 Adams-Bashforth-Moulton 预测-校正法以及扩散表示法（指数模叠加），以确保结果的稳健性。

核心贡献与结果

1. 时钟协变性的失效： 本文证明了对于 $0 < \alpha < 1$，分数阶系统的轨道在状态相关时钟下不是不变的。

   • 时间律差异（$t^{2\alpha}$）： 在第一个非零阶处，时钟解与重参数化参考解之间的差异与轨道相切。该差异受普适函数控制：
     $$c(\alpha) = \frac{1}{\Gamma(1+2\alpha)} - \frac{\alpha}{\Gamma(1+\alpha)\Gamma(2+\alpha)}$$
     该函数仅在 $\alpha=1$ 时消失（恢复整数阶定理），并在 $\alpha \to 0$ 时趋于 1。
   • 几何差异（$t^{3\alpha}$）： 在下一阶，会出现一个横截分量，证明轨道是几何上不同的曲线。这由第二个普适函数 $d(\alpha)$ 控制，且对于所有 $\alpha \in (0, 1)$，$d(\alpha)$ 均不为零。
   • 数值确认： 通过数值计算，系数 $c(\alpha)$ 和 $d(\alpha)$ 的精度得到了高度验证（$c$ 的误差 < 0.1%，$d$ 的误差 < 1%），符合理论预测。

2. 局部稳定性的不变性： 作者证明了局部性质保持不变。平衡点、线性稳定性以及 Matignon 临界阶（决定稳定性丧失的阈值 $\alpha^*$，由特征值的辐角决定）不受时钟影响。这是因为 Matignon 判据取决于特征值的辐角，而乘以正标量 $g(x^*)$ 不会改变辐角。

3. 全局动力学的不变性失效： 虽然局部稳定性得以保留，但全局结构发生了形变。

   • 极限环： 在分数阶范德波尔振子中，改变时钟会使极限环的形状和振幅发生形变，在中间阶数时变化幅度高达 16%，这在整数阶中是不可能发生的现象。
   • 分岔： 在分数阶 Rössler 系统中，时钟会移动局部极大值的值，但不会改变倍周期分岔的序列或混沌发生的临界阶。
   • 混沌与周期性： 在分数阶 Chen 系统中，一个状态相关的 Arrhenius 时钟（依赖于类温度变量 $z$）可以在相同的阶数 $\alpha$ 和相同的初始条件下，将周期机制转化为混沌机制（或反之）。具体而言，在 $\alpha=0.94$ 时，使用平凡时钟时系统是周期的，而使用特定的 Arrhenius 时钟时则是混沌的。

4. 通用规则： 本文确立了分岔移动的规则：

   • 局部分岔（平衡点稳定性、Hopf 分岔）不会移动。
   • 全局分岔（倍周期、危机、周期窗口的产生/消亡）随时钟强度的增加而线性偏移。
   • 系统对时钟变化的“脆弱性”取决于其与全局分岔的接近程度。

意义与主张

本文声称，对于分数阶系统，“时钟”不仅仅是一个数学上的便利，它是一个模型的物理组成部分。

• 物理阐释： 在整数阶模型中，通过一个与状态相关的因子（例如随温度变化的速率）对场进行乘法运算，等同于重新定义时间单位。在分数阶模型中，这种选择是一个截然不同的物理假设：

  • $D^\alpha x = g(x)f(x)$ 意味着记忆是在物理时间中累积的。
  • $D^\alpha X = f(X)$ 在内在时间 $s$ 中运行，意味着记忆是在内在/材料时间中累积的（这是粘弹性学和玻璃态动力学中的一个概念）。
    本文断言这是两种不同的模型，具有不同的预测，且其差异由普适函数 $c(\alpha)$ 和 $d(\alpha)$ 进行量化。

• 对混沌理论的启示： 作者得出结论，“混沌不是分数阶系统的不变量”。最小混沌阶、分岔图以及奇异吸引子的形状都是特定的场-时钟对的属性，而非系统本身的属性。因此，基于 Matignon 判据的混沌必要条件是稳健的，但混沌的实际发生取决于动力学因子的具体形式。

• 实际后果： 使用分数阶模型的研究人员必须明确声明记忆是在物理时间中累积，还是在内在的材料时间中累积。分数阶吸引子在状态空间中的绘制，是向量场与运行它的时钟共同作用的结果；不同的时钟会产生不同的吸引子，甚至可能导致性质的变化（如从周期性变为混沌）。