FACULTAD DE INGENIERÍA

DEPARTAMENTO DE INGENIERÍA ELÉCTRICA

VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

Matriz Ybus: teoría y ejercicios

Elaborado por: Christian Efraín Cardona Morales

La matriz Y bus es una matriz de números complejos que proporciona información de la

topología de un sistema de potencia. Se tiene entonces que la parte real de Y b us es G, lo cual
corresponde a la conductancia, es una matriz y se representa como en la ecuación 1, y la
parte imaginaria de Y bus es B, que corresponde a la susceptancia, es una matriz y se
representa como en la ecuación 2, por ende, la matriz Y bus es como se expresa en la
ecuación 3.
R { Y bus }=G (1)

I m {Y bus }=B (2)

Y bus =G+ jB (3)

Para llegar a la matriz Y bus se pueden utilizar dos métodos:

 Bloques constructivos.
 Por observación.
Método por bloque constructivo:
Lo primero que se debe analizar es como está el diagrama, si en impedancias o admitancias,
se recuerda que una bobina en impedancia es con signo positivo, mientras que en
admitancia es con signo negativo y viceversa para un condensador, esto con la razón de que
para este método lo que necesitamos es admitancias. El número de bloques dependerá de
que tantas ramas hallan, incluyendo las ramas que están conectadas a tierra, en donde la
tierra será el nodo cero, y el resto se enumeraran de 1 hasta tantos nodos hallan en el
diagrama diferentes de tierra.
Para ver cómo funciona el método se tiene el sistema de la figura 1, al cual se le sacara el
bloque constructivo, que se representara en la ecuación 4.

m + - n
Va
Ia
Im In
+ +
Ya = 1/Za
Vm

Vn

- -
Nodo de referencia

Figura 1. Diagrama del sistema.

Bloque constructivo de la matriz Y bus :

m n
m (4)
Ya −Ya Vm = ℑ
[−Ya Ya Vn ][ ]
¿ [ ] n
 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

A continuación, se desarrollará el ejemplo visto en clase por medio del método de bloque
constructivo:
Ejemplo 1:

-0,2j
Ŀ ° 0,2j
0,125j
2
1

0,1j 0,1j

0,2j

3 4

-0,1j 0,1j

0,2j Ŀ °

Figura 2. Diagrama en impedancias.

Solución:
Se puede observar que el diagrama de la figura 2 está en impedancias, por lo tanto, se debe
pasar los elementos de las ramas a admitancias y luego crear los bloques constructivos para
cada rama.
Se comienza con la rama 1-0:
En la rama 1-0 de la figura 2, se puede observar que hay una bobina de valor 2j y una fuente
de corriente que está inyectando 0,9∟0° en el nodo 1. Utilizando la ecuación 4 y pasando la
impedancia de la rama a admitancia con la ecuación 5 se tiene:
1 (5)
Y=
Z

1
Y= =−5 j
0,2 j

1 0
−5 j 5 j V 1 = 0,9∟ 0 ° 1
[ ][ ] [
5 j −5 j V 0 −0,9 ∟0 ° 0 ]
 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

La corriente de inyección al nodo cero puede ser 0 o -0,9∟0°, no importa, ya que la fila y
columna en cero se eliminan, debido a que no se tiene en cuenta en la matriz Y bus, por lo
tanto, se tiene al final la siguiente matriz:

[−55 j j 5 j V 1 = 0,9∟ 0 °
−5 j V 0][ ] [
−0,9 ∟0 ° ]
1
[ −5 j ] [ V 1 ] =[ 0,9 ∟0 ° ] 1
Siempre que se tenga conectada la rama al nodo cero se elimina la columna y fila de este
nodo.
Se prueba con otra rama conectada a tierra, la rama 2-0:
Se tiene entonces que para esta rama hay un condensador de valor -2j y no hay inyección de
corriente, utilizando las ecuaciones 4 y 5, y sabiendo que la fila y columna del nodo cero se
eliminan, quedara la siguiente expresión:
1
Y= =5 j
−0,2 j

2
[ 5 j ] [ V 2 ] =[ 0 ] 2

