[论文解读] Discrete averaging for discrete time dynamical systems

Vassili Gelfreich, Arturo Vieiro|arXiv (Cornell University)|Mar 11, 2026

Numerical methods for differential equations被引用 0

一句话总结

该论文为映射开发离散平均化以近似由自治流驱动的动力学，从而实现显式的绝热不变量和统一误差界。它将该方法应用于近单位映射和共振固定点，在数值与解析方面相对于经典平均化具有优势。

ABSTRACT

In this paper we develop the theory of discrete averaging designed to study discrete time dynamical systems defined by iterates of a map. The discrete averaging uses weighted averages over a segment of trajectory to find an autonomous vector field that approximates the original map. The method provides a simple and effective tool for finding adiabatic invariants, both numerically and analytically. It is capable of strengthening various theorems of the classical averaging theory because it eliminates two intermediate steps used in the classical averaging: the suspension procedure that assigns a rapidly oscillating flow to the map and time-dependent coordinate changes that eliminate the dependence on time. We discuss two applications of the discrete averaging - to the dynamics of a near-identity map and to the dynamics of a map in a neighbourhood of a resonant fixed point. We show that the discrete averaging provides explicit uniform bounds for approximation errors. We also show that the discrete averaging can be used to establish domain of validity of adiabatic approximations in numerical experiments.

研究动机与目标

• 引入离散平均化作为工具，以研究由映射迭代定义的离散时间动力学系统。
• 从轨迹片段加权平均中构造内插向量场。
• 获得将映射近似为时间单位流的显式统一误差界。
• 展示在近单位映射和共振邻域的应用。
• 展示离散平均化如何通过避免悬挂和时变坐标变换来加强经典平均化。

提出的方法

• 从 F 的加权轨迹片段平均中定义内插向量场 X_n。
• 将 X_n 表示为 F^k 的线性组合，系数 p_{nk} 由内插多项式导出。
• 证明对于距单位距离为 epsilon 的切向单位映射，F = Phi^1_{X_n} + O(epsilon^{n+1})。
• 建立统一分析界，表明近似误差取决于 analytic diffeomorphisms 的 epsilon/delta 比值。
• 将更高阶的内插量与不需要正则形式转换的绝热不变量联系起来。
• 将该方法应用于共轭保守的 Hénon 映射在共振附近以从原坐标提取一个绝热不变量。

Figure 1: Iterates of the Henon map ( 7 ) (left) and level lines of the Taylor polynomial of the adiabatic invariant $h_{2}$ of order 6 in $x,y$ and one in $\varepsilon$ (right) for $c=1+\varepsilon$ with $\varepsilon=10^{-3}$ .

实验结果

研究问题

• RQ1 离散平均化是否能够在原始坐标下产生一个自治向量场，以紧密近似给定的近单位映射？
• RQ2 将映射近似为内插向量场的时间单位流的显式误差界是什么？
• RQ3 在共振固定点附近离散平均化的表现如何，能否产生可靠的绝热不变量？
• RQ4 离散平均化是否能够在数值实验中确定绝热近似的有效性域？
• RQ5 与经典平均化相比，离散方法在计算效率与解析控制方面的表现如何？

主要发现

• 离散平均化提供显式的内插向量场 X_n，使得 F = Phi^1_{X_n} + O(epsilon^{n+1})。
• 对于解析映射，存在统一的近似误差界，由 epsilon/delta 比值控制，并随 n 的增大而改善。
• 高阶内插量得到的绝热不变量在更高阶误差下得到保持，在 Hénon 映射示例中得到证明。
• 该方法在原坐标中直接构造绝热不变量，无需调用正则形式转换。
• 离散平均化的有效域可以通过在某一域内考察最优插值阶数来数值评估。
• 在某些情形下获得了指数级精度界，与 Neishtadt 定理的 refined 版本相吻合。

Figure 2: Domain of validity of discrete averaging for the Hénon map with $\varepsilon=-10^{-3}$ (left column) and $\varepsilon=10^{-3}$ (right column). In the top row, the colors indicate the optimal $n$ values that minimize $G_{(x,y)}(n)=|X_{2n}(x,y)-X_{2n+2}(x,y)|_{2}$ . The bottom row plots exhi