Método de Newton-Raphson

Método de Newton-Raphson
Ejemplo:
Sea f(x) una función polinomial, determine una aproximación para una raíz real mediante el método
de Newton-Raphson.
1. Obtenga una expresión para la derivada de la función f(x).
2. Suponga un valor inicial de la variable, xk, y defina una magnitud del error ε > 0.
Planteando como una ecuación polinomial
2x4 - 3x3 + 2x2 - 5x – 5 = 0
Iniciando con x0 = 0
Y hasta que l xk+1 – xk l < ε = 0.01
Solución:
Primera iteración
Dada la función polinomial f(x) = 2x4 - 3x3 + 2x2 - 5x - 5
Su derivada f' (x) = 8x3 - 9x2 + 4x - 5
El valor inicial de x (k = 0) x0 = 0
Y el tamaño del error ε = 0.01
3. Evalué la función f(x) y su derivada f' (x) en el punto x = xk
Como se indica en la secuencia de cálculo, se debe evaluar la función y su derivada en x0:
f(x0 = 0) = 2(0)4 - 3(0)3 + 2(0)2 - 5(0) – 5 = -5
f' (x0 = 0) = 8(0)3 - 9(0)2 + 4(0) -5 = -5
4. Compare el valor de la función con el error definido en el paso 2.
4.1. Si l f( xk ) l < ε, la raíz aproximada es xk
4.2. Si l f( xk ) l ≥ ε, continúe en el paso 5.
l f(0) = -5 l ≥ ε
5. Observe el valor de la derivada.
5.1. Si f' (xk) = 0, el método no converge en ese punto. Intente otro valor inicial xk
5.2. Si f' (xk) ≠ 0, continué con el paso 6.
f' (0) = -5 es ≠ 0, se procede con el paso siguiente.

𝒇(xk )
6. Calcule la siguiente aproximación xk+1 con la relación xk+1= xk – 𝒇'(x
k)
𝑓(x0 ) −5
x1 = x0 – 𝑓'(x =0- =-1
0) −5
7. Criterio de convergencia: compare el valor de xk+1 con el de xk
7.1. Si l xk+1 – xk l ≥ ε, siga al paso 8.
7.2. Si l xk+1 – xk l < ε, la raíz es xk+1
De donde l x1 - x0 l = l - 1 - 0 l = 1 > ε
8. Asigne al valor supuesto xk el valor calculado xk+1
Como xk+1 = -1 entonces el valor de xk será xk = -1 para la segunda iteración.
9. Vuelva al paso 3.
Ya se ha calculado un valor aproximado a la raíz en esta primera iteración, sin embargo, como se
han cumplido los criterios y aún es posible obtener un valor más cercano a la raíz podemos continuar
con la siguiente iteración.
Segunda iteración
f(-1) = 2(-1)4 - 3(-1)3 + 2(-1)2 - 5(-1) – 5 = 7
f'(-1) = 8(-1)3 - 9(-1)2 + 4(-1) – 5 = -26
Nuevamente, como l f(-1) l ≥ ε y f'(-1) ≠ 0, se procede a obtener el valor de x2

7
x2 = - 1 - = -0.7308
(−26)
En donde l x2 - x1 l = l - 0.7308 - (-1) l = 0.2692 > ε.
De acuerdo con la información que nos brinda la segunda repetición del método, es posible
calcular una razón aún más cercana.
Tercera iteración
f(-0.7308) = 2(-0.7308)4 - 3(-0.7308)3 + 2(-0.7308)2 - 5(-0.7308) – 5 = 1.4630
f'(-0.7308) = 8(-0.7308)3 - 9(-0.7308)2 + 4(-0.7308) – 5 = -15.8513
Con lo anterior l f(-0.7308) l ≥ ε y f'(-0.7308) ≠ 0, se procede a obtener el valor de x3

1.4630
x3 = -0.7308 - (−15.8513) = -0.6385
En donde l x3 – x2 l = l -0.6385 - (-0.7308) l = 0.0923 > ε.
Como el valor absoluto de la diferencia de la última anterior con su antecesora sigue siendo mayor
a la épsilon dada, aún se puede calcular otra aproximación.
Cuarta iteración
f(-0.6385) = 2(-0.6385)4 - 3(-0.6385)3 + 2(-0.6385)2 - 5(-0.6385) – 5 = 0.1208

