バーンスタイン多項式

提供: testwiki
ナビゲーションに移動 検索に移動

テンプレート:簡易区別 テンプレート:読み仮名バーンスタイン基底関数線形結合で与えられる多項式である。テンプレート:読み仮名[1]ベルンシュタイン多項式とも。

概要

バーンスタイン多項式はバーンスタイン基底関数線形結合で与えられる多項式である(⇒#定義)。テンプレート:Math 次のバーンスタイン多項式は任意の[[多項式#高々 n 次の多項式|高々 テンプレート:Mvar 次の多項式]]を表現できるため[2]、バーンスタイン多項式は多項式の別形式での表現であるといえる(⇒#特性[1]

バーンスタイン形式の数値的に安定な手法は、テンプレート:Illとして知られている。

バーンスタイン多項式はセルゲイ・ベルンシュテインが確率論を用いてワイエルシュトラスの近似定理の別証明をする際に初めて導入された(テンプレート:Harvnb)。のちに彼の名を取ってバーンスタイン多項式と呼ばれるようになった。

コンピュータ・グラフィックスの出現により、 テンプレート:Math の範囲におけるバーンスタイン多項式は、ベジェ曲線の重要な要素となった。

定義

バーンスタイン基底関数

n 次のテンプレート:読み仮名bν,n(x) は以下で定義される[3][注 1]

bν,n(x)=(nν)xν(1x)nν,ν=0,,n.

n 次のバーンスタイン基底関数は、[[多項式#高々 n 次の多項式|高々 テンプレート:Mvar 次の多項式]]からなるベクトル空間の基底をなす[2]

バーンスタイン多項式

n 次のバーンスタイン多項式 Bn(x) は以下で定義される:

Bn(x)=ν=0nβνbν,n(x)

すなわちバーンスタイン基底関数の線形結合であり、その係数 βνバーンスタイン係数あるいはベジェ係数と呼ばれる。

バーンスタイン基底関数は以下のような式となる。

b0,0(x)=1b0,1(x)=1xb1,1(x)=xb0,2(x)=(1x)2b1,2(x)=2x(1x)b2,2(x)=x2

特性

バーンスタイン基底関数の特性

バーンスタイン基底関数は以下のような特性を持つ。

b'ν,n(x)=n(bν1,n1(x)bν,n1(x))
  • 高次のバーンスタイン基底関数の和としても記述可能
bν,n1(x)=nνnbν,n(x)+ν+1nbν+1,n(x).

二項分布と区間同値

バーンスタイン基底関数は閉区間 [0,1]単位区間)において二項分布確率質量関数 (nk)pk(1p)nk と同値である。

区間正値

バーンスタイン基底関数は開区間 (0,1) において正の値のみをとる。すなわち bν,n(x)>0 for x(0,1)

区間非負

バーンスタイン基底関数は閉区間 [0,1]単位区間)において非負の値のみをとる。すなわち bν,n(x)0 for x[0,1]

これは 0x, 0(1x) for x[0,1] かつ二項係数が非負より明らかである。またこの関数がこの閉区間において二項分布と同値であることからも明らかである(⇒#二項分布と区間同値)。

1の分割

テンプレート:Math 次のバーンスタイン基底関数は1の分割をなす特性をもつ[4]

この特性は以下で示される:

ν=0nbν,n(x)=ν=0n(nν)xν(1x)nν=(x+(1x))n=1.

また、バーンスタイン基底関数が二項分布の確率質量関数と同じ形であることからもこの特性がわかる(⇒#二項分布と区間同値)。

バーンスタイン多項式の特性

バーンスタイン多項式は以下のような特性を持つ。

高々 n 次の多項式と同値

テンプレート:Math 次のバーンスタイン多項式は[[多項式#高々 n 次の多項式|高々 テンプレート:Mvar 次の多項式]]と同値である。言い換えれば、高々 テンプレート:Math 次の多項式からなるベクトル空間の任意の元を表現できる[2]

テンプレート:Math を例にとると、1次バーンスタイン多項式 B(x) は次のように展開できる:

B(x)=ν=01βνbν,1(x)=β0(1x)+β1x=(β1β0)x+β0

ここで傾き a・切片 b の高々1次の多項式 p(x)=ax+bB(x) が表現できるか考える。β0=b, β1=a+b とすると B(x) は、

B(x)=(a+bb)x+b=ax+b=p(x)

となる。ゆえに1次バーンスタイン多項式は高々1次の多項式と同値であるといえる。

この特性のためバーンスタイン多項式は多項式のバーンスタイン形式とも呼ばれる[1]

連続関数の近似

[0, 1] の範囲において連続な関数 テンプレート:Math を用いたバーンスタイン多項式

Bn(f)(x)=ν=0nf(νn)bν,n(x).

は、[0, 1] の範囲で以下のように、一様に収束する。

limnBn(f)(x)=f(x)

このことは、各点収束するが一様収束はしないという命題に比べ、より強い命題である。この一様収束は、以下のように明確に示される。

limnsup0x1|f(x)Bn(f)(x)|=0.

上述のように、バーンスタイン多項式はワイエルシュトラスの近似定理の証明にも用いられる。

また、より一般的に、連続な k 次導関数についても、

Bn(f)(k)(n)knkf(k) and f(k)Bn(f)(k)0,

であることが示せる。ここで (n)knk=(10n)(11n)(1k1n)Bn固有値である。

limnBn(f)(x)=f(x)であることの初等的な説明

bν,n(x)=(nν)xν(1x)nν

は確率 テンプレート:Math で事象 テンプレート:Math が起こる試行を テンプレート:Math 回繰り返したとき、事象 テンプレート:Math がちょうどテンプレート:Math回起こる確率を表す。試行をn回繰り返す場合において、テンプレート:Mathテンプレート:Math 回起こったときに得られる確率変数をテンプレート:Mathとすると、

Bn(f)(x)=ν=0nf(νn)bν,n(x).

は期待値を表す。 一方、テンプレート:Math 回試行を繰り返す場合、事象 テンプレート:Math が起こる回数は平均して テンプレート:Math である。よって、平均して得られる確率変数、すなわち期待値は テンプレート:Math であると考えられる。今 テンプレート:Mathテンプレート:Math は整数で、テンプレート:Mathテンプレート:Math を分母とする有理数とは限らないので テンプレート:Mathテンプレート:Math の誤差も テンプレート:Math とは限らないが、テンプレート:Math を大きくしていくと両者の誤差は テンプレート:Math に近づいていくと考えられるので、

limnBn(f)(x)=f(x)

が成り立つ。

脚注

注釈

テンプレート:Reflist

出典

テンプレート:Reflist

参考文献

関連項目

外部リンク

テンプレート:PlanetMath attribution

  1. 1.0 1.1 1.2 "In computer aided geometric design, polynomials are usually expressed in Bernstein form. ... Let p(t) ... be a polynomial in the Bernstein form" p.744,746 より引用。テンプレート:Cite journal
  2. 2.0 2.1 2.2 テンプレート:Cite book
  3. "バーンスタイン基底関数 xin=nCiti(1t)ni n は次数を表す" テンプレート:Harvnb より引用。
  4. "P(t)=i=0nBinPiBin=nCiti(1t)ni ある比率で各制御点の座標を混ぜ合わせる ... 混合比(和は 1 になる) 混合比を関数で表したものを「基底関数」とよぶ" テンプレート:Harvnb より引用。


引用エラー: 「注」という名前のグループの <ref> タグがありますが、対応する <references group="注"/> タグが見つかりません