Teorema de Taylor

De testwiki
Ir a la navegación Ir a la búsqueda
La función exponencial y=ex (línea roja continua) y su aproximación mediante un polinomio de Taylor alrededor del origen de coordenadas (línea verde discontinua)

En cálculo diferencial, el teorema de Taylor, recibe su nombre del matemático británico, Brook Taylor, quien lo enunció con mayor generalidad en 1712, aunque previamente James Gregory lo había descubierto en 1671.[1] Este teorema permite obtener aproximaciones polinómicas de una función en un entorno de cierto punto en el que la función sea diferenciable. Además el teorema permite acotar el error obtenido mediante dicha estimación.

Caso de una variable

Enunciado del teorema

La versión más básica del teorema es como sigue:

Plantilla:Teorema

El polinomio que aparece en el teorema de Taylor,

Pk(x)=f(a)+f(a)(xa)+f(a)2!(xa)2++f(k)(a)k!(xa)k,

se denomina polinomio de Taylor de orden k de la función f en el punto a. El polinomio de Taylor es el único polinomio que «mejor aproxima la función en forma asintótica», en el sentido de que si existe una función hk: y un polinomio p de orden k tal que

f(x)=p(x)+hk(x)(xa)k,limxahk(x)=0,

entonces p=Pk. El teorema de Taylor describe el comportamiento asintótico del término del resto

 Rk(x)=f(x)Pk(x),

el cual es el error de aproximación cuando se aproxima f con su polinomio de Taylor. Utilizando la notación o el teorema de Taylor se puede expresar de la siguiente forma:

Rk(x)=o(|xa|k),xa.

Fórmulas explícitas para el resto

Existen diferentes formas de expresar Rk(x) que se mencionan a continuación: Plantilla:Definición

Usualmente, esta refinación del teorema de Taylor, se demuestra con el teorema del valor medio, de ahí su nombre. También se pueden hallar expresiones similares. Por ejemplo, si G(t) es continua en el intervalo cerrado y diferenciable con derivadas no nulas en el intervalo abierto entre a y x, entonces

Rk(x)=f(k+1)(ξ)k!(xξ)kG(x)G(a)G(ξ)

para algún número ξ entre a y x. Esta versión generaliza las formas de Lagrange y Cauchy del resto, que son tomadas como casos especiales, y se demuestran usando el teorema del valor medio de Cauchy.

En el caso de la forma integral del resto, se requieren conceptos de la teoría integral de Lebesgue para una generalidad completa. Sin embargo, se mantiene el concepto que provee la integral de Riemann donde la derivada (k+1)-ésima de f es continua en el intervalo cerrado [a,x]. Plantilla:Definición

Debido a la continuidad absoluta de f(k) en el intervalo cerrado entre a y x su derivada f(k+1) existe como una función L1, y el resultado puede probarse con un cálculo formal usando el teorema fundamental del cálculo e integración por partes.

Para algunas funciones f(x), se puede probar que el resto, Rn(f), se aproxima a cero cuando n se acerca al infinito; dichas funciones pueden ser expresadas como series de Taylor en un entorno reducido alrededor de un punto a y son denominadas funciones analíticas.

El teorema de Taylor con Rn(f) expresado de la segunda forma es también válido si la función f tiene números complejos o valores vectoriales. Además existe una variación del teorema de Taylor para funciones con múltiples variables.

Acotación del resto

Suele ser muy útil en la práctica acotar el término del resto de la aproximación de Taylor, en lugar de tener la fórmula exacta de este. Suponiendo que f es continuamente diferenciable k+1 veces en un intervalo I que contiene a a. Suponemos que hay constantes q y Q tal que

qf(k+1)(x)Q

en el intervalo I. Entonces el término del resto satisface la desigualdad[2]

q(xa)k+1(k+1)!Rk(x)Q(xa)k+1(k+1)!.