Ya se sabe cómo trabajar con las ramas conectadas a tierra, por lo que para las ramas 4-0,
3-0 y 5-0 se hace de la misma forma.
Ahora se continua con las ramas que están conectadas entre nodos diferentes de tierra:
Rama 3-5:
En esta rama se puede observar que hay una bobina con valor de 0,2j, y hay una inyección
de corriente en el nodo 5 de 1∟10°, utilizando las ecuaciones 4 y 5 se tiene:
1
Y= =−5 j
0,2 j

3 5
3
[−55 j j 5 j V3 =
−5 j V 5
0
][ ] [
1 ∟10 ° ] 5

Si se quisiera tomar la rama de manera que fuera 5-3, se podría hacer, pero teniendo en
cuenta que el bloque 3-5 no se haya hecho, ya que, cada rama representa un bloque, por lo
que la rama 5-3 y 3-5 son lo mismo, con la diferencia que en la matriz se cambie la posición
de los voltajes y corrientes.
Rama 4-5:
Como se dijo anteriormente, la rama se pudo haber cogido 5-4 pero teniendo en cuenta la
posición de los voltajes y corrientes para la ecuación 4. En esta rama se puede observar que
hay una bobina de valor 0,25j, y hay inyección de corriente en el nodo 5 de 1∟10°, teniendo
la rama analizada y utilizando las ecuaciones 4 y 5 se tiene:
1
Y= =−4 j
0,25 j

4 5
4
5
 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

[−44 j j 4j V 4 =
][ ] [
−4 j V 5
0
1 ∟10° ]
Rama 1-3:
Se puede observar que en la rama 1-3 de la figura 2, hay una bobina de valor 0,1j, y hay una
inyección de corriente de 0,9∟0° en el nodo 1, utilizando las ecuaciones 4 y 5 se tiene:
1
Y= =−10 j
0,1 j

1 3
j 10 j V 1 = 0,9∟ 0°1
[−10
10 j ][ ] [
−10 j V 3 0 3 ]
Rama 1-2:
Se puede observar que en la rama 1-2 de la figura 2, hay una bobina de valor 0,125j, y hay
inyección de corriente de 0,9∟0° en el nodo 1, utilizando las ecuaciones 4 y 5 se tiene:
1
Y= =−8 j
0,125 j

1 2
1
[−88 j j 8 j V 1 = 0,9∟ 0 °
][ ] [
−8 j V 2 0 ] 2

Luego de realizar el proceso para cada rama se procede a calcular el valor de cada posición
de la matriz Y bus. Para obtener los valores de cada posición de la matriz Y bus , se toma los
valores en cada bloque constructivo que contenga esta posición y se suman, como se
muestra en la ecuación 6, a continuación, se mostrará como hallar la posición (1,1) de la
matriz Y bus:

 Posición (1,1) de la rama 1-0 = -5j

 Posición (1,1) rama 1-3 = -10j
 Posición (1,1) rama 1-2 = -8j

Y bus(1,1) =Rama [1−0](1,1) + Rama [1−3](1,1 )+ Rama [1−2](1,1 ) (6)

Y bus(1,1) =−5 j−10 j−8 j

Y bus(1,1) =−23 j

Luego de realizar el proceso para cada posición, sumando los valores de los bloques
constructivos que contengan la misma posición, se llega a la matriz Y bus de 5x5:

Figura 3. Matriz Ybus 5x5.

 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

Método por observación:

Este método consiste en observar el diagrama y sacar la matriz Y bus, teniendo en cuenta que
para la diagonal es la suma de las ramas que se conectan en cada nodo y para las otras
posiciones es la admitancia con signo contrario, para el ejemplo 1 se tiene:

Para la posición (3,3) de la diagonal de la matriz Y bus:

De la figura 2 se tiene que al nodo 3 le llegan la rama 1-3, 4-3, 5-3, y 3-0, con valores de 0,1j,
0,2j, 0,2j y -0,1j respectivamente, que luego de pasarlos a admitancias queda -10j, -5j y -5j y
10j y al sumarlos da:
Y bus(3,3) =−10 j−5 j−5 j+10 j

Y bus(3,3) =−10 j

Ahora para la posición (3,4) y (4,3) de la matriz Y bus :

