UNIVERSIDAD NACIONAL DE SAN AGUSTIN DE

AREQUIPA
FACULTAD DE INGENIERIA DE PROCESOS
ESC. PROFESIONAL DE INGENIERIA DE INDUSTRIAS
ALIMENTARIAS

CURSO:
MÉTODOS NUMÉRICOS
DOCENTE:
PERCY TICONA CENTENO
ALUMNO:
FERNÁNDEZ QUEA JAHIR ARNALDO
CUI:
20150867

2018
AREQUIPA - PERÚ
 1
TRABAJO ENCARGADO

INTRODUCCIÓN

En ingeniería hay procesos que son modelados con ecuaciones

diferenciales ordinarias, cuya solución es imposible determinar por
métodos analíticos es allí la utilidad de los métodos numéricos que
calcula una solución aproximada por medio de un número finito de
iteraciones que mejora su eficiencia de manera rápida, al utilizar un
software adecuado.

Se define como un procedimiento mediante el cual se obtiene, casi

siempre de manera aproximada, la solución de ciertos problemas
realizando cálculos puramente aritméticos y lógicos como operaciones
aritméticas elementales, cálculo de funciones, consulta de una tabla de
valores, cálculo preposicional, etc. El procedimiento consiste de una lista
finita de instrucciones precisas que especifican una secuencia de
operaciones algebraicas y lógicas conocida como algoritmo, que
producen o bien una aproximación de la solución del problema o bien
un mensaje. La eficiencia en el cálculo de dicha aproximación depende,
en parte, de la facilidad de implementación del algoritmo y de las
características especiales y limitaciones de los instrumentos de cálculo.
En general, al emplear estos instrumentos de cálculo se introducen errores
llamados de redondeo.

METODOS NUMERICOS Mg. Percy Ticona Centeno

 2
TRABAJO ENCARGADO

MARCO TEÓRICO

Método De La Bisección

Supongamos que se tiene una función continua en el intervalo [a,b] y de

tal manera que f(a)*f(b) < 0 esto quiere decir que f(x) tiene un cero en
el intervalo abierto (a,b). Por la razón que el producto del valor de la
función en a y b es negativo es decir cambia de signo en el intervalo
[a,b], lo que afirma es una consecuencia del teorema del valor medio.

Pues el método en análisis explota el hecho anterior para su

fundamento, pues dicho método determina c = (a+b)/2 y averigua si
f(a) f(c) <0 si esto resulta siendo cierto entonces f(x) tiene una raíz en el
intervalo [a,c]. En seguida tomamos el valor de c como b y realizamos
el mismo análisis anterior.

Si ocurriera que f(a) f(c) >0 entonces f(c) f(b) < 0 en este caso redefinimos
a c =a En ambos caso a sucedido que se a determinado un nuevo
intervalo que contiene una raíz de la función y el proceso puede
repetirse.

Si f(a)* f(c) = 0, o f(c)* f(b) = 0; entonces f(c) = 0 y con esto se ha

determinado una raíz del polinomio, pero vale aclarar que este caso no
sucede en general puesto que los redondeos en una computadora difícil
es cero. Por esta razón es que para concluir se debe realizar con una
tolerancia de 10-3.

Este método también se le conoce con el nombre de método de la

bipartición, pero debemos destacar este método es el más sólido y
seguro que los otros métodos para encontrar una raíz en un intervalo.

INTERPRETACIÓN GRÁFICA

f(a)
f(b)

a c b

a c b

f(c) f(b) f(a)

caso a) f(a)*f(c)<0¸ Caso b) f*(c)f(b)<0

METODOS NUMERICOS Mg. Percy Ticona Centeno

 3
TRABAJO ENCARGADO

Matlab: Software Matemático

MATLAB es un entorno de cálculo técnico de altas prestaciones para

cálculo numérico y visualización. Integra:

 Análisis numérico
 Cálculo matricial
 Procesamiento de señales
 Gráficos

en un entorno fácil de usar, donde los

problemas y las soluciones son
expresados como se escriben
matemáticamente, sin la programación
tradicional. El nombre MATLAB proviene
de ``MATrix LABoratory'' (Laboratorio de
Matrices). MATLAB fue escrito
originalmente para proporcionar un
acceso sencillo al software matricial
desarrollado por los
proyectos LINPACK y EISPACK, que juntos representan lo más avanzado
en programas de cálculo matricial. MATLAB es un sistema interactivo
cuyo elemento básico de datos es una matriz que no requiere
dimensionamiento. Esto permite resolver muchos problemas numéricos en
una fracción del tiempo que llevaría hacerlo en lenguajes
como C, BASIC o FORTRAN. MATLAB ha evolucionado en los últimos años
a partir de la colaboración de muchos usuarios. En entornos universitarios
se ha convertido en la herramienta de enseñanza estándar para cursos
de introducción en álgebra lineal aplicada, así como cursos avanzados
en otras áreas.

