[论文解读] Weak error expansion of the implicit Euler scheme

Omar Aboura|arXiv (Cornell University)|Apr 25, 2013

Advanced Numerical Methods in Computational Mathematics参考文献 6被引用 3

一句话总结

本文为伊藤随机微分方程应用隐式欧拉格式建立了弱误差展开，将塔拉伊-图巴罗框架扩展至隐式格式。证明了该格式具有弱阶1的收敛性，并推导出弱误差的一阶展开，其中涉及一个新函数 $\psi_i$，其与显式格式的 $\psi_e$ 不同，差异项包含 $b^2\Delta u$ 和 $\sigma^2\Delta(b\partial u)$，其余项有 $O(h^2)$ 的界。该结果依赖于柯尔莫哥洛夫方程、伊藤公式与克拉克-翁孔公式，以及插值技术。

ABSTRACT

In this paper, we extend the Talay Tubaro theorem to the implicit Euler scheme.

研究动机与目标

• 将塔拉伊-图巴罗弱误差展开框架从显式格式扩展至随机微分方程的隐式欧拉格式。
• 分析在光滑有界系数下隐式欧拉格式的弱收敛行为。
• 推导一个精确的一阶弱误差展开，其中涉及一个表征主导误差项的新函数 $\psi_i$。
• 为研究随机偏微分方程中的弱误差提供理论基础，因为在随机偏微分方程中隐式格式更为自然。

提出的方法

• 利用与SDE相关的柯尔莫哥洛夫后向方程推导弱误差展开。
• 引入离散隐式欧拉格式的连续插值，以支持弱误差的分析。
• 应用伊藤公式与克拉克-翁孔公式，对误差项进行分解与估计。
• 使用函数 $S_h(x) = 1/(1 - h b'(x))$ 建模隐式格式的稳定性与校正项。
• 通过求和引理建立离散到连续的逼近，证明 $h\sum E[v(t_{k+1}, X^{N}_{t_{k+1}}) w(t_k, X^{N}_{t_k})] \to E\int v w(t, X_t) dt$，误差为 $O(h)$。
• 比较 $\psi_{ih}$（$\psi_i$ 的离散逼近）与连续形式 $\psi_i$，证明其期望差值为 $O(h)$。

实验结果

研究问题

• RQ1 塔拉伊-图巴罗弱误差展开框架能否被扩展至随机微分方程的隐式欧拉格式？
• RQ2 隐式欧拉格式弱误差展开的主导项是什么？它与显式情况有何不同？
• RQ3 隐式格式的弱误差在标准 $O(h)$ 收敛之外，如何依赖于漂移与扩散系数？
• RQ4 函数 $S_h(x) = 1/(1 - h b'(x))$ 在校正隐式格式弱误差中起什么作用？
• RQ5 隐式欧拉格式的弱误差是否仍为阶1？能否推导出更高阶展开？

主要发现

• 隐式欧拉格式具有弱阶1，满足 $|\mathbb{E}[f(X_T^N)] - \mathbb{E}[f(X_T)]| \leq Ch$，其中常数 $C$ 与 $N$ 无关。
• 弱误差具有如下一阶展开：$\mathbb{E}[f(X_T^N)] - \mathbb{E}[f(X_T)] = h\mathbb{E}\int_0^T \psi_i(t, X_t) dt + O(h^2)$，其中 $\psi_i$ 显式地由 $b$、$\sigma$ 及柯尔莫哥洛夫PDE解 $u$ 的导数定义。
• $\psi_i$ 与显式格式的 $\psi_e$ 的差异项包含 $b^2\Delta u$、$\sigma^2\Delta(b\partial u)$ 及更高阶导数，反映了格式的隐式特性。
• 离散逼近 $\psi_{ih}$ 与连续形式 $\psi_i$ 的差异在期望意义下有 $O(h)$ 的界，确保展开在极限下仍有效。
• 证明依赖于柯尔莫哥洛夫方程、伊藤公式与克拉克-翁孔公式，以及控制离散-连续逼近误差的求和引理。
• 当 $b=0$ 时，$\psi_i = \psi_e = \frac{1}{8}\sigma^4 \partial^4 u - \frac{1}{8}\sigma^2 \Delta(\sigma^2 \Delta u)$，确认与显式情况的一致性。