03 Apunte Método de Bisección y Regula Falsi

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

Departamento de Matemática

Unidad II: Resolución de ecuaciones no lineales

Apunte: Métodos de Bisección y Regula Falsi.

Objetivo:

- Estudiar los métodos de bisección y Regula Falsi para hallar raı́ces reales
de una ecuación no lineal en una variable

Introducción

Las asignaturas de matemática estudiadas hasta ahora se han dedicado a


desarrollar métodos que nos proporcionen la solución exacta de un problema.

Por ejemplo, calcular la solución de una ecuación del tipo f (x) = 0 rea-
lizando operaciones elementales sobre la misma para conseguir despejar la
incógnita x.
Rb
O bien, calcular el valor de una integral definida a f (x)dx calculando una
primitiva de f (x).

Infortunadamente, en la gran mayorı́a de los casos que se presentan en la


práctica, estos métodos no son de aplicación. Ello puede deberse a que el
método para calcular la solución exacta sea muy complicado, a que no se
conozca un método adecuado, o incluso a que no exista un método que nos
permita, mediante cálculos elementales, encontrar la solución.

En estos casos es necesario recurrir a métodos numéricos, denomina-


dos ası́ porque, usualmente, consisten en realizar una sucesión más o menos
larga de operaciones numéricas (normalmente mediante la ayuda de un orde-
nador), al cabo de las cuales encontramos un valor numérico que, si bien

1
Departamento de Matemática

no es la solución exacta del problema, se le parece mucho, es decir, aproxima


la solución buscada con una precisión razonablemente buena.

Definición de raı́z de una ecuación

Definición: Sea f : A 7−→ B, A, B ⊆ R una función dada. Un número


α ∈ A se dice que es raı́z de la ecuación f (x) = 0, o un cero de la función f
si f (α) = 0.

Un método simple para obtener una aproximación a la raı́z de la ecuación


f (x) = 0 consiste en graficar la función y observar dónde cruza el eje ”x” (ver
Figura 1). La técnica gráfica tiene un valor práctico limitado, ya que no es
precisa. Sin embargo, las aproximaciones de la raı́z obtenidas por el método
gráfico se pueden usar como valores iniciales en los métodos numéricos que
se estudiaran en ésta unidad.

Los métodos numéricos que en cada paso dan un intervalo cerrado donde
se encuentra la raı́z buscada, son llamados métodos cerrados. Aquı́ estudia-
remos dos de tales métodos: el método de Bisección y el método de Regula
Falsi (falsa posición).

Método de Bisección: Supongamos que f (x) es una función continua


en el intervalo cerrado [a, b] y f (a) · f (b) < 0. Entonces, por el Teorema
de Bolzano (conocido también como el teorema del valor intemedio) existe
al menos un número α ∈ (a, b) tal que f (α) = 0. Asumiremos en lo que
sigue que la raı́z en éste intervalo es única (en el desarrollo de los ejemplos se
demostrará la unicidad de dicha raı́z).

El método de Bisección aplicada a la función f (x) para aproximar


la raı́z α ∈ (a, b), consiste en dividir sucesivamente el intervalo [a, b] por
la mitad, hasta que la longitud del subintervalo que contiene a la raı́z α sea
menor que alguna tolerancia especificada .

2
Departamento de Matemática

Figura 1. Método gráfico para hallar la raı́z de una ecuación

Entonces, tomando a1 = a, b1 = b y x1 el punto medio de [a1, b1], o sea


a1 + b1 b1 − a1
x1 = = a1 +
2 2
es la primera aproximación de la raı́z α.

b1 − a1
Si f (x1) = 0 o < , entonces α = x1 y el proceso termina.
2
En el caso que no suceda lo anterior:

Si f (a1) · f (x1) < 0 entonces α ∈ (a1, x1) y tomamos a2 = a1, b2 = x1;


en caso contrario f (x1) · f (b1) < 0 tomamos a2 = x1, b2 = b1.

En general, después de (n-1)-iteración, la raı́z α ∈ (an, xn) y tomamos

3
Departamento de Matemática

