组合积分器的升阶 | Codeboy's Blog

形式陈述

已有一个对称二阶方法，怎样只调整子步长度就升到四阶？对自治系统，设 $S_h$ 是一致、对称、足够光滑的二阶一步法（非自治问题须把时间并入状态），例如Strang 分裂。构造

\Psi_h=S_{ah}\circ S_{bh}\circ S_{ah},\qquad 2a+b=1,\quad2a^3+b^3=0.

第一式保证推进的净时间是一整步，第二式消去领先三阶局部误差。解为

a=\frac1{2-2^{1/3}}\approx1.351207, \qquad b=-\frac{2^{1/3}}{2-2^{1/3}}\approx-1.702414.

在各正负子步均有定义、误差展开有效且稳定的条件下， $\Psi_h$ 具有固定终点上的四阶全局离散化精度，即误差为 $O(h^4)$ ；这不是相邻迭代误差的 Q-阶。如果 $S_h$ 保辛，三个子步的复合仍保辛。

直觉

对称方法的形式对数可写为 $\log S_h=hF+h^3E_3+h^5E_5+\cdots$ 。三个子步中，净流项系数为 $2a+b$ ，领先误差系数为 $2a^3+b^3$ 。对称排列继续消除偶次对数项，所以消掉三阶后，下一个局部误差从五阶开始。

负子步是这次消差的必要部分：若 $a,b$ 都非负且净时间为一，立方和就不可能为零。算法先向前超过一整步，再向后，最后向前结束；中间状态可能越出最终一步所在的局部区域。

例子与边界

对 $y'=y$ ，隐式中点的放大因子是 $r(h)=(1+h/2)/(1-h/2)$ ，其对数展开为

\log r(h)=h+\frac{h^3}{12}+\frac{h^5}{80}+O(h^7).

组合后放大因子为 $r(ah)^2r(bh)$ ，取对数得到

h+\frac{2a^3+b^3}{12}h^3+\frac{2a^5+b^5}{80}h^5+O(h^7).

三阶项为零，五阶系数一般非零，因此这是四阶而非无限阶。用它积分到固定时间，减半步长后全局误差在渐近区应约缩为 $1/16$ ；普通中点约缩为 $1/4$ 。

对热方程的高频模态，正时间子流乘 $e^{-k^2h}$ ，负时间子流却乘 $e^{k^2|h|}$ 。高频噪声会被指数放大，反向扩散在连续问题中本来就不适定。Hamilton 轨道可逆，不代表任意耗散子系统也适合这种负步组合。

同样，带碰撞、终止事件或只在正时间定义的子流程，可能没有可用逆流。即使所有子流存在，较大的中间步系数也可能缩小稳定区或放大舍入。

推论与应用

实现先核验系数的一阶与三阶恒等式，再测试时间对称性和误差阶数。基本方法每步成本为 $C$ 时，这个三重组合约需 $3C$ ，相邻相同子流可合并的情况下略省。比较效率应看同一误差目标下的总成本，而非只比较阶数。

更高阶组合继续消去高阶误差项，通常引入更多子步与更复杂的交换子条件。低阶公式的推导机制可复用，具体系数需要逐阶验证。

参考资料

Haruo Yoshida, Construction of higher order symplectic integrators, Physics Letters A 150, 1990, pp. 262–268，对称组合的升阶构造。
Hairer, Lubich and Wanner, Geometric Numerical Integration, 2nd ed., Chapter II，对称方法的组合与负子步。