Polinomio de Bernstein
Los polinomios de Bernstein o polinomios en la base de Bernstein son una clase particular de polinomios en el campo de los números reales, que son utilizados dentro del ámbito del análisis numérico. El nombre hace referencia al matemático ucraniano Sergei Natanovich Bernstein.
El algoritmo de evaluación más numéricamente estable es el de
De Casteljau.
Índice
1 Definición
2 Propiedades
3 Ejemplo
4 Aplicaciones
5 Véase también
Definición
Un polinomio de Bernstein P(x){displaystyle P(x),} de orden n aproxima una función f(x){displaystyle f(x),} en un intervalo, mejor cuanto mayor sea n, a partir de esta fórmula:
- P(x)=∑i=0nciBin(x){displaystyle P(x)=sum _{i=0}^{n}{c_{i}B_{i}^{n}(x)}}
donde los Bin(x){displaystyle B_{i}^{n}(x)} son elementos de la distribución binomial respecto de la variable x{displaystyle x,} y los ci{displaystyle c_{i},} son valores de la función que queremos aproximar.
Para aproximar la función en el intervalo [0,1]{displaystyle [0,1],} estos elementos toman los siguientes valores:
- ci=f(in)yBin(x)=(ni)xi(1−x)n−i{displaystyle c_{i}=fleft({frac {i}{n}}right)qquad {text{y}}qquad B_{i}^{n}(x)={n choose i}x^{i}(1-x)^{n-i}}
(aquí (ni){displaystyle {n choose i}} es el coeficiente binomial).
y más en general transformando las ecuaciones para un intervalo [a,b]{displaystyle [a,b],}, los Bin(x)[a,b]{displaystyle B_{i}^{n}(x)_{[a,b]}} se convierten en polinomios de la base de Bernstein:
- ci=f(ib−an+a)yBin(x)[a,b]=(ni)(x−a)i(b−x)n−i(b−a)n{displaystyle c_{i}=fleft(i,{frac {b-a}{n}}+aright)qquad {text{y}}qquad B_{i}^{n}(x)_{[a,b]}={n choose i}{(x-a)^{i}(b-x)^{n-i} over (b-a)^{n}}}
Así, la fórmula general desarrollada es:
- P(x)=∑i=0nf(ib−an+a)n!i!(n−i)!(x−a)i(b−x)n−i(b−a)n{displaystyle P(x)=sum _{i=0}^{n}{fleft(i,{frac {b-a}{n}}+aright){frac {n!}{i!(n-i)!}}{frac {(x-a)^{i}(b-x)^{n-i}}{(b-a)^{n}}}}}
Propiedades
Para un grado n, existen n+1 polinomios de Bernstein B0n,…,Bnn{displaystyle B_{0}^{n},dots ,B_{n}^{n}} definidos sobre el intervalo [a,b]{displaystyle [a,b],} , por
- Bin(x)[a,b]=(ni)(x−a)i(b−x)n−i(b−a)n{displaystyle B_{i}^{n}(x)_{[a,b]}={n choose i}{(x-a)^{i}(b-x)^{n-i} over (b-a)^{n}}}
Estos polinomios presentan estas propiedades importantes, que cumplen para cualquier valor de x{displaystyle x,} en el intervalo [a,b]{displaystyle [a,b],}
Partición de la unidad : ∑i=0nBin(x)=1{displaystyle qquad sum _{i=0}^{n}B_{i}^{n}(x)=1}
Positividad : Bin(x)≥0,∀i∈0…n{displaystyle B_{i}^{n}(x)geq 0,qquad forall iin 0dots n}
Simetría : Bin(x)=Bn−in(1−x),∀i∈0…n{displaystyle B_{i}^{n}(x)=B_{n-i}^{n}(1-x),qquad forall iin 0dots n}
Las dos primeras propiedades nos indican que forman una combinación convexa.
La modificación por escala y traslación de intervalo no influye sobre los coeficientes del polinomio en cuestión.
Se ha de notar la gran semejanza de estos polinomios con la distribución binomial.
Para el intervalo [0,1]{displaystyle [0,1],} existe esta fórmula de recurrencia:
Bin(x)={(1−x)Bin−1(x),i=0(1−x)Bin−1(x)+xBi−1n−1(x),i=1…n−1xBi−1n−1(x),i=n{displaystyle B_{i}^{n}(x)={begin{cases}(1-x)B_{i}^{n-1}(x),&i=0\(1-x)B_{i}^{n-1}(x)+xB_{i-1}^{n-1}(x),&i=1dots n-1\xB_{i-1}^{n-1}(x),&i=nend{cases}}}.
Ejemplo
En el caso de un polinomio de orden 2{displaystyle 2} la base en [0,1]{displaystyle [0,1],} está compuesta de:
- B02(x)=(20)x0(1−x)2−0=(1−x)2{displaystyle B_{0}^{2}(x)={2 choose 0}x^{0}(1-x)^{2-0}=(1-x)^{2}}
- B12(x)=(21)x1(1−x)2−1=2x(1−x){displaystyle B_{1}^{2}(x)={2 choose 1}x^{1}(1-x)^{2-1}=2x(1-x)}
- B22(x)=(22)x2(1−x)2−2=x2{displaystyle B_{2}^{2}(x)={2 choose 2}x^{2}(1-x)^{2-2}=x^{2}}
Un polinomio expresado en esta base tendría entonces la forma:
- P(x)=c0B02(x)+c1B12(x)+c2B22(x)=f(0)(1−x)2+2f(12)x(1−x)+f(1)x2{displaystyle P(x)=c_{0}B_{0}^{2}(x)+c_{1}B_{1}^{2}(x)+c_{2}B_{2}^{2}(x)=f(0)(1-x)^{2}+2fleft({frac {1}{2}}right)x(1-x)+f(1)x^{2}}
Si aproximamos f1(x)=x{displaystyle f_{1}(x)=x,} obtenemos el mismo polinomio: P1(x)=x{displaystyle P_{1}(x)=x,}
si evaluamos f2(x)=x2{displaystyle f_{2}(x)=x^{2},} aproxima a: P2(x)=x2+x2{displaystyle P_{2}(x)={frac {x^{2}+x}{2}},}
y probando con f3(x)=ex{displaystyle f_{3}(x)=e^{x},} resulta: P3(x)=(1−x)2+2ex(1−x)+ex2≈ 0.421x2+1.29x+1{displaystyle P_{3}(x)=(1-x)^{2}+2{sqrt {e}},x(1-x)+ex^{2}approx 0.421x^{2}+1.29x+1,}
Aplicaciones
Los polinomios de Bernstein son utilizados para demostrar el teorema de aproximación de Weierstrass y por esto son también utilizados para efectuar aproximaciones e interpolaciones de funciones como, por ejemplo, la curva de Beziér, así como para la estimación de las funciones de densidad de probabilidad:
Para n que tiende al infinito, el polinomio converge uniformamente hacia la función f (x), o sea
- |Bn(x)−f(x)|≤5/4 ω(f,1/n){displaystyle |B_{n}(x)-f(x)|leq 5/4 omega (f,1/{sqrt {n}})}
donde
ω(f,δ)=sup|h|≤δ|f(x+h)−f(x)|{displaystyle omega (f,delta )=sup _{|h|leq delta }|f(x+h)-f(x)|}, llamado módulo de continuidad.
Véase también
- Algoritmo de "De Casteljau"
- Curva de Beziér
Aproximación de Bernstein, permite aproximar uniformemente funciones continuas.