[论文解读] Recurrence relations for orthogonal polynomials on a triangle

Sheehan Olver, Alex Townsend+1 more|arXiv (Cornell University)|Jan 27, 2018

Aerospace Engineering and Control Systems参考文献 3被引用 7

一句话总结

本文提出了一套系统性框架，利用提升算符推导三角形上正交多项式（特别是Koornwinder IV类多项式）的稀疏递推关系。通过映射多项式在参数和次数上增减±1，该方法使在三角形区域上求解PDE的谱方法得以高效实现，借助稀疏线性代数，类似于超球谱方法在ODE中的应用。

ABSTRACT

This paper derives sparse recurrence relations between orthogonal polynomials on a triangle and their partial derivatives, which are analogous to recurrence relations for Jacobi polynomials. We derive these recurrences in a systematic fashion by introducing ladder operators that map an orthogonal polynomial to another by incrementing or decrementing its associated parameters by one.

研究动机与目标

• 开发一种系统性方法，推导三角形上正交多项式的递推关系，扩展雅可比多项式递推的结构。
• 解决在三角形区域上偏微分方程（PDE）的稀疏谱离散化问题，其中标准方法缺乏稀疏性。
• 将提升算符技术从单变量推广至双变量正交多项式，特别是三角形上的Koornwinder多项式。
• 通过推导乘法、微分及多项式基转换的递推关系，实现PDE的高效求解，借助稀疏线性代数。

提出的方法

• 引入24个提升算符，用于将双变量Koornwinder多项式的参数（a, b, c, d）和次数（n, k）分别增减±1。
• 利用这些提升算符，以移位的Koornwinder多项式表示x、y和z（其中z = 1 - x - y）乘法的递推关系。
• 通过组合提升算符构造微分特征值问题，证明多项式是二阶PDE的特征函数。
• 通过涉及移位雅可比多项式导数的链式法则恒等式，表达多项式的偏导数。
• 利用递推关系的稀疏性，实现三角形区域上PDE的高效谱方法，类比于超球谱方法在ODE中的应用。
• 通过正交多项式的仿射不变性，将结果从直角三角形推广至任意三角形。

实验结果

研究问题

• RQ1 如何利用基于算符的方法系统性地推导出三角形上双变量正交多项式的递推关系？
• RQ2 以移位参数的Koornwinder多项式表示时，x、y和z乘法的递推关系结构是什么？
• RQ3 能否利用提升算符推导出这些多项式所满足的微分特征值问题？
• RQ4 所推导的递推关系如何实现三角形区域上PDE的稀疏谱离散化？
• RQ5 第四参数d在系统性推导此前未报告的递推关系中起到什么作用？

主要发现

• 本文推导出12个雅可比多项式的提升算符和24个双变量Koornwinder多项式的提升算符，实现了参数与次数的系统性平移。
• 推导出x、y和z乘法的稀疏递推关系，仅涉及少数具有移位参数的邻近多项式项。
• x和y的递推关系分别涉及最多四项，参数为(a-1,b,c)和(a,b-1,c)，并给出了显式系数。
• 证明多项式是两个二阶微分算符的特征函数，对应的特征值为 -k(k+b+c+1) 和 -n(n+a+b+c+2)。
• 通过链式法则恒等式，用移位雅可比多项式的导数表达双变量多项式的偏导数。
• 该框架使利用稀疏线性代数在三角形上高效求解PDE成为可能，将超球谱方法推广至三角形区域。