2025-11-19T12:16:13.706625

Uniform Asymptotic approximation and numerical evaluation of the Reverse Generalized Bessel Polynomial zeros

Dunster, Gil, Ruiz-Antolin et al.
Uniform asymptotic expansions are derived for the zeros of the reverse generalized Bessel polynomials of large degree $n$ and real parameter $a$. It is assumed that $-Δ_{1} n+\frac{3}{2} \leq a \leq Δ_{2} n$ for fixed arbitrary $Δ_{1} \in (0,1)$ and bounded positive $Δ_{2}$. For this parameter range at most one of the zeros is real, with the rest being complex conjugates. The new expansions are uniformly valid for all the zeros, and are shown to be highly accurate for moderate or large values of $n$. They are consequently used as initial values in a very efficient numerical algorithm designed to obtain the remaining complex zeros using Taylor series.
academic

Uniform Asymptotic approximation and numerical evaluation of the Reverse Generalized Bessel Polynomial zeros

基本信息

  • 论文ID: 2510.14295
  • 标题: Uniform Asymptotic approximation and numerical evaluation of the Reverse Generalized Bessel Polynomial zeros
  • 作者: T. M. Dunster, A. Gil, D. Ruiz-Antolin, J. Segura
  • 分类: math.CA (Classical Analysis and ODEs)
  • 发表时间: 2025年10月16日
  • 论文链接: https://arxiv.org/abs/2510.14295

摘要

本文推导了大度数nn和实参数aa的逆广义贝塞尔多项式零点的一致渐近展开。假设参数范围为−Δ1n+32≤a≤Δ2n-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n，其中Δ1∈(0,1)\Delta_1 \in (0,1)为固定任意常数，Δ2\Delta_2为有界正常数。在此参数范围内，至多有一个零点为实数，其余均为复共轭对。新的展开对所有零点一致有效，在中等或大的nn值下表现出高精度。这些展开被用作高效数值算法的初值，该算法利用泰勒级数计算剩余的复零点。

研究背景与动机

问题描述

广义贝塞尔多项式定义为： yn(z;a)=∑k=0n(nk)(n+a−1)k(12z)ky_n(z; a) = \sum_{k=0}^{n} \binom{n}{k} (n+a-1)_k \left(\frac{1}{2z}\right)^k

其中(α)k=Γ(α+k)/Γ(α)(α)_k = Γ(α + k)/Γ(α)是Pochhammer符号。这些多项式的零点通常为复值，在应用数学和工程学中有多种应用。

研究动机

  1. 实际应用需求：广义贝塞尔多项式的零点在应用数学和工程领域有重要应用
  2. 计算挑战：现有计算这些零点的方法在效率和精度上存在局限
  3. 理论完善：需要建立一致有效的渐近理论来描述所有零点的行为
  4. 数值算法优化：需要高效的数值方法来计算复零点

现有方法局限性

  • 以往工作2, 13在计算精度和效率上有限制
  • 缺乏对所有零点一致有效的渐近展开
  • 数值算法的初值选择不够优化

核心贡献

  1. 推导了一致渐近展开：为逆广义贝塞尔多项式的零点建立了对所有零点一致有效的渐近展开式
  2. 建立了转折点理论：基于WKB理论和Liouville-Green方法，处理了具有转折点的微分方程
  3. 构造了高效数值算法：结合渐近展开和泰勒级数，开发了计算复零点的高效算法
  4. 提供了高精度近似：展开式在中等或大的nn值下达到10−1510^{-15}量级的相对误差

方法详解

任务定义

寻找逆广义贝塞尔多项式θn(z;a)=znyn(z−1;a)θ_n(z; a) = z^n y_n(z^{-1}; a)的零点tm(u,a)t_m(u, a)，其中u=n+12u = n + \frac{1}{2}，在参数范围−Δ1n+32≤a≤Δ2n-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n下。

核心理论框架

1. 微分方程变换