an + bn bn − an
xn = = an +
2 2
como la n-ésima aproximación de la raı́z α, además
bn − an b − a
|α − xn| ≤ = n ≤ (1)
2 2
.
b−a
Nota 1: Como lı́m = 0, entonces lı́m xn = α, es decir, la su-
n→+∞ 2n n→+∞
cesión {xn}n converge a la raı́z α; lo que significa que el método de Bisección
siempre converge.

Nota 2: El número de iteraciones necesarias para obtener una aproxima-


ción de la raı́z α menor o igual a la tolerancia  se obtiene despejando n en
(1), es decir,

 
b−a
n ≥ log2 (2)


A modo de resumen podemos decir que el método de bisección consiste en


lo siguiente:

1. Subdividir en dos partes el intervalo en que se sabe que la función cambia


de signo y tiene una sola raı́z.
2. Averiguar, utilizando el Teorema de Bolzano, en cual de las dos mitades
se encuentra la raiz y descartar la otra mitad del intervalo.
3. Reiniciar este proceso con el subintervalo elegido.
4. Continuar con este proceso hasta que el subintervalo elegido tenga una
longitud lo suficientemente pequeña como para que cualquiera de sus
puntos sea una aproximación aceptable de la solución. La elección óptima
como aproximación es, entonces, el punto medio del subintervalo.

4
Departamento de Matemática

Ejemplo: Encuentre una aproximación de la solución de la siguiente ecua-


ción dado el intervalo y la tolerancia,

x − e−x = 0, [0, 1],  = 0. 1

Primer paso: Identificar la función f , en éste ejemplo tenemos que f (x) =


x − e−x.

Segundo paso: Comprobar las hipótesis del Teorema de Bolzano, es decir


ver si la función es continua en el intervalo dado y si cambia de signo en los
extremos de dicho intervalo.

- De estudios anteriores no es difı́cil ver que la función f no se indefine en


ningún punto, por tanto es continua en todo R y por tanto es continua en el
intervalo [0, 1].

- Como f (0) = 0 − e0 = −1 y f (1) = 1 − e−1 ≈ 0.6321, entonces la


función f cambia de signo en los extremos del intervalo.

Quedando ası́ comprobada las hipótesis del Teorema de Bolzano, por lo que
podemos decir que existe al menos una raı́z de la ecuación en el intervalo.
Para garantizar que existe una única raı́z basta con demostrar que la función
f es estrictamente monótona en [0, 1], es decir, que en dicho intervalo o bien
la función es siempre decreciente o bien es siempre creciente.

Es conocido de cursos anteriores que podemos obtener la monotonı́a de


una función con ayuda de la primera derivada de dicha función, en éste caso
tenemos que,

f 0(x) = 1 + e−x,

es decir la primera derivada de f siempre es positiva, por tanto la función


es estrictamente monótona creciente, lo cual demuestra que la ecuación tiene

5
Departamento de Matemática

una única raı́z en el intervalo [0, 1].

Note que utilizando la desigualdad (2), obtenemos que para éste ejemplo
se necesitan 4 iteraciones para obtener un error menor a 0.1, debido a que,

n > log2 10 ≈ 3.3219

En la siguiente tabla se resume los cálculos realizados para obtener la apro-


ximación de la raı́z de la ecuación con la tolerancia indicada.

b−a
n an bn xn f (an) f (bn) f (xn) 2n
1 0 1 0.5 -1 0.6321 -0.1065 0.5
2 0.5 1 0.75 -0.1065 0.6321 0.2776 0.25
3 0.5 0.75 0.625 -0.1065 0.2776 0.0897 0.125
4 0.5 0.625 0.5625 -0.1065 0.0897 -0.0073 0.0625

Finalmente, se ha obtenido la aproximación de la raı́z de la ecuación igual


a α = 0.5625 para una tolerancia menor o igual a 0.1.

Método de la Falsa Posición (o Regula Falsi)

Aun cuando la bisección es una técnica perfectamente válida para deter-


