Me Todos
Me Todos
Me Todos
Si f es una función continua sobre el intervalo [a;b] y si f(a).f(b)<0, entonces f
debe tener un cero en (a;b). Dado que f(a).f(b)<0, la función cambia de signo en el
intervalo [a;b] y, por lo tanto, tiene por lo menos un cero en el intervalo. Ésta es una
consecuencia del teorema del valor intermedio para funciones continuas.
El método de bisección explota esta idea ya que los extremos del intervalo se van
a ir acercando sistemáticamente hasta obtener un intervalo de longitud
suficientemente pequeña en el que se localiza un cero. El proceso de decisión para
subdividir el intervalo consiste en tomar el punto medio del intervalo c = (a+b)/2 y
luego analizar las tres posibilidades que pueden darse:
Si f(a) y f(c) tienen signos opuestos, entonces hay un cero en [a;c].
Si f(c) y f(b) tienen signos opuestos, entonces hay un cero en [c;b].
Si f(c)=0, entonces c es un cero.
Si ocurre la primera opción, o bien la segunda, entonces se ha encontrado un
intervalo de mitad de ancho que el original que contiene una raíz. Para continuar el
proceso, se renombra el nuevo intervalo más pequeño como [a;b] y se repite el
proceso hasta que el intervalo sea tan pequeño como se desee. Puesto que el proceso
de bisección genera una sucesión de intervalos encajados, con sus correspondientes
puntos medios, se usará la siguiente notación para tener registro de los detalles del
proceso:
[a1;b1] es el intervalo de partida y c1=(a1+b1)/2 es su punto medio.
[a2;b2] es el segundo intervalo y c2=(a2+b2)/2 es su punto medio; el
intervalo [a2;b2] es la mitad de ancho que [a1;b1].
Después de llegar al intervalo [an;bn], en el que también se localiza un
cero y cuyo punto medio está dado por c n=(an+bn)/2, se construye el intervalo
[an+1;bn+1], en el que también se localiza un cero, y que mide la mitad que
[an;bn]
Animación del método
Análisis del error
Con el objeto de analizar el método de la bisección, se llamará [a i;bi], con i = 1,
2, 3, …, a los intervalos que surgen en el proceso. Se pueden hacer las siguientes
observaciones sobre estos números:
(1) Dado que r (la raíz buscada) está en una de las dos mitades del intervalo [a n,
bn], la distancia del punto medio del intervalo a r debe ser menor o igual que la
mitad de la longitud del intervalo.
(2) Por la igualdad demostrada en (*)
Ventajas y desventajas
Una gran ventaja de este método es que siempre converge para funciones
continuas f(x). Además, proporciona el tamaño exacto del intervalo en cada iteración
(en ausencia de errores de redondeo). Para aclarar esto, se puede observar que en
este método después de cada iteración el tamaño del intervalo se reduce a la mitad;
después de n iteraciones, el intervalo original se habrá reducido 2 n veces. Por lo
anterior, si el intervalo original es de tamaño M y el criterio de convergencia aplicado al
valor absoluto de la diferencia de dos aproximaciones sucesivas es ε, se puede saber
de antemano el número de iteraciones que se requieren:
Así, el proceso iterativo producirá puntos Pi= (xi; f(xi)) hasta que el último punto
Pn se encuentre en la banda horizontal comprendida entre las rectas de ecuaciones y=ε
e y=-ε.
Otro criterio de parada involucra las abscisas. Se trata de determinar si la
sucesión {xn} converge y, para ello, se dibujan dos rectas verticales de ecuaciones
x=P+β y x=P-β a cada lado de x=P. De esta manera, el proceso se detendrá cuando P n
está entre ambas rectas.
El método de Newton
Se toma una aproximación r0 del cero p en el intervalo [a, b], tal que f´(r 0) ≠ 0 y
|r0- r| es pequeño. Entonces, se puede construir el primer polinomio de Taylor de f(x)
alrededor de r0:
Por el teorema de Taylor, se sabe que:
con ζn entre rn y r. Despejando de esta ecuación, resulta:
Los métodos cuyo orden de convergencia es mayor que uno pero menor que dos
reciben el nombre de superlineales.
El método de Regula - Falsi
Otro algoritmo popular para resolver ecuaciones no lineales es el método de
Regula - Falsi o método de la posición falsa.
Al igual que en el método de bisección, se supondrá que f(a) y f(b) tienen distinto
signo. Si bien en el método de bisección se usa el punto medio del intervalo [a;b] para
llevar a cabo el siguiente paso, es posible conseguir una mejor aproximación usando el
punto (c;0) en el que la recta secante L que pasa por los puntos (a;f(a)) y (b;f(b))
cruza el eje OX.
Para hallar el punto c, se igualan las dos fórmulas para la pendiente m de la recta
L:
que resulta de usar los puntos (c;0) y (b;f(b)). Igualando las dos expresiones
anteriores resulta:
Ventajas y desventajas
La ventaja del método de Regula - Falsi, al igual que el de bisección, es que es
siempre convergente para funciones continuas f(x). Aunque en general, converge más
rápidamente que el método de la bisección, su velocidad de convergencia es baja.
Orden de convergencia
Se puede demostrar que bajo ciertas condiciones el método de Regula - Falsi
tiene orden de convergencia lineal. Esto sucede cuando la función f(x) es cóncava
entre a y b y el punto a es siempre uno de los dos puntos usados para la siguiente
iteración. Lo mismo ocurriría si fuese convexa en las inmediaciones de la raíz.