clc
variasejecuciones=1;
while variasejecuciones==1
clear
% ---------------------------------------------------
fprintf('PC54\n');
fprintf('FC54.m\n\n');
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,'\n');
fprintf(fid,'PC54\n');
if salida==2
fprintf(fid,'%s\n',cadena1);
end
fprintf(fid,'Esquema Linea Doble SR, Bifurcando a X,Y\n');
fprintf(fid,'5 Tramos\n');
fprintf(fid,'Unominal kV: 400,220,132,66\n');
fprintf(fid,'Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf(fid,'Eleccion: kE,L1,L2,L3,L4,L5,Sccs,Sccx,Sccr,Sccy\n');
fprintf(fid,' Kx0,Kr0,Ky0,k0=zs0/zr0\n');
fprintf(fid,'Compensacion 21: Si,No\n');
fprintf(fid,'Falta Simple m, o Doble m,p\n');
fprintf(fid,'Barra m,p a elegir entre a,b,c,d,e,f,g,h,q,r\n');
fprintf(fid,'Para cada Tramo: Barra m=1, Barra p=2\n');
% ---------------------------------------------------
%PC54
datos=input('datos: standard=1,input=2: ');
if datos==1
Un=220;
kE=1;
L1=100;
L2=100;
L3=100;
L4=100;
L5=100;
compens21=2;
Sccs=5000;
Sccx=5000;
Sccr=5000;
Sccy=5000;
Kx0=1;
Kr0=1;
Ky0=1;
% ---------------------------------------------------
else
Un=input('Unominal kV:400,220,132,66: ');
kE=input('kE=E/Enominal: ');
L1=input('L1(km)= ');
L2=input('L2(km)= ');
L3=input('L3(km)= ');
L4=input('L4(km)= ');
L5=input('L5(km)= ');
compens21=input('Compensacion 21:Si=1,No=2: ');
Sccs=input('Sccs(MVA)= ');
Sccx=input('Sccx(MVA)= ');
Sccr=input('Sccr(MVA)= ');
Sccy=input('Sccy(MVA)= ');
Kx0=input('Kx0= ');
Kr0=input('Kr0= ');
Ky0=input('Ky0= ');
end
% ---------------------------------------------------
k0=input('k0=zs0/zr0= ');
TipoFalta=input('Falta Simple m=1,Doble m,p=2: ');
fprintf(fid,'\n');
% ---------------------------------------------------
% Esquema
% 4 Barras (S,X,R,Y), Linea Doble SR, Falta Multiple (m,p)
% Tramo 1: S-Entronque 1 = L1 Km
% Tramo 2:Entronque 1-Entronque 2 = L2 Km
% Tramo 3:Entronque 2-R = L3 Km
% Tramo 4:Bifurcacion de SR1 en Linea Doble Entronque 1 - X = L4 Km
% Tramo 5:Bifurcacion de SR2 en Linea Doble Entronque 2 - Y = L5 Km
% Fuentes S,R,X,Y.
% Red previa a la Falta en Vacio
% ------------------------------------------------------------------
% DATOS
% ------------------------------------------------------------------
Ubase=Un;
Vbase=Ubase/(sqrt(3));
MVAbase=100;
Ibase=MVAbase/(3*Vbase);
zbase=Vbase/Ibase;
if Un==400
zl1u=(0.0346+0.3208i)/zbase; zl2u=zl1u; zl0u=(0.1828+1.1482i)/zbase;
zlM0u=(0.1482+0.7383i)/zbase;
elseif Un==220
zl1u=(0.0829+0.4135i)/zbase; zl2u=zl1u; zl0u=(0.2311+1.3085i)/zbase;
zlM0u=(0.1482+0.8031i)/zbase;
elseif Un==132
zl1u=(0.1342+0.4042i)/zbase; zl2u=zl1u; zl0u=(0.2824+1.3718i)/zbase;
zlM0u=(0.1482+0.8947i)/zbase;
elseif Un==66
zl1u=(0.2156+0.3952i)/zbase; zl2u=zl1u; zl0u=(0.3636+1.4278i)/zbase;
zlM0u=(0.1482+0.9288i)/zbase;
end
% ----------------------------------------------------------------------
a=0.5;b=0.5;c=0.5;d=0.5;e=0.5;f=0.5;g=0.5;h=0.5;q=0.5;r=0.5;
% ----------------------------------------------------------------------
%CALCULOS PRELIMINARES
% ----------------------------------------------------------------------
U=kE*Un; % kV
E=U/sqrt(3)+0i; E=E/Vbase; % E pu
k=(zl0u-zl1u)/(3*zl1u);
Km=zlM0u/zl1u;
zs1=(100/Sccs)*i; zx1=(100/Sccx)*i; zr1=(100/Sccr)*i; zy1=(100/Sccy)*i; % pu
zs2=zs1; zx2=zx1; zr2=zr1; zy2=zy1; % pu
zx0=Kx0*zx1; zr0=Kr0*zr1; zy0=Ky0*zy1; % pu
zs0=k0*zr0;
% ----------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'Un =%-4.0f(kV)\n',Un);
fprintf(fid,'kE=E/Enominal =%-4.0f\n',kE);
fprintf(fid,'Tipo de Falta, Simple(1), Doble(2) =%-1.0f\n',TipoFalta);
fprintf(fid,'Compensacion 21, Si(1), No(2) =%-1.0f\n',compens21);
fprintf(fid,'L1=%-5.1f(km)\tL2=%-5.1f(km)\tL3=%-5.1f(km)',L1,L2,L3);
fprintf(fid,'\tL4=%-5.1f(km)\tL5=%-5.1f(km)\n',L4,L5);
fprintf(fid,'Sccs=%-5.0f(MVA)\tSccx=%-5.0f(MVA)\tSccr=%-5.0f(MVA)\tSccy=%-5.0f(MVA)\n',Sccs,Sccx,Sccr,Sccy);
fprintf(fid,'Kx0=%-6.2f\tKr0=%-6.2f\tKy0=%-6.2f\tk0=%-6.2f\n',Kx0,Kr0,Ky0,k0);
fprintf(fid,'\n');
fprintf(fid,'Tramo con Falta m: T1sup=11,T4izq=12,T4der=13,T2sup=14,T3sup=15,T3inf=16\n');
fprintf(fid,'T5der=17,T5izq=18,T2inf=19,T1inf=20:');
if salida==1
pfaltam=input('Tramo con Falta m: ');
elseif salida==2
fprintf('Tramo con Falta m: T1sup=11,T4izq=12,T4der=13,T2sup=14,T3sup=15,T3inf=16\n');
fprintf('T5der=17,T5izq=18,T2inf=19,T1inf=20 \n');
pfaltam=input('Tramo con Falta m: ');
fprintf(fid,'%3.0f\n',pfaltam);
end
if pfaltam==11
fprintf(fid,'Falta m en Linea SX, entre S y E1\n');
a=input('m =(0...1)=');
m=a;
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Sm=zl1u*m*L1; Xm=zl1u*(L4+(1-m)*L1); % pu
Sm=Sm*zbase; Xm=Xm*zbase; % Ohm
fprintf(fid,'Sm =%-5.2f+i(%-5.2f)Ohm\n',real(Sm),imag(Sm));
fprintf(fid,'Xm =%-5.2f+i(%-5.2f)Ohm\n',real(Xm),imag(Xm));
fprintf(fid,'\n');
elseif pfaltam==12
fprintf(fid,'Falta m en Linea SX, entre E1 y X\n');
b=input('m=(0...1)=');
m=b;
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Sm=zl1u*(L1+m*L4); Xm=zl1u*(1-m)*L4;
Sm=Sm*zbase; Xm=Xm*zbase;
fprintf(fid,'Sm =%-5.2f+i(%-5.2f)Ohm\n',real(Sm),imag(Sm));
fprintf(fid,'Xm =%-5.2f+i(%-5.2f)Ohm\n',real(Xm),imag(Xm));
fprintf(fid,'\n');
elseif pfaltam==13
fprintf(fid,'Falta m en Linea RX, entre E1 y X\n');
c=input('m=(0...1)=');
m=c;
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Xm=zl1u*(1-m)*L4; Rm=zl1u*(L3+L2+m*L4);
Xm=Xm*zbase; Rm=Rm*zbase;
fprintf(fid,'Xm =%-5.2f+i(%-5.2f)Ohm',real(Xm),imag(Xm));
fprintf(fid,'Rm =%-5.2f+i(%-5.2f)Ohm\n',real(Rm),imag(Rm));
fprintf(fid,'\n');
elseif pfaltam==14
fprintf(fid,'Falta m en Linea RX, entre E1 y E2\n');
d=input('m=(0...1)=');
m=d;
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Xm=zl1u*(L4+m*L2); Rm=zl1u*(L3+(1-m)*L2);
Xm=Xm*zbase; Rm=Rm*zbase;
fprintf(fid,'Xm =%-5.2f+i(%-5.2f)Ohm\n',real(Xm),imag(Xm));
fprintf(fid,'Rm =%-5.2f+i(%-5.2f)Ohm\n',real(Rm),imag(Rm));
fprintf('\n');
elseif pfaltam==15
fprintf(fid,'Falta m en Linea RX, entre E2 y R\n');
e=input('m=(0...1)=');
m=e;
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Xm=zl1u*(L4+L2+m*L3); Rm=zl1u*(1-m)*L3;
Xm=Xm*zbase; Rm=Rm*zbase;
fprintf(fid,'Xm =%-5.2f+i(%-5.2f)Ohm\n',real(Xm),imag(Xm));
fprintf(fid,'Rm =%-5.2f+i(%-5.2f)Ohm\n',real(Rm),imag(Rm));
fprintf(fid,'\n');
elseif pfaltam==16
fprintf(fid,'Falta m en Linea RY, entre E2 y Y\n');
f=input('m=(0...1)=');
m=f; if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Rm=zl1u*(1-m)*L3; Ym=zl1u*(L5+m*L3);
Rm=Rm*zbase; Ym=Ym*zbase;
fprintf(fid,'Rm =%-5.2f+i(%-5.2f)Ohm',real(Rm),imag(Rm));
fprintf(fid,'Ym =%-5.2f+i(%-5.2f)Ohm\n',real(Ym),imag(Ym));
fprintf(fid,'\n');
elseif pfaltam==17
fprintf(fid,'Falta m en Linea RY, entre E2 y Y\n');
g=input('m=(0...1)=');
m=g;
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Rm=zl1u*(L3+m*L5); Ym=zl1u*(1-m)*L5;
Rm=Rm*zbase; Ym=Ym*zbase;
fprintf(fid,'Rm =%-5.2f+i(%-5.2f)Ohm\n',real(Rm),imag(Rm));
fprintf(fid,'Ym =%-5.2f+i(%-5.2f)Ohm\n',real(Ym),imag(Ym));
fprintf(fid,'\n');
elseif pfaltam==18
fprintf(fid,'Falta m en Linea SY, entre E2 y Y\n');
h=input('m=(0...1)=');
m=h;
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Sm=zl1u*(L1+L2+m*L5); Ym=zl1u*(1-m)*L5;
Sm=Sm*zbase; Ym=Ym*zbase;
fprintf(fid,'Sm =%-5.2f+i(%-5.2f)Ohm\n',real(Sm),imag(Sm));
fprintf(fid,'Ym =%-5.2f+i(%-5.2f)Ohm\n',real(Ym),imag(Ym));
fprintf(fid,'\n');
elseif pfaltam==19
fprintf(fid,'Falta m en Linea SY, entre E1 y E2\n');
q=input('m=(0...1)=');
m=q;
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Sm=zl1u*(L1+m*L2); Ym=zl1u*(L5+(1-m)*L2);
Sm=Sm*zbase; Ym=Ym*zbase;
fprintf(fid,'Sm =%-5.2f+i(%-5.2f)Ohm\n',real(Sm),imag(Sm));
fprintf(fid,'Ym =%-5.2f+i(%-5.2f)Ohm\n',real(Ym),imag(Ym));
fprintf('\n');
elseif pfaltam==20
fprintf(fid,'Falta m en Linea SY, entre S y E1\n');
r=input('m=(0...1)=');
m=r;
if fid==1
else
fprintf(fid,'m=%4.2f\n',m);
end
Sm=zl1u*m*L1; Ym=zl1u*(L5+L2+(1-m)*L1);
Sm=Sm*zbase; Ym=Ym*zbase;
fprintf(fid,'Sm =%-5.2f+i(%-5.2f)Ohm\n',real(Sm),imag(Sm));
fprintf(fid,'Ym =%-5.2f+i(%-5.2f)Ohm\n',real(Ym),imag(Ym));
fprintf(fid,'\n');
end
if TipoFalta==2
fprintf(fid,'Tramo con Falta p: T1sup=11,T4izq=12,T4der=13,T2sup=14,T3sup=15,T3inf=16\n');
fprintf(fid,'T5der=17,T5izq=18,T2inf=19,T1inf=20: ');
if salida==1
pfaltap=input('Tramo con Falta p: ');
elseif salida==2
fprintf('Tramo con Falta p: T1sup=11,T4izq=12,T4der=13,T2sup=14,T3sup=15,T3inf=16\n');
fprintf('T5der=17,T5izq=18,T2inf=19,T1inf=20 \n');
pfaltap=input('Tramo con Falta p: ');
fprintf(fid,'%3.0f\n',pfaltap);
end
if pfaltap==11
fprintf(fid,'Falta p en Linea SX, entre S y E1\n');
a=input('p =(0...1)=');
p=a;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Sp=zl1u*p*L1; Xp=zl1u*(L4+(1-p)*L1);
Sp=Sp*zbase; Xp=Xp*zbase;
fprintf(fid,'Sp =%-5.2f+i(%-5.2f)Ohm\n',real(Sp),imag(Sp));
fprintf(fid,'Xp =%-5.2f+i(%-5.2f)Ohm\n',real(Xp),imag(Xp));
fprintf(fid,'\n');
elseif pfaltap==12
fprintf(fid,'Falta p en Linea SX, entre E1 y X\n');
b=input('p=(0...1)=');
p=b;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Sp=zl1u*(L1+p*L4); Xp=zl1u*(1-p)*L4;
Sp=Sp*zbase; Xp=Xp*zbase;
fprintf(fid,'Sp =%-5.2f+i(%-5.2f)Ohm\n',real(Sp),imag(Sp));
fprintf(fid,'Xp =%-5.2f+i(%-5.2f)Ohm\n',real(Xp),imag(Xp));
fprintf(fid,'\n');
elseif pfaltap==13
fprintf(fid,'Falta p en Linea RX, entre E1 y X\n');
c=input('p=(0...1)=');
p=c;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Xp=zl1u*(1-p)*L4; Rp=zl1u*(L3+L2+p*L4);
Xp=Xp*zbase; Rp=Rp*zbase;
fprintf(fid,'Xp =%-5.2f+i(%-5.2f)Ohm\n',real(Xp),imag(Xp));
fprintf(fid,'Rp =%-5.2f+i(%-5.2f)Ohm\n',real(Rp),imag(Rp));
fprintf(fid,'\n');
elseif pfaltap==14
fprintf(fid,'Falta p en Linea RX, entre E1 y E2\n');
d=input('p=(0...1)=');
p=d;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Xp=zl1u*(L4+p*L2); Rp=zl1u*(L3+(1-p)*L2);
Xp=Xp*zbase; Rp=Rp*zbase;
fprintf(fid,'Xp =%-5.2f+i(%-5.2f)Ohm\n',real(Xp),imag(Xp));
fprintf(fid,'Rp =%-5.2f+i(%-5.2f)Ohm\n',real(Rp),imag(Rp));
fprintf(fid,'\n');
elseif pfaltap==15
fprintf(fid,'Falta p en Linea RX, entre E2 y R\n');
e=input('p=(0...1)=');
p=e;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Xp=zl1u*(L4+L2+p*L3); Rp=zl1u*(1-p)*L3;
Xp=Xp*zbase; Rp=Rp*zbase;
fprintf(fid,'Xp =%-5.2f+i(%-5.2f)Ohm\n',real(Xp),imag(Xp));
fprintf(fid,'Rp =%-5.2f+i(%-5.2f)Ohm\n',real(Rp),imag(Rp));
fprintf(fid,'\n');
elseif pfaltap==16
fprintf(fid,'Falta p en Linea RY, entre E2 y R\n');
f=input('p=(0...1)=');
p=f;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Rp=zl1u*(1-p)*L3; Yp=zl1u*(L5+p*L3);
Rp=Rp*zbase; Yp=Yp*zbase;
fprintf(fid,'Rp =%-5.2f+i(%-5.2f)Ohm',real(Rp),imag(Rp));
fprintf(fid,'Yp =%-5.2f+i(%-5.2f)Ohm\n',real(Yp),imag(Yp));
fprintf(fid,'\n');
elseif pfaltap==17
fprintf(fid,'Falta p en Linea RY, entre E2 y Y\n');
g=input('p=(0...1)=');
p=g;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Rp=zl1u*(L3+p*L5); Yp=zl1u*(1-p)*L5;
Rp=Rp*zbase; Yp=Yp*zbase;
fprintf(fid,'Rp =%-5.2f+i(%-5.2f)Ohm\n',real(Rp),imag(Rp));
fprintf(fid,'\tYp =%-5.2f+i(%-5.2f)Ohm\n',real(Yp),imag(Yp));
fprintf(fid,'\n');
elseif pfaltap==18
fprintf(fid,'Falta p en Linea SY, entre E2 y Y\n');
h=input('p=(0...1)=');
p=h;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Sp=zl1u*(L1+L2+p*L5); Yp=zl1u*(1-p)*L5;
Sp=Sp*zbase; Yp=Yp*zbase;
fprintf(fid,'Sp =%-5.2f+i(%-5.2f)Ohm\n',real(Sp),imag(Sp));
fprintf(fid,'Yp =%-5.2f+i(%-5.2f)Ohm\n',real(Yp),imag(Yp));
fprintf(fid,'\n');
elseif pfaltap==19
fprintf(fid,'Falta p en Linea SY, entre E1 y E2\n');
q=input('p=(0...1)=');
p=q;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Sp=zl1u*(L1+p*L2); Yp=zl1u*(L5+(1-p)*L2);
Sp=Sp*zbase; Yp=Yp*zbase;
fprintf(fid,'Sp =%-5.2f+i(%-5.2f)Ohm\n',real(Sp),imag(Sp));
fprintf(fid,'Yp =%-5.2f+i(%-5.2f)Ohm\n',real(Yp),imag(Yp));
fprintf(fid,'\n');
elseif pfaltap==20
fprintf(fid,'Falta p en Linea SY, entre S y E1\n');
r=input('p=(0...1)=');
p=r;
if fid==1
else
fprintf(fid,'p=%4.2f\n',p);
end
Sp=zl1u*p*L1; Yp=zl1u*(L5+L2+(1-p)*L1);
Sp=Sp*zbase; Yp=Yp*zbase;
fprintf(fid,'Sp = %5.2f+i(%-5.2f)Ohm\n',real(Sp),imag(Sp));
fprintf(fid,'Yp = %5.2f+i(%-5.2f)Ohm\n',real(Yp),imag(Yp));
fprintf(fid,'\n');
end
end % Tipo de Falta ==2
% ----------------------------------------------------------------------
% CALCULOS
% ----------------------------------------------------------------------
[Minc]=subMinc(1);
[Minca]=subMinca(1);
[z0]=subz0(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r); % pu
[z1]=subz1(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r); % pu
[z2]=subz2(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r); % pu
% 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
if TipoFalta==1
% CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p 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 % TipoFalta==2
fprintf(fid,'\n');
fprintf(fid,'m =%-3.2f\n',m);
fprintf(fid,'Falta:%s\n',cfaltam);
fprintf(fid,'Rfm:%-10.2f Ohm\n',Rfm*zbase);
if (faltam==7)|(faltam==8)|(faltam==9)|(faltam==11)
fprintf(fid,'Rgm:%-10.2f Ohm\n',Rgm*zbase);
end
if TipoFalta==2
fprintf(fid,'\n');
fprintf(fid,'p =%-3.2f\n',p);
fprintf(fid,'Falta:%s\n',cfaltap);
fprintf(fid,'Rfp:%-10.2f Ohm\n',Rfp*zbase);
if (faltap==7)|(faltap==8)|(faltap==9)|(faltap==11)
fprintf(fid,'Rgp:%-10.2f Ohm\n',Rgp*zbase);
end
end
% IBUS(0,1,2)
I0=zeros(20,1); I1=zeros(20,1); I2=zeros(20,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;
% I(0,1,2) DE LOS ELEMENTOS
[Ie0]=subIe0(Minc,V0,z0);
[Ie1]=subIe1(Minca,V1a,z1);
[Ie2]=subIe2(Minc,V2,z2);
Vbase=Vbase*1000;
Ibase=Ibase*1000;
% ASIGNACIONES V,I
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;
Vx0=V0(2);Vx1=V1(2);Vx2=V2(2); VVx0=Vx0*Vbase; VVx1=Vx1*Vbase; VVx2=Vx2*Vbase;
Vr0=V0(3);Vr1=V1(3);Vr2=V2(3); VVr0=Vr0*Vbase; VVr1=Vr1*Vbase; VVr2=Vr2*Vbase;
Vy0=V0(4);Vy1=V1(4);Vy2=V2(4); VVy0=Vy0*Vbase; VVy1=Vy1*Vbase; VVy2=Vy2*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;
Ix0=Ie0(2);Ix1=Ie1(2);Ix2=Ie2(2); IIx0=Ix0*Ibase; IIx1=Ix1*Ibase; IIx2=Ix2*Ibase;
Ir0=Ie0(3);Ir1=Ie1(3);Ir2=Ie2(3); IIr0=Ir0*Ibase; IIr1=Ir1*Ibase; IIr2=Ir2*Ibase;
Iy0=Ie0(4);Iy1=Ie1(4);Iy2=Ie2(4); IIy0=Iy0*Ibase; IIy1=Iy1*Ibase; IIy2=Iy2*Ibase;
Isx0=Ie0(5);Isx1=Ie1(5);Isx2=Ie2(5); IIsx0=Isx0*Ibase; IIsx1=Isx1*Ibase; IIsx2=Isx2*Ibase;
Isy0=Ie0(6);Isy1=Ie1(6);Isy2=Ie2(6); IIsy0=Isy0*Ibase; IIsy1=Isy1*Ibase; IIsy2=Isy2*Ibase;
Ixs0=Ie0(7);Ixs1=Ie1(7);Ixs2=Ie2(7); IIxs0=Ixs0*Ibase; IIxs1=Ixs1*Ibase; IIxs2=Ixs2*Ibase;
Ixr0=Ie0(8);Ixr1=Ie1(8);Ixr2=Ie2(8); IIxr0=Ixr0*Ibase; IIxr1=Ixr1*Ibase; IIxr2=Ixr2*Ibase;
Irx0=Ie0(9);Irx1=Ie1(9);Irx2=Ie2(9); IIrx0=Irx0*Ibase; IIrx1=Irx1*Ibase; IIrx2=Irx2*Ibase;
Iry0=Ie0(10);Iry1=Ie1(10);Iry2=Ie2(10); IIry0=Iry0*Ibase; IIry1=Iry1*Ibase; IIry2=Iry2*Ibase;
Iys0=Ie0(11);Iys1=Ie1(11);Iys2=Ie2(11); IIys0=Iys0*Ibase; IIys1=Iys1*Ibase; IIys2=Iys2*Ibase;
Iyr0=Ie0(12);Iyr1=Ie1(12);Iyr2=Ie2(12); IIyr0=Iyr0*Ibase; IIyr1=Iyr1*Ibase; IIyr2=Iyr2*Ibase;
% CALCULO DE LAS COMPONENTES DE FASE DE I,V EN LOS PUNTOS m,p DE FALTA
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);
if TipoFalta==2
[Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2);
end
[Isa,Isb,Isc]=abc(Is0,Is1,Is2);
[Ixa,Ixb,Ixc]=abc(Ix0,Ix1,Ix2);
[Ira,Irb,Irc]=abc(Ir0,Ir1,Ir2);
[Iya,Iyb,Iyc]=abc(Iy0,Iy1,Iy2);
[Isxa,Isxb,Isxc]=abc(Isx0,Isx1,Isx2);
[Isya,Isyb,Isyc]=abc(Isy0,Isy1,Isy2);
[Ixsa,Ixsb,Ixsc]=abc(Ixs0,Ixs1,Ixs2);
[Ixra,Ixrb,Ixrc]=abc(Ixr0,Ixr1,Ixr2);
[Irxa,Irxb,Irxc]=abc(Irx0,Irx1,Irx2);
[Irya,Iryb,Iryc]=abc(Iry0,Iry1,Iry2);
[Iyra,Iyrb,Iyrc]=abc(Iyr0,Iyr1,Iyr2);
[Iysa,Iysb,Iysc]=abc(Iys0,Iys1,Iys2);
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
if TipoFalta==2
[Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2);
end
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vxa,Vxb,Vxc]=abc(Vx0,Vx1,Vx2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[Vya,Vyb,Vyc]=abc(Vy0,Vy1,Vy2);
[VVma,VVmb,VVmc]=abc(VVm0,VVm1,VVm2);
if TipoFalta==2
[VVpa,VVpb,VVpc]=abc(VVp0,VVp1,VVp2);
end
[VVsa,VVsb,VVsc]=abc(VVs0,VVs1,VVs2);
[VVxa,VVxb,VVxc]=abc(VVx0,VVx1,VVx2);
[VVra,VVrb,VVrc]=abc(VVr0,VVr1,VVr2);
[VVya,VVyb,VVyc]=abc(VVy0,VVy1,VVy2);
[IIma,IImb,IImc]=abc(IIm0,IIm1,IIm2);
if TipoFalta==2
[IIpa,IIpb,IIpc]=abc(IIp0,IIp1,IIp2);
end
[IIsa,IIsb,IIsc]=abc(IIs0,IIs1,IIs2);
[IIxa,IIxb,IIxc]=abc(IIx0,IIx1,IIx2);
[IIra,IIrb,IIrc]=abc(IIr0,IIr1,IIr2);
[IIya,IIyb,IIyc]=abc(IIy0,IIy1,IIy2);
[IIsxa,IIsxb,IIsxc]=abc(IIsx0,IIsx1,IIsx2);
[IIsya,IIsyb,IIsyc]=abc(IIsy0,IIsy1,IIsy2);
[IIxsa,IIxsb,IIxsc]=abc(IIxs0,IIxs1,IIxs2);
[IIxra,IIxrb,IIxrc]=abc(IIxr0,IIxr1,IIxr2);
[IIrxa,IIrxb,IIrxc]=abc(IIrx0,IIrx1,IIrx2);
[IIrya,IIryb,IIryc]=abc(IIry0,IIry1,IIry2);
[IIyra,IIyrb,IIyrc]=abc(IIyr0,IIyr1,IIyr2);
[IIysa,IIysb,IIysc]=abc(IIys0,IIys1,IIys2);
% PROTECCION DE DISTANCIA (21)
switch(compens21)
case{1},
% COMPENSADA
[Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isxa,Isxb,Isxc,Isx0,k,Km,Isy0);
[Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isya,Isyb,Isyc,Isy0,k,Km,Isx0);
[Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG]=...
PDistComp(Vxa,Vxb,Vxc,Ixsa,Ixsb,Ixsc,Ixs0,k,Km,Ixr0);
[Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG]=...
PDistComp(Vxa,Vxb,Vxc,Ixra,Ixrb,Ixrc,Ixr0,k,Km,Ixs0);
[Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irxa,Irxb,Irxc,Irx0,k,Km,Iry0);
[Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irya,Iryb,Iryc,Iry0,k,Km,Irx0);
[Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG]=...
PDistComp(Vya,Vyb,Vyc,Iyra,Iyrb,Iyrc,Iyr0,k,Km,Iys0);
[Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG]=...
PDistComp(Vya,Vyb,Vyc,Iysa,Iysb,Iysc,Iys0,k,Km,Iyr0);
case{2},
% NO COMPENSADA
[Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isxa,Isxb,Isxc,Isx0,k);
[Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isya,Isyb,Isyc,Isy0,k);
[Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG]=...
PDist(Vxa,Vxb,Vxc,Ixsa,Ixsb,Ixsc,Ixs0,k);
[Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG]=...
PDist(Vxa,Vxb,Vxc,Ixra,Ixrb,Ixrc,Ixr0,k);
[Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG]=...
PDist(Vra,Vrb,Vrc,Irxa,Irxb,Irxc,Irx0,k);
[Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG]=...
PDist(Vra,Vrb,Vrc,Irya,Iryb,Iryc,Iry0,k);
[Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG]=...
PDist(Vya,Vyb,Vyc,Iyra,Iyrb,Iyrc,Iyr0,k);
[Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG]=...
PDist(Vya,Vyb,Vyc,Iysa,Iysb,Iysc,Iys0,k);
end
% 67N
[Vspol,Isxop,AngDirSx]=PDire(Vs0,Isx0);
[Vspol,Isyop,AngDirSy]=PDire(Vs0,Isy0);
[Vxpol,Ixsop,AngDirXs]=PDire(Vx0,Ixs0);
[Vxpol,Ixrop,AngDirXr]=PDire(Vx0,Ixr0);
[Vrpol,Irxop,AngDirRx]=PDire(Vr0,Irx0);
[Vrpol,Iryop,AngDirRy]=PDire(Vr0,Iry0);
[Vypol,Iyrop,AngDirYr]=PDire(Vy0,Iyr0);
[Vypol,Iysop,AngDirYs]=PDire(Vy0,Iys0);
[VVspol,IIsxop,AngDirSx]=PDire(VVs0,IIsx0);
[VVspol,IIsyop,AngDirSy]=PDire(VVs0,IIsy0);
[VVxpol,IIxsop,AngDirXs]=PDire(VVx0,IIxs0);
[VVxpol,IIxrop,AngDirXr]=PDire(VVx0,IIxr0);
[VVrpol,IIrxop,AngDirRx]=PDire(VVr0,IIrx0);
[VVrpol,IIryop,AngDirRy]=PDire(VVr0,IIry0);
[VVypol,IIyrop,AngDirYr]=PDire(VVy0,IIyr0);
[VVypol,IIysop,AngDirYs]=PDire(VVy0,IIys0);
% ---------------------------------------------------
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 67N, pu, cartesiana *****************
% -----------------------------------------------------------------------
% Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys
% -----------------------------------------------------------------------
% -----------------------------------------------------------------------
% ** C1 ************** SALIDA V(0,1,2),I(0,1,2) pu cartesiana ***********
% -----------------------------------------------------------------------
% Vm,Vp_2,Vs,Vx,Vr,Vy
% Im,Ip_2,Is,Ix,Ir,Iy,Isx,Isy,Ixs,Ixr,Irx,Iry,Ir,Iys
% -----------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'Vm0 :%9.2f+i(%-9.2f)pu',real(Vm0),imag(Vm0));
fprintf(fid,'\tVm1 :%9.2f+i(%-9.2f)pu',real(Vm1),imag(Vm1));
fprintf(fid,'\tVm2 :%9.2f+i(%-9.2f)pu',real(Vm2),imag(Vm2));
fprintf(fid,'\n');
if TipoFalta==2
fprintf(fid,'Vp0 :%9.2f+i(%-9.2f)pu',real(Vp0),imag(Vp0));
fprintf(fid,'\tVp1 :%9.2f+i(%-9.2f)pu',real(Vp1),imag(Vp1));
fprintf(fid,'\tVp2 :%9.2f+i(%-9.2f)pu',real(Vp2),imag(Vp2));
fprintf(fid,'\n');
end
fprintf(fid,'Vs0 :%9.2f+i(%-9.2f)pu',real(Vs0),imag(Vs0));
fprintf(fid,'\tVs1 :%9.2f+i(%-9.2f)pu',real(Vs1),imag(Vs1));
fprintf(fid,'\tVs2 :%9.2f+i(%-9.2f)pu',real(Vs2),imag(Vs2));
fprintf(fid,'\n');
fprintf(fid,'Vx0 :%9.2f+i(%-9.2f)pu',real(Vx0),imag(Vx0));
fprintf(fid,'\tVx1 :%9.2f+i(%-9.2f)pu',real(Vx1),imag(Vx1));
fprintf(fid,'\tVx2 :%9.2f+i(%-9.2f)pu',real(Vx2),imag(Vx2));
fprintf(fid,'\n');
fprintf(fid,'Vr0 :%9.2f+i(%-9.2f)pu',real(Vr0),imag(Vr0));
fprintf(fid,'\tVr1 :%9.2f+i(%-9.2f)pu',real(Vr1),imag(Vr1));
fprintf(fid,'\tVr2 :%9.2f+i(%-9.2f)pu',real(Vr2),imag(Vr2));
fprintf(fid,'\n');
fprintf(fid,'Vy0 :%9.2f+i(%-9.2f)pu',real(Vy0),imag(Vy0));
fprintf(fid,'\tVy1 :%9.2f+i(%-9.2f)pu',real(Vy1),imag(Vy1));
fprintf(fid,'\tVy2 :%9.2f+i(%-9.2f)pu',real(Vy2),imag(Vy2));
fprintf(fid,'\n\n');
% ----------------------------------------------------------------------
fprintf(fid,'Im0 :%9.2f+i(%-9.2f)pu',real(Im0),imag(Im0));
fprintf(fid,'\tIm1 :%9.2f+i(%-9.2f)pu',real(Im1),imag(Im1));
fprintf(fid,'\tIm2 :%9.2f+i(%-9.2f)pu',real(Im2),imag(Im2));
fprintf(fid,'\n');
if TipoFalta==2
fprintf(fid,'Ip0 :%9.2f+i(%-9.2f)pu',real(Ip0),imag(Ip0));
fprintf(fid,'\tIp1 :%9.2f+i(%-9.2f)pu',real(Ip1),imag(Ip1));
fprintf(fid,'\tIp2 :%9.2f+i(%-9.2f)pu',real(Ip2),imag(Ip2));
fprintf(fid,'\n');
end
fprintf(fid,'Is0 :%9.2f+i(%-9.2f)pu',real(Is0),imag(Is0));
fprintf(fid,'\tIs1 :%9.2f+i(%-9.2f)pu',real(Is1),imag(Is1));
fprintf(fid,'\tIs2 :%9.2f+i(%-9.2f)pu',real(Is2),imag(Is2));
fprintf(fid,'\n');
fprintf(fid,'Ix0 :%9.2f+i(%-9.2f)pu',real(Ix0),imag(Ix0));
fprintf(fid,'\tIx1 :%9.2f+i(%-9.2f)pu',real(Ix1),imag(Ix1));
fprintf(fid,'\tIx2 :%9.2f+i(%-9.2f)pu',real(Ix2),imag(Ix2));
fprintf(fid,'\n');
fprintf(fid,'Ir0 :%9.2f+i(%-9.2f)pu',real(Ir0),imag(Ir0));
fprintf(fid,'\tIr1 :%9.2f+i(%-9.2f)pu',real(Ir1),imag(Ir1));
fprintf(fid,'\tIr2 :%9.2f+i(%-9.2f)pu',real(Ir2),imag(Ir2));
fprintf(fid,'\n');
fprintf(fid,'Iy0 :%9.2f+i(%-9.2f)pu',real(Iy0),imag(Iy0));
fprintf(fid,'\tIy1 :%9.2f+i(%-9.2f)pu',real(Iy1),imag(Iy1));
fprintf(fid,'\tIy2 :%9.2f+i(%-9.2f)pu',real(Iy2),imag(Iy2));
fprintf(fid,'\n');
fprintf(fid,'Isx0:%9.2f+i(%-9.2f)pu',real(Isx0),imag(Isx0));
fprintf(fid,'\tIsx1:%9.2f+i(%-9.2f)pu',real(Isx1),imag(Isx1));
fprintf(fid,'\tIsx2:%9.2f+i(%-9.2f)pu',real(Isx2),imag(Isx2));
fprintf(fid,'\n');
fprintf(fid,'Isy0:%9.2f+i(%-9.2f)pu',real(Isy0),imag(Isy0));
fprintf(fid,'\tIsy1:%9.2f+i(%-9.2f)pu',real(Isy1),imag(Isy1));
fprintf(fid,'\tIsy2:%9.2f+i(%-9.2f)pu',real(Isy2),imag(Isy2));
fprintf(fid,'\n');
fprintf(fid,'Ixs0:%9.2f+i(%-9.2f)pu',real(Ixs0),imag(Ixs0));
fprintf(fid,'\tIxs1:%9.2f+i(%-9.2f)pu',real(Ixs1),imag(Ixs1));
fprintf(fid,'\tIxs2:%9.2f+i(%-9.2f)pu',real(Ixs2),imag(Ixs2));
fprintf(fid,'\n');
fprintf(fid,'Ixr0:%9.2f+i(%-9.2f)pu',real(Ixr0),imag(Ixr0));
fprintf(fid,'\tIxr1:%9.2f+i(%-9.2f)pu',real(Ixr1),imag(Ixr1));
fprintf(fid,'\tIxr2:%9.2f+i(%-9.2f)pu',real(Ixr2),imag(Ixr2));
fprintf(fid,'\n');
fprintf(fid,'Irx0:%9.2f+i(%-9.2f)pu',real(Irx0),imag(Irx0));
fprintf(fid,'\tIrx1:%9.2f+i(%-9.2f)pu',real(Irx1),imag(Irx1));
fprintf(fid,'\tIrx2:%9.2f+i(%-9.2f)pu',real(Irx2),imag(Irx2));
fprintf(fid,'\n');
fprintf(fid,'Iry0:%9.2f+i(%-9.2f)pu',real(Iry0),imag(Iry0));
fprintf(fid,'\tIry1:%9.2f+i(%-9.2f)pu',real(Iry1),imag(Iry1));
fprintf(fid,'\tIry2:%9.2f+i(%-9.2f)pu',real(Iry2),imag(Iry2));
fprintf(fid,'\n');
fprintf(fid,'Iyr0:%9.2f+i(%-9.2f)pu',real(Iyr0),imag(Iyr0));
fprintf(fid,'\tIyr1:%9.2f+i(%-9.2f)pu',real(Iyr1),imag(Iyr1));
fprintf(fid,'\tIyr2:%9.2f+i(%-9.2f)pu',real(Iyr2),imag(Iyr2));
fprintf(fid,'\n');
fprintf(fid,'Iys0:%9.2f+i(%-9.2f)pu',real(Iys0),imag(Iys0));
fprintf(fid,'\tIys1:%9.2f+i(%-9.2f)pu',real(Iys1),imag(Iys1));
fprintf(fid,'\tIys2:%9.2f+i(%-9.2f)pu',real(Iys2),imag(Iys2));
fprintf(fid,'\n');
fprintf(fid,'\n');
% ***********************************************************************************
% ------------------------------------------------------------------
% ** C1 *********** SALIDA V(a,b,c),I(a,b,c) pu cartesiana *********
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vx,Vr,Vy
% Im,Ip_2,Is,Ix,Ir,Iy,Isx,Isy,Ixs,Ixr,Irx,Iry,Ir,Iys
% ------------------------------------------------------------------
% ***********************************************************************************
fprintf(fid,'Vma :%9.2f+i(%-9.2f)pu',real(Vma),imag(Vma));
fprintf(fid,'\tVmb :%9.2f+i(%-9.2f)pu',real(Vmb),imag(Vmb));
fprintf(fid,'\tVmc :%9.2f+i(%-9.2f)pu',real(Vmc),imag(Vmc));
fprintf(fid,'\n');
if TipoFalta==2
fprintf(fid,'Vpa :%9.2f+i(%-9.2f)pu',real(Vpa),imag(Vpa));
fprintf(fid,'\tVpb :%9.2f+i(%-9.2f)pu',real(Vpb),imag(Vpb));
fprintf(fid,'\tVpc :%9.2f+i(%-9.2f)pu',real(Vpc),imag(Vpc));
fprintf(fid,'\n');
end
fprintf(fid,'Vsa :%9.2f+i(%-9.2f)pu',real(Vsa),imag(Vsa));
fprintf(fid,'\tVsb :%9.2f+i(%-9.2f)pu',real(Vsb),imag(Vsb));
fprintf(fid,'\tVsc :%9.2f+i(%-9.2f)pu',real(Vsc),imag(Vsc));
fprintf(fid,'\n');
fprintf(fid,'Vxa :%9.2f+i(%-9.2f)pu',real(Vxa),imag(Vxa));
fprintf(fid,'\tVxb :%9.2f+i(%-9.2f)pu',real(Vxb),imag(Vxb));
fprintf(fid,'\tVxc :%9.2f+i(%-9.2f)pu',real(Vxc),imag(Vxc));
fprintf(fid,'\n');
fprintf(fid,'Vra :%9.2f+i(%-9.2f)pu',real(Vra),imag(Vra));
fprintf(fid,'\tVrb :%9.2f+i(%-9.2f)pu',real(Vrb),imag(Vrb));
fprintf(fid,'\tVrc :%9.2f+i(%-9.2f)pu',real(Vrc),imag(Vrc));
fprintf(fid,'\n');
fprintf(fid,'Vya :%9.2f+i(%-9.2f)pu',real(Vya),imag(Vya));
fprintf(fid,'\tVyb :%9.2f+i(%-9.2f)pu',real(Vyb),imag(Vyb));
fprintf(fid,'\tVyc :%9.2f+i(%-9.2f)pu',real(Vyc),imag(Vyc));
fprintf(fid,'\n\n');
% ----------------------------------------------------------------------
fprintf(fid,'Ima :%9.2f+i(%-9.2f)pu',real(Ima),imag(Ima));
fprintf(fid,'\tImb :%9.2f+i(%-9.2f)pu',real(Imb),imag(Imb));
fprintf(fid,'\tImc :%9.2f+i(%-9.2f)pu',real(Imc),imag(Imc));
fprintf(fid,'\n');
if TipoFalta==2
fprintf(fid,'Ipa :%9.2f+i(%-9.2f)pu',real(Ipa),imag(Ipa));
fprintf(fid,'\tIpb :%9.2f+i(%-9.2f)pu',real(Ipb),imag(Ipb));
fprintf(fid,'\tIpc :%9.2f+i(%-9.2f)pu',real(Ipc),imag(Ipc));
fprintf(fid,'\n');
end
fprintf(fid,'Isa :%9.2f+i(%-9.2f)pu',real(Isa),imag(Isa));
fprintf(fid,'\tIsb :%9.2f+i(%-9.2f)pu',real(Isb),imag(Isb));
fprintf(fid,'\tIsc :%9.2f+i(%-9.2f)pu',real(Isc),imag(Isc));
fprintf(fid,'\n');
fprintf(fid,'Ixa :%9.2f+i(%-9.2f)pu',real(Ixa),imag(Ixa));
fprintf(fid,'\tIxb :%9.2f+i(%-9.2f)pu',real(Ixb),imag(Ixb));
fprintf(fid,'\tIxc :%9.2f+i(%-9.2f)pu',real(Ixc),imag(Ixc));
fprintf(fid,'\n');
fprintf(fid,'Ira :%9.2f+i(%-9.2f)pu',real(Ira),imag(Ira));
fprintf(fid,'\tIrb :%9.2f+i(%-9.2f)pu',real(Irb),imag(Irb));
fprintf(fid,'\tIrc :%9.2f+i(%-9.2f)pu',real(Irc),imag(Irc));
fprintf(fid,'\n');
fprintf(fid,'Iya :%9.2f+i(%-9.2f)pu',real(Iya),imag(Iya));
fprintf(fid,'\tIyb :%9.2f+i(%-9.2f)pu',real(Iyb),imag(Iyb));
fprintf(fid,'\tIyc :%9.2f+i(%-9.2f)pu',real(Iyc),imag(Iyc));
fprintf(fid,'\n\n');
fprintf(fid,'Isxa:%9.2f+i(%-9.2f)pu',real(Isxa),imag(Isxa));
fprintf(fid,'\tIsxb:%9.2f+i(%-9.2f)pu',real(Isxb),imag(Isxb));
fprintf(fid,'\tIsxc:%9.2f+i(%-9.2f)pu',real(Isxc),imag(Isxc));
fprintf(fid,'\n');
fprintf(fid,'Isya:%9.2f+i(%-9.2f)pu',real(Isya),imag(Isya));
fprintf(fid,'\tIsyb:%9.2f+i(%-9.2f)pu',real(Isyb),imag(Isyb));
fprintf(fid,'\tIsyc:%9.2f+i(%-9.2f)pu',real(Isyc),imag(Isyc));
fprintf(fid,'\n');
fprintf(fid,'Ixsa:%9.2f+i(%-9.2f)pu',real(Ixsa),imag(Ixsa));
fprintf(fid,'\tIxsb:%9.2f+i(%-9.2f)pu',real(Ixsb),imag(Ixsb));
fprintf(fid,'\tIxsc:%9.2f+i(%-9.2f)pu',real(Ixsc),imag(Ixsc));
fprintf(fid,'\n');
fprintf(fid,'Ixra:%9.2f+i(%-9.2f)pu',real(Ixra),imag(Ixra));
fprintf(fid,'\tIxrb:%9.2f+i(%-9.2f)pu',real(Ixrb),imag(Ixrb));
fprintf(fid,'\tIxrc:%9.2f+i(%-9.2f)pu',real(Ixrc),imag(Ixrc));
fprintf(fid,'\n');
fprintf(fid,'Irxa:%9.2f+i(%-9.2f)pu',real(Irxa),imag(Irxa));
fprintf(fid,'\tIrxb:%9.2f+i(%-9.2f)pu',real(Irxb),imag(Irxb));
fprintf(fid,'\tIrxc:%9.2f+i(%-9.2f)pu',real(Irxc),imag(Irxc));
fprintf(fid,'\n');
fprintf(fid,'Irya:%9.2f+i(%-9.2f)pu',real(Irya),imag(Irya));
fprintf(fid,'\tIryb:%9.2f+i(%-9.2f)pu',real(Iryb),imag(Iryb));
fprintf(fid,'\tIryc:%9.2f+i(%-9.2f)pu',real(Iryc),imag(Iryc));
fprintf(fid,'\n');
fprintf(fid,'Iyra:%9.2f+i(%-9.2f)pu',real(Iyra),imag(Iyra));
fprintf(fid,'\tIyrb:%9.2f+i(%-9.2f)pu',real(Iyrb),imag(Iyrb));
fprintf(fid,'\tIyrc:%9.2f+i(%-9.2f)pu',real(Iyrc),imag(Iyrc));
fprintf(fid,'\n');
fprintf(fid,'Iysa:%9.2f+i(%-9.2f)pu',real(Iysa),imag(Iysa));
fprintf(fid,'\tIysb:%9.2f+i(%-9.2f)pu',real(Iysb),imag(Iysb));
fprintf(fid,'\tIysc:%9.2f+i(%-9.2f)pu',real(Iysc),imag(Iysc));
fprintf(fid,'\n\n');
% **********************************************************************
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 ********************
% ------------------------------------------------------------------
% Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'Sx_67N:\tVspol=%-8.4f pu\tIsxop=%-8.4f pu\tAngDirSx=%-5.1f\n',Vspol,Isxop,AngDirSx);
fprintf(fid,'Sy_67N:\tVspol=%-8.4f pu\tIsyop=%-8.4f pu\tAngDirSy=%-5.1f\n',Vspol,Isyop,AngDirSy);
fprintf(fid,'Xs_67N:\tVxpol=%-8.4f pu\tIxsop=%-8.4f pu\tAngDirXs=%-5.1f\n',Vxpol,Ixsop,AngDirXs);
fprintf(fid,'Xr_67N:\tVxpol=%-8.4f pu\tIxrop=%-8.4f pu\tAngDirXr=%-5.1f\n',Vxpol,Ixrop,AngDirXr);
fprintf(fid,'Rx_67N:\tVrpol=%-8.4f pu\tIrxop=%-8.4f pu\tAngDirRx=%-5.1f\n',Vrpol,Irxop,AngDirRx);
fprintf(fid,'Ry_67N:\tVrpol=%-8.4f pu\tIryop=%-8.4f pu\tAngDirRy=%-5.1f\n',Vrpol,Iryop,AngDirRy);
fprintf(fid,'Yr_67N:\tVypol=%-8.4f pu\tIyrop=%-8.4f pu\tAngDirYr=%-5.1f\n',Vypol,Iyrop,AngDirYr);
fprintf(fid,'Ys_67N:\tVypol=%-8.4f pu\tIysop=%-8.4f pu\tAngDirYs=%-5.1f\n',Vypol,Iysop,AngDirYs);
fprintf(fid,'\n');
% ------------------------------------------------------------------
% ** C2 ************ SALIDA V(0,1,2), I(0,1,2) pu Polar ************
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vx,Vr,Vy
% Im,Ip_2,Is,Ix,Ir,Iy,Isx,Isy,Ixs,Ixr,Irx,Iry,Ir,Iys
% ------------------------------------------------------------------
% ------------------------------------------------------------------
% ** C2 ************ SALIDA V(a,b,c), I(a,b,c) pu Polar ************
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vx,Vr,Vy
% Im,Ip_2,Is,Ix,Ir,Iy,Isx,Isy,Ixs,Ixr,Irx,Iry,Ir,Iys
% ------------------------------------------------------------------
% **********************************************************************
case{3}, % OVA(cartesiana): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
% --
% **********************************************************************
% -----------------------------------------------------------------------
% ** C3 ******************** SALIDA 21, Ohm, cartesiana *****************
% -----------------------------------------------------------------------
% Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys
% -----------------------------------------------------------------------
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');
[Sx_ZABre,Sx_ZABim]=zmed(Sx_ZAB,zbase); fprintf(fid,'Sx\tAB:%8.2f+i(%-8.2f)Ohm',Sx_ZABre,Sx_ZABim);
[Sx_ZBCre,Sx_ZBCim]=zmed(Sx_ZBC,zbase); fprintf(fid,'\tBC:%8.2f+i(%-8.2f)Ohm',Sx_ZBCre,Sx_ZBCim);
[Sx_ZCAre,Sx_ZCAim]=zmed(Sx_ZCA,zbase); fprintf(fid,'\tCA:%8.2f+i(%-8.2f)Ohm\n',Sx_ZCAre,Sx_ZCAim);
[Sx_ZAGre,Sx_ZAGim]=zmed(Sx_ZAG,zbase); fprintf(fid,'\tAG:%8.2f+i(%-8.2f)Ohm',Sx_ZAGre,Sx_ZAGim);
[Sx_ZBGre,Sx_ZBGim]=zmed(Sx_ZBG,zbase); fprintf(fid,'\tBG:%8.2f+i(%-8.2f)Ohm',Sx_ZBGre,Sx_ZBGim);
[Sx_ZCGre,Sx_ZCGim]=zmed(Sx_ZCG,zbase); fprintf(fid,'\tCG:%8.2f+i(%-8.2f)Ohm\n',Sx_ZCGre,Sx_ZCGim);
fprintf(fid,'\n');
[Sy_ZABre,Sy_ZABim]=zmed(Sy_ZAB,zbase);fprintf(fid,'Sy\tAB:%8.2f+i(%-8.2f)Ohm',Sy_ZABre,Sy_ZABim);
[Sy_ZBCre,Sy_ZBCim]=zmed(Sy_ZBC,zbase);fprintf(fid,'\tBC:%8.2f+i(%-8.2f)Ohm',Sy_ZBCre,Sy_ZBCim);
[Sy_ZCAre,Sy_ZCAim]=zmed(Sy_ZCA,zbase);fprintf(fid,'\tCA:%8.2f+i(%-8.2f)Ohm\n',Sy_ZCAre,Sy_ZCAim);
[Sy_ZAGre,Sy_ZAGim]=zmed(Sy_ZAG,zbase);fprintf(fid,'\tAG:%8.2f+i(%-8.2f)Ohm',Sy_ZAGre,Sy_ZAGim);
[Sy_ZBGre,Sy_ZBGim]=zmed(Sy_ZBG,zbase);fprintf(fid,'\tBG:%8.2f+i(%-8.2f)Ohm',Sy_ZBGre,Sy_ZBGim);
[Sy_ZCGre,Sy_ZCGim]=zmed(Sy_ZCG,zbase);fprintf(fid,'\tCG:%8.2f+i(%-8.2f)Ohm\n',Sy_ZCGre,Sy_ZCGim);
fprintf(fid,'\n');
[Xs_ZABre,Xs_ZABim]=zmed(Xs_ZAB,zbase);fprintf(fid,'Xs\tAB:%8.2f+i(%-8.2f)Ohm',Xs_ZABre,Xs_ZABim);
[Xs_ZBCre,Xs_ZBCim]=zmed(Xs_ZBC,zbase);fprintf(fid,'\tBC:%8.2f+i(%-8.2f)Ohm',Xs_ZBCre,Xs_ZBCim);
[Xs_ZCAre,Xs_ZCAim]=zmed(Xs_ZCA,zbase);fprintf(fid,'\tCA:%8.2f+i(%-8.2f)Ohm\n',Xs_ZCAre,Xs_ZCAim);
[Xs_ZAGre,Xs_ZAGim]=zmed(Xs_ZAG,zbase);fprintf(fid,'\tAG:%8.2f+i(%-8.2f)Ohm',Xs_ZAGre,Xs_ZAGim);
[Xs_ZBGre,Xs_ZBGim]=zmed(Xs_ZBG,zbase);fprintf(fid,'\tBG:%8.2f+i(%-8.2f)Ohm',Xs_ZBGre,Xs_ZBGim);
[Xs_ZCGre,Xs_ZCGim]=zmed(Xs_ZCG,zbase);fprintf(fid,'\tCG:%8.2f+i(%-8.2f)Ohm\n',Xs_ZCGre,Xs_ZCGim);
fprintf(fid,'\n');
[Xr_ZABre,Xr_ZABim]=zmed(Xr_ZAB,zbase);fprintf(fid,'Xr\tAB:%8.2f+i(%-8.2f)Ohm',Xr_ZABre,Xr_ZABim);
[Xr_ZBCre,Xr_ZBCim]=zmed(Xr_ZBC,zbase);fprintf(fid,'\tBC:%8.2f+i(%-8.2f)Ohm',Xr_ZBCre,Xr_ZBCim);
[Xr_ZCAre,Xr_ZCAim]=zmed(Xr_ZCA,zbase);fprintf(fid,'\tCA:%8.2f+i(%-8.2f)Ohm\n',Xr_ZCAre,Xr_ZCAim);
[Xr_ZAGre,Xr_ZAGim]=zmed(Xr_ZAG,zbase);fprintf(fid,'\tAG:%8.2f+i(%-8.2f)Ohm',Xr_ZAGre,Xr_ZAGim);
[Xr_ZBGre,Xr_ZBGim]=zmed(Xr_ZBG,zbase);fprintf(fid,'\tBG:%8.2f+i(%-8.2f)Ohm',Xr_ZBGre,Xr_ZBGim);
[Xr_ZCGre,Xr_ZCGim]=zmed(Xr_ZCG,zbase);fprintf(fid,'\tCG:%8.2f+i(%-8.2f)Ohm\n',Xr_ZCGre,Xr_ZCGim);
fprintf(fid,'\n');
[Rx_ZABre,Rx_ZABim]=zmed(Rx_ZAB,zbase);fprintf(fid,'Rx\tAB:%8.2f+i(%-8.2f)Ohm',Rx_ZABre,Rx_ZABim);
[Rx_ZBCre,Rx_ZBCim]=zmed(Rx_ZBC,zbase);fprintf(fid,'\tBC:%8.2f+i(%-8.2f)Ohm',Rx_ZBCre,Rx_ZBCim);
[Rx_ZCAre,Rx_ZCAim]=zmed(Rx_ZCA,zbase);fprintf(fid,'\tCA:%8.2f+i(%-8.2f)Ohm\n',Rx_ZCAre,Rx_ZCAim);
[Rx_ZAGre,Rx_ZAGim]=zmed(Rx_ZAG,zbase);fprintf(fid,'\tAG:%8.2f+i(%-8.2f)Ohm',Rx_ZAGre,Rx_ZAGim);
[Rx_ZBGre,Rx_ZBGim]=zmed(Rx_ZBG,zbase);fprintf(fid,'\tBG:%8.2f+i(%-8.2f)Ohm',Rx_ZBGre,Rx_ZBGim);
[Rx_ZCGre,Rx_ZCGim]=zmed(Rx_ZCG,zbase);fprintf(fid,'\tCG:%8.2f+i(%-8.2f)Ohm\n',Rx_ZCGre,Rx_ZCGim);
fprintf(fid,'\n');
[Ry_ZABre,Ry_ZABim]=zmed(Ry_ZAB,zbase);fprintf(fid,'Ry\tAB:%8.2f+i(%-8.2f)Ohm',Ry_ZABre,Ry_ZABim);
[Ry_ZBCre,Ry_ZBCim]=zmed(Ry_ZBC,zbase);fprintf(fid,'\tBC:%8.2f+i(%-8.2f)Ohm',Ry_ZBCre,Ry_ZBCim);
[Ry_ZCAre,Ry_ZCAim]=zmed(Ry_ZCA,zbase);fprintf(fid,'\tCA:%8.2f+i(%-8.2f)Ohm\n',Ry_ZCAre,Ry_ZCAim);
[Ry_ZAGre,Ry_ZAGim]=zmed(Ry_ZAG,zbase);fprintf(fid,'\tAG:%8.2f+i(%-8.2f)Ohm',Ry_ZAGre,Ry_ZAGim);
[Ry_ZBGre,Ry_ZBGim]=zmed(Ry_ZBG,zbase);fprintf(fid,'\tBG:%8.2f+i(%-8.2f)Ohm',Ry_ZBGre,Ry_ZBGim);
[Ry_ZCGre,Ry_ZCGim]=zmed(Ry_ZCG,zbase);fprintf(fid,'\tCG:%8.2f+i(%-8.2f)Ohm\n',Ry_ZCGre,Ry_ZCGim);
fprintf(fid,'\n');
[Yr_ZABre,Yr_ZABim]=zmed(Yr_ZAB,zbase);fprintf(fid,'Yr\tAB:%8.2f+i(%-8.2f)Ohm',Yr_ZABre,Yr_ZABim);
[Yr_ZBCre,Yr_ZBCim]=zmed(Yr_ZBC,zbase);fprintf(fid,'\tBC:%8.2f+i(%-8.2f)Ohm',Yr_ZBCre,Yr_ZBCim);
[Yr_ZCAre,Yr_ZCAim]=zmed(Yr_ZCA,zbase);fprintf(fid,'\tCA:%8.2f+i(%-8.2f)Ohm\n',Yr_ZCAre,Yr_ZCAim);
[Yr_ZAGre,Yr_ZAGim]=zmed(Yr_ZAG,zbase);fprintf(fid,'\tAG:%8.2f+i(%-8.2f)Ohm',Yr_ZAGre,Yr_ZAGim);
[Yr_ZBGre,Yr_ZBGim]=zmed(Yr_ZBG,zbase);fprintf(fid,'\tBG:%8.2f+i(%-8.2f)Ohm',Yr_ZBGre,Yr_ZBGim);
[Yr_ZCGre,Yr_ZCGim]=zmed(Yr_ZCG,zbase);fprintf(fid,'\tCG:%8.2f+i(%-8.2f)Ohm\n',Yr_ZCGre,Yr_ZCGim);
fprintf(fid,'\n');
[Ys_ZABre,Ys_ZABim]=zmed(Ys_ZAB,zbase);fprintf(fid,'Ys\tAB:%8.2f+i(%-8.2f)Ohm',Ys_ZABre,Ys_ZABim);
[Ys_ZBCre,Ys_ZBCim]=zmed(Ys_ZBC,zbase);fprintf(fid,'\tBC:%8.2f+i(%-8.2f)Ohm',Ys_ZBCre,Ys_ZBCim);
[Ys_ZCAre,Ys_ZCAim]=zmed(Ys_ZCA,zbase);fprintf(fid,'\tCA:%8.2f+i(%-8.2f)Ohm\n',Ys_ZCAre,Ys_ZCAim);
[Ys_ZAGre,Ys_ZAGim]=zmed(Ys_ZAG,zbase);fprintf(fid,'\tAG:%8.2f+i(%-8.2f)Ohm',Ys_ZAGre,Ys_ZAGim);
[Ys_ZBGre,Ys_ZBGim]=zmed(Ys_ZBG,zbase);fprintf(fid,'\tBG:%8.2f+i(%-8.2f)Ohm',Ys_ZBGre,Ys_ZBGim);
[Ys_ZCGre,Ys_ZCGim]=zmed(Ys_ZCG,zbase);fprintf(fid,'\tCG:%8.2f+i(%-8.2f)Ohm\n',Ys_ZCGre,Ys_ZCGim);
fprintf(fid,'\n');
% ------------------------------------------------------------------
% ** C3 ******** SALIDA V(0,1,2),I(0,1,2) Volt-Amp cartesiana ******
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vx,Vr,Vy
% Im,Ip_2,Is,Ix,Ir,Iy,Isx,Isy,Ixs,Ixr,Irx,Iry,Ir,Iys
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'Vm0 :%9.2f+i(%-9.2f)(V)',real(VVm0),imag(VVm0));
fprintf(fid,'\tVm1 :%9.2f+i(%-9.2f)(V)',real(VVm1),imag(VVm1));
fprintf(fid,'\tVm2 :%9.2f+i(%-9.2f)(V)',real(VVm2),imag(VVm2));
fprintf(fid,'\n');
if TipoFalta==2
fprintf(fid,'Vp0 :%9.2f+i(%-9.2f)(V)',real(VVp0),imag(VVp0));
fprintf(fid,'\tVp1 :%9.2f+i(%-9.2f)(V)',real(VVp1),imag(VVp1));
fprintf(fid,'\tVp2 :%9.2f+i(%-9.2f)(V)',real(VVp2),imag(VVp2));
fprintf(fid,'\n');
end
fprintf(fid,'Vs0 :%9.2f+i(%-9.2f)(V)',real(VVs0),imag(VVs0));
fprintf(fid,'\tVs1 :%9.2f+i(%-9.2f)(V)',real(VVs1),imag(VVs1));
fprintf(fid,'\tVs2 :%9.2f+i(%-9.2f)(V)',real(VVs2),imag(VVs2));
fprintf(fid,'\n');
fprintf(fid,'Vx0 :%9.2f+i(%-9.2f)(V)',real(VVx0),imag(VVx0));
fprintf(fid,'\tVx1 :%9.2f+i(%-9.2f)(V)',real(VVx1),imag(VVx1));
fprintf(fid,'\tVx2 :%9.2f+i(%-9.2f)(V)',real(VVx2),imag(VVx2));
fprintf(fid,'\n');
fprintf(fid,'Vr0 :%9.2f+i(%-9.2f)(V)',real(VVr0),imag(VVr0));
fprintf(fid,'\tVr1 :%9.2f+i(%-9.2f)(V)',real(VVr1),imag(VVr1));
fprintf(fid,'\tVr2 :%9.2f+i(%-9.2f)(V)',real(VVr2),imag(VVr2));
fprintf(fid,'\n');
fprintf(fid,'Vy0 :%9.2f+i(%-9.2f)(V)',real(VVy0),imag(VVy0));
fprintf(fid,'\tVy1 :%9.2f+i(%-9.2f)(V)',real(VVy1),imag(VVy1));
fprintf(fid,'\tVy2 :%9.2f+i(%-9.2f)(V)',real(VVy2),imag(VVy2));
fprintf(fid,'\n\n');
% ----------------------------------------------------------------------
fprintf(fid,'Im0 :%9.2f+i(%-9.2f)(A)',real(IIm0),imag(IIm0));
fprintf(fid,'\tIm1 :%9.2f+i(%-9.2f)(A)',real(IIm1),imag(IIm1));
fprintf(fid,'\tIm2 :%9.2f+i(%-9.2f)(A)',real(IIm2),imag(IIm2));
fprintf(fid,'\n');
if TipoFalta==2
fprintf(fid,'Ip0 :%9.2f+i(%-9.2f)(A)',real(IIp0),imag(IIp0));
fprintf(fid,'\tIp1 :%9.2f+i(%-9.2f)(A)',real(IIp1),imag(IIp1));
fprintf(fid,'\tIp2 :%9.2f+i(%-9.2f)(A)',real(IIp2),imag(IIp2));
fprintf(fid,'\n');
end
fprintf(fid,'Is0 :%9.2f+i(%-9.2f)(A)',real(IIs0),imag(IIs0));
fprintf(fid,'\tIs1 :%9.2f+i(%-9.2f)(A)',real(IIs1),imag(IIs1));
fprintf(fid,'\tIs2 :%9.2f+i(%-9.2f)(A)',real(IIs2),imag(IIs2));
fprintf(fid,'\n');
fprintf(fid,'Ix0 :%9.2f+i(%-9.2f)(A)',real(IIx0),imag(IIx0));
fprintf(fid,'\tIx1 :%9.2f+i(%-9.2f)(A)',real(IIx1),imag(IIx1));
fprintf(fid,'\tIx2 :%9.2f+i(%-9.2f)(A)',real(IIx2),imag(IIx2));
fprintf(fid,'\n');
fprintf(fid,'Ir0 :%9.2f+i(%-9.2f)(A)',real(IIr0),imag(IIr0));
fprintf(fid,'\tIr1 :%9.2f+i(%-9.2f)(A)',real(IIr1),imag(IIr1));
fprintf(fid,'\tIr2 :%9.2f+i(%-9.2f)(A)',real(IIr2),imag(IIr2));
fprintf(fid,'\n');
fprintf(fid,'Iy0 :%9.2f+i(%-9.2f)(A)',real(IIy0),imag(IIy0));
fprintf(fid,'\tIy1 :%9.2f+i(%-9.2f)(A)',real(IIy1),imag(IIy1));
fprintf(fid,'\tIy2 :%9.2f+i(%-9.2f)(A)',real(IIy2),imag(IIy2));
fprintf(fid,'\n');
fprintf(fid,'Isx0:%9.2f+i(%-9.2f)(A)',real(IIsx0),imag(IIsx0));
fprintf(fid,'\tIsx1:%9.2f+i(%-9.2f)(A)',real(IIsx1),imag(IIsx1));
fprintf(fid,'\tIsx2:%9.2f+i(%-9.2f)(A)',real(IIsx2),imag(IIsx2));
fprintf(fid,'\n');
fprintf(fid,'Isy0:%9.2f+i(%-9.2f)(A)',real(IIsy0),imag(IIsy0));
fprintf(fid,'\tIsy1:%9.2f+i(%-9.2f)(A)',real(IIsy1),imag(IIsy1));
fprintf(fid,'\tIsy2:%9.2f+i(%-9.2f)(A)',real(IIsy2),imag(IIsy2));
fprintf(fid,'\n');
fprintf(fid,'Ixs0:%9.2f+i(%-9.2f)(A)',real(IIxs0),imag(IIxs0));
fprintf(fid,'\tIxs1:%9.2f+i(%-9.2f)(A)',real(IIxs1),imag(IIxs1));
fprintf(fid,'\tIxs2:%9.2f+i(%-9.2f)(A)',real(IIxs2),imag(IIxs2));
fprintf(fid,'\n');
fprintf(fid,'Ixr0:%9.2f+i(%-9.2f)(A)',real(IIxr0),imag(IIxr0));
fprintf(fid,'\tIxr1:%9.2f+i(%-9.2f)(A)',real(IIxr1),imag(IIxr1));
fprintf(fid,'\tIxr2:%9.2f+i(%-9.2f)(A)',real(IIxr2),imag(IIxr2));
fprintf(fid,'\n');
fprintf(fid,'Irx0:%9.2f+i(%-9.2f)(A)',real(IIrx0),imag(IIrx0));
fprintf(fid,'\tIrx1:%9.2f+i(%-9.2f)(A)',real(IIrx1),imag(IIrx1));
fprintf(fid,'\tIrx2:%9.2f+i(%-9.2f)(A)',real(IIrx2),imag(IIrx2));
fprintf(fid,'\n');
fprintf(fid,'Iry0:%9.2f+i(%-9.2f)(A)',real(IIry0),imag(IIry0));
fprintf(fid,'\tIry1:%9.2f+i(%-9.2f)(A)',real(IIry1),imag(IIry1));
fprintf(fid,'\tIry2:%9.2f+i(%-9.2f)(A)',real(IIry2),imag(IIry2));
fprintf(fid,'\n');
fprintf(fid,'Iyr0:%9.2f+i(%-9.2f)(A)',real(IIyr0),imag(IIyr0));
fprintf(fid,'\tIyr1:%9.2f+i(%-9.2f)(A)',real(IIyr1),imag(IIyr1));
fprintf(fid,'\tIyr2:%9.2f+i(%-9.2f)(A)',real(IIyr2),imag(IIyr2));
fprintf(fid,'\n');
fprintf(fid,'Iys0:%9.2f+i(%-9.2f)(A)',real(IIys0),imag(IIys0));
fprintf(fid,'\tIys1:%9.2f+i(%-9.2f)(A)',real(IIys1),imag(IIys1));
fprintf(fid,'\tIys2:%9.2f+i(%-9.2f)(A)',real(IIys2),imag(IIys2));
fprintf(fid,'\n');
fprintf(fid,'\n');
% ***********************************************************************************
% ------------------------------------------------------------------
% ** C3 ******* SALIDA V(a,b,c),I(a,b,c) Volt,Amp cartesiana *******
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vx,Vr,Vy
% Im,Ip_2,Is,Ix,Ir,Iy,Isx,Isy,Ixs,Ixr,Irx,Iry,Ir,Iys
% ------------------------------------------------------------------
% **********************************************************************
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 ******************
% ------------------------------------------------------------------
% Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys
% ------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA SALIDA 21, Ohm , polar\n');
fprintf(fid,'\n');
[mod,arg]=mod_argz(Sx_ZAB*zbase);fprintf(fid,'Sx_ZAB=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Sx_ZBC*zbase);fprintf(fid,'Sx_ZBC=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Sx_ZCA*zbase);fprintf(fid,'Sx_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sx_ZAG*zbase);fprintf(fid,'Sx_ZAG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Sx_ZBG*zbase);fprintf(fid,'Sx_ZBG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Sx_ZCG*zbase);fprintf(fid,'Sx_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_argz(Sy_ZAB*zbase);fprintf(fid,'Sy_ZAB=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Sy_ZBC*zbase);fprintf(fid,'Sy_ZBC=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Sy_ZCA*zbase);fprintf(fid,'Sy_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Sy_ZAG*zbase);fprintf(fid,'Sy_ZAG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Sy_ZBG*zbase);fprintf(fid,'Sy_ZBG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Sy_ZCG*zbase);fprintf(fid,'Sy_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_argz(Xs_ZAB*zbase);fprintf(fid,'Xs_ZAB=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Xs_ZBC*zbase);fprintf(fid,'Xs_ZBC=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Xs_ZCA*zbase);fprintf(fid,'Xs_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Xs_ZAG*zbase);fprintf(fid,'Xs_ZAG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Xs_ZBG*zbase);fprintf(fid,'Xs_ZBG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Xs_ZCG*zbase);fprintf(fid,'Xs_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_argz(Xr_ZAB*zbase);fprintf(fid,'Xr_ZAB=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Xr_ZBC*zbase);fprintf(fid,'Xr_ZBC=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Xr_ZCA*zbase);fprintf(fid,'Xr_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Xr_ZAG*zbase);fprintf(fid,'Xr_ZAG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Xr_ZBG*zbase);fprintf(fid,'Xr_ZBG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Xr_ZCG*zbase);fprintf(fid,'Xr_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_argz(Rx_ZAB*zbase);fprintf(fid,'Rx_ZAB=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Rx_ZBC*zbase);fprintf(fid,'Rx_ZBC=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Rx_ZCA*zbase);fprintf(fid,'Rx_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Rx_ZAG*zbase);fprintf(fid,'Rx_ZAG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Rx_ZBG*zbase);fprintf(fid,'Rx_ZBG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Rx_ZCG*zbase);fprintf(fid,'Rx_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_argz(Ry_ZAB*zbase);fprintf(fid,'Ry_ZAB=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Ry_ZBC*zbase);fprintf(fid,'Ry_ZBC=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Ry_ZCA*zbase);fprintf(fid,'Ry_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Ry_ZAG*zbase);fprintf(fid,'Ry_ZAG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Ry_ZBG*zbase);fprintf(fid,'Ry_ZBG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Ry_ZCG*zbase);fprintf(fid,'Ry_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_argz(Yr_ZAB*zbase);fprintf(fid,'Yr_ZAB=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Yr_ZBC*zbase);fprintf(fid,'Yr_ZBC=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Yr_ZCA*zbase);fprintf(fid,'Yr_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Yr_ZAG*zbase);fprintf(fid,'Yr_ZAG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Yr_ZBG*zbase);fprintf(fid,'Yr_ZBG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Yr_ZCG*zbase);fprintf(fid,'Yr_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_argz(Ys_ZAB*zbase);fprintf(fid,'Ys_ZAB=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Ys_ZBC*zbase);fprintf(fid,'Ys_ZBC=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Ys_ZCA*zbase);fprintf(fid,'Ys_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
[mod,arg]=mod_argz(Ys_ZAG*zbase);fprintf(fid,'Ys_ZAG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Ys_ZBG*zbase);fprintf(fid,'Ys_ZBG=%7.2f(%-5.1f)(Ohm)\t',mod,arg);
[mod,arg]=mod_argz(Ys_ZCG*zbase);fprintf(fid,'Ys_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
% -----------------------------------------------------------------------
% ** C4 ****************** SALIDA 67N, Volt, Amp, Polar *****************
% -----------------------------------------------------------------------
% Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys
% -----------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'SALIDA SALIDA 67N, Volt, Amp, Polar\n');
fprintf(fid,'\n');
fprintf(fid,'Sx_67N:\tVspol=%-8.0f(V)\tIsxop=%-8.0f(A)\tAngDirSx=%-5.1f\n',VVspol,IIsxop,AngDirSx);
fprintf(fid,'Sy_67N:\tVspol=%-8.0f(V)\tIsyop=%-8.0f(A)\tAngDirSy=%-5.1f\n',VVspol,IIsyop,AngDirSy);
fprintf(fid,'Xs_67N:\tVxpol=%-8.0f(V)\tIxsop=%-8.0f(A)\tAngDirXs=%-5.1f\n',VVxpol,IIxsop,AngDirXs);
fprintf(fid,'Xr_67N:\tVxpol=%-8.0f(V)\tIxrop=%-8.0f(A)\tAngDirXr=%-5.1f\n',VVxpol,IIxrop,AngDirXr);
fprintf(fid,'Rx_67N:\tVrpol=%-8.0f(V)\tIrxop=%-8.0f(A)\tAngDirRx=%-5.1f\n',VVrpol,IIrxop,AngDirRx);
fprintf(fid,'Ry_67N:\tVrpol=%-8.0f(V)\tIryop=%-8.0f(A)\tAngDirRy=%-5.1f\n',VVrpol,IIryop,AngDirRy);
fprintf(fid,'Yr_67N:\tVypol=%-8.0f(V)\tIyrop=%-8.0f(A)\tAngDirYr=%-5.1f\n',VVypol,IIyrop,AngDirYr);
fprintf(fid,'Ys_67N:\tVypol=%-8.0f(V)\tIysop=%-8.0f(A)\tAngDirYs=%-5.1f\n',VVypol,IIysop,AngDirYs);
fprintf(fid,'\n');
% ------------------------------------------------------------------
% ** C4 ********* SALIDA V(0,1,2), I(0,1,2) Volt,Amp polar ********
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vx,Vr,Vy
% Im,Ip_2,Is,Ix,Ir,Iy,Isx,Isy,Ixs,Ixr,Irx,Iry,Ir,Iys
% ------------------------------------------------------------------
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\t',mod,arg);
[mod,arg]=mod_arg(VVm1);fprintf(fid,'Vm1 =%8.0f(%-5.1f)(V)\t\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\t',mod,arg);
[mod,arg]=mod_arg(VVp1);fprintf(fid,'Vp1 =%8.0f(%-5.1f)(V)\t\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\t',mod,arg);
[mod,arg]=mod_arg(VVs1);fprintf(fid,'Vs1 =%8.0f(%-5.1f)(V)\t\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(VVx0);fprintf(fid,'Vx0 =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVx1);fprintf(fid,'Vx1 =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVx2);fprintf(fid,'Vx2 =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVr0);fprintf(fid,'Vr0 =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVr1);fprintf(fid,'Vr1 =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVr2);fprintf(fid,'Vr2 =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVy0);fprintf(fid,'Vy0 =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVy1);fprintf(fid,'Vy1 =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVy2);fprintf(fid,'Vy2 =%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\t',mod,arg);
[mod,arg]=mod_arg(IIm1);fprintf(fid,'Im1 =%8.0f(%-5.1f)(A)\t\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\t',mod,arg);
[mod,arg]=mod_arg(IIp1);fprintf(fid,'Ip1 =%8.0f(%-5.1f)(A)\t\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\t',mod,arg);
[mod,arg]=mod_arg(IIs1);fprintf(fid,'Is1 =%8.0f(%-5.1f)(A)\t\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(IIx0);fprintf(fid,'Ix0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIx1);fprintf(fid,'Ix1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIx2);fprintf(fid,'Ix2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIr0);fprintf(fid,'Ir0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIr1);fprintf(fid,'Ir1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIr2);fprintf(fid,'Ir2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIy0);fprintf(fid,'Iy0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIy1);fprintf(fid,'Iy1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIy2);fprintf(fid,'Iy2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIsx0);fprintf(fid,'Isx0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIsx1);fprintf(fid,'Isx1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIsx2);fprintf(fid,'Isx2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIsy0);fprintf(fid,'Isy0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIsy1);fprintf(fid,'Isy1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIsy2);fprintf(fid,'Isy2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIxs0);fprintf(fid,'Ixs0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxs1);fprintf(fid,'Ixs1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxs2);fprintf(fid,'Ixs2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIxr0);fprintf(fid,'Ixr0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxr1);fprintf(fid,'Ixr1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxr2);fprintf(fid,'Ixr2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIrx0);fprintf(fid,'Irx0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIrx1);fprintf(fid,'Irx1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIrx2);fprintf(fid,'Irx2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIry0);fprintf(fid,'Iry0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIry1);fprintf(fid,'Iry1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIry2);fprintf(fid,'Iry2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIyr0);fprintf(fid,'Iyr0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIyr1);fprintf(fid,'Iyr1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIyr2);fprintf(fid,'Iyr2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIys0);fprintf(fid,'Iys0 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIys1);fprintf(fid,'Iys1 =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIys2);fprintf(fid,'Iys2 =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
% ------------------------------------------------------------------
% ** C4 ******** SALIDA V(a,b,c), I(a,b,c) Volt,Amp polar *********
% ------------------------------------------------------------------
% Vm,Vp_2,Vs,Vx,Vr,Vy
% Im,Ip_2,Is,Ix,Ir,Iy,Isx,Isy,Ixs,Ixr,Irx,Iry,Ir,Iys
% ------------------------------------------------------------------
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\t',mod,arg);
[mod,arg]=mod_arg(VVmb);fprintf(fid,'Vmb =%8.0f(%-5.1f)(V)\t\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\t',mod,arg);
[mod,arg]=mod_arg(VVpb);fprintf(fid,'Vpb =%8.0f(%-5.1f)(V)\t\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\t',mod,arg);
[mod,arg]=mod_arg(VVsb);fprintf(fid,'Vsb =%8.0f(%-5.1f)(V)\t\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(VVxa);fprintf(fid,'Vxa =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVxb);fprintf(fid,'Vxb =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVxc);fprintf(fid,'Vxc =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVra);fprintf(fid,'Vra =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVrb);fprintf(fid,'Vrb =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVrc);fprintf(fid,'Vrc =%8.0f(%-5.1f)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(VVya);fprintf(fid,'Vya =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVyb);fprintf(fid,'Vyb =%8.0f(%-5.1f)(V)\t\t',mod,arg);
[mod,arg]=mod_arg(VVyc);fprintf(fid,'Vyc =%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\t',mod,arg);
[mod,arg]=mod_arg(IImb);fprintf(fid,'Imb =%8.0f(%-5.1f)(A)\t\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\t',mod,arg);
[mod,arg]=mod_arg(IIpb);fprintf(fid,'Ipb =%8.0f(%-5.1f)(A)\t\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\t',mod,arg);
[mod,arg]=mod_arg(IIsb);fprintf(fid,'Isb =%8.0f(%-5.1f)(A)\t\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(IIxa);fprintf(fid,'Ixa =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxb);fprintf(fid,'Ixb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxc);fprintf(fid,'Ixc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIra);fprintf(fid,'Ira =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIrb);fprintf(fid,'Irb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIrc);fprintf(fid,'Irc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIya);fprintf(fid,'Iya =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIyb);fprintf(fid,'Iyb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIyc);fprintf(fid,'Iyc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIsxa);fprintf(fid,'Isxa =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIsxb);fprintf(fid,'Isxb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIsxc);fprintf(fid,'Isxc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIsya);fprintf(fid,'Isya =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIsyb);fprintf(fid,'Isyb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIsyc);fprintf(fid,'Isyc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIxsa);fprintf(fid,'Ixsa =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxsb);fprintf(fid,'Ixsb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxsc);fprintf(fid,'Ixsc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIxra);fprintf(fid,'Ixra =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxrb);fprintf(fid,'Ixrb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIxrc);fprintf(fid,'Ixrc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIrxa);fprintf(fid,'Irxa =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIrxb);fprintf(fid,'Irxb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIrxc);fprintf(fid,'Irxc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIrya);fprintf(fid,'Irya =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIryb);fprintf(fid,'Iryb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIryc);fprintf(fid,'Iryc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIyra);fprintf(fid,'Iyra =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIyrb);fprintf(fid,'Iyrb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIyrc);fprintf(fid,'Iyrc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(IIysa);fprintf(fid,'Iysa =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIysb);fprintf(fid,'Iysb =%8.0f(%-5.1f)(A)\t\t',mod,arg);
[mod,arg]=mod_arg(IIysc);fprintf(fid,'Iysc =%8.0f(%-5.1f)(A)',mod,arg);
fprintf(fid,'\n');
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.