Si x>a, y similarmente si x<a. Esta es una consecuencia simple de la forma de Lagrange del resto. En particular, si

|f(k+1)(x)|M

sobre un intervalo I=(ar,a+r) con algún r>0 , entonces

|Rk(x)|M|xa|k+1(k+1)!Mrk+1(k+1)!

para todo x(ar,a+r). A la segunda desigualdad se la llama acotación uniforme, porque permanece uniformemente para todo x sobre el intervalo (ar,a+r).

Ejemplo

Aproximación de ex (azul) por su polinomio de Taylor Pk de orden Plantilla:Esd centrado en x = 0 (rojo)

Suponiendo que se desea aproximar la función f(x)=ex en el intervalo [1,1] con un error no mayor a 105. Este ejemplo necesita que se conozcan las siguientes propiedades de la función exponencial:

Plantilla:Ecuación

De estas propiedades se tiene que f(k)(x)=ex para todo k, y en particular, f(k)(0)=1. Entonces el polinomio de Taylor de orden k de f en 0 y su resto bajo la forma de Lagrange son:

Pk(x)=1+x+x22!++xkk!,Rk(x)=eξ(k+1)!xk+1,

donde ξ es algún número entre 0 y x. Ya que ex es creciente Plantilla:Eqnref, podemos usar simplemente que ex1 para x[1,0] para acotar el resto sobre el subintervalo [1,0]. Para obtener una cota superior para el resto en [1,0], usamos la propiedad eξ<ex para 0<ξ<x para acotar

ex=1+x+eξ2x2<1+x+ex2x2,0<x1

usando la expansión de Taylor de segundo orden. Entonces resolvemos ex para deducir que

ex1+x1x22=21+x2x24,0x1

simplemente maximizando el numerador y minimizando el denominador. Combinando estas acotaciones para ex vemos que

|Rk(x)|4|x|k+1(k+1)!4(k+1)!,1x1,

así se alcanza la precisión requerida, donde

4(k+1)!<1054105<(k+1)!k9.

(ver factorial o calcular manualmente los valores 9!=362 880 y 10!=3 628 800). Como conclusión, el teorema de Taylor permite la aproximación

ex=1+x+x22!++x99!+R9(x),|R9(x)|<105,1x1.

Luego, esta aproximación nos da la expresión decimal e ≈ 2,71828, correcta hasta cinco dígitos decimales.

Demostración

Sea[3]

