Contenido Explicación teórica
Sistemas Eléctricos de Potencia
Formulación método
Newton-Raphson
Canal de YouTube: Aprende con Romero
M. C. David Romero Rojas
Lección: Flujos de Potencia
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
1 Explicación teórica
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
Explicación teórica: Método Newton-Raphson
Newton-Raphson es un método iterativo para la solución de
ecuaciones no lineales. Se formula a partir de dos funciones g1 &
g2 con dos variables x1 & x2 que se representan mediante:
g1 (x1 , x2 , u) = h1 (x1 , x2 , u) − b1 = 0 (1)
g2 (x1 , x2 , u) = h2 (x1 , x2 , u) − b2 = 0 (2)
donde u representa una variable de control (se puede considerar
como una constante) y b1 & b2 se consideran constantes
conocidas.
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
Explicación teórica: Método Newton-Raphson
Los valores iniciales del método iterativo se representan mediante
(0) (0) (0) (0)
x1 & x2 , que se suman a los incrementos ∆x1 & ∆x2 para
(∗) (∗)
obtener las soluciones reales x1 & x2 .
(∗) (0) (0)
x1 = x1 + ∆x1
(∗) (0) (0)
x2 = x2 + ∆x2
Donde los incrementos se formulan mas adelante.
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
Explicación teórica: Método Newton-Raphson
Por lo tanto, al sustituir las expresiones anteriores en las funciones
(1) y (2)
(∗) (∗) (0) (0) (0) (0)
g1 x1 , x2 , u = h1 x1 + ∆x1 , x2 + ∆x2 , u − b1 = 0 (3)
(∗) (∗) (0) (0) (0) (0)
g2 x1 , x2 , u = h2 x1 + ∆x1 , x2 + ∆x2 , u − b2 = 0 (4)
Se obtienen las expresiones que contienen los valores iniciales
(0) (1)
x1,2 y los incrementos iterativos ∆x1,2 .
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
Explicación teórica: Método Newton-Raphson
Al realizar la expansión de las Series de Taylor para (3) y (4)
anteriores se obtiene
0 0
(0) ∂g1 (0) ∂g1
(∗) (∗) (0) (0)
g1 x1 , x2 , u = g1 x1 , x2 , u + ∆x1 + ∆x2 +···=0 (5)
∂x1 ∂x2
0 0
(0) ∂g2 (0) ∂g2
(∗) (∗) (0) (0)
g2 x1 , x2 , u = g2 x1 , x2 , u + ∆x1 + ∆x2 +···=0 (6)
∂x1 ∂x2
Donde solo se considera la primera expansión de Taylor debido a
que es suficiente para obtener resultados óptimos.
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
Explicación teórica: Método Newton-Raphson
Realizando el acomodo de las expresiones anteriores en forma
matricial y sustituyendo (1) y (2) se obtiene lo siguiente
" ∂g ∂g1
#0 " #
(0) (0)
(0) (0)
1 (0) 0 − g1 x1 , x2 , u b1 − h1 x1 , x2 , u
∂x1 ∂x2 ∆x1
(0) =
∂g2 ∂g2
= (7)
(0) (0) (0) (0)
∂x1 ∂x2 ∆x2 0 − g2 x1 , x2 , u b2 − h2 x1 , x2 , u
donde la matriz que contiene las derivadas parciales se conoce
como Matriz Jacobiana, que se expresa mediante
" #0
∂g1 ∂g1
J0 ∈ Rn×n = ∂x1
∂g2
∂x2
∂g2
∂x1 ∂x2
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
Explicación teórica: Método Newton-Raphson
El (7) se busca determinar el valor de los incrementos por medio
de
" #
(0)
∆x1
(0) ,
∆x2
y finalmente, las ecuaciones de error se expresan por medio de
" #
(0) (0)
(0) (0)
(0)
∆g1 0 − g1 x1 , x2 , u b1 − h1 x1 , x2 , u
(0) =
=
(0) (0) (0) (0)
∆g2 0 − g2 x1 , x2 , u b2 − h2 x1 , x2 , u
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
Explicación teórica: Método Newton-Raphson
Por lo tanto, el método se puede expresar de forma general
" # " #
(0) (0)
∆x1 ∆g1
J(0) (0) = (0) (8)
∆x2 ∆g2
(0)
donde J(0) representa la matriz Jacobiana evaluada con x1,2 ,
(0) (0)
mientras que ∆x1,2 es el incremento de la iteración y ∆g1,2
(0)
expresa las funciones de error evaluada con x1,2 .
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
Explicación teórica: Método Newton-Raphson
Los valores para la siguiente iteración se expresan como
(1) (0) (0)
x1 = x1 + ∆x1
(1) (0) (0)
x2 = x2 + ∆x2 .
Ahora se deben usar los nuevos valores obtenidos en (7) para
determinar los incrementos que corresponde a la siguiente
iteración.
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
Contenido Explicación teórica
Explicación teórica: Método Newton-Raphson
Por lo que finalmente la siguiente forma de (8) es
" # " #
(1) (1)
(1) ∆x1 ∆g1
J (1) = (1)
∆x2 ∆g2
donde J(1) es la matriz Jacobiana evaluada con el nuevo valor de
(1) (1) (0)
x1,2 , mientras que ∆x1,2 es la nueva incógnita y ∆g1,2 expresa
(1)
las nuevas funciones de error evaluadas con x1,2 .
NOTA: Se debe repetir el proceso hasta que se cumpla con error
establecido o hasta que la función de error tenga valores muy
pequeños.
Canal de YouTube: Aprende con Romero Universidad de Guadalajara