Debido a que la posición (3,4) y (4,3) comparten el mismo elemento estas dos tendrán el
mismo valor, que se halla obteniendo la admitancia que comparte esta rama la cual es 3-4 o
4-3, esta admitancia es de valor -5j y con signo contrario sería 5j.
Cuando no hay conexión entre nodos el valor es cero, como en el caso de la posición (1,4),
en donde se puede ver en la figura 2 que no hay conexión alguna.
Observando la figura 3, se puede concluir que el método por observación arrojo los mismos
resultados que el método por bloques constructivos.
Ejemplo 2:

Teniendo una matriz Y bus 4 x4 , sacar el diagrama de admitancias.

Figura 4. Matriz Ybus 4x4.

Solución:
 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

-4j 5j
-j
-2j

1 2

-4j -8j

4
3

-4j 13j 0,9

Figura 5. Diagrama de admitancias.

Ramas acopladas mutuamente en Ybarra

Cuando se tiene las ramas acopladas ya no es posible utilizar el método por observación, por
lo tanto, se debe hacer otro tipo de bloque constructivo, el cual es el que se muestra en la
figura 6.

Figura 6. Matriz de nodos.

La figura 6 es una demostración que sale de la figura 7.

m p
Ym
Ip
+
+
Vb

Ya Yb
Va

In - -
Iq

n q

Figura 7. Ramas mutuamente acopladas: parámetros de admitancia.

Donde Y a , Y M y Y b se calculan con la inversa de la matriz de impedancias, que corresponde a

la ecuación 7.

Za Z M
−1
1 Z b −Z M Y YM (7)
[ Z M Zb ] = 2
[
Z a Z b−Z M −Z M Za
= a
YM Yb ][ ]
 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

Ejemplo 3:
0,1j

2
1

0,2j

3
Za

Zm
4 5

-0,2j
-0,5j

Figura 8. Diagrama de impedancias.

Se tiene que:

 Za = 0,5j
 Zb = 0,3j
 Zm = 0,2j
De la figura 7 y 8 se puede observar que el nodo 4 es m y p, el nodo 1 es n y el nodo 3 es q.
Se continua entonces con el cálculo del bloque constructivo de la rama acoplada, con
respecto a la figura 6:
4 1 4 3
4
[−11 −11 ]Y [−11 −11 ] Y
[ a

[−11 −11 ] Y [−11 −11 ] Y

M
M

b ] 1
4
3

Se procede a calcular Y a , Y M y Y b con la ecuación 7:

−1

[ 0,5 j 0,2 j
0,2 j 0,3 j ] = 0,5 j∗0,31j−(0,2 j) [−0,2
0,3 j
j 2
−0,2 j = −2,7273 j 1,8182 j
0,5 j ][
1,8182 j −4,5455 j ]
Así, el bloque de la rama acoplada queda de la siguiente manera:

4 1 4 3
4
1 −1 −2,7273 j 1 1 −1 1,8182 j
[ ] [ ]
[[ −1 1
1 −1
−1 1 ]
1,8182 j
4
31
[
−1 1
−1
−1 1 ]
−4,5455 j ]
Se procede entonces a calcular los bloques constructivos faltantes:
Rama 1-2:
De la figura 8 se puede observar que para esta rama hay un inductor de valor 0,1j, y no hay
inyección de corriente en ninguno de los dos nodos, así, utilizando las ecuaciones 4 y 5 se
tiene:
1
Y= =−10 j
0,1 j
 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

1 2
j 10 j V 1 = 0 1
[−10
10 j −10 j V 2 ][ ] [ ]
0 2

Rama 2-3:
De la figura 8 se puede observar que para esta rama hay un inductor de valor 0,2j, y no hay
inyección de corriente en ninguno de los dos nodos, así, utilizando las ecuaciones 4 y 5 se
tiene:
1
Y= =−5 j
0,2 j

2 3
−5 j 5 j V 2 = 0 2
[ 5 j −5 j V 3 0][ ] [ ] 3

Rama 3-5:
De la figura 8 se puede observar que para esta rama hay un inductor de valor 0,1j, y no hay
inyección de corriente en ninguno de los dos nodos, así, utilizando las ecuaciones 4 y 5 se
tiene:
1
Y= =−10 j
0,1 j

3 5
j 10 j V 3 = 0 3
[−10
10 j −10 j V 5 ][ ] [ ]
0 5

