形式陈述
能否先离散作用量,再由驻点条件导出更新,而不直接离散运动方程?以下在构形的实局部坐标 中取光滑的 与 ,并固定非零步长 。设 Lagrangian 为 ,连续运动对应其 Euler–Lagrange 方程公理库常微分方程Ordinary differential equation · ODE未知函数及其单一自变量导数组成的方程。。离散 Lagrangian 近似一小段作用量,定义
记 为对第 个端点的多元导数公理库多元函数导数Derivative in several variables · Jacobian derivativeFréchet 导数在有限维 Euclidean 坐标中的 Jacobian、梯度与方向导数表示。。固定两个端点,对内部 求变分,每个 只出现在相邻两项中,得到离散 Euler–Lagrange 方程
若混合导数 非奇异,隐函数定理公理库隐函数定理Implicit function theorem当相关偏导块可逆时,方程组局部可把部分变量表示为其余变量的函数。可在所选解附近局部解出下一点。定义离散动量 、,得到相空间的一步映射。
直觉
每段离散轨道贡献一个作用量。内部节点的微小移动同时改变左段和右段;驻点要求两侧端点动量刚好接上。动力学更新因此来自一项全局变分原则的局部一致性。
对单段有 。再取一次外微分公理库外微分Exterior derivative把 k 形式映为 k+1 形式且满足 d²=0 与分次 Leibniz 规则的算子。, 给出端点辛二形式一致,从而解释正则离散 Legendre 变换下的局部保辛映射公理库辛映射Symplectic map以 Jacobian 的 J 恒等式定义辛映射,证明 Hamilton 流保持它,并区分辛结构、体积与能量。。这里需要非退化条件把关系真正变成局部映射。
例子与边界
最简单的机械离散作用量
取 ,用
离散 Euler–Lagrange 方程给出
定义上述离散动量后,有 、,即
这正是先动量后位置的辛 Euler公理库辛 Euler 方法Symplectic Euler method以更新后的动量推动位置,构造一阶辛方法,并在谐振子上精确计算稳定区和离散不变量。。二阶位置差分的外观并不意味着这套端点相空间变量自动获得二阶一步精度;离散 Lagrangian 的求积选择及动量定义共同决定方法阶数。
若 对平移不变,。对 求导得 ,所以 ,给出离散平移动量守恒。更一般对称性也产生相应离散动量守恒律。
推论与应用
构造流程是选择作用量近似、求两端导数、检查混合导数可逆性、解离散动量方程、核验边界匹配与对称性。若离散作用量保留某种空间对称性,就有机会保留对应动量;粗糙求积或不对称处理可能破坏它。
固定步长的变分积分器通常不精确保留原 Hamiltonian,时间离散打破了连续时间平移的直接结构。保辛、保动量和保能量应分别陈述。退化 Lagrangian、约束和外力也需要相应扩展,不能直接套用正则无力公式。
参考资料
- Marsden and West, Discrete mechanics and variational integrators, Acta Numerica 10, 2001, pp. 357–514, Part 1,离散作用量、Legendre 变换及离散 Noether 定理。
- Hairer, Lubich and Wanner, Geometric Numerical Integration, 2nd ed., Chapter VI,变分形式与辛方法。