[论文解读] Orthogonal polynomials in and on a quadratic surface of revolution

Sheehan Olver, Yuan Xu|arXiv (Cornell University)|Jun 25, 2019

Electromagnetic Scattering and Analysis参考文献 8被引用 4

一句话总结

本文为 R^{d+1} 中的二次旋转曲面（如圆锥、双曲面和抛物面）上的函数构造了显式的正交多项式基——包括曲面及其内部区域，推广了球谐函数与球内正交多项式。通过快速变换与积分求积规则，实现了在这些非标准区域上的高效高精度函数逼近与谱方法。

ABSTRACT

We present explicit constructions of orthogonal polynomials inside quadratic bodies of revolution, including cones, hyperboloids, and paraboloids. We also construct orthogonal polynomials on the surface of quadratic surfaces of revolution, generalizing spherical harmonics to the surface of a cone, hyperboloid, and paraboloid. We use this construction to develop cubature and fast approximation methods.

研究动机与目标

• 将球谐函数与球内经典正交多项式推广至 R^{d+1} 中的二次旋转曲面（d ≥ 2）。
• 为这些曲面及其边界所围成的体积区域上的函数构造显式的正交多项式基。
• 利用谱方法实现对这些区域上光滑函数的高精度逼近。
• 构建高效的积分求积规则与正交展开的快速算法，适用于这些非典型几何结构。
• 将近期快速变换技术（如 Slevinsky 的方法）拓展至非球形二次域，以支持实际计算。

提出的方法

• 在二次旋转曲面上定义权函数，并通过分离变量法与特殊函数理论构造正交多项式基。
• 在柱坐标 (r, θ, t) 下使用张量积基表示圆锥等曲面上的函数，利用雅可比多项式与切比雪夫多项式。
• 应用 Slevinsky 的快速正交多项式变换算法（适用于圆盘与球面），实现函数的高效重展开，复杂度接近最优。
• 基于高斯-雅可比与高斯-切比雪夫点的张量积网格构造积分求积规则。
• 利用线性代数方法求解张量积基与正交多项式基之间的变换矩阵的逆，确保多项式展开的精确恢复。
• 利用函数的解析性特征预测正交展开中系数的衰减速率，以奇点位置作为收敛速度的预测指标。

实验结果

研究问题

• RQ1 如何显式构造在圆锥、双曲面等二次旋转曲面及其内部区域上的正交多项式？
• RQ2 在这些曲面上，对于一般权函数，正交多项式基的结构是什么？其如何推广球谐函数？
• RQ3 能否利用现有的谱变换技术，将正交展开的快速算法适配至非球形二次域？
• RQ4 在这些曲面上，解析函数的展开系数衰减速率如何？其与奇点位置有何关联？
• RQ5 如何基于正交多项式基构造高效积分求积规则，以实现对二次旋转体的积分与逼近？

主要发现

• 为 R^{d+1} 中圆锥、双曲面与抛物面的曲面及其内部区域上的函数构造了显式的正交多项式基。
• 圆锥曲面上的正交多项式推广了球谐函数，其构造基于连带勒让德多项式与雅可比多项式。
• 对于圆锥上的光滑函数，其正交展开系数呈指数衰减，衰减速率由最近奇点的位置决定。
• 在圆锥上关于 x, y, t 的整函数表现出超指数级系数衰减，表明具有谱收敛性。
• 基于 Slevinsky 算法的快速变换方法实现了接近最优的复杂度，支持正交展开的高效计算。
• 基于张量积积分网格的积分求积规则被推导并验证在这些区域上的积分与逼近中具有高效性。