Bernsteinpolynoom

De bernsteinpolynomen (genoemd naar Sergej Natanovitsj Bernstein) zijn een familie speciale reële polynomen met geheeltallige coëfficiënten.

Betekenis en geschiedenis

De bernsteinpolynomen vinden hun oorsprong in de approximatietheorie. Hun bedenker, Bernstein, kon met behulp van deze polynomen een constructief bewijs leveren voor de stelling van Stone-Weierstrass]. Aan het einde van de jaren vijftig werden voor het eerst pogingen ondernomen methoden op basis van bernsteinpolynomen te gebruiken bij het ontwerpen van krommen en oppervlakken. Paul de Faget de Casteljau bij Citroën en Pierre Bézier bij Renault gebruikten de polynomen bij hun ontwikkeling van de béziercurven en legden zo de basis van het huidige Computer Aided Design (CAD).

Definitie

Voor n N 0 {\displaystyle n\in \mathbb {N} _{0}} zijn de n + 1 {\displaystyle n+1} bernsteinpolynomen van graad n {\displaystyle n} de reële polynomen

B k , n ( t ) = ( n k ) t k ( 1 t ) n k {\displaystyle B_{k,n}(t)={n \choose k}\,t^{k}\,(1-t)^{n-k}}

met 0 k n {\displaystyle 0\leq k\leq n} en t [ 0 , 1 ] {\displaystyle t\in [0,1]} .

Door affiene transformaties van het interval [ 0 , 1 ] {\displaystyle [0,1]} naar een interval [ a , b ] {\displaystyle [a,b]} ontstaan de gegeneraliseerde bernsteinpolynomen:

B k , n [ a , b ] ( t ) = 1 ( b a ) n ( n k ) ( t a ) k ( b t ) n k {\displaystyle B_{k,n}^{[a,b]}(t)={\frac {1}{(b-a)^{n}}}{n \choose k}(t-a)^{k}\,(b-t)^{n-k}}

voor t [ a , b ] {\displaystyle t\in [a,b]}

Hierin is ( n k ) {\displaystyle {n \choose k}} een binomiaalcoëfficiënt.

De eerste bernsteinpolynomen B k , n ( t ) {\displaystyle B_{k,n}(t)}
k {\displaystyle k} 0 {\displaystyle 0} 1 {\displaystyle 1} 2 {\displaystyle 2} 3 {\displaystyle 3}
n {\displaystyle n}
0 1 {\displaystyle 1}
1 1 t {\displaystyle 1-t} t {\displaystyle t}
2 ( 1 t ) 2 {\displaystyle (1-t)^{2}} 2 t ( 1 t ) {\displaystyle 2t(1-t)} t 2 {\displaystyle t^{2}}
3 ( 1 t ) 3 {\displaystyle (1-t)^{3}} 3 t ( 1 t ) 2 {\displaystyle 3t(1-t)^{2}} 3 t 2 ( 1 t ) {\displaystyle 3t^{2}(1-t)} t 3 {\displaystyle t^{3}}

Voorbeeld

De bernsteinpolynomen B_ {k, 4}

De afbeelding toont de bernsteinpolynomen B i , 4 {\displaystyle B_{i,4}} , 0 k 4 {\displaystyle 0\leq k\leq 4} van graad 4 {\displaystyle 4} op het interval [ 0 , 1 ] {\displaystyle [0,1]} .

Eigenschappen

De bernsteinpolynomen op het interval [ 0 , 1 ] {\displaystyle [0,1]} hebben de volgende eigenschappen:

  • Basiseigenschap: De bernsteinpolynomen { B k , n : 0 k n } {\displaystyle \{B_{k,n}:0\leq k\leq n\}} zijn lineair onafhankelijk en vormen een basis van Π n {\displaystyle \Pi _{n}} , de ruimte van de polynomen van graad kleiner of gelijk aan n {\displaystyle n} .
  • Positiviteit: B k , n ( t ) > 0 {\displaystyle B_{k,n}(t)>0} voor alle t ( 0 , 1 ) {\displaystyle t\in (0,1)} .
  • Extrema: De polynoom B k , n {\displaystyle B_{k,n}} heeft precies één (absoluut) maximum op het interval [ 0 , 1 ] {\displaystyle [0,1]} in het punt t = k / n {\displaystyle t=k/n} . In het bijzonder is dus:
B 0 , n ( 0 ) = B n , n ( 1 ) = 1 {\displaystyle B_{0,n}(0)=B_{n,n}(1)=1}
  • Opdeling van 1: De bernsteinpolynomen van graad n {\displaystyle n} zijn de termen in de binomiale ontwikkeling:
