Stability Analysis of Pantograph Delay Differential Equations

技术摘要： pantograph 时滞微分方程的稳定性分析

问题陈述
本文研究标量 pantograph 时滞微分方程（DDE）的稳定性性质：
$$\dot{x}(t) = ax(t) + bx(qt), \quad 0 < q < 1$$
其中时滞与当前时间成正比，表示为 $x(qt) = x(t - (1-q)t)$。与常时滞 DDE 不同，pantograph 方程无法通过标准指数试探解（$x(t) = e^{\lambda t}$）导出简单的特征拟多项式，致使经典谱稳定性理论不再适用。此外，由于显式依赖于 $t$ 和 $qt$，无法直接约化为特征方程，且处理比例时滞的标准常微分方程数值方法需要进行重大修改。

方法论
本研究采用结合解析技术与数值模拟的混合方法：

1. 级数表示：作者利用该方程已知的幂级数解，推导出基本的不稳定性充分条件。
2. Lyapunov 方法：
   • Lyapunov–Krasovskii 理论：构造泛函以推导渐近稳定性的充分条件。
   • Lyapunov–Chetaev 理论：定义函数以建立不稳定性的充分条件。
3. 数值模拟：实施预测 - 校正方法（详见附录），以验证理论发现并探索解析结果尚不明确的参数平面区域（$a, b, q$）。当延迟项 $x(qt)$ 的自变量落在离散时间步长之间时，该方法利用线性插值进行处理。
4. 向非线性系统的扩展：将稳定性框架应用于 Mackey–Glass 方程的 $q$-模拟，以研究分岔和混沌行为。

主要贡献与结果

• 解析稳定性判据：

  • 不稳定性：若 $a > 0$ 且 $b > 0$，或若 $a > -b > 0$，则系统被证明是不稳定的。此外，若 $a + b > 0$，则对于所有 $q \in (0, 1)$，系统均不稳定。
  • 渐近稳定性：推导出了渐近稳定性的严格充分条件：$a < 0$，$a < b < -a$，且 $b^2/a^2 < q < 1$。作者指出，该条件是充分的但非必要的。
  • 参数平面划分：研究将参数平面映射为三个截然不同的区域：与时滞无关的不稳定区、与时滞无关的稳定区以及时滞依赖的稳定区。

• 数值观察与改进：

  • 精细化的稳定区域：数值计算表明，由 Lyapunov–Krasovskii 泛函导出的充分条件较为保守。在由 $a < 0$ 和 $a < b < -a$ 定义的更宽区域内，系统似乎对于 所有 $q \in (0, 1)$ 均呈现渐近稳定性。
  • 阈值曲线（$b^*$）：在 $a < 0$ 且 $b < a$ 的区域中，模拟表明存在一个阈值 $b^*(a)$，它将导致对所有 $q$ 均不稳定的参数与表现出时滞依赖稳定性的参数区分开来。该阈值通过二次拟合经验近似为：
    $$b^*(a) \approx 0.014455a^2 + 1.391065a - 0.03994898$$
  • 动力学行为：本文在边界情况下识别出“超慢”收敛或发散现象，这通常伴随着振荡，且由于时滞 $\tau(t) = (1-q)t$ 的无界增长，振荡周期随时间增加。

• Mackey–Glass 模拟：

  • 作者构建了 Mackey–Glass 方程的比例时滞版本：$\dot{x}(t) = -\beta x(t) + \frac{\alpha x(qt)}{1 + x(qt)^c}$。
  • 通过在非零平衡点附近线性化，该问题简化为具有特定系数的 pantograph 方程。
  • 确定了参数 $\alpha$ 的两个分岔阈值：
    $$\alpha^*_1 = \frac{\beta c}{c - 2}, \quad \alpha^*_2 = \frac{\beta^2 c}{b^* + \beta(c - 1)}$$
  • 这些阈值划分了与时滞无关的稳定区、时滞依赖的稳定区以及不稳定区。当 $\alpha > \alpha^*_2$ 时，系统表现出混沌振荡和共存混沌吸引子。

意义与局限性
本文的主要意义在于解决了比例时滞带来的解析挑战，这些挑战排除了标准特征方程方法的应用。通过结合级数解、Lyapunov 理论和数值验证，该工作为线性 pantograph 方程提供了全面的稳定性图谱，并将这些见解扩展至非线性混沌系统。

作者明确承认其关于稳定性边界的声明具有适度性。他们指出，关于整个区域 $a < 0, a < b < -a$ 对所有 $q$ 均渐近稳定的猜想需要严格的证明。同样，阈值曲线 $b^*(a)$ 目前是基于经验的；从解析上推导它或证明精确界限仍是一个未决问题。这项工作为 pantograph 型时滞系统更完整的理论奠定了基础，突出了比例时滞所固有的独特动力学特征（如超慢发散和随时间变化的振荡周期）。