f'(-0.6385) = 8(-0.6385)3 - 9(-0.6385)2 + 4(-0.6385) – 5 = -13.3049
Nuevamente, como l f((-0.6385) l ≥ ε y f'((-0.6385) ≠ 0, se procede a obtener el valor de x3

0.1208
x4 = -0.6385- (−13.3049) = -0.6294
En donde l x4 – x3 l = l -0.6294- (-0.6385) l = 0.0091 < ε.
Se recomienda la construcción de una tabla, tanto para simplificar el procedimiento, como para
observar la disminución del error en cada iteración, además que en algunos polinomios basta con
pocas iteraciones para obtener una raíz.
El método de Newton-Raphson converge a la solución de una manera más rápida con respecto a
otros métodos. El tamaño del error se reduce con pocas iteraciones.
Por lo tanto, el resultado final que se obtiene es:
x4 = - 0.6294, con un error aproximado de 0.0091
Ejercicio:
Use el método de Newton-Raphson para aproximar la raíz de la función trascendente
f(x) = x2 – 5x – ex con l xk+1 – xk l < ε = 0.01

Método de Newton-Raphson

Cargado por

Copyright:

Formatos disponibles

Método de Newton-Raphson

Cargado por

Información del documento

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Método de Newton-Raphson

Cargado por

Copyright:

Formatos disponibles

Método de Newton-Raphson

1. Obtenga una expresión para la derivada de la función f(x).

Planteando como una ecuación polinomial

2x4 - 3x3 + 2x2 - 5x – 5 = 0

Y hasta que l xk+1 – xk l < ε = 0.01

Dada la función polinomial f(x) = 2x4 - 3x3 + 2x2 - 5x - 5

Su derivada f' (x) = 8x3 - 9x2 + 4x - 5

Y el tamaño del error ε = 0.01

3. Evalué la función f(x) y su derivada f' (x) en el punto x = xk

Como se indica en la secuencia de cálculo, se debe evaluar la función y su derivada en x0:

f(x0 = 0) = 2(0)4 - 3(0)3 + 2(0)2 - 5(0) – 5 = -5

f' (x0 = 0) = 8(0)3 - 9(0)2 + 4(0) -5 = -5

4. Compare el valor de la función con el error definido en el paso 2.

4.1. Si l f( xk ) l < ε, la raíz aproximada es xk

4.2. Si l f( xk ) l ≥ ε, continúe en el paso 5.

5. Observe el valor de la derivada.

5.2. Si f' (xk) ≠ 0, continué con el paso 6.

f' (0) = -5 es ≠ 0, se procede con el paso siguiente.

7. Criterio de convergencia: compare el valor de xk+1 con el de xk

7.1. Si l xk+1 – xk l ≥ ε, siga al paso 8.

7.2. Si l xk+1 – xk l < ε, la raíz es xk+1

8. Asigne al valor supuesto xk el valor calculado xk+1

Como xk+1 = -1 entonces el valor de xk será xk = -1 para la segunda iteración.

f(-1) = 2(-1)4 - 3(-1)3 + 2(-1)2 - 5(-1) – 5 = 7

f'(-1) = 8(-1)3 - 9(-1)2 + 4(-1) – 5 = -26

Nuevamente, como l f(-1) l ≥ ε y f'(-1) ≠ 0, se procede a obtener el valor de x2

En donde l x2 - x1 l = l - 0.7308 - (-1) l = 0.2692 > ε.

f(-0.7308) = 2(-0.7308)4 - 3(-0.7308)3 + 2(-0.7308)2 - 5(-0.7308) – 5 = 1.4630

f'(-0.7308) = 8(-0.7308)3 - 9(-0.7308)2 + 4(-0.7308) – 5 = -15.8513

Con lo anterior l f(-0.7308) l ≥ ε y f'(-0.7308) ≠ 0, se procede a obtener el valor de x3

En donde l x3 – x2 l = l -0.6385 - (-0.7308) l = 0.0923 > ε.

f(-0.6385) = 2(-0.6385)4 - 3(-0.6385)3 + 2(-0.6385)2 - 5(-0.6385) – 5 = 0.1208

Nuevamente, como l f((-0.6385) l ≥ ε y f'((-0.6385) ≠ 0, se procede a obtener el valor de x3

En donde l x4 – x3 l = l -0.6294- (-0.6385) l = 0.0091 < ε.

Por lo tanto, el resultado final que se obtiene es:

x4 = - 0.6294, con un error aproximado de 0.0091

Use el método de Newton-Raphson para aproximar la raíz de la función trascendente

f(x) = x2 – 5x – ex con l xk+1 – xk l < ε = 0.01

También podría gustarte