Skip to content

定义Definition

Bernoulli 多项式与周期余项核

Bernoulli polynomials · Bernoulli numbers · 伯努利多项式 · 伯努利数

由生成函数构造Bernoulli多项式,证明导数、差分和端点规范,区分B₁符号约定,并以周期延拓构造求和误差核。

求和与积分相差的,不只是首尾两个点。要把更细的差别分层记录,需要一组在相邻整数处恰好接得上的多项式。Bernoulli多项式把导数、单位步差分和端点修正组织在一起。

形式陈述 ​

在 t=0 把可去奇点补齐,以幂级数定义

(1)textet−1=∑n=0∞Bn(x)tnn!,|t|<2π.

系数 Bn(x) 是关于 x 的 n 次多项式,称Bernoulli多项式;Bn:=Bn(0) 称Bernoulli数。本页固定 B1=−1/2 约定。

最初几项为

B0(x)=1,B1(x)=x−12,B2(x)=x2−x+16,B3(x)=x3−32x2+12x,B4(x)=x4−2x3+x2−130.

它们满足

(2)Bn′(x)=nBn−1(x),Bn(x+1)−Bn(x)=nxn−1(n≥1),(3)∫01Bn(x)dx=0(n≥1).

导数关系加上零均值和 B0=1 也唯一确定这一族:先积分得到下一阶,再由均值条件确定积分常数。

直觉

普通积分只确定原函数到一个常数。式(3)每次都把这个常数调到“一个单元内平均为零”,这样周期复制后不会凭空加入固定偏置。式(2)则说明它对导数和单位步差分各保留一种简洁结构。

从同一个生成函数推导三条规则 ​

对式(1)求 x 导数,右端系数整体乘 t,比较同次项即得导数式。把 x 换成 x+1 后相减,左端变成 text,比较系数得到差分式。

在 x∈[0,1] 上逐项积分,左端为1,故除常数项外的均值全为零。这里固定 |t|<2π 的较小闭圆盘后,级数在紧区间上一致收敛,才可逐项积分。

将 x 换成 1−x、t 换成 −t,得到

Bn(1−x)=(−1)nBn(x).

差分式在 x=0 又给 Bn(1)=Bn(0)(n≠1)。两者合并,所有奇数 n≥3 的 Bn 都为0;n=1 例外,两个端点分别为 −1/2 与 1/2。

例子与边界

精确生成系数,而非记一张表 ​

由 [t/(et−1)](et−1)=t 比较系数,得到

(4)B0=1,∑k=0n(n+1k)Bk=0(n≥1).

因此可逐项算出 B1=−1/2,B2=1/6,B3=0,B4=−1/30,再用

Bn(x)=∑k=0n(nk)Bkxn−k

恢复多项式。计算到阶 N 需要 O(N2) 次有理数乘加;分子分母的位数会增长,这个计数不是固定机器字的位复杂度保证。

周期化只复制单元,不能抹去接缝 ​

记 {x}=x−⌊x⌋,定义 B~n(x)=Bn({x})。当 n≥2 时,两个端点值相等,所以周期函数连续;B~1 则在整数处跳跃。每个开整数单元内仍满足导数规则,跨接缝分部积分时必须逐单元记边界项。

特别地

B4(u)=u2(1−u)2−130(0≤u≤1),

所以

(5)−130≤B4(u)≤7240,−116≤B4−B4(u)=−u2(1−u)2≤0.

这些是整段解析界,因为 0≤u(1−u)≤1/4;不是在有限网格上观察到的最大值。它们会直接控制求和余项的大小和符号。

另一些文献采用第一Bernoulli数 +1/2 的记号。不能仅替换一个符号而保留式(1)和相同求和端点;本页的生成函数已经固定了约定。

推论与应用

差分式使 Bp+1(x)/(p+1) 成为 xp 的离散原函数,从而

∑k=0N−1kp=Bp+1(N)−Bp+1(0)p+1.

例如 p=1 给 N(N−1)/2,核对的是左闭右开求和;改成1到N必须改端点。旧有限差分页用二项式基解决一般多项式求和,这里不另提供重复的幂和算法。

本页新的接口在周期核。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₄的全区间因式分解界。
关系图谱7 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

  1. 前置三跳
  2. 前置二跳
  3. 前置一跳
  4. 当前条目
  5. 后续一跳
  6. 后续二跳
  7. 后续三跳
文字版关系按与当前条目的最短距离分组
类型化关系