En la industria, MATLAB se utiliza para investigación y para resolver

problemas prácticos de ingeniería y matemáticas, con un gran énfasis en
aplicaciones de control y procesamiento de señales. MATLAB también
proporciona una serie de soluciones específicas
denominadas TOOLBOXES. Estas son muy importantes para la mayoría de
los usuarios de MATLAB y son conjuntos de funciones MATLAB que
extienden el entorno MATLAB para resolver clases particulares de
problemas como:

 Procesamiento de señales
 Diseño de sistemas de control
 Simulación de sistemas dinámicos
 Identificación de sistemas
 Redes neuronales y otros.

METODOS NUMERICOS Mg. Percy Ticona Centeno

 4
TRABAJO ENCARGADO

Probablemente la característica más importante de MATLAB es su

capacidad de crecimiento. Esto permite convertir al usuario en un autor
contribuyente, creando sus propias aplicaciones. En resumen, las
prestaciones más importantes de MATLAB son:

 Escritura del programa en lenguaje matemático.

 Implementación de las matrices como elemento básico del
lenguaje, lo que permite una gran reducción del código, al no
necesitar implementar el cálculo matricial.
 Implementación de aritmética compleja.
 Un gran contenido de órdenes específicas, agrupadas
en TOOLBOXES.
 Posibilidad de ampliar y adaptar el lenguaje, mediantes ficheros
de script y funciones

METODOS NUMERICOS Mg. Percy Ticona Centeno

 5
TRABAJO ENCARGADO

EJERCICIOS PROPUESTOS

Resolución
1.- f (a) f (b) ≤ 0

Se plantea en esta ecuación donde se presencia el caso donde NO

EXISTE UNA RAIZ ENTRE a y b por tanto f (a) y f (b) tienen el mismo SIGNO.
2.- f (a) f (b) ≥ 0

Se plantea en esta ecuación donde se presencia el caso donde EXISTE

AL MENOS UNA RAIZ ENTRE a y b por tanto f(a) y f(b) tienen un SIGNO
completamente DIFERENTE.

1.- En la ventana del programa Matlab se usaron los siguientes

COMANDOS para poder realizar el grafico correspondiente:
>> x=[-2:0.01:1];
>> y=1000*sin(x.^3+1)./log10(5+x.^2);
>> plot(x,y)
Dando como resultado de la aplicación de estos la siguiente gráfica:

Interpretación:

En la gráfica se observa que

existe dos raíces entre -2 y -
1.5 asi como también se
observa que existe una raíz
entre -1.75 y -1.

METODOS NUMERICOS Mg. Percy Ticona Centeno

 6
TRABAJO ENCARGADO

2.- En la ventana del programa Matlab se usaron los siguientes

COMANDOS para poder realizar el grafico correspondiente:
>> x=[-10:0.01:10];
>> y=log(x.^2+1)-200*(x+10).^3+9*x.^2+5;
>> plot(x,y)
Dando como resultado de la aplicación de estos la siguiente gráfica:

Interpretación:

En La Gráfica Se Observa Que

Existe Una Raíz Entre -10 Y -8.

1.- En la ventana del programa Matlab se usaron los siguientes

COMANDOS para poder realizar el grafico correspondiente:
>> x=[0:0.01:20];
>> y=x-5*log(x);
>> plot(x,y)
Dando como resultado de la aplicación de estos la siguiente gráfica:

Interpretación:

En la gráfica se observa que

existe una raiz entre 1 y 2 asi
como también se observa que
existe una raiz entre 12 y 13

METODOS NUMERICOS Mg. Percy Ticona Centeno

 7
TRABAJO ENCARGADO

COMANDOS APLICADOS EN EL METODO DE LA BISECCION:

function biseccion()

xai=input('Ingrese LIMITE INFERIOR ');

xbi=input('Ingrese LIMITE SUPERIOR ');
tol=input(' error: ');
syms x;
f=input(' funciòn: ');
i=1;

f1=subs(f,x,xai);
f2=subs(f,x,xbi);
%GRAFICA

error(i)=1;
if f1*f2 < 0
xa(i)=xai; f1=subs(f,x,xa(i));

xb(i)=xbi; f2=subs(f,x,xb(i));

xr(i)=(xa(i)+xb(i))/2;
f3=subs(f,x,xr(i));

