修正方程与后向误差分析

形式陈述 ​

离散轨迹究竟在近似哪一个邻近微分方程？对自治方程 $y^{\prime }=f(y)$ 和一步方法公理库常微分方程时间步进框架Time-stepping method for ODE · ODE integrator · Numerical time integration将常微分方程初值问题离散为时间网格上的数值状态与一步或多步更新，并分离流映射、误差层和退出状态。 $y_{n+1}={\mathrm{\Phi }}_h(y_n)$，先假设 $f$ 与 $\mathrm{\Phi }(h,y)={\mathrm{\Phi }}_h(y)$ 在 $h=0$ 及所考察状态区域的共同邻域内为 $C^{\mathrm{\infty }}$，并满足 ${\mathrm{\Phi }}_0=\mathrm{id}$、${\partial }_h\mathrm{\Phi }(0,y)=f(y)$。后向误差分析寻找形式修正向量场

$$f_h=f+h{f}_2+h^2{f}_3+\cdots ,$$

使其形式时间 $h$ 流与 ${\mathrm{\Phi }}_h$ 逐阶相同。这里“形式”表示按 $h$ 的幂次比较系数，尚未断言无穷级数收敛。

实际误差界使用有限截断 $f_h^{[r]}$，其中保留至 $h^r$ 项。对固定整数 $r\ge 0$，一个足够的有限阶条件是 $f$ 与 $\mathrm{\Phi }$ 联合具有 $C^{2r+3}$ 正则性，所需导数在共同邻域内有界，并且对足够小的 $h$，截断向量场的时间 $h$ 流在所考察的共同区域内存在。逐阶匹配及 Taylor 余项估计于是给出 ${\mathrm{\Phi }}_h(y)-{\phi }_h^{[r]}(y)=O(h^{r+2})$，常数由这些导数界决定。这是一个方便的充分条件，并非最少导数数目的断言。

当原向量场来自Hamilton 系统公理库Hamilton 系统Hamiltonian system · Hamilton 系统由 Hamiltonian 的梯度经标准辛矩阵生成相空间流的常微分方程系统。且一步映射保辛配对公理库辛映射Symplectic map以 Jacobian 的 J 恒等式定义辛映射，证明 Hamilton 流保持它，并区分辛结构、体积与能量。时，形式修正场在适当局部 Hamilton 表示中仍为 Hamilton 场。一类标准的解析结论适用于解析 Hamiltonian 的固定步长辛（分区）RK 方法：在轨迹所处紧区域的共同复邻域中有统一解析界，步长相对于相关频率足够小，并选择与 $1/h$ 同阶的截断阶，可得到指数小的单步缺陷和长时间近似能量保持。更一般的方法须另行核验相应的解析系数与截断估计。一般结论是“截断后极接近”，并不是存在一个收敛的修正 Hamiltonian 让任意离散轨迹永久精确等于它的流。

直觉

前向分析固定原方程，问数值点离真轨迹多远；后向分析稍改方程，问这串点可以由怎样的连续运动解释。修改的是动力学，形式上类似把算法的系统性偏差写入向量场。

对显式 Euler，${\mathrm{\Phi }}_h(y)=y+hf(y)$。假设修正场为 $f+hg+O(h^2)$，其时间 $h$ 流由Taylor 展开公理库多元 Taylor 定理Multivariable Taylor theorem用各阶导数给出多元函数的局部多项式展开及余项控制。给出

$$y+hf+h^2\{g+\frac{1}{2}{f}^{\prime }f\}+O(h^3).$$

为匹配 Euler 缺少的二阶项，应取 $g=-\frac{1}{2}{f}^{\prime }f$。这个计算揭示算法带来的首阶动力学扰动，且不需要声称整个修正级数收敛。

例子与边界

线性标量的特殊精确嵌入 ​

对实标量方程 $y^{\prime }=\lambda y$，Euler 一步乘 $1+h\lambda$。当 $h\ne 0$ 且 $1+h\lambda >0$ 时，恰好是常系数方程

$$y^{\prime }=\frac{\log (1+h\lambda )}{h}y=(\lambda -\frac{h}{2}{\lambda }^2+\frac{h^2}{3}{\lambda }^3-\cdots )y$$

的一步精确流。幂级数在 $|h\lambda |<1$ 收敛，这比对数系数作为实数存在的范围更窄。线性例子确实可精确嵌入，但不能把它推广成所有非线性方法的一般事实。

若取 $1+h\lambda =-1$，一步映射为 $y\mapsto -y$。一维光滑自治实微分方程的局部流保持次序，不能实现这种反向映射，展示了步长与嵌入条件的真实限制。

能量界的正确读法 ​

对阶数 $p$ 的辛方法，修正 Hamiltonian 通常从 $H+O(h^p)$ 开始。在上述解析、有界、固定小步长条件下，典型估计具有 $|H(y_n)-H(y_0)|\le C{h}^p+C{t}_n{e}^{-c/h}$ 的结构。这解释较长时间内的小幅能量误差，同时保留截断缺陷累积项。

步长接近高频稳定边界、轨迹接近奇点，或采用一般状态依赖步长时，统一解析界或保辛结构可能失效。长期结论应在实际参数范围内重新检查。

推论与应用

辛 Euler公理库辛 Euler 方法Symplectic Euler method以更新后的动量推动位置，构造一阶辛方法，并在谐振子上精确计算稳定区和离散不变量。在谐振子上有可显式验证的精确二次不变量；一般非线性系统则依赖本页的形式与截断分析。两种陈述的强度不同，报告时应保留这个区分。

后向误差可以解释能量表现，却不保证长时间位置与相位都准确。轨道可能保持在正确附近的能量曲线上，同时逐渐积累相位误差。

参考资料

• Ernst Hairer, Challenges in Geometric Numerical Integration, Trends in Contemporary Mathematics, Springer INdAM Series 8, 2014, pp. 125–135；§2 中的 “Backward Error Analysis — General Situation”，截断修正方程与误差。
• Hairer, Lubich and Wanner, Geometric Numerical Integration, 2nd ed., Chapter IX，形式级数、最优截断和长时间能量估计。