UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
UNIVERSIDAD NACIONAL DE INGENIERÍA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
Escuela Profesional de Ingeniería Química
INFORME FINAL
REGLA FALSA: UNA FAMILIA DE MÉTODOS DE BÚSQUEDA DE UNA RAÍZ
(A FAMILY OF REGULA FALSI ROOT FINDING METHODS)
Integrantes:
Jiménez Ferro, Rita ………………………………. 20172181H
Méjico Gutiérrez, Raúl Enrique……………….. ….20171505D
Toro Isla, Iván ………………………………………20170370H
Navarro La Rosa, Juan Eduardo …………….......20082175I
Docentes:
Dávila
LIMA-PERU (2019-11-30)
2019-2
1
UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
I. RESUMEN
En la ingeniería química se presentan muchos problemas donde para su
solución es necesario utilizar ecuaciones no lineales. Este tipo de problemas
representa en la actualidad uno de los temas básicos y necesarios dentro del
análisis numérico.
El presente informe, está referido a la aplicación de ecuaciones no lineales en
problemas aplicados en ingeniería química haciendo uso de una técnica
numérica, al mismo tiempo se han resuelto problemas con diversos softwares
matemáticos. Introduciendo conceptos de química relacionados a los
problemas a tratar, se estudia a profundidad una familia de métodos de
búsqueda de una raíz: el método de regula falsi, también llamado método de
la falsa posición, en la solución de ecuaciones no lineales justificando su
eficiencia mediante su orden de convergencia o, en algunos casos, sobre la
base de resultados numéricos de tal manera que se define con precisión los
resultados obtenidos.
Por tanto, el trabajo comienza con una descripción del método de falsa
posición para resolver ecuaciones no lineales. Se discuten sus ventajas y
desventajas, y se refieren algunas estrategias para acelerar su convergencia.
Además de ello se detallan diversos problemas de ingeniería química
donde intervienen directamente ecuaciones no lineales y que para hallar su
solución es necesario recurrir algún método numérico que detalle en forma
simple y precisa la solución del mismo.
Finalmente, estos resultados son asistidos y comprobados mediante diversos
softwares matemáticos utilizados en el presente informe.
2
UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
II. INTRODUCCIÓN
En muchas áreas de las matemáticas y la ingeniería existen problemas de
modelado y simulación de procesos que implican en algún momento resolver
una ecuación no lineal.
El problema de la resolución de ecuaciones no lineales figura entre los más
importantes en la teoría y la práctica, no solo de las matemáticas aplicadas,
sino también de muchas ramas de las ciencias, la ingeniería, la física, la
informática, la astronomía, las finanzas, etc.
En una situación práctica, un problema matemático no es más que la
modelización de un fenómeno, en nuestro caso químico. Este modelo
matemático es la imagen matemática del fenómeno de estudio y se
representa mediante una gran variedad de ecuaciones: transcendentes, en
derivadas ordinarias o en derivadas parciales, sujetos a condiciones iniciales
o de frontera (o combinación de ambas), ecuaciones integrales, todo ellos en
dimensión finita o infinita, en su mayoría no lineales.
En los casos en los que no somos capaces de encontrar la solución analítica
o esta es muy costosa, optamos por obtener una aproximación de dicha
solución mediante procesos numéricos que aproximan de manera eficiente y
con suficiente precisión, las soluciones de problemas expresados
matemáticamente. El modelo matemático implementado en algoritmos lógico-
numéricos en un ordenador, permite estudiar las cualidades del proceso
original uniendo las ventajas de la teoría y del experimento.
3
UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
Al trabajar con un modelo matemático y no con el fenómeno de estudio, se
pueden estudiar y pronosticar sus propiedades de estado (ventaja teórica) en
forma relativamente rápida y a bajo costo. Al mismo tiempo los algoritmos
numéricos permiten, apoyándose en la potencia de calculo de los
ordenadores, verificar las cualidades del fenómeno de estudio en una forma
no accesible para los enfoques teóricos (ventaja del experimento).
En el campo de ecuaciones lineales se ha llevado a cabo una intensa labor
investigadora, siendo considerable el volumen de publicaciones sobre este
tema que se han generado en los últimos años. Sin embargo, el número de
publicaciones para la solución de ecuaciones no lineales es menor y mas aun
si esta relacionado a su aplicación en la Ingeniería Química.
La metodología empleada en el presente trabajo consiste básicamente en
detallar la teoría referente a la parte matemática y química, para después
hacer la resolución practica en algunos ejemplos y problemas planteados
para facilitar la comprensión de lo estudiado, de tal forma que se desarrolle
en forma clara y precisa la aplicación de ecuaciones no lineales en ingeniería
química. Además, se utilizan diversos softwares matemáticos para hallar sus
resultados en forma rápida y precisa.
4
UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
III. MARCO TEORICO
3.1. ANÁLISIS NUMÉRICO
El análisis numérico tiene que ver con el desarrollo y evaluación de
métodos para calcular los resultados numéricos requeridos a partir de datos
numéricos (Scheid Francis, 1972).
Los análisis numéricos han sido desarrollados con el objeto de resolver
problemas matemáticos cuya solución es difícil o en ocasiones imposible de
resolver por medio de Jos análisis tradicionales.
Las soluciones que ofrecen los análisis numéricos son aproximaciones
de los valores reales y, por tanto, se tendrá un cierto grado de error que será
conveniente determinar.
Son herramientas muy poderosas utilizadas para la solución de
problemas. Tienen la capacidad de manejar sistemas de ecuaciones grandes,
no linealidades y geometrías complicadas las cuales son bastante comunes
en la práctica de la ingeniería y que, en ocasiones, son imposibles de
resolverse de una manera analítica. Por lo tanto, los análisis numéricos
amplían la habilidad de quien los estudia para la solución de problemas.
La gran mayoría de los análisis numéricos son procesos cíclicos o
iterativos, en los cuales se repite una serie de pasos y estos se basan en las
denominadas ecuaciones o fórmulas de recurrencia, las cuales relacionan
dos o más elementos consecutivos de una sucesión de números, funciones,
matrices, etc.
5
UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
3.2. ECUACIONES NO LINEALES
Algunos conceptos generales sobre las ecuaciones no lineales de una sola
variable independiente, así como algunos métodos numéricos para
resolverlas se mostrarán a continuación.
3.2.1 Generalidades
El problema de resolver una ecuación no lineal de la forma:
f(x)=0 (3.2.1)
consiste en encontrar el o los valores x que la anulan. A éstos se les
denomina raíces, ceros o soluciones de la función f(x).
Tomando en cuenta que en general las raíces de las ecuaciones
pueden ser reales o complejas, los problemas de búsqueda de raíces
pueden dividirse en dos (Chapra & Canale, 1999) :
Determinación de las raíces reales de ecuaciones algebraicas y
trascendentes. Los métodos numéricos utilizados para resolver
este tipo de problemas sólo permiten determinar el valor de una
sola raíz.
Determinación de todas las raíces reales y complejas de
polinomios. Los métodos numéricos en este rubro determinan
sistemáticamente todas las raíces del polinomio.
Por otra parte, los métodos de resolución de las ecuaciones no lineales
pueden clasificarse en los siguientes rubros:
6
UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
Método gráfico. Consiste en graficar la función y observar dónde
cruza el eje de las x. Este valor es donde f(x)=0, por lo tanto, es
la raíz. Este método es poco preciso, pero puede utilizarse
como punto de partida para otros métodos más sofisticados.
Método de prueba y error. Consiste en elegir un valor de x y
evaluar si f(x) es cero. Si no es así se elige otro valor y se
evalúa nuevamente. El proceso se repite hasta que se obtenga
un valor que proporcione una f(x) igual o cercano a cero bajo
determinada tolerancia prestablecida.
Métodos iterativos: Con base en un algoritmo específico y un
valor inicial supuesto para la raíz, 𝑥0 ,se construye una sucesión
de números reales x=bbbbb convergente a la solución
de la ecuación. Si el método es adecuado se cumplirá que:
en cuyo caso se dice que el método converge a la solución; de
lo contrario el método diverge.
Métodos directos: Proporcionan la solución mediante un número
finito de operaciones elementales. Por ejemplo, la conocida
ecuación para calcular las raíces de un polinomio de segundo
grado:
Ccccccc
En la siguiente sección se describen el método iterativo de Regla falsa
para resolver ecuaciones no lineales de una sola variable
independiente.
7
UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
3.2.2 Métodos iterativos
Los métodos iterativos consisten en aplicar una regla
determinada para generar una sucesión de valores cada vez más
próximos a alguna de las raíces de la ecuación no lineal que se
trata de resolver. Dependiendo de las características de la
ecuación, de la aproximación inicial y de la regla utilizada, los
métodos iterativos pueden o no converger a la solución.
Métodos cerrados y métodos abiertos
Los métodos para resolver una ecuación no lineal se denominan
cerrados si requieren un conocimiento previo del intervalo que
contiene la raíz. Si esto no es necesario, se dice que el método
es abierto, y basta una aproximación inicial para comenzar el
proceso iterativo.
Para identificar el intervalo que contiene la raíz en los métodos
cerrados o establecer la aproximación inicial de un método
abierto, debe tenerse un buen entendimiento del problema físico
que se está modelando con la ecuación y adquirir la mayor
cantidad de información posible respecto al mismo; la utilidad de
graficar la función de interés resulta más que evidente.
8
UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
En los casos más complicados, como por ejemplo cuando la
función a resolver se define durante la ejecución de un programa
computacional, una opción es incorporar una búsqueda
incremental. Ésta consiste en realizar evaluaciones de la función
en dos puntos. Si la función cambia de signo, entonces al menos
existe una raíz en ese intervalo. Si esto no es así, se amplía el
intervalo de búsqueda incrementando el valor de uno de los
puntos y realizando la evaluación correspondiente de la función;
este proceso termina hasta identificar el cambio de signo referido,
o al sobrepasar un intervalo máximo prestablecido, en cuyo caso
no ha sido posible identificar la existencia de la raíz en el mismo.
Un problema con la búsqueda incremental radica en la elección
de la longitud del incremento. De esta manera, la búsqueda será
muy tardada si la longitud es muy pequeña. Por otro lado, si la
longitud es muy grande, entonces las raíces cercanas entre sí
pueden pasar inadvertidas. En este último caso, si existen raíces
múltiples, se puede calcular la primera derivada de la función
f(x) al inicio y al final de cada intervalo. Cuando la derivada
cambia de signo, puede existir un máximo o un mínimo en ese
intervalo, lo que representa una búsqueda más detallada para
detectar la posibilidad de una raíz.
3.2.3 Método de la falsa posición
9
UNIVERSIDAD NACIONAL DE INGENÍERIA
FACULTAD DE INGENIERÍA QUÍMICA Y TEXTIL
DEPARTAMENTO ACADÉMICO DE CIENCIAS BÁSICAS
3.2.4
10