Clase 4 2022 2 Newton Raphson
Clase 4 2022 2 Newton Raphson
Clase 4 2022 2 Newton Raphson
Métodos Numéricos
2
𝑥𝑥 − 𝑥𝑥
𝑓𝑓 𝑥𝑥 = 𝑓𝑓 𝑥𝑥 + 𝑥𝑥 − 𝑥𝑥 𝑓𝑓 ′ 𝑥𝑥 + 𝑓𝑓 ′′ 𝜉𝜉 𝑥𝑥
2!
ó
2
1
𝑥𝑥 − 𝑥𝑥 2
𝑓𝑓 𝑥𝑥 = 𝑓𝑓 𝑥𝑥 + 𝑥𝑥 − 𝑥𝑥 𝑓𝑓 𝑥𝑥 + 𝑓𝑓 𝜉𝜉 𝑥𝑥
2!
𝑓𝑓 𝑥𝑥
− 1
≅ 𝑝𝑝 − 𝑥𝑥
𝑓𝑓 𝑥𝑥
𝑓𝑓 𝑥𝑥 𝑓𝑓 𝑝𝑝 𝑛𝑛−1
∴ 𝑝𝑝 ≅ 𝑥𝑥 − 1
⇒ 𝑝𝑝 𝑛𝑛 ≅ 𝑝𝑝 𝑛𝑛−1 − 1 , 𝑛𝑛 ≥ 1
𝑓𝑓 𝑥𝑥 𝑓𝑓 𝑝𝑝 𝑛𝑛−1
Tenemos que:
𝒇𝒇 𝒙𝒙𝒏𝒏
𝒙𝒙 𝒏𝒏+𝟏𝟏 = 𝒙𝒙 𝒏𝒏 − ′ , 𝒏𝒏 ≥ 𝟎𝟎
𝒇𝒇 𝒙𝒙𝒏𝒏
Nota
1) 𝑓𝑓 ′ 𝑥𝑥 ≠ 0 , pues la recta tangente que pasa por el punto (𝑥𝑥0 , 𝑓𝑓(𝑥𝑥0 )) , no debe ser paralelo al eje de
las X.
2) 𝐶𝐶 2 [𝑎𝑎, 𝑏𝑏].= Conjunto de todas las funciones que tienen segunda derivada, continua en x
ALGORITMO DE NEWTON RAPHSON
𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 1
𝑥𝑥1 = 6.099630
S𝑒𝑒𝑒𝑒 𝑛𝑛 = 2
𝑥𝑥2 = 6.093552
𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 3
𝑥𝑥3 = 6.093529
n 𝒙𝒙𝒏𝒏 𝝃𝝃𝒔𝒔 = ⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |
0 𝑥𝑥0 = 6
1 6.099630 0.099630
2 6.093552 0.006078
3 6.093529 0.000023
EJEMPLO
(4 dígitos)
Sea
𝑥𝑥0 = 23
𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 1
𝑥𝑥1 = 24.8567
INTRODUCCIÓN
S𝑒𝑒𝑒𝑒 𝑛𝑛 = 2
𝑥𝑥2 = 24.3840
𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 3
𝑥𝑥3 = 24.5206
n 𝒙𝒙𝒏𝒏 𝝃𝝃𝒔𝒔 = ⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |
0 𝑥𝑥0 = 23
1 24.8567 1.8567
2 24.3840 0.4726
3 24.5206 0.1366
4 24.5307 0.0101
5 24.5308 0.0001
6 24.5308 0.0000
EJEMPLO:
https://picazo.eltiempo.es/2015-el-ano-en-el-que-superamos-las-400-ppm-de-c02/
𝑥𝑥
𝑓𝑓 𝑥𝑥 = 3𝑥𝑥−5 − 𝑐𝑐𝑐𝑐𝑐𝑐 −6
2
′ 𝑥𝑥−5
1 𝑥𝑥
𝑓𝑓 𝑥𝑥 = 3 𝑙𝑙𝑙𝑙 3 + 𝑠𝑠𝑠𝑠𝑠𝑠
2 2
𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 1
𝑥𝑥1 = 7.260567
S𝑒𝑒𝑒𝑒 𝑛𝑛 = 2
𝑥𝑥2 = 6.729556
𝑆𝑆𝑆𝑆𝑆𝑆 𝑛𝑛 = 3
𝑥𝑥3 = 6.499873
n 𝒙𝒙𝒏𝒏 𝐸𝐸𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎𝑎
0 𝑥𝑥0 = 8
1 7.260567 0.101842
2 6.729556 0.078907
3 6.499873 0.035337
4 6.466344 0.005185
5 6.465731 0.000095 ≤ 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎
𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 − 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨
𝑬𝑬𝒂𝒂𝒂𝒂𝒂𝒂𝒂𝒂𝒂𝒂 = 𝑬𝑬𝑬𝑬𝑬𝑬𝑬𝑬𝑬𝑬 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 = ⃒ ⃒
𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨 𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨𝑨
EJEMPLO:
Tenemos: x0 = 0.5
0 0.5000
1 0.3211
2 0.3168
3 0.3168
4 0.3168
5 0.3168
6 0.3168
n Xn
f(x)= 0 0.5000 -0.0109
f'(X)= -0.0607
1 0.3211 -0.0003
-0.0579
2 0.3168 0.0000
-0.0578
3 0.3168 0.0000
-0.0578
4 0.3168 0.0000
-0.0578
5 0.3168 0.0000
-0.0578
6 0.3168 0.0000
-0.0578
EJEMPLO:
2
Resolver 𝑓𝑓 𝑥𝑥 = 𝑥𝑥 − 2 − 𝑙𝑙𝑙𝑙 𝑥𝑥 utilizando el método de
Tenemos:
⇒ 𝒙𝒙𝟏𝟏 = 1.33333333
0 1.00000000
1 1.33333333 0.33333333
2 1.40857927 0.07524594
3 1.41238156 0.00380229
4 1.41239117 0.00000961
5 1.41239117 0.00000000
n Xi
f (x) = 0 1.00000000 1.00000000
f '(x)= -3.00000000
1 1.33333333 0.15676237
-2.08333333
2 1.40857927 0.00719689
-1.89277665
3 1.41238156 0.00001809
-1.88326084
4 1.41239117 0.00000000
-1.88323681
5 1.41239117 -6.17219078
-1.88323681
EJEMPLO:
Tenemos:
0 𝑥𝑥0 = 5
1 3.424657534 1.575342466
2 2.924356997 0.500300537
3 2.861146976 0.063210021
4 2.860104689 0.001042287
5 2.860104406 0.000000283 < 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
EJEMPLO:
𝑓𝑓 1 𝑓𝑓 3 < 0 ⇒ ∃ 𝑥𝑥 ∗ ∈ 1, 3
𝑓𝑓 1 = 0.841470985 𝑓𝑓 3 = −1.858879992
𝑓𝑓 ′′ 1 = −0.841470985 𝑓𝑓 ′′ 3 = −0.141120008
Tenemos:
𝑠𝑠𝑠𝑠𝑠𝑠 3−3+1
⇒ 𝑥𝑥1 = 3 −
cos 3−1
⇒ 𝑥𝑥1 = 2.065885929
n Xn
0
3.000000000
1
2.065885929
2
1.939820941
3
1.934572691
4
1.934563211
5
1.934563211
OBSERVACION
3) 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 ó está el valor inicial muy alejado de la
raíz buscada ó se entrapa.
𝑓𝑓 𝑥𝑥 = 𝑥𝑥 3 − 5𝑥𝑥 2 + 7𝑥𝑥 − 3
𝑓𝑓 𝑥𝑥 = 𝑥𝑥 − 3 𝑥𝑥 − 1 𝑥𝑥 − 1
Ejemplo 2.-
𝑢𝑢 𝑥𝑥𝑖𝑖
𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 − 2
𝑢𝑢′ 𝑥𝑥𝑖𝑖
2
𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓′ 𝑥𝑥𝑖𝑖
𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 −
𝑓𝑓′ 𝑥𝑥𝑖𝑖 𝑓𝑓′ 𝑥𝑥𝑖𝑖 2 − 𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓′′ 𝑥𝑥𝑖𝑖 ]
Por lo tanto
𝒇𝒇 𝒙𝒙𝒊𝒊 𝒇𝒇′ 𝒙𝒙𝒊𝒊
𝒙𝒙 𝒊𝒊+𝟏𝟏 = 𝒙𝒙 𝒊𝒊 −
𝒇𝒇′ 𝒙𝒙𝒊𝒊 𝟐𝟐 − 𝒇𝒇 𝒙𝒙𝒊𝒊 𝒇𝒇′′ 𝒙𝒙𝒊𝒊
EJEMPLO:
𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′ 𝑥𝑥𝑖𝑖
𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 − ′
𝑓𝑓 𝑥𝑥𝑖𝑖 2 − 𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′′ 𝑥𝑥𝑖𝑖
𝟏𝟏𝟏𝟏−𝟓𝟓 = 0.0001
n 𝒙𝒙𝒏𝒏
⃒𝒙𝒙𝒏𝒏+𝟏𝟏 − 𝒙𝒙𝒏𝒏 |
0 0.500000000
1 0.568013734
0.068013734
2 0.567380913
0.000632821
3 0.567208171
0.000172742
4 0.567129245
0.000078926
5 0.567146333
0.000017087
6 0.567142632 0.000003701 < 𝟏𝟏𝟏𝟏−𝟓𝟓 = 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎
EJEMPLO:
𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′ 𝑥𝑥𝑖𝑖
𝑥𝑥 𝑖𝑖+1 = 𝑥𝑥 𝑖𝑖 − ′
𝑓𝑓 𝑥𝑥𝑖𝑖 2 − 𝑓𝑓 𝑥𝑥𝑖𝑖 𝑓𝑓 ′′ 𝑥𝑥𝑖𝑖
𝑥𝑥0 =1.3
PRIMERA ITERACION
𝑓𝑓(𝑥𝑥0 ) 𝑓𝑓𝑓(𝑥𝑥0 )
𝑥𝑥1 = 𝑥𝑥0 − = 1.003083062
𝑓𝑓𝑓(𝑥𝑥0 ) 2 − 𝑓𝑓(𝑥𝑥0 ) 𝑓𝑓𝑓𝑓(𝑥𝑥0 )
𝒙𝒙𝟏𝟏 = 1.003083062
𝑓𝑓(𝑥𝑥1 ) 𝑓𝑓𝑓(𝑥𝑥1 )
𝑥𝑥2 = 𝑥𝑥1 − = 1.042758770
𝑓𝑓𝑓(𝑥𝑥1 ) 2 − 𝑓𝑓(𝑥𝑥1 ) 𝑓𝑓𝑓𝑓(𝑥𝑥1 )
𝒙𝒙𝟐𝟐 = 1.042758770
2
𝑓𝑓 ′(1.042758770) = 9 1.003083062 + 10 1.003083062 − 𝑐𝑐𝑐𝑐𝑐𝑐 1.003083062 = 19.70986122
𝑓𝑓(𝑥𝑥2 ) 𝑓𝑓′(𝑥𝑥2 )
𝑥𝑥3 = 𝑥𝑥2 − = 1.044052552
𝑓𝑓′(𝑥𝑥2 ) 2 − 𝑓𝑓(𝑥𝑥2 ) 𝑓𝑓′′(𝑥𝑥2 )
𝒙𝒙𝟑𝟑 = 1.044052552
0 𝑥𝑥0 = 𝜋𝜋
1 3.490658504 0349065850
2 3.476054580 0.014603924
La velocidad de producción de gas cloro en una celda electrolítica viene dada por
la siguiente función :
𝑥𝑥
𝑓𝑓 𝑥𝑥 = 2𝑐𝑐𝑐𝑐𝑐𝑐 𝑥𝑥 − 𝑡𝑡𝑡𝑡𝑡𝑡( ) donde utiliza Cloruro de Sodio (NaCl) acuoso
2
′′
1 𝑥𝑥 2 𝑥𝑥
𝑓𝑓 𝑥𝑥 = −2 cos 𝑥𝑥 − sec 𝑡𝑡𝑡𝑡𝑡𝑡
2 2 2
𝑥𝑥0 =0.008727
PRIMERA ITERACION
0.008727
𝑓𝑓 0.008727 = 2𝑐𝑐𝑐𝑐𝑐𝑐 0.008727 − 𝑡𝑡𝑡𝑡𝑡𝑡 =1.9956
2
2
1 0.008727
𝑓𝑓 ′ 0.008727 = −2𝑠𝑠𝑠𝑠𝑠𝑠 0.008727 − sec = −0.5175
2 2
2
′′
1 0.008727 0.008727
𝑓𝑓 0.008727 = −2 cos 0.008727 − sec 𝑡𝑡𝑡𝑡𝑡𝑡 = −2.0021
2 2 2
𝑓𝑓(𝑥𝑥0 ) 𝑓𝑓𝑓(𝑥𝑥0 )
𝑥𝑥1 = 𝑥𝑥0 − = 0.2510
𝑓𝑓𝑓(𝑥𝑥0 ) 2 − 𝑓𝑓(𝑥𝑥0 ) 𝑓𝑓𝑓𝑓(𝑥𝑥0 )
𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
𝑓𝑓 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = 2𝑐𝑐𝑐𝑐𝑐𝑐 𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎𝟎 − 𝑡𝑡𝑡𝑡𝑡𝑡 = 1.8112
2
2
1 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
𝑓𝑓 ′(𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐) = −2𝑠𝑠𝑠𝑠𝑠𝑠 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 − 𝑠𝑠𝑠𝑠𝑠𝑠 = −1.0046
2 2
1 𝟎𝟎.𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 2 𝟎𝟎.𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐
′′
𝑓𝑓 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 = −2 𝑐𝑐𝑐𝑐𝑐𝑐 𝟎𝟎. 𝟐𝟐𝟐𝟐𝟐𝟐𝟐𝟐 − 𝑠𝑠𝑠𝑠𝑠𝑠 𝑡𝑡𝑡𝑡𝑡𝑡 = −2.0014
2 2 2
𝑓𝑓(𝑥𝑥1 ) 𝑓𝑓𝑓(𝑥𝑥1 )
𝑥𝑥2 = 𝑥𝑥1 − = 0.6437
𝑓𝑓𝑓(𝑥𝑥1 ) 2 − 𝑓𝑓(𝑥𝑥1 ) 𝑓𝑓𝑓𝑓(𝑥𝑥1 )
𝒙𝒙𝟐𝟐 = 0.6436
0.6437
𝑓𝑓 0.6437 = 2𝑐𝑐𝑐𝑐𝑐𝑐 0.6437 − 𝑡𝑡𝑡𝑡𝑡𝑡 = 1.2665
2
2
1 0.6437
𝑓𝑓 ′(0.6437) = −2𝑠𝑠𝑠𝑠𝑠𝑠 0.6437 − sec = −1.7557
2 2
2
′′
1 0.6437 0.6437
𝑓𝑓 0.6437 = −2 cos 0.6437 − sec 𝑡𝑡𝑡𝑡𝑡𝑡 = −1.7851
2 2 2
𝑓𝑓(𝑥𝑥2 ) 𝑓𝑓′(𝑥𝑥2 )
𝑥𝑥3 = 𝑥𝑥2 − 2
= 0.3461
𝑓𝑓′(𝑥𝑥2 ) − 𝑓𝑓(𝑥𝑥2 ) 𝑓𝑓′′(𝑥𝑥2 )
𝒙𝒙𝟑𝟑 = 1.0597
0 0.008727
1 0.2510 0.2422
2 0.6436 0.3926
3 1.0597 0.4161
4 1.2085 0.1487
5 1.2155 0.0071
6 1.2156 0 < 𝟎𝟎. 𝟎𝟎𝟎𝟎𝟎𝟎
EJEMPLO:
Para un proceso minero metalúrgico se evalúa la etapa de conminución,
empleando el Método de Newton Raphson Modificado en la función
𝒇𝒇 𝒙𝒙 = 𝒙𝒙𝟑𝟑 − 𝟓𝟓𝒙𝒙𝟐𝟐 + 𝟕𝟕𝒙𝒙 − 𝟑𝟑, para ello se desea saber sobre el rendimiento de las
chancadoras primarias en función de su capacidad esperando obtener un
resultado con error relativo porcentual pequeño, con un valor inicial de 𝒙𝒙𝟎𝟎 = 𝟎𝟎 y el
valor verdadero es 1. (6 dígitos)
0 𝑥𝑥0 = 0 100
1 1.105263158 10.5263158
2 1.003081664 0.3081664
3 1.000002402 0.0002395
n 𝒙𝒙𝒏𝒏 𝑬𝑬𝑹𝑹%
0 𝒙𝒙𝟎𝟎 = 0
1 0.42857 57.143
2 0.68571 31.429
3 0.83286 16.714
4 0.91332 8.668
5 0.95577 4.423
6 0.97764 2.236
7 0.98875 1.125