Polinomios de Bernstein

Por Lizbeth Fernández Villegas

MATERIAL EN REVISIÓN

Introducción

Es sabido que existen funciones que no es tan sencillo evaluar en todos los puntos de su dominio. Sin embargo, cuando la función f es nveces derivable en un punto a podemos definir polinomios de Taylor Tn,a (ver entrada Cálculo Diferencial e Integral I: Polinomios de Taylor (Parte 1)). Conforme aumentamos el grado del polinomio más nos acercamos al valor real de f, incluso tenemos resultados en relación al residuo de Taylor (ver Cálculo Diferencial e Integral I: Polinomios de Taylor (Parte 2). Pero, ¿qué podemos hacer ante una función que no es derivable en ningún punto? En esta sección presentaremos una forma alternativa de estimar funciones continuas, útil para aquellas en las que no podemos identificar los polinomios de Taylor, pues recordemos que la continuidad no es una condición suficiente para que una función sea derivable. Por ejemplo, la función de Weierstrass dada por f(x)=n=0ancos(bnπx)   con 0<a<1 y ab>1+32π,b>1 entero impar,
es continua en R pero no es diferenciable en ningún punto.

La función de Weierstrass, es continua y no diferenciable en ningún punto.

Las funciones con las que nos vamos a aproximar son conocidas como polinomios de Bernstein. Una forma de entenderlos es a través de la probabilidad. A continuación presentamos ideas tomadas del artículo de la Dra. Ana Meda cuya lectura sugerimos:
Meda A. (2005) Interpolar con volados, o los Polinomios de Bernstein. Miscelánea Matemática, 41, 1-12.

Sea f:[0,1]R una función continua y nN. Supón además que conocemos los valores que toma esta función en los puntos 0n,1n,2n,3n, y sea x[0,1]. Para saber cuánto vale f en x vamos a tomar una moneda cuya probabilidad de arrojar águila al lanzarse sea precisamente x. Ahora la vamos a lanzar n veces mientras registramos las veces en que salió águila. Sea k ese número de resultados. Identifiquemos el valor

(1)f(kn).

Definimos φk:[0,1]R, como
(2)φk(x)=(nk)xk(1x)nk.

Con
(nk)=n!k!(nk)!.

Corresponde a la probabilidad de que k de los lanzamientos sean águila con esa moneda cargada de probabilidad x.

Con (1) y (2) definimos

(3)Bn(f,x)=k=0nf(kn)(nk)xk(1x)nk,

que es la esperanza de la variable aleatoria que acabamos de describir, debido a que corresponde a la suma de los valores que toma esta variable ponderada por la probabilidad de que los tome.

A lo largo de esta entrada mostraremos formalmente que esta estimación funciona. Comenzamos diciendo qué es «acercarse mucho» a una función continua. Presentamos la definición con polinomios, pues son funciones continuas y derivables, lo cual facilita su manejo.

Definición. Función continua aproximada por polinomios. Sea fC0([0,1],R). Si para cada ε>0 existe un polinomio Pε:[0,1]R tal que

fPε=Supx[0,1]{|f(x)Pε(x)|}<ε,

diremos que la función f puede aproximarse por polinomios. Demostraremos que toda función continua en C0([0,1],R) tiene esta propiedad. Específicamente, los polinomios que usaremos en esa aproximación están dados por la siguiente:

Definición. El nésimo polinomio de Bernstein. Sea fC0([0,1],R). El nésimo polinomio de Bernstein de f de grado a lo más n es:
Bn(f,x)=k=0nf(kn)(nk)xk(1x)nk.

Que es la igualdad (3) definida arriba.

Mostremos el nésimo polinomio de Bernstein para tres funciones particulares.

El nésimo polinomio de Bernstein para la función constante, f(x)=1.

Del teorema del binomio sabemos que
(s+t)n=k=0n(nk)sktnk

Haciendo s=x y t=1x tenemos que
1=k=0n(nk)xk(1x)nk

Si consideramos f(x)=1 entonces f(kn)=1 y así
1=k=0n1(nk)xk(1x)nk=k=0nf(kn)(nk)xk(1x)nk

Por lo tanto
Bn(f(x)=1,x)=1.

El nésimo polinomio de Bernstein para la función identidad, f(x)=x.

Partiendo de
1=k=0n(nk)xk(1x)nk

Reemplacemos n por n1 y k por j para obtener
(4)1=j=0n1(n1j)xj(1x)n(1+j)

El siguiente resultado será usado en el cálculo del polinomio:

(n1k1)=(n1)!(k1)!(nk)!=nnkk(n1)!(k1)!(nk)!=knn(n1)!k(k1)!(nk)!=knn!k!(nk)!=kn(nk)

Por lo tanto

(5)(n1k1)=kn(nk)

Multipliquemos por x ambos lados de la igualdad (4) y apliquemos la igualdad (5) en el coeficiente para tener

x=j=0n1j+1n(nj+1)xj+1(1x)n(1+j)
Haciendo k=j+1 y usando que f(x)=x entonces f(kn)=kn. Se sigue:
(6)x=k=1nkn(nk)xk(1x)nk=k=0nf(kn)(nk)xk(1x)nk

Por lo tanto
Bn(f(x)=x,x)=x.

El nésimo polinomio de Bernstein para la función identidad, f(x)=x2.

Partiendo de
1=k=0n(nk)xk(1x)nk

Reemplacemos n por n2 y k por j para obtener

(7)1=j=0n2(n2j)xj(1x)n(2+j)

El siguiente resultado será usado en el cálculo del polinomio:

(n2k2)=(n2)!(k2)!(nk)!=nnn1n1kkk1k1(n2)!(k2)!(nk)!=k(k1)n(n1)n(n1)(n2)!k(k1)(k2)!(nk)!=k(k1)n(n1)n!k!(nk)!=k(k1)n(n1)(nk)

Por lo tanto

(8)(n2k2)=k(k1)n(n1)(nk)

Multipliquemos por x2 ambos lados de la igualdad (7) y apliquemos la igualdad (8) en el coeficiente para tener
x2=j=0n2(j+2)(j+1)n(n1)(nj+2)xj+2(1x)n(j+2)
Haciendo k=j+2 se sigue que:
x2=k=2nk(k1)n(n1)(nk)xk(1x)nk=k=0nk(k1)n(n1)(nk)xk(1x)nk

Entonces
(9)x2(n2n)=k=0nk(k1)(nk)xk(1x)nk

Al dividir ambos lados de la igualdad por n2 se obtiene:
(11n)x2=k=0nk(k1)n2(nk)xk(1x)nk=k=0n(k2n2kn2)(nk)xk(1x)nk=k=0nk2n2(nk)xk(1x)nkk=0nkn2(nk)xk(1x)nk=k=0n(kn)2(nk)xk(1x)nk1nk=0nkn(nk)xk(1x)nk=k=0n(kn)2(nk)xk(1x)nk1nx

De modo que
(11n)x2+1nx=k=0n(kn)2(nk)xk(1x)nk=k=0nf(kn)(nk)xk(1x)nk

Por lo tanto
Bn(x,f(x)=x2)=(11n)x2+1nx.

Ahora daremos paso al

Teorema de aproximación de Bernstein. Sea f:[0,1]R continua. Entonces la sucesión de polinomios de Bernstein para f converge uniformemente a f en [0,1].

Demostración:
Partimos de
1=k=0n(nk)xk(1x)nk
Al multiplicar ambos lados de la igualdad por f(x) obtenemos

f(x)=k=0nf(x)(nk)xk(1x)nk.

Por lo tanto

f(x)Bn(f,x)=k=0n(f(x)f(kn))(nk)xk(1x)nk.

Usando la desigualdad del triángulo y el hecho de que x[0,1] implica que (nk)xk(1x)nk>0, se sigue que

(10)|f(x)Bn(f,x)|k=0n|(f(x)f(kn))|(nk)xk(1x)nk

Nuestra intención ahora será mostrar que cuando n tiende a infinito esta distancia se hace muy pequeña.

Sea ε>0. Como f es continua en [0,1] entonces es uniformemente continua, así existe δ>0 tal que
(11) si |st|<δ entonces |f(s)f(t)|<ε2.

Sea x[0,1]. Considera nN tal que
n máx {1δ4,f2ε2}.

Separemos los elementos de N{0} en dos conjuntos como sigue:

(12)I1={kN{0}:0kn,|knx|<(1n)14},(13)I2={kN{0}:0kn,kI1},

Como n1δ4 entonces
n141δδ(1n)14

De modo que si kI1 satisface |knx|<(1n)14δ y por la desigualdad (11) se cumple que

kI1|(f(x)f(kn))|(nk)xk(1x)nkkI1ε2(nk)xk(1x)nkε2k=0(nk)xk(1x)nk(14)=ε2.

En cuanto a los ksI2 tenemos lo siguiente:

Si kI2 entonces
(1n)14|xkn|1n|xkn|4(xkn)4n(xkn)2n.

Usaremos la última desigualdad al final del siguiente grupo de ecuaciones.

kI2|(f(x)f(kn))|(nk)xk(1x)nkkI22f(nk)xk(1x)nk=2fkI2(xkn)2(xkn)2(nk)xk(1x)nk(15)2fnkI2(xkn)2(nk)xk(1x)nk

Ahora, partiendo de 1=k=0n(nk)xk(1x)nk, tenemos:

(x2)1=k=0n(x2)(nk)xk(1x)nk(16)x2=k=0nx2(nk)xk(1x)nk

Por la igualdad (9) tenemos:

(17)x2(n2n)=k=0nk(k1)(nk)xk(1x)nk

Partiendo de (6) se sigue:

(2x)x=k=0n(2x)kn(nk)xk(1x)nk(18)2x2=k=0n2knx(nk)xk(1x)nk

Y también que

(n)x=k=0n(n)kn(nk)xk(1x)nk(19)nx=k=0nk(nk)xk(1x)nk

Sumando (17) con (19).

x2(n2n)+nx=k=0n(k(k1)+k)(nk)xk(1x)nkx2(n2n)+nx=k=0nk2(nk)xk(1x)nk(1n2)(x2(n2n)+nx)=k=0n(1n2)k2(nk)xk(1x)nk(20)(11n)x2+1nx=k=0nk2n2(nk)xk(1x)nk

Ahora sumemos (16), (18) y (20) para obtener:

(21)1n(xx2)=k=0n(xkn)2(nk)xk(1x)nk

Dado que máx x[0,1](xx2)=14 podemos concluir que

(22)k=0n(xkn)2(nk)xk(1x)nk=1n(xx2)(1n)(14)=14n

Ya que nf2ε2 podemos concluir que εnf. Se sigue de (15) y (22) que:

kI2|(f(x)f(kn))|(nk)xk(1x)nk2fnkI2(xkn)2(nk)xk(1x)nk2εnn14n(23)=ε2.

Ahora, de (14) y (23) tenemos, finalmente que

|f(x)Bn(f,x)|k=0n|(f(x)f(kn))|(nk)xk(1x)nk=kI1|(f(x)f(kn))|(nk)xk(1x)nk+kI2|(f(x)f(kn))|(nk)xk(1x)nkε2+ε2(24)=ε.

Lo cual demuestra que la sucesión de polinomios de Bernstein converge uniformemente a f.

Si bien, el teorema anterior aplica para funciones con dominio en [0,1] se puede generalizar a cualquier intervalo cerrado en R según enunciamos a continuación.

Teorema. De aproximación de Weierstrass. Sea f:[a,b]R una función continua. Entonces existe una sucesión de polinomios que converge uniformemente a f en [a,b].

Demostración:
Definimos ρ:[0,1][a,b] donde para cada x[0,1],ρ(x):=(1x)a+xb. En consecuencia, la función g:=fρ es una función continua en [0,1] y por el teorema anterior sabemos que la sucesión de polinomios de Bernstein de g converge en g, es decir, Bn(g,x)g.

Ahora definimos pn:=Bn(g,x)ρ1. Nota que es un polinomio. Ocurre que

pnf:=máxt[a,b]|pn(t)f(t)|=máxx[0,1]|Bn(g,x)(x)g(x)|=Bn(g,x)g.

Por lo tanto se da la convergencia uniforme pnf en C0[a,b].

Más adelante…

Conoceremos condiciones bajo las que es posible acercarnos a funciones continuas que tienen su dominio en conjuntos más generales que un intervalo cerrado: los espacios compactos. Lo haremos a través del teorema de Stone-Weierstrass que enunciaremos en la siguiente sección.

Tarea moral

  1. Considera la función φk definida en (2). Demuestra que φk alcanza su máximo en el punto kn.
  2. Desarrolla las funciones φk para k=0,1,2 para n=2.
  3. Considera la función f:[0,1]R dada por f(x)=sen(πx). Grafica los polinomios de Bernstein de f para n=1,2,4,8. Visualiza la gráfica para f.

Enlaces

Deja una respuesta

Tu dirección de correo electrónico no será publicada. Los campos obligatorios están marcados con *

Este sitio usa Akismet para reducir el spam. Aprende cómo se procesan los datos de tus comentarios.