Volterra Integral Equations

考虑如下形式的Volterra积分方程：

其中，是三角形区域上的一个给定核函数。下面介绍几种常用来求解此种问题的数值方法。

Collocation Method

REFERENCES
Shen, Jie, Tao Tang, and Li-Lian Wang. Spectral methods: algorithms, analysis and applications. (2011) Section 5.1

配点法的思路非常直接：选定一组配点，要求找到一个多项式使得在这些配点处满足方程，即

然而由于方程中的积分是变上限的，这样的形式不利于构建离散的方程，为了对于近似点值建立线性系统进行求解，我们需要先使用变量代换来将上式的变上限积分转化为固定的上的积分。令

其中而，于是原方程变为

对应的配点方程为

下一步我们需要建立插值矩阵，使用配点网格上的值进行插值来近似，从而将上式中的积分转化为一个矩阵-向量乘积。设

是以配点网格上的Lagrange基函数为基底的插值多项式，则

接下来将上面的表达式代入配点方程中得到

我们可以把这一表达式写为矩阵形式：

其中是第个标准基向量，，。将对应的方程合起来，我们就得到了一个线性系统

其中

一般需要通过数值积分来计算，常用的数值积分方法包括高斯-勒让德积分和Clenshaw-Curtis积分等，因此矩阵的构建复杂度依赖于核函数的光滑性以及数值积分方法的效率，通常复杂度为，其中是计算积分使用的点的数量(相比于一般无法忽略)。

通常实际计算中为了兼顾计算效率与数值稳定性会将配点选为Gauss点，此时基函数可被相应的正交多项式表示，例如当使用Legendre-Gauss配点时，可以将上述积分中的Lagrange基函数表示成Legendre多项式展开的形式：

其中当使用Legendre-Gauss与Legendre-Gauss-Radau配点时，，而当使用Legendre-Gauss-Lobatto配点时，，对于。

最后我们可以通过求解上述线性系统来得到近似解。由于矩阵通常是满矩阵，这类方法的总体复杂度为，其中是配点的数量。

Chebfun目前使用此类方法来求解Volterra积分方程，具体实现细节可以参考Chebfun的文档和源代码。

Olver’s Method

REFERENCES
Gutleb, Timon S., and Sheehan Olver. A sparse spectral method for Volterra integral equations using orthogonal polynomials on the triangle. SINUM (2020).