[論文レビュー] Derivatives with respect to the order of the Bessel function of the first kind

J. Sesma|arXiv (Cornell University)|Jan 20, 2014

Differential Equations and Boundary Problems参考文献 1被引用数 8

ひとこと要約

本稿では、第一種ベッセル関数 $ J_\nu(z) $ の順序 $ \nu $ に関する $ k $ 階微分について、閉形式の式を導出する。この導出には、逆ガンマ関数の微分、ポッホハマー記号、およびそれらの一般化形態であるスターリング数と一般化ベルヌーイ多項式が用いられ、$ \nu $ が正または負の整数である場合の正確な計算が可能になる。主な貢献は、低階微分の逐次評価を回避する体系的な手法を提供することであり、高階微分の計算を直接的に行えるようにする。

ABSTRACT

An explicit expression of the k-th derivative of the Bessel function $J_ν(z)$, with respect to its order $ν$, is given. Particularizations for the cases of positive or negative $ν$ are considered.

研究の動機と目的

• 実数および複素数 $ \nu $ に対して有効な、$ J_\nu(z) $ の順序 $ \nu $ に関する $ k $ 階微分の閉形式表現を導出すること。
• 系列展開 $ J_\nu(z) $ に不可欠な $ \Gamma^{-1}(\nu+1+m) $、$ (\nu+1+m)_m $、および $ (\nu+1+m)^{-1}_m $ の微分に対する明示的公式を提供すること。
• 低階微分の逐次評価を回避し、直接的な高階微分計算を可能にすることで、数値的効率を向上させること。
• 特殊関数の恒等式と級数展開を用いて、$ \nu $ が正または負の整数である場合への結果の拡張を図ること。

提案手法

• $ J_\nu(z) = (z/2)^\nu \sum_{m=0}^\infty \frac{(-z^2/4)^m}{m! \Gamma(\nu+1+m)} $ の級数定義を、$ \nu $ に関する微分の出発点として用いる。
• ガンマ関数およびポッホハマー関数の微分を扱うために、補助関数 $ \mathcal{G}^{(k)}(t) = \frac{d^k}{dt^k} \Gamma^{-1}(t) $、$ \mathcal{P}_m^{(k)}(t) = \frac{1}{k!} \frac{d^k}{dt^k} (t)_m $、および $ \mathcal{Q}_m^{(k)}(t) = \frac{1}{k!} \frac{d^k}{dt^k} (t)_m^{-1} $ を導入する。
• スターリング数を用いて、$ \mathcal{P}_m^{(k)}(t) $ の明示的表現を導出する：$ \mathcal{P}_m^{(k)}(t) = (-1)^{m-k} \sum_{l=0}^{m-k} (-1)^l \binom{m}{l} s(m-l,k) (t)_l $（$ m \geq k $ の場合）。
• 母関数と級数展開を用いて、$ |\varepsilon| < 1 $ の範囲で $ \mathcal{G}^{(k)}(1+\varepsilon) $ を $ \varepsilon $ のべき級数として表現する。
• 数値的安定性を高めるために、$ \mathcal{Q}_m^{(k)}(t) $ の再帰関係 $ \mathcal{Q}_{m+1}^{(k)}(t) = \left( \mathcal{Q}_m^{(k)}(t) - \mathcal{Q}_{m+1}^{(k-1)}(t) \right)/(t+m) $ を適用する。
• $ \partial^k J_\nu(z)/\partial \nu^k $ の最終的表現を、級数再展開により導出し、$ \nu = N + \varepsilon $（$ N \in \mathbb{Z} $、$ |\Re \varepsilon| \leq 1/2 $）のケースに応じて場合分けを行う。

実験結果

リサーチクエスチョン

• RQ1 低階微分の逐次評価を回避して、$ J_\nu(z) $ の順序 $ \nu $ に関する $ k $ 階微分の閉形式表現を導出可能か？
• RQ2 逆ガンマ関数およびポッホハマー記号の微分は、スターリング数や一般化ベルヌーイ多項式といった既知の特殊関数を用いてどのように表現できるか？
• RQ3 $ \nu $ が正または負の整数である場合の $ \partial^k J_\nu(z)/\partial \nu^k $ の正確な表現は何か？
• RQ4 本手法は複素数 $ \nu $ に対しても適用可能か？収束性および数値的安定性の観点から、その限界は何か？

主な発見

• $ \partial^k J_\nu(z)/\partial \nu^k $ の $ k $ 階微分は、$ \mathcal{G}^{(k_2)}(1+\varepsilon) $、$ \mathcal{P}_{-N-m}^{(k-k_1-k_2)}(-\varepsilon) $、および $ \mathcal{Q}_{m+N}^{(k-k_1-k_2)}(1+\varepsilon) $ を含む有限和として表現され、直接計算が可能になる。
• $ \nu = n \in \mathbb{Z}_{\geq 0} $ の場合、最初の微分は $ \left. \frac{\partial}{\partial \nu} J_\nu(z) \right|_{\nu=n} = \left( \ln(z/2) + \gamma \right) J_n(z) - (z/2)^n \sum_{m=0}^\infty \frac{(-z^2/4)^m}{m! (m+n)!} \hat{H}_{m+n}^{(1)} $ となる。ここで $ \gamma $ はオイラー・マスケローニ定数である。
• $ \nu = -n $（$ n > 0 $）の場合、最初の微分は $ \left. \frac{\partial}{\partial \nu} J_\nu(z) \right|_{\nu=-n} = \left( \ln(z/2) + \gamma \right) J_{-n}(z) - (z/2)^{-n} \left[ (-1)^n \sum_{m=0}^{n-1} \frac{(z^2/4)^m}{m! (n-m-1)!} + \sum_{m=n}^\infty \frac{(-z^2/4)^m}{m! (m-n)!} \hat{H}_{m-n}^{(1)} \right] $ となる。
• 複素数 $ \nu $ に対しても、$ |\Im \nu| \lesssim 1/2 $ の範囲では本手法が有効であるが、$ |t| $ が大きくなると $ \mathcal{G}^{(k)}(t) $ の級数の収束性が劣化する。
• $ \mathcal{P}_m^{(k)}(t) $ の明示的表現はスターリング数を用いて導出され、$ \mathcal{P}_m^{(k)}(t) = (-1)^{m-k} \sum_{l=0}^{m-k} (-1)^l \binom{m}{l} s(m-l,k) (t)_l $（$ m \geq k $）となる。
• 一般化ベルヌーイ多項式による表現 $ \mathcal{P}_m^{(k)}(t) = (-1)^{m-k} \binom{m}{k} B_{m-k}^{(m+1)}(1-t) $ が、代替的な閉形式表現を提供する。