[論文レビュー] Alternative Jacobi Polynomials and Orthogonal Exponentials

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

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

ひとこと要約

本稿では、ヤコビの『逆に、常に逆に』の原則にインspiredされた双方向直交化アルゴリズムを用いて、[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.

研究の動機と目的

• 近壁乱流のためのスペクトル法にインspiredされ、区間 [0,1] 上で古典的ヤコビ多項式の代替となる新たな直交多項式のクラスを開発すること。
• 半直線 [0,∞) 上で (α,β)-パrameter化された直交指数多項式系を構築し、指数的減衰を示す関数の近似に適すること。
• 直交指数多項式および単位関数から導かれる、[0,1] 上の2つのパrameter化された離散的ほぼ直交関数系を構築し、hp-近似および数値積分に用いること。
• パrameterシフトの下で、代替ヤコビ多項式と古典的ヤコビ多項式との間のアルゴリズム的同等性を確立し、数値解析における有用性を示すこと。
• 特異直交性の例外的ケースを探索し、カウチの公式による積分表現を提供することで、理論的および計算的適用性を高めること。

提案手法

• ロドリゲス型公式を用いて代替ヤコビ多項式を定義する：$\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)$。
• 部分積分を用いて、重み $x^\alpha(1-x)^\beta$ に関して [0,1] 上での直交性を証明し、$\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)$ を導入し、ヤコビ系とアルゴリズム的に同等になるように保証する。
• スケーリングされた指数関数および単位関数を用いて、2つの離散的ほぼ直交関数系 $\boldsymbol{\widetilde{\mathcal{Z}}}_n^{(\omega_n)}(t)$ および $\boldsymbol{\mathcal{Z}}_n^{(\omega_n)}(t)$ を構築し、$\gamma_n$ が直交性の質を制御する。

実験結果

リサーチクエスチョン

• RQ1 双方向直交化アルゴリズムは、[0,1] 上で古典的ヤコビ多項式と同等の性質を持つ代替直交多項式をどのように生成できるか？
• RQ2 重み $x^\alpha(1-x)^\beta$ に関して、代替ヤコビ多項式 $\mathcal{P}_{nk}^{(\alpha,\beta)}(x)$ の直交構造とノルムは何か？
• RQ3 [0,∞) 上で、ヤコビ多項式系 [-1,1] におけるアルゴリズム的効率と同等となる直交指数多項式を体系的に構築できるか？
• RQ4 パラメータのシフトとスケーリングは、これらの代替系の直交性および再帰構造にどのように影響するか？
• RQ5 パラメータ $\gamma_n$ は、有限区間における hp-近似に用いる離散的ほぼ直交関数の構築において、果たす役割は何か？

主な発見

• 代替ヤコビ多項式 $\mathcal{P}_{nk}^{(\alpha,\beta)}(x)$ は、重み $x^\alpha(1-x)^\beta$ に関して [0,1] 上で直交し、$\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)$ であり、異なるパラメータ領域を結びつける。
• 直交指数多項式 $\mathcal{E}_{nk}^{(\alpha,\beta)}(t)$ は、双方向直交化を用いて構築され、ヤコビ系とアルゴリズム的に同等である。
• 2つの離散的ほぼ直交関数系 $\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$ が導出された。