时间分数阶微分方程的数值解法

---

文档主要内容

论文 | 适用人群：数学、物理、金融、材料科学等领域的研究人员、工程技术人员及高校相关专业学生。

本文聚焦于时间分数阶微分方程的数值解法，旨在解决该类方程解析解难以直接获取的实际问题。随着计算工具的持续升级，分数阶微分方程在金融经济学、材料科学、遗传学、物理学等多个研究领域的应用日益广泛。然而，由于解析解通常不具备实用性，研究者转而寻求高效、准确的数值求解方法。本文的核心目标正是介绍一种简单且精确的数值算法，并通过MATLAB验证其精度，最终完成误差分析。

全文围绕该算法展开，共分为三个部分。第一部分系统阐述了分数阶微积分方程的定义及相关定理，重点介绍了Caputo导数、Riemann-Liouville积分以及微积分算子的基本性质。第二部分详细推导了数值算法的实施过程，包含两个关键步骤：首先利用微积分算子的相关定理将原方程转化为等价形式，随后采用离散近似方法对等价方程进行离散化处理，从而建立可计算的数值格式。第三部分通过数值实验求解具体的时间分数阶微分方程，以检验算法的有效性，并同步进行误差分析。

数值实验的结果表明，近似解与精确解之间的误差值极小，从而充分验证了该数值方法的有效性与可靠性。这一结论为分数阶微分方程的实际应用提供了坚实的数值基础，尤其适用于需要高精度模拟的复杂系统建模场景。文档中涉及的关键技术包括Caputo导数、Riemann-Liouville积分、数值逼近与线性化，这些内容对于理解分数阶微分方程的数值求解原理具有重要参考价值。

通过本文，读者可以快速掌握一种经过验证的数值解法流程，并直接借鉴其MATLAB实现思路，用于自身研究或工程计算中的误差评估与算法优化。