clc, clear all Va1=Vfalla-Zposit*Ia1;

%FALLA LINEA - LINEA fprintf('Va1= %9.4f angulo %9.4f ° pu

disp('FALLA LINEA- LINEA') \n',abs(Va1),radtodeg(angle(Va1)))
Vbase=220; Va2=0-Znegat*Ia2;
Sb=25; fprintf('Va2= %9.4f angulo %9.4f ° pu
Zb=Vbase^2/Sb; \n',abs(Va2),radtodeg(angle(Va2)))
Zf=4; fprintf('Va0= %9.4f angulo %9.4f ° pu
Zfpu=Zf/Zb; \n',abs(Va0),radtodeg(angle(Va0)))
Xg1=j*0.125*(25/25)*(11/11)^2; Va=Va1+Va2+Va0;
Xt1=j*0.105*(25/30)*(220/220)^2; fprintf('Va= %9.4f angulo %9.4f ° pu \n',abs(Va),radtodeg(angle(Va)))
XL=(2+6*j)/Zb; Vb=a2*Va1+a1*Va2;
Xt2=j*0.10*(25/30)*(220/220)^2; fprintf('Vb= %9.4f angulo %9.4f ° pu \n',abs(Vb),radtodeg(angle(Vb)))
zpc=(10.8/11)^2/conj((10+8*j)/Sb); Vc=a1*Va1+a2*Va2;
disp('-----------------DATOS EN PU----------------------') fprintf('Vc = %9.4f angulo %9.4f ° pu \n',abs(Vc),radtodeg(angle(Vc)))
fprintf('Xgenerador 1= %9.4fj pu \n',imag(Xg1)) fprintf('\n')
fprintf('X transf. 1= %9.4fj pu \n',imag(Xt1)) Vb0=Va0; Vc0=Va0;
fprintf('X linea= %9.4f + %9.4fj pu \n',real(XL),imag(XL)) Vab=Va-Vb;
fprintf('X transf. 2= %9.4f j pu \n',imag(Xt2)) fprintf('Vab= %9.4f angulo %9.4f ° pu
fprintf('Zcarga= %9.4f + %9.4fj pu \n',real(zpc),imag(zpc)) \n',abs(Vab),radtodeg(angle(Vab)))
fprintf('\n') Vbc=Vb-Vc;
disp('-----------------RED DE SECUENCIA POSITIVA--------') fprintf('Vbc= %9.4f angulo %9.4f ° pu
Z1=(Xg1+Xt1+XL); Z2=(Xt2+zpc); \n',abs(Vbc),radtodeg(angle(Vbc)))
Zposit=Z1*Z2/(Z1+Z2); Vca=Vc-Va;
fprintf('Zpositivo= %9.6f + %9.4fj pu \n',real(Zposit),imag(Zposit)) fprintf('Vca= %9.4f angulo %9.4f ° pu
fprintf('\n') \n',abs(Vca),radtodeg(angle(Vca)))
disp('-----------------RED DE SECUENCIA NEGATIVA--------') disp('COMPROBAMOS Vb-Vc=Ib*Zf ')
Znegat=Zposit; Vbcfallaz=Zfpu*Ib;
fprintf('Znegativo= %9.6f + %9.4fj pu \n',real(Znegat),imag(Znegat)) fprintf('Vbcfalla con Zf= %9.4f angulo %9.4f ° pu
fprintf('\n') \n',abs(Vbcfallaz),radtodeg(angle(Vbcfallaz)))
disp('-----------------RED DE SECUENCIA CERO------------') FALLA LINEA- LINEA
Z0linea=3*XL; -----------------DATOS EN PU----------------------
Z0carga=3*zpc; Z3=Xt1+Z0linea; Z4=Xt2; Xgenerador 1= 0.1250j pu
Zcero=Z3*Z4/(Z3+Z4); X transf. 1= 0.0875j pu
fprintf('Zcero= %9.6f + %9.4fj pu \n',real(Zcero),imag(Zcero)) X linea= 0.0010 + 0.0031j pu
disp('Nuestra red de secuencia cero para una falla linea-linea') X transf. 2= 0.0833 j pu
disp('no interviene en el analisis') Zcarga= 1.4695 + 1.1756j pu
fprintf('\n') -----------------RED DE SECUENCIA POSITIVA-----------------
disp('-----------------I carga y V falla---------------') Zpositivo= 0.016643 + 0.1996j pu
Icarga=1/zpc; -----------------RED DE SECUENCIA NEGATIVA-----------------
fprintf('Icarga= %9.4f %9.4f j pu \n',real(Icarga),imag(Icarga)) Znegativo= 0.016643 + 0.1996j pu
%Voltage de falla en la barra 3 -----------------RED DE SECUENCIA CERO---------------------
Vfalla=1+Icarga*Xt2; Zcero= 0.000663 + 0.0448j pu
fprintf('Vfalla= %9.4f angulo %9.4f ° pu Nuestra red de secuencia cero para una falla linea-linea
\n',abs(Vfalla),radtodeg(angle(Vfalla))) no interviene en el analisis
Zfallapu=3*Zfpu; -----------------I carga y V falla----------------------
%% CORRIENTES DE SECUENCIA POSITIVA Icarga= 0.4150 -0.3320 j pu
a1=-1/2+i*sqrt(3)/2; Vfalla= 1.0282 angulo 1.9272 ° pu
a2=-1/2-i*sqrt(3)/2; -------------------CORRIENTES Ia1
A=[1 1 1;1 a2 a1;1 a1 a2]; Ia1= 2.5652 angulo -83.0130 ° pu
disp('CORRIENTES Ia1') CORRIENTES DE FALLA
Vprefalla=Vfalla; Ia falla= 0.0000 pu
Ia1=Vprefalla/(Zposit+Znegat+Zfpu); Ib falla= 4.4430 angulo -173.0130 ° pu
fprintf('Ia1= %9.4f angulo %9.4f ° pu Ib amperes= 291.4966 A
\n',abs(Ia1),radtodeg(angle(Ia1))) Ic falla= 4.4430 angulo 6.9870 ° pu
Iabc=A*[0;Ia1;-Ia1]; Ic amperes= 291.4966 A
disp('CORRIENTES DE FALLA') Ia2= 2.5652 angulo 96.9870 ° pu
Ia=0; LOS VOLTAGES DE SECUENCIA SERÁN
fprintf('Ia falla= %9.4f pu \n',abs(Ia)) Va1= 0.5144 angulo 1.6332 ° pu
Ib=(a2-a1)*Ia1; Va2= 0.5139 angulo 2.2215 ° pu
fprintf('Ib falla= %9.4f angulo %9.4f ° pu Va0= 0.0000 angulo 0.0000 ° pu
\n',abs(Ib),radtodeg(angle(Ib))) LOS TENSIONES EN LAS LINEAS SERÁN
Ibase=Sb*10^3/(sqrt(3)*Vbase); Va= 1.0282 angulo 1.9272 ° pu
IbA=Ibase*abs(Ib); Vb= 0.5187 angulo -178.0281 ° pu
fprintf(' Ib amperes= %9.4f A \n',abs(IbA)) Vc = 0.5096 angulo -178.1183 ° pu
Ic=(a1-a2)*Ia1; LOS TENSIONES DE LINEA-LINEA SOBRE UNA BASE DE VOLTAGE DE FASE SON:
fprintf('Ic falla= %9.4f angulo %9.4f ° pu Vab= 1.5469 angulo 1.9422 ° pu
\n',abs(Ic),radtodeg(angle(Ic))) Vbc= 0.0092 angulo -173.0130 ° pu
IcA=Ibase*abs(Ic); Vca= 1.5378 angulo -178.0879 ° pu
fprintf(' Ic amperes= %9.4f A \n',abs(IcA)) COMPROBAMOS Vb-Vc=Ib*Zf
Ia2=-Ia1; Vbcfalla con Zf= 0.0092 angulo -173.0130 ° pu
fprintf('Ia2= %9.4f angulo %9.4f ° pu
\n',abs(Ia2),radtodeg(angle(Ia2)))
clc, clear all fprintf('Va2= %9.4f angulo %9.4f ° pu
%FALLA LINEA - LINEA \n',abs(Va2),radtodeg(angle(Va2)))
% Ocurre una falla 2LG con Zf=2.5 Ohms y Zg=7 Ohms. La falla ocurre en la fprintf('\n')
barra C. %%
disp('FALLA DOBLE LINEA- TIERRA') disp('LOS VOLTAGES DE FALLA SERÁN')
Vbase=220;
Sb=25; Va=Va0+Va1+Va2;
Zb=Vbase^2/Sb; fprintf('Va= %9.4f angulo %9.4f ° pu \n',abs(Va),radtodeg(angle(Va)))
Zf=2.5;Zfpu=Zf/Zb; Vb=Va0+a2*Va1+a1*Va2;
Zg=7;Zgpu=Zg/Zb; fprintf('Vb= %9.4f angulo %9.4f ° pu \n',abs(Vb),radtodeg(angle(Vb)))
Vc=Va0+a1*Va1+a2*Va2;
Xg1=j*0.125*(25/25)*(11/11)^2; fprintf('Vc = %9.4f angulo %9.4f ° pu \n',abs(Vc),radtodeg(angle(Vc)))
Xt1=j*0.105*(25/30)*(220/220)^2; fprintf('\n')
XL=(2+6*j)/Zb;
Xt2=j*0.10*(25/30)*(220/220)^2; disp('CAIDAS DE TENSION EN LA RAMA EN PARALELO')
zpc=(10.8/11)^2/conj((10+8*j)/Sb); deltaV=Va1-Zfpu*Ia1;
fprintf('deltaV= %9.4f angulo %9.4f ° pu
disp('-----------------DATOS EN PU----------------------') \n',abs(deltaV),radtodeg(angle(deltaV)))
fprintf('Xgenerador 1= %9.4fj pu \n',imag(Xg1)) fprintf('\n')
fprintf('X transf. 1= %9.4fj pu \n',imag(Xt1)) %%
fprintf('X linea= %9.4f + %9.4fj pu \n',real(XL),imag(XL)) disp('COMPROBANDO LAS TENSIONES DE FASE B Y C')
fprintf('X transf. 2= %9.4f j pu \n',imag(Xt2)) %formula 4.13 y 4.15 respectivamente
fprintf('Zcarga= %9.4f + %9.4fj pu \n',real(zpc),imag(zpc)) Vbk1=(Zfpu+Zgpu)*Ib+Zgpu*Ic;
fprintf('\n') fprintf('Vb= %9.4f angulo %9.4f ° pu
disp('-----------------RED DE SECUENCIA POSITIVA----------------------') \n',abs(Vbk1),radtodeg(angle(Vbk1)))
Vck1=(Zfpu+Zgpu)*Ic+Zgpu*Ib;
Z1=(Xg1+Xt1+XL); Z2=(Xt2+zpc); fprintf('Vc= %9.4f angulo %9.4f ° pu
Zposit=Z1*Z2/(Z1+Z2); \n',abs(Vck1),radtodeg(angle(Vck1)))
fprintf('Zpositivo= %9.6f + %9.4fj pu \n',real(Zposit),imag(Zposit)) FALLA DOBLE LINEA- TIERRA
fprintf('\n') -----------------DATOS EN PU----------------------
disp('-----------------RED DE SECUENCIA NEGATIVA----------------------')
Znegat=Zposit; Xgenerador 1= 0.1250j pu
fprintf('Znegativo= %9.6f + %9.4fj pu \n',real(Znegat),imag(Znegat)) X transf. 1= 0.0875j pu
fprintf('\n') X linea= 0.0010 + 0.0031j pu
disp('-----------------RED DE SECUENCIA CERO----------------------')
X transf. 2= 0.0833 j pu
Z0linea=3*XL; Zcarga= 1.4695 + 1.1756j pu
Z0carga=3*zpc; Z3=Xt1+Z0linea; Z4=Xt2; -----------------RED DE SECUENCIA POSITIVA----------
Zcero=Z3*Z4/(Z3+Z4); Zpositivo= 0.016643 + 0.1996j pu
fprintf('Zcero= %9.6f + %9.4fj pu \n',real(Zcero),imag(Zcero))
fprintf('\n')
-----------------RED DE SECUENCIA NEGATIVA----------
disp('-----------------I carga y V falla----------------------') Znegativo= 0.016643 + 0.1996j pu
Icarga=1/zpc;
fprintf('Icarga= %9.4f %9.4f j pu \n',real(Icarga),imag(Icarga))
%Voltage de falla en la barra 3 -----------------RED DE SECUENCIA CERO--------------
Vfalla=1+Icarga*Xt2; Zcero= 0.000663 + 0.0448j pu
fprintf('Vfalla= %9.4f angulo %9.4f ° pu
\n',abs(Vfalla),radtodeg(angle(Vfalla)))
fprintf('\n') -----------------I carga y V falla-----------------
disp('------------------------------------------------------------') Icarga= 0.4150 -0.3320 j pu
%% CORRIENTES DE SECUENCIA POSITIVA Vfalla= 1.0282 angulo 1.9272 ° pu
a1=-1/2+i*sqrt(3)/2;
a2=-1/2-i*sqrt(3)/2; ----------------------------------------------------
A=[1 1 1;1 a2 a1;1 a1 a2]; CORRIENTES DE SECUENCIA
disp('CORRIENTES DE SECUENCIA') Ia1= 4.3208 angulo -81.5459 ° pu
Vprefalla=Vfalla;
Zt=Zposit+Zfpu+((Znegat+Zfpu)*(Zcero+Zfpu+3*Zgpu)/(Zcero+Znegat+2*Zfpu+3* Ia0= 3.5155 angulo 100.4877 ° pu
Zgpu)); Ia2= 0.8170 angulo 89.6713 ° pu
Ia1=Vprefalla/Zt;
fprintf('Ia1= %9.4f angulo %9.4f ° pu
\n',abs(Ia1),radtodeg(angle(Ia1))) CORRIENTES DE FALLA
Ia0=-Ia1*((Zposit+Zfpu)/(Znegat+2*Zfpu+Zcero+3*Zgpu)); Ia falla= 0.0000 angulo 180.0000 pu
fprintf('Ia0= %9.4f angulo %9.4f ° pu Ib falla= 7.0954 angulo 139.1694 ° pu
\n',abs(Ia0),radtodeg(angle(Ia0)))
Ia2=-Ia1*((Zcero+Zfpu+3*Zgpu)/(Zcero+2*Zfpu+3*Zgpu+Znegat)); Ic falla= 6.6890 angulo 58.9614 ° pu
fprintf('Ia2= %9.4f angulo %9.4f ° pu
\n',abs(Ia2),radtodeg(angle(Ia2))) COMPROBANDO LAS CORRIENTES
fprintf('\n')
%% Ifalla1 (Ib+Ic)= 10.5466 angulo 100.4877 ° pu
disp('CORRIENTES DE FALLA') Ifalla2 (3*Ia0)= 10.5466 angulo 100.4877 ° pu
Ia=Ia0+Ia1+Ia2;
fprintf('Ia falla= %9.4f angulo %9.4f pu
\n',abs(Ia),radtodeg(angle(Ia))) LOS VOLTAGES DE SECUENCIA SERÁN
Ib=Ia0+a2*Ia1+a1*Ia2; Va0= 0.1575 angulo 9.6396 ° pu
fprintf('Ib falla= %9.4f angulo %9.4f ° pu Va1= 0.1652 angulo -7.3406 ° pu
\n',abs(Ib),radtodeg(angle(Ib)))
Ibase=Sb*10^3/(sqrt(3)*Vbase); Va2= 0.1637 angulo -5.0942 ° pu
IbA=Ibase*abs(Ib);
%fprintf(' Ib amperes= %9.4f A \n',abs(IbA)) LOS VOLTAGES DE FALLA SERÁN
Ic=Ia0+a1*Ia1+a2*Ia2;
fprintf('Ic falla= %9.4f angulo %9.4f ° pu Va= 0.4822 angulo -1.1015 ° pu
\n',abs(Ic),radtodeg(angle(Ic))) Vb= 0.0456 angulo 107.6946 ° pu
IcA=Ibase*abs(Ic); Vc = 0.0450 angulo 93.1712 ° pu
%fprintf(' Ic amperes= %9.4f A \n',abs(IcA))
fprintf('\n')
%% CAIDAS DE TENSION EN LA RAMA EN PARALELO
disp('COMPROBANDO LAS CORRIENTES') deltaV= 0.1638 angulo -5.4620 ° pu
Ifalla1=Ib+Ic; Ifalla2=3*Ia0;
fprintf('Ifalla1 (Ib+Ic)= %9.4f angulo %9.4f ° pu
\n',abs(Ifalla1),radtodeg(angle(Ifalla1))) COMPROBANDO LAS TENSIONES DE FASE B Y C
fprintf('Ifalla2 (3*Ia0)= %9.4f angulo %9.4f ° pu Vb= 0.0456 angulo 107.6946 ° pu
\n',abs(Ifalla2),radtodeg(angle(Ifalla2)))
fprintf('\n') Vc= 0.0450 angulo 93.1712 ° pu
%%
disp('LOS VOLTAGES DE SECUENCIA SERÁN')
Va0=-Zcero*Ia0;
fprintf('Va0= %9.4f angulo %9.4f ° pu
\n',abs(Va0),radtodeg(angle(Va0)))
Va1=Vprefalla-Zposit*Ia1;
fprintf('Va1= %9.4f angulo %9.4f ° pu
\n',abs(Va1),radtodeg(angle(Va1)))
Va2=-Znegat*Ia2;
clc;clear all;format short; Znegat=(g1+g4)*(g2+g5)/(g1+g2+g4+g5);
%% DATOS EN PU
SB=100; fprintf('Znegativo= %9.6f + %9.4fj pu
VB=138; \n',real(Znegat),imag(Znegat))
ZB=VB^2/SB; fprintf('\n')
a=cosd(120)+j*sind(120);
disp('-----------------RED DE SECUENCIA CERO---------------')
Xg1_1=j*.275*(SB/27.5)*(10/10)^2; eq1=Xs_T1/2+XL1_0/2;
Xg1_2=j*.165*SB/27.5; Zcero=eq1*(Xs_T2/2)/(eq1+Xs_T2/2);
Xg1_0=j*.11*SB/27.5; fprintf('Zcero= %9.6f + %9.4fj pu
\n',real(Zcero),imag(Zcero))
XT1=j*0.0966*SB/27.6; fprintf('\n')
Xp_T1=XT1/5;
Xs_T1=4*XT1/5; %%
a1=-1/2+i*sqrt(3)/2;
XT2=j*0.09625*SB/27.5; a2=-1/2-i*sqrt(3)/2;
Xp_T2=XT2/5; A=[1 1 1;1 a2 a1;1 a1 a2];
Xs_T2=4*XT2/5;
disp('CORRIENTES DE SECUENCIA EN NODO 4')
XL1=((j*0.484)*69.96)/ZB; Vprefalla=V4;
XL1_0=3*XL1; Zt=Zposit+Znegat+Zcero;
Ia0=Vprefalla/Zt;
Z1=(13.8/13.8)^2/(conj(20+j*10)/SB); fprintf('Ia 0= %9.4f angulo %9.4f ° pu
\n',abs(Ia0),radtodeg(angle(Ia0)))
disp('-----------------DATOS EN PU----------------------') Ia1=Ia0; Ia2=Ia0;
fprintf('Xgenerador 1 (positivo)= %9.4fj pu \n',imag(Xg1_1)) fprintf('Ia 1= %9.4f angulo %9.4f ° pu
fprintf('Xgenerador 1(negativo)= %9.4fj pu \n',imag(Xg1_2)) \n',abs(Ia1),radtodeg(angle(Ia1)))
fprintf('Xgenerador 1(cero)= %9.4fj pu \n',imag(Xg1_0)) fprintf('Ia 2= %9.4f angulo %9.4f ° pu
fprintf('\n') \n',abs(Ia2),radtodeg(angle(Ia2)))
fprintf('X transf. 1= %9.4fj pu \n',imag(XT1))
fprintf('X transf. 1(primario)= %9.4fj pu \n',imag(Xp_T1)) I1=[Ia0;Ia1;Ia2];
fprintf('X transf. 1(secundario)= %9.4fj pu \n',imag(Xs_T1)) II=A*I1;
fprintf('\n') disp('CORRIENTES DE FALLA FASICA EN NODO 4')
fprintf('X transf. 2= %9.4fj pu \n',imag(XT2)) fprintf('Ia= %9.4f angulo %9.4f ° pu
fprintf('X transf. 2(primario)= %9.4fj pu \n',imag(Xp_T2)) \n',abs(II(1)),radtodeg(angle(II(1))))
fprintf('X transf. 2(secundario)= %9.4fj pu \n',imag(Xs_T2)) fprintf('Ib= %9.4f angulo %9.4f ° pu
fprintf('\n') \n',abs(II(2)),radtodeg(angle(II(2))))
fprintf('X linea(positivo)= %9.4f + %9.4fj pu fprintf('Ic= %9.4f angulo %9.4f ° pu
\n',real(XL1),imag(XL1)) \n',abs(II(3)),radtodeg(angle(II(3))))
fprintf('X linea (cero)= %9.4f + %9.4fj pu
\n',real(XL1_0),imag(XL1_0)) disp('CORRIENTES de SECUENCIA DE FALLA EN G1')
fprintf('\n') Xizq=z4+z1;
fprintf('Z carga= %9.4f + %9.4fj pu \n',real(Z1),imag(Z1)) Xder=z5+z2;
%% B)DETERMINAR LA TENSION DE LA BARRA DEL GENERADOR 2 I1g1=(Ia1*(Xder)/(Xder+Xizq))/2;
fprintf('I1g1= %9.4f angulo %9.4f ° pu
V5=1; \n',abs(I1g1),radtodeg(angle(I1g1)))
S5=20+j*10; Xizq_2=g4+g1;
I5=(conj(S5)/100)/conj(V5); Xder_2=g5+g2;
I5p=[abs(I5) angle(I5)*180/pi];
IG1_2=(Ia1*Xder_2/(Xder_2+Xizq_2))/2;
V4=V5+XT2*I5/2; IG1_2p=[abs(IG1_2) angle(IG1_2)*180/pi];
V4p=[abs(V4) angle(V4)*180/pi];
IG1_0=0;
V3=V4+XL1*I5/2; fprintf('I2g1= %9.4f angulo %9.4f ° pu
V3p=[abs(V3) angle(V3)*180/pi]; \n',abs(IG1_2),radtodeg(angle(IG1_2)))
fprintf('I0g1= %9.4f angulo %9.4f ° pu
V1=V3+XT1*I5/2; \n',abs(IG1_0),radtodeg(angle(IG1_0)))
V1p=[abs(V1) angle(V1)*180/pi]; fprintf('\n')
V2=V1; disp('CORRIENTES FASICAS DE FALLA EN G1')
fprintf('\n') IG1_012=[IG1_0;I1g1;IG1_2];
disp('------------TENSIONES EN CADA BARRA-----------------')
fprintf('V5= %9.4f angulo %9.4f ° pu IG1_A=IG1_0+I1g1+IG1_2;
\n',abs(V5),radtodeg(angle(V5))) IG1_Ap=[abs(IG1_A) angle(IG1_A)*180/pi];
fprintf('V4= %9.4f angulo %9.4f ° pu
\n',abs(V4),radtodeg(angle(V4))) IG1_B=IG1_0+a^2*I1g1+a*IG1_2;
fprintf('V3= %9.4f angulo %9.4f ° pu IG1_Bp=[abs(IG1_B) angle(IG1_B)*180/pi];
\n',abs(V3),radtodeg(angle(V3)))
fprintf('V2= %9.4f angulo %9.4f ° pu IG1_C=IG1_0+a*I1g1+a^2*IG1_2;
\n',abs(V2),radtodeg(angle(V2))) IG1_Cp=[abs(IG1_C) angle(IG1_C)*180/pi];
fprintf('V1= %9.4f angulo %9.4f ° pu
\n',abs(V1),radtodeg(angle(V1))) IG1_ABC=[IG1_Ap;IG1_Bp;IG1_Cp]
%% c)CALCULAR LA CORRIENTE QUE CIRCULA EN CADA FASE DEL G1 fprintf('\n')
%% CALCULO DE POTENCIA DE CORTOCIRCUITO
disp('-----------------RED DE SECUENCIA POSITIVA------------')
disp('CALCULO DE POTENCIA DE CORTOCIRCUITO PARA CADA NODO')
z1=((1j+0.35j)*(1j+0.35j))/(1j+0.35j+1j+0.35j); Z1th5=(z1+z4+z5)*z2/(z1+z4+z5+z2);
z2=4+2j; Z2th5=(g1+g4+g5)*g2/(g1+g4+g5+g2);
z4=(0.1778j*0.1778j)/(0.1778j+0.1778j); %% BARRA 1
z5=(0.35j*0.35j)/(0.35j+0.35j); Z1th1_eq=(z2+z5+z4)*(1j+0.35j)/(z2+z5+z4+1j+0.35j)+j*0.35;
Zposit=(z1+z4)*(z2+z5)/(z1+z2+z4+z5); Z1th1=Z1th1_eq*1j/(Z1th1_eq+1j);

fprintf('Zpositivo= %9.6f + %9.4fj pu Z2th1_eq=(g2+g5+g4)*(Xg1_2+0.35j)/(g2+g5+g4+Xg1_2+0.35j)+j*0.3

\n',real(Zposit),imag(Zposit)) 5;
fprintf('\n') Z2th1=Z2th1_eq*Xg1_2/(Z2th1_eq+Xg1_2);
disp('-----------------RED DE SECUENCIA NEGATIVA-----------') Z0th1=Xg1_0;
Zt=Z0th1+Z1th1+Z2th1;
g1=((0.6j+0.35j)*(0.6j+0.35j))/(0.6j+0.35j+0.6j+0.35j);
g2=4+2j; Ia0_1=V1/Zt;
g4=(0.1778j*0.1778j)/(0.1778j+0.1778j); Ia1_1=Ia0_1; Ia2_1=Ia0_1;
g5=(0.35j*0.35j)/(0.35j+0.35j); I1=[Ia0_1;Ia1_1;Ia2_1];
 II=A*I1; RESULTADOS EJERCICIO4.4
PCC=3*abs(II(1))*SB; -----------------DATOS EN PU----------------------
fprintf(' PCC BARRA 1= %.4f MVA \n',PCC) Xgenerador 1 (positivo)= 1.0000j pu
%% BARRA 2 Xgenerador 1(negativo)= 0.6000j pu
Z2th1_eq=(z2+z5+z4)*(1j+0.35j)/(z2+z5+z4+1j+0.35j)+j*0.35; Xgenerador 1(cero)= 0.4000j pu
Z2th1=Z1th1_eq*1j/(Z1th1_eq+1j); X transf. 1= 0.3500j pu
X transf. 1(primario)= 0.0700j pu
X transf. 1(secundario)= 0.2800j pu
Z2th1_eq=(g2+g5+g4)*(Xg1_2+0.35j)/(g2+g5+g4+Xg1_2+0.35j)+j*0.3 X transf. 2= 0.3500j pu
5; X transf. 2(primario)= 0.0700j pu
Z2th1=Z2th1_eq*Xg1_2/(Z2th1_eq+Xg1_2); X transf. 2(secundario)= 0.2800j pu
Z0th1=Xg1_0; X linea(positivo)= 0.0000 + 0.1778j pu
Zt=Z0th1+Z1th1+Z2th1; X linea (cero)= 0.0000 + 0.5334j pu
Z carga= 4.0000 + 2.0000j pu
Ia0_1=V2/Zt;
Ia1_1=Ia0_1; Ia2_1=Ia0_1; ------------TENSIONES EN CADA BARRA----------------------
I1=[Ia0_1;Ia1_1;Ia2_1]; V5= 1.0000 angulo 0.0000 ° pu
II=A*I1; V4= 1.0181 angulo 1.9701 ° pu
%fprintf('I1_2= %9.4f \n',abs(II(1))); V3= 1.0277 angulo 2.9437 ° pu
PCC=3*abs(II(1))*SB; V2= 1.0476 angulo 4.8067 ° pu
fprintf(' PCC BARRA 2= %.4f MVA \n',PCC) V1= 1.0476 angulo 4.8067 ° pu
%% BARRA 3 -----------------RED DE SECUENCIA POSITIVA----------------
Zpositivo= 0.094742 + 0.6943j pu
Z1th3=(z2+z4+z5)*z1/(z2+z4+z5+z1);
Z2th3=(g2+g4+g5)*g1/(g2+g4+g5+g1); -----------------RED DE SECUENCIA NEGATIVA--------------
Z0th3=(XL1_0/2+Xs_T2/2)*(Xs_T1/2)/(XL1_0/2+Xs_T2/2+Xs_T1/2); Znegativo= 0.054121 + 0.5268j pu

Zt=Z0th3+Z1th3+Z2th3; -----------------RED DE SECUENCIA CERO----------------------

Zcero= 0.000000 + 0.1041j pu
Ia0_1=V3/Zt;
Ia1_1=Ia0_1; Ia2_1=Ia0_1; CORRIENTES DE SECUENCIA EN NODO 4
I1a=[Ia0_1;Ia1_1;Ia2_1]; Ia 0= 0.7634 angulo -81.6210 ° pu
II=A*I1a; Ia 1= 0.7634 angulo -81.6210 ° pu
PCC=3*abs(II(1))*SB; Ia 2= 0.7634 angulo -81.6210 ° pu
fprintf(' PCC BARRA 3= %.4f MVA \n',PCC) CORRIENTES DE FALLA FASICA EN NODO 4
%% BARRA 4 Ia= 2.2902 angulo -81.6210 ° pu
Zt=Zposit+Znegat+Zcero; Ib= 0.0000 angulo 0.0000 ° pu
Ic= 0.0000 angulo 0.0000 ° pu
Ia0_1=V4/Zt; CORRIENTES de SECUENCIA DE FALLA EN G1
Ia1_1=Ia0_1; Ia2_1=Ia0_1; I1g1= 0.3501 angulo -89.3915 ° pu
I1a=[Ia0_1;Ia1_1;Ia2_1]; I2g1= 0.3585 angulo -87.4863 ° pu
II=A*I1a; I0g1= 0.0000 angulo 0.0000 ° pu
PCC=3*abs(II(1))*SB;
fprintf(' PCC BARRA 4= %.4f MVA \n',PCC) CORRIENTES FASICAS DE FALLA EN G1
%% BARRA 5 IG1_ABC =
Z1th3=(z1+z4+z5)*z2/(z2+z4+z5+z1); 0.7085 -88.4276
Z2th3=(g1+g4+g5)*g2/(g2+g4+g5+g1); 0.3645 90.4343
Z0th3=Z1; 0.3441 92.7779
Zt=Z0th3+Z1th3+Z2th3;
CALCULO DE POTENCIA DE CORTOCIRCUITO PARA CADA NODO
Ia0_1=V5/Zt; PCC BARRA 1= 672.9785 MVA
Ia1_1=Ia0_1; Ia2_1=Ia0_1; PCC BARRA 2= 672.9785 MVA
I1a=[Ia0_1;Ia1_1;Ia2_1]; PCC BARRA 3= 784.6243 MVA
II=A*I1a; PCC BARRA 4= 687.0732 MVA
PCC=3*abs(II(1))*SB; PCC BARRA 5= 163.6166 MVA
fprintf(' PCC BARRA 5= %.4f MVA \n',PCC) >>