Error Taylor

PRCTICA
POLINOMIOS DE TAYLOR. ACOTACIN DEL RESTO
Prcticas Matlab
Prctica 4
Objetivos
Analizar cmo incide en la aproximacin mediante polinomios de Taylor el considerar

distintas funciones.
Mostrar cmo acotar el error que se comente al aproximar el valor de una funcin en
un punto por su polinomio de Taylor utilizando la expresin del resto de Lagrange.
Comandos
Los comandos que se utilizarn ya se han visto en prcticas anteriores.
Ejemplos resueltos
Consideremos f (x ) = x , g (x ) = (1.5) y supongamos que

x
deseamos obtener de forma aproximada el valor de 1.5 utilizando

polinomios de Taylor. Observamos que f (1.5) = g (0.5) = 1.5
Se pide:
(a) Obtener los polinomios de Taylor de rdenes 2, 3, 4 en a=1 de la
funcin f(x) y de g(x).
(b) Rellenar la siguiente tabla sustituyendo los polinomios en los

puntos que se indican
Polinomios de
f (x ) = x a = 1
T2
T3
T4
T2
T3
T4
x = 1.5
Polinomios de
g (x ) = (1.5) a = 1
x
x = 0.5
Qu polinomio de Taylor aproxima mejor al valor 1.5 ?
PGINA 2
MATLAB: POLINOMIOS DE TAYLOR. ACOTACIN DEL RESTO
Una aproximacin de
1.5 con 20 cifras decimales es:
1.2247448713915890491. Puedes utilizar Matlab para calcular este
valor tecleando:
>>vpa(sqrt(1.5),20)
Solucin:
(a) Cdigo Matlab
syms x
% Funciones a utilizar
f=sqrt(x);
g=(1.5)^x;
% Punto en el que se desarrolla el polinomio
a=1.0;
% Polinomios de Taylor de grados 3, 4 y 5
polinomiosf=[taylor(f,3,a) taylor(f,4,a) taylor(f,5,a)]
polinomiosg=[taylor(g,3,a) taylor(g,4,a) taylor(g,5,a)]
Cdigo Matlab
Al cdigo del apartado anterior habr que aadir:
(b)
% Puntos a analizar
punf=1.5;
pung=0.5;
% Sustitucin de los polinomios de Taylor de f en puntf
valorPolinomiof=subs(polinomiosf,x,punf);
% Sustitucin de los polinomos de Taylor de f en puntg
valorPolinomiog=subs(polinomiosg,x,pung);
% Mostramos los valores de las funciones y de los polinomios
disp('---------------Valores funciones--------------------')
disp('
x
f(x)')
disp([punf
subs(f,x,punf)])
disp('
x
g(x)')
disp([pung
subs(g,x,pung)])
disp('---------------Valores polinomios--------------------')
disp('
grado
Taylor_f
Taylor_g')
disp([[2 3 4]' valorPolinomiof' valorPolinomiog'])
Se considera la funcin f (x ) =
1
1+x
(a) Calcula para esta funcin y y la diferencial en a = 0 para
PGINA 3
distintos valores de x (considera x = 0.5 , x = 0.1 , x = 0.01 )

y = f (0 + x ) f (0)
dy = f ' (0) x
y dy
x = 0.5
x = 0.1
x = 0.01
Nota: Como puedes observar en la tabla anterior se tiene que y dy

para valores pequeos de x .
(b) Podras obtener un valor aproximado de f (0 + x ) realizando
nicamente operaciones aritmticas elementales y conociendo el valor
de f (0) y el de f ' (0) ? Podras estimar el error de esta aproximacin?
(c) Aplica el apartado (b) al caso en el que x = 0.5 y escribe una

cota del error que se comente al aproximar f (0.5) =
f (0 ) +
2
por
3
f ' (0 ) 1
. Da un intervalo en el que puedas asegurar que est el
1!
2
valor exacto de
2
.
3
Solucin:
(a) Cdigo Matlab
syms x
f=1/sqrt(1+x);
a=0;
n=1;
%Clculo a+incremento
incrementox=[0.5 0.1
0.01];
x=a+incrementox;
% Calculo del incremento y de la diferencial para los valores de h
incrementoy=subs(f,x)-subs(f,a);
derivada=diff(f,1);
der=subs(derivada,a);
diferencialy=der*incrementox;
%Escribimos en forma de tabla
disp([incrementox' incrementoy' diferencialy'])
PGINA 4
(b) S, el valor de f (0 + x ) se podra obtener de forma aproximada como

f (0) + f ' (0) x .
La diferencia entre estos dos valores es el resto de orden 1:

f '' (t )
2
R1 =
x )
(
2!
siendo t un punto intermedio a 0 y x .
Si calculamos un valor M de manera que
f '' (t ) M
para t comprendido entre 0 y x
se tendra que una cota del error de la aproximacin sera:

error = R1 =
En
el caso particular
de
que
f '' (t )
2!
(x )
f (x ) =
2
M
x )
(
2!
1
1+x
la
derivada
segunda
es
5/2
3
que al ser una funcin decreciente y positiva toma su valor ms
1 + x)
(
4
grande en valor absoluto en el punto 0 , por lo que bastara tomar M = f '' (0)
f '' (x ) =
(c) Observamos que
f ' (0 )
2
= f (0.5) y que T1 (0.5) = f (0) +
0.5 = 0.75 .
3
1!
Calculamos una cota del error que se comete al decir que
2
0.75 .
3
Cdigo Matlab
syms x
f=1/sqrt(1+x);
a=0;
n=1;
puntox=0.5;
%Para acotar el resto de Lagrange dibujamos la derivada
%de orden n+1 y vemos donde alcanza su valor mximo
dominio=a:0.01:puntox;
figure(1)
plot(dominio,subs(diff(f,n+1),dominio))
legend('Derivada n+1')
%Acotamos la expresin del resto acotando la
%derivada n+1 en t por M. Para este ejemplo
%la derivada n+1 toma el valor mayor, en valor
%absoluto, en el punto 0
M=abs(subs(diff(f,n+1),0));
%Escribimos una cota del Resto de Lagrange en valor absoluto
cotaResto=M/factorial(n+1)*(puntox-a).^(n+1);
PGINA 5
%Calculamos el valor de la funcin y de su polinomio en el punto

%para ver la bondad de la cota del error calculada
valorFuncion=subs(f,puntox);
valorPolinomio=subs(taylor(f,n+1),puntox);
error=abs(valorFuncion-valorPolinomio);
disp('-------------------------------------------------')
disp([ '
x
f(x)
T(x)
error
cotaResto'])
disp([puntox' valorFuncion'
valorPolinomio' error'
cotaResto'])
Una cota del error sera 0.0938. El intervalo en el que se encuentra el valor exacto de
2 / 3 es el siguiente:
0.75 0.0938, 0.75 + 0.0938 = 0.6562, 0.8438 .
(a) De forma anloga al apartado anterior, podras realizar una

estimacin del error de la aproximacin de f (x ) =
1
1+x
por su
polinomio de Taylor de grado 2 en el punto a = 0 cuando x pertenece

al intervalo 0 x
1
?
2
(b) Cul sera el grado del polinomio de Taylor que se necesitara

utilizar para aproximar
2
con un error menor que 108 ?
3
(a) Bastara sustituir en el cdigo del apartado (c) del ejercicio anterior puntox por un
1
vector con valores 0 x . Por ejemplo: puntox=0:0.1:0.5;
2
Observa que como la derivada tercera es creciente pero tomando valores negativos el
valor ms grande en valor absoluto se encuentra en el punto 0 . Es decir, M = f ''' (0)
Si adems incluyes el siguiente cdigo puedes representar el error en valor absoluto

entre el valor de la funcin y su polinomio que da el ordenador y el valor que hemos
calculado como cota de dicho error.
figure(2)
plot(puntox,error,'b-',puntox,cotaResto,'ro-')
legend('|f(x)-T(x)|','Cota del resto')
PGINA 6
(d) El grado del polinomio es 23. Basta ir aumentando el valor de n hasta comprobar
que cotaResto es menor que 108 .
Cdigo Matlab
syms x
f=1/sqrt(1+x);
a=0;
punto=0.5;
n=23;
polinomio=taylor(f,n+1,a);
derivadaN1=diff(f,n+1);
M=abs(subs(derivadaN1,x,0));
cotaResto=M/factorial(n+1)*(punto-a)^(n+1);
format long
disp([n cotaResto])
Importante: Observa que la cota de la derivada n+1 que se considera en la acotacin

del resto de Lagrange siempre es la derivada n+1 en el 0 ya que
(1) (2n 1) !!
n
(n +1
(t ) =
(1 + t )
(n +1)
2
con
0 t 0.5
f (n +1 (t ) f (n +1 (0) = M
Otra forma:
syms x
f=1/sqrt(1+x);
a=0;
punto=0.5;
n=0;
cotaResto=1;
while cotaResto>(10^(-8))
n=n+1;
polinomio=taylor(f,n+1,a);
derivadaN1=diff(f,n+1);
M=abs(subs(derivadaN1,x,0));
cotaResto=M/factorial(n+1)*(punto-a)^(n+1);
end
disp('-------------------------------------------------')
disp(['El grado del polinomio es ' num2str(n)])
PGINA 7
Ejercicios propuestos
(1) Dibujar las grficas de las siguientes funciones en 0,1
( )
(a )
f (x ) = s en e x
(b )
f (x ) = cos log x
3
2
(2) Determinar el error que se comete al aproximar el valor de f (1) al

considerar el polinomio de Taylor de grado 3 en el punto a = 0.5 para
las dos funciones anteriores.
Determinar el grado mnimo del polinomio de Taylor alrededor del

punto a = 0.5 para cada una de las siguientes funciones de tal
manera que el error que se cometa en la aproximacin de f (1) sea de
4 cifras significativas
( )
(a )
f (x ) = s en e x
(b )
f (x ) = cos log x
3
2

Error Taylor

Cargado por

Copyright:

Formatos disponibles

Error Taylor

Cargado por

Información del documento

Descripción original:

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Error Taylor

Cargado por

Copyright:

Formatos disponibles

PRCTICA

POLINOMIOS DE TAYLOR. ACOTACIN DEL RESTO

Analizar cmo incide en la aproximacin mediante polinomios de Taylor el considerar

Consideremos f (x ) = x , g (x ) = (1.5) y supongamos que

deseamos obtener de forma aproximada el valor de 1.5 utilizando

(b) Rellenar la siguiente tabla sustituyendo los polinomios en los

MATLAB: POLINOMIOS DE TAYLOR. ACOTACIN DEL RESTO

(a) Calcula para esta funcin y y la diferencial en a = 0 para

MATLAB: POLINOMIOS DE TAYLOR. ACOTACIN DEL RESTO

distintos valores de x (considera x = 0.5 , x = 0.1 , x = 0.01 )

Nota: Como puedes observar en la tabla anterior se tiene que y dy

(c) Aplica el apartado (b) al caso en el que x = 0.5 y escribe una

MATLAB: POLINOMIOS DE TAYLOR. ACOTACIN DEL RESTO

(b) S, el valor de f (0 + x ) se podra obtener de forma aproximada como

La diferencia entre estos dos valores es el resto de orden 1:

para t comprendido entre 0 y x

se tendra que una cota del error de la aproximacin sera:

(c) Observamos que

Calculamos una cota del error que se comete al decir que

MATLAB: POLINOMIOS DE TAYLOR. ACOTACIN DEL RESTO

%Calculamos el valor de la funcin y de su polinomio en el punto

(a) De forma anloga al apartado anterior, podras realizar una

polinomio de Taylor de grado 2 en el punto a = 0 cuando x pertenece

(b) Cul sera el grado del polinomio de Taylor que se necesitara

Si adems incluyes el siguiente cdigo puedes representar el error en valor absoluto

MATLAB: POLINOMIOS DE TAYLOR. ACOTACIN DEL RESTO

Importante: Observa que la cota de la derivada n+1 que se considera en la acotacin

MATLAB: POLINOMIOS DE TAYLOR. ACOTACIN DEL RESTO

(1) Dibujar las grficas de las siguientes funciones en 0,1

(2) Determinar el error que se comete al aproximar el valor de f (1) al

Determinar el grado mnimo del polinomio de Taylor alrededor del

También podría gustarte