UNIVERSIDAD CATOLICA DE SANTA MARIA Página:1/7
FACULTAD DE CIENCIAS E INGENIERÍAS FISICAS Y FORMALES
ESCUELA PROFESIONAL Jefes de Prácticas:
INGENIERIA MECÁNICA, MECÁNICA-ELECTRÍCA Y MECATRÓNICA Ing. Juan Carlos Cuadros
Ing. Fernando Siles
Laboratorio de Métodos Numéricos para Ingenieros
Código: 4K05107
Tema: Solución de Ecuaciones No Lineales f(x)=0
Semestre: V
Método de la Bisección Grupo: FECHA:
Apellidos y Nombres (1):____________________________________________
Lab. Nº: 03 08/abr/2025
Apellidos y Nombres (2):____________________________________________
I. OBJETIVOS
• Entender y aplicar los métodos iterativos de intervalo, específicamente el método de la bisección en la solución de
ecuaciones no lineales algebraicas y trascendentes.
• Comprender y aplicar los métodos iterativos abiertos, específicamente el método de Newton-Raphson en la solución de
ecuaciones no lineales, algebraicas y trascendentes.
• Implementar programas en MATLAB que solucionen las raíces de ecuaciones algebraicas y trascendentes por medio del
método de la bisección.
• Implementar programas en MATLAB que solucionen las raíces de ecuaciones algebraicas y trascendentes por medio del
método de Newton-Raphson.
II. MARCO TEORICO
2.1 Método de la Bisección.
En esta práctica, estudiaremos dos de los métodos para la solución numérica de ecuaciones algébricas y trascendentes
no lineales, esto es, ecuaciones que se puedan escribir en la forma f(x) = 0, donde f es una función real de variable real.
Todos los valores s que anulan f, esto es, tales que f(s) = 0, se llaman ceros o raíces de la función f o solución de la
ecuación f(x)=0.
Figura 1 à Ceros de una función- Diferentes casos
Para una ecuación del tipo f(x) = 0, antes de intentar aplicar cualquier método de resolución, es importante garantizar
que la ecuación tenga solución, o sea, que existe un real; s tal que f(s)=0. Muchas veces importa también determinar si
la solución es única, o si existen diferentes soluciones y, en este caso, saber cual o cuales importan determinar.
2.2 Teorema de Bolzano
Dado un intervalo cerrado [a,b] y una función continua f(x), existe al menos una solución a la ecuación algebraica no
lineal o trascendente f(x)=0 si f(a)*f(b) 0
Figura 2 → Aplicación del Teorema de Bolzano
Laboratorio de Métodos Numéricos para Ingenieros Página: 2/7
Tema: Solución de Ecuaciones No Lineales f(x)=0
LAB N° 3 JCC - FSN
Método de la Bisección
El teorema de Bolzano garantiza la existencia de una raíz si existe un cambio de signo en el intervalo [a,b], pero la antítesis
es falsa, sino existe un cambio de signo, también puede existir una raíz en [a,b].
Los métodos numéricos que están basados en el teorema de Bolzano se denominan “métodos cerrados”, ya que exigen
como argumento un intervalo cerrado donde la función experimente un cambio de signo.
2.3. Método de las Bisecciones Sucesivas
En la resolución de ecuaciones no lineales se utilizan, salvo soluciones analíticas simples, métodos iterativos que generan
una sucesión de valores que tienden al valor de la raíz. Este método presenta la ventaja de acotar no sólo el valor de la
función, sino también el intervalo al que pertenece la raíz. Para su aplicación es necesario que verifique las condiciones
del Teorema de Bolzano, esto es, la función debe ser continua y cambiar de signo en sus extremos. Por ejemplo: Resolver
en los números reales las ecuaciones:
x 3 − 1 = 0, sen( x) − x = 0 ó tan( x) − x = 0 (*)
Las gráficas siguientes ilustran mejor la situación. Indican que un problema con enunciado tan simple puede ser difícil de
resolver:
Figura 3 à Gráficas de las ecuaciones (*)
La primera gráfica sugiere que hay un único cero, precisamente donde ese cero está y las otras dos indican que hay alguno
en cada caso, pero nada más. En realidad, sen(x) y x se encuentran en un único punto, x = 0; pero en cambio tan(x) y x se
encuentran en un número infinito de puntos. Ciertamente se requiere disponer de herramientas teóricas (teoremas) y
prácticas (algoritmos) para resolver problemas de esta clase.
El método de las bisecciones sucesivas parte del intervalo inicial [a, b] que se sabe contiene un cero de f, supuestamente
único. En cada iteración se produce la reducción del intervalo a la mitad del intervalo actual. Para lo cual, se divide el
intervalo actual escogiéndose el intervalo izquierdo o derecho de forma que la función tenga signo diferente en los
extremos del sub-intervalo escogido. O sea,
Laboratorio de Métodos Numéricos para Ingenieros Página: 3/7
Tema: Solución de Ecuaciones No Lineales f(x)=0
LAB N° 3 JCC - FSN
Método de la Bisección
Figura 4 → Bisecciones Sucesivas
1
siendo [an, bn] el intervalo de la iteración n, se calcula 𝑥𝑛+1 = (𝑎𝑛 + 𝑏𝑛 ). El valor 𝑥𝑛+1 sustituye a an o bn de acuerdo si
2
se cumple que 𝑓(𝑥𝑛+1) 𝑓(𝑏𝑛 ) < 0 o 𝑓(𝑥𝑛+1) 𝑓(𝑎𝑛 ) < 0. De esta forma, se asegura que s [an,bn] en cualquier iteración.
2.5. Algoritmo de Bisección
Método de las bisecciones sucesivas
Inicialización a0 , b0 = a, b
Repetir a n + bn
1. xn+1 =
2
2. Si f ( x n +1 ) f (a n ) 0
Entonces a n +1 = a n ; bn +1 = x n +1 ;
Sino a n +1 = x n +1 ; bn +1 = bn ;
Hasta que Verificar criterio de parada
2.5. Error del Algoritmo de Bisección
Figura 5 à Error del Método de la Bisección
Teorema 1: Sea f una función continua en el intervalo [a,b] tal que 𝑓(𝑎)𝑓(𝑏) < 0. El método de la bisección genera una
sucesión {pn} que converge a p con la propiedad:
𝑏−𝑎
|𝑝 − 𝑝𝑛 | ≤ 𝑛≥1
2𝑛
Demostración:
Para cada 𝑛 ≥ 1, tenemos:
1
𝑏𝑛 − 𝑎𝑛 = (𝑏 − 𝑎) y 𝑝 ∈ (𝑎𝑛 , 𝑏𝑛 )
2𝑛−1
Laboratorio de Métodos Numéricos para Ingenieros Página: 4/7
Tema: Solución de Ecuaciones No Lineales f(x)=0
LAB N° 3 JCC - FSN
Método de la Bisección
Puesto que 𝑝𝑛 = (𝑎𝑛 + 𝑏𝑛 )/2, para cada 𝑛 ≥ 1, se tiene:
1 𝑏−𝑎
|𝑝 − 𝑝𝑛 | ≤ (𝑏𝑛 − 𝑎𝑛 ) = 𝑛
2 2
y esta desigualdad implica que lim 𝑝𝑛 = 𝑝
𝑛→∞
Si llamamos δ al valor absoluto del error, donde 𝛿 = |𝑝 − 𝑝𝑛 |, se puede entonces afirmar que el error absoluto del
estimado pn está acotado por:
𝑏−𝑎
2𝑛
𝑏−𝑎
El número de iteraciones suficientes para garantizar un error absoluto no superior a se puede calcular haciendo ≤
2𝑛
𝑏−𝑎
𝛿 obteniéndose el valor: 𝑛 ≥ 𝑙𝑜𝑔2 ( 𝛿 )
Otra sugerencia sería finalizar el cálculo cuando el error verdadero se encuentra por debajo de algún nivel prefijado.
Por lo tanto se requiere estimar el error de forma tal que no se necesite el conocimiento previo de la raíz de la ecuación.
Se puede calcular el error relativo porcentual 𝜀𝑎 de la siguiente manera:
𝑥𝑟𝑛𝑢𝑒𝑣𝑜 − 𝑥𝑟𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟
𝜀𝑎 = | | 100%
𝑥𝑟𝑛𝑢𝑒𝑣𝑜
donde 𝑥𝑟𝑛𝑢𝑒𝑣𝑜 es la raíz en la iteración actual y 𝑥𝑟𝑎𝑛𝑡𝑒𝑟𝑖𝑜𝑟 es el valor de la raíz en la iteración anterior. Se utiliza el valor
absoluto, ya que solo importa la magnitud del error. Cuando 𝜀𝑎 es menor que un valor previamente fijado 𝜀𝑠 , termina el
cálculo.
El siguiente ejemplo ilustra la aplicación de este método.
Ejemplo 1:
Determinar una aproximación con un error absoluto inferior a 7x10-2 de la (única) solución de la ecuación 1 + 𝑥 + 𝑒 𝑥 = 0
que se sabe esta en el intervalo [-2, -1].
Solución
Paso 1: Verificación de las condiciones de convergencia
La función f(x)=1+x+ex es monótona, f(-2)=-1+e-2 <0, y f(-1)= e-1 >0.
Paso 2: Determinación del número de iteraciones
Como se desea alcanzar una precisión de 7x10-2se debe escoger n tal que:
− 1 − ( −2 )
n log2 n 3.8365
7 x10− 2
Efectuando 4 iteraciones a partir de [-2, -1] se tiene un error máximo absoluto de:
1
4
6 x10− 2
2
Laboratorio de Métodos Numéricos para Ingenieros Página: 5/7
Tema: Solución de Ecuaciones No Lineales f(x)=0
LAB N° 3 JCC - FSN
Método de la Bisección
Iteraciones:
n an f(an) bn f(bn) xn+1 f(xn+1)
1 -2.000 -0.865 -1.000 0.3679 -1.500 -0.2769
Podemos realizar una comprobación gráfica de la solución del problema, entonces, graficamos la función 1 + 𝑥 + 𝑒 𝑥 = 0
y vemos donde tiene un cruce por el eje x.
Figura 6 → Grafico de la Ec: f(x)=1+x+ex
La solución de la ecuación será x = , o sea, x [ , ]
III. MATERIAL Y EQUIPO
• Una PC con Sistema Operativo (S.O.) Windows X o superior.
• MATLAB 2016 o superior
IV. PROCEDIMIENTO
1. Localice gráficamente las raíces de 𝑓(𝑥) = 0, siendo 𝒇(𝒙) = |𝒙| − 𝒆𝒙 .
2. Usando el método de la bisección sucesivas, determine un valor aproximado para los ceros de 𝒇(𝒙) = |𝒙| − 𝒆𝒙 con
un error que no exceda a 0.015
3. Confeccionar en MATLAB un programa que:
3.1. Calcule la raíz de una ecuación algebraica no lineal o ecuación trascendente.
3.2. Adjunte su algoritmo y diagrama de flujo
Laboratorio de Métodos Numéricos para Ingenieros Página: 6/7
Tema: Solución de Ecuaciones No Lineales f(x)=0
LAB N° 3 JCC - FSN
Método de la Bisección
% Algoritmo de Solución % Diagrama de flujo
% Inicio de programa
% Fin de programa
V. CUESTIONARIO FINAL
1. Use cuatro iteraciones del método de la Bisección para determinar las raíces de 𝑒 2𝑥 − 6𝑥 = 0 en el intervalo [0, 0.5].
¿Cuántas iteraciones son necesarias para obtener la aproximación a la raíz redondeada a 5 cifras decimales?
2. Hacer una prueba de datos del programa realizado en el procedimiento con los datos de la pregunta anterior.
Laboratorio de Métodos Numéricos para Ingenieros Página: 7/7
Tema: Solución de Ecuaciones No Lineales f(x)=0
LAB N° 3 JCC - FSN
Método de la Bisección
VI. OBSERVACIONES Y CONCLUSIONES
Haga sus observaciones y emita al menos cinco conclusiones en torno al trabajo realizado
• ________________________________________________________________________________________________
• ________________________________________________________________________________________________
• ________________________________________________________________________________________________
• ________________________________________________________________________________________________
• ________________________________________________________________________________________________