Analisis Numerico Metodos Abiertos
Analisis Numerico Metodos Abiertos
Analisis Numerico Metodos Abiertos
2019
Resumen
Aproximaciones sucesivas es un método abierto, es decir, no necesita de un intervalo que atrape una
raı́z, sino que requiere de un valor x0 que representa una aproximación a la raı́z; de la cercanı́a de
ésta a la raı́z dependerá la velocidad en que se cumpla con una tolerancia preestablecida.
Una forma sencilla de definir un método de aproximaciones sucesivas consiste en despejar de una
ecuación a la variable independiente; esto se aplica particularmente en ecuaciones que por su forma
*
Profesores de la División de Ciencias Básicas de la Facultad de Ingenierı́a de la UNAM.
1
Análisis numérico 2
ex − x2
x2 + 7x − ex = 0 =⇒ x =
7
exi − x2i
xi+1 =
7
donde xi es un valor inicial y xi+1 es un valor corregido que, en un escenario favorable, tendrá
una cantidad de error menor con respecto a la raı́z de la ecuación. El proceso iterativo se detendrá
cuando entre dos aproximaciones sucesivas (de aquı́ el nombre del método) se satisfaga la tolerancia
preestablecida.
No obstante la aparente facilidad que se muestra, el método no es ciento por ciento eficaz en todas
las ecuaciones, como se verá posteriormente.
Una mamera de obtener una ecuación de recurrencia general es la siguiente:
Sea f (x) una función algebraica o trascendente:
f (x) = 0 (1)
f (x) + x = x (2)
Definiendo al término:
G(x) = f (x) + x (3)
Sustituyendo (3) en (2) se tiene:
G(x) = x (4)
La ecuación (4) representa el método de Aproximaciones sucesivas, por lo cual debe expresarse en
forma iterativa:
xi+1 = G(xi ) (5)
La aplicación del método consiste en proporcionar una aproximación inicial a la raı́z de la ecuación
(que puede obtenerse por medios gráficos o al detectar un cambio de signo en la función tabular) y
sustituirla en la ecuación (5), obteniéndose una nueva aproximación. De nuevo deberá sustituirse esta
última hasta que la diferencia entre dos aproximaciones sucesivas satisfaga determinada tolerancia
preestablecida. Es importante aclarar que aún cuando se utilice la ecuación de recurrencia (5) las
raı́ces corresponden a la función original f (x).
Análisis numérico 3
Sea f (x) una ecuación algebraica o trascendente que tiene como raı́z real al número a, sustituyendo
en la ecuación (4):
a = G(a) (6)
Restando (5) a la ecuación (6):
a − xi = G(a) − G(xi−1 ) (7)
Multiplicando el segundo miembro de (7) por el factor unitario:
a − xi−1
a − xi−1
G(a) − G(xi−1 )
a − xi = · (a − xi ) (8)
a − xi
Aplicando el teorema del valor medio del cálculo diferencial:
G(a) − G(xi−1 )
G0 (τ ) = ; xi−1 ≤ τ ≤ a (9)
a − xi−1
Sustituyendo la ecuación (9) en (8):
Despejando G0 (τ ):
a − xi
G0 (τ ) = ; xi−1 ≤ τ ≤ a (11)
a − xi−1
En el segundo miembro de la ecuación (11) puede observarse que su denominador debe ser mayor
que el numerador, toda vez que xi−1 posee un mayor error que xi ya que es una aproximación
previa. Esto implica que:
| a − xi |
| G0 (τ ) |= <1 (12)
| a − xi−1 |
En consecuencia, el método convergerá si se cumple que:
Este criterio de convergencia debe probarse para cada una de las aproximaciones a cada raı́z.
Los casos de convergencia y divergencia se explican en las siguientes figuras; en todas ellas, el éxito
o fracaso del método depende del valor de G0 (τ ).
Consideremos como ejemplo una función sencilla que nos permita verificar resultados fácilmente
(Olivera Salazar, s.f.) (Garcı́a B., 2017). Se propone f (x) = x2 − 0,5. Se percibe que este polinomio
de segundo grado representa a una parábola que abre hacia arriba; naturalmente, posee dos raı́ces
√
cuyos valores son ± 0,5.
Ahora bien, suponiendo desconocida esta información, se realizará la exploración de la función para
encontrar sus raı́ces. El paso más recomendado es graficar la función.
Se propone como valores iniciales x0 = −1 para la raı́z negativa y x0 = 1 para la raı́z positiva.
Sea f (x) = x2 − 0,5 y en consecuencia G(xi = x2 + x − 0,5 y la ecuación de recurrencia:
xi+1 = x2i + x − i − 0,5
Las iteraciones para la obtención de la aproximación a la raı́z negativa se detallan en el cuadro 1.
La aproximación a la raı́z es −0,70701 con un error absoluto de 0,00033 después de 8 iteraciones.
Las iteraciones para la obtención de la aproximación a la raı́z positiva se detallan en el cuadro 2.
Es evidente que para la raı́z positiva el método no converge
La explicación se obtiene a partir del criterio de convergencia:
Para la raı́z negativa
1.4. Conclusiones
2. Método de Newton-Raphson
El método Newton Raphson (N-R) es, junto con Bisección, uno de los más utilizados. Su preferencia
radica en su robustez y velocidad para encontrar la raı́z cuando la aproximación cumple con su
criterio de convergencia. Se aplica a ecuaciones algebraicas y trascendentes y proporciona raı́ces
reales y complejas.
El nombre original del método N-R es de las tangentes. La tangente es una recta que intersecta a
una curva en un sólo punto; en consecuencia, es perpendicular a su radio. A partir de la figura 5, se
plantea que en un valor x0 que represente una aproximación a la raı́z de la ecuación, se trace una
tangente en el punto f (x0 ).
Esta recta tangente deberá cortar al eje horizontal y el punto donde esto ocurra será la nueva aproxi-
mación x1 , de tal forma que en el punto f (x1 ) se trace una nueva tangente. Este proceso se repetirá
hasta que el corte de la tangente en el eje horizontal coincida con la raı́z de la ecuación, o bien,
cuando la diferencia entre dos aproximaciones sucesivas cumpla con una tolerancia preestablecida.
De nuevo a partir de la figura 6, con base en las dos primeras iteraciones, se define la siguiente
relación entre el triángulo formado por la recta tangente y el ángulo
theta:
f (x0 )
tan(θ) = (17)
x0 − x1
Por otra parte, se conoce que:
f 0 (x0 ) = tan(θ) (18)
Análisis numérico 8
f (xi )
xi+1 = xi − (20)
f 0 (xi )
Por ser un método de punto fijo, el criterio de convergencia que deberá cumplirse es:
Para adaptar la ecuación (21) al método N-R, se sustituye la aproximación τ por xi+1 , de acuerdo
a lo siguiente:
f (xi )
| G(xi ) |= xi − 0 (22)
f (xi )
De tal forma que debe obtenerse la derivada de la ecuación (22):
f (xi ) · f 00 (xi )
G0 (xi ) = (23)
[f 0 (xi )]2
Análisis numérico 9
Encontrar una raı́z de la ecuación f (x) = sen(x) · e−x + 1. La gráfica de la ecuación es:
En la figura 7 se observa que la función tiene tres raı́ces reales en los intervalos [−8, −6], [−4, −2] y
[−2, 0]. Calculemos cada una de ellas, considerando el criterio de equivalencia en cada una de ellas
de acuerdo a las siguientes expresiones:
Las ecuaciones (25), (26) y (27) deben sustituirse cada una en la ecuación (24) tomando como
primera aproximación a la raı́z el punto medio de cada uno de los intervalos 1 y, si el resultado
cumple con el criterio de convergencia, utilizar esta aproximación en la ecuación de recurrencia:
sen(x) · e−x + 1
xi+1 = xi − (28)
e−x · [cos(x) − sen(x)]
1
Este puede ser un criterio cómoda para localizar una primera aproximación, aunque no es necesario contar con el
intervalo
Análisis numérico 10
1. x = −6,2850492734
2. x = −3,0963639324
3. x = −0,5885327440
Figura 10: Cuatro raı́ces múltiples Figura 11: Cinco raı́ces múltiples
Por otra parte, en las funciones mostradas, conforme crece el número de raı́ces múltiples, en números
pares (figuras 8 y 10) y en nones (figuras 9 y 11), la pendiente de la tangente tiende a cero; aunque
la raı́z no sea 0, el método fracasará.
Esto no quiere decir que no se pueda resolver una función con raı́ces múltiples, existe una versión
modificada del método Newton-Raphson que contempla esta situación.
3.1. Conclusiones
Como lo muestran las soluciones anteriores, el método de N-R resulta una herramienta ágil y robusta
en el cálculo de raı́ces. Por este motivo es el algoritmo preferido de los fabricantes de calculadoras
programables. Para los casos en que no resulta convergente, su complemento ideal resulta del método
de Bisección.
De esta forma, se considera que se pueden obtener las raı́ces reales de prácticamente cualquier
ecuación algebraica o trascendente.
Notas
1
Las figuras y gráficas incluidas en este trabajo fueron elaboradas por los autores
Referencias
Borras, H., Duran, R., y Iriarte, R. (1984). Apuntes de métodos numéricos (F. de Ingenierı́a UNAM,
Ed.).
Burden, R., y Faires, D. (2011). Análisis numérico (C. Learning, Ed.).
Chapra, S., y Canale, R. (2015). Métodos numéricos para ingenieros (M. Hill, Ed.).
Garcı́a B., S. (2017). Métodos numéricos.
Gerald, C. (1991). Análisis numérico (Alfaomega, Ed.).
Luthe, R., Olivera, A., y Schutz, F. (1985). Métodos numéricos.
Olivera Salazar, A. (s.f.). Métodos numéricos (Limusa, Ed.).
Sandoval, H. (2017). Métodos numéricos.