De la Vallée Poussin type approximation for solving some Fredholm integral equations

技术摘要：用于求解某些 Fredholm 积分方程的 de la Vallée Poussin 型逼近法

问题陈述
本文研究了第二类 Fredholm 积分方程（FIEs）的数值解，其形式为：
$$f(y) - \nu \int_{-1}^{1} f(x)k(x, y)w(x) dx = g(y), \quad y \in (-1, 1),$$
其中 $w(x)$ 是 Jacobi 权重，$g$ 和 $k$ 是已知函数，$f$ 是未知解。本研究重点关注已知函数可能表现出由 Jacobi 权重 $u$ 控制的端点奇异性的情形，该权重 $u$ 可能与定义积分算子的权重 $w$ 不同。目标是在局部连续函数空间 $C_u$（配备加权一致范数）中求解这些方程，并处理涉及弱奇异核（代数型或对数型）以及高度振荡核的情况。

方法论
作者提出了一种名为 VP-method 的数值方法，该方法基于在 Jacobi 零点处的 de la Vallée Poussin (VP) 多项式逼近。与使用单一节点集构造唯一插值多项式的经典 Lagrange 插值不同，VP 算子 $V_n^m(v_{\alpha,\beta})$ 是一个准投影算子（quasi-projector），它依赖于一个额外的参数 $m$（其中 $1 < m < n$）。

该方法步骤如下：

1. 逼近： 使用与 Jacobi 权重 $v_{\alpha,\beta}$ 相关的 VP 算子 $V_n^m$ 对积分算子 $K$ 和源函数 $g$ 进行逼近。该权重 $v_{\alpha,\beta}$ 的选择独立于积分权重 $w$，以优化逼近性质。
2. 离散化： 将连续方程替换为有限维方程 $(I - \nu K_n^m)f_n^m = g_n^m$，其中 $f_n^m$ 是特定子空间 $S_n^m$ 中的 $n$ 维多项式。
3. 线性系统： 将问题简化为求解一个 $n$ 阶线性方程组。系统矩阵涉及“修正矩”（即核函数对 Jacobi 多项式的积分）。作者提供了高效计算这些矩的递推关系，适用于弱奇异、对数型及高度振荡型核。
4. 比较： 将该方法与 L-method 进行比较，后者是基于文献 [5] 中提出的相同节点处的 Lagrange 插值的投影法。

主要贡献

• 理论稳定性与收敛性： 本文证明，在对权重和核函数做出适当假设的前提下，VP-method 在加权一致空间 $C_u$ 中是稳定且收敛的。至关重要的是，VP 算子的 Lebesgue 常数相对于 $n$ 是一致有界的，而 Lagrange 插值的 Lebesgue 常数至少以 $\log n$ 的速度增长。
• 收敛阶： VP-method 的误差估计被证明与最佳多项式逼近阶一致，即 $O(E_{n-m}(f)_u)$。相比之下，L-method 要求解属于 Hölder–Zygmund 子空间，并在误差估计中引入额外的 $\log n$ 因子以确保收敛。
• 泛用性： VP-method 允许逼近权重 $v_{\alpha,\beta}$ 不同于积分权重 $w$。这种灵活性显著扩大了适用权重函数和解行为的范围，优于通常要求 $v_{\alpha,\beta} = w$ 且对权重施加更严格条件的 L-method。
• 处理奇异性与振荡： 作者确立了在加权空间中，对于弱奇异和高度振荡核，积分算子具有紧性的充分条件。他们还推导了算法所需的特定修正矩递推关系，以应对这些挑战性情况。
• 条件数： VP-method 所产生的线性系统表现出良好的条件特性，其条件数随着离散维度 $n$ 的增加保持有界。

结果
针对以下各类测试用例进行了数值实验：

• 弱奇异核（代数型和对数型）。
• 高度振荡核。
• 具有端点奇异性和剧烈变化的函数。

实验结果证实了理论误差估计。主要发现包括：

• 全局精度： 当两者均适用时，VP-method 和 L-method 的最大绝对误差通常处于同一数量级。
• 局部精度： VP-method 始终比 L-method 获得更高的局部（逐点）精度。这归功于 VP 算子减轻 Gibbs 现象的能力，特别是在函数变化剧烈或存在奇异性的区域附近。
• 参数敏感性： VP-method 的性能取决于参数 $\theta = m/n$ 的选择。虽然没有提供关于最优 $\theta$ 的通用理论准则，但数值测试表明，接近 1 的值（趋向于 Fejér 算子）往往会降低性能，而中间值通常能产生最佳结果。
• 鲁棒性： 在 L-method 无法满足收敛假设的情形下（例如特定的权重 $u$ 和 $w$ 组合），VP-method 仍能成功收敛。

意义与主张
本文声称，VP-method 为求解 Fredholm 积分方程提供了一种优于经典基于 Lagrange 投影方法的稳健替代方案。其主要意义在于：

1. 消除对数增长： 通过利用一致有界的 Lebesgue 常数，该方法避免了与 Lagrange 插值相关的 $\log n$ 误差估计退化和稳定性问题。
2. 灵活性： 将逼近权重与积分权重解耦的能力，使其能够处理更广泛的问题类，包括那些经典方法可能失效的带有奇异性的问题。
3. 卓越的逐点行为： 该方法通过有效抑制 Gibbs 型振荡，在关键区域（如奇异点或快速变化处）提供更精确的解重构。
4. 计算效率： 该方法产生的线性系统条件良好，并利用高效的矩计算递推关系，使其适用于大规模计算。

作者总结道，VP-method 结合了优异的逼近性质、稳定性和计算效率，使其成为求解第二类 Fredholm 积分方程（特别是在存在奇异性和振荡行为时）的有力工具。