1 = ( t + ( 1 t ) ) n = k = 0 n B k , n ( t ) {\displaystyle 1=(t+(1-t))^{n}=\sum _{k=0}^{n}B_{k,n}(t)}
  • Symmetrie: B k , n ( t ) = B n k , n ( 1 t ) {\displaystyle B_{k,n}(t)=B_{n-k,n}(1-t)}
  • Recursie:
B k , n ( t ) = ( 1 t ) B k , n 1 ( t ) + t B k 1 , n 1 ( t ) {\displaystyle B_{k,n}(t)=(1-t)\,B_{k,n-1}(t)+t\,B_{k-1,n-1}(t)}
B k , n = 0 {\displaystyle B_{k,n}=0} voor k < 0 {\displaystyle k<0} en k > n {\displaystyle k>n}
B 0 , 0 = 1 {\displaystyle B_{0,0}=1}
  • Teruglopend:
B k , n ( t ) = k + 1 n + 1 B k + 1 , n + 1 ( t ) + n + 1 k n + 1 B k , n + 1 ( t ) {\displaystyle B_{k,n}(t)={\frac {k+1}{n+1}}B_{k+1,n+1}(t)+{\frac {n+1-k}{n+1}}B_{k,n+1}(t)}
  • Afgeleiden:
B k , n ( t ) = n ( B k 1 , n 1 ( t ) B k , n 1 ( t ) ) {\displaystyle B'_{k,n}(t)=n(B_{k-1,n-1}(t)-B_{k,n-1}(t))}
B 1 , n 1 = B n , n 1 = 0 {\displaystyle B_{-1,n-1}=B_{n,n-1}=0}
  • Primitieve:
B k , n ( t ) d t = 1 n + 1 i = k + 1 n + 1 B i , n + 1 ( t ) {\displaystyle \int B_{k,n}(t)\,\mathrm {d} t={\frac {1}{n+1}}\sum \limits _{i=k+1}^{n+1}B_{i,n+1}(t)}

Benadering door bernsteinpolynomen

Voor een functie f : [ 0 , 1 ] R {\displaystyle f\colon [0,1]\to \mathbb {R} } heet de polynoom B n ( f ) {\displaystyle B_{n}(f)} gedefinieerd door

B n ( f ) ( t ) = k = 0 n B k , n ( t ) f ( k n ) {\displaystyle B_{n}(f)(t)=\sum _{k=0}^{n}B_{k,n}(t)\cdot f\left({\frac {k}{n}}\right)}

de n {\displaystyle n} -de bernsteinpolynoom van f {\displaystyle f} .

Als f {\displaystyle f} een continue functie is op het interval [ 0 , 1 ] {\displaystyle [0,1]} , convergeert de rij van zijn bernsteinpolynomen B n ( f ) {\displaystyle B_{n}(f)} uniform naar f {\displaystyle f} .

Het bewijs hiervan kan onder andere geleverd wordenmet behulp van de zwakke wet van de grote getallen.

Voorbeeld

Benaderingen van de functie f {\displaystyle f} (rood) door bernsteinpolymomen van graad 4 (blauw) en van graad 10 (geel)

De benadering van de functie

f ( t ) = { 2 t  voor  0 t 1 2 2 ( 1 t )  voor  1 2 t 1 {\displaystyle f(t)={\begin{cases}2t&{\text{ voor }}0\leq t\leq {\tfrac {1}{2}}\\2(1-t)&{\text{ voor }}{\tfrac {1}{2}}\leq t\leq 1\end{cases}}}

door bernsteinpolynomen van de graad 4 is de bernsteinpolynoom van f {\displaystyle f} :

B 4 ( f ) ( t ) = k = 0 n B k , 4 ( t ) f ( k 4 ) = {\displaystyle B_{4}(f)(t)=\sum _{k=0}^{n}B_{k,4}(t)\cdot f\left({\frac {k}{4}}\right)=}
= 1 2 ( 4 t ( 1 t ) 3 + 4 t 3 ( 1 t ) ) + 6 t 2 ( 1 t ) 2 = {\displaystyle ={\tfrac {1}{2}}\left(4t(1-t)^{3}+4t^{3}(1-t)\right)+6t^{2}(1-t)^{2}=}
= 2 t 4 t 3 + 2 t 4 {\displaystyle =2t-4t^{3}+2t^{4}}

In de figuur staat de functie f {\displaystyle f} en de benaderingen voor n = 4 {\displaystyle n=4} en n = 10 {\displaystyle n=10} .

Literatuur

  • Bernstein, S.N., Démonstration du théorème de Weierstrass fondée sur le calcul des probabilités, Commun. soc. Math. Charkov, deel 12, nr. 2, blz. 1-2, 1912/1913.

Weblinks

  • Bernsteinpolynomen (applet)