INSTITUTO POLITCNICO NACIONAL
UNIDAD PROFESIONAL INTERDISCIPLINARIA DE BIOTECNOLOGA
INGENIERA EN ALIMENTOS
MTODOS NUMRICOS (TALLER)
GRUPO: 4AM1
TAREA 2
SERIE DE TAYLOR, MTODO DE BISECCIN, MTODO DEL PUNTO FIJO Y
MTODO DE NEWTON - RAPHSON
REYES RODRGUEZ ALEJANDRO
RUIZ HERNNDEZ BRENDA ANGLICA
PROFESOR:
RAMREZ BARRIOS MIGUEL LUIS
MXICO D.F. A 30 DE ABRIL DE 2015
Ejercicio 1
4
Sea f ( x)=x sin(2 x)
a) Con la ayuda de Matlab realice una grfica de f(x) desde x 0=0 hasta xf=1
CODIFICACIN EN MATLAB.
clc
clear all
close all
x=0:0.1:1
y=x.^4-sin(2*x)
plot(x,y)
grid on
title('Grfica de f(x)=x^4sin(2x) desde x0=0 hasta
xf=1')
xlabel('Eje X')
ylabel('Eje Y')
Figura 1. Grfica de la funcin en el intervalo de x0=0 hasta xf=1
.
b) Escriba analticamente una aproximacin polinomial de f(x) mediante la serie de Taylor de
4to Orden.
Teorema de Taylor: Si la funcin f y sus primeras n+1 derivada son continuas en un intervalo
que contiene a y x, entonces el valor de la funcin de x est dado por:
f ( x )=f ( a ) +f ' ( a ) ( xa )+
''
''
(n )
f (a)
f (a)
2 f ' (a)
3
n
( xa ) +
( x a ) ++
( x a ) + Rn
2!
3!
n!
Frecuentemente es conveniente simplificar la Serie de Taylor definiendo un paso
y expresado en la serie se tiene entonces:
''
f ( x i+1 ) =f ( x i ) + f ' ( x i ) h+
''
(n )
f ( xi ) 2 f ' ( xi ) 3
f ( xi ) n
h+
h ++
h + Rn
2!
3!
n!
Entonces, dada la funcin, se obtienen las derivadas hasta el orden 4.
f ' ( x )=4 x 3 2cos ( 2 x)
f ' ' ( x )=12 x 2+ 4 sin ( 2 x)
f ' ' ' ( x )=24 x +8 cos ( 2 x)
f IV ( x )=2416 sin ( 2 x)
h=xi +1x 1
c) Desarrolle un programa que calcule la serie de Taylor de orden 4 con h=0.1 de x 0=0 hasta
xf=1 y junto con ayuda del inciso a) llene la siguiente tabla.
x
0,2
0,5
0,6
0,7
Tabla1. Aproximacin de la serie de Taylor.
f(x)
Aproximacin de 4to Orden
-0,1986
1,8836
-0,6918
0,5605
-0,7790
0,2747
-0,8024
0,1062
d) Realice una grfica que compare f(x) con la aproximacin de 4to orden.
Figura 2. Grfica de las Aproximaciones de la funcin original.
CODIFICACIN:
clc
clear all
close all
h=0.1
j=1
y_aux=0
f_0_aux=0;
for y=0:h:1
f_real(j)=y_aux^4-sin(2*y_aux);
f_0(j)=y_aux^4-sin(2*y_aux);
f_1(j)=f_0(j)+(4*y_aux^3-2*cos(2*y_aux))*h;
f_2(j)=f_0(j)+(12*y_aux^2+4*sin(2*y_aux))*h;
f_3(j)=f_0(j)+(24*y_aux+8*cos(2*y_aux))*h;
f_4(j)=f_0(j)+(24-16*sin(2*y_aux))*h;
Ea_0(j)=(f_0(j)-f_0_aux)*100/f_0(j);
Ea_1(j)=(f_1(j)-f_0(j))*100/f_1(j);
Ea_2(j)=(f_2(j)-f_1(j))*100/f_2(j);
Ea_3(j)=(f_3(j)-f_2(j))*100/f_3(j);
Ea_4(j)=(f_4(j)-f_3(j))*100/f_4(j);
f_0_aux=f_0(j);
y_aux=y;
x(j)=y;
j=j+1;
end
figure(1)
plot(x,f_real,'b')
hold on
plot(x,f_1,'r')
plot(x,f_2,'g')
plot(x,f_3,'c')
plot(x,f_4,'m')
title('Grfica de f(x)junto con sus Aproximaciones')
xlabel('Eje X')
ylabel('Eje Y')
legend('Real','Primer Orden','Segundo Orden','Tercer Orden','Cuarto Orden')
grid on
Ejercicio 2
c
La velocidad de cada de un paracaista est dada por:
V ( t )=
t
gm
(1e m ) , conociendo los
c
siguientes valores:
t=20 s
g=9.81m/s 2
m=60 Kg
v =36 m/ s
Realice un programa para calcular el valor de la constante c con un error de aproximacin del
Ea=0.015 por medio del mtodo de Biseccin.
a) Describa el problema como un problema para encontrar races es decir, escriba f(c), poner
f en funcin de c.
Obtener el valor del coeficiente de resistencia (c).
t
gm
f ( c )=
1e m v=0
c
b) Elabore el programa, Cul es el valor de c?
La raz es
xr =
16.2769
c) Llene la siguiente tabla con los datos obtenidos en el inciso anterior, donde n es el nmero
de iteraciones, (xa:xb) el intervalo inicial del mtodo y x r es la raz por calcular que en este
caso es c.
Tabla 2. Obtencin de raz por el mtodo de biseccin.
n
0
1
2
3
4
5
6
xa
10
15
15
16,25
16,25
16,25
16,25
Xb
20
20
17,5
17,5
16,875
16,5625
16,40625
16,25
16,328125
xr
15
17,5
16,25
16,875
16,5625
16,40625
16,328125
16,289062
5
Ea
14,2857143
7,69230769
3,7037037
1,88679245
0,95238095
0,4784689
0,23980815
16,25
16,269531
3
16,269531
3
16,274414
1
9
10
11
16,289062
5
16,289062
5
16,279296
9
16,279296
9
16,269531
3
16,279296
9
16,274414
1
16,276855
5
0,12004802
0,059988
0,030003
0,01499925
d) Realice una grfica del nmero de iteraciones contra el valor absoluto del error
aproximado.
Nmero de Iteraciones VS Error Aproximado
16
14
12
10
Error Aproximado 8
6
4
2
0
0 2 4 6 8 10 12
Nmero de Iteraciones
Figura 3. Numero de iteraciones vs Error Aprox.
e) Realice una grfica del nmero de iteraciones contra xr
Nmero de Iteraciones VS Xr
18
17
16
Xr
15
14
13
0
Nmero de Iteraciones
Figura 4. Nmero de iteraciones vs Xr.
clc
clear all
close all
v=36;
m=60;
t=20;
g=9.81;
Es=0.015;
x_a=10;
x_b=20;
f_a=g*m*(x_a^-1)*(1-exp(-x_a*t/m))-v
f_b=g*m*(x_b^-1)*(1-exp(-x_b*t/m))-v
if f_a*f_b>=0
disp('La Raz no esta en el intervalo')
else
Ea=1;
CODIFICACIN
i=1;
while abs(Ea)>Es
x_r=(x_a+x_b)/2;
f_r=g*m*(x_r^-1)*(1-exp(-x_r*t/m))-v;
if f_a*f_r<0
x_b=x_r;% esta a la izquierda
end
if f_a*f_r>0
x_a=x_r;% esta a la derecha
end
if f_a*f_r==0
Ea=Es;
end
xr=(x_a+x_b)/2;
Ea=(xr-x_r)*100/xr
end
E_a(i)=Ea;
Es1(i)=Es;%auxiliar
n(i)=i;
i=1+i;
10
12
Ejercicio 3
2 x
Sea f ( x )=e
cos(0.192 x) . Mediante el Mtodo del Punto Fijo, calcule la raz de f en el
intervalo (2,15) con un error de aproximacin de Ea=0.016.
a) Escriba g(x).
xi+ 1=g( xi)
2 x
x=e
cos ( 0.192 x ) + x
xi+ 1=e xi cos ( 0.192 xi )+ x
xi+ 1=g( xi)
b) Cul es la raz de f en (2,15)?
La raz es
xr = 8.181
4.- Calcule la raz de la funcin del problema 3 por el Mtodo de Newton Raphson.
a) Calcule la derivada de f.
f ' ( x )=2e2 x +0.192sin ( 0.192 x)
b) Elabore un programa que calcule la raz de f Cul es la raz de f en (2,15)?
xr = 8.181
c) Llene la siguiente tabla donde xr es tal que f(xr)=0
Tabla 3. Obtencion de la raz por Newton-Raphson.
n
xr
Ea
1
8.1035
23.4033
2
8.1812
0.9501
3
8.1812
7,05902799122066e-05
d) Realice una grfica del nmero de iteraciones contra el valor absoluto del error
aproximado.
Figura 5. Nmero de Iteraciones vs Error Aprox.
e) Realice una grfica del nmero de iteraciones contra Xr.
Figura 6. Nmero de Iteraciones vs Xr.
f) Cul Mtodo, Punto Fijo o Newton- Raphson, converge en menos iteraciones?
Mtodo de Newton-Raphson
CONCLUSIONES
El mtodo de Newton-Raphson es el mtodo programado en las calculadoras y computadoras,
que el resultado es visto en menor nmero de iteraciones.
