[论文解读] Alternative Jacobi Polynomials and Orthogonal Exponentials

Vladimir S. Chelyshkov|arXiv (Cornell University)|May 10, 2011

Mathematical functions and polynomials参考文献 3被引用 9

一句话总结

本文引入了在区间 [0,1] 上正交的替代雅可比多项式，以及在 [0,∞) 上正交的指数多项式，采用受雅可比的“倒置，始终倒置”原则启发的双向正交化算法。主要贡献在于建立了递推关系、积分表示以及显式权函数下的正交性，表明这些系统在光谱方法（如流体动力学和逼近理论）中的算法实用性与经典雅可比多项式相当。

ABSTRACT

Sequences of orthogonal polynomials that are alternative to the Jacobi polynomials on the interval $[0,1]$ are defined and their properties are established. An $(α,β)$-parameterized system of orthogonal polynomials of the exponential function on the semi-axis $[0,\infty)$ is presented. Two subsystems of the alternative Jacobi polynomials, as well as orthogonal exponential polynomials are described. Two parameterized systems of discretely almost orthogonal functions on the interval $[0,1]$ are introduced.

研究动机与目标

• 开发一种在区间 [0,1] 上替代经典雅可比多项式的新型正交多项式类，动机源于近壁湍流的光谱方法。
• 构建一个 (α,β)-参数化的半轴 [0,∞) 上的正交指数多项式系统，适用于逼近具有指数衰减特性的函数。
• 引入两组基于正交指数多项式与常数项构造的、在 [0,1] 上离散几乎正交的参数化函数系统，用于 hp-逼近与数值积分。
• 在参数平移下，建立替代雅可比多项式与经典雅可比多项式之间的算法等价性，证明其在数值分析中的实用性。
• 探讨奇异正交性的特殊情况，并通过柯西公式提供积分表示，增强其理论与计算适用性。

提出的方法

• 通过罗德里格斯型公式定义替代雅可比多项式：$\mathcal{P}_{nk}^{(\alpha,\beta)}(x) = \frac{x^{-\alpha-k-1}(1-x)^{-\beta}}{(n-k)!} \frac{d^{n-k}}{dx^{n-k}} \left( x^{\alpha+n+k+1}(1-x)^{\beta+n-k} \right)$。
• 通过分部积分法证明在 [0,1] 上关于权函数 $x^\alpha(1-x)^\beta$ 的正交性，得到 $\int_0^1 x^\alpha(1-x)^\beta \mathcal{P}_{nk}^{(\alpha,\beta)}(x) \mathcal{P}_{nl}^{(\alpha,\beta)}(x) dx = h_{nk}^{(\alpha,\beta)} \delta_{kl}$。
• 通过重复使用分部积分法，推导出范数表达式 $h_{nk}^{(\alpha,\beta)} = \frac{1}{\alpha+2k+1} \frac{\Gamma(\alpha+n+k+2)\Gamma(\beta+n-k+1)}{(n-k)!\Gamma(\alpha+\beta+n+k+2)}$。
• 建立互惠关系：$x^p \mathcal{P}_{nk}^{(\alpha,\beta)}(x) = \mathcal{P}_{n+p,k+p}^{(\alpha-2p,\beta)}(x)$，将不同参数区间的系统联系起来。
• 通过指数函数 $\{e^{-kt}\}$ 的双向正交化，引入 [0,∞) 上的正交指数多项式 $\mathcal{E}_{nk}^{(\alpha,\beta)}(t)$，确保其与雅可比系统的算法等价性。
• 利用缩放指数函数与常数项，构造两组离散几乎正交函数系统 $\boldsymbol{\widetilde{\mathcal{Z}}}_n^{(\omega_n)}(t)$ 和 $\boldsymbol{\mathcal{Z}}_n^{(\omega_n)}(t)$，其中 $\gamma_n$ 控制正交性的质量。

实验结果

研究问题

• RQ1 如何通过双向正交化算法在 [0,1] 上生成与经典雅可比多项式性质相匹配的替代正交多项式？
• RQ2 替代雅可比多项式 $\mathcal{P}_{nk}^{(\alpha,\beta)}(x)$ 关于权函数 $x^\alpha(1-x)^\beta$ 的正交结构与范数为何？
• RQ3 能否系统性地构造出在 [0,∞) 上的正交指数多项式，使其在算法效率上与 [-1,1] 上的雅可比多项式相当？
• RQ4 参数平移与缩放如何影响这些替代系统的正交性与递推结构？
• RQ5 参数 $\gamma_n$ 在构造用于有限区间上 hp-逼近的离散几乎正交函数中的作用是什么？

主要发现

• 替代雅可比多项式 $\mathcal{P}_{nk}^{(\alpha,\beta)}(x)$ 在 [0,1] 上关于权函数 $x^\alpha(1-x)^\beta$ 正交，满足 $\int_0^1 x^\alpha(1-x)^\beta \mathcal{P}_{nk}^{(\alpha,\beta)}(x) \mathcal{P}_{nl}^{(\alpha,\beta)}(x) dx = h_{nk}^{(\alpha,\beta)} \delta_{kl}$。
• 范数显式给出为 $h_{nk}^{(\alpha,\beta)} = \frac{1}{\alpha+2k+1} \frac{\Gamma(\alpha+n+k+2)\Gamma(\beta+n-k+1)}{(n-k)!\Gamma(\alpha+\beta+n+k+2)}$，其中 $\alpha > -1$，$\beta > -1$。
• 该系统满足互惠性质：$x^p \mathcal{P}_{nk}^{(\alpha,\beta)}(x) = \mathcal{P}_{n+p,k+p}^{(\alpha-2p,\beta)}(x)$，将不同参数区域的系统联系起来。
• 在 [0,∞) 上的正交指数多项式 $\mathcal{E}_{nk}^{(\alpha,\beta)}(t)$ 通过双向正交化构造，实现了与雅可比系统的算法等价性。
• 引入两组离散几乎正交函数系统 $\boldsymbol{\widetilde{\mathcal{Z}}}_n^{(\omega_n)}(t)$ 和 $\boldsymbol{\mathcal{Z}}_n^{(\omega_n)}(t)$，其中 $\gamma_n$ 控制正交性的程度。
• 利用柯西积分公式，推导出积分表示 $\mathcal{P}_{nk}^{(\alpha,\beta)}(x) = \frac{1}{2\pi i} \frac{(x^{-1})^{\alpha+\beta+n+2}}{(1-x^{-1})^\beta} \int_C \frac{z^{-(\alpha+\beta+n+k+2)}(1-z)^{\beta+n-k}}{(z-x^{-1})^{n-k+1}} dz$。