fprintf('It. Xa Xr Xb
Error \n');
fprintf('%2d \t %11.7f \t %11.7f \t %11.7f
\n',i,xa(i),xr(i),xb(i));
while abs(error(i)) >= tol,
if f1*f3<0
xa(i+1)=xa(i);f1=subs(f,x,xa(i+1));
xb(i+1)=xr(i);f2=subs(f,x,xb(i+1));
end
if f1*f3> 0
xa(i+1)=xr(i);f1=subs(f,x,xa(i+1));
xb(i+1)=xb(i);f2=subs(f,x,xb(i+1));
end
xr(i+1)=(xa(i+1)+xb(i+1))/2; f3=subs(f,x,xr(i+1));
error(i+1)=abs((xr(i+1)-xr(i))/(xr(i+1)));
fprintf('%2d \t %11.7f \t %11.7f \t %11.7f \t %7.3f
\n',i+1,xa(i+1),xr(i+1),xb(i+1),error(i+1));
i=i+1;
end

end

METODOS NUMERICOS Mg. Percy Ticona Centeno

 8
TRABAJO ENCARGADO

PARA LA FUNCION: y=1000*sin(x.^3+1)./log10(5+x.^2);

>> biseccion()

Ingrese LIMITE INFERIOR -2

Ingrese LIMITE SUPERIOR -1.75

error: 0.000000001

funciòn: 1000*sin(x^3+1)/log10(5+x^2)

It. Xa Xr Xb Error

1 -2.0000000 -1.8750000 -1.7500000

2 -2.0000000 -1.9375000 -1.8750000 0.032

3 -2.0000000 -1.9687500 -1.9375000 0.016

4 -1.9687500 -1.9531250 -1.9375000 0.008

5 -1.9531250 -1.9453125 -1.9375000 0.004

6 -1.9453125 -1.9414063 -1.9375000 0.002

7 -1.9414063 -1.9394531 -1.9375000 0.001

8 -1.9394531 -1.9384766 -1.9375000 0.001

9 -1.9384766 -1.9379883 -1.9375000 0.000

10 -1.9384766 -1.9382324 -1.9379883 0.000

11 -1.9384766 -1.9383545 -1.9382324 0.000

12 -1.9384766 -1.9384155 -1.9383545 0.000

13 -1.9384155 -1.9383850 -1.9383545 0.000

14 -1.9384155 -1.9384003 -1.9383850 0.000

15 -1.9384003 -1.9383926 -1.9383850 0.000

16 -1.9383926 -1.9383888 -1.9383850 0.000

17 -1.9383888 -1.9383869 -1.9383850 0.000

18 -1.9383869 -1.9383860 -1.9383850 0.000

19 -1.9383869 -1.9383864 -1.9383860 0.000

20 -1.9383869 -1.9383867 -1.9383864 0.000

21 -1.9383869 -1.9383868 -1.9383867 0.000

22 -1.9383869 -1.9383869 -1.9383868 0.000

23 -1.9383869 -1.9383868 -1.9383868 0.000

METODOS NUMERICOS Mg. Percy Ticona Centeno

 9
TRABAJO ENCARGADO

24 -1.9383869 -1.9383868 -1.9383868 0.000

25 -1.9383869 -1.9383869 -1.9383868 0.000

26 -1.9383869 -1.9383868 -1.9383868 0.000

27 -1.9383868 -1.9383868 -1.9383868 0.000

LA RAIZ ES -1.9383868

PARA LA FUNCION: log(x.^2+1)-200*(x+10).^3+9*x.^2+5

>> biseccion()

Ingrese LIMITE INFERIOR -10

Ingrese LIMITE SUPERIOR -8

error: 0.000000001

funciòn: log(x^2+1)-200*(x+10)^3+9*x^2+5

It. Xa Xr Xb Error

1 -10.0000000 -9.0000000 -8.0000000

2 -9.0000000 -8.5000000 -8.0000000 0.059

3 -9.0000000 -8.7500000 -8.5000000 0.029

4 -8.7500000 -8.6250000 -8.5000000 0.014

5 -8.6250000 -8.5625000 -8.5000000 0.007

6 -8.5625000 -8.5312500 -8.5000000 0.004

7 -8.5312500 -8.5156250 -8.5000000 0.002

8 -8.5156250 -8.5078125 -8.5000000 0.001

9 -8.5156250 -8.5117188 -8.5078125 0.000

10 -8.5117188 -8.5097656 -8.5078125 0.000

11 -8.5117188 -8.5107422 -8.5097656 0.000

12 -8.5107422 -8.5102539 -8.5097656 0.000

13 -8.5107422 -8.5104980 -8.5102539 0.000

14 -8.5104980 -8.5103760 -8.5102539 0.000

15 -8.5103760 -8.5103149 -8.5102539 0.000

16 -8.5103760 -8.5103455 -8.5103149 0.000

17 -8.5103455 -8.5103302 -8.5103149 0.000

METODOS NUMERICOS Mg. Percy Ticona Centeno

 10
TRABAJO ENCARGADO

18 -8.5103455 -8.5103378 -8.5103302 0.000

19 -8.5103455 -8.5103416 -8.5103378 0.000

20 -8.5103455 -8.5103436 -8.5103416 0.000

21 -8.5103455 -8.5103445 -8.5103436 0.000

22 -8.5103455 -8.5103450 -8.5103445 0.000

23 -8.5103455 -8.5103452 -8.5103450 0.000

24 -8.5103452 -8.5103451 -8.5103450 0.000

25 -8.5103452 -8.5103452 -8.5103451 0.000

26 -8.5103452 -8.5103452 -8.5103452 0.000

27 -8.5103452 -8.5103452 -8.5103452 0.000

28 -8.5103452 -8.5103452 -8.5103452 0.000

LA RAIZ ES: -8.5103452

PARA LA FUNCION: X-5LN(X)

>> biseccion()

Ingrese LIMITE INFERIOR 1

Ingrese LIMITE SUPERIOR 2

error: 0.000000001

funciòn: x-5*log(x)

It. Xa Xr Xb Error

1 1.0000000 1.5000000 2.0000000

2 1.0000000 1.2500000 1.5000000 0.200

3 1.2500000 1.3750000 1.5000000 0.091

4 1.2500000 1.3125000 1.3750000 0.048

5 1.2500000 1.2812500 1.3125000 0.024

6 1.2812500 1.2968750 1.3125000 0.012

7 1.2812500 1.2890625 1.2968750 0.006

8 1.2890625 1.2929688 1.2968750 0.003

9 1.2929688 1.2949219 1.2968750 0.002

METODOS NUMERICOS Mg. Percy Ticona Centeno

 11
TRABAJO ENCARGADO

10 1.2949219 1.2958984 1.2968750 0.001

11 1.2949219 1.2954102 1.2958984 0.000

12 1.2954102 1.2956543 1.2958984 0.000

13 1.2956543 1.2957764 1.2958984 0.000

14 1.2957764 1.2958374 1.2958984 0.000

15 1.2958374 1.2958679 1.2958984 0.000

16 1.2958374 1.2958527 1.2958679 0.000

17 1.2958527 1.2958603 1.2958679 0.000

18 1.2958527 1.2958565 1.2958603 0.000

19 1.2958527 1.2958546 1.2958565 0.000

20 1.2958546 1.2958555 1.2958565 0.000

21 1.2958546 1.2958550 1.2958555 0.000

22 1.2958550 1.2958553 1.2958555 0.000

23 1.2958553 1.2958554 1.2958555 0.000

24 1.2958554 1.2958555 1.2958555 0.000

25 1.2958555 1.2958555 1.2958555 0.000

26 1.2958555 1.2958555 1.2958555 0.000

27 1.2958555 1.2958555 1.2958555 0.000

28 1.2958555 1.2958555 1.2958555 0.000

29 1.2958555 1.2958555 1.2958555 0.000

30 1.2958555 1.2958555 1.2958555 0.000

LA RAIZ: 1.2958555

METODOS NUMERICOS Mg. Percy Ticona Centeno

 12
TRABAJO ENCARGADO

BIBLIOGRAFIA

 “TEXTO: METODOS NUMERICOS PARA ECUACIONES DIFERENCIALES

ORDINARIAS CON MATLAB Recuperado de:”
https://unac.edu.pe/images/documentos/organizacion/vri/cdcitr
a/Informes_Finales_Investigacion/2011/Junio/IF_COLLANTE_HUANT
O_FIME.pdf
 Unidad I: Introducción A Los Métodos Numéricos. Recuperado de:
https://sites.google.com/site/mecanalisis/home/unidad-i-
introduccion-a-los-metodos-numericos

METODOS NUMERICOS Mg. Percy Ticona Centeno

 13
TRABAJO ENCARGADO

ANEXOS
 INTRODUCCION AL PROGRAMA MATLAB

METODOS NUMERICOS Mg. Percy Ticona Centeno

 14
TRABAJO ENCARGADO

METODOS NUMERICOS Mg. Percy Ticona Centeno

 15
TRABAJO ENCARGADO

METODOS NUMERICOS Mg. Percy Ticona Centeno