A numerical Bernstein splines approach for nonlinear initial value problems with Hilfer fractional derivative

技术摘要：一种用于非线性 Hilfer 分数阶初值问题的数值 Bernstein 样条法

问题陈述
本文研究了涉及 Hilfer 分数阶导数的非线性分数阶初值问题（FIVPs）的数值近似。Hilfer 导数由阶数 $\alpha \in (0,1)$ 和类型 $\beta \in [0,1]$ 参数化，它在 Riemann-Liouville（$\beta=0$）和 Caputo（$\beta=1$）导数之间进行了插值。一个主要的挑战在于，当 $\beta \neq 1$ 且初始条件为非齐次时，解在 $t \downarrow 0$ 时表现出奇异行为（具体而言，它们是非连续的，且其行为类似于 $t^{\gamma-1}$，其中 $\gamma = \alpha + \beta(1-\alpha)$）。现有的数值方法，例如广泛使用的 Adams–Bashforth–Moulton 预测器-校正器方案，通常是为 Caputo 情况量身定制的，或者要求解具有连续可微性，这限制了它们在处理一般 Hilfer 情况下的奇异解时的适用性。此外，现有文献缺乏能够处理非线性 Hilfer FIVP 全范围奇异解且具备理论误差界的有效方法。

方法论
作者提出了一种基于 Bernstein 样条 的数值方法，用于近似求解如下非线性 FIVP：
$$\begin{cases} {}_{0}D^{\alpha,\beta}_{t} y(t) = f(t, y(t)), & 0 < t \leq T \\ {}_{0}I^{1-\gamma}_{t} y(0) = \tilde{y}_0 \end{cases}$$
该方法论包含以下几个关键组成部分：

1. 加权空间变换： 为了处理 $t=0$ 处的奇异性，问题被转化为加权连续函数空间 $C^{1-\gamma}[0, T]$。此外，引入了一个时间域偏移 $\varepsilon > 0$，在 $[\varepsilon, T]$ 上求解该问题（此时解是连续的），然后取 $\varepsilon \downarrow 0$ 的极限。
2. Bernstein 样条近似： 解使用定义在节点集上的分段 Bernstein 多项式（样条）进行近似。该方法采用了局部“逐节点”（knot-for-knot）迭代策略。在每个子区间上，通过将 Hilfer 导数定义的等效积分方程应用的 Picard 迭代来构造解。
3. 算子矩阵实现： 作者利用不完全 Beta 函数推导出了 Bernstein 样条分数阶积分的闭式表达式。这使得分数阶积分算子能够对样条基函数进行解析应用，从而避免了与解的历史相关的数值求积误差。
4. 向量化与并行化计算： 实现过程利用张量积和向量化技术（具体使用 Python 中的 numpy.einsum）来高效计算分数阶积分。这种方法允许针对多项式阶数 $q$ 进行并行化计算，从而减轻了与分数阶导数“记忆效应”相关的计算成本。

核心贡献

• 理论框架： 本文在放宽的收敛假设下，建立了样条近似的存在性与唯一性结果（定理 3.1）。文中提供了关于近似误差相对于多项式阶数 $q$（$O(q^{-\alpha/2})$）、节点大小 $h$（$O(h^\alpha)$）以及时间偏移参数 $\varepsilon$（$O(\varepsilon^\alpha)$）的显式收敛速率（定理 3.2）。
• 处理奇异性： 不同于以往受限于连续解或 Caputo 情况的研究，本方法通过利用加权空间变换和 $\varepsilon$-偏移，专门设计用于捕捉 $\beta \in [0, 1)$ 时解的奇异行为。
• 高效实现： 作者提出了一种向量化算法，在多项式基函数评估方面将计算复杂度与节点大小解耦，这为相比于传统的逐步法实现加速提供了一条极具前景的途径。
• 非线性应用： 该方法成功应用于非线性系统——分数阶 Van der Pol 振荡器，证明了其处理非线性项和变化的 $\beta$ 值的能力。

数值结果

• 收敛性验证： 在多项式测试问题上的数值实验证实了理论收敛速率。随着节点大小 $h$ 和偏移 $\varepsilon$ 的减小，以及多项式阶数 $q$ 的增加（直到达到数值稳定性限制进一步提升的程度），误差随之降低。
• 与预测器-校正器对比： 在线性 Caputo 情况（$\beta=1$）下与 Adams–Bashforth–Moulton 预测器-校正器方法（Diethelm 等人）进行对比时，Bernstein 样条法在相同的节点大小时展现出显著更高的精度（较低的平均误差和上确界误差）。在测试设置中，预测器-校正器方法表现出持续的偏差，并且在更细的网格下并未改善其上确界误差。
• Van der Pol 振荡器： 对 Hilfer 分数阶 Van der Pol 振荡器的模拟表明，无论 Hilfer 类型 $\beta$ 如何，该方法都能收敛到稳定的极限环。结果表明，该系统表现出对唯一极限环的强吸引性，这与 Caputo 情况下的半解析发现一致。

意义与主张
本文声称提供了第一种能够获得非线性 Hilfer 分数阶初值问题全范围解（包括奇异解）的高效数值方法，并具备理论误差界。作者强调，通过使用局部迭代，他们的方法消除了其他方法中存在的对时间域区间大小的严格约束。虽然作者承认样条法目前比预测器-校正器法具有更高的计算需求，但他们强调了其卓越的精度，以及处理标准方法失效或缺乏理论保证的奇异非线性系统的能力。这项工作为未来扩展到边值问题和分数阶偏微分方程（FPDEs）奠定了基础。