0% encontró este documento útil (0 votos)
207 vistas12 páginas

Método de Newton Raphson

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1/ 12

INTRODUCCIN

En el presente escrito, se pretende abordar una de las diversas alternativas


existentes para la resolucin de ecuaciones, segn menciona Nieves Hurtado y
Domnguez Snchez (2006) de la forma f(x) = 0, donde f(x) es una funcin real de
una variable x, como un polinomio en x:
5

f ( x )=4 x + x 8 x +2
o una funcin transcendente
f ( x )=e x senx + ln 3 x + x 3

En general, se conoce que no es posible determinar los ceros de una funcin, es


decir, valores x* tal que f (x*) = 0, en un nmero finito de pasos.
En casos como este, tenemos que usar mtodos de aproximacin los
cuales son usualmente iterativos y tienen la forma: iniciando con una
aproximacin inicial x0, se calculan aproximaciones sucesivas x 1,x2 y
elegimos xn como aproximacin de x* cuando se cumpla un criterio de
parada dado. A los ceros de un polinomio se les conoce tambin como
races.
Los mtodos iterativos ms usuales son: biseccin, regula falsi, punto fijo,
Newton y el mtodo de la secante. En este caso, nos enfocaremos en
analizar el mtodo de Newton Raphson, sus aplicaciones, as como
observar las ventajas y desventajas que podemos hallar al hacer uso de l.
Este mtodo, tambin llamado Mtodo de Newton, es uno de los mtodos
que muestra mejor velocidad de convergencia llegando (bajo ciertas
condiciones) a duplicar, en cada iteracin, los decimales exactos.
En la literatura inglesa se acostumbra llamar al mtodo de Newton, mtodo
de Newton-Raphson.
El mtodo de Newton aparece en su libro Method of Fluxions escrito en
1671 pero publicado hasta 1736. El mtodo fue publicado por primera vez
en un libro de J. Raphson en 1690.
Se dice que Raphson tuvo acceso al manuscrito de Newton. En todo caso
ni Newton ni Raphson mencionan las derivadas. El primero que dio una
descripcin del mtodo de Newton usando derivadas fue T. Simpson en
1740.
(Mora F. 2016, p.108)

FUNDAMENTOS TERICOS
MTODO DE NEWTON RAPHSON
Este es un mtodo de segundo orden de convergencia cuando se trata de races
reales no repetidas. Consiste en un procedimiento que lleva la ecuacin f ( x )=0
a la forma

x=g(x ) , de modo que

g' ( x i ) =0 . Su deduccin se presenta en

seguida:
En la siguiente figura (Fig.1) se tiene la grfica
es una raz real

f ( x ) cuyo cruce con el eje

x .

Fig. 1

Vamos a suponer un valor inicial

x0

una tangente a la curva en el punto

que se sita en el eje horizontal. Si se traza


x
( 0 , f ( x0 ) ) y a partir de ese punto se sigue

x ; el punto de corte

por la tangente hasta su interseccin con el eje


una nueva aproximacin a
f ( x ) por su tangente en

as

x
( 0 , f ( x0 ) ) .

sucesivamente,

|f (x i)| 1 ,|x i +1x|<

es

