2.3. Método Intervalo Medio

Descargar como pdf o txt
Descargar como pdf o txt
Está en la página 1de 4

Métodos Numéricos para Ingenieros

Raúl Barahona

Métodos que utilizan intervalos

Al buscar valores de raíces que hacen a una función igual a cero, se analizan
técnicas numéricas que aprovechan el hecho de que una función, típicamente,
cambia de signo en la vecindad de una de esas raíces. Estas técnicas se conocen
como métodos que utilizan intervalos, pues, su funcionamiento requiere de dos
valores iniciales. Como su nombre lo indica, estos valores deben “delimitar" o estar
uno de cada lado de la raíz (denominados límite inferior y superior).

En esta sección estudiaremos los métodos que son clasificados como métodos que
utilizan intervalos para obtener raíces, pues, se requiere de dos valores iniciales
para obtener una aproximación a la raíz. Estos son: intervalo medio o Bisección y
Regula Falsi o Falsa Posición.

Antes de iniciar el proceso se requiere comprobar que el valor raíz se encuentre


dentro del intervalo a evaluar, de lo contrario no se podrá encontrar la raíz de forma
rápida y se tendrán que realizar muchas iteraciones hasta encontrar el valor. Sin
embargo, existe otro tipo de métodos que no requieren intervalos para el proceso
de búsqueda de la raíz. Mas bien, utilizan un valor cualquiera para iniciar el proceso
y que, debido al comportamiento de estos se obtendrá una aproximación bastante
exacta. Estos métodos se conocen como métodos abiertos; métodos de Newton-
Raphson y Secante.

2.3. Método de intervalo medio (bisección)

Este método consiste en dividir el intervalo en dos partes iguales o punto medio
reteniendo la mitad en donde f(x) cambia de signo, para conservar al menos una
raíz y repetir el proceso varias veces.

Una situación típica se describe en la siguiente ilustración, en donde el punto medio


del intervalo (xi, xd), se denota por xm = (xi + xd)/2 y f(xm) es positivo. Ya que el
producto f(xi)*f(xm)<0, debe conservar (xi, xm) como sub-íntervalo que contiene al
menos una raíz. El siguiente paso es evaluar f(xm) y un nuevo intervalo (xi, xm), etc.

intervalo para iniciar (xi, xd) o lo que el mismo (a, b)


xm = (xi, + xd) / 2
siguiente intervalo (xi, xm)

1
Métodos Numéricos para Ingenieros

Raúl Barahona

Intervalo inicial
(xi=a, xd=b)

Sub intervalo Sub intervalo


(xi, xm) (xm, xd )

a o xi xm xm b o xd

