[论文解读] Non-recursive expressions for even-index Bernoulli numbers: A remarkable sequence of determinants

Renaat Van Malderen|ArXiv.org|May 20, 2005

Neural Networks and Applications被引用 7

一句话总结

本文提出了偶数索引伯努利数 $ B_{2p} $ 的非递归显式公式，利用一组行列式 $ D_k $，这些行列式满足递推关系，并与生成函数 $ \frac{z}{\sin z} $ 相关联。主要贡献在于推导出 $ B_{2p} $ 的三种等价闭式表达式，其中行列式序列 $ D_p $ 的计算结果为 $ D_p = \frac{2}{\pi^{2p}} \left(1 - \frac{1}{2^{2p-1}} \right) \zeta(2p) $，从而将伯努利数与黎曼 zeta 函数联系起来。

ABSTRACT

Bernoulli numbers are usually expressed in terms of their lower index numbers (recursive). This paper gives explicit formulas for Bernoulli numbers of even index. The formulas contain a remarkable sequence of determinants. The value of these determinants for variable dimension is investigated.

研究动机与目标

• 推导偶数索引伯努利数 $ B_{2p} $ 的显式非递归表达式，避免递归计算。
• 研究支撑新公式的行列式序列 $ D_k $ 的结构及其计算方法。
• 建立行列式序列 $ D_k $、生成函数 $ \frac{z}{\sin z} $ 与黎曼 zeta 函数 $ \zeta(2p) $ 之间的联系。
• 利用变换技术与行列式恒等式，提供 $ B_{2p} $ 的多种等价闭式表达式。

提出的方法

• 定义一组大小递增的行列式 $ D_k $，其中 $ D_0 = 1 $，每个 $ D_k $ 是一个类似三对角的行列式，涉及奇数阶乘倒数。
• 通过从已知的伯努利数递推恒等式导出的线性系统，利用矩阵求逆与克莱姆法则，推导出第一个非递归公式 (2.2)。
• 应用生成函数技术与几何变换（通过 $ D(z) = \sum D_k z^{2k} $）证明公式 (2.2)、(2.3) 和 (2.4) 的等价性。
• 通过证明 $ D(z) $ 与 $ \frac{\sin z}{z} $ 的卷积结果为单位脉冲，建立生成函数恒等式 $ D(z) = \frac{z}{\sin z} $。
• 利用复变围道积分计算 $ D_p $ 作为 $ \frac{1}{z^{2p} \sin z} $ 的留数，从而导出积分表示与级数展开。
• 通过将 $ \frac{1}{\sin z} $ 展开为部分分式并提取留数，推导出闭式表达式 $ D_p = \frac{2}{\pi^{2p}} \left(1 - \frac{1}{2^{2p-1}} \right) \zeta(2p) $。

实验结果

研究问题

• RQ1 能否利用一组行列式推导出偶数索引伯努利数的显式非递归公式？
• RQ2 行列式序列 $ D_k $ 的闭式计算结果是什么？它与已知特殊函数有何关联？
• RQ3 通过变换技术，行列式序列与伯努利数的生成函数之间有何关系？
• RQ4 能否利用生成函数恒等式推导出多个 $ B_{2p} $ 的等价表达式，并证明其等价性？
• RQ5 $ D_p $ 的渐近行为如何？它与黎曼 zeta 函数有何联系？

主要发现

• 行列式序列 $ D_k $ 满足递推关系 $ D_k = \sum_{l=1}^k \frac{(-1)^{l-1} D_{k-l}}{(2l+1)!} $，其中 $ D_0 = 1 $，且对 $ k \geq 1 $ 有 $ \sum_{l=0}^k \frac{(-1)^l D_{k-l}}{(2l+1)!} = 0 $。
• 序列 $ D_k $ 的生成函数为 $ D(z) = \frac{z}{\sin z} $，这是系数 $ D_p $ 的已知生成函数。
• 显式公式 $ B_{2p} = \frac{(-1)^{p+1}(2p)! D_p}{2(2^{2p-1} - 1)} $ 提供了以 $ D_p $ 表示的 $ B_{2p} $ 的非递归表达式。
• 行列式 $ D_p $ 的计算结果为 $ D_p = \frac{2}{\pi^{2p}} \left(1 - \frac{1}{2^{2p-1}} \right) \zeta(2p) $，直接将 $ D_p $ 与偶数点处的黎曼 zeta 函数联系起来。
• 渐近地，$ D_p \sim \frac{2}{\pi^{2p}} $，且 $ \lim_{p \to \infty} \frac{\pi^{2p} D_p}{2} = 1 $，表明其模长迅速衰减。
• 通过生成函数卷积与留数分析，已证明所有三个公式 (2.2)、(2.3) 和 (2.4) 等价，确认了不同表示形式之间的一致性。