Programa Matlab Para Resolver Ejercicio 13.3 Del Libro Analisis De Sistemas De Potencia John J. Grainger & William D. Stevenson Jr.

ANALISIS DE SISTEMAS DE POTENCIA John J. Grainger & William D. Stevenson Jr.
Ejemplo 13.3. El sistema de cuatro barras de la figura 13.5 tiene los datos de líneas y barras dados en la tabla 13.2. Calcular los coeficientes B del sistema y demuestre que las pérdidas de transmisión calculadas con la formula de perdidas coinciden con los resultados del flujo de potencia.

Figura 13.5 Sistemade cuatro barras
Datos de Líneas | Datos de Barras |
| Z serie | Y paralelo | Generación | Carga |
De barra a barra | R | X | B | Barra | P | voltaje | ángulo | P | Q |
línea 1-4 | 0.00744 | 0.0372 | 0.0775 | 1 | | 1 | 0 | 0 | 0 |
línea 1-3 | 0.01008 | 0.0504 | 0.1025 | 2 | 3.18 | 1 | | 0 | 0 |
línea 2-3 | 0.00744 | 0.0372 | 0.0775 | 3 | 0 | 0 | 0 | 2.2 | 1.3634 |
línea 2-4 |0.01272 | 0.0636 | 0.1275 | 4 | 0 | 0 | 0 | 2.8 | 1.7352 |
Tabla 13.2
Solución de flujos de potencia |
| Generación | Voltaje |
Barra | P | Q | Magnitud | Angulo |
1 | 1.913152 | 1.87224 | 1 | 0 |
2 | 3.18 | 1.325439 | 1 | 2.43995 |
3 | 0 | 0 | 1.96051 | -1.07932 |
4 | 0 | 0 | 0.94304 | -2.62658 |

Algoritmo computacional que calcula las pérdidas del sistema |
clcclear all%%DatosNL=4;%Numero de Lineas de la redNb=4;%Numero de nodosNc=2;%Numero de cargasNg=2;%Numero de fuentes o generadoresslack=1;%Numero de nodo al que se encuentra conectado el generador pivote o slack en el sistema %Caracteristicas de la impedancia de la linea y su admitancia en derivacion%[Nodo salida: nodo llegada:Resistencia: Reactancia: Admitancia en derivacion]Lines=[1 4 0.00744 0.0372 0.0775/21 3 0.01008 0.0504 0.1025/2 2 3 0.00744 0.0372 0.0775/2 2 4 0.01272 0.0636 0.1275/2]; %Caracteristicas de la carga% [Nodo en que esta conectado: Potencia real demandada: Potencia reactiva demandada]Loads=[3 2.20 1.3634 4 2.80 1.7352]; %Caracteristicas de los generadores% [Nodo al que se encuentra conectado: Nodo al que se encuentra conectado: Potencia real generada:Potencia reactiva generada]Gen=[1 1 1.913152 1.872240 2 2 3.180000 1.325439 ]; %Voltajes nodales provenientes de una corrida de flujos% [Magnitud: Angulo en grados]V=[ 1.0 0 1.0 2.43995 0.96051 -1.07932 0.94304 -2.62658]; %% Calculo de la Ybus Ybus=zeros(Nb,Nb); %Crea una matriz de ceros de NbxNb for k=1:NL sen=Lines(k,1); %guarda el nodo en quese encuntra conectado la linea rece=Lines(k,2);%guarda el nodo en que se encuntra conectado la linea Z=Lines(k,3)+j*Lines(k,4); %Obtiene la impedancia de cada elemento de la matriz lines juntando la resistencia con la reactancia Y=(1/Z);%Saca el inverso de la matriz de impedancias Ybus(sen,rece)=Ybus(sen,rece)-Y;%Multiplica por un signo (-) todos los elementos de la matriz fuera de ladiagonal Ybus(rece,sen)=Ybus(rece,sen)-Y; Y=Y+j*Lines(k,5);%Agrega la admitancia en derivacion a la nueva matriz Y Ybus(sen,sen)= Ybus(sen,sen)+Y; Ybus(rece,rece)= Ybus(rece,rece)+Y;%Obtiene la Ybus del sistema end %% Calculo de Zbus Zbus=inv(Ybus); %Calcula la matriz de impedancias Rbus=real(Zbus);%Calculo de la matriz de la parte real de las impedancias de ZbusXbus=imag(Zbus);%Calculo de la matriz de la parte imaginaria de las impedancias de Zbus %% Calculo de corrientes de cargas R=V(:,1);%LLena un vector R con los voltajes nodales TH=V(:,2)*pi/180;%Convierte el angulo de grados a radianes [Vre,Vim] = pol2cart(TH,R)% Covierte de polar a rectangular Vaux=Vre+j*Vim;%Los valores rectangulares se guardan en este vector for k=1:Ncsen=Loads(k,1);%pregunta por el nodo en que se encuentran las cargas Ic(k,1)=(Loads(k,2)-j*Loads(k,3))/conj(Vaux(sen))%Dependiendo del nodo se calculan las corrientes de carga end %% calculo de la ID ID=ones(1,Nc)*Ic;%Crea un vector fila de 1xNc y lo multiplica por el vector columna Ic para calcular las corrientes totales de carga %% Calculo de las fracciones de carga D=Ic*(1/ID)%...