Metos Numericos Ejercicios Matlab PARTE 2

LATEX
Ingenieria Mecatronica
UNIVERSIDAD DE LAS FUERZAS ARMADAS - ESPE

DEPARTAMENTO DE CIENCIAS EXACTAS
METODOS
NUMERICOS
DEBER 2
Nombre: Andr
es Enrquez
CURSO : H - 203
Profesor: Ing. Patricio Pugarn
Fecha: Lunes 18 de Mayo de 2015
SANGOLQUI - ECUADOR
1.
Ejercicios de Deber
1.1.
Ejercicio 1
En cada uno de los casos siguientes, halle el error absoluto Ez y el error relativo Rz y determine el numero de cifras
significativas de la aproximacion.
a) x = 2, 71828182, x
= 2, 7182
b) y = 98350, y = 98000
c) z = 0, 000068, z = 0, 00006
1. Solucion Literal a.
x = 2, 71828182, x
= 2, 7182
E x = x x
Ex = 2,71828182 2,7182
Ex = 8,182 105
Rx = xx
x
Rx =
8,182105
2,71828182
LATEX
Rx = 3,009989597 105
Cifras signiticativas:
Rx <
10d
2
Si d = 1
1
0,00003009989597 < 102
0,00006019979194 < 101
0,00006019979194 < 0,1
Si d = 2
0,00006019979194 < 102
0,00006019979194 < 0,01
Si d = 3
0,00006019979194 < 103
0,00006019979194 < 0,001
Si d = 4
0,00006019979194 < 104
0,00006019979194 < 0,0001
Si d = 5
0,00006019979194 < 105
0,00006019979194 < 0,00001
Con d=5 No satisface la desigualdad entonces tiene 4 cifras significativas.

2. Solucion Literal b.
y = 98350, y = 98000
Ey = y y
Ey = 98350 98000
Ey = 350
Ry = yy
y
350
Ry = 98350
Ry = 0,003558718861
Ry <
10d
2
Si d = 1
1
0,003558718861 < 102
0,007117437722 < 101
0,007117437722 < 0,1
LATEX
Si d = 2
0,007117437722 < 102
0,007117437722 < 0,01
Si d = 3
0,007117437722 < 103
0,007117437722 < 0,001
Para d = 3; no cumple entonces tiene 2 cifras significativas.

