[论文解读] 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 阶导数的闭式表达式。该方法利用倒伽马函数的导数、Pochhammer 符号及其广义形式，结合斯特林数与广义贝尔努利多项式，实现了对 $ \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) $ 的 k 阶导数的闭式表达式，适用于实数与复数 $ \nu $。
• 提供 $ \Gamma^{-1}(\nu+1+m) $、$ (\nu+1+m)_m $ 与 $ (\nu+1+m)^{-1}_m $ 的导数的显式公式，这些在 $ J_\nu(z) $ 的级数展开中至关重要。
• 实现高阶导数的直接计算，避免逐次计算低阶导数，从而提升数值效率。
• 将结果扩展至整数 $ \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} $，以处理伽马函数与 Pochhammer 符号的导数。
• 利用第一类斯特林数推导 $ \mathcal{P}_m^{(k)}(t) $ 的显式表达式：当 $ m \geq k $ 时，有 $ \mathcal{P}_m^{(k)}(t) = (-1)^{m-k} \sum_{l=0}^{m-k} (-1)^l \binom{m}{l} s(m-l,k) (t)_l $。
• 利用生成函数与级数展开，将 $ \mathcal{G}^{(k)}(1+\varepsilon) $ 表示为 $ \varepsilon $ 的幂级数，适用于 $ |\varepsilon| < 1 $。
• 应用 $ \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 能否推导出关于 $ \nu $ 的 $ J_\nu(z) $ 的 k 阶导数的闭式表达式，从而避免递归计算低阶导数？
• RQ2 如何将倒伽马函数与 Pochhammer 符号的导数用已知特殊函数（如斯特林数与广义贝尔努利多项式）表示？
• RQ3 当 $ \nu $ 为正整数或负整数时，$ \partial^k J_\nu(z)/\partial \nu^k $ 的精确表达式是什么？
• RQ4 该方法如何推广至复数 $ \nu $？其收敛性与数值稳定性存在哪些限制？

主要发现

• 关于 $ \nu $ 的 $ J_\nu(z) $ 的 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] $。
• 只要 $ |\Im \nu| \lesssim 1/2 $，该方法对复数 $ \nu $ 依然有效，但当 $ |t| $ 较大时，$ \mathcal{G}^{(k)}(t) $ 的级数收敛性会变差。
• 利用斯特林数推导出 $ \mathcal{P}_m^{(k)}(t) $ 的显式表达式：当 $ m \geq k $ 时，有 $ \mathcal{P}_m^{(k)}(t) = (-1)^{m-k} \sum_{l=0}^{m-k} (-1)^l \binom{m}{l} s(m-l,k) (t)_l $。
• 广义贝尔努利多项式表示 $ \mathcal{P}_m^{(k)}(t) = (-1)^{m-k} \binom{m}{k} B_{m-k}^{(m+1)}(1-t) $ 提供了另一种闭式表达形式。