(x =x
xm = (xi + xd) /2

El método es algo lento, sin embargo, en cada paso o iteración se reduce la longitud
del intervalo a la mitad. Es común utilizar el método de intervalo medio (Bisección)
para localizar una aproximación a la raíz dentro de un intervalo mayor o igual a
0.001 (llamado épsilon Ɛ). Lo ideal es que f(x) = 0. Con el objeto de obtener
intervalos finales más pequeños. Se podrían emplear otros métodos para reducir la
longitud del intervalo de una manera más eficiente.

Algoritmo del método de intervalo medio (bisección)


1. Dado un intervalo inicial (xi, xd) de tal forma que la función cambie de signo
dentro del mismo, lo que se puede verificar con la expresión: f(xi)*f(xd)<0, o
de otra forma se debe corregir el intervalo inicial
2. Obtener el valor de referencia Ɛ que condiciona detener el proceso.
3. Primera aproximación a la raíz xm, mediante: xm = (xi + xd) / 2, ir al paso 6
4. Para continuar el proceso, se deben verificar los siguientes criterios, a fin,
de determinar el siguiente sub-íntervalo.
• Si f(xi)*f(xm) > 0, continuamos en el sub-intervalo izquierdo, para ello
hacemos, xi = xm, es decir (xi, xd)
• Si f(xi)*f(xm) < 0, continuamos en el sub-intervalo derecho, para ello
hacemos xd = xm, es decir (xi, xd)
5. Con el nuevo intervalo (xi, xd) Obtener otra aproximación: xm = (xi + xd)/2
6. Si |f(xm)| ≤ Ɛ, (lo ideal sería un Ɛ = 0), xm es la raíz (aproximada), de otra forma
ir al paso 4

Los ejemplos a continuación demuestran el uso del método. En cada paso observen
los valores que determinan en siguiente intervalo.

2
Métodos Numéricos para Ingenieros

Raúl Barahona

Ejemplo: obtener una aproximación a la raíz para la ecuación f(x) = x3 - x - 1, que


esté entre [0, 3], con una aproximación, no menor que Ɛ = 2x10-2 (0.02), y cuatro
cifras después del punto decimal.
i a b xm f(a) f(xm) f(a)*f(xm)
0 0 3 1.5 -1 0.875 -0.8750
1 0 1.5 0.75 -1 -1.3281 1.3281
2 0.75 1.5 1.1250 -1.3281 -0.7012 0.9312
3 1.125 1.5 1.3125 -0.7012 -0.0515 0.0361
4 1.3125 1.5 1.4063 -0.0515 0.3747 -0.0193
5 1.3125 1.4063 1.3594 -0.0515 0.1526 -0.0079
6 1.3125 1.3594 1.3359 -0.0515 0.0483 -0.0025
7 1.3125 1.3359 1.3242 -0.0515 -0.0021 0.0001
8 1.3242 1.3359 1.3301 -0.0021 0.0230 0.0000
9 1.3242 1.3301 1.3271 -0.0021 0.0104 0.0000
10 1.3242 1.3271 1.3257 -0.0021 0.0041 0.0000
11 1.3242 1.3257 1.3250 -0.0021 0.0010 0.0000
12 1.3242 1.3250 1.3246 -0.0021 -0.0006 0.0000
13 1.3246 1.3250 1.3248 -0.0006 0.0002 0.0000
14 1.3246 1.3248 1.3247 -0.0006 -0.0002 0.0000
15 1.3247 1.3248 1.3247 -0.0002 0.0000 0.0000

La raíz aproximada se obtuvo en la iteración 11, pues, f(xm) es menor a 0.02 o sea
< Ɛ. Sin embargo, si cotinuamos se obtiene f(xm) = 0, por lo tanto es mas exacta en
4 decimales.

Ejemplo: dada la función f(x) = x3 + x2 - 3x - 3. Encuentre la raíz más próxima que


esté entre x = 1 y x = 2. Con una estimación de Ɛ = 0.02

i xi o a xd o b xm f(a) f(x m) f(a)*f(x m)


0 1 2 1.5 -4 -1.875 7.5000
1 1.5000 2.0000 1.7500 -1.8750 0.1719 -0.3223
2 1.5000 1.7500 1.6250 -1.8750 -0.9434 1.7688
3 1.6250 1.7500 1.6875 -0.9434 -0.4094 0.3862
4 1.6875 1.7500 1.7188 -0.4094 -0.1248 0.0511
5 1.7188 1.7500 1.7344 -0.1248 0.0220 -0.0027
6 1.7188 1.7344 1.7266 -0.1248 -0.0518 0.0065
7 1.7266 1.7344 1.7305 -0.0518 -0.0150 0.0008
8 1.7305 1.7344 1.7324 -0.0150 0.0035 -0.0001
9 1.7305 1.7324 1.7314 -0.0150 -0.0057 0.0001
10 1.7314 1.7324 1.7319 -0.0057 -0.0011 0.0000
11 1.7319 1.7324 1.7322 -0.0011 0.0012 0.0000
12 1.7319 1.7322 1.7321 -0.0011 0.0000 0.0000

La raíz aproximada se obtuvo en la iteración 12, pues, f(xm) es menor a 0.02 o sea
< Ɛ. Siendo esta, la mas exacta en 4 decimales.

3
Métodos Numéricos para Ingenieros

Raúl Barahona

Ejemplo: encuentre una aproximación a la raíz dentro del intervalo (-0.56, 0.57),
para f(x) = ex - 1 = 0. No más de 3 cifras posteriores al punto decimal.

i xi o a xd o b xm f(a) f(x m) f(a)*f(x m)


0 -0.560 0.570 0.005 -0.429 0.005 -0.002

1 -0.560 0.005 -0.278 -0.429 -0.242 0.104


2 -0.278 0.005 -0.136 -0.242 -0.127 0.031
3 -0.136 0.005 -0.066 -0.127 -0.064 0.008
4 -0.066 0.005 -0.030 -0.064 -0.030 0.002
5 -0.030 0.005 -0.013 -0.030 -0.013 0.000
6 -0.013 0.005 -0.004 -0.013 -0.004 0.000
7 -0.004 0.005 0.001 -0.004 0.001 0.000
8 -0.004 0.001 -0.002 -0.004 -0.002 0.000
9 -0.002 0.001 -0.001 -0.002 -0.001 0.000
10 -0.001 0.001 0.000 -0.001 0.000 0.000

La raíz aproximada se obtuvo en la iteración 10, pues, f(xm) es menor a 0.000 o sea
< Ɛ. Siendo esta, la mas exacta en 3 decimales.

También podría gustarte