clc, clear all fprintf('Vcfalla = %9.4f angulo %9.

4f ° pu
Vb=220; \n',abs(Vcfalla),radtodeg(angle(Vcfalla)))
Sb=25;
Zb=Vb^2/Sb; Ib1=Ia1*a2;
Zf=4; Ic1=Ia1*a1;
Zfpu=Zf/Zb; %% graficas
figure
Xg1=j*0.125*(25/25)*(11/11)^2; compass(Va0,'b');
Xt1=j*0.105*(25/30)*(220/220)^2; hold on
XL=(2+6*j)/Zb; compass(Ia1,'r');
Xt2=j*0.10*(25/30)*(220/220)^2; hold on
zpc=(10.8/11)^2/conj((10+8*j)/Sb); compass(Ib1,'r');
hold on
disp('-----------------DATOS EN PU----------------------') compass(Ic1,'r');
fprintf('Xgenerador 1= %9.4fj pu \n',imag(Xg1))
fprintf('X transf. 1= %9.4fj pu \n',imag(Xt1)) hold on
fprintf('X linea= %9.4f + %9.4fj pu \n',real(XL),imag(XL)) compass(Vafalla,'y');
fprintf('X transf. 2= %9.4f j pu \n',imag(Xt2)) hold on
fprintf('Zcarga= %9.4f + %9.4fj pu \n',real(zpc),imag(zpc)) compass(Vbfalla,'y');
fprintf('\n') hold on
disp('-----------------RED DE SECUENCIA POSITIVA----------- compass(Vcfalla,'y');
-----------')
hold on
Z1=(Xg1+Xt1+XL); Z2=(Xt2+zpc); Vb1=a2*Va1;
Zposit=Z1*Z2/(Z1+Z2); Vc1=a1*Va1;
fprintf('Zpositivo= %9.6f + %9.4fj pu compass(Va1,'g');
\n',real(Zposit),imag(Zposit)) hold on
fprintf('\n') compass(Vb1,'g');
disp('-----------------RED DE SECUENCIA NEGATIVA----------- hold on
-----------') compass(Vc1,'g');
Znegat=Zposit;
fprintf('Znegativo= %9.6f + %9.4fj pu hold on
\n',real(Znegat),imag(Znegat)) Vb2=a2*Va2;
fprintf('\n') Vc2=a1*Va2;
disp('-----------------RED DE SECUENCIA CERO---------------
-------') compass(Va2,'m');
hold on
Z0linea=3*XL; compass(Vb2,'m');
Z0carga=3*zpc; Z3=Xt1+Z0linea; Z4=Xt2; hold on
Zcero=Z3*Z4/(Z3+Z4); compass(Vc2,'m');
fprintf('Zcero= %9.6f + %9.4fj pu
\n',real(Zcero),imag(Zcero)) title('tensiones y corrientes')
fprintf('\n') legend('Va0','Ia1','Ib1','Ic1','Vafalla','Vbfalla','Vcfalla
','Va1','Vb1','Vc1','Va2','Vb2','Vc2')
disp('-----------------I carga y V falla------------------ hold off
----')
Icarga=1/((11/11)^2/conj((10+8*j)/Sb)); -----------------DATOS EN PU----------------------
fprintf('Icarga= %9.4f %9.4f j pu
\n',real(Icarga),imag(Icarga)) Xgenerador 1= 0.1250j pu
%Voltage de falla en la barra 3 X transf. 1= 0.0875j pu
disp('-----------------------------------------------------
-------') X linea= 0.0010 + 0.0031j pu
Vfalla=1+Icarga*Xt2; X transf. 2= 0.0833 j pu
fprintf('Vfalla= %9.4f angulo %9.4f ° pu
\n',abs(Vfalla),radtodeg(angle(Vfalla))) Zcarga= 1.4695 + 1.1756j pu
Zfallapu=3*Zfpu;
Vprefalla=Vfalla; -----------------RED DE SECUENCIA POSITIVA----------------------
Ia1=Vprefalla/(Zposit+Znegat+Zcero+Zfallapu); Zpositivo= 0.016643 + 0.1996j pu
fprintf('Ia 1= %9.4f angulo %9.4f ° pu
\n',abs(Ia1),radtodeg(angle(Ia1)))
-----------------RED DE SECUENCIA NEGATIVA----------------------
Ifallapu=3*Ia1;
fprintf('Ifalla pu= %9.4f angulo %9.4f ° pu Znegativo= 0.016643 + 0.1996j pu
\n',abs(Ifallapu),radtodeg(angle(Ifallapu)))
Va1=Vfalla-Zposit*Ia1;
fprintf('Va1= %9.4f angulo %9.4f ° pu -----------------RED DE SECUENCIA CERO----------------------
\n',abs(Va1),radtodeg(angle(Va1))) Zcero= 0.000663 + 0.0448j pu
Va2=0-Znegat*Ia1;
fprintf('Va2= %9.4f angulo %9.4f ° pu -----------------I carga y V falla----------------------
\n',abs(Va2),radtodeg(angle(Va2))) Icarga= 0.4000 -0.3200 j pu
Va0=0-Zcero*Ia1;
fprintf('Va3= %9.4f angulo %9.4f ° pu ------------------------------------------------------------
\n',abs(Va0),radtodeg(angle(Va0))) Vfalla= 1.0272 angulo 1.8596 ° pu
Vafalla=Va1+Va2+Va0; Ia 1= 2.3037 angulo -82.9744 ° pu
fprintf('Vafalla= %9.4f angulo %9.4f ° pu Ifalla pu= 6.9112 angulo -82.9744 ° pu
\n',abs(Vafalla),radtodeg(angle(Vafalla)))
Va1= 0.5657 angulo 1.5329 ° pu
Vafallaz=Zfpu*Ifallapu; Va2= 0.4615 angulo -177.7399 ° pu
fprintf('Vafalla con Zf= %9.4f angulo %9.4f ° pu
\n',abs(Vafallaz),radtodeg(angle(Vafallaz))) Va3= 0.1032 angulo -173.8226 ° pu
fprintf('\n') Vafalla= 0.0143 angulo -82.9744 ° pu
a1=-1/2+i*sqrt(3)/2;
a2=-1/2-i*sqrt(3)/2; Vafalla con Zf= 0.0143 angulo -82.9744 ° pu
Vbfalla=a2*Va1+a1*Va2+Va0;
fprintf('Vbfalla = %9.4f angulo %9.4f ° pu
\n',abs(Vbfalla),radtodeg(angle(Vbfalla))) Vbfalla = 0.9075 angulo -97.9753 ° pu
Vb0=Va0; Vc0=Va0; Vcfalla = 0.8985 angulo 101.7935 ° pu
Vcfalla=a1*Va1+a2*Va2+Va0;
>>