Tarea 3 Metodo de Newton Raphson

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

INSTITUTO POLITECNICO

NACIONAL
ESCUELA SUPERIOR DE INGENIERIA
MECANICA Y ELECTRICA

INVESTIGACION DE
METODO DE NEWTON
RAPHSON

Alumno:

 Garay López Aldo Arturo.


Profesor: Yáñez López José Luis.

Unidad de aprendizaje: Análisis Numérico.

Boleta: 2019302548.
Grupo: 4CV1.
1.1 METODO DE NEWTON RAPHSON

Este método es uno de los más utilizados para localizar raíces ya


que en general es muy eficiente y siempre converge para una
función polinomial.
Se requiere que las funciones sean diferenciables, y por tanto,
continuas, para poder aplicar este método.
Se debe partir de un valor inicial para la raíz: xi , este puede ser
cualquier valor, el método convergirá a la raíz más cercana.
Si se extiende una tangente desde el punto b, el punto donde esta
tangente cruza al eje x representa una aproximación mejorada de la
raíz.

La fórmula de Newton-Raphson se deduce a partir de la fórmula de la


pendiente de una recta.

Pendiente de una recta:


Hay que determinar un número máximo de iteraciones

Normalmente esto se hace considerando una “tolerancia”, esto es:

El valor absoluto de la diferencia de la debe ser menor que la


tolerancia o el resultado de alguna fórmula de error debe ser menor que
la tolerancia dada.

Una de las fórmulas de error más útiles es la del error relativo porcentual
aproximado:

100 %

El método de Newton-Raphson es convergente cuadráticamente, es decir,


el error es aproximadamente al cuadrado del error anterior.

Esto significa que el número de cifras decimales correctas se duplica


aproximadamente en cada interacción.

Cuando el método de Newton-Raphson converge, se obtienen resultados


en relativamente pocas interacciones, ya que para raíces no repetidas este
método converge con orden 2 y el error Ei+1 es proporcional al cuadrado
del resultado anterior Ei

Supóngase que el error en una iteración es 10-n el error en la siguiente,


(que es proporcional al cuadrado del error anterior) es entonces
aproximadamente 10-2n, el que sigue será aproximadamente 10-4n etc.

De esto puede afirmarse que de cada iteración duplica aproximadamente


el número de dígitos correctos.

Sin embargo el método de Newton-Raphson algunas veces no converge,


sino que oscila. Esto ocurre si no hay raíz real, si la raíz es un punto de
inflexión o si el valor inicial está muy alejado de la raíz buscada y alguna
otra parte de la función “atrapa” la iteración.

Este método, el cual es un método iterativo, es uno de los más usados y


efectivos. A diferencia de los métodos anteriores, el método de Newton-
Raphson no trabaja sobre un intervalo sino que basa su fórmula en un
proceso iterativo.
Supongamos que tenemos la aproximación a la raíz de ,

Trazamos la recta tangente a la curva en el punto ; ésta cruza al


eje en un punto que será nuestra siguiente aproximación a la
raíz .
Para calcular el punto , calculamos primero la ecuación de la recta
tangente. Sabemos que tiene pendiente
Y por lo tanto la ecuación de la recta tangente es:

Hacemos :

Y despejamos :

Que es la fórmula iterativa de Newton-Raphson para calcular la siguiente


aproximación:

, si

Note que el método de Newton-Raphson no trabaja con intervalos donde


nos asegure que encontraremos la raíz, y de hecho no tenemos ninguna
garantía de que nos aproximaremos a dicha raíz. Desde luego, existen
ejemplos donde este método no converge a la raíz, en cuyo caso se dice
que el método diverge. Sin embargo, en los casos donde si converge a la
raíz lo hace con una rapidez impresionante, por lo cual es uno de los
métodos preferidos por excelencia.
También observe que en el caso de que , el método no se puede
aplicar. De hecho, vemos geométricamente que esto significa que la recta
tangente es horizontal y por lo tanto no intersecta al eje en ningún
punto, a menos que coincida con éste, en cuyo caso mismo es una raíz
de !
Ejemplo1
Usar el método de Newton-Raphson, para aproximar la raíz
de , comenzando con y hasta que .
Solución
En este caso, tenemos que
De aquí tenemos que:

Comenzamos con y obtenemos:

En este caso, el error aproximado es,

