[论文解读] Approximating linear response by nonintrusive shadowing algorithms

Angxiu Ni|arXiv (Cornell University)|Mar 22, 2020

Theoretical and Computational Physics参考文献 40被引用 4

一句话总结

该论文证明，非侵入式影子法仅计算混沌系统线性响应中的“影子贡献”，而非完整导数。它证明了剩余的“不稳定贡献”是影子法的系统性误差，并表明当不稳定维数较低时，该误差较小，从而解释了影子法在流体动力学等高维系统中成功的经验现象。

ABSTRACT

Nonintrusive shadowing algorithms efficiently compute $v$, the difference between shadowing trajectories, then use $v$ to compute derivatives of averaged objectives of chaos with respect to parameters of the dynamical system. However, previous proofs of shadowing methods wrongly assume that shadowing trajectories are representative. In contrast, the linear response formula is proved rigorously, but is more difficult to compute. We prove that $v$ gives only a part, called the shadowing contribution, of the linear response; hence, the other part, the unstable contribution, is the systematic error of shadowing methods. For systems with a small ratio of unstable dimensions, with some further statistical assumptions, we show that the unstable contribution is small. We also briefly describe an algorithm for the unstable contribution, which is simpler to derive but less efficient than the fast linear response algorithm. Moreover, we prove the convergence of the nonintrusive shadowing algorithm, the fastest shadowing algorithm, to $v$ and to the shadowing contribution.

研究动机与目标

• 严格分析非侵入式影子算法在近似混沌系统线性响应时的误差。
• 将影子法的系统性误差识别为“不稳定贡献”——即完整线性响应中缺失的部分。
• 证明非侵入式影子算法收敛于影子贡献，而非完整线性响应。
• 提供不稳定贡献较小的条件，解释影子法在高维系统中的经验成功。
• 通过将响应分解为影子贡献与不稳定贡献，为快速线性响应算法奠定理论基础。

提出的方法

• 将线性响应分解为两部分：影子贡献（由影子法计算）与不稳定贡献（即误差）。
• 应用非侵入式影子算法，该算法仅使用 $ m $ 个切线方程的解来计算影子方向 $ v $，其中 $ m $ 为不稳定维数。
• 采用基于采样的误差分析，通过相关性衰减假设来界定采样误差。
• 引入假设3，即 $ \Phi_u v $ 的相关性呈指数衰减，从而对采样误差 $ \|\tilde{e}^S\| $ 给出定量界。
• 利用霍尔德连续性与一致双曲系统中的相关性衰减，推导非侵入式影子算法的收敛界。
• 在后续工作中提出一种快速线性响应算法，利用二阶切线方程高效计算不稳定贡献。

实验结果

研究问题

• RQ1 非侵入式影子算法实际计算了线性响应的哪一部分？
• RQ2 尽管存在已知的理论局限性，为何影子法在流体动力学等高维混沌系统中仍能成功？
• RQ3 影子法的系统性误差是什么？在何种条件下该误差较小？
• RQ4 能否严格证明非侵入式影子算法收敛于影子贡献？
• RQ5 如何通过结合影子法与不稳定贡献的高效计算，恢复完整的线性响应？

主要发现

• 非侵入式影子算法仅计算线性响应中的“影子贡献”，而非完整导数。
• 影子法的系统性误差即为“不稳定贡献”，即线性响应中缺失的部分。
• 在不稳定维数较低且相关性呈指数衰减的系统中，不稳定贡献较小，解释了影子法的经验成功。
• 影子贡献的采样误差满足 $ \|\tilde{e}^S\| \leq \sqrt{\frac{2C_3}{K(1 - \kappa_3)}} \|v^A\| $，在假设3下随 $ K^{-0.5} $ 衰减。
• 非侵入式影子算法收敛于影子贡献的速率与以往影子方法相同，证实其对 $ v $ 的准确性。
• 近似线性响应的总误差由不稳定贡献、采样误差与数值误差之和界定，后两者在 $ K \to \infty $ 时趋于零。