x (hay que observar que se ha reemplazado la curva

El proceso se repite comenzando con


x2

x1

x1

hasta

, se obtiene una nueva aproximacin


que

un

x1

valor

satisfaga

o ambos. Si lo anterior no se cumpliera en un mximo de

iteraciones debe reiniciarse con un nuevo valor

x0

La ecuacin central del algoritmo se obtiene as:

La pendiente de la tangente a la curva en el punto

x
( 0 , f ( x0 ) ) es:

as que:

y sustituyendo:

en general:

(Ecuacin 1)

Se considera que este mtodo es de orden 2 porque

g' ( x )=0

g' ' ( x ) 0.

Ejemplo 1
Encuentre una raz real de la ecuacin
mtodo de Newton Raphson,

x 0=1

f ( x )=x 3 +2 x 2 +10 x20

, con =10-3 aplicado

|x i+1x|

mediante el
.

Solucin:
La primera derivada de la funcin es

f ' ( x )=3 xi2 + 4 x i+10

Se sustituyen f(x) y f(x) en la Ecuacin 1


3

x i+1=x i

xi +2 x i +10 x i20
2
3 x i + 4 x i +10

Primera iteracin:

( 1 )3+ 2 ( 1 )2 +10 ( 1 )20


x 1=1
=1.41176
3 ( 1 )2 +4 (1 )+10
Como

x 1 x0

, se calcula

x2

Segunda iteracin:
x 1=1.41176

( 1.41176 )3 +2 ( 1.41176 )2+ 10 (1.41176 )20


=1.36934
3 (1.41176 )2 +4 (1.41176 ) +10

Con este proceso se obtiene la siguiente tabla:


i

xi

0
1
2
3
4

1
1.41176
1.36934
1.36881
1.36881

|x i+1x i|

g ' ( x i )

0.41176
0.04243
0.00053
0.0000

0.24221
0.02446
0.00031
4.6774 X 10-8
9.992 X 10-16

Se requirieron slo tres iteraciones para satisfacer el criterio de convergencia.

(Nieves Hurtado y Domnguez Snchez, 2006, pp.44-45)


Ejemplo 2
Utilice el mtodo de Newton-Raphson para calcular la raz de f(x)=e xx
empleando como valor inicial x0 = 0.
Solucin:
La primera derivada de la funcin es:

f ' ( x )=e x 1

que se sustituye, junto con

la funcin original en nuestra frmula para tener:


xi

x i+1=x i

e x i
xi

e 1

Empezando con un valor inicial x 0 = 0, se aplica esta ecuacin iterativa para


calcular
i

xi

|x i+1x i|

0
1
2
3
4

0
0.5
0.566311003
0.567143165
0.567143290

0.5
0.066311
0.00083216
1.25E-07

Error relativo
porcentual
100.0000
11.7093
0.1467
0.0000

As, el mtodo converge rpidamente a la raz verdadera.


Criterio de terminacin y estimacin de errores
Como en los otros mtodos para localizar races, la ecuacin:

Se utiliza como un criterio de terminacin. No obstante, el desarrollo del mtodo


con base en la serie de Taylor (cuadro 6.2), proporciona una comprensin terica
respecto a la velocidad de convergencia expresada por Ei+1 = O(Ei2).
De esta forma, el error debe ser proporcional al cuadrado del error anterior. En
otras palabras, el nmero de cifras significativas de precisin aproximadamente se
duplica en cada iteracin. Dicho comportamiento se examina en el siguiente
ejemplo.

Ejemplo 3
Sabiendo que el mtodo de Newton-Raphson es convergente en forma cuadrtica.
Es decir, el error es proporcional al cuadrado del error anterior:
Ecuacin 2

Examine esta frmula y observe si concuerda con los resultados del ejemplo 2.
La primera derivada de de f(x)=exx

es

f ' ( x )=e x 1

que se evala en

xr=0.56714329 para dar (0.56714329) = 1.56714329.


''
x
La segunda derivada es: f ( x )=e
la cual se evala como (0.56714329)=
0.56714329.
Estos resultados se sustituyen en la Ecuacin 2
Et ,i+1

0.56714329
2
2
Et , i=0.18095 Et ,i
2 ( 1.56714329 )

En el ejemplo 2, el error inicial fue E t,1 0.56714329, el cual se sustituye en la


ecuacin de error que predice:
2

Et ,1 0.18095(0.56714329) =0.0582
que es cercano al error verdadero de 0.06714329. En la siguiente iteracin,
2

Et ,2 0.18095(0.06714329) =0.0008158
que tambin se compara de manera favorable con el error verdadero 0.0008323.
Para la tercera iteracin,
Et ,3 0.18095(0.0008323)2=0.000000125
que es el error obtenido en el ejemplo anterior (2). As, la estimacin del error
mejora, ya que conforme nos acercamos a la raz, x y se aproximan mejor
mediante xr. Finalmente:
Et ,4 0.18095(0.000000125)2=2.8310 x 1015
As, este ejemplo ilustra que el error en el mtodo de Newton-Raphson para este
caso es, de hecho, proporcional (por un factor de 0.18095) al cuadrado del error
en la iteracin anterior.

Desventajas del mtodo de Newton-Raphson


Aunque en general el mtodo de Newton-Raphson es muy eficiente, hay
situaciones donde se comporta de manera deficiente. Por ejemplo en el caso
especial de races mltiples que se analizar ms adelante en este captulo. Sin
embargo, tambin cuando se trata de races simples, se encuentran dificultades,
como en el siguiente ejemplo:
Ejemplo 4

Ejemplo de una funcin que converge lentamente con el mtodo de NewtonRaphson.


Determine la raz positiva de f(x) = x 101 usando el mtodo de Newton-Raphson y
un valor inicial x = 0.5.
Solucin:
La frmula de Newton-Raphson en este caso es:
x10
1
x i+1=x i 1 9
10 x i
que se utiliza para calcular:
Iteracin
0
1
2
3
4
5
.
.
.

x
0.5
51.65
46.485
41.8365
37.65285
33.887565

1.000000

De esta forma, despus de la primera prediccin deficiente, la tcnica converge a


la raz verdadera, 1, pero muy lentamente.
Adems de la convergencia lenta debido a la naturaleza de la funcin, es posible
que se presenten otras dificultades, como se ilustra en la figura 2 (fig2).
Por ejemplo el inciso a) muestra el caso donde un punto de inflexin [esto es,
(x)= 0] ocurre en la vecindad de una raz.
Se observa que las iteraciones que empiezan con x 0 divergen progresivamente de
la raz.
En el inciso b) se ilustra la tendencia del mtodo de Newton-Raphson a oscilar
alrededor de un mnimo o mximo local. Tales oscilaciones pueden persistir o
alcanzar una pendiente cercana a cero, despus de lo cual la solucin se aleja del
rea de inters. En el inciso c) se muestra cmo un valor inicial cercano a una raz
salta a una posicin varias races ms lejos. Esta tendencia a alejarse del rea de
inters se debe a que se encuentran pendientes cercanas a cero. En efecto, una
pendiente cero [(x) = 0] es un verdadero desastre, ya que causa una divisin
entre cero en la frmula de Newton-Raphson.
En forma grfica (inciso d), esto significa que la solucin se dispara
horizontalmente y jams toca al eje x.
De manera que no hay un criterio general de convergencia para el mtodo de
Newton- Raphson. Su convergencia depende de la naturaleza de la funcin y de la