hk(x)={f(x)P(x)(xa)kx=a0x=a

donde, como dice en el enunciado del teorema de Taylor,

P(x)=f(a)+f(a)(xa)+f(a)2!(xa)2++f(k)(a)k!(xa)k.

Es suficiente mostrar que

limxahk(x)=0.

La demostración de Plantilla:Eqnref se basa en la aplicación repetida de la regla de L'Hôpital. Se observa que, para cada j=0,1,...,k1, f(j)(a)=P(j)(a). De aquí que cada una de los primeras k1 derivadas del numerador en hk(x) se anula en x=a, y lo mismo sucede con el denominador. También, ya que la condición de que la función f sea k veces diferenciable en un punto requiere diferenciabilidad de orden k1 en un entorno de dicho punto (esto es así, porque la diferenciabilidad requiere una función definida en un entorno del punto), el numerador y sus k2 derivadas son diferenciables en un entorno de a. Claramente, el denominador también satisface dicha condición, y adicionalmente, no se anula a menos que x=a, por lo tanto se satisfacen todas las condiciones para la regla de L'Hôpital, y así se justifica su utilización. Por lo tanto,

limxaf(x)P(x)(xa)k=limxaddx(f(x)P(x))ddx(xa)k==limxadk1dxk1(f(x)P(x))dk1dxk1(xa)k=1k!limxaf(k1)(x)P(k1)(x)xa=1k!(f(k)(a)f(k)(a))=0

donde queda la anteúltima igualdad por la definición de la derivada en x=a.

Obtención de la forma de valor medio del resto

Sea G una función real, continua sobre un intervalo cerrado entre a y x y diferenciable con derivadas no nulas sobre el intervalo abierto entre a y x, y la función que se define como

F(t)=f(t)+f(t)(xt)+f(t)2!(xt)2++f(k)(t)k!(xt)k.

Entonces, por el teorema del valor medio de Cauchy,

(*)F(ξ)G(ξ)=F(x)F(a)G(x)G(a)

para algún ξ sobre el intervalo abierto entre a y x. Se observa que el numerador F(x)F(a)=Rk(x) es exactamente el resto del polinomio de Taylor para f(x). Calculando

F(t)=f(t)+(f(t)(xt)f(t))+(f(3)(t)2!(xt)2f(2)(t)1!(xt))++(f(k+1)(t)k!(xt)kf(k)(t)(k1)!(xt)k1)=f(k+1)(t)k!(xt)k,

reemplazando en (*) y reorganizando los términos para hallar que

Rk(x)=f(k+1)(ξ)k!(xξ)kG(x)G(a)G(ξ).

Esta es la forma del término que mencionamos como «resto», después enunciamos el teorema de Taylor con el resto bajo la forma del valor medio. La forma de Lagrange del resto puede obtenerse haciendo G(t)=(tx)k+1 y la forma de Cauchy haciendo G(t)=ta.

Observación. Usando este método también se puede recurrir a la forma integral del resto haciendo

G(t)=atf(k+1)(s)k!(xs)kds,

pero los requerimientos de f necesitados para usar el teorema del valor medio son más fuertes, si se tiene el objetivo de probar el caso en que f(k) es únicamente continua absolutamente. Sin embargo, si se usa la integral de Riemann en vez de la integral de Lebesgue, los requerimientos no pueden ser tan débiles.

Obtención de la forma integral del resto

Debido a la continuidad absoluta de f(k) sobre el intervalo cerrado entre a y x su derivada f(k+1) existe como una función L1, y se usa el teorema fundamental del cálculo y la integración por partes. Esta misma demostración se aplica para la integral de Riemann teniendo en cuenta que f(k) es continua sobre el intervalo cerrado y diferenciable sobre el intervalo abierto entre a y x, y esto permite llegar al mismo resultado que cuando se utilizó el teorema del valor medio.

El teorema fundamental del cálculo dice que

f(x)=f(a)+axf(t)dt.

A partir de aquí se usa la integración por partes y se usa una vez más el teorema fundamental del cálculo para ver que

f(x)=f(a)+(xf(x)af(a))axtf(t)dt=f(a)+x(f(a)+axf(t)dt)af(a)axtf(t)dt=f(a)+(xa)f(a)+ax(xt)f(t)dt,

que es exactamente el teorema de Taylor con resto en la forma integral para el caso k=1. La enunciación general se demuestra usando la inducción. Suponiendo que

(*)f(x)=f(a)+f(a)1!(xa)++f(k)(a)k!(xa)k+axf(k+1)(t)k!(xt)kdt.

Integrando el término del resto por partes se llega a que

axf(k+1)(t)k!(xt)kdt=[f(k+1)(t)(k+1)k!(xt)k+1]ax+axf(k+2)(t)(k+1)k!(xt)k+1dt= f(k+1)(a)(k+1)!(xa)k+1+axf(k+2)(t)(k+1)!(xt)k+1dt.

Substituyendo esto en la fórmula (*) se muestra que si se tiene para el valor k, debe obtenerse también para el valor k+1. Por lo tanto, ya que se tiene para k=1, se tiene para cualquier valor entero positivo k.

Caso de varias variables

El teorema de Taylor anterior Plantilla:Eqnref puede generalizarse al caso de varias variables como se explica a continuación. Sea B una bola en n centrada en a, y f una función real definida sobre la clausura B¯ cuyas derivadas parciales hasta de orden n+1 son todas continuas en cada punto de la bola. El teorema de Taylor establece que para cualquier xB: Plantilla:Ecuación donde la suma se extiende sobre los multi-índices α (esta fórmula usa la notación multi-índice). El resto satisface la desigualdad: Plantilla:Ecuación para todo α con |α|=n+1. Tal como sucede en el caso de una variable, el resto puede expresarse explícitamente en términos de derivadas superiores (véase la demostración para los detalles).

Demostración

Para demostrar el teorema de Taylor para el caso multidimensional, considérese un función f:n o campo escalar, que suponemos continuo y, para simplificar lo expuesto (aunque una generalización es trivial), de clase C. Sea 𝐫(t) una función vectorial que va de n, y definámosla como 𝐫(t)=𝐚+𝐮t (de ahora en adelante, se omitirán las flechas de los vectores). Pongamos 𝐫(t)=𝐲 Ahora hagamos g(t)=f[𝐫(t)] y recordemos que g(t)=f(𝐲)𝐫(t). Notemos ahora que: Plantilla:Ecuación Ahora, derivando sucesivas veces, encontramos que podemos poner de forma muy cómoda: Plantilla:Ecuación donde el producto escalar del vector 𝐮 con el gradiente Plantilla:Ecuación representa la derivada direccional, y el exponente n sobre ella es entendido como las sucesivas veces que la hacemos sobre la función; es decir, hacemos la derivada direccional n veces. Ahora, empleando el teorema de Taylor para una variable real, expandimos g(t) en su serie de McLaurin: Plantilla:Ecuación y haciendo t=1 y sustituyendo las derivadas por las expresiones antes hallada se evidencia que: Plantilla:Ecuación Obsérvese que el primer término aparece el gradiente y en el segundo la matriz hessiana, pero escrito con esta notación particular que resulta más cómoda y compacta. La expresión obtenida es equivalente a la expresada más arriba mediante la notación multiíndice.

Conjunto Ox,αn(h) de operadores fraccionales

El cálculo fraccional de conjuntos (Fractional Calculus of Sets (FCS)), mencionado por primera vez en el artículo titulado "Sets of Fractional Operators and Numerical Estimation of the Order of Convergence of a Family of Fractional Fixed-Point Methods",[4] es una metodología derivada del cálculo fraccional.[5] El concepto principal detrás del FCS es la caracterización de los elementos del cálculo fraccional utilizando conjuntos debido a la gran cantidad de operadores fraccionales disponibles.[6][7][8] Esta metodología se originó a partir del desarrollo del método de Newton-Raphson fraccional [9] y trabajos relacionados posteriores.[10][11][12]

Ilustración de algunas líneas generadas por el método de Newton–Raphson fraccional para la misma condición inicial x0 pero con diferentes órdenes α del operador fraccional implementado. Fuente: Applied Mathematics and Computation

El cálculo fraccional, una rama de las matemáticas que trata con derivadas de orden no entero, surgió casi simultáneamente con el cálculo tradicional. Esta emergencia fue en parte debido a la notación de Leibniz para derivadas de orden entero: dndxn. Gracias a esta notación, L'Hopital pudo preguntar en una carta a Leibniz sobre la interpretación de tomar n=12 en una derivada. En ese momento, Leibniz no pudo proporcionar una interpretación física o geométrica para esta pregunta, por lo que simplemente respondió a L'Hopital en una carta que «... es una aparente paradoja de la cual, algún día, se derivarán consecuencias útiles».

El nombre «cálculo fraccional» se origina a partir de una pregunta histórica, ya que esta rama del análisis matemático estudia derivadas e integrales de un cierto orden α. Actualmente, el cálculo fraccional carece de una definición unificada de lo que constituye una derivada fraccional. En consecuencia, cuando no es necesario especificar explícitamente la forma de una derivada fraccional, típicamente se denota de la siguiente manera:

dαdxα.

Los operadores fraccionales tienen varias representaciones, pero una de sus propiedades fundamentales es que recuperan los resultados del cálculo tradicional a medida que αn. Considerando una función escalar h:m y la base canónica de m denotada por {e^k}k1, el siguiente operador fraccional de orden α se define utilizando notación de Einstein:[13]

oxαh(x):=e^kokαh(x).

Denotando kn como la derivada parcial de orden n con respecto al componente k-ésimo del vector x, se define el siguiente conjunto de operadores fraccionales:

Ox,αn(h):={oxα:okαh(x) y limαnokαh(x)=knh(x) k1},

cuyo complemento es:

Ox,αn,c(h):={oxα:okαh(x) k1 y limαnokαh(x)knh(x) para al menos un k1}.

Como consecuencia, se define el siguiente conjunto:

Ox,αn,u(h):=Ox,αn(h)Ox,αn,c(h).

Extensión a funciones vectoriales

Para una función h:Ωmm, el conjunto se define como:

mOx,αn,u(h):={oxα:oxαOx,αn,u([h]k) km},

donde [h]k:Ωm denota el k-ésimo componente de la función h.

Conjunto mSx,αn,γ(h) de operadores fraccionales

Sea 0 el conjunto {0}. Si γ0m y xm, entonces es posible definir la siguiente notación multi-índice:

{γ!:=k=1m[γ]k!,|γ|:=k=1m[γ]k,xγ:=k=1m[x]k[γ]kγxγ:=[γ]1[x]1[γ]2[x]2[γ]m[x]m.

Entonces, considerando una función h:Ωm y el operador fraccional:

sxαγ(oxα):=o1α[γ]1o2α[γ]2omα[γ]m,

se define el siguiente conjunto de operadores fraccionales:

Sx,αn,γ(h):={sxαγ=sxαγ(oxα) : sxαγh(x)  con  oxαOx,αs(h) sn2  y  limαksxαγh(x)=kγxkγh(x) α,|γ|n}.

De donde se obtienen los siguientes resultados:

Si sxαγSx,αn,γ(h)  {limα0sxαγh(x)=o10o20om0h(x)=h(x)limα1sxαγh(x)=o1[γ]1o2[γ]2om[γ]mh(x)=γxγh(x) |γ|nlimαqsxαγh(x)=o1q[γ]1o2q[γ]2omq[γ]mh(x)=qγxqγh(x) q|γ|qnlimαnsxαγh(x)=o1n[γ]1o2n[γ]2omn[γ]mh(x)=nγxnγh(x) n|γ|n2.

Como consecuencia, considerando una función h:Ωmm, se define el siguiente conjunto de operadores fraccionales:

mSx,αn,γ(h):={sxαγ : sxαγSx,αn,γ([h]k) km}.

Conjunto mTx,α,γ(a,h) de operadores fraccionales

Considerando una función h:Ωmm y el siguiente conjunto de operadores fraccionales:

mSx,α,γ(h):=limnmSx,αn,γ(h).

Entonces, tomando una bola B(a;δ)Ω, es posible definir el siguiente conjunto de operadores fraccionales:

mTx,α,γ(a,h):={txα,=txα,(sxαγ) : sxαγmSx,α,γ(h)  y  txα,h(x):=|γ|=01γ!e^jsxαγ[h]j(a)(xa)γ},

el cual permite generalizar la expansión en serie de Taylor de una función vectorial en notación multi-índice. Como consecuencia, es posible obtener el siguiente resultado:

Si txα,mTx,α,γ(a,h){txα,h(x)=e^j[h]j(a)+|γ|=11γ!e^jsxαγ[h]j(a)(xa)γ+|γ|=21γ!e^jsxαγ[h]j(a)(xa)γ=h(a)+k=1ne^jokα[h]j(a)[(xa)]k+|γ|=21γ!e^jsxαγ[h]j(a)(xa)γ.

Referencias

Plantilla:Listaref

Bibliografía

Plantilla:Control de autoridades