形式陈述
只改变 Euler 更新的先后顺序,为什么会改变长期轨迹?对可分Hamiltonian公理库Hamilton 系统Hamiltonian system · Hamilton 系统由 Hamiltonian 的梯度经标准辛矩阵生成相空间流的常微分方程系统。 ,先更新动量再更新位置的辛 Euler 方法是
两步都显式可算。一般不可分 Hamiltonian 对应的半隐式形式为
第一式通常需要求解。在 足够光滑且关于新动量的块 非奇异时,隐函数定理公理库隐函数定理Implicit function theorem当相关偏导块可逆时,方程组局部可把部分变量表示为其余变量的函数。给出所选解支的局部可解性。作为一步积分法公理库常微分方程时间步进框架Time-stepping method for ODE · ODE integrator · Numerical time integration将常微分方程初值问题离散为时间网格上的数值状态与一步或多步更新,并分离流映射、误差层和退出状态。,在光滑性、轨迹区域和稳定传播条件下,它有全局 误差与局部 缺陷。
对可分系统,动量更新和位置更新各自是 Hamilton 子流,这是按该顺序的Lie–Trotter 分裂公理库Lie–Trotter 与 Strang 分裂Lie-Trotter splitting · Strang splitting把可解子流按一阶或对称二阶次序复合,用非交换误差解释分裂的精度与结构保持。;辛映射公理库辛映射Symplectic map以 Jacobian 的 J 恒等式定义辛映射,证明 Hamilton 流保持它,并区分辛结构、体积与能量。对复合封闭,所以组合保辛。一般形式也可通过生成函数证明保辛。把顺序反过来得到另一种辛 Euler;两者互为伴随,单独一种通常不时间对称。
直觉
显式 Euler 用旧位置推动动量,也用旧动量推动位置,两个变化没有在同一步内相互回应。辛 Euler 让第二次更新使用刚得到的值,形成一个保持面积配对的剪切组合。
这种改变针对几何;方法仍为全局一阶,局部截断误差为 。它适合展示“误差阶数相同,长期行为仍可能不同”。
例子与边界
谐振子的映射与精确离散不变量
对 ,方法写成
矩阵行列式为一,所以在二维中保辛。直接展开还能验证
在每一步精确保持。 时这个二次型正定,因而把轨迹限制在一个椭圆上。它是本线性例子的精确不变量;一般非线性问题的修正 Hamiltonian 需按形式级数和截断误差理解。
取 、初值 ,得到 。真实能量从 变成 ,而 仍为 。显式 Euler 则给出 ,能量为 ,每步能量乘 ,会持续增长。
稳定区与可逆性
辛 Euler 矩阵的迹为 ,仅在 时具有不同的单位圆特征值;边界 可产生线性增长,超过二则有指数不稳定方向。保辛没有取消步长限制。
用同一种辛 Euler 走一步 再走一步 ,一般不回到原点。精确反演需要交换子步顺序。对称组合的构造见 Lie–Trotter 与 Strang 分裂。
推论与应用
可分形式每步各需一次位置力和动量导数评值,向量更新为 ,主要成本由梯度计算决定。不可分形式还要计入隐式求解成本,并控制求解容差。
离散梯度方法公理库离散梯度保能量方法Discrete gradient method用两点链式法则与反对称结构构造精确保能量更新,并在四次势能上解出一个完整隐式步。选择另一项结构目标:在隐式方程精确求解时逐步保持给定 Hamiltonian,却一般不保辛。该页的四次势能步给出了保能量而面积不保的具体对照;辛 Euler 则保辛而未必保持原能量。
比较方法时同时检查状态误差、辛缺陷、往返误差和真实能量轨迹。对一般系统,固定小步长下的结构保持与长时间精度要结合具体频率和轨迹区域判断。
参考资料