Continuamos el proceso hasta reducir el error aproximado hasta donde se


pidió.
Resumimos los resultados en la siguiente tabla:

Aprox. a la raíz Error aprox.


1
1.268941421 21.19%
1.309108403 3.06%
1.309799389 0.052%

De lo cual concluimos que , la cual es correcta en todos


sus dígitos!
La misma idea puede aplicarse para crear algoritmos que aproximen
raíces -ésimas de números reales positivos.
Observe que cuando el método de Newton-Raphson converge a la raíz, lo
hace de una forma muy rápida y de hecho, observamos que el error
aproximado disminuye a pasos agigantados en cada paso del proceso.
Aunque no es nuestro objetivo establecer formalmente las cotas para los
errores en cada uno de los métodos que hemos estudiado, cabe
mencionar que si existen estas cotas que miden con mayor precisión la
rapidez o lentitud del método en estudio.

Veremos a continuación un ejemplo del método de Newton Raphson, con


la siguiente ecuación:

# Fxn Dfxn Nuevo Xm


1 18 4 -3.5
2 -30.375 37.75 -2.6953642384106
3 -6.2771541041392 22.794965133108 -2.419989651633
4 -0.59229583988115 18.569049742033 -2.3880927130115
5 -0.0073539466744812 18.108960417816 -2.3876866186524
6 -1.1814129692311E-6 18.103142166676 -2.3876865533923

Hemos terminado de analizar el método de la Newton Raphson, en este


ejemplo con un error de 0.0001; se encuentra la última raíz (Xm): -
2.3876865533923 con 6 iteraciones.

ALGORITMO
Tres son las formas principales por las que tradicionalmente se ha
obtenido el algoritmo de Newton-Raphson.
La primera de ellas es una simple interpretación geométrica. En efecto,
atendiendo al desarrollo geométrico del método de la secante, podría
pensarse en que si los puntos de iteración están lo suficientemente cerca
(a una distancia infinitesimal), entonces la secante se sustituye por la
tangente a la curva en el punto. Así pues, si por un punto de iteración
trazamos la tangente a la curva, por extensión con el método de la
secante, el nuevo punto de iteración se tomará como la abscisa en el
origen de la tangente (punto de corte de la tangente con el eje X). Esto es
equivalente a linealizar la función, es decir, f se reemplaza por una recta
tal que contiene al punto ( , ( )) y cuya pendiente coincide con la
derivada de la función en el punto, . La nueva aproximación a la
raíz, , se logra la intersección de la función lineal con el eje X de
abscisas. Matemáticamente:

Una forma alternativa de obtener el algoritmo es desarrollando la función


f(x) en serie de Taylor, para un entorno del punto : En la ilustración
adjunta del método de Newton se puede ver que es una mejor
aproximación que para el cero (x) de la función f.

Si se trunca el desarrollo a partir del término de grado 2, y


evaluamos en :

Si además se acepta que tiende a la raíz, se ha de cumplir que


, luego, sustituyendo en la expresión anterior, obtenemos el
algoritmo.
Finalmente, hay que indicar que el método de Newton-Raphson puede
interpretarse como un método de iteración de punto fijo. Así, dada la
ecuación , se puede considerar el siguiente método de iteración
de punto fijo:

Se escoge h (x) de manera que g'(r)=0 (r es la raíz buscada).


Dado que g'(r) es:

Entonces:

Como h (x) no tiene que ser única, se escoge de la


forma más sencilla:
Por tanto, imponiendo subíndices:

Expresión que coincide con la del algoritmo de Newton-Raphson

Teorema de Convergencia Global del Método de Newton

Sea verificando1 :

para todo

para todo

Entonces existe un único tal que por lo que la sucesión


converge a s.
ESTIMACIÓN DEL ERROR
Se puede demostrar que el método de Newton-Raphson tiene
convergencia cuadrática: si es raíz, entonces:

para una cierta constante . Esto significa que si en algún momento el


error es menor o igual a 0,1, a cada nueva iteración doblamos
(aproximadamente) el número de decimales exactos. En la práctica puede
servir para hacer una estimación aproximada del error:
Error relativo entre dos aproximaciones sucesivas:

Con lo cual se toma el error relativo como si la última aproximación fuera


el valor exacto. Se detiene el proceso iterativo cuando este error relativo
es aproximadamente menor que una cantidad fijada previamente.

También podría gustarte