Rama 5-0:
De la figura 8 se puede observar que para esta rama hay un condensador de valor -0,5j, y no
hay inyección de corriente, además, hay una conexión en el nodo 0, el cual corresponde a la
tierra, sabiendo lo anterior y utilizando las ecuaciones 4 y 5 se tiene:
1
Y= =2 j
−0,5 j

5
[ 2 j ][ V 5 ] =[ 0 ] 5

Rama 4-0:
De la figura 8 se puede observar que para esta rama hay un condensador de valor -0,2j, y no
hay inyección de corriente, además, hay una conexión en el nodo 0, el cual corresponde a la
tierra, sabiendo lo anterior y utilizando las ecuaciones 4 y 5 se tiene:
1
Y= =5 j
−0,2 j

4
[ 5 j ] [ V 4 ] =[ 0 ] 4

Teniendo los bloques constructivos, se continua con calcular cada posición de la matriz Y bus ,
que para este caso será una matriz 5x5 y se hará como se explicó en la ecuación 6:
 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

Se comenzará con la diagonal:

Y bus(1,1) =−2,7273 j−10 j=−12,7273 j

Y bus(2,2) =−10 j−5 j=−15 j

Y bus(3,3) =−4,5455 j−5 j−10 j=−19,5455 j

Y bus(4,4 )=−2,7273 j+1,8182 j +1,8182 j−4,5455 j+5 j=1,3636 j

Y bus(5,5) =−10 j+2 j=−8 j

Ahora se continuará con las posiciones faltantes:

Y bus (1,2 )=10 j

Y bus ( 1,3 )=1,8182 j

Y bus (1,4 )=2,7273 j−1,8182 j=0,9091 j

Y bus (1,5 )=0

Y bus (2,3 )=5 j

Y bus (2,4 )=0

Y bus (2,5 )=0

Y bus (3,4 )=−1,8182+4,5455 j=2,7273 j

Y bus ( 3,5 )=10 j

Y bus (4,5 )=0

Teniendo los valores de las posiciones se procede a hacer la matriz Y bus 5 x 5

Figura 9. Matrix Ybus 5x5, por método de bloques constructivos.

Red de admitancias equivalente

Esta red equivalente me ayuda a obtener un circuito con el cual puedo eliminar el acople que
hay en las ramas, y utilizar el método por observación para hallar la matriz Y bus, en la figura
10 se muestra el circuito del cual se parte para obtener la red equivalente de la figura 11.
 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

m p

Im Zm
Ip
+
+

Vb
Za Zb
Va

In - -
Iq

n q

Figura 10. Ramas mutuamente acopladas: parámetros de impedancias

-Ym p
Im m Ip

Ym
Ya as Ym Yb

In
Iq
n -Ym q

Figura 11. Red equivalente en admitancias.

Tomando el ejemplo 3, se procede a sacar el circuito equivalente con respecto a la figura 11.

m=p

n q
-Ym

Figura 12. Circuito equivalente del ejemplo 3.

Realizando el paralelo entre Ya y Ym, y Yb y Ym se tiene:

m=p

-0,9091j -2,7273j

n q

-1,8182j

Figura 13. Circuito equivalente reducido.

Ahora se procede a colocarlo en el diagrama de la figura 6, con todos los elementos en

admitancias, para así poder aplicar el método por observación.
 FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA ELÉCTRICA
VISIÓN: "El programa de Ingeniería Eléctrica de la Universidad de Antioquia será reconocido por sus
aportes al MANEJO EFICIENTE DE LA ENERGÍA ELÉCTRICA en la región y el país, manifestados
en la formación de ingenieros emprendedores y en proyectos de INVESTIGACIÓN, INNOVACIÓN Y
DESARROLLO."

-10j

2
1
-1,8182j
-5j

-0,9091j 3

-2,7273j

4 5

5j
2j

Figura 14. Diagrama en admitancias.

Aplicando el método por observación se obtiene:

Figura 15. Matriz Ybus 5x5, método por observación.

Observando la figura 9 que es la matriz Y bus por método de bloques constructivos y la figura
15 que es la matriz Y bus por el método de por observación, se concluye que la red equivalente
cumple perfectamente y facilita el trabajo, ya que, se crea un modelo con el cual se puede
sacar la matriz Y bus más rápidamente.