通过缩放函数： wn(0)(z;a)=2−n−a+1z1−n−a/2e−zθn(z;a)w_n^{(0)}(z; a) = 2^{-n-a+1} z^{1-n-a/2} e^{-z} θ_n(z; a)

将问题转化为求解微分方程： d2wdz2={1+a−2z+(2n+a)(2n+a−2)4z2}w\frac{d^2w}{dz^2} = \left\{1 + \frac{a-2}{z} + \frac{(2n+a)(2n+a-2)}{4z^2}\right\} w

2. Liouville-Green理论

定义参数： u=n+12,α=a−2uu = n + \frac{1}{2}, \quad α = \frac{a-2}{u}

微分方程重写为： d2wdz2={u2f(α,z)+g(z)}w\frac{d^2w}{dz^2} = \{u^2 f(α, z) + g(z)\} w

其中： f(α,z)=(z+12α)2+1+αz2,g(z)=−14z2f(α, z) = \frac{(z + \frac{1}{2}α)^2 + 1 + α}{z^2}, \quad g(z) = -\frac{1}{4z^2}

3. 转折点分析

转折点位于： z1,2(α)=±iσ−12α,σ=1+αz_{1,2}(α) = ±iσ - \frac{1}{2}α, \quad σ = \sqrt{1 + α}

通过Liouville变量ζζ和LG变量ξξ： 23ζ3/2=ξ=∫z1(α)zf1/2(α,t)dt\frac{2}{3}ζ^{3/2} = ξ = \int_{z_1(α)}^z f^{1/2}(α, t) dt

零点的渐近展开

零点满足隐式方程： Z(u,a,u−1tm(u,a))=u−2/3amZ(u, a, u^{-1}t_m(u, a)) = u^{-2/3}a_m

其中ama_m是Airy函数Ai(x)Ai(x)的第mm个负零点。

展开形式为： tm(u,a)∼u∑s=0∞τm,s(α)u2s(u→∞)t_m(u, a) \sim u \sum_{s=0}^{\infty} \frac{τ_{m,s}(α)}{u^{2s}} \quad (u → ∞)

主要系数计算

  • 领导项τm,0τ_{m,0}：由复杂的隐式方程(3.16)确定
  • 后续项：通过递推关系(3.20)-(3.23)计算：

