[论文解读] Approximations for the Caputo derivative (II)

Yuri Dimitrov, Venelin Todorov+1 more|arXiv (Cornell University)|Apr 25, 2016

Fractional Differential Equations Solutions参考文献 20被引用 6

一句话总结

本论文利用多 polylogarithm 展开和中点求积法，为 Caputo 分数阶导数开发了高阶数值逼近方法，推导了对 Beta 函数和幂函数 Caputo 导数的 Riemann 和逼近的渐近误差展开。主要贡献是提出了一类新型的移位 Grünwald 型逼近格式，显著提升了收敛速率，通过求解阶数高达 2−α 和 2−2α 的分数阶微分方程进行了验证，在测试案例中实现了二阶精度。

ABSTRACT

In the present paper we use the expansion formula of the polylogarithm function to construct approximations of the Caputo derivative which are related to the midpoint approximation of the integral in the definition of the Caputo derivative. The asymptotic expansion formula of the Riemann sum approximation of the beta function and the first terms of the expansion formulas of the approximations of the Caputo derivative of the power function are obtained in the paper. The induced shifted approximations of the Grünwald formula and the approximations of the Caputo derivative studied in the first part of the paper are constructed and applied for numerical solution of fractional differential equations.

研究动机与目标

• 开发超越标准 L1 格式的 Caputo 分数阶导数更高阶数值逼近方法。
• 通过构建移位 Grünwald 型公式，提升求解分数阶微分方程的收敛速率。
• 推导 Beta 函数和幂函数 Caputo 导数的 Riemann 和逼近的渐近展开公式。
• 通过已知解析解的两 Term 和三 Term 分数阶微分方程的数值解，验证所提逼近方法。
• 利用最大 l∞ 范数分析新逼近方法的收敛行为与误差阶。

提出的方法

• 利用多 polylogarithm 函数展开，基于积分定义的中点求积法构造 Caputo 导数的逼近格式。
• 推导 Beta 函数 Riemann 和逼近的渐近展开公式，从而实现对数值格式的误差分析。
• 通过修改标准权重，构造 Grünwald 公式的诱导移位逼近格式，以提升收敛阶。
• 将推导出的逼近方法应用于数值求解两 Term 和三 Term 分数阶微分方程，采用迭代格式并进行稳定性与收敛性分析。
• 利用 Euler-Maclaurin 公式和 zeta 函数恒等式，推导涉及解的导数与伯努利数的展开项。
• 通过已知解析形式的分数阶微分方程的精确解验证方法，计算最大 l∞ 误差与收敛阶。

实验结果

研究问题

• RQ1 基于多 polylogarithm 的展开是否能相对于标准 L1 方法，为 Caputo 导数提供更高阶的逼近？
• RQ2 对 Beta 函数和幂函数 Caputo 导数的 Riemann 和逼近的渐近误差结构为何？
• RQ3 移位 Grünwald 型逼近在求解分数阶微分方程时如何提升收敛速率？
• RQ4 当将新数值格式应用于两 Term 和三 Term 分数阶微分方程时，其观测收敛阶是多少？
• RQ5 所推导的展开公式在分数阶导数逼近背景下，与 zeta 函数和伯努利数有何关联？

主要发现

• 所提出的移位逼近格式在两 Term 分数阶微分方程中实现了二阶收敛（收敛阶为 2−α），数值测试中观测收敛阶接近 1.99。
• 对于 α=0.25 的两 Term 方程，当 h=0.003125 时，新格式的最大误差为 0.7875×10⁻⁶，收敛阶达 1.9998。
• 当 α=0.4 时，三 Term 方程的数值解在 h=0.003125 处达到 1.1727 的收敛阶，表明其具有更高的稳定性和精度。
• 当 α=0.6 时，二阶格式实现了 1.9292 的收敛阶，验证了理论上的 2−α 收敛速率。
• 误差分析表明，涉及 Riemann zeta 函数和解的二阶导数的渐近展开项对主导误差项有显著贡献。
• 该方法成功扩展了 L1 逼近框架，通过引入基于多 polylogarithm 和 zeta 函数恒等式的高阶校正项。