[论文解读] A note on parameter derivatives of classical orthogonal polynomials

Radosław Szmytkowski|arXiv (Cornell University)|Jan 17, 2009

Mathematical functions and polynomials参考文献 5被引用 15

一句话总结

本文提出了一种计算经典正交多项式——勒让德、盖根堡和广义拉盖尔多项式——参数导数的新方法，通过推导如下形式的展开式：∂Pₙ(λ;z)/∂λ = Σ aₙₖ(λ)Pₖ(λ;z)。利用微分方程和克里斯托夫-达布瓦公式，作者重新推导了勒让德和盖根堡多项式已知的结果，并首次以新方法推导出广义拉盖尔多项式简单展开式，其方法与弗罗利希和科夫此前的工作有本质不同。

ABSTRACT

Coefficients in the expansions of the form $\partial P_{n}(λ;z)}/\partialλ=\sum_{k=0}^{n}a_{nk}(λ)P_{k}(λ;z)$, where $P_{n}(λ;z)$ is the $n$th classical (the generalized Laguerre, Gegenbauer or Jacobi) orthogonal polynomial of variable $z$ and $λ$ is a parameter, are evaluated. A method we adopt in the present paper differs from that used by Fröhlich [Integral Transforms Spec. Funct. 2 (1994) 253] for the Jacobi polynomials and by Koepf [Integral Transforms Spec. Funct. 5 (1997) 69] for the generalized Laguerre and the Gegenbauer polynomials.

研究动机与目标

• 通过一种新方法重新推导经典正交多项式参数导数的已知展开式。
• 为广义拉盖尔多项式提供一种简单形式 ∂Lₙ^(λ)(z)/∂λ = Σ_{k=0}^{n-1} 1/(n−k) Lₖ^(λ)(z) 的替代推导。
• 通过其微分方程的参数微分，建立适用于勒让德、盖根堡和拉盖尔多项式的系统化方法。
• 提出一种独立于弗罗利希的超几何函数方法和科夫的极限过程的方法，基于幂级数展开中系数的比较。

提出的方法

• 通过对勒让德多项式的二阶微分方程关于参数 α 求导，推导其参数导数。
• 利用在 z′ = 1 处求值的克里斯托夫-达布瓦公式，将导数表示为低次勒让德多项式的线性组合。
• 应用已知的微分恒等式 (2n+α+β)(z²−1)Pₙ′(z) = n[β−α+(2n+α+β)z]Pₙ(z) − 2(n+α)(n+β)Pₙ₋₁(z)，以简化导数表达式。
• 通过比较 (z−z₀)ⁿ 展开中幂级数系数，建立系数 aₙₖ(λ)，并利用一个将 aₙₙ(λ) 与首项系数 kₙ(λ) 的对数导数联系起来的引理。
• 通过利用其与勒让德多项式的关系及极限形式，将相同方法应用于盖根堡和广义拉盖尔多项式。
• 利用恒等式 Lₙ^(α)(z) = lim_{β→∞} Pₙ^(α,β)(1−2z/β) 推导拉盖尔情形下的简单形式。

实验结果

研究问题

• RQ1 能否开发一种新方法，无需依赖超几何函数或极限过程，来计算经典正交多项式参数导数？
• RQ2 参数导数 ∂Pₙ(λ;z)/∂λ 在相同正交多项式基下的展开系数 aₙₖ(λ) 的显式形式是什么？
• RQ3 克里斯托夫-达布瓦恒等式如何用于将勒让德多项式的导数表示为低次多项式的线性组合？
• RQ4 为何广义拉盖尔多项式导数简化为 ∂Lₙ^(λ)(z)/∂λ = Σ_{k=0}^{n-1} 1/(n−k) Lₖ^(λ)(z)，且能否直接从微分方程推导出此结果？
• RQ5 是否存在一种统一方法，适用于勒让德、盖根堡和拉盖尔多项式的参数导数，且避免依赖于特殊函数的方法？

主要发现

• 在展开式 ∂Pₙ(λ;z)/∂λ = Σₖ aₙₖ(λ)Pₖ(λ;z) 中，系数 aₙₙ(λ) 由首项系数的对数导数给出：aₙₙ(λ) = ∂/∂λ ln kₙ(λ)。
• 对于勒让德多项式，关于 α 的导数为：∂Pₙ^(α,β)(z)/∂α = [ψ(2n+α+β+1)−ψ(n+α+β+1)]Pₙ^(α,β)(z) + Σ_{k=0}^{n−1} [ (2k+α+β+1)Γ(n+β+1)Γ(k+α+β+1) / ((n−k)(k+n+α+β+1)Γ(k+β+1)Γ(n+α+β+1)) ] Pₖ^(α,β)(z)。
• 盖根堡多项式 Cₙ^(λ)(z) 关于 λ 的导数为：∂Cₙ^(λ)(z)/∂λ = [ψ(n+λ)−ψ(λ)]Cₙ^(λ)(z) + 4Σ_{k=0}^{n−1} [1+(−1)^{k+n}]/2 × [ (k+λ) / ((n−k)(k+n+2λ)) ] Cₖ^(λ)(z)。
• 广义拉盖尔多项式导数简化为 ∂Lₙ^(λ)(z)/∂λ = Σ_{k=0}^{n−1} 1/(n−k) Lₖ^(λ)(z)，该结果通过与勒让德多项式的极限关系及新方法推导得出。
• 该方法在不使用超几何函数或组合恒等式的情况下，对已知结果提供了自洽的推导，与弗罗利希和科夫的方法有本质不同。
• 克里斯托夫-达布瓦公式在 z′=1 处求值时，可将导数表达式转化为低次正交多项式的线性组合。