3. Solucion Literal c.
z = 0, 000068, z = 0, 00006
Ez = z z
Ez = 0,000068 0,00006
Ez = 0,000008
Rz = zz
z
Rz = 0,000008
0,000068
Rz = 0,1176470588
Rz <
10d
2
Si d = 1
1
0,1176470588 < 102
0,2352941176 < 101
0,2352941176 < 0,1
No cumple con ningun d entonces no tiene cifras significativas.
1.2.
Ejercicio 2
Usando Aritmetica de punto flotante con tres cifras significativas y redondeo, calcule las siguientes
sumas (sumando en el orden que se indica).
a)
b)
P6
1
k=1 3k
P6
1
k=1 37k
a) Solucion Literal a.
P6
1
k=1 3k
En forma fraccionaria.
LATEX
1
P6
= 3k =
P6
k=1
k=1
364
729
1
3
1
9
1
27
1
81
1
243
1
739
= 0,499
En forma decimal.
1
k
P6
=3
P6
= 0,498
k=1
k=1
= 0,333 + 0,111 + 0,037 + 0,012 + 0,004 + 0,001
b) Solucion literal b.
P6
1
k=1 37k
En forma fraccionaria.
P6
k=1
P6
=3
k=1
1
7k
364
729
1
729
1
243
1
81
1
27
1
9
1
3
= 0,499
En forma decimal.
1
7k
P6
=3
P6
= 0,498
k=1
k=1
1.3.
= 0,001 + 0,004 + 0,012 + 0,037 + 0,111 + 0,333
Ejercicio 3
Dados los desarrollos de Taylor:
1
= 1 + h + h2 + h3 + O(h4 )
1h
y
h2
h4
+
+ O(h6 )
2!
4!
Determine el orden de aproximacion de su suma y de su producto.
cos(h) = 1
Suma
1
+ cos(h) =
1h
h2
h4
+
+ O(h6 ))
2!
4!
h2
h2
([1 + h + h2 + h3 ] + O(h4 ) + ([1
+
] + O(h6 )))
2!
4!
1
h2
+ cos(h) = 2 + h +
+ h3 + O(h4 )
1h
2
(1 + h + h2 + h3 + O(h4 )) + (1
LATEX
Orden de aproximacion es O(h4 )

Producto
1
cos(h) =
1h
h2
h4
+
+ O(h6 ))
2!
4!
h2
h2
+
] + O(h6 )))
([1 + h + h2 + h3 ] + O(h4 ) + ([1
2!
4!
(1 + h + h2 + h3 + O(h4 )) (1
X = 1 + h + h2 + h3
4
2
Y = 1 h2! + h4!
1
cos(h) = X Y + X O(h6 ) + Y O(h4 ) + O(h4 ) O(h6 )
1h
h2
h3
1
cos(h) = 1 + h +
+
+ O(h4 )
1h
2!
2!
1.4.
Ejercicio 4
Encuentre una raiz positiva, con los metodos de la biseccion y Newton, a mano, de la ecuacion:
x2 4x sin(x) + (2 sin(x))2 = 0;
que se exacta hasta la segunda cifra significativa, solo con el empleo de una calculadora.
Metodo de Newton:
f (x) = x2 4xsen(x) + (2sen(x))2 ;
0
f (x) = 2x 4sen(x) 4xcos(x) + 4sen(x)cos(x);

f (x) = xn
f (xn )
f 0 (xn )
x1 = 2,34
x2 = 2,15
x3 = 2,04
x4 = 1,97
x5 = 1,94
x6 = 1,92
x7 = 1,92
La raiz es -1.92
LATEX
Metodo de la Biseccion:
f (x) = x2 4xsen(x) + (2sen(x))2 ;
a = 0, b = 2
x = 1,92
1.5.
Ejercicio 5
Con los programas: metodos de la biseccion, newton y metodo de la secante, encuentre una raiz de:
f (x) = x tan(x),
en el intervalo [1,2]. Compare y comente los resultados obtenidos.
disp(CALCULAR LA RAIZ)
x0=input(Ingrese el valor:
y=@(x) (x-tan(x));
a=diff(sym(x-tan(x)));
b=inline(a);
);
%derivada de la funcion
%Hacemos funcion a la derivada
x1=x0-(y(x0)/b(x0));
disp(x1)
while abs(x1-x0)>0.00001
x2=x1-(y(x1)/b(x1));
x0=x1;
x1=x2;
end
disp(x1)
1.6.
Ejercicio 6
Encontrar mediante el metodo de Newton las raices de la ecuacion..
3 sin(x) = cos(x)
disp(CALCULO DE RAICES -METODO NEWTON)
LATEX
x0=input(Ingrese el valor:
);
y=@(x) (3^(1/2)*sin(x)-cos(x));
der=diff(sym(3^(1/2)*sin(x)-cos(x))); %derivada de la funcion
%der
der1=inline(der);
%Hacemos funcion a la derivada
x1=x0-(y(x0)/der1(x0));
while abs(x1-x0)>0.00001
x2=x1-(y(x1)/der1(x1));
x0=x1;
x1=x2;
end
disp(x1)
1.7.
Ejercicio 7
Considerese la funcion f (x) = cosh(x) + cos(x) , con = 1, 2, 3. Encuentre un intervalo que contenga
un cero de f para cada valor de y calcule dicho cero el mEtodo de la biseccion.
Grafico de la funcion para = 1,2 y 3.
disp(Funcion: cosh(x)+cos(x)-g )
g=input(Elija la que numero desea para g(1,2,3) : );
a=input(Ingresar el limite a: );
b=input(Ingresar el limite b: );
fx=@(x) (cosh(x)+cos(x)-g);
dif=abs(b-a);
while dif>0.00001
if fx(a)*fx(b)<0
c=((0.5)*(a+b));
if (fx(a)*fx(c)<0)
b=c;
else
a=c;
end
else
disp(No hay raiz en la funcion)
return
end
dif=abs(b-a);
end
disp(La raiz de la funcion es: )
disp(c)
LATEX
Figura 1: Cuando = 1.
LATEX
1.8.
Ejercicio 8
Un objeto esta situado en un plano cuya pendiente varia a una tasa constante w. La posicion del objeto,
al instante t, esta dada por la formula.
s(t, w) =
g
[sinh(wt) sin(wt)],
2w2
donde g = 9, 8m/s2 es la aceleracion de la gravedad. Asumiendo que el objeto se ha desplazado 1 metro

en 1 segundo, calcule el valor de w, usando el metodo de la biseccion, con una tolerancia de 105 .
Cuantas iteraciones se requieren para alcanzar la tolerancia indicada?
function ejercicio27(a,b,t,f)
if f(a)*f(b) > 0
disp(No hay ra~
Az)
return
end
n=0;
while abs(a - b) > t
if f(a)*f(b)<0
c = (a+b)/2;
else
f(c)==0
x = c;
return
end
if f(a)*f(c)<0
b = c;
else
a = c;
end
n = n+1;
end
raiz = c
iteraciones = n
w = -2:0.0001:3;
y = w.^2 - 4.9.*sinh(w) + 4.9.*sin(w);
plot(w,y)
axis([-2 2 -2 2])
grid on
LATEX
1.9.
10
Ejercicio 9
Escriba e implemente un programa en Matlab para calcular la raiz cuadrada de un numero positivo a,
basado en el metodo de Newton.
function [ ] = ejercicio29(a)
%UNTITLED5 Summary of this function goes here
%
Detailed explanation goes here
if a<0
fprintf(La raiz cuadrada es imaginaria\n)
return
end
x0 = a+1;
f = @(x) x^2-a;
g = @(x) 2*x;
x(1) = x0;
x(2) = x(1)-f(x(1))/g(x(1));
n = 1;
while abs(x(n+1)-x(n))> sqrt(eps)
n = n+1;
x(n+1) = x(n)-f(x(n))/g(x(n));
end
resp1 = x(n+1);
resp2 = -x(n+1);
disp(La primera raiz es: )
disp(resp1)
disp(La segunda raiz des:)
disp(resp2)
end

Metos Numericos Ejercicios Matlab PARTE 2

Cargado por

Copyright:

Formatos disponibles

Metos Numericos Ejercicios Matlab PARTE 2

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

Metos Numericos Ejercicios Matlab PARTE 2

Cargado por

Copyright:

Formatos disponibles

LATEX

UNIVERSIDAD DE LAS FUERZAS ARMADAS - ESPE

Con d=5 No satisface la desigualdad entonces tiene 4 cifras significativas.

Para d = 3; no cumple entonces tiene 2 cifras significativas.

No cumple con ningun d entonces no tiene cifras significativas.

= 0,333 + 0,111 + 0,037 + 0,012 + 0,004 + 0,001

= 0,001 + 0,004 + 0,012 + 0,037 + 0,111 + 0,333

Dados los desarrollos de Taylor:

Orden de aproximacion es O(h4 )

f (x) = 2x 4sen(x) 4xcos(x) + 4sen(x)cos(x);

%Hacemos funcion a la derivada

Encontrar mediante el metodo de Newton las raices de la ecuacion..

disp(CALCULO DE RAICES -METODO NEWTON)

Grafico de la funcion para = 1,2 y 3.

donde g = 9, 8m/s2 es la aceleracion de la gravedad. Asumiendo que el objeto se ha desplazado 1 metro

También podría gustarte