Scribd
Cargar
19 vistas3 páginas

Método Newton-Raphson Multivariable

El documento presenta el problema de encontrar las raíces del sistema de ecuaciones no lineales mediante el método de Newton-Raphson multivariable. Se evalúan las ecuaciones y sus derivadas parciales para el punto inicial (-1,3) y se realizan 10 iteraciones para encontrar una solución cercana a (0.2248, 5.7576) donde las funciones son menores a 0.0001.

Cargado por

Abel Ramos
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd
19 vistas3 páginas

Método Newton-Raphson Multivariable

El documento presenta el problema de encontrar las raíces del sistema de ecuaciones no lineales mediante el método de Newton-Raphson multivariable. Se evalúan las ecuaciones y sus derivadas parciales para el punto inicial (-1,3) y se realizan 10 iteraciones para encontrar una solución cercana a (0.2248, 5.7576) donde las funciones son menores a 0.0001.

Cargado por

Abel Ramos
Derechos de autor
© All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como DOCX, PDF, TXT o lee en línea desde Scribd

Newton Raphson Multivariable

Realizar 10 iteraciones con el Método de Newton – Raphson Multivariable o


hasta que |fi| ≤ 0.0001 para resolver el siguiente sistema de ecuaciones no
lineales:
x(0) = (-1,3)

x
f 1 ( x , y ) =7
( )
ln ( y )
−2 x 2 +1=0

f 2 ( x , y ) =9 x 2 y−5 x 2 y 2 + y =0

Derivadas parciales
∂f 1 7
= −4 x
∂ x ln ( y )
∂f 1 −7 x
=
∂ y y ( ln ( y ) )2

∂f 2
=18 xy−10 x y 2
∂x
∂f 2
=9 x 2−10 x2 y +1
∂y
Evaluando f1 y f2 con los valores iniciales

f 1 ( x , y ) =7
( ln−1( 3) )−2(−1) +1=−7.3716
2

f 2 ( x , y ) =9 (−1 )2 ( 3 )−5 (−1 )2 ( 3 )2 + ( 3 ) =−15

Evaluamos con las derivadas parciales


∂f 1 7
= −4 (−1 ) =10.3716
∂ x ln (3 )

∂ f 1 −7(−1)
= =1.9332
∂ y 3 ( ln ( 3 ) )2

∂f 2
=18 (−1 )( 3 )−10 (−1 )( 3 )2 =36
∂x
∂f 2
=9(−1)2−10 (−1 )2 (3 )+ 1=−20
∂y

El jacobiano esta dado por

∂f 1 ∂f 1
10.3716 1.9332
∂x ∂y
J= =
-20
∂f 2 ∂f 2
36
∂x ∂y

J=−207.432−(69.5952)=−277.0272
Cálculo de ∆ x y ∆ y

∆ x=
( f 2) ( ∂∂fy1 )−( f 1) ( ∂∂f y2 ) = (−15) ( 1.9332)−(−7.3716)(−20) =0.6 368
J −277.0272

∆ y=
( f 1) ( ∂∂fx2 )−( f 2 )( ∂∂fx1 ) = (−7.3716 )( 36) −(−15)(10.3716) =0. 3963
J −277. 0272

Nuevos valores
(1)
x =−1+ 0.6368=−0. 3632
(1)
y =3+0. 3963=3. 3963
Interac 0 1 2 3 4 5 6 7 8 9 10
x -1 -0.3631 -0.2428 -0.2346 -0.2248 -0.2248 -0.2248 -0.2248 -0.2248 -0.2248 -0.2248
y 3 3.3964 4.3526 5.6077 5.7468 5.7576 5.7576 5.7576 5.7576 5.7576 5.7576
f1 -7.3717 -1.3427 -0.2733 -0.0626 -0.0011 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
f2 -15.0000 -0.1784 1.0787 -0.2685 0.0143 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
 ∂f1/ ∂x 10.3717 7.1775 5.7304 4.9984 4.9024 4.8980 4.8980 4.8980 4.8980 4.8980 4.8980
 ∂f1/ ∂y 1.9332 0.5006 0.1805 0.0985 0.0896 0.0892 0.0892 0.0892 0.0892 0.0892 0.0892
 ∂f2/ ∂x 36 19.6884 26.9719 50.0944 50.9944 51.2246 51.2245 51.2245 51.2245 51.2245 51.2245
 ∂f2/ ∂y -20 -2.2918 -1.0347 -1.5911 -1.4499 -1.4548 -1.4548 -1.4548 -1.4548 -1.4548 -1.4548
J -277.0305 -26.3057 -10.7975 -12.8881 -11.6753 -11.6944 -11.6944 -11.6944 -11.6944 -11.6944 -11.6944
Δx 0.6369 0.1204 0.0082 0.0098 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000
Δy 0.3964 0.9562 1.2551 0.1391 0.0108 0.0000 0.0000 0.0000 0.0000 0.0000 0.0000

También podría gustarte