exactitud del valor inicial. La nica solucin en estos casos es tener un valor inicial
que sea suficientemente cercano a la raz.
Cabe sealar, que para algunas funciones ningn valor inicial funcionar. Los
buenos valores iniciales por lo comn se predicen con un conocimiento del
problema fsico o mediante el uso de recursos alternativos, tales como las
grficas, que proporcionan mayor claridad en el comportamiento de la solucin.
Ante la falta de un criterio general de convergencia se sugiere el diseo de
programas computacionales eficientes que reconozcan la convergencia lenta o la
divergencia.
Fig2

(Chapra y Canale,2007, pp.149-152)

CONCLUSIONES
Habiendo analizado los fundamentos, as como aplicacin, ventajas y desventajas
del mtodo de Newton Raphson, podemos concluir con los ejemplos vistos, que
presenta bastante practicidad cuando se trata de funciones no muy complejas, y
que la rapidez con la que se puede llegar a un valor deseado nos da la ventaja de
realizar el menor nmero de pasos iterativos y de alcanzar la solucin pronta a
nuestro problema en cuestin.
Por otra parte, observamos que a pesar de su rapidez el mtodo de Newton
Raphson necesita el clculo de la derivada y se sabe que no todas las funciones

derivables tienen derivadas que se pueden expresar en trminos de funciones


elementales o funciones especiales y podra colapsar si la derivada toma valores
muy pequeos en el proceso.
Tambin en las funciones obtenidas por interpolacin corre el riesgo de caer en un
ciclo. Por las razones citadas anteriormente, el mtodo presenta cierto grado de
inexactitud debido a su rapidez para converger.
De esta manera, se puede afirmar, que se debe poner especial atencin a la
funcin con la que se trabaje puesto que adems del mtodo de Newton Raphson,
existen distintos algoritmos para encontrar las races de las funciones pero
ninguno, incluido Newton Raphson, es general, es decir, no hay un algoritmo que
funcione con todas las ecuaciones.

BIBLIOGRAFA
Nieves H.A. & Domnguez S. F. (2006) Mtodos numricos aplicados a la
ingeniera. 2Ed. Editorial CECSA. Mxico.
Mora F.W., (2016) Introduccin a los Mtodos Numricos. Costa Rica.
Recuperado de http://www.tec-digital.itcr.ac.cr/revistamatematica/
Chapa S., Canale R. (2007) Mtodos numricos para ingenieros. 5 Ed. Editorial
McGraw-Hill. Mxico

También podría gustarte

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy