clc
variasejecuciones=1;
while variasejecuciones==1
clear
% --------------------------------------------------
salida=input('Pantalla=1, Archivo=2 :');
if salida==2
cadena=input('Poner NombreArchivo.txt:','s');
cadena1=strcat('Salida\',cadena);
end
if salida==1
fid=1;
elseif salida==2
fid=fopen(cadena1,'w+');
end
% --------------------------------------------------
fprintf(fid,'PC124\n');
fprintf(fid,'Esquema Linea Doble SR\n');
fprintf(fid,'Unominal kV: 400,220,132,66\n');
fprintf(fid,'Impedancia de Transferencia: Si,No\n');
fprintf(fid,'Compensacion 21: Si,No\n');
fprintf(fid,'Falta Simple m, o Doble m,p\n');
fprintf(fid,'Barras Ficticias Iniciales m,p\n');
fprintf(fid,'Esquema Normal, Sm fuera, Rm fuera\n');
fprintf(fid,'\n');
% --------------------------------------------------
%PC24
datos=input('datos: standard=1,input=2: ');
if datos==1
Unominal=220;
kE=1;
L=220;
imptransfer=2;
Compens21=2;
Sccs=5000;
Sccr=5000;
Kr0=1;
esquema=1;
% -----------------------------------------------------------------
else
Unominal=input('Unominal kV: 400,220,132,66 :');
kE=input('kE=E/Enominal: ');
L=input('L km: ');
Sccs=input('Sccs(MVA): ');
Sccr=input('Sccr(MVA): ');
Kr0=input('Kr0: ');
imptransfer=input('zsr Si=1, No=2 :');
if imptransfer==1
zsr1=input('zsr1=zsr2(Complejo r+xi en Ohms): ');
zsr2=zsr1;
Ksr0=input('Ksr0=zsr0/zsr1: ');
zsr0=Ksr0*zsr1;
end
Compens21=input('Compensacion 21: Si=1, No=2 :');
esquema=input('Esquema: Normal=1, Sm Abierta=2, Rm Abierta=3 :');
end
% -----------------------------------------------------------------
m=0.5; p=0.5;
% -----------------------------------------------------------------
k0=input('k0=zs0/zr0: ');
TipoFalta=input('Falta Simple m=1, Doble m,p=2 :');
fprintf(fid,'Falta m en Linea SR1, entre S y R\n');
m=input('m =(0...1)= :');
fprintf(fid,'\n');
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
% -----------------------------------------------------------------
Ubase=Unominal;
Vbase=Ubase/(sqrt(3));
MVAbase=100;
Ibase=MVAbase/(3*Vbase);
zbase=Vbase/Ibase;
if Unominal==400
zl1u=(0.0346+0.3208i)/zbase; zl2u=zl1u; zl0u=(0.1828+1.1482i)/zbase;
zlM0u=(0.1482+0.7383i)/zbase;
elseif Unominal==220
zl1u=(0.0829+0.4135i)/zbase; zl2u=zl1u; zl0u=(0.2311+1.3085i)/zbase;
zlM0u=(0.1482+0.8031i)/zbase;
elseif Unominal==132
zl1u=(0.1342+0.4042i)/zbase; zl2u=zl1u; zl0u=(0.2824+1.3718i)/zbase;
zlM0u=(0.1482+0.8947i)/zbase;
elseif Unominal==66
zl1u=(0.2156+0.3952i)/zbase; zl2u=zl1u; zl0u=(0.3636+1.4278i)/zbase;
zlM0u=(0.1482+0.9288i)/zbase;
end
if imptransfer==1
zsr1=zsr1/zbase;
zsr2=zsr1;
zsr0=zsr0/zbase;
end
% -----------------------------------------------------------------
Sm=zl1u*m*L; Rm=zl1u*(1-m)*L; % pu
SSm=Sm*zbase; RRm=Rm*zbase; % Ohm
if TipoFalta==2
fprintf(fid,'Falta p en Linea SR2, entre S y R\n');
p=input('p =(0...1)= :');
fprintf(fid,'\n');
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
fprintf(fid,'\n');
end
Sp=zl1u*p*L; Rp=zl1u*(1-p)*L; % pu
SSp=Sp*zbase; RRp=Rp*zbase; % Ohm
end
% -----------------------------------------------------------------
U=kE*Unominal; % kV
E=U/sqrt(3)+0i; % kV
E=E/Vbase; % E pu
k=(zl0u-zl1u)/(3*zl1u);
Km=zlM0u/zl1u;
zs1=(MVAbase/Sccs)*i; % pu
zr1=(MVAbase/Sccr)*i; % pu
zs2=zs1; % pu
zr2=zr1; % pu
zr0=Kr0*zr1; % pu
zs0=k0*zr0; % pu
zl0=zl0u*L; zl1=zl1u*L; zl2=zl2u*L; zlM0=zlM0u*L; % pu
% -----------------------------------------------------------------
% CALCULOS
% -----------------------------------------------------------------
if imptransfer==2&esquema==1
[Minc]=subMinc1(1);
[Minca]=subMinca1(1);
[z0]=subz0_1(zs0,zr0,zl0,zlM0,m,p);
[z1]=subz1_1(zs1,zr1,zl1,m,p);
[z2]=subz2_1(zs2,zr2,zl2,m,p);
elseif imptransfer==2&esquema==2
[Minc]=subMinc2(1);
[Minca]=subMinca2(1);
[z1]=subz1_2(zs1,zr1,zl1,m,p);
[z2]=subz2_2(zs2,zr2,zl2,m,p);
[z0]=subz0_2(zs0,zr0,zl0,zlM0,m,p);
elseif imptransfer==2&esquema==3
[Minc]=subMinc3(1)
[Minca]=subMinca3(1)
[z1]=subz1_3(zs1,zr1,zl1,m,p)
[z2]=subz2_3(zs2,zr2,zl2,m,p)
[z0]=subz0_3(zs0,zr0,zl0,zlM0,m,p)
elseif imptransfer==1&esquema==1
[Minc]=subMinc1sr(1);
[Minca]=subMinca1sr(1);
[z0]=subz0_1sr(zs0,zsr0,zr0,zl0,zlM0,m,p);
[z1]=subz1_1sr(zs1,zsr1,zr1,zl1,m,p);
[z2]=subz2_1sr(zs2,zsr2,zr2,zl2,m,p);
elseif imptransfer==1&esquema==2
[Minc]=subMinc2sr(1);
[Minca]=subMinca2sr(1);
[z1]=subz1_2sr(zs1,zsr1,zr1,zl1,m,p);
[z2]=subz2_2sr(zs2,zsr2,zr2,zl2,m,p);
[z0]=subz0_2sr(zs0,zsr0,zr0,zl0,zlM0,m,p)
elseif imptransfer==1&esquema==3
[Minc]=subMinc3sr(1);
[Minca]=subMinca3sr(1);
[z1]=subz1_3sr(zs1,zsr1,zr1,zl1,m,p);
[z2]=subz2_3sr(zs2,zsr2,zr2,zl2,m,p);
[z0]=subz0_3sr(zs0,zsr0,zr0,zl0,zlM0,m,p);
end
% MATRIZ DE ADMITANCIA DE BARRAS
[Y0]=subY0(Minc,z0); % pu
[Y1]=subY1(Minc,z1); % pu
[Y2]=subY2(Minc,z2); % pu
% ZBUS(0,1,2)
Z0=inv(Y0); % pu
Z1=inv(Y1); % pu
Z2=inv(Y2); % pu
pfaltam=3;
pfaltap=4;
if TipoFalta==1
% CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN EL PUNTO m DE FALTA
[Im0,Im1,Im2,Vm0,Vm1,Vm2,faltam,Rfm,Rgm]=ma2m(Z0,Z1,Z2,E,zbase,pfaltam); % pu
elseif TipoFalta==2
% CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p DE FALTA
[Im0,Im1,Im2,Vm0,Vm1,Vm2,Ip0,Ip1,Ip2,Vp0,Vp1,Vp2,faltam,faltap,Rfm,Rgm,Rfp,Rgp]=...
ma4mp(Z0,Z1,Z2,E,zbase,pfaltam,pfaltap); % pu
end
if faltam==1 cfaltam='ag';elseif faltam==2 cfaltam='bg';elseif faltam==3 cfaltam='cg';
elseif faltam==4 cfaltam='ab';elseif faltam==5 cfaltam='bc';elseif faltam==6 cfaltam='ca';
elseif faltam==7 cfaltam='abg';elseif faltam==8 cfaltam='bcg';elseif faltam==9 cfaltam='cag';
elseif faltam==10 cfaltam='abc';elseif faltam==11 cfaltam='abcg';
end
if TipoFalta==2
if faltap==1 cfaltap='ag';elseif faltap==2 cfaltap='bg';elseif faltap==3 cfaltap='cg';
elseif faltap==4 cfaltap='ab';elseif faltap==5 cfaltap='bc';elseif faltap==6 cfaltap='ca';
elseif faltap==7 cfaltap='abg';elseif faltap==8 cfaltap='bcg';elseif faltap==9 cfaltap='cag';
elseif faltap==10 cfaltap='abc';elseif faltap==11 cfaltap='abcg';
end
end
fprintf(fid,'\n');
% IBUS(0,1,2)
I0=zeros(4,1);
I1=zeros(4,1);
I2=zeros(4,1);
I0(pfaltam,1)=-Im0;
I1(pfaltam,1)=-Im1;
I2(pfaltam,1)=-Im2;
if TipoFalta==2
I0(pfaltap,1)=-Ip0;
I1(pfaltap,1)=-Ip1;
I2(pfaltap,1)=-Ip2;
end
% VBUS(0,1,2)
V0=Z0*I0;
V1=Z1*I1+E; V1a=[V1;E];
V2=Z2*I2;
% Ie(0,1,2) DE LOS ELEMENTOS
[Ie0]=subIe0(Minc,V0,z0);
[Ie1]=subIe1(Minca,V1a,z1);
[Ie2]=subIe2(Minc,V2,z2);
% ASIGNACIONES V,I
Vbase=Vbase*1000;
Ibase=Ibase*1000;
VVm0=Vm0*Vbase;VVm1=Vm1*Vbase;VVm2=Vm2*Vbase;
if TipoFalta==2
VVp0=Vp0*Vbase;VVp1=Vp1*Vbase;VVp2=Vp2*Vbase;
end
Vs0=V0(1);Vs1=V1(1);Vs2=V2(1); VVs0=Vs0*Vbase;VVs1=Vs1*Vbase;VVs2=Vs2*Vbase;
Vr0=V0(2);Vr1=V1(2);Vr2=V2(2); VVr0=Vr0*Vbase;VVr1=Vr1*Vbase;VVr2=Vr2*Vbase;
IIm0=Im0*Ibase;IIm1=Im1*Ibase;IIm2=Im2*Ibase;
if TipoFalta==2
IIp0=Ip0*Ibase;IIp1=Ip1*Ibase;IIp2=Ip2*Ibase;
end
Is0=Ie0(1);Is1=Ie1(1);Is2=Ie2(1); IIs0=Is0*Ibase;IIs1=Is1*Ibase;IIs2=Is2*Ibase;
Ir0=Ie0(2);Ir1=Ie1(2);Ir2=Ie2(2); IIr0=Ir0*Ibase;IIr1=Ir1*Ibase;IIr2=Ir2*Ibase;
Ism0=Ie0(3);Ism1=Ie1(3);Ism2=Ie2(3); IIsm0=Ism0*Ibase;IIsm1=Ism1*Ibase;IIsm2=Ism2*Ibase;
Isp0=Ie0(4);Isp1=Ie1(4);Isp2=Ie2(4); IIsp0=Isp0*Ibase;IIsp1=Isp1*Ibase;IIsp2=Isp2*Ibase;
Irm0=Ie0(5);Irm1=Ie1(5);Irm2=Ie2(5); IIrm0=Irm0*Ibase;IIrm1=Irm1*Ibase;IIrm2=Irm2*Ibase;
Irp0=Ie0(6);Irp1=Ie1(6);Irp2=Ie2(6); IIrp0=Irp0*Ibase;IIrp1=Irp1*Ibase;IIrp2=Irp2*Ibase;
if imptransfer==1
if esquema==1
Izsr0=Ie0(7);Izsr1=Ie1(7);Izsr2=Ie2(7); IIzsr0=Izsr0*Ibase;IIzsr1=Izsr1*Ibase;IIzsr2=Izsr2*Ibase;
elseif esquema==2|esquema==3
Izsr0=Ie0(8);Izsr1=Ie1(8);Izsr2=Ie2(8); IIzsr0=Izsr0*Ibase;IIzsr1=Izsr1*Ibase;IIzsr2=Izsr2*Ibase;
end
end
[Ima,Imb,Imc]=abc(Im0,Im1,Im2); IIma=Ima*Ibase;IImb=Imb*Ibase;IImc=Imc*Ibase;
if TipoFalta==2
[Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2); IIpa=Ipa*Ibase;IIpb=Ipb*Ibase;IIpc=Ipc*Ibase;
end
[Isa,Isb,Isc]=abc(Is0,Is1,Is2); IIsa=Isa*Ibase;IIsb=Isb*Ibase;IIsc=Isc*Ibase;
[Ira,Irb,Irc]=abc(Ir0,Ir1,Ir2); IIra=Ira*Ibase;IIrb=Irb*Ibase;IIrc=Irc*Ibase;
[Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2); IIsma=Isma*Ibase;IIsmb=Ismb*Ibase;IIsmc=Ismc*Ibase;
[Ispa,Ispb,Ispc]=abc(Isp0,Isp1,Isp2); IIspa=Ispa*Ibase;IIspb=Ispb*Ibase;IIspc=Ispc*Ibase;
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2); IIrma=Irma*Ibase;IIrmb=Irmb*Ibase;IIrmc=Irmc*Ibase;
[Irpa,Irpb,Irpc]=abc(Irp0,Irp1,Irp2); IIrpa=Irpa*Ibase;IIrpb=Irpb*Ibase;IIrpc=Irpc*Ibase;
if imptransfer==1
[Izsra,Izsrb,Izsrc]=abc(Izsr0,Izsr1,Izsr2); IIzsra=Izsra*Ibase;IIzsrb=Izsrb*Ibase;IIzsrc=Izsrc*Ibase;
end
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2); VVma=Vma*Vbase;VVmb=Vmb*Vbase;VVmc=Vmc*Vbase;
if TipoFalta==2
[Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2); VVpa=Vpa*Vbase;VVpb=Vpb*Vbase;VVpc=Vpc*Vbase;
end
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2); VVsa=Vsa*Vbase;VVsb=Vsb*Vbase;VVsc=Vsc*Vbase;
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2); VVra=Vra*Vbase;VVrb=Vrb*Vbase;VVrc=Vrc*Vbase;
% PROTECCION DE DISTANCIA (21)
switch(Compens21)
case{1},
% COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k,Km,Isp0);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k,Km,Ism0);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k,Km,Irp0);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k,Km,Irm0);
case{2},
% NO COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
PDist(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
PDist(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k);
end
% 67N
[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0); VVspol=Vspol*Vbase;IIsmop=Ismop*Ibase;
[Vspol,Ispop,AngDirSp]=PDire(Vs0,Isp0); IIspop=Ispop*Ibase;
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0); VVrpol=Vrpol*Vbase;IIrmop=Irmop*Ibase;
[Vrpol,Irpop,AngDirRp]=PDire(Vr0,Irp0); IIrpop=Irpop*Ibase;
% ******************************************************************
% *************************** SALIDA ******************************
% ******************************************************************
fprintf(fid,'\n');
fprintf(fid,'***************************************************\n');
fprintf(fid,'Unominal =%-4.0f(kV)\n',Unominal);
fprintf(fid,'kE=E/Enominal =%-4.0f\n',kE);
fprintf(fid,'L=%-5.1f(km)\n',L);
fprintf(fid,'Ibase=%-8.2f(kA)\n',Ibase*1000);
fprintf(fid,'Tipo de Falta, Simple(1), Doble(2) =%-1.0f\n',TipoFalta);
fprintf(fid,'Esquema, Normal(1), Sm fuera(2), Rm fuera(3) =%-1.0f\n',esquema);
fprintf(fid,'Impedancia de Transferencia, Si(1), No(2) =%-1.0f\n',imptransfer);
fprintf(fid,'Compensacion 21, Si(1), No(2) =%-1.0f\n',Compens21);
if imptransfer==1
fprintf('zsr1=zsr2=%-7.2f+i(%-7.2f)Ohm',real(zsr1),imag(zsr1));
fprintf('\tzsr0=%-7.2f+i(%-7.2f)Ohm\n',real(zsr0),imag(zsr0));
end
fprintf(fid,'Sccs =%-8.0f(MVA)\n',Sccs);
fprintf(fid,'Sccr =%-8.0f(MVA)\n',Sccr);
fprintf(fid,'Kr0 =%-4.1f\n',Kr0);
fprintf(fid,'k0 =%-4.1f\n',k0);
fprintf(fid,'\n');
fprintf(fid,'m =%-3.2f\t',m);
fprintf(fid,'Falta:%s\t',cfaltam);
fprintf(fid,'Rfm(Ohm):%-10.2f',Rfm*zbase);
if (faltam==7)|(faltam==8)|(faltam==9)|(faltam==11)
fprintf(fid,'\tRgm(Ohm):%-10.2f',Rgm*zbase);
end
if TipoFalta==2
fprintf(fid,'\n');
fprintf(fid,'p =%-3.2f\t',p);
fprintf(fid,'Falta:%s\t',cfaltap);
fprintf(fid,'Rfp(Ohm):%-10.2f',Rfp*zbase);
if (faltap==7)|(faltap==8)|(faltap==9)|(faltap==11)
fprintf(fid,'\tRgp(Ohm):%-10.2f',Rgp*zbase);
end
end
fprintf(fid,'\n');
fprintf(fid,'***************************************************\n');
variassalidas=1;
while variassalidas==1
fprintf(fid,'\n');
output=input('Salida: pu(cart)=1, pu(polar)=2, Ohm-Volt-Amp(cart)=3, Ohm-Volt-Amp(polar)=4 :');
if salida==2
fprintf(fid,'Salida: pu(cart)=1, pu(polar)=2, Ohm-Volt-Amp(cart)=3, Ohm-Volt-Amp(polar)=4 : %1.0f',output);
fprintf(fid,'\n');
end
switch(output)
% OPCIONES DE SALIDA
% --------------------------------
% 21 67N V(0,1,2) V(a,b,c)
% I(0,1,2) I(a,b,c)
% 1: pu(cartesiana) *** ***
% 2: pu(polar) *** *** ***
% 3: OVA(cartesiana) *** *** ***
% 4: OVA(polar) *** *** *** ***
% **********************************************************************
case{1}, % pu(cartesiana): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
% ----------------- -----------------
% **********************************************************************
% -----------------------------------------------------------------------
% ** C1 ************** SALIDA V(0,1,2),I(0,1,2) pu cartesiana ***********
% -----------------------------------------------------------------------
% Vm,Vp_2,Vs,Vr
% Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr
% -----------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA V(0,1,2),I(0,1,2) pu cartesiana\n\n');
fprintf(fid,'\n');
fprintf(fid,'Vm0 =%8.4f+i(%-8.4f)pu\t',real(Vm0),imag(Vm0));
fprintf(fid,'Vm1 =%8.4f+i(%-8.4f)pu\t',real(Vm1),imag(Vm1));
fprintf(fid,'Vm2 =%8.4f+i(%-8.4f)pu\t',real(Vm2),imag(Vm2));
fprintf(fid,'\n');
if TipoFalta==2
fprintf(fid,'Vp0 =%8.4f+i(%-8.4f)pu\t',real(Vp0),imag(Vp0));
fprintf(fid,'Vp1 =%8.4f+i(%-8.4f)pu\t',real(Vp1),imag(Vp1));
fprintf(fid,'Vp2 =%8.4f+i(%-8.4f)pu\t',real(Vp2),imag(Vp2));
fprintf(fid,'\n');
end
fprintf(fid,'Vs0 =%8.4f+i(%-8.4f)pu\t',real(Vs0),imag(Vs0));
fprintf(fid,'Vs1 =%8.4f+i(%-8.4f)pu\t',real(Vs1),imag(Vs1));
fprintf(fid,'Vs2 =%8.4f+i(%-8.4f)pu\t',real(Vs2),imag(Vs2));
fprintf(fid,'\n');
fprintf(fid,'Vr0 =%8.4f+i(%-8.4f)pu\t',real(Vr0),imag(Vr0));
fprintf(fid,'Vr1 =%8.4f+i(%-8.4f)pu\t',real(Vr1),imag(Vr1));
fprintf(fid,'Vr2 =%8.4f+i(%-8.4f)pu\t',real(Vr2),imag(Vr2));
fprintf(fid,'\n\n');
% ----------------------------------------------------------
fprintf(fid,'Im0 =%8.4f+i(%-8.4f)pu\t',real(Im0),imag(Im0));
fprintf(fid,'Im1 =%8.4f+i(%-8.4f)pu\t',real(Im1),imag(Im1));
fprintf(fid,'Im2 =%8.4f+i(%-8.4f)pu\t',real(Im2),imag(Im2));
fprintf(fid,'\n');
if TipoFalta==2
fprintf(fid,'Ip0 =%8.4f+i(%-8.4f)pu\t',real(Ip0),imag(Ip0));
fprintf(fid,'Ip1 =%8.4f+i(%-8.4f)pu\t',real(Ip1),imag(Ip1));
fprintf(fid,'Ip2 =%8.4f+i(%-8.4f)pu\t',real(Ip2),imag(Ip2));
fprintf(fid,'\n');
end
fprintf(fid,'Is0 =%8.4f+i(%-8.4f)pu\t',real(Is0),imag(Is0));
fprintf(fid,'Is1 =%8.4f+i(%-8.4f)pu\t',real(Is1),imag(Is1));
fprintf(fid,'Is2 =%8.4f+i(%-8.4f)pu\t',real(Is2),imag(Is2));
fprintf(fid,'\n');
fprintf(fid,'Ir0 =%8.4f+i(%-8.4f)pu\t',real(Ir0),imag(Ir0));
fprintf(fid,'Ir1 =%8.4f+i(%-8.4f)pu\t',real(Ir1),imag(Ir1));
fprintf(fid,'Ir2 =%8.4f+i(%-8.4f)pu\t',real(Ir2),imag(Ir2));
fprintf(fid,'\n');
if esquema==1|esquema==3
fprintf(fid,'Ism0 =%8.4f+i(%-8.4f)pu\t',real(Ism0),imag(Ism0));
fprintf(fid,'Ism1 =%8.4f+i(%-8.4f)pu\t',real(Ism1),imag(Ism1));
fprintf(fid,'Ism2 =%8.4f+i(%-8.4f)pu\t',real(Ism2),imag(Ism2));
fprintf(fid,'\n');
end
if TipoFalta==2
fprintf(fid,'Isp0 =%8.4f+i(%-8.4f)pu\t',real(Isp0),imag(Isp0));
fprintf(fid,'Isp1 =%8.4f+i(%-8.4f)pu\t',real(Isp1),imag(Isp1));
fprintf(fid,'Isp2 =%8.4f+i(%-8.4f)pu\t',real(Isp2),imag(Isp2));
fprintf(fid,'\n');
elseif TipoFalta==1
fprintf(fid,'Isr0 =%8.4f+i(%-8.4f)pu\t',real(Isp0),imag(Isp0));
fprintf(fid,'Isr1 =%8.4f+i(%-8.4f)pu\t',real(Isp1),imag(Isp1));
fprintf(fid,'Isr2 =%8.4f+i(%-8.4f)pu\t',real(Isp2),imag(Isp2));
fprintf(fid,'\n');
end
if esquema==1|esquema==2
fprintf(fid,'Irm0 =%8.4f+i(%-8.4f)pu\t',real(Irm0),imag(Irm0));
fprintf(fid,'Irm1 =%8.4f+i(%-8.4f)pu\t',real(Irm1),imag(Irm1));
fprintf(fid,'Irm2 =%8.4f+i(%-8.4f)pu\t',real(Irm2),imag(Irm2));
fprintf(fid,'\n');
end
if TipoFalta==2
fprintf(fid,'Irp0 =%8.4f+i(%-8.4f)pu\t',real(Irp0),imag(Irp0));
fprintf(fid,'Irp1 =%8.4f+i(%-8.4f)pu\t',real(Irp1),imag(Irp1));
fprintf(fid,'Irp2 =%8.4f+i(%-8.4f)pu\t',real(Irp2),imag(Irp2));
fprintf(fid,'\n');
elseif TipoFalta==1
fprintf(fid,'Irs0 =%8.4f+i(%-8.4f)pu\t',real(Irp0),imag(Irp0));
fprintf(fid,'Irs1 =%8.4f+i(%-8.4f)pu\t',real(Irp1),imag(Irp1));
fprintf(fid,'Irs2 =%8.4f+i(%-8.4f)pu\t',real(Irp2),imag(Irp2));
fprintf(fid,'\n');
end
if imptransfer==1
fprintf(fid,'Izsr0=%8.4f+i(%-8.4f)pu\t',real(Izsr0),imag(Izsr0));
fprintf(fid,'Izsr1=%8.4f+i(%-8.4f)pu\t',real(Izsr1),imag(Izsr1));
fprintf(fid,'Izsr2=%8.4f+i(%-8.4f)pu\t',real(Izsr2),imag(Izsr2));
fprintf(fid,'\n');
end
% ------------------------------------------------------------------
% ** C1 *********** SALIDA V(a,b,c),I(a,b,c) pu cartesiana *********
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vr
% Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA V(a,b,c),I(a,b,c) pu cartesiana\n\n');
fprintf(fid,'\n');
fprintf(fid,'Vma =%8.4f+i(%-8.4f)pu\t',real(Vma),imag(Vma));
fprintf(fid,'Vmb =%8.4f+i(%-8.4f)pu\t',real(Vmb),imag(Vmb));
fprintf(fid,'Vmc =%8.4f+i(%-8.4f)pu\t',real(Vmc),imag(Vmc));
if TipoFalta==2
fprintf(fid,'\n');
fprintf(fid,'Vpa =%8.4f+i(%-8.4f)pu\t',real(Vpa),imag(Vpa));
fprintf(fid,'Vpb =%8.4f+i(%-8.4f)pu\t',real(Vpb),imag(Vpb));
fprintf(fid,'Vpc =%8.4f+i(%-8.4f)pu\t',real(Vpc),imag(Vpc));
end
fprintf(fid,'\n');
fprintf(fid,'Vsa =%8.4f+i(%-8.4f)pu\t',real(Vsa),imag(Vsa));
fprintf(fid,'Vsb =%8.4f+i(%-8.4f)pu\t',real(Vsb),imag(Vsb));
fprintf(fid,'Vsc =%8.4f+i(%-8.4f)pu\t',real(Vsc),imag(Vsc));
fprintf(fid,'\n');
fprintf(fid,'Vra =%8.4f+i(%-8.4f)pu\t',real(Vra),imag(Vra));
fprintf(fid,'Vrb =%8.4f+i(%-8.4f)pu\t',real(Vrb),imag(Vrb));
fprintf(fid,'Vrc =%8.4f+i(%-8.4f)pu\t',real(Vrc),imag(Vrc));
fprintf(fid,'\n\n');
fprintf(fid,'Ima =%8.4f+i(%-8.4f)pu\t',real(Ima),imag(Ima));
fprintf(fid,'Imb =%8.4f+i(%-8.4f)pu\t',real(Imb),imag(Imb));
fprintf(fid,'Imc =%8.4f+i(%-8.4f)pu\t',real(Imc),imag(Imc));
if TipoFalta==2
fprintf(fid,'\n');
fprintf(fid,'Ipa =%8.4f+i(%-8.4f)pu\t',real(Ipa),imag(Ipa));
fprintf(fid,'Ipb =%8.4f+i(%-8.4f)pu\t',real(Ipb),imag(Ipb));
fprintf(fid,'Ipc =%8.4f+i(%-8.4f)pu\t',real(Ipc),imag(Ipc));
end
fprintf(fid,'\n');
fprintf(fid,'Isa =%8.4f+i(%-8.4f)pu\t',real(Isa),imag(Isa));
fprintf(fid,'Isb =%8.4f+i(%-8.4f)pu\t',real(Isb),imag(Isb));
fprintf(fid,'Isc =%8.4f+i(%-8.4f)pu\t',real(Isc),imag(Isc));
fprintf(fid,'\n');
fprintf(fid,'Ira =%8.4f+i(%-8.4f)pu\t',real(Ira),imag(Ira));
fprintf(fid,'Irb =%8.4f+i(%-8.4f)pu\t',real(Irb),imag(Irb));
fprintf(fid,'Irc =%8.4f+i(%-8.4f)pu\t',real(Irc),imag(Irc));
fprintf(fid,'\n');
if esquema==1|esquema==3
fprintf(fid,'Isma =%8.4f+i(%-8.4f)pu\t',real(Isma),imag(Isma));
fprintf(fid,'Ismb =%8.4f+i(%-8.4f)pu\t',real(Ismb),imag(Ismb));
fprintf(fid,'Ismc =%8.4f+i(%-8.4f)pu\t',real(Ismc),imag(Ismc));
fprintf(fid,'\n');
end
if TipoFalta==2
fprintf(fid,'Ispa =%8.4f+i(%-8.4f)pu\t',real(Ispa),imag(Ispa));
fprintf(fid,'Ispb =%8.4f+i(%-8.4f)pu\t',real(Ispb),imag(Ispb));
fprintf(fid,'Ispc =%8.4f+i(%-8.4f)pu\t',real(Ispc),imag(Ispc));
fprintf(fid,'\n');
elseif TipoFalta==1
fprintf(fid,'Isra =%8.4f+i(%-8.4f)pu\t',real(Ispa),imag(Ispa));
fprintf(fid,'Isrb =%8.4f+i(%-8.4f)pu\t',real(Ispb),imag(Ispb));
fprintf(fid,'Isrc =%8.4f+i(%-8.4f)pu\t',real(Ispc),imag(Ispc));
fprintf(fid,'\n');
end
if esquema==1|esquema==2
fprintf(fid,'Irma =%8.4f+i(%-8.4f)pu\t',real(Irma),imag(Irma));
fprintf(fid,'Irmb =%8.4f+i(%-8.4f)pu\t',real(Irmb),imag(Irmb));
fprintf(fid,'Irmc =%8.4f+i(%-8.4f)pu\t',real(Irmc),imag(Irmc));
fprintf(fid,'\n');
end
if TipoFalta==2
fprintf(fid,'Irpa =%8.4f+i(%-8.4f)pu\t',real(Irpa),imag(Irpa));
fprintf(fid,'Irpb =%8.4f+i(%-8.4f)pu\t',real(Irpb),imag(Irpb));
fprintf(fid,'Irpc =%8.4f+i(%-8.4f)pu\t',real(Irpc),imag(Irpc));
fprintf(fid,'\n');
elseif TipoFalta==1
fprintf(fid,'Irsa =%8.4f+i(%-8.4f)pu\t',real(Irpa),imag(Irpa));
fprintf(fid,'Irsb =%8.4f+i(%-8.4f)pu\t',real(Irpb),imag(Irpb));
fprintf(fid,'Irsc =%8.4f+i(%-8.4f)pu\t',real(Irpc),imag(Irpc));
fprintf(fid,'\n');
end
if imptransfer==1
fprintf(fid,'Izsra=%8.4f+i(%-8.4f)pu\t',real(Izsra),imag(Izsra));
fprintf(fid,'Izsrb=%8.4f+i(%-8.4f)pu\t',real(Izsrb),imag(Izsrb));
fprintf(fid,'Izsrc=%8.4f+i(%-8.4f)pu\t',real(Izsrc),imag(Izsrc));
fprintf(fid,'\n\n');
end
% ******************************************************************
case{2}, % pu(polar): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
% --- ----------------- ----------------
% ******************************************************************
% ------------------------------------------------------------------
% ** C2 ******************* SALIDA 67N pu polar ********************
% ------------------------------------------------------------------
% Sm(e1,e3),Sp_2,Sr_1,Rm(e1,e2),Rp_2,Rs_1
% ------------------------------------------------------------------
fprintf(fid,'SALIDA 67N, pu, polar\n');
fprintf(fid,'\n');
if (esquema==1)|(esquema==3)
fprintf(fid,'Sm_67N:\tVspol=%8.4f\tIsmop=%8.4f\tAngDirSm=%5.1f\n',Vspol,Ismop,AngDirSm);
end
if TipoFalta==2
fprintf(fid,'Sp_67N:\tVspol=%8.4f\tIspop=%8.4f\tAngDirSp=%5.1f\n',Vspol,Ispop,AngDirSp);
end
if TipoFalta==1
fprintf(fid,'Sr_67N:\tVspol=%8.4f\tIsrop=%8.4f\tAngDirSr=%5.1f\n',Vspol,Ispop,AngDirSp);
end
if (esquema==1)|(esquema==2)
fprintf(fid,'Rm_67N:\tVrpol=%8.4f\tIrmop=%8.4f\tAngDirRm=%5.1f\n',Vrpol,Irmop,AngDirRm);
end
if TipoFalta==2
fprintf(fid,'Rp_67N:\tVrpol=%8.4f\tIrpop=%8.4f\tAngDirRp=%5.1f\n',Vrpol,Irpop,AngDirRp);
end
if TipoFalta==1
fprintf(fid,'Rs_67N:\tVrpol=%8.4f\tIrsop=%8.4f\tAngDirRs=%5.1f\n',Vrpol,Irpop,AngDirRp);
end
fprintf(fid,'\n');
% ------------------------------------------------------------------
% ** C2 ************* SALIDA V(0,1,2), I(0,1,2) pu polar ***********
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vr
% Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr
% ------------------------------------------------------------------
fprintf(fid,'SALIDA V(0,1,2), I(0,1,2) pu polar\n');
fprintf(fid,'\n');
[mod,arg]=mod_arg(Vm0);fprintf(fid,'Vm0 =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vm1);fprintf(fid,'Vm1 =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vm2);fprintf(fid,'Vm2 =%8.4f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
if TipoFalta==2
fprintf(fid,'SALIDA V(0,1,2), I(0,1,2) pu Polar\n');
[mod,arg]=mod_arg(Vp0);fprintf(fid,'Vp0 =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vp1);fprintf(fid,'Vp1 =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vp2);fprintf(fid,'Vp2 =%8.4f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
end
[mod,arg]=mod_arg(Vs0);fprintf(fid,'Vs0 =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vs1);fprintf(fid,'Vs1 =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vs2);fprintf(fid,'Vs2 =%8.4f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Vr0);fprintf(fid,'Vr0 =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vr1);fprintf(fid,'Vr1 =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vr2);fprintf(fid,'Vr2 =%8.4f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
[mod,arg]=mod_arg(Im0);fprintf(fid,'Im0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Im1);fprintf(fid,'Im1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Im2);fprintf(fid,'Im2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
if TipoFalta==2
[mod,arg]=mod_arg(Ip0);fprintf(fid,'Ip0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ip1);fprintf(fid,'Ip1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ip2);fprintf(fid,'Ip2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
[mod,arg]=mod_arg(Is0);fprintf(fid,'Is0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Is1);fprintf(fid,'Is1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Is2);fprintf(fid,'Is2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Ir0);fprintf(fid,'Ir0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ir1);fprintf(fid,'Ir1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ir2);fprintf(fid,'Ir2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
if (esquema==1)|(esquema==3)
[mod,arg]=mod_arg(Ism0);fprintf(fid,'Ism0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ism1);fprintf(fid,'Ism1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ism2);fprintf(fid,'Ism2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_arg(Isp0);fprintf(fid,'Isp0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isp1);fprintf(fid,'Isp1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isp2);fprintf(fid,'Isp2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_arg(Isp0);fprintf(fid,'Isr0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isp1);fprintf(fid,'Isr1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isp2);fprintf(fid,'Isr2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if (esquema==1)|(esquema==2)
[mod,arg]=mod_arg(Irm0);fprintf(fid,'Irm0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irm1);fprintf(fid,'Irm1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irm2);fprintf(fid,'Irm2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_arg(Irp0);fprintf(fid,'Irp0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irp1);fprintf(fid,'Irp1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irp2);fprintf(fid,'Irp2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_arg(Irp0);fprintf(fid,'Irs0 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irp1);fprintf(fid,'Irs1 =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irp2);fprintf(fid,'Irs2 =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if imptransfer==1
[mod,arg]=mod_arg(Izsr0);fprintf(fid,'Izsr0=%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Izsr1);fprintf(fid,'Izsr1=%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Izsr2);fprintf(fid,'Izsr2=%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
fprintf(fid,'\n');
% ------------------------------------------------------------------
% ** C2 ************ SALIDA V(a,b,c), I(a,b,c) pu Polar ************
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vr
% Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr
% ------------------------------------------------------------------
fprintf(fid,'SALIDA V(a,b,c), I(a,b,c) pu Polar\n');
fprintf(fid,'\n');
[mod,arg]=mod_arg(Vma);fprintf(fid,'Vma =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vmb);fprintf(fid,'Vmb =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vmc);fprintf(fid,'Vmc =%8.4f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
if TipoFalta==2
[mod,arg]=mod_arg(Vpa);fprintf(fid,'Vpa =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vpb);fprintf(fid,'Vpb =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vpc);fprintf(fid,'Vpc =%8.4f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
end
[mod,arg]=mod_arg(Vsa);fprintf(fid,'Vsa =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vsb);fprintf(fid,'Vsb =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vsc);fprintf(fid,'Vsc =%8.4f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Vra);fprintf(fid,'Vra =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vrb);fprintf(fid,'Vrb =%8.4f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vrc);fprintf(fid,'Vrc =%8.4f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
[mod,arg]=mod_arg(Ima);fprintf(fid,'Ima =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Imb);fprintf(fid,'Imb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Imc);fprintf(fid,'Imc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
if TipoFalta==2
[mod,arg]=mod_arg(Ipa);fprintf(fid,'Ipa =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ipb);fprintf(fid,'Ipb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ipc);fprintf(fid,'Ipc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
[mod,arg]=mod_arg(Isa);fprintf(fid,'Isa =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isb);fprintf(fid,'Isb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isc);fprintf(fid,'Isc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Ira);fprintf(fid,'Ira =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irb);fprintf(fid,'Irb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irc);fprintf(fid,'Irc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n\n');
if (esquema==1)|(esquema==3)
[mod,arg]=mod_arg(Isma);fprintf(fid,'Isma =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ismb);fprintf(fid,'Ismb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ismc);fprintf(fid,'Ismc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_arg(Ispa);fprintf(fid,'Ispa =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ispb);fprintf(fid,'Ispb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ispc);fprintf(fid,'Ispc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_arg(Ispa);fprintf(fid,'Isra =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ispb);fprintf(fid,'Isrb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ispc);fprintf(fid,'Isrc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if (esquema==1)|(esquema==2)
[mod,arg]=mod_arg(Irma);fprintf(fid,'Irma =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irmb);fprintf(fid,'Irmb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irmc);fprintf(fid,'Irmc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_arg(Irpa);fprintf(fid,'Irpa =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irpb);fprintf(fid,'Irpb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irpc);fprintf(fid,'Irpc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_arg(Irpa);fprintf(fid,'Irsa =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irpb);fprintf(fid,'Irsb =%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irpc);fprintf(fid,'Irsc =%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if imptransfer==1
[mod,arg]=mod_arg(Izsra);fprintf(fid,'Izsra=%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Izsrb);fprintf(fid,'Izsrb=%8.4f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(Izsrc);fprintf(fid,'Izsrc=%8.4f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
% ***********************************************************************
case{3}, % OVA cartesiana: 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
% --- ----------------- -----------------
% ***********************************************************************
% -----------------------------------------------------------------------
% C3.1 SALIDA 21, Ohm, cartesiana
% -----------------------------------------------------------------------
% Sm(e1,e3),Sp_2,Sr_1,Rm(e1,e2),Rp_2,Rs_1
% -----------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA 21,Ohm, cartesiana \n');
fprintf(fid,'\n');
fprintf(fid,'Nota: Componentes de Z >=10000 Ohm se representan por Inf\n');
fprintf(fid,'Nota: Componentes de Z <=-10000 Ohm se representan por -Inf\n');
fprintf(fid,'\n');
if (esquema==1)|(esquema==3)
fprintf(fid,'Sm =%-7.2f+i(%-7.2f)Ohm\n',real(SSm),imag(SSm));
end
if (esquema==1)|(esquema==2)
fprintf(fid,'Rm =%-7.2f+i(%-7.2f)Ohm\n',real(RRm),imag(RRm));
end
if TipoFalta==2
fprintf('Sp =%-7.2f+i(%-7.2f)Ohm\n',real(SSp),imag(SSp));
fprintf('Rp =%-7.2f+i(%-7.2f)Ohm\n',real(RRp),imag(RRp));
end
fprintf(fid,'\n');
BarraInt=input('Supone una Barra Intermedia A, en Linea 1, un poco antes de m?: Si=1,No=2 : ');
if BarraInt==1
SApu=input('SA/SR = ');
if salida ==2
fprintf(fid,'Suponer Barra Intermedia A, Linea 1: SA/SR =%4.2f\n',SApu);
end
SA=SApu*zl1*zbase;
fprintf(fid,'SA =%-7.2f+i(%-7.2f)Ohm\n',real(SA),imag(SA));
SAzona1=0.8*SA;
fprintf(fid,'SAL1(zona1) =%-7.2f+i(%-7.2f)Ohm\n',real(SAzona1),imag(SAzona1));
fprintf(fid,'\n');
end
if esquema==1|esquema==3
[Sm_ZABre,Sm_ZABim]=zmed(Sm_ZAB,zbase); fprintf(fid,'Sm_ZAB=%7.2f+i(%-7.2f)Ohm\n',Sm_ZABre,Sm_ZABim);
[Sm_ZBCre,Sm_ZBCim]=zmed(Sm_ZBC,zbase); fprintf(fid,'Sm_ZBC=%7.2f+i(%-7.2f)Ohm\n',Sm_ZBCre,Sm_ZBCim);
[Sm_ZCAre,Sm_ZCAim]=zmed(Sm_ZCA,zbase); fprintf(fid,'Sm_ZCA=%7.2f+i(%-7.2f)Ohm\n',Sm_ZCAre,Sm_ZCAim);
[Sm_ZAGre,Sm_ZAGim]=zmed(Sm_ZAG,zbase); fprintf(fid,'Sm_ZAG=%7.2f+i(%-7.2f)Ohm\n',Sm_ZAGre,Sm_ZAGim);
[Sm_ZBGre,Sm_ZBGim]=zmed(Sm_ZBG,zbase); fprintf(fid,'Sm_ZBG=%7.2f+i(%-7.2f)Ohm\n',Sm_ZBGre,Sm_ZBGim);
[Sm_ZCGre,Sm_ZCGim]=zmed(Sm_ZCG,zbase); fprintf(fid,'Sm_ZCG=%7.2f+i(%-7.2f)Ohm\n',Sm_ZCGre,Sm_ZCGim);
fprintf(fid,'\n');
end
if TipoFalta==2
[Sp_ZABre,Sp_ZABim]=zmed(Sp_ZAB,zbase);fprintf(fid,'Sp_ZAB=%7.2f+i(%-7.2f)Ohm\n',Sp_ZABre,Sp_ZABim);
[Sp_ZBCre,Sp_ZBCim]=zmed(Sp_ZBC,zbase);fprintf(fid,'Sp_ZBC=%7.2f+i(%-7.2f)Ohm\n',Sp_ZBCre,Sp_ZBCim);
[Sp_ZCAre,Sp_ZCAim]=zmed(Sp_ZCA,zbase);fprintf(fid,'Sp_ZCA=%7.2f+i(%-7.2f)Ohm\n',Sp_ZCAre,Sp_ZCAim);
[Sp_ZAGre,Sp_ZAGim]=zmed(Sp_ZAG,zbase);fprintf(fid,'Sp_ZAG=%7.2f+i(%-7.2f)Ohm\n',Sp_ZAGre,Sp_ZAGim);
[Sp_ZBGre,Sp_ZBGim]=zmed(Sp_ZBG,zbase);fprintf(fid,'Sp_ZBG=%7.2f+i(%-7.2f)Ohm\n',Sp_ZBGre,Sp_ZBGim);
[Sp_ZCGre,Sp_ZCGim]=zmed(Sp_ZCG,zbase);fprintf(fid,'Sp_ZCG=%7.2f+i(%-7.2f)Ohm\n',Sp_ZCGre,Sp_ZCGim);
fprintf(fid,'\n');
elseif TipoFalta==1
[Sp_ZABre,Sp_ZABim]=zmed(Sp_ZAB,zbase);fprintf(fid,'Sr_ZAB=%7.2f+i(%-7.2f)Ohm\n',Sp_ZABre,Sp_ZABim);
[Sp_ZBCre,Sp_ZBCim]=zmed(Sp_ZBC,zbase);fprintf(fid,'Sr_ZBC=%7.2f+i(%-7.2f)Ohm\n',Sp_ZBCre,Sp_ZBCim);
[Sp_ZCAre,Sp_ZCAim]=zmed(Sp_ZCA,zbase);fprintf(fid,'Sr_ZCA=%7.2f+i(%-7.2f)Ohm\n',Sp_ZCAre,Sp_ZCAim);
[Sp_ZAGre,Sp_ZAGim]=zmed(Sp_ZAG,zbase);fprintf(fid,'Sr_ZAG=%7.2f+i(%-7.2f)Ohm\n',Sp_ZAGre,Sp_ZAGim);
[Sp_ZBGre,Sp_ZBGim]=zmed(Sp_ZBG,zbase);fprintf(fid,'Sr_ZBG=%7.2f+i(%-7.2f)Ohm\n',Sp_ZBGre,Sp_ZBGim);
[Sp_ZCGre,Sp_ZCGim]=zmed(Sp_ZCG,zbase);fprintf(fid,'Sr_ZCG=%7.2f+i(%-7.2f)Ohm\n',Sp_ZCGre,Sp_ZCGim);
fprintf(fid,'\n');
end
if esquema==1|esquema==2
[Rm_ZABre,Rm_ZABim]=zmed(Rm_ZAB,zbase);fprintf(fid,'Rm_ZAB=%7.2f+i(%-7.2f)Ohm\n',Rm_ZABre,Rm_ZABim);
[Rm_ZBCre,Rm_ZBCim]=zmed(Rm_ZBC,zbase);fprintf(fid,'Rm_ZBC=%7.2f+i(%-7.2f)Ohm\n',Rm_ZBCre,Rm_ZBCim);
[Rm_ZCAre,Rm_ZCAim]=zmed(Rm_ZCA,zbase);fprintf(fid,'Rm_ZCA=%7.2f+i(%-7.2f)Ohm\n',Rm_ZCAre,Rm_ZCAim);
[Rm_ZAGre,Rm_ZAGim]=zmed(Rm_ZAG,zbase);fprintf(fid,'Rm_ZAG=%7.2f+i(%-7.2f)Ohm\n',Rm_ZAGre,Rm_ZAGim);
[Rm_ZBGre,Rm_ZBGim]=zmed(Rm_ZBG,zbase);fprintf(fid,'Rm_ZBG=%7.2f+i(%-7.2f)Ohm\n',Rm_ZBGre,Rm_ZBGim);
[Rm_ZCGre,Rm_ZCGim]=zmed(Rm_ZCG,zbase);fprintf(fid,'Rm_ZCG=%7.2f+i(%-7.2f)Ohm\n',Rm_ZCGre,Rm_ZCGim);
fprintf(fid,'\n');
end
if TipoFalta==2
[Rp_ZABre,Rp_ZABim]=zmed(Rp_ZAB,zbase);fprintf(fid,'Rp_ZAB=%7.2f+i(%-7.2f)Ohm\n',Rp_ZABre,Rp_ZABim);
[Rp_ZBCre,Rp_ZBCim]=zmed(Rp_ZBC,zbase);fprintf(fid,'Rp_ZBC=%7.2f+i(%-7.2f)Ohm\n',Rp_ZBCre,Rp_ZBCim);
[Rp_ZCAre,Rp_ZCAim]=zmed(Rp_ZCA,zbase);fprintf(fid,'Rp_ZCA=%7.2f+i(%-7.2f)Ohm\n',Rp_ZCAre,Rp_ZCAim);
[Rp_ZAGre,Rp_ZAGim]=zmed(Rp_ZAG,zbase);fprintf(fid,'Rp_ZAG=%7.2f+i(%-7.2f)Ohm\n',Rp_ZAGre,Rp_ZAGim);
[Rp_ZBGre,Rp_ZBGim]=zmed(Rp_ZBG,zbase);fprintf(fid,'Rp_ZBG=%7.2f+i(%-7.2f)Ohm\n',Rp_ZBGre,Rp_ZBGim);
[Rp_ZCGre,Rp_ZCGim]=zmed(Rp_ZCG,zbase);fprintf(fid,'Rp_ZCG=%7.2f+i(%-7.2f)Ohm\n',Rp_ZCGre,Rp_ZCGim);
fprintf(fid,'\n\n');
elseif TipoFalta==1
[Rp_ZABre,Rp_ZABim]=zmed(Rp_ZAB,zbase);fprintf(fid,'Rs_ZAB=%7.2f+i(%-7.2f)Ohm\n',Rp_ZABre,Rp_ZABim);
[Rp_ZBCre,Rp_ZBCim]=zmed(Rp_ZBC,zbase);fprintf(fid,'Rs_ZBC=%7.2f+i(%-7.2f)Ohm\n',Rp_ZBCre,Rp_ZBCim);
[Rp_ZCAre,Rp_ZCAim]=zmed(Rp_ZCA,zbase);fprintf(fid,'Rs_ZCA=%7.2f+i(%-7.2f)Ohm\n',Rp_ZCAre,Rp_ZCAim);
[Rp_ZAGre,Rp_ZAGim]=zmed(Rp_ZAG,zbase);fprintf(fid,'Rs_ZAG=%7.2f+i(%-7.2f)Ohm\n',Rp_ZAGre,Rp_ZAGim);
[Rp_ZBGre,Rp_ZBGim]=zmed(Rp_ZBG,zbase);fprintf(fid,'Rs_ZBG=%7.2f+i(%-7.2f)Ohm\n',Rp_ZBGre,Rp_ZBGim);
[Rp_ZCGre,Rp_ZCGim]=zmed(Rp_ZCG,zbase);fprintf(fid,'Rs_ZCG=%7.2f+i(%-7.2f)Ohm\n',Rp_ZCGre,Rp_ZCGim);
fprintf(fid,'\n\n');
end
% ------------------------------------------------------------------
% C3.3 SALIDA V(0,1,2),I(0,1,2) Volt-Amp cartesiana
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vr
% Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr
% ------------------------------------------------------------------
fprintf(fid,'SALIDA V(0,1,2),I(0,1,2) Volt, Amp, cartesiana \n');
fprintf(fid,'\n');
fprintf(fid,'Vm0 =%8.0f+i(%-8.0f)V\t',real(VVm0),imag(VVm0));
fprintf(fid,'Vm1 =%8.0f+i(%-8.0f)V\t',real(VVm1),imag(VVm1));
fprintf(fid,'Vm2 =%8.0f+i(%-8.0f)V\t',real(VVm2),imag(VVm2));
if TipoFalta==2
fprintf(fid,'\n');
fprintf(fid,'Vp0 =%8.0f+i(%-8.0f)V\t',real(VVp0),imag(VVp0));
fprintf(fid,'Vp1 =%8.0f+i(%-8.0f)V\t',real(VVp1),imag(VVp1));
fprintf(fid,'Vp2 =%8.0f+i(%-8.0f)V\t',real(VVp2),imag(VVp2));
end
fprintf(fid,'\n');
fprintf(fid,'Vs0 =%8.0f+i(%-8.0f)V\t',real(VVs0),imag(VVs0));
fprintf(fid,'Vs1 =%8.0f+i(%-8.0f)V\t',real(VVs1),imag(VVs1));
fprintf(fid,'Vs2 =%8.0f+i(%-8.0f)V\t',real(VVs2),imag(VVs2));
fprintf(fid,'\n');
fprintf(fid,'Vr0 =%8.0f+i(%-8.0f)V\t',real(VVr0),imag(VVr0));
fprintf(fid,'Vr1 =%8.0f+i(%-8.0f)V\t',real(VVr1),imag(VVr1));
fprintf(fid,'Vr2 =%8.0f+i(%-8.0f)V\t',real(VVr2),imag(VVr2));
fprintf(fid,'\n\n');
% ---------------------------------------------------------------
fprintf(fid,'Im0 =%8.0f+i(%-8.0f)A\t',real(IIm0),imag(IIm0));
fprintf(fid,'Im1 =%8.0f+i(%-8.0f)A\t',real(IIm1),imag(IIm1));
fprintf(fid,'Im2 =%8.0f+i(%-8.0f)A\t',real(IIm2),imag(IIm2));
if TipoFalta==2
fprintf(fid,'\n');
fprintf(fid,'Ip0 =%8.0f+i(%-8.0f)A\t',real(IIp0),imag(IIp0));
fprintf(fid,'Ip1 =%8.0f+i(%-8.0f)A\t',real(IIp1),imag(IIp1));
fprintf(fid,'Ip2 =%8.0f+i(%-8.0f)A\t',real(IIp2),imag(IIp2));
end
fprintf(fid,'\n');
fprintf(fid,'Is0 =%8.0f+i(%-8.0f)A\t',real(IIs0),imag(IIs0));
fprintf(fid,'Is1 =%8.0f+i(%-8.0f)A\t',real(IIs1),imag(IIs1));
fprintf(fid,'Is2 =%8.0f+i(%-8.0f)A\t',real(IIs2),imag(IIs2));
fprintf(fid,'\n');
fprintf(fid,'Ir0 =%8.0f+i(%-8.0f)A\t',real(IIr0),imag(IIr0));
fprintf(fid,'Ir1 =%8.0f+i(%-8.0f)A\t',real(IIr1),imag(IIr1));
fprintf(fid,'Ir2 =%8.0f+i(%-8.0f)A\t',real(IIr2),imag(IIr2));
fprintf(fid,'\n');
if esquema==1|esquema==3
fprintf(fid,'Ism0 =%8.0f+i(%-8.0f)A\t',real(IIsm0),imag(IIsm0));
fprintf(fid,'Ism1 =%8.0f+i(%-8.0f)A\t',real(IIsm1),imag(IIsm1));
fprintf(fid,'Ism2 =%8.0f+i(%-8.0f)A\t',real(IIsm2),imag(IIsm2));
fprintf(fid,'\n');
end
if TipoFalta==2
fprintf(fid,'Isp0 =%8.0f+i(%-8.0f)A\t',real(IIsp0),imag(IIsp0));
fprintf(fid,'Isp1 =%8.0f+i(%-8.0f)A\t',real(IIsp1),imag(IIsp1));
fprintf(fid,'Isp2 =%8.0f+i(%-8.0f)A\t',real(IIsp2),imag(IIsp2));
fprintf(fid,'\n');
elseif TipoFalta==1
fprintf(fid,'Isr0 =%8.0f+i(%-8.0f)A\t',real(IIsp0),imag(IIsp0));
fprintf(fid,'Isr1 =%8.0f+i(%-8.0f)A\t',real(IIsp1),imag(IIsp1));
fprintf(fid,'Isr2 =%8.0f+i(%-8.0f)A\t',real(IIsp2),imag(IIsp2));
fprintf(fid,'\n');
end
if esquema==1|esquema==2
fprintf(fid,'Irm0 =%8.0f+i(%-8.0f)A\t',real(IIrm0),imag(IIrm0));
fprintf(fid,'Irm1 =%8.0f+i(%-8.0f)A\t',real(IIrm1),imag(IIrm1));
fprintf(fid,'Irm2 =%8.0f+i(%-8.0f)A\t',real(IIrm2),imag(IIrm2));
fprintf(fid,'\n');
end
if TipoFalta==2
fprintf(fid,'Irp0 =:%8.0f+i(%-8.0f)A\t',real(IIrp0),imag(IIrp0));fprintf(fid,'Irp1 :%8.0f+i(%-8.0f)A\t',real(IIrp1),imag(IIrp1));
fprintf(fid,'Irp2 =%8.0f+i(%-8.0f)A\t',real(IIrp2),imag(IIrp2));
fprintf(fid,'\n\n');
elseif TipoFalta==1
fprintf(fid,'Irs0 =%8.0f+i(%-8.0f)A\t',real(IIrp0),imag(IIrp0));
fprintf(fid,'Irs1 =%8.0f+i(%-8.0f)A\t',real(IIrp1),imag(IIrp1));
fprintf(fid,'Irs2 =%8.0f+i(%-8.0f)A\t',real(IIrp2),imag(IIrp2));
fprintf(fid,'\n');
end
if imptransfer==1
fprintf(fid,'Izsr0=%8.0f+i(%-8.0f)A\t',real(IIzsr0),imag(IIzsr0));
fprintf(fid,'Izsr1=%8.0f+i(%-8.0f)A\t',real(IIzsr1),imag(IIzsr1));
fprintf(fid,'Izsr2=%8.0f+i(%-8.0f)A\t',real(IIzsr2),imag(IIzsr2));
fprintf(fid,'\n');
end
% ------------------------------------------------------------------
% C3.4 SALIDA V(a,b,c),I(a,b,c) Volt,Amp cartesiana
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vr
% Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA V(a,b,c),I(a,b,c) Volt,Amp cartesiana\n');
fprintf(fid,'\n');
fprintf(fid,'Vma =%8.0f+i(%-8.0f)V\t',real(VVma),imag(VVma));
fprintf(fid,'Vmb =%8.0f+i(%-8.0f)V\t',real(VVmb),imag(VVmb));
fprintf(fid,'Vmc =%8.0f+i(%-8.0f)V\t',real(VVmc),imag(VVmc));
if TipoFalta==2
fprintf(fid,'\n');
fprintf(fid,'Vpa =%8.0f+i(%-8.0f)V\t',real(VVpa),imag(VVpa));
fprintf(fid,'Vpb =%8.0f+i(%-8.0f)V\t',real(VVpb),imag(VVpb));
fprintf(fid,'Vpc =%8.0f+i(%-8.0f)V\t',real(VVpc),imag(VVpc));
end
fprintf(fid,'\n');
fprintf(fid,'Vsa =%8.0f+i(%-8.0f)V\t',real(VVsa),imag(VVsa));
fprintf(fid,'Vsb =%8.0f+i(%-8.0f)V\t',real(VVsb),imag(VVsb));
fprintf(fid,'Vsc =%8.0f+i(%-8.0f)V\t',real(VVsc),imag(VVsc));
fprintf(fid,'\n');
fprintf(fid,'Vra =%8.0f+i(%-8.0f)V\t',real(VVra),imag(VVra));
fprintf(fid,'Vrb =%8.0f+i(%-8.0f)V\t',real(VVrb),imag(VVrb));
fprintf(fid,'Vrc =%8.0f+i(%-8.0f)V\t',real(VVrc),imag(VVrc));
fprintf(fid,'\n\n');
fprintf(fid,'Ima =%8.0f+i(%-8.0f)A\t',real(IIma),imag(IIma));
fprintf(fid,'Imb =%8.0f+i(%-8.0f)A\t',real(IImb),imag(IImb));
fprintf(fid,'Imc =%8.0f+i(%-8.0f)A\t',real(IImc),imag(IImc));
if TipoFalta==2
fprintf(fid,'\n');
fprintf(fid,'Ipa =%8.0f+i(%-8.0f)A\t',real(IIpa),imag(IIpa));
fprintf(fid,'Ipb =%8.0f+i(%-8.0f)A\t',real(IIpb),imag(IIpb));
fprintf(fid,'Ipc =%8.0f+i(%-8.0f)A\t',real(IIpc),imag(IIpc));
end
fprintf(fid,'\n');
fprintf(fid,'Isa =%8.0f+i(%-8.0f)A\t',real(IIsa),imag(IIsa));
fprintf(fid,'Isb =%8.0f+i(%-8.0f)A\t',real(IIsb),imag(IIsb));
fprintf(fid,'Isc =%8.0f+i(%-8.0f)A\t',real(IIsc),imag(IIsc));
fprintf(fid,'\n');
fprintf(fid,'Ira =%8.0f+i(%-8.0f)A\t',real(IIra),imag(IIra));
fprintf(fid,'Irb =%8.0f+i(%-8.0f)A\t',real(IIrb),imag(IIrb));
fprintf(fid,'Irc =%8.0f+i(%-8.0f)A\t',real(IIrc),imag(IIrc));
fprintf(fid,'\n');
if esquema==1|esquema==3
fprintf(fid,'Isma =%8.0f+i(%-8.0f)A\t',real(IIsma),imag(IIsma));
fprintf(fid,'Ismb =%8.0f+i(%-8.0f)A\t',real(IIsmb),imag(IIsmb));
fprintf(fid,'Ismc =%8.0f+i(%-8.0f)A\t',real(IIsmc),imag(IIsmc));
fprintf(fid,'\n');
end
if TipoFalta==2
fprintf(fid,'Ispa =%8.0f+i(%-8.0f)A\t',real(IIspa),imag(IIspa));
fprintf(fid,'Ispb =%8.0f+i(%-8.0f)A\t',real(IIspb),imag(IIspb));
fprintf(fid,'Ispc =%8.0f+i(%-8.0f)A\t',real(IIspc),imag(IIspc));
fprintf(fid,'\n');
elseif TipoFalta==1
fprintf(fid,'Isra =%8.0f+i(%-8.0f)A\t',real(IIspa),imag(IIspa));
fprintf(fid,'Isrb =%8.0f+i(%-8.0f)A\t',real(IIspb),imag(IIspb));
fprintf(fid,'Isrc =%8.0f+i(%-8.0f)A\t',real(IIspc),imag(IIspc));
fprintf(fid,'\n');
end
if esquema==1|esquema==2
fprintf(fid,'Irma =%8.0f+i(%-8.0f)A\t',real(IIrma),imag(IIrma));
fprintf(fid,'Irmb =%8.0f+i(%-8.0f)A\t',real(IIrmb),imag(IIrmb));
fprintf(fid,'Irmc =%8.0f+i(%-8.0f)A\t',real(IIrmc),imag(IIrmc));
fprintf(fid,'\n');
end
if TipoFalta==2
fprintf(fid,'Irpa =%8.0f+i(%-8.0f)A\t',real(IIrpa),imag(IIrpa));
fprintf(fid,'Irpb =%8.0f+i(%-8.0f)A\t',real(IIrpb),imag(IIrpb));
fprintf(fid,'Irpc =%8.0f+i(%-8.0f)A\t',real(IIrpc),imag(IIrpc));
fprintf(fid,'\n');
elseif TipoFalta==1
fprintf(fid,'Irsa =%8.0f+i(%-8.0f)A\t',real(IIrpa),imag(IIrpa));
fprintf(fid,'Irsb =%8.0f+i(%-8.0f)A\t',real(IIrpb),imag(IIrpb));
fprintf(fid,'Irsc =%8.0f+i(%-8.0f)A\t',real(IIrpc),imag(IIrpc));
fprintf(fid,'\n');
end
if imptransfer==1
fprintf(fid,'Izsra=%8.0f+i(%-8.0f)A\t',real(IIzsra),imag(IIzsra));
fprintf(fid,'Izsrb=%8.0f+i(%-8.0f)A\t',real(IIzsrb),imag(IIzsrb));
fprintf(fid,'Izsrc=%8.0f+i(%-8.0f)A\t',real(IIzsrc),imag(IIzsrc));
fprintf(fid,'\n');
end
% ******************************************************************
case{4}, % OVA polar: 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
% --- ---
% ******************************************************************
% ------------------------------------------------------------------
% ** C4 ****************** SALIDA 21, Ohm , polar ******************
% ------------------------------------------------------------------
% Sm(e1,e3),Sp_2,Sr_1,Rm(e1,e2),Rp_2,Rs_1
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA 21, Ohm , polar \n');
fprintf(fid,'\n');
fprintf(fid,'Nota: Componentes de Z >=10000 Ohm se representan por Inf\n');
fprintf(fid,'Nota: Componentes de Z <=-10000 Ohm se representan por -Inf\n');
fprintf(fid,'\n');
if (esquema==1)|(esquema==3)
[mod,arg]=mod_arg(SSm);fprintf(fid,'Sm=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
end
if (esquema==1)|(esquema==2)
[mod,arg]=mod_arg(RRm);fprintf(fid,'Rm=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
end
if TipoFalta==2
[mod,arg]=mod_arg(SSp);fprintf(fid,'Sp=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_arg(RRp);fprintf(fid,'Rp=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
end
fprintf(fid,'\n');
if (esquema==1)|(esquema==3)
[mod,arg]=mod_argz(Sm_ZAB*zbase);fprintf(fid,'Sm_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sm_ZBC*zbase);fprintf(fid,'Sm_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sm_ZCA*zbase);fprintf(fid,'Sm_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sm_ZAG*zbase);fprintf(fid,'Sm_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sm_ZBG*zbase);fprintf(fid,'Sm_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sm_ZCG*zbase);fprintf(fid,'Sm_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_argz(Sp_ZAB*zbase);fprintf(fid,'Sp_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZBC*zbase);fprintf(fid,'Sp_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZCA*zbase);fprintf(fid,'Sp_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZAG*zbase);fprintf(fid,'Sp_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZBG*zbase);fprintf(fid,'Sp_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZCG*zbase);fprintf(fid,'Sp_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_argz(Sp_ZAB*zbase);fprintf(fid,'Sr_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZBC*zbase);fprintf(fid,'Sr_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZCA*zbase);fprintf(fid,'Sr_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZAG*zbase);fprintf(fid,'Sr_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZBG*zbase);fprintf(fid,'Sr_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sp_ZCG*zbase);fprintf(fid,'Sr_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
end
if (esquema==1)|(esquema==2)
[mod,arg]=mod_argz(Rm_ZAB*zbase);fprintf(fid,'Rm_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rm_ZBC*zbase);fprintf(fid,'Rm_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rm_ZCA*zbase);fprintf(fid,'Rm_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rm_ZAG*zbase);fprintf(fid,'Rm_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rm_ZBG*zbase);fprintf(fid,'Rm_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rm_ZCG*zbase);fprintf(fid,'Rm_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_argz(Rp_ZAB*zbase);fprintf(fid,'Rp_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZBC*zbase);fprintf(fid,'Rp_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZCA*zbase);fprintf(fid,'Rp_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZAG*zbase);fprintf(fid,'Rp_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZBG*zbase);fprintf(fid,'Rp_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZCG*zbase);fprintf(fid,'Rp_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_argz(Rp_ZAB*zbase);fprintf(fid,'Rs_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZBC*zbase);fprintf(fid,'Rs_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZCA*zbase);fprintf(fid,'Rs_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZAG*zbase);fprintf(fid,'Rs_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZBG*zbase);fprintf(fid,'Rs_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rp_ZCG*zbase);fprintf(fid,'Rs_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
end
% -----------------------------------------------------------------------
% ** C4 ****************** SALIDA 67N, Volt, Amp, Polar *****************
% -----------------------------------------------------------------------
% Sm(e1,e3),Sp_2,Sr_1,Rm(e1,e2),Rp_2,Rs_1
% -----------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA 67N, Volt, Amp, polar\n');
fprintf(fid,'\n');
if (esquema==1)|(esquema==3)
fprintf(fid,'Sm_67N:\tVspol=%-8.0f(V)\tIsmop=%-8.0f(A)\tAngDirSm=%-5.1f\n',VVspol,IIsmop,AngDirSm);
end
if TipoFalta==2
fprintf(fid,'Sp_67N:\tVspol=%-8.0f(V)\tIspop=%-8.0f(A)\tAngDirSp=%-5.1f\n',VVspol,IIspop,AngDirSp);
elseif TipoFalta==1
fprintf(fid,'Sr_67N:\tVspol=%-8.0f(V)\tIsrop=%-8.0f(A)\tAngDirSr=%-5.1f\n',VVspol,IIspop,AngDirSp);
end
if (esquema==1)|(esquema==2)
fprintf(fid,'Rm_67N:\tVrpol=%-8.0f(V)\tIrmop=%-8.0f(A)\tAngDirRm=%-5.1f\n',VVrpol,IIrmop,AngDirRm);
end
if TipoFalta==2
fprintf(fid,'Rp_67N:\tVrpol=%-8.0f(V)\tIrpop=%-8.0f(A)\tAngDirRp=%-5.1f\n',VVrpol,IIrpop,AngDirRp);
elseif TipoFalta==1
fprintf(fid,'Rs_67N:\tVrpol=%-8.0f(V)\tIrsop=%-8.0f(A)\tAngDirRs=%-5.1f\n',VVrpol,IIrpop,AngDirRp);
end
fprintf(fid,'\n');
% ------------------------------------------------------------------
% ** C4 ********* SALIDA V(0,1,2), I(0,1,2) Volt,Amp polar ********
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vr
% Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA V(0,1,2), I(0,1,2) Volt-Amp polar\n');
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVm0);fprintf(fid,'Vm0 =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVm1);fprintf(fid,'Vm1 =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVm2);fprintf(fid,'Vm2 =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
if TipoFalta==2
[mod,arg]=mod_arg(VVp0);fprintf(fid,'Vp0 =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVp1);fprintf(fid,'Vp1 =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVp2);fprintf(fid,'Vp2 =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
end
[mod,arg]=mod_arg(VVs0);fprintf(fid,'Vs0 =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVs1);fprintf(fid,'Vs1 =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVs2);fprintf(fid,'Vs2 =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVr0);fprintf(fid,'Vr0 =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVr1);fprintf(fid,'Vr1 =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVr2);fprintf(fid,'Vr2 =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
[mod,arg]=mod_arg(IIm0);fprintf(fid,'Im0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIm1);fprintf(fid,'Im1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIm2);fprintf(fid,'Im2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
if TipoFalta==2
[mod,arg]=mod_arg(IIp0);fprintf(fid,'Ip0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIp1);fprintf(fid,'Ip1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIp2);fprintf(fid,'Ip2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
[mod,arg]=mod_arg(IIs0);fprintf(fid,'Is0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIs1);fprintf(fid,'Is1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIs2);fprintf(fid,'Is2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIr0);fprintf(fid,'Ir0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIr1);fprintf(fid,'Ir1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIr2);fprintf(fid,'Ir2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
if (esquema==1)|(esquema==3)
[mod,arg]=mod_arg(IIsm0);fprintf(fid,'Ism0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsm1);fprintf(fid,'Ism1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsm2);fprintf(fid,'Ism2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_arg(IIsp0);fprintf(fid,'Isp0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsp1);fprintf(fid,'Isp1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsp2);fprintf(fid,'Isp2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_arg(IIsp0);fprintf(fid,'Isr0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsp1);fprintf(fid,'Isr1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsp2);fprintf(fid,'Isr2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if (esquema==1)|(esquema==2)
[mod,arg]=mod_arg(IIrm0);fprintf(fid,'Irm0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrm1);fprintf(fid,'Irm1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrm2);fprintf(fid,'Irm2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_arg(IIrp0);fprintf(fid,'Irp0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrp1);fprintf(fid,'Irp1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrp2);fprintf(fid,'Irp2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_arg(IIrp0);fprintf(fid,'Irs0 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrp1);fprintf(fid,'Irs1 =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrp2);fprintf(fid,'Irs2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if imptransfer==1
[mod,arg]=mod_arg(IIzsr0);fprintf(fid,'Izsr0=%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIzsr1);fprintf(fid,'Izsr1=%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIzsr2);fprintf(fid,'Izsr2=%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
% ------------------------------------------------------------------
% ** C4 ******** SALIDA V(a,b,c), I(a,b,c) Volt,Amp polar *********
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vr
% Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA V(a,b,c), I(a,b,c) Volt,Amp Polar\n');
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVma);fprintf(fid,'Vma =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVmb);fprintf(fid,'Vmb =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVmc);fprintf(fid,'Vmc =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
if TipoFalta==2
[mod,arg]=mod_arg(VVpa);fprintf(fid,'Vpa =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVpb);fprintf(fid,'Vpb =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVpc);fprintf(fid,'Vpc =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
end
[mod,arg]=mod_arg(VVsa);fprintf(fid,'Vsa =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVsb);fprintf(fid,'Vsb =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVsc);fprintf(fid,'Vsc =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVra);fprintf(fid,'Vra =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVrb);fprintf(fid,'Vrb =%8.0f(%-5.1f)(V)\t',mod,arg);
[mod,arg]=mod_arg(VVrc);fprintf(fid,'Vrc =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
[mod,arg]=mod_arg(IIma);fprintf(fid,'Ima =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IImb);fprintf(fid,'Imb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IImc);fprintf(fid,'Imc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
if TipoFalta==2
[mod,arg]=mod_arg(IIpa);fprintf(fid,'Ipa =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIpb);fprintf(fid,'Ipb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIpc);fprintf(fid,'Ipc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
[mod,arg]=mod_arg(IIsa);fprintf(fid,'Isa =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsb);fprintf(fid,'Isb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsc);fprintf(fid,'Isc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIra);fprintf(fid,'Ira =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrb);fprintf(fid,'Irb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrc);fprintf(fid,'Irc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n\n');
if (esquema==1)|(esquema==3)
[mod,arg]=mod_arg(IIsma);fprintf(fid,'Isma =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsmb);fprintf(fid,'Ismb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIsmc);fprintf(fid,'Ismc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_arg(IIspa);fprintf(fid,'Ispa =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIspb);fprintf(fid,'Ispb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIspc);fprintf(fid,'Ispc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_arg(IIspa);fprintf(fid,'Isra =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIspb);fprintf(fid,'Isrb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIspc);fprintf(fid,'Isrc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if (esquema==1)|(esquema==2)
[mod,arg]=mod_arg(IIrma);fprintf(fid,'Irma =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrmb);fprintf(fid,'Irmb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrmc);fprintf(fid,'Irmc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==2
[mod,arg]=mod_arg(IIrpa);fprintf(fid,'Irpa =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrpb);fprintf(fid,'Irpb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrpc);fprintf(fid,'Irpc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if TipoFalta==1
[mod,arg]=mod_arg(IIrpa);fprintf(fid,'Irsa =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrpb);fprintf(fid,'Irsb =%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIrpc);fprintf(fid,'Irsc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
if imptransfer==1
[mod,arg]=mod_arg(IIzsra);fprintf(fid,'Izsra=%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIzsrb);fprintf(fid,'Izsrb=%8.0f(%-5.1f)(A)\t',mod,arg);
[mod,arg]=mod_arg(IIzsrc);fprintf(fid,'Izsrc=%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
end
end % end del switch(output)
fprintf(fid,'\n');
variassalidas=input('Mas Salidas=1,Acabar=2 :');
fprintf(fid,'\n');
end % end del while variassalidas
if salida==2
fclose(fid);
end
fprintf(1,'\n');
variasejecuciones=input('Mas Ejecuciones=1,Acabar=2 :');
end % end del while variasejecuciones
fprintf(1,'Pulse Enter\n');
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
