[论文解读] Quantifying the computability of the Lorenz system

Benjamin Kehlet, Anders Logg|arXiv (Cornell University)|Jun 12, 2013

Cellular Automata and Applications参考文献 18被引用 7

一句话总结

本文通过分析舍入误差随时间对数值精度的限制，量化了洛伦兹系统可计算性的程度。利用基于伴随的后验误差分析，推导出一个模型，表明可计算性与机器精度呈线性关系，解的可计算时间 T ≈ 2.5 × n_mach，其中 n_mach 为有效数字位数；例如，420 位精度可使计算时间达到 T ≈ 1050。

ABSTRACT

It is well known that the computation of accurate trajectories of the Lorenz system is a difficult problem. Computed solutions are very sensitive to the discretization error determined by the time step size and polynomial order of the method, as well as round-off errors. In this work, we show how round-off errors limit the computability of the Lorenz system and quantify exactly the length of intervals over which solutions can be computed, expressed in terms of the floating point precision. Using adjoint-based \emph{a posteriori} error analysis techniques, we estimate the stability of computations with respect to initial data, discretization, and round-off errors, respectively. The analysis is verified by computing an accurate solution on the time interval $[0, 1000]$ using a very high order (order 200) finite element method and very high floating point precision 400 digits).

研究动机与目标

• 解决一个基本问题：对于给定的机器精度，洛伦兹系统能被计算多远？
• 量化舍入误差在混沌系统（如洛伦兹吸引子）数值解中的限制作用。
• 证明标准的先验误差估计过于悲观，无法捕捉真实的可计算性极限。
• 通过使用高阶（q=200）有限元方法与超高的 400 位精度算术，验证理论模型。

提出的方法

• 采用基于伴随的后验误差分析方法，估算对初值、离散化和舍入误差的敏感性。
• 推导出综合误差模型：E ≈ [0.5‖U(0)−u(0)‖ + 0.001Δt^{2q} + (0.002 + 0.0005q)Δt^{-1/2}ε_mach] × 10^{0.388T}，涵盖所有主要误差来源。
• 采用高阶连续伽辽金（cG(q)）有限元方法，多项式阶数 q=200，在区间 [0,1000] 上计算参考解。
• 采用 400 位浮点精度以最小化舍入误差，并验证理论误差模型。
• 通过数值实验，从经验上拟合常数 C₂^{[q]} 和 C₃^{[q]}，以校准误差模型。
• 通过平衡离散化误差与舍入误差的贡献，推导出最优时间步长 Δt ≈ ε_mach^{1/(2q+0.5)}。

实验结果

研究问题

• RQ1 对于给定的机器精度，洛伦兹系统能被准确计算的最大时间区间 T 是多少？
• RQ2 与离散化误差相比，舍入误差如何限制混沌系统（如洛伦兹系统）的可计算性？
• RQ3 后验误差分析能否提供比标准先验界更准确、更真实的可计算性估计？
• RQ4 随着机器精度（有效数字位数）的提高，可计算性如何变化？

主要发现

• 洛伦兹系统的可计算性与机器精度呈线性关系：T ≈ 2.5 × n_mach，其中 n_mach 为有效数字位数。
• 在 16 位精度（标准双精度）下，可计算性被限制在 T ≈ 40，远低于通常假设的 T ≈ 1.1（来自先验估计）。
• 在 420 位精度下，该方法可计算出在 T = 1000 时仍保持准确的解，与理论预测的 T ≈ 1050 一致。
• 误差模型 E ≈ 0.002 × ε_mach × 10^{0.4T} 准确预测了误差增长，当 ε_mach = 10^{-16} 时，E ≈ 0.001 在 T ≈ 50 处成立。
• 最优时间步长为 Δt ≈ ε_mach^{1/(2q+0.5)}，当 q=100 且 ε_mach=10^{-420} 时，Δt ≈ 0.008，接近实际使用的 Δt=0.0037。
• 在 T=1000 时，稳定性因子 S_C(T) ≈ 2 × 10^{388}，确认了误差增长率为 10^{0.388T}，验证了误差模型中指数项的正确性。