τm,1=−Υm,1ζm,0′τ_{m,1} = -\frac{Υ_{m,1}}{ζ'_{m,0}}

τm,2=−12ζm,0′{τm,12ζm,0′′+2τm,1Υm,1′+2Υm,2}τ_{m,2} = -\frac{1}{2ζ'_{m,0}}\{τ_{m,1}^2 ζ''_{m,0} + 2τ_{m,1} Υ'_{m,1} + 2Υ_{m,2}\}

数值算法设计

1. 迭代函数

Tn(a,z)=z−1Ωn(a,z)arctan⁡(Ωn(a,z)wn(0)(z;a)∂wn(0)(z;a)/∂z)T_n(a, z) = z - \frac{1}{\sqrt{Ω_n(a, z)}} \arctan\left(\frac{\sqrt{Ω_n(a, z)}w_n^{(0)}(z; a)}{\partial w_n^{(0)}(z; a)/\partial z}\right)

2. 步长函数

Hn(a,z)=z+πΩn(a,z)H_n(a, z) = z + \frac{π}{\sqrt{Ω_n(a, z)}}

3. 泰勒级数展开

利用递推关系(4.7)计算高阶导数，避免直接计算函数值。

实验设置

测试参数

  • 度数范围：n=15,30,50,200,500,1000,2000n = 15, 30, 50, 200, 500, 1000, 2000
  • 参数值：a=1.01,20.2a = 1.01, 20.2等
  • 零点索引：m=1,3,10,15,25m = 1, 3, 10, 15, 25等

评价指标

  • 相对误差：与Maple高精度实现的比较
  • 计算时间：CPU执行时间
  • 收敛性：迭代算法的收敛精度设为ε=10−12ε = 10^{-12}

实现细节

  • 渐近展开使用前5项系数
  • Maple实现使用60位精度
  • Matlab算法在Dell Latitude 7410上测试

实验结果

主要结果

1. 高精度近似

表1和表2显示的相对误差：

  • 对于a=1.01a = 1.01：相对误差在10−1810^{-18}到10−1510^{-15}量级
  • 对于a=20.2a = 20.2：相对误差在10−1810^{-18}到10−1510^{-15}量级
  • 所有测试案例的相对误差都接近或优于10−1510^{-15}

2. 参数依赖性分析

图1显示随着参数aa增加，相对误差递减，最大相对误差小于4×10−154×10^{-15}。

3. 算法效率

表3的计算时间结果：

  • n=30n = 30（15个零点）：2.8×10−32.8×10^{-3}秒
  • n=2000n = 2000（1000个零点）：3.7×10−23.7×10^{-2}秒
  • 零点数量增加67倍，计算时间仅增加13倍，显示算法的高效性

收敛性分析

图4显示即使计算大量零点，相对误差仍保持良好控制，算法具有优秀的数值稳定性。

案例验证

图2展示了F(w)=0F(w) = 0方程求解的可视化，Newton方法能准确定位零点。

相关工作

历史发展

  1. Carpenter (1992)2：早期的广义贝塞尔多项式零点渐近分析
  2. Pasquini (2000)13：精确计算广义贝塞尔多项式零点的数值方法
  3. Segura (2013)14：特殊函数复零点的通用计算方法
  4. de Bruin等 (1981)3,4：广义贝塞尔多项式零点的分布域研究

本文优势

  • 一致性：展开对所有零点一致有效
  • 精度：达到机器精度水平
  • 效率：算法复杂度显著优于现有方法
  • 适用性：覆盖更广的参数范围

结论与讨论

主要结论

  1. 理论贡献：建立了逆广义贝塞尔多项式零点的完整渐近理论
  2. 数值成就：开发了高精度、高效率的零点计算算法
  3. 实用价值：为相关应用提供了可靠的计算工具

局限性

  1. 参数限制：方法适用于−Δ1n+32≤a≤Δ2n-\Delta_1 n + \frac{3}{2} \leq a \leq \Delta_2 n的参数范围
  2. 复杂性：渐近展开的推导过程数学上较为复杂
  3. 实现要求：需要高精度算术运算支持

未来方向

  1. 扩展到更一般的参数范围
  2. 研究其他特殊函数零点的类似方法
  3. 优化算法的计算效率
  4. 探索在实际工程问题中的应用

深度评价

优点

  1. 理论严谨：基于转折点理论和WKB方法的严格数学推导
  2. 方法创新：将渐近展开与数值算法巧妙结合
  3. 精度卓越：达到10−1510^{-15}量级的相对误差
  4. 效率突出：算法复杂度优于传统方法
  5. 实用性强：提供了完整的实现和开源代码

不足

  1. 理论复杂：方法的数学基础较为深奥，理解门槛较高
  2. 参数依赖：方法的有效性依赖于特定的参数范围假设
  3. 实现复杂：需要计算多个复杂的系数和导数

影响力

  1. 学术价值：为特殊函数理论提供了重要的理论工具
  2. 应用前景：在信号处理、控制理论等领域有潜在应用
  3. 方法论意义：为类似问题的研究提供了范例

适用场景

  1. 科学计算：需要高精度计算贝塞尔多项式零点的场合
  2. 工程应用：滤波器设计、控制系统分析等
  3. 理论研究：特殊函数和渐近分析的进一步研究

参考文献

本文引用了14篇重要参考文献，涵盖了渐近分析、特殊函数、数值算法等相关领域的经典和前沿工作，为研究提供了坚实的理论基础。

总评：这是一篇高质量的数学分析论文，在理论推导和数值实现方面都达到了很高的水准。方法的创新性、结果的精确性以及算法的高效性都值得称赞，为相关领域的研究和应用提供了有价值的工具。