[论文解读] A Lohner-type algorithm for control systems and ordinary differential inclusions

Tomasz Kapela, Piotr Zgliczyński|ArXiv.org|Dec 6, 2007

Numerical methods for differential equations参考文献 9被引用 4

一句话总结

本文提出了一种Lohner型算法，用于通过将$C^0$-Lohner方法扩展至处理常微分方程的扰动，实现对控制系统和常微分包含中可达集的严格计算。该方法利用分量估计与对数范数，对不确定条件下的解进行界约束，从而实现严格Poincaré映射计算，并证明了在$ uf 0.032, 0.127]$范围内，Kuramoto-Sivashinsky PDE中存在多个周期轨道。该方法已集成至CAPD库中，在有界扰动下表现出优于标准方法的精度。

ABSTRACT

We describe a Lohner-type algorithm for the computation of rigorous upper bounds for reachable set for control systems, solutions of ordinary differential inclusions and perturbations of ODEs.

研究动机与目标

• 开发一种严格数值方法，用于计算具有有界不确定性的控制系统的可达集及微分包含的解。
• 将$C^0$-Lohner算法扩展至处理具有保证误差界之扰动ODE与非自治系统。
• 实现对非自治ODE及小延迟系统在截面处Poincaré映射的精确计算。
• 为耗散型PDE（如Kuramoto-Sivashinsky方程）中周期轨道的存在性提供计算框架。
• 在基准问题上对比并验证两种估计策略——对数范数（LN）与分量估计（CW）方法——的性能。

提出的方法

• 该算法使用区间算术与泰勒模型，对形如$x'(t) = f(x(t), y(t))$的扰动ODE解进行严格包围，其中$|y(t)| < \epsilon$。
• 结合基于单边Lipschitz条件的分量估计与对数范数，推导出在不确定性下解增长的紧致界。
• 关键创新在于提出一个新不等式，将对数范数与单边Lipschitz常数关联，从而提升界估计的紧致性。
• 通过区间牛顿型校正与区间凸包传播，在时间步长上实现严格包围的积分。
• 通过追踪轨迹直至其穿过截面$\alpha(x) = 0$，支持Poincaré映射的计算。
• 实现已集成至CAPD库中，支持$C^0$-Lohner与扩展的分量估计策略。

实验结果

研究问题

• RQ1 能否将Lohner型算法严格扩展至处理具有有界扰动的控制系统与微分包含？
• RQ2 在扰动ODE中，基于分量估计（CW）与对数范数（LN）的估计方法在精度与计算成本方面如何比较？
• RQ3 该算法能否对非自治ODE与小延迟系统实现严格Poincaré映射计算？
• RQ4 该算法在证明Kuramoto-Sivashinsky PDE中存在多个周期轨道方面的性能如何？
• RQ5 当扰动超过向量场分量大小时，该方法是否仍能提供可靠的包围？

主要发现

• 对于$ uf = 0.1$的扰动谐振子，分量估计（CW）方法产生的包围比对数范数（LN）方法更紧致，时间步长从100提升至100,000时，包围直径由1.620降至0.820。
• 该算法成功计算了具有$[-10^{-4}, 10^{-4}]^3$扰动的Rössler系统Poincaré映射，在截面$x=0, x'>0$上得到直径为(0.619, 0.322, 0.003)的严格包围。
• 对于Kuramoto-Sivashinsky PDE，该方法证明了在$ uf 0.032, 0.127)$范围内存在多个周期轨道，包括具有一个不稳定方向的吸引与不稳定轨道。
• 即使扰动大到$ uf = 10$，该算法仍能生成有效包围，其直径随时间成比例增长，但保持可计算性与严格性。
• 在所有测试问题中，CW方法在精度上始终优于LN方法，尤其在较大时间步长与扰动规模下表现更优。
• 对于谐振子在$T = 2\pi$处的解直径估计，该算法在精细时间离散化下收敛至约1.620（LN）与0.820（CW），表现出高精度。