[论文解读] Solving generalized Abel's integral equations of the first and second kinds via Taylor-collocation method

Eisa Zarei, Samad Noeiaghdam|arXiv (Cornell University)|Apr 23, 2018

Fractional Differential Equations Solutions参考文献 28被引用 11

一句话总结

本文提出了一种新颖的数值方法，结合泰勒多项式与配点法，求解广义阿贝尔积分方程的第一类和第二类。通过使用泰勒展开近似未知函数，并在离散点上施加配点条件，该方法将积分方程转化为线性系统。数值结果表明，该方法具有高精度和快速收敛性，在测试实例中当 n=9 时，最大误差降至 1.4×10⁻¹³ 以下。

ABSTRACT

The aim of this paper is to present an efficient numerical procedure to approximate the generalized Abel's integral equations of the first and second kinds. For this reason, the Taylor polynomials and the collocation method are applied. Also, the error analysis of presented method is illustrated. Several examples are approximated and the numerical results show the accuracy and efficiency of this method.

研究动机与目标

• 开发一种高效数值格式，用于求解广义阿贝尔积分方程的第一类和第二类。
• 解决具有分数阶奇异性特征的伏尔泰拉型积分方程中奇异核的挑战。
• 通过混合泰勒-配点方法，相较于现有方法，提升求解精度与收敛速度。
• 为所提出的方法提供严格的误差界与收敛性分析。
• 在具有已知精确解的基准问题上验证该方法的有效性。

提出的方法

• 使用以点 z 为中心的 n 次泰勒多项式近似未知函数 Φ(x)。
• 通过选择 n+1 个等距点 xi = a + (b−a)i/n（i=0 到 n）应用配点法。
• 将泰勒展开式代入积分方程，并在每个配点处强制满足方程。
• 构建线性系统 AX = G，其中 A 包含 (t−z)^j / (ϕ(xi)−ϕ(t))^α 的积分，X 包含泰勒系数 Φ^(j)(z)，G 包含 g(xi)。
• 求解线性系统以确定泰勒系数，并重构近似解。
• 利用所得多项式计算不同 x 和 n 值下的数值解与绝对误差。

实验结果

研究问题

• RQ1 泰勒-配点法是否能在求解广义阿贝尔积分方程的第一类和第二类时实现高精度？
• RQ2 所提出方法的收敛速率如何依赖于泰勒多项式次数 n？
• RQ3 对于不同的 n 值和测试问题，最大绝对误差的量级是多少？
• RQ4 该方法在具有已知精确解的问题上的表现如何，特别是在数值一致性方面？
• RQ5 该方法是否能有效处理具有 0 < α < 1 的奇异核，同时保持稳定性和精度？

主要发现

• 泰勒-配点法实现了高精度，当 n=9 时，示例1的最大绝对误差降至 1.37683×10⁻¹³。
• 对于示例2（Φ(x) = cos x），当 n=9 时最大误差降低至 3.06992×10⁻¹⁰，表明收敛速度极快。
• 在示例3（Φ(x) = e^x）中，当 n=9 时最大误差达到 1.28571×10⁻⁸，显示出随着 n 增大，误差持续减小。
• 对于示例4（Φ(x) = x²），当 n=5 时得到精确解；对于示例5（Φ(x) = 1−2x），当 n=3 时也得到精确解，证实该方法对某些多项式解具有精确性。
• 误差分析表明该方法具有稳定性，且误差传播受控，随着 n 增大，误差显著减小。
• 所有示例的数值结果均表明该方法具有鲁棒性、高效性，并且相较于标准方法表现更优。