“设 $a<b$ 为整数,$p\ge1$ 为整数,$f\in C^{2p}[a,b]$。以Bernoulli数与周期多项式的 $B 1= 1/2$ 约定,有 $$ \sum {k=a}^b f…”
求和与积分相差的,不只是首尾两个点。要把更细的差别分层记录,需要一组在相邻整数处恰好接得上的多项式。Bernoulli多项式把导数、单位步差分和端点修正组织在一起。
形式陈述
在
系数
最初几项为
它们满足
导数关系加上零均值和
直觉
普通积分只确定原函数到一个常数。式(3)每次都把这个常数调到“一个单元内平均为零”,这样周期复制后不会凭空加入固定偏置。式(2)则说明它对导数和单位步差分各保留一种简洁结构。
从同一个生成函数推导三条规则
对式(1)求
在
将
差分式在
例子与边界
精确生成系数,而非记一张表
由
因此可逐项算出
恢复多项式。计算到阶
周期化只复制单元,不能抹去接缝
记
特别地
所以
这些是整段解析界,因为
另一些文献采用第一Bernoulli数
推论与应用
差分式使
例如
本页新的接口在周期核。Euler–Maclaurin公式逐单元使用式(2),让内部端点抵消,剩下有限端点导数差与一个准确积分余项。B₁的跳跃贡献采样值,较高阶的端点匹配让过程继续,式(5)再把余项转成证书。
参考资料
- NIST DLMF,§24.2(i)、(iii),式24.2.1、24.2.3及24.2.11–12:数、多项式生成函数与周期延拓。
- NIST DLMF,§24.4,式24.4.1、24.4.3、24.4.7及导数小节:差分、对称和幂和。本文从生成函数重建这些关系,并单独给出B₄的全区间因式分解界。