?-Caputo分数阶非线性微分方程在分级网格上的高阶数值方法

《Mathematics》：Higher-Order Numerical Methods for Solving ?-Caputo Fractional Nonlinear Differential Equation with Graded Meshes

　　

该研究针对一类含关于另一函数的Caputo分数阶导数（?-Caputo fractional derivative）的非线性分数阶微分方程，构造了高阶数值方法。通过对分数阶微分算子作变量变换，将原控制问题转化为弱奇异Volterra积分方程。在作变量代换x=?(t)后，关于?的分数次积分项被重新表述为变换变量下的标准分数次积分。研究人员随后在分级网格（graded meshes）上针对分数次积分算子构造二次与三次Lagrange插值逼近，并在初值点附近及最后一个子区间引入特殊处理，使所得格式可显式实现。在允许弱奇异的适当正则性假设下，推导了详细的误差估计。理论结果表明，选取合适的分级参数后，所提方法可恢复期望收敛阶??(N?(3+α))与??(N?4)，其中α∈(0,1)为分数阶导数阶数。若干数值算例验证了理论误差估计并展示了所构格式的有效性。

研究背景方面，分数阶微分方程因能刻画异常扩散、黏弹性、输运过程等复杂动力学系统中的记忆与遗传效应而被广泛采用。相比整数阶模型，分数阶模型含弱奇异卷积核，对非局部时间行为建模更灵活。经典Caputo分数阶非线性方程的标准策略是化为弱奇异核的第二类Volterra积分方程，再施以Adams型或配置型离散；但在均匀网格上，解在初值附近的弱奇异会导致收敛阶下降。分级网格结合高阶插值是补偿初值奇异性、恢复最优收敛阶的有效途径。现有文献对经典Caputo情形已有较系统的高阶格式与误差分析，但对?-Caputo（即导数关于严格增函数?(t)定义）这一更一般算子的非线性初值问题，可用数值方法仍较少，且直接离散系数复杂、单调性难以验证。因此，研究人员开展本研究，旨在为?-Caputo分数阶非线性微分方程建立基于分级网格与Lagrange插值的高阶显式格式并完成严格误差分析。论文发表于《Mathematics》。

关键技术方法上，研究人员先将原?-Caputo初值问题等价转化为关于变换变量x=?(t)的标准弱奇异Volterra积分方程，使核不再显含?；随后在x变量下构造分级网格t_n=?^?1((n/N)^rT)（r为分级参数），对积分项分别采用分段二次与三次Lagrange插值，并在首段用常数/线性起步、末子区间用向后外推以保证显式性；误差分析依赖插值余项估计、?分级网格下的步长估计、以及弱奇异离散Gronwall不等式（基于Mittag–Leffler函数）。数值验证选取?(t)=t（经典Caputo）与?(t)=t+?t2两类算例，在双精度下计算最大节点误差与实验收敛阶。

研究结果部分，第二节“Problem Formulation”给出?-Caputo导数与积分定义，导出等价Volterra方程（式4），引入x=?(t)变换与分级网格（式5），定义变换后非线性项G，并给出解在变换变量下的正则性假设（Assumption 1）。第三节“Quadratic Interpolation Method”构造二次Lagrange插值格式：首点常数、次点线性、其后二次，末区间奇数/偶数n分别处理并给出显式权重；通过Lemma 1证得局部截断误差界??(N^?(3+α))，Lemma 2界住历史系数，Lemma 3给出?分级网格上的离散Gronwall不等式，最终Theorem 1得全局误差??(N^?(3+α))。第四节“Cubic Interpolation Method”将插值升为三次（首段线性、次段二次、其后三次，末段向后三次外推），Lemma 4证截断误差??(N^?4)，Lemma 5界住系数，Theorem 2得全局误差??(N^?4)（α<1时），α→1临界情形带对数因子。第五节“Numerical Simulations”分5.1经典?(t)=t与5.2非平凡?(t)=t+?t2两组：前者取解析解y(t)=t2E_α,3(t^α)型构造，后者取y(t)=E_α(?(t)^α)；权重中分数矩用不完全Beta函数I_z(a,b)计算以避免近消失去。表1–4显示二次法EOC≈3+α、三次法EOC≈4，与理论吻合。第六节“Conclusions”总结：二次法收敛阶??(N^?(3+α))，三次法在α<1时??(N^?4)、α=1临界带对数因子；分级网格有效补偿初值弱奇异，方法不限于经典Caputo。讨论指出当前限制为标量初值、?严格增、f Lipschitz及变换后源项G的正则假设；边界值问题与α≥1情形需另述可解性与稳定性，留待未来工作。

研究结论翻译：本文针对含?-Caputo分数阶导数的非线性微分方程，在分级网格上发展了二次与三次Lagrange插值方法。通过将原问题转化为等价弱奇异Volterra积分方程，在低正则性假设下建立了严格收敛估计。二次法收敛阶为??(N^?(3+α))，三次法在α<1时恢复四阶收敛，临界α=1时带对数因子。对经典?(t)=t与非平凡?(t)=t+?t2的数值实验均与理论估计相符。结果表明所提方法的高阶精度不局限于经典Caputo情形，分级网格能有效补偿解在初值处的弱奇异性。本分析限于标量初值问题、容许的严格增?、Lipschitz非线性及变换后源项G的正则假设；分数阶边值问题与α≥1情形将导致不同积分表述，需补充可解性、稳定性与误差分析，扩展至边值及?-Caputo分数阶偏微分方程为后续工作。