[论文解读] Eigenvalues for Two-Lag Linear Delay Differential Equations

David M. Bortz|arXiv (Cornell University)|Jun 27, 2012

Numerical methods for differential equations参考文献 33被引用 6

一句话总结

本文提出了一种级数展开方法，用于计算双时滞线性延迟微分方程（DDEs）的特征值，将此前仅适用于单时滞系统的Lambert W函数方法扩展至双时滞系统，通过推导非线性特征方程根的渐近展开式实现。该方法实现了精确的稳定性分析，通过在表现出Hopf分岔的生态学果蝇模型中的应用得到验证。

ABSTRACT

We consider the class of two-lag linear delay differential equations and develop a series expansion to solve for the roots of the nonlinear characteristic equation. Supporting numerical results are presented along with application of our method to study the stability of a two-lag model from ecology.

研究动机与目标

• 将此前仅适用于单时滞DDEs的Lambert W函数框架扩展至具有多个时滞的双时滞线性DDEs。
• 推导出双时滞DDEs中非线性特征方程根的计算高效渐近级数展开式。
• 对级数展开的精度与收敛性进行数值验证，以支持实际应用。
• 将该方法应用于生态学中的一类双时滞模型，具体为澳大利亚果蝇种群模型，分析其稳定性与Hopf分岔行为。

提出的方法

• 针对形式为 $ \dot{x}(t) = \alpha x(t) + \beta x(t - \tau_1) + \gamma x(t - \tau_2) $ 的双时滞线性DDE，推导出非线性特征方程，进而转化为指数多项式特征值问题。
• 引入一种变换，将特征值表示为广义Lambert W型函数的形式，利用复对数函数及分支相关的展开式。
• 提出一种包含广义Bell多项式与Pochhammer符号的双重求和级数，以渐近方式计算特征值根。
• 通过递归计算Bell多项式 $ \mathrm{BellY} $ 来处理展开项的多变量结构。
• 利用复对数函数与Lambert W函数的主支及高阶分支，捕获复平面上所有无穷多个特征值。
• 在Mathematica中实现该方法，对级数进行符号与数值计算，截断至有限阶次以支持实际计算。

实验结果

研究问题

• RQ1 能否将Lambert W函数框架从单时滞系统扩展至双时滞线性DDEs，以实现特征值的精确计算？
• RQ2 双时滞DDEs的特征方程具有何种结构？其根能否以收敛级数形式系统展开？
• RQ3 与非线性特征值问题的数值根求解方法相比，所提出的级数展开在精度与收敛性方面表现如何？
• RQ4 当延迟参数变化时，双时滞生态学模型（如果蝇种群模型）的稳定性行为如何？
• RQ5 双时滞DDEs中渐近展开所需的系数约束在实际模型中是否通常成立？

主要发现

• 所提出的级数展开能精确计算双时滞线性DDEs的特征值，数值结果表明在中等截断阶次下即实现收敛与高精度。
• 该方法成功识别出主特征值 $ s_0 $ 及其实部 $ \mathrm{Re}(s_0) $，后者决定系统的稳定性，且可作为延迟参数 $ \tau_1 $ 与 $ \tau_2 $ 的函数。
• 在双时滞果蝇模型中，分析揭示在 $ \tau_2 \approx 1.3 $ 处发生Hopf分岔，系统由稳定转为不稳定振荡。
• 该方法正确捕捉了系统由稳定到不稳定动力学的转变，经由初始函数 $ \phi(t) = 1 $（稳定）与 $ \phi(t) = 0.1 $（不稳定）的数值模拟得到验证。
• 推导展开式所需的系数约束在许多实际模型中均成立，表明该方法具有广泛适用性，超越理论框架。
• Mathematica实现支持多分支特征值的高效计算，有效支持稳定性分析与分岔检测。