“Año de la recuperación y consolidación de la economía peruana.
UNIVERSIDAD NACIONAL DE SAN AGUSTÍN
FACULTAD DE INGENIERÍA DE PROCESOS
ESCUELA PROFESIONAL DE INGENIERÍA DE MATERIALES
ASIGNATURA: METODOS NUMERICOS APLICADO A LA INGENIERIA
DOCENTE: Ing. LUIS ANTONIO LAZO ALARCON
TEMA: MÉTODO DE LA BISECCIÓN EN MATLAB
INTEGRANTE:
Ramos Ccolqque Jostin Ryan Matteo
Semestre 2025 – A
Arequipa – Perú
Introducción
Encontrar la raíz de una función, es decir, el valor de x para el cual f ( x )=0, es un problema muy
común en matemáticas, ingeniería y otras áreas científicas. Aunque en algunos casos es posible
obtener una solución exacta, muchas funciones no se pueden resolver de manera directa, por lo
que se recurre a métodos numéricos para obtener aproximaciones confiables.
Uno de los métodos más conocidos y utilizados para este propósito es el método de la bisección.
Se trata de un procedimiento simple pero muy efectivo, que parte de un intervalo donde la
función cambia de signo. A partir de ese intervalo, el método va dividiéndolo en mitades y
evaluando hasta acercarse progresivamente a la raíz. Aunque no es el método más rápido, sí es
uno de los más seguros y fáciles de entender, lo que lo convierte en una excelente herramienta
para aprender los fundamentos de la aproximación numérica.
En este trabajo, he implementado el método de la bisección en el programa MATLAB, creando
un código que permite al usuario ingresar una función, los extremos del intervalo, y el margen de
error deseado. El programa realiza las iteraciones necesarias, muestra los resultados paso a paso
y genera una gráfica de la función, lo que facilita visualizar cómo se comporta y cómo se
encuentra la raíz.
A través de este proyecto, busco no solo aplicar lo aprendido en el curso de Métodos Numéricos,
sino también mejorar mis habilidades en programación matemática y comprensión de algoritmos
que son fundamentales en el análisis computacional.
Objetivos
1. Comprender y aplicar el método de la bisección para encontrar raíces aproximadas de
funciones no lineales mediante la implementación del algoritmo en MATLAB.
2. Desarrollar competencias en la programación y análisis numérico, mediante la
implementación de un script en MATLAB que permita ingresar datos de forma
interactiva, evaluar funciones y visualizar gráficamente el proceso de aproximación hacia
la raíz.
Marco Teórico
El método de la bisección es una técnica numérica utilizada para encontrar raíces reales de
funciones continuas. Se basa en el teorema del valor intermedio, el cual establece que si una
función continua f (x) cambia de signo en un intervalo cerrado [ a , b ], es decir, si f (a). f (b)<0 ,
entonces existe al menos un númeroc ∈(a , b) tal que f ( c ) =0.
Este algoritmo pertenece a la categoría de métodos de búsqueda de raíces por intervalos, también
conocidos como métodos brutos, en los que se garantiza convergencia siempre que se cumpla la
condición inicial del cambio de signo. A diferencia de otros métodos más sofisticados como
Newton-Raphson o la secante, la bisección no requiere derivadas ni aproximaciones iniciales
precisas, lo que lo convierte en una opción robusta y segura, especialmente en etapas iniciales
del análisis numérico. (Chapra, 2015)
El procedimiento consiste en:
1. Verificar que f ( x 1 ) . f (x 2)<0, lo cual garantiza la existencia de una raíz en el intervalo
[ x 1 , x 2 ].
x 1+ x 2
2. Calcular el punto medio xp= .
2
3. Evaluar f ( xp ) :
Si f ( xp )=0 , se ha encontrado la raíz exacta.
Si f ( x 1 ) . f (xp)<0, la raíz está en [ x 1 , xp ] .
Si f ( x 2 ) . f (xp)<0 , la raíz está en [ xp , x 2 ] .
4. Repetir el proceso en el subintervalo que contiene la raíz, hasta que la longitud del
intervalo sea menor que una tolerancia predefinida ε .
Este proceso garantiza una convergencia segura, aunque no necesariamente rápida, ya que reduce
el tamaño del intervalo a la mitad en cada iteración.
Algunas ventajas son:
Fácil de entender e implementar.
Garantiza convergencia si se cumple la condición inicial.
No depende del valor de la derivada ni de una buena estimación inicial.
El método de la bisección se utiliza ampliamente en diversos campos de la ciencia e ingeniería
para resolver ecuaciones que surgen en problemas como la estática de estructuras, el análisis de
circuitos eléctricos, el ajuste de parámetros en modelos matemáticos, o la resolución de
ecuaciones de transferencia de calor y masa. (Burden, 2011)
Procedimiento experimental
1. Apertura del entorno de trabajo:
Se abrió el software MATLAB desde el escritorio del computador
2. Creación de un nuevo script:
En el entorno de MATLAB, se seleccionó la opción "New Script" en la pestaña Home
para crear un archivo, donde se escribirá el código del método de la bisección.
3. Impresión de un mensaje inicial
Se utiliza fprintf para mostrar un mensaje en la consola que identifica el programa. Es útil
para organizar la salida y brindar contexto al usuario.
4. Ingreso de la función y graficación inicial.
Input (...,'s') permite que el usuario ingrese una función como texto, por ejemplo:
2
x −4 .
Ezplot (Fx) grafica la función ingresada en un intervalo automático.
Grid activa la cuadrícula en el gráfico para mejor visualización.
5. Ingreso de intervalo inicial y tolerancia
Se solicita al usuario ingresar los extremos del intervalo [ x 1 , x 2 ] y el error de tolerancia e
, que será el criterio de paro del algoritmo.
6. Evaluación de los extremos iniciales
Aquí se evalúa la función en x 1 y x 2 .
Como Fx es una cadena de texto, se necesita asignar el valor a x y luego usar eval
(Fx) para obtener el valor numérico de la función en ese punto.
7. Encabezado de la tabla de iteraciones
Se imprime un encabezado con títulos de columnas que se utilizarán para mostrar los
datos de cada iteración: los extremos, el punto medio y los valores de la función.
8. Bucle principal: método de la bisección
El ciclo while continúa mientras la longitud del intervalo sea mayor que el error.
Se calcula el punto medio xp, se evalúa la función en ese punto y se almacena en
Fxp.
Se imprimen los valores de cada iteración con formato de 4 decimales.
Se aplica el criterio de cambio de signo:
Si el producto Fx 1∗Fxp ≤ 0, la raíz está en [ x 1 , xp ] , por lo x2 se actualiza
a xp.
Si no, se actualiza x 1=xp, y la raíz está en [ xp , x 2 ] .
9. Fin del ciclo e impresión del resultado
Al terminar el ciclo, se imprime la raíz aproximada encontrada con el formato
correspondiente.
10. Graficación final de la función
Se vuelve a graficar la función para mostrar visualmente dónde se encuentra la raíz.
Resultados
Para comprobar que el método de la bisección funciona correctamente, se realizaron dos pruebas
con la función f ( x )=x 3−2 x +1, utilizando el código implementado en MATLAB. Este código
fue diseñado para recibir los valores del intervalo, la función y el error tolerado, y luego mostrar
el proceso iterativo en pantalla junto con una gráfica de la función.
Primer caso: Intervalo [ 0 , 2 ]
Se ingresaron los siguientes datos:
Función: x 3−2 x+ 1
x 1=0
x 2=2
Error: 0.0001
Al ejecutar el programa, el primer paso fue graficar la función.
En la gráfica se observa que la curva cruza el eje X entre 0 y 1, lo cual ya sugiere visualmente
que existe una raíz en ese intervalo. El método, sin embargo, fue aplicado desde 0 hasta 2, lo
cual sigue siendo válido porque el cambio de signo entre los extremos garantiza la presencia de
una raíz.
Durante la ejecución, el código fue calculando el punto medio xp en cada iteración, evaluando la
función en ese punto, y decidiendo si debía actualizar el extremo izquierdo o derecho del
intervalo. Todo este proceso se repitió 15 veces, como se puede ver en la tabla de resultados
impresa por MATLAB.
Durante la ejecución, el código fue calculando el punto medio xp en cada iteración, evaluando la
función en ese punto, y decidiendo si debía actualizar el extremo izquierdo o derecho del
intervalo. Todo este proceso se repitió 15 veces, como se puede ver en la tabla de resultados
impresa por MATLAB.
En cada fila de la tabla, se muestran los valores actuales de x 1 , xp , x 2 , f ( x 1 ) , f ( xp ) y f (x 2), lo
que permite seguir paso a paso cómo se va reduciendo el intervalo y cómo se acerca el valor de
f (xp) a cero.
El programa detuvo las iteraciones cuando la diferencia entre los extremos fue menor al error
especificado. La raíz aproximada encontrada fue:
x=0.6180
Segundo caso: Intervalo [ −2 , 0 ]
Se ingresaron los siguientes datos:
x 1=−2
x 2=0
Error: 0.0001
El algoritmo se ejecutó con éxito, realizando nuevamente 15 iteraciones, donde el intervalo se
fue acortando y el punto medio fue acercándose a la raíz. La tabla generada muestra cómo el
valor de f (xp) se acerca a cero tras iteración.
La raíz encontrada fue:
x=−1.6180
Análisis general de los gráficos
En ambos casos, el gráfico de la función proporcionado por ezplot en MATLAB muestra con
claridad el comportamiento de la curva en los intervalos elegidos. Las raíces encontradas por el
método de la bisección coinciden visualmente con los puntos donde la función se anula, lo que
permite confirmar que la implementación del algoritmo es correcta.
Analisis final de los resultados:
Estos resultados demuestran que el método de la bisección, aunque simple, es altamente
confiable cuando se aplica a funciones continuas con cambio de signo en un intervalo. El hecho
de que el código proporcione tanto el seguimiento numérico como la visualización gráfica ayuda
a entender mejor el proceso de convergencia, algo muy valioso desde el punto de vista educativo.
Además, se evidencia que la implementación en MATLAB permite automatizar el proceso,
solicitar los datos de entrada de forma clara y mostrar el resultado con precisión, lo que refuerza
el aprendizaje tanto del algoritmo como del uso del software en aplicaciones matemáticas.
Cuestionario
1. ¿Qué es el método de la bisección y para qué se utiliza?
El método de la bisección es un procedimiento numérico que se utiliza para encontrar
raíces reales de funciones continuas. Se basa en el Teorema del Valor Intermedio, el cual
garantiza que si una función continua cambia de signo en un intervalo [ a , b ], entonces
existe al menos una raíz en ese intervalo. El método consiste en dividir el intervalo en dos
mitades sucesivas, seleccionando aquella que contenga el cambio de signo, hasta que la
solución se aproxime a la raíz con el error deseado.
2. ¿Qué condiciones deben cumplirse para aplicar correctamente el método de la bisección?
Para que el método de la bisección se aplique correctamente, deben cumplirse las
siguientes condiciones:
La función f (x) debe ser continua en el intervalo cerrado [ a , b ].
Debe existir un cambio de signo, es decir, f ( a ) . f (b)<0 . Esto asegura que hay al
menos una raíz en el intervalo.
Se debe establecer un criterio de paro, como un número máximo de iteraciones o
una tolerancia de error ( ∣b−a ∣<ε ).
3. ¿Qué ventajas tiene el método de la bisección frente a otros métodos numéricos?
Las principales ventajas del método de la bisección son:
Simplicidad: Es fácil de entender y programar.
Fiabilidad: Siempre converge si se cumplen las condiciones del cambio de signo.
Control del error: Permite controlar de forma directa el error máximo en la
estimación de la raíz.
Sin embargo, su convergencia es más lenta que la de otros métodos como el de Newton
Raphson o el de la secante.
4. ¿Cómo se calcula el número máximo de iteraciones necesarias para alcanzar un error
determinado?
El número máximo de iteraciones n necesarias para que el método de la bisección alcance
un error menor que una tolerancia ε , puede estimarse con la siguiente fórmula:
log ( b−a )−log ( ε )
n≥
log 2
5. ¿Por qué el método de la bisección puede fallar o no ser eficiente en algunos casos?
Aunque el método de la bisección es robusto, puede presentar algunas limitaciones:
Si la función no cambia de signo en el intervalo elegido, el método no puede
aplicarse.
Es lento, ya que su convergencia es lineal, lo que significa que reduce el error a la
mitad en cada iteración. Métodos como el de Newton tienen convergencia
cuadrática, mucho más rápida.
Solo encuentra una raíz por vez. Si hay varias raíces en el intervalo, solo detectará
una de ellas, dependiendo del punto medio y el cambio de signo.
Ejercicios Aplicativos (se utilizó el mismo código)
1. Control de temperatura en un reactor químico
En un reactor, se requiere mantener una temperatura donde la velocidad de reacción es
óptima. Esa temperatura se puede estimar con la siguiente función:
3 2
f ( x )=x −6 x +11 x−6
Intervalo [ 1.5 , 2.5 ]
Error: 0.0001
La raíz encontrada es X ≈ 1.9999 ºC, lo que indica la temperatura óptima de operación
ese rango.
Resultados
Interpretación de la Tabla de Iteraciones:
a. Primera iteración
Dividió el intervalo [ 1.5 , 2.5 ] en xp=2.0
f ( 1.5 )=0.375 (positivo).
f ( 2.0 )=0.000 (la raíz exacta).
f ( 2.5 )=−0.375 (negativo).
Seleccionó el subintervalo [ 1.5 , 2.5 ] porque f ( 1.5 )∗f (2.5)≤ 0
b. Iteraciones posteriores:
El método continuó dividiendo el intervalo por la mitad.
En cada paso, la aproximación se acercaba más a la raíz verdadera ( x=2.0 )
Los valores de f (xp) disminuían sistemáticamente hacia cero.
c. Convergencia:
El método requirió 14 iteraciones para alcanzar la precisión deseada.
La convergencia es monótona y garantizada por el método de bisección.
La aproximación final (1.9999) tiene un error absoluto de 0.0001 respecto a la
solución exacta (2.0).
Gráfico
Una curva cúbica que cruza el eje x en x=1, x=2 y x=3
La raíz encontrada ( x=2.0 ) marcada en el intervalo [ 1.5 , 2.5 ]
Comportamiento de la función que confirma el teorema de Bolzano (cambio de
signo)
2. Optimización de velocidad en una turbina hidráulica
En una planta hidroeléctrica, se desea encontrar la velocidad en la que la potencia
generada por una turbina es máxima. La función que modela este comportamiento se
puede representar así:
3 2
f ( x )=x −4 x + x+ 6
Intervalo [ 1.5 , 2.5 ]
Error: 0.0001
La velocidad crítica estimada es v ≈ 1.9999 m/s.
Resultados
Interpretación de los resultados:
a. Convergencia rápida:
El método encontró la raíz exacta ( x=2.0 ) en pocas iteraciones.
Esto ocurre porque 2 es exactamente el punto medio del intervalo inicial.
b. Verificación de la solución:
f ( 2 )=8−16 +2+6=0
El error es prácticamente cero.
Grafico
La función cruza el eje x exactamente en x=2
Comportamiento en el intervalo:
Creciente en [ 1.5 , 2.5 ]
Decreciente después de x=2.0
Bibliografía
Burden, R. L., & Faires, J. D. (2011). Análisis numérico (9.ª ed.). Cengage Learning.
Chapra, S. C., & Canale, R. P. (2015). Métodos numéricos para ingenieros (7.ª ed.).
McGraw-Hill Education.
