We derive two distinct asymptotic expansions for the zeros $j_{ν,k}^{(n)}$ of the $n$-th derivative of Bessel function $J_ν^{(n)}(x)$. The first is a McMahon-type expansion for the case when $k \to \infty$ with fixed $ν$, for which we also establish an explicit error bound. The second addresses the case when $ν\to \infty$ with fixed $k$ and it involves the zeros of Airy functions and their derivatives. These results extend and refine the classical work of Wong, Lang, and Olver on the zeros of Bessel functions. In the course of obtaining our main results, we also generalize several auxiliary results, which in turn provide a broader framework for the study of zeros of special functions.
论文ID : 2510.12353标题 : Asymptotic behavior of zeros of Bessel function derivatives作者 : Árpád Baricz, Pranav Kumar, Saminathan Ponnusamy分类 : math.CA (Classical Analysis and ODEs)发表时间 : 2025年10月14日论文链接 : https://arxiv.org/abs/2510.12353v1 本文推导了Bessel函数的n阶导数 J ν ( n ) ( x ) J_\nu^{(n)}(x) J ν ( n ) ( x ) j ν , k ( n ) j_{\nu,k}^{(n)} j ν , k ( n ) k → ∞ k \to \infty k → ∞ ν \nu ν ν → ∞ \nu \to \infty ν → ∞ k k k
Bessel函数零点在应用数学和数学物理的各个问题中发挥着重要作用，特别是在量子力学、散射理论、波传播等相关领域。由于其广泛的适用性，理解Bessel函数及其零点的渐近行为受到了相当的关注。
经典McMahon展开的限制 : 早期McMahon对Bessel函数 J ν ( x ) J_\nu(x) J ν ( x ) Wong和Lang的工作 : 虽然他们扩展了这些思想来获得 J ν ′ ′ ( x ) J''_\nu(x) J ν ′′ ( x ) 研究空白 : 尽管在理解n阶Bessel函数导数零点性质方面取得了重大进展，但这些零点的渐近展开在文献中尚未被探索。本文旨在填补这一空白，通过呈现McMahon型渐近展开并分析 ν → ∞ \nu \to \infty ν → ∞ J ν ( n ) ( x ) J_\nu^{(n)}(x) J ν ( n ) ( x ) j ν , k ( n ) j_{\nu,k}^{(n)} j ν , k ( n ) k → ∞ k \to \infty k → ∞
建立了McMahon型渐近展开 : 对于Bessel函数n阶导数的零点，当 k → ∞ k \to \infty k → ∞ ν \nu ν 提供了显式误差界 : 为McMahon型展开建立了明确的误差估计，量化了近似的精度。推导了大参数渐近展开 : 处理 ν → ∞ \nu \to \infty ν → ∞ k k k 推广了辅助结果 : 在获得主要结果的过程中，推广了若干基础理论结果，为特殊函数零点研究提供了更广泛的框架。扩展了经典理论 : 将Wong、Lang和Olver的经典工作从Bessel函数本身扩展到其任意阶导数。研究Bessel函数第一类 J ν ( x ) J_\nu(x) J ν ( x ) J ν ( n ) ( x ) J_\nu^{(n)}(x) J ν ( n ) ( x ) j ν , k ( n ) j_{\nu,k}^{(n)} j ν , k ( n )
ν > 0 \nu > 0 ν > 0 n ≥ 1 n \geq 1 n ≥ 1 k k k 定理2.1 : 对于 n ∈ N 0 n \in \mathbb{N}_0 n ∈ N 0 x x x
1 2 π x J ν ( 2 n ) ( x ) = cos ( x − ν π 2 − π 4 ) τ ν ( 2 n ) ( x ) + sin ( x − ν π 2 − π 4 ) θ ν ( 2 n ) ( x ) \sqrt{\frac{1}{2}\pi x} J_\nu^{(2n)}(x) = \cos\left(x - \frac{\nu\pi}{2} - \frac{\pi}{4}\right)\tau_\nu^{(2n)}(x) + \sin\left(x - \frac{\nu\pi}{2} - \frac{\pi}{4}\right)\theta_\nu^{(2n)}(x) 2 1 π x J ν ( 2 n ) ( x ) = cos ( x − 2 ν π − 4 π ) τ ν ( 2 n ) ( x ) + sin ( x − 2 ν π − 4 π ) θ ν ( 2 n ) ( x )
其中：
τ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ α 2 m , 2 n x 2 m , θ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ α 2 m + 1 , 2 n x 2 m + 1 \tau_\nu^{(2n)}(x) = \sum_{m=0}^{\infty} \frac{\alpha_{2m,2n}}{x^{2m}}, \quad \theta_\nu^{(2n)}(x) = \sum_{m=0}^{\infty} \frac{\alpha_{2m+1,2n}}{x^{2m+1}} τ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ x 2 m α 2 m , 2 n , θ ν ( 2 n ) ( x ) = ∑ m = 0 ∞ x 2 m + 1 α 2 m + 1 , 2 n
系数满足递推关系：
α 2 m , 2 n = α 2 m , 2 n − 1 − 4 m − 1 2 α 2 m − 1 , 2 n − 1 \alpha_{2m,2n} = \alpha_{2m,2n-1} - \frac{4m-1}{2}\alpha_{2m-1,2n-1} α 2 m , 2 n = α 2 m , 2 n − 1 − 2 4 m − 1 α 2 m − 1 , 2 n − 1 α 2 m + 1 , 2 n = − α 2 m − 1 , 2 n − 1 − 4 m + 1 2 α 2 m , 2 n − 1 \alpha_{2m+1,2n} = -\alpha_{2m-1,2n-1} - \frac{4m+1}{2}\alpha_{2m,2n-1} α 2 m + 1 , 2 n = − α 2 m − 1 , 2 n − 1 − 2 4 m + 1 α 2 m , 2 n − 1
通过反演技术，得到零点的McMahon型展开：
偶数阶导数 (J ν ( 2 m ) J_\nu^{(2m)} J ν ( 2 m ) j ν , k + 1 ( 2 m ) = α + α 1 , 2 m α 0 , 2 m ⋅ 1 α + O ( 1 α 3 ) j_{\nu,k+1}^{(2m)} = \alpha + \frac{\alpha_{1,2m}}{\alpha_{0,2m}} \cdot \frac{1}{\alpha} + O\left(\frac{1}{\alpha^3}\right) j ν , k + 1 ( 2 m ) = α + α 0 , 2 m α 1 , 2 m ⋅ α 1 + O ( α 3 1 )
其中 α = k π + ν π 2 − π 4 \alpha = k\pi + \frac{\nu\pi}{2} - \frac{\pi}{4} α = kπ + 2 ν π − 4 π
奇数阶导数 (J ν ( 2 m + 1 ) J_\nu^{(2m+1)} J ν ( 2 m + 1 ) j ν , k ( 2 m + 1 ) = β − α 1 , 2 m + 1 α 0 , 2 m + 1 ⋅ 1 β + O ( 1 β 3 ) j_{\nu,k}^{(2m+1)} = \beta - \frac{\alpha_{1,2m+1}}{\alpha_{0,2m+1}} \cdot \frac{1}{\beta} + O\left(\frac{1}{\beta^3}\right) j ν , k ( 2 m + 1 ) = β − α 0 , 2 m + 1 α 1 , 2 m + 1 ⋅ β 1 + O ( β 3 1 )
其中 β = k π + ν π 2 − 3 π 4 \beta = k\pi + \frac{\nu\pi}{2} - \frac{3\pi}{4} β = kπ + 2 ν π − 4 3 π
定理2.2 : 对于渐近展开的误差项，有界估计：
对于2n阶导数，当 ν ≥ − 2 n + 1 2 \nu \geq -2n + \frac{1}{2} ν ≥ − 2 n + 2 1 n ∈ N n \in \mathbb{N} n ∈ N x x x ∣ δ 2 n ( ν , x ) ∣ ≤ 4 ( ν + 2 n ) 2 − 1 4 x exp { 4 ( ν + 2 n ) 2 − 1 4 x } |\delta_{2n}(\nu,x)| \leq \frac{4(\nu+2n)^2-1}{4x} \exp\left\{\frac{4(\nu+2n)^2-1}{4x}\right\} ∣ δ 2 n ( ν , x ) ∣ ≤ 4 x 4 ( ν + 2 n ) 2 − 1 exp { 4 x 4 ( ν + 2 n ) 2 − 1 }
当 ν → ∞ \nu \to \infty ν → ∞ k k k
对于偶数阶导数：
j ν , k ( 2 n ) = ν − a k 2 1 / 3 ν 1 / 3 + 2 1 / 3 a k 2 10 ν − 1 / 3 + O ( ν − 1 ) j_{\nu,k}^{(2n)} = \nu - \frac{a_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}a_k^2}{10}\nu^{-1/3} + O(\nu^{-1}) j ν , k ( 2 n ) = ν − 2 1/3 a k ν 1/3 + 10 2 1/3 a k 2 ν − 1/3 + O ( ν − 1 )
对于奇数阶导数：
j ν , k ( 2 n + 1 ) = ν − a k ′ 2 1 / 3 ν 1 / 3 + 2 1 / 3 ( a k ′ ) 2 10 ν − 1 / 3 + O ( ν − 1 ) j_{\nu,k}^{(2n+1)} = \nu - \frac{a'_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}(a'_k)^2}{10}\nu^{-1/3} + O(\nu^{-1}) j ν , k ( 2 n + 1 ) = ν − 2 1/3 a k ′ ν 1/3 + 10 2 1/3 ( a k ′ ) 2 ν − 1/3 + O ( ν − 1 )
其中 a k a_k a k a k ′ a'_k a k ′ Ai ( x ) \text{Ai}(x) Ai ( x ) Ai ′ ( x ) \text{Ai}'(x) Ai ′ ( x )
统一的递推框架 : 建立了所有阶导数系数的统一递推关系，避免逐个计算每阶导数。改进的误差分析 : 通过修改Hethcote的方法，得到了比经典结果更精确的误差界。扩展的适用范围 : 将经典理论从特定的低阶导数扩展到任意阶导数。Uniform渐近理论的应用 : 巧妙地将Olver的uniform渐近展开理论应用到高阶导数情况。由于这是纯数学理论研究，"实验"主要指理论验证和数值检验：
数学归纳法验证 : 使用数学归纳法证明递推关系的正确性渐近分析 : 通过极限分析验证渐近展开的有效性误差界检验 : 通过具体计算验证误差界的紧致性McMahon的原始结果（针对 J ν ( x ) J_\nu(x) J ν ( x ) Wong和Lang的结果（针对 J ν ′ ′ ( x ) J''_\nu(x) J ν ′′ ( x ) Olver的uniform渐近展开理论 McMahon型展开的建立 : 成功推导出任意阶Bessel函数导数零点的McMahon型渐近展开，具有明确的系数表达式。误差界的改进 : 相比于已有结果，新的误差界更加精确，特别是对于高阶导数的情况。大参数展开的完整性 : 对于 ν → ∞ \nu \to \infty ν → ∞ 对于三阶和四阶导数的零点，得到：
j ν , k ( 3 ) = ν − a k 2 1 / 3 ν 1 / 3 + 2 1 / 3 a k 2 10 ν − 1 / 3 + a k 3 700 ν − 1 + ⋯ j_{\nu,k}^{(3)} = \nu - \frac{a_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}a_k^2}{10}\nu^{-1/3} + \frac{a_k^3}{700}\nu^{-1} + \cdots j ν , k ( 3 ) = ν − 2 1/3 a k ν 1/3 + 10 2 1/3 a k 2 ν − 1/3 + 700 a k 3 ν − 1 + ⋯
j ν , k ( 4 ) = ν − a k ′ 2 1 / 3 ν 1 / 3 + 2 1 / 3 ( a k ′ ) 2 10 ν − 1 / 3 + ( a k ′ ) 3 700 ν − 1 + ⋯ j_{\nu,k}^{(4)} = \nu - \frac{a'_k}{2^{1/3}}\nu^{1/3} + \frac{2^{1/3}(a'_k)^2}{10}\nu^{-1/3} + \frac{(a'_k)^3}{700}\nu^{-1} + \cdots j ν , k ( 4 ) = ν − 2 1/3 a k ′ ν 1/3 + 10 2 1/3 ( a k ′ ) 2 ν − 1/3 + 700 ( a k ′ ) 3 ν − 1 + ⋯
对于McMahon型展开，当满足条件时：
偶数阶导数：∣ j ν , k + 1 ( 2 n ) − ( k + ν 2 − 1 4 ) π ∣ ≤ 0.26 ( 4 ( ν + 2 n ) 2 − 1 ) π ( k + ν 2 − 0.314 ) |j_{\nu,k+1}^{(2n)} - (k + \frac{\nu}{2} - \frac{1}{4})\pi| \leq \frac{0.26(4(\nu+2n)^2-1)}{\pi(k + \frac{\nu}{2} - 0.314)} ∣ j ν , k + 1 ( 2 n ) − ( k + 2 ν − 4 1 ) π ∣ ≤ π ( k + 2 ν − 0.314 ) 0.26 ( 4 ( ν + 2 n ) 2 − 1 ) 奇数阶导数：∣ j ν , k ( 2 n + 1 ) − ( k + ν 2 − 1 4 ) π ∣ ≤ 0.26 ( 4 ( ν + 2 n + 1 ) 2 − 1 ) π ( k + ν 2 − 0.314 ) |j_{\nu,k}^{(2n+1)} - (k + \frac{\nu}{2} - \frac{1}{4})\pi| \leq \frac{0.26(4(\nu+2n+1)^2-1)}{\pi(k + \frac{\nu}{2} - 0.314)} ∣ j ν , k ( 2 n + 1 ) − ( k + 2 ν − 4 1 ) π ∣ ≤ π ( k + 2 ν − 0.314 ) 0.26 ( 4 ( ν + 2 n + 1 ) 2 − 1 ) McMahon (1895) : 建立了Bessel函数零点的经典渐近展开Watson (1944) : 完善了理论基础并给出了更严格的证明Olver (1954) : 发展了uniform渐近展开理论Wong & Lang (1990-1991) : 将理论扩展到Bessel函数的二阶导数Baricz等 (2018) : 证明了n阶导数零点的实性和正性本文在这一发展脉络中的贡献是将理论从特定的低阶导数情况完全推广到任意阶导数，并提供了更精确的误差分析。这是该领域的一个重要理论完善。
完整性 : 建立了Bessel函数任意阶导数零点的完整渐近理论精确性 : 提供了比已有结果更精确的误差界统一性 : 给出了统一的理论框架，适用于所有阶的导数填补空白 : 完成了Bessel函数导数零点渐近理论的最后一块拼图方法创新 : 提供了处理高阶导数渐近问题的一般性方法应用潜力 : 为相关的物理和工程问题提供了理论工具适用条件 : 某些结果需要满足特定的参数范围条件计算复杂性 : 高阶项的计算可能变得相当复杂数值验证 : 理论结果需要更多的数值验证来确认其实用性推广到修正Bessel函数 : 将类似方法应用到修正Bessel函数及其导数其他特殊函数 : 扩展到Wright函数、超Bessel函数等数值算法 : 基于理论结果开发高效的零点计算算法理论完整性 : 提供了该问题的完整理论解决方案方法严谨性 : 使用了成熟的数学分析方法，证明过程严谨结果的一般性 : 不仅解决了特定问题，还提供了一般性的理论框架误差分析 : 提供了明确的误差界，增强了结果的实用性缺少数值验证 : 作为纯理论研究，缺少充分的数值实验验证应用示例不足 : 没有提供具体的应用场景和实例计算复杂度 : 对于非常高阶的导数，计算可能变得不切实际学术影响 : 完善了特殊函数理论的重要分支实用价值 : 为需要Bessel函数导数零点的应用提供了理论基础方法论价值 : 提供的技术方法可能适用于其他类似问题数学物理 : 波动方程、热传导等问题的解析解工程应用 : 信号处理、振动分析等领域数值计算 : 高精度零点计算算法的理论基础论文引用了该领域的关键文献，包括：
McMahon的开创性工作 (1895) Watson的经典教科书 (1944) Olver的渐近理论 (1954, 1962, 1974) Wong & Lang的扩展工作 (1990, 1991) 近期相关研究 (Baricz等, 2018; Dimitrov & Lun, 2025) 这些引用体现了研究的历史连续性和理论完整性。
