Skip to content

方法Method

Bernstein 多项式与构造型 Weierstrass 逼近

Bernstein polynomial approximation · Bernstein polynomials · Weierstrass approximation theorem · Weierstrass 逼近定理

用非负二项式权重构造连续函数的多项式逼近,并给出连续模、Lipschitz及二阶光滑情形的一致误差证书。

形式陈述 ​

给定 f∈C([0,1],R) 或 C([0,1],C),整数 n≥1,定义

bn,k(x)=(nk)xk(1−x)n−k,Bnf(x)=∑k=0nf(k/n)bn,k(x).

端点值按多项式连续延拓解释,故 Bnf(0)=f(0)、Bnf(1)=f(1)。Bnf 是次数至多 n 的多项式,其系数直接来自 n+1 个等距样本。

构造型 Weierstrass 定理:Bnf 一致收敛到 f。因此闭区间上任意连续实值或复值函数都可由同标量域的多项式一致逼近。对任意 δ>0,有可计算的通用界

‖Bnf−f‖∞≤ωf(δ)+‖f‖∞2nδ2,

其中 ωf 是连续模。特别地,若 f 为 L-Lipschitz,则

‖Bnf−f‖∞≤L2n.

这些是充分误差上界,不宣称 Bernstein 构造在给定次数下达到最佳一致误差。

直觉

固定 x 时,权重 bn,k(x) 非负且和为一。Bnf(x) 是若干采样值的加权平均,主要权重集中在 k/n 接近 x 的位置。连续模限制近处样本与 f(x) 的差,权重的二阶偏差则控制远处样本总共能占多少分量。

若用概率语言,Bnf(x)=E[f(Sn/n)],其中 Sn 服从参数 (n,x) 的二项分布。但下面只用有限和恒等式,构造和证明都不要求先学习概率论。

三个权重恒等式 ​

二项式展开给出 ∑kbn,k=1。利用 k(nk)=n(n−1k−1) 得到

∑kknbn,k(x)=x.

同样由 k(k−1)(nk)=n(n−1)(n−2k−2),得到

∑k(kn−x)2bn,k(x)=x(1−x)n≤14n.

n=1 时第二阶下降阶乘和为零,同一公式仍成立。三个恒等式分别表达总质量、一阶位置和离散程度。

由近远分解证明一致逼近 ​

在误差和 ∑k[f(k/n)−f(x)]bn,k(x) 中,将索引按 |k/n−x|≤δ 与大于 δ 分组。近组的函数差至多为 ωf(δ);远组的差至多为 2‖f‖∞,而远组总权重满足

∑|k/n−x|>δbn,k(x)≤δ−2∑k(k/n−x)2bn,k(x)≤14nδ2.

这就给出形式陈述中的误差界。为使误差小于给定 ε,先取 δ 使连续模小于 ε/2,再取 n 使第二项小于 ε/2,两次选择都与 x 无关。

若有 Lipschitz 常数,使用有限加权和的Cauchy–Schwarz 不等式可避免近远分组:

|Bnf(x)−f(x)|≤L∑k|k/n−x|bn,k(x)≤L[∑k(k/n−x)2bn,k(x)]1/2≤L2n.

平方根来自二阶偏差,常数 1/2 来自 x(1−x) 的最大值。

例子与边界

给尖角一个完整多项式证书 ​

取 v(x)=|x−1/2|。反三角不等式给它的 Lipschitz 常数为一,故

pn(x)=∑k=0n|kn−12|(nk)xk(1−x)n−k,‖pn−v‖∞≤12n.

要误差至多 10−2,取 n=2500 即得一个有明确全部系数的多项式证书。若要求严格小于该值,可取 n>2500。不需要把2501项展开成容易误读的单项式系数表。

低次情形可完全手算。n=4 时采样值依次为 1/2,1/4,0,1/4,1/2,因此

p4(x)=12−x+2x3−x4.

在 0≤x≤1/2,误差为 p4−v=2x3−x4,其导数 6x2−4x3≥0;对称性给另一半区间的误差。因此实际最大误差在 x=1/2 取得,恰为 3/16,小于通用上界 1/4。

采样不等于插值 ​

对 f(x)=x2,上述矩公式直接给

Bnf(x)=x2+x(1−x)n.

内点通常不满足 Bnf(k/n)=f(k/n),尽管构造使用了这些样本。多项式插值要求精确穿过所有节点,本构造通过非负平均取得统一稳定与收敛保证,两者不是同一运算。

连续性也不能取消。跳跃函数若被连续多项式一致逼近,一致极限应连续,产生矛盾。对一般连续函数,次数还不能只由 ‖f‖∞ 与容差决定:固定 n 后,可在相邻网点之间放高为一而所有采样值为零的连续尖峰,此时 Bnf=0 而误差为一;其连续模随着尖峰变窄而变差。

推论与应用

按容差构造并求值 ​

已知 L 和目标误差 ε>0 时,可取

n=max{1,⌈L2/(4ε2)⌉},

采样 ak=f(k/n),并以这 n+1 个数保存 Bernstein 形式。一次安全避免显式大二项系数的精确算术求值方式是逐层凸组合:初始 ak(0)=ak,对 r=1,…,n 更新

ak(r)=(1−x)ak(r−1)+xak+1(r−1),0≤k≤n−r.

输出 a0(n)。归纳不变量为

ak(r)=∑j=0r(rj)(1−x)r−jxjak+j,

由二项系数递推即可从第 r−1 层推出第 r 层,最终正是 Bnf(x)。共有 n(n+1)/2 次组合,成本 O(n2) 算术运算、O(n) 存储;从小到大的下标原地更新不会覆盖尚需读取的右邻值。

若每个样本误差至多 η,非负且和为一的权重使最终采样误差也至多 η。因此总证书可加上 η。浮点舍入误差仍需另行预算,数学一致误差界没有自动包含机器算术误差。

其他正则性与区间 ​

连续模的区间缩放性质配合同一二阶矩给出

‖Bnf−f‖∞≤32ωf(n−1/2).

确实,ωf(t)≤(1+nt)ωf(n−1/2),取加权和并使用上面的平均距离界即得。若 f∈C2[0,1] 且 ‖f″‖∞≤A,在 x 处作带二阶余项的展开,一阶项因权重均值等于 x 而抵消,从而误差至多 A/(8n)。

在非退化闭区间 [a,b] 上,先对 F(t)=f(a+(b−a)t) 构造 BnF,再代入 t=(x−a)/(b−a),仍为 x 的多项式。原函数若有 Lipschitz 常数 L,新的证书为 L(b−a)/(2n)。Stone–Weierstrass 定理进一步刻画一般紧空间上的函数代数何时稠密;Korovkin 定理则抽出这里的三个权重矩,形成正线性逼近的通用检验。

参考资料
关系图谱15 个相邻概念 · 2 类关系

拖动节点调整位置。

显示关系

显示:依赖

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