[论文解读] DDE-BIFTOOL Manual - Bifurcation analysis of delay differential equations

Jan Sieber, Koen Engelborghs+3 more|arXiv (Cornell University)|Jun 27, 2014

Numerical methods for differential equations参考文献 25被引用 72

一句话总结

本手册详述了 DDE-BIFTOOL，这是一个基于 MATLAB 的软件包，用于延迟微分方程（DDEs）的数值分歧分析，支持常数延迟和状态依赖延迟。该工具可使用配点法与延拓技术计算平衡点、周期轨道、稳定性以及一阶和二阶分歧，还扩展支持正规形计算与对称性感知的动力学分析。

ABSTRACT

DDEBIFTOOL is a collection of Matlab routines for numerical bifurcation analysis of systems of delay differential equations with discrete constant and state-dependent delays. The package supports continuation and stability analysis of steady state solutions and periodic solutions. Further one can compute and continue several local and global bifurcations: fold and Hopf bifurcations of steady states; folds, period doublings and torus bifurcations of periodic orbits; and connecting orbits between equilibria. To analyse the stability of steady state solutions, approximations are computed to the rightmost, stability-determining roots of the characteristic equation which can subsequently be used as starting values in a Newton procedure. For periodic solutions, approximations to the Floquet multipliers are computed. The manual describes the structure of the package, its routines, and its data and method parameter structures.

研究动机与目标

• 为延迟微分方程（DDEs）的数值分歧分析提供一个可移植、用户友好的工具，尤其适用于稳态和周期解。
• 支持常数延迟与状态依赖延迟系统，包括稳定性分析与解分支的延拓。
• 通过正规形计算，将数值能力扩展至二阶分歧，如霍普夫分歧、广义霍普夫分歧与鲍廷分歧。
• 通过专用扩展支持对称系统分析，利用相对平衡点与周期轨道计算对称性相关的动力学行为。
• 通过引用核心 DDE-BIFTOOL 发表文献及相关理论扩展，确保结果可复现并正确归因。

提出的方法

• 使用分段多项式逼近与自适应网格选择的配点法，实现 DDE 的高精度数值求解。
• 采用伪弧长延拓法，追踪参数变化下平衡点与周期轨道的解分支。
• 应用数值方法计算最右特征根与 Floquet 乘子，以实现稳定性分析。
• 实现右端项与延迟项（状态依赖情形）的雅可比矩阵，以提升收敛性与计算精度。
• 通过模块化设计与外部延拓框架集成，支持扩展功能，如正规形计算与对称性处理。
• 通过函数句柄与结构化数据类型支持用户自定义系统结构，用于问题定义、求解点与分支的组织。

实验结果

研究问题

• RQ1 如何有效应用数值分歧分析方法于具有常数延迟与状态依赖延迟的 DDE？
• RQ2 哪些数值方法可确保 DDE 中稳定性与周期解计算的鲁棒性与高精度？
• RQ3 如何通过正规形系数检测与分析二阶分歧，如霍普夫与广义霍普夫分歧？
• RQ4 分析 DDE 时面临哪些计算挑战？如何通过自适应网格与雅可比矩阵使用加以缓解？
• RQ5 如何利用 DDE 系统中的对称性与旋转不变性，计算相对平衡点与周期轨道？

主要发现

• DDE-BIFTOOL 成功实现了基于配点法与延拓法的 DDE 数值分歧分析，支持常数延迟与状态依赖延迟。
• 该软件包可通过 ddebiftool_extra_nmfm 扩展，实现对一阶分歧（折叠、倍周期、环面）与二阶分歧（霍普夫、广义霍普夫、鲍廷）的检测与计算。
• 基于 Kuznetsov、Wage 与 Bosschaert 的理论工作，通过扩展实现霍普夫及其相关分歧的正规形系数计算。
• 软件支持强制周期性达 $2\pi$ 的周期解计算，如 phase_oscillator 演示所示。
• 通过自动步长选择与自适应网格技术增强了鲁棒性，但用户应测试离散化效应，因系统具有固有的数值敏感性。
• 该软件包采用 BSD 2-条款许可证发布，可在 SourceForge 获取，版本 3.1 及以上支持正规形计算与改进的稳定性处理。