[论文解读] Shifted Chebyshev polynomials for Solving Three-Dimensional Volterra Integral Equations of the second kind

Doaa Mohamed|arXiv (Cornell University)|Sep 27, 2016

Fractional Differential Equations Solutions参考文献 9被引用 5

一句话总结

本文提出了一种基于平移切比雪夫多项式的数值方法，用于求解三维第二类伏尔泰拉积分方程。通过将解近似为张量积平移切比雪夫基函数的截断级数，并在高斯-切比雪夫-洛巴托点处进行配点，该方法将积分方程转化为代数方程组，在测试的线性和非线性情形下均实现了高精度的精确解。

ABSTRACT

In this paper, an efficient method is presented for solving three dimensional Volterra integral equations of the second kind with continuous kernel. Shifted Chebyshev polynomial is applied to approximate a solution for these integral equations. This method transforms the integral equation to algebraic equations with unknown Chebyshev coefficients. Numerical results are calculated and the estimated error in each example is computed using Maple 17.

研究动机与目标

• 开发一种高效数值方法，用于求解具有连续核的三维第二类伏尔泰拉积分方程。
• 克服在应用数学与工程中获得多维积分方程解析解的困难。
• 应用平移切比雪夫多项式作为谱方法，对未知解进行高精度逼近。
• 在具有已知精确解的线性和非线性测试问题上验证该方法。
• 通过数值示例与误差分析，展示该方法的收敛性与精度。

提出的方法

• 将未知解 $ u(x,y,z) $ 近似为 $ x $、$ y $ 和 $ z $ 的平移切比雪夫多项式截断三重级数。
• 使用高斯-切比雪夫-洛巴托配点对定义域进行离散化，并在这些点处强制使残差方程为零。
• 将多项式逼近代入原积分方程，推导出涉及核函数与基函数积分的残差方程。
• 将残差方程转化为关于未知切比雪夫系数 $ a_{i,j,k} $ 的大型线性或非线性代数方程组。
• 使用标准线性代数求解器数值求解所得方程组，以确定系数。
• 利用计算得到的系数与基函数重构近似解。

实验结果

研究问题

• RQ1 平移切比雪夫多项式能否有效逼近三维第二类伏尔泰拉积分方程的解？
• RQ2 基于平移切比雪夫多项式的配点方法在求解线性和非线性三维伏尔泰拉积分方程时的精度如何？
• RQ3 该方法在不同测试问题中的收敛行为与误差量级如何？
• RQ4 当已知精确解时，该方法能否准确复现，特别是对于多项式或指数形式的精确解？
• RQ5 随着多项式逼近次数 $ N $ 的增加，该方法在计算精度方面表现如何？

主要发现

• 在例3.1中，当 $ N=1 $ 时，该方法成功恢复了精确解 $ u(x,y,z) = x+y+z $，所有计算系数均与解析值完全一致。
• 在例3.2中，当 $ N=2 $ 时，该方法精确复现了解 $ u(x,y,z) = x^2y + yz^2 + xyz $，证实了其高精度。
• 在非线性例3.3中，该方法在 $ N=1 $ 时即达到精确解 $ u(x,y,z) = xyz $，表明其对非线性核具有鲁棒性。
• 在例3.4中，最大绝对误差在 $ (0.1,0.1,0.1) $ 处为 $ 2.299 \times 10^{-4} $，在 $ (0.001,0.001,0.001) $ 处减小至 $ 2.99 \times 10^{-8} $，表明随着网格细化，解具有收敛性。
• 在例3.5中，绝对误差从 $ (0.1,0.1,0.1) $ 处的 $ 0.033 $ 降低至 $ (0.001,0.001,0.001) $ 处的 $ 0.00037 $，表明具有良好的收敛性与稳定性。
• 对于具有已知精确解的问题，该方法始终达到机器精度，验证了其可靠性和高效性。