minar raı́ces, su método de aproximación por “fuerza bruta” es relativamente
ineficiente. La falsa posición es una alternativa basada en una visualización
gráfica.

Método de la Falsa Posición (o Regula Falsi): Consideremos una


función f continua en un intervalo [a, b] y tal que f (a)f (b) < 0. El método
de la Falsa Posición, para encontrar una aproximación de una raı́z α ∈ (a, b)
de f (x) = 0, es similar al método de Bisección en el sentido de que se generan
subintervalos [an, bn] que encierran a la raı́z α, pero ésta vez xn no es el punto
medio de [an, bn], sino el punto de interseccion de la recta que pasa por los

6
Departamento de Matemática

punto (an, f (an)), (bn, f (bn)) con el eje x (ver la Figura 2)

Figura 2.

Empezamos tomando a1 = a, b1 = b y encontramos la primera aproxima-


ción de la raı́z, x1, como la intersección con el eje x, de la recta secante a f
que pasa por los puntos (a1, f (a1)), (b2, f (b2)), es decir:

(b1 − a1) f (a1) a1f (b1) − b1f (a1)


x1 = a1 − =
f (b1) − f (a1) f (b1) − f (a1)

Si f (x1) = 0, entonces α = x1 y el proceso termina.

Si f (a1)f (x1) < 0 entonces α ∈ (a1, x1) y tomamos a2 = a1, b2 = x1, de


lo contrario tomamos a2 = x1, b2 = b1.

Aplicamos nuevamente el proceso anterior al intervalo [a2, b2], es decir,


hacemos

7
Departamento de Matemática

(b2 − a2) f (a2)


x2 = a2 −
f (b2) − f (a2)
Después de la (n − 1)-ésima iteración, tenemos α ∈ (an, bn) y tomamos a

(bn − an) f (an)


xn = an −
f (bn) − f (an)

Éste método tiene la desventaja, con respecto al método de Bisección, que


la longitud del subintervalo que contiene a la raı́z en general no tiende a cero,
porque las gráficas de las funciones son cóncavas (hacia arriba o hacia abajo)
en la vecindad de la raı́z, lo que hace que uno de los extremos de los subin-
tervalos se aproxime a la raı́z, mientras el otro permanece fijo (ver la Figura 2).

Por lo anterior, la longitud del subintervalo [an, bn] no puede tomarse como
un criterio de aproximación a la raı́z; se requiere una tolerancia en el valor
de la función en la aproximación xn, es decir |f (xn)| <  o |xn − xn−1| < 
para alguna tolerancia  > 0 previamente escogida. El procedimiento ter-
mina cuando se alcance esta tolerancia o un número máximo de iteraciones
previamente establecido.

Ejemplo: Utilizando el método de de Regula Falsi y el ejemplo anterior,


obtenga una aproximación de la solución de la siguiente ecuación dado el
intervalo y la tolerancia,

x − e−x = 0, [0, 1],  = 0.1

Usando Regula Falsi, los resultados se resumen en la siguiente tabla.

n an bn xn f (an ) f (bn ) |f (xn )| < 0,1


1 0 1 0.6127 -1 0.6321 0.0708

8
Departamento de Matemática

Note que solo se necesitó una iteración del método de la falsa posición para
obtener una aproximación de la raı́z con una tolerancia  = 0.1, mientras que
se necesitaron 4 iteraciones usando el método de Bisección. Los estudiantes
podrán verificar que la aproximación x3 = 0.5677 obtenida con el metodo
de Regular Falsi es obtenida con una tolerancia menor a 10−3. Aunque un
método como el de la Falsa Posición casi siempre es superior al de Biseccion,
hay algunos casos que violan esta conclusión general.

Ejercicios que pueden ser desarrollados en la clase:

Utilizando los métodos de Bisección y Regula Falsi, encuentre una aproxi-


mación de la solución para las siguientes ecuaciones en los intervalos indicados
y considerando las tolerancias dadas,

a) x3 + 4x2 − 12 = 0, [1, 2],  = 0.05

b) ln(x) − x + 2 = 0, [0,1, 1],  = 0.05

También podría gustarte