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,'PC162\n');
fprintf(fid,'Esquema Linea Doble SR\n');
fprintf(fid,'Una de las Lineas abierta y puesta a Tierra en ambos extremos\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,'Impedancia de Transferencia: Si,No\n');
fprintf(fid,'Compensacion 21: Si,No\n');
fprintf(fid,'Eleccion: kE,L,Sccs,Sccr,Kr0,k0=zs0/zr0\n');
fprintf(fid,'Falta Simple\n');
fprintf(fid,'\n');
% -------------------------------------------------------------------
%PC162
datos=input('datos: standard=1,input=2: ');
if datos==1 
Un=220;
kE=1;
L=220;
imptransfer=2;
compens21=2;
Sccs=5000;
Sccr=5000;
Kr0=1;
% -------------------------------------------------------------------
else
Un=input('U(Fase-Fase) kV:400,220,132,66: ');
kE=input('kE=E/Enominal: ');
L=input('Longitud Linea km:  ');
imptransfer=input('imptransfer: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 :');
Sccs=input('Sccs(MVA): ');
Sccr=input('Sccr(MVA): ');
Kr0=input('Kr0: ');
end
%-------------------------------------------------------------------
k0=input('k0:  ');
m=input('m:  ');
%-------------------------------------------------------------------
pfaltam=1;
%-------------------------------------------------------------------
En=Un/sqrt(3);
%U=kE*Un; 
E=kE*En;
if Un==400
   zl1u=0.0346+0.3208i; zl2u=zl1u; zl0u=0.1828+1.1482i; zlM0u=0.1482+0.7383i;
elseif Un==220
   zl1u=0.0829+0.4135i; zl2u=zl1u; zl0u=0.2311+1.3085i; zlM0u=0.1482+0.8031i;
elseif Un==132
   zl1u=0.1342+0.4042i; zl2u=zl1u; zl0u=0.2824+1.3718i; zlM0u=0.1482+0.8947i;
elseif Un==66
   zl1u=0.2156+0.3952i; zl2u=zl1u; zl0u=0.3636+1.4278i; zlM0u=0.1482+0.9288i; 
end
k=(zl0u-zl1u)/(3*zl1u);
zl1=zl1u*L; zl2=zl2u*L; zl0=zl0u*L; zlM0=zlM0u*L;
zs1=(Un^2/Sccs)*i;
zs2=zs1;
zr1=(Un^2/Sccr)*i;
zr2=zr1;
zr0=Kr0*zr1;
zs0=k0*zr0;
Km=zlM0u/zl1u;
Sm=zl1u*m*L;
Rm=zl1u*(1-m)*L;

% CALCULOS
% -------------------------------------------------------------------
switch(imptransfer)
case{1},
[Minc]=subMincT(1);
[Minca]=subMincaT(1);
[z0,P0,Q0]=subz0T(zs0,zr0,zsr0,zl0,zlM0,m);
[z1]=subz1T(zs1,zr1,zsr1,zl1,m);
[z2]=subz2T(zs2,zr2,zsr2,zl2,m);
case(2)
[Minc]=subMinc(1);
[Minca]=subMinca(1);
[z0,P0,Q0]=subz0(zs0,zr0,zl0,zlM0,m);
[z1]=subz1(zs1,zr1,zl1,m);
[z2]=subz2(zs2,zr2,zl2,m);
end 

% MATRIZ DE ADMITANCIA DE BARRAS
[Y0]=subY0(Minc,z0);
[Y1]=subY1(Minc,z1);
[Y2]=subY2(Minc,z2);

% ZBUS(0,1,2)
Z0=inv(Y0);
Z1=inv(Y1);
Z2=inv(Y2);

zbase=1;
% V(0,1,2), I(0,1,2) EN EL PUNTO m DE LA FALTA
[Im0,Im1,Im2,Vm0,Vm1,Vm2,faltam,Rfm,Rgm]=ma2m(Z0,Z1,Z2,E,zbase,pfaltam);

% IBUS(0,1,2)
I0=[-Im0;0;0];
I1=[-Im1;0;0];
I2=[-Im2;0;0];

% 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);

% ASIGNACIONES V(0,1,2)
Vs0=V0(2); Vs1=V1(2); Vs2=V2(2);
Vr0=V0(3); Vr1=V1(3); Vr2=V2(3);

% ASIGNACIONES I(0,1,2)
Is0=Ie0(1); Is1=Ie1(1); Is2=Ie2(1);
Ir0=Ie0(2); Ir1=Ie1(2); Ir2=Ie2(2);

Ism0=Ie0(3); Ism1=Ie1(3); Ism2=Ie2(3);
Irm0=Ie0(4); Irm1=Ie1(4); Irm2=Ie2(4);
if imptransfer==1
Isrm0=Ie0(5); Isrm1=Ie1(5); Isrm2=Ie2(5);
end

% CALCULO DE LA CORRIENTE EN LA LINEA PUESTA A TIERRA EN AMBOS EXTREMOS
%[It0]=subIt0(Ism0,Irm0,zl0,zlM0,m);
It0=-inv(Q0)*P0*Ie0;

% V(a,b,c), I(a,b,c) EN EL PUNTO m DE LA FALTA
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);

% V(a,b,c), I(a,b,c) EN LOS EXTREMOS DE LINEA
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[Isa,Isb,Isc]=abc(Is0,Is1,Is2);
[Ira,Irb,Irc]=abc(Ir0,Ir1,Ir2);
[Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2);
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2);
if imptransfer==1
[Isrma,Isrmb,Isrmc]=abc(Isrm0,Isrm1,Isrm2);
end

% 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,-It0);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
   PDistComp(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k,Km,It0);
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);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
   PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k);
end

% CALCULO DE LAS MAGNITUDES DE ALIMETACION DE LA PROTECCION
% DIRECCIONAL DE TIERRA EN LOS EXTREMOS DE LINEA
[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0);
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0);
%
% PREPARACION SALIDAS
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
%-----------------------------------------------------
Vm0=Vm0*1000; Vm1=Vm1*1000; Vm2=Vm2*1000;
Vs0=Vs0*1000; Vs1=Vs1*1000; Vs2=Vs2*1000;
Vr0=Vr0*1000; Vr1=Vr1*1000; Vr2=Vr2*1000;

Im0=Im0*1000; Im1=Im1*1000; Im2=Im2*1000;
Is0=Is0*1000; Is1=Is1*1000; Is2=Is2*1000;
Ir0=Ir0*1000; Ir1=Ir1*1000; Ir2=Ir2*1000;

Ism0=Ism0*1000; Ism1=Ism1*1000; Ism2=Ism2*1000;
Irm0=Irm0*1000; Irm1=Irm1*1000; Irm2=Irm2*1000;
if imptransfer==1
Isrm0=Isrm0*1000; Isrm1=Isrm1*1000; Isrm2=Isrm2*1000;
end

It0=It0*1000;

[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);

[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[Isa,Isb,Isc]=abc(Is0,Is1,Is2);
[Ira,Irb,Irc]=abc(Ir0,Ir1,Ir2);
[Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2);
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2);
if imptransfer==1
[Isrma,Isrmb,Isrmc]=abc(Isrm0,Isrm1,Isrm2);
end

[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0);
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0);


% -------------------------------------------------------------------
% SALIDAS
% -------------------------------------------------------------------
fprintf(fid,'\n');
fprintf(fid,'***************************************************\n');
fprintf(fid,'Un=%-4.0f(kV)\n',Un);
fprintf(fid,'kE=E/En=%-4.0f\n',kE);
fprintf(fid,'L=%-5.1f(km)\n',L);
fprintf(fid,'Compensacion 21, Si(1), No(2) =%-1.0f\n',compens21);
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=zs0/zr0=%-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
fprintf(fid,'\n');
fprintf(fid,'***************************************************\n');
variassalidas=1;
while variassalidas==1
fprintf(fid,'\n');
output=input('Salida: Ohm-Volt-Amp(cart)=1, Ohm-Volt-Amp(polar)=2 :');
if salida==2
fprintf(fid,'Salida: Ohm-Volt-Amp(cart)=1, Ohm-Volt-Amp(polar)=2 : %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:OVA(cartesiana)   ***          ***        ***
   % 2:OVA(polar)        ***  ***     ***        ***
   
% **********************************************************************
% ***********************************************************************
case{1}, % 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 
% -----------------------------------------------------------------------
% -----------------------------------------------------------------------
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');
fprintf(fid,'Sm =%-7.2f+i(%-7.2f)Ohm\n',real(Sm),imag(Sm));
fprintf(fid,'Rm =%-7.2f+i(%-7.2f)Ohm\n',real(Rm),imag(Rm));
fprintf(fid,'\n');
[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');
[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');
   
% ------------------------------------------------------------------
%   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(Vm0),imag(Vm0));
fprintf(fid,'Vm1  =%8.0f+i(%-8.0f)V\t',real(Vm1),imag(Vm1));
fprintf(fid,'Vm2  =%8.0f+i(%-8.0f)V\t',real(Vm2),imag(Vm2));
fprintf(fid,'\n');
fprintf(fid,'Vs0  =%8.0f+i(%-8.0f)V\t',real(Vs0),imag(Vs0));
fprintf(fid,'Vs1  =%8.0f+i(%-8.0f)V\t',real(Vs1),imag(Vs1));
fprintf(fid,'Vs2  =%8.0f+i(%-8.0f)V\t',real(Vs2),imag(Vs2));
fprintf(fid,'\n');
fprintf(fid,'Vr0  =%8.0f+i(%-8.0f)V\t',real(Vr0),imag(Vr0));
fprintf(fid,'Vr1  =%8.0f+i(%-8.0f)V\t',real(Vr1),imag(Vr1));
fprintf(fid,'Vr2  =%8.0f+i(%-8.0f)V\t',real(Vr2),imag(Vr2));
fprintf(fid,'\n');
% ---------------------------------------------------------------
fprintf(fid,'Im0  =%8.0f+i(%-8.0f)A\t',real(Im0),imag(Im0));
fprintf(fid,'Im1  =%8.0f+i(%-8.0f)A\t',real(Im1),imag(Im1));
fprintf(fid,'Im2  =%8.0f+i(%-8.0f)A\t',real(Im2),imag(Im2));
fprintf(fid,'\n');
fprintf(fid,'Is0  =%8.0f+i(%-8.0f)A\t',real(Is0),imag(Is0));
fprintf(fid,'Is1  =%8.0f+i(%-8.0f)A\t',real(Is1),imag(Is1));
fprintf(fid,'Is2  =%8.0f+i(%-8.0f)A\t',real(Is2),imag(Is2));
fprintf(fid,'\n');
fprintf(fid,'Ir0  =%8.0f+i(%-8.0f)A\t',real(Ir0),imag(Ir0));
fprintf(fid,'Ir1  =%8.0f+i(%-8.0f)A\t',real(Ir1),imag(Ir1));
fprintf(fid,'Ir2  =%8.0f+i(%-8.0f)A\t',real(Ir2),imag(Ir2));
fprintf(fid,'\n');
fprintf(fid,'Ism0 =%8.0f+i(%-8.0f)A\t',real(Ism0),imag(Ism0));
fprintf(fid,'Ism1 =%8.0f+i(%-8.0f)A\t',real(Ism1),imag(Ism1));
fprintf(fid,'Ism2 =%8.0f+i(%-8.0f)A\t',real(Ism2),imag(Ism2));
fprintf(fid,'\n');
fprintf(fid,'Irm0 =%8.0f+i(%-8.0f)A\t',real(Irm0),imag(Irm0));
fprintf(fid,'Irm1 =%8.0f+i(%-8.0f)A\t',real(Irm1),imag(Irm1));
fprintf(fid,'Irm2 =%8.0f+i(%-8.0f)A\t',real(Irm2),imag(Irm2));
fprintf(fid,'\n');
if imptransfer==1
fprintf(fid,'Isrm0 =%8.0f+i(%-8.0f)A\t',real(Isrm0),imag(Isrm0));
fprintf(fid,'Isrm1 =%8.0f+i(%-8.0f)A\t',real(Isrm1),imag(Isrm1));
fprintf(fid,'Isrm2 =%8.0f+i(%-8.0f)A\t',real(Isrm2),imag(Isrm2));
fprintf(fid,'\n');
end
fprintf(fid,'It0 =%8.0f+i(%-8.0f)A\t',real(It0),imag(It0));
fprintf(fid,'\n');
% ------------------------------------------------------------------
%    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(Vma),imag(Vma));
fprintf(fid,'Vmb  =%8.0f+i(%-8.0f)V\t',real(Vmb),imag(Vmb));
fprintf(fid,'Vmc  =%8.0f+i(%-8.0f)V\t',real(Vmc),imag(Vmc));
fprintf(fid,'\n');
fprintf(fid,'Vsa  =%8.0f+i(%-8.0f)V\t',real(Vsa),imag(Vsa));
fprintf(fid,'Vsb  =%8.0f+i(%-8.0f)V\t',real(Vsb),imag(Vsb));
fprintf(fid,'Vsc  =%8.0f+i(%-8.0f)V\t',real(Vsc),imag(Vsc));
fprintf(fid,'\n');
fprintf(fid,'Vra  =%8.0f+i(%-8.0f)V\t',real(Vra),imag(Vra));
fprintf(fid,'Vrb  =%8.0f+i(%-8.0f)V\t',real(Vrb),imag(Vrb));
fprintf(fid,'Vrc  =%8.0f+i(%-8.0f)V\t',real(Vrc),imag(Vrc));
fprintf(fid,'\n');
fprintf(fid,'Ima  =%8.0f+i(%-8.0f)A\t',real(Ima),imag(Ima));
fprintf(fid,'Imb  =%8.0f+i(%-8.0f)A\t',real(Imb),imag(Imb));
fprintf(fid,'Imc  =%8.0f+i(%-8.0f)A\t',real(Imc),imag(Imc));
fprintf(fid,'\n');
fprintf(fid,'Isa  =%8.0f+i(%-8.0f)A\t',real(Isa),imag(Isa));
fprintf(fid,'Isb  =%8.0f+i(%-8.0f)A\t',real(Isb),imag(Isb));
fprintf(fid,'Isc  =%8.0f+i(%-8.0f)A\t',real(Isc),imag(Isc));
fprintf(fid,'\n');
fprintf(fid,'Ira  =%8.0f+i(%-8.0f)A\t',real(Ira),imag(Ira));
fprintf(fid,'Irb  =%8.0f+i(%-8.0f)A\t',real(Irb),imag(Irb));
fprintf(fid,'Irc  =%8.0f+i(%-8.0f)A\t',real(Irc),imag(Irc));
fprintf(fid,'\n');
fprintf(fid,'Isma =%8.0f+i(%-8.0f)A\t',real(Isma),imag(Isma));
fprintf(fid,'Ismb =%8.0f+i(%-8.0f)A\t',real(Ismb),imag(Ismb));
fprintf(fid,'Ismc =%8.0f+i(%-8.0f)A\t',real(Ismc),imag(Ismc));
fprintf(fid,'\n');
fprintf(fid,'Irma =%8.0f+i(%-8.0f)A\t',real(Irma),imag(Irma));
fprintf(fid,'Irmb =%8.0f+i(%-8.0f)A\t',real(Irmb),imag(Irmb));
fprintf(fid,'Irmc =%8.0f+i(%-8.0f)A\t',real(Irmc),imag(Irmc));
fprintf(fid,'\n');
if imptransfer==1
fprintf(fid,'Isrma =%8.0f+i(%-8.0f)A\t',real(Isrma),imag(Isrma));
fprintf(fid,'Isrmb =%8.0f+i(%-8.0f)A\t',real(Isrmb),imag(Isrmb));
fprintf(fid,'Isrmc =%8.0f+i(%-8.0f)A\t',real(Isrmc),imag(Isrmc));
fprintf(fid,'\n');
end
% ******************************************************************
case{2}, % OVA polar: 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c)
%                    ---  --- 
% ******************************************************************
% ------------------------------------------------------------------
% ** C4.1 ****************** SALIDA 21, Ohm , polar ******************
% ------------------------------------------------------------------
% ------------------------------------------------------------------
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');
   [mod,arg]=mod_arg(Sm);fprintf(fid,'Sm=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
   [mod,arg]=mod_arg(Rm);fprintf(fid,'Rm=%7.2f(%-5.1fº)(Ohm)\n',mod,arg);
fprintf(fid,'\n');
[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');
[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');
% -----------------------------------------------------------------------
% ** C4.2 ****************** 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');
fprintf(fid,'Sm_67N:\tVspol=%-8.0f(V)\tIsmop=%-8.0f(A)\tAngDirSm=%-5.1fº\n',Vspol,Ismop,AngDirSm);
fprintf(fid,'Rm_67N:\tVrpol=%-8.0f(V)\tIrmop=%-8.0f(A)\tAngDirRm=%-5.1fº\n',Vrpol,Irmop,AngDirRm);
fprintf(fid,'\n');
% ------------------------------------------------------------------
% ** C4.3 ********* 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(Vm0);fprintf(fid,'Vm0  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vm1);fprintf(fid,'Vm1  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vm2);fprintf(fid,'Vm2  =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Vs0);fprintf(fid,'Vs0  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vs1);fprintf(fid,'Vs1  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vs2);fprintf(fid,'Vs2  =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Vr0);fprintf(fid,'Vr0  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vr1);fprintf(fid,'Vr1  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vr2);fprintf(fid,'Vr2  =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
fprintf(fid,'\n');
%-------------------------------------------------------------------------------
[mod,arg]=mod_arg(Im0);fprintf(fid,'Im0  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Im1);fprintf(fid,'Im1  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Im2);fprintf(fid,'Im2  =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Is0);fprintf(fid,'Is0  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Is1);fprintf(fid,'Is1  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Is2);fprintf(fid,'Is2  =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Ir0);fprintf(fid,'Ir0  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ir1);fprintf(fid,'Ir1  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ir2);fprintf(fid,'Ir2  =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Ism0);fprintf(fid,'Ism0 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ism1);fprintf(fid,'Ism1 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ism2);fprintf(fid,'Ism2 =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Irm0);fprintf(fid,'Irm0 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irm1);fprintf(fid,'Irm1 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irm2);fprintf(fid,'Irm2 =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
if imptransfer==1
[mod,arg]=mod_arg(Isrm0);fprintf(fid,'Isrm0 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isrm1);fprintf(fid,'Isrm1 =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isrm2);fprintf(fid,'Isrm2 =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
end
[mod,arg]=mod_arg(It0);fprintf(fid,'It0  =%8.0f(%-5.1fº)(A)\t',mod,arg);
fprintf(fid,'\n');

% ------------------------------------------------------------------
% ** C4.4 ******** 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(Vma);fprintf(fid,'Vma  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vmb);fprintf(fid,'Vmb  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vmc);fprintf(fid,'Vmc  =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Vsa);fprintf(fid,'Vsa  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vsb);fprintf(fid,'Vsb  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vsc);fprintf(fid,'Vsc  =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Vra);fprintf(fid,'Vra  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vrb);fprintf(fid,'Vrb  =%8.0f(%-5.1fº)(V)\t',mod,arg);
[mod,arg]=mod_arg(Vrc);fprintf(fid,'Vrc  =%8.0f(%-5.1fº)(V)',mod,arg);
fprintf(fid,'\n');
%------------------------------------------------------------------------
fprintf(fid,'\n');
%------------------------------------------------------------------------
[mod,arg]=mod_arg(Ima);fprintf(fid,'Ima  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Imb);fprintf(fid,'Imb  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Imc);fprintf(fid,'Imc  =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Isa);fprintf(fid,'Isa  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isb);fprintf(fid,'Isb  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isc);fprintf(fid,'Isc  =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Ira);fprintf(fid,'Ira  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irb);fprintf(fid,'Irb  =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irc);fprintf(fid,'Irc  =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Isma);fprintf(fid,'Isma =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ismb);fprintf(fid,'Ismb =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Ismc);fprintf(fid,'Ismc =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
[mod,arg]=mod_arg(Irma);fprintf(fid,'Irma =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irmb);fprintf(fid,'Irmb =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Irmc);fprintf(fid,'Irmc =%8.0f(%-5.1fº)(A)',mod,arg);
fprintf(fid,'\n');
if imptransfer==1
[mod,arg]=mod_arg(Isrma);fprintf(fid,'Isrma =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isrmb);fprintf(fid,'Isrmb =%8.0f(%-5.1fº)(A)\t',mod,arg);
[mod,arg]=mod_arg(Isrmc);fprintf(fid,'Isrmc =%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');
<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd">
<html xmlns="http://www.w3.org/1999/xhtml" xml:lang="es">
	<head>
            <title>e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.</title>

                <meta http-equiv="Content-Type" content="text/html; charset=utf-8"/>
<meta name="description" content="Buscador de proyectos ETSI" />
<meta name="author" content="Daniel Callejo"/>
<meta name="author" content="Federico L&oacute;pez"/>
<meta name="author" content="Juli&aacute;n P&eacute;rez"/>
<meta name="author" content="Jos&eacute; Mar&iacute;a Vidal Vidal"/>
<base href="https://biblus.us.es/bibing/proyectos/" /><link rel="shortcut icon" href="./favicon.ico" />
<link rel="stylesheet" href="./publico/css/estructura.css" type="text/css" />
<link rel="stylesheet" href="./publico/css/bootstrap.css" type="text/css" media="screen"/>
<link rel="stylesheet" href="./publico/css/enlaces-capas.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/pestanas.css" type="text/css"  /> 
<link rel="stylesheet" href="./publico/css/formularios.css" type="text/css" />
<link rel="stylesheet" href="./publico/css/lightbox.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/autocompleter.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/resultados.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/milista.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/explorador.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/ultimos.css" type="text/css"  />
<link rel="stylesheet" href="./publico/css/contactar.css" type="text/css"  />
<!-- <link rel="stylesheet" href="./css/print.css" type="text/css" media="print" /> -->
<script src="./publico/javascript/jquery-1.9.1.min.js" type="text/javascript"></script>
<script src="./publico/javascript/bootstrap.js" type="text/javascript"></script>

<script src="./publico/javascript/mootools-1.2.4-core.js" type="text/javascript"></script>
<script src="./publico/javascript/mootools-1-1.2.4.2-more.js" type="text/javascript"></script>
<script src="./publico/javascript/targetbl.js" type="text/javascript"></script>
<script src="./publico/javascript/mensajes.js" type="text/javascript"></script>
<script src="./publico/javascript/form.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Observer.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Autocompleter.js" type="text/javascript"></script>
<script src="./publico/javascript/autocompleter/Autocompleter.Request.js" type="text/javascript"></script>
<script src="./publico/javascript/domready.js" type="text/javascript"></script>

<!-- Google tag (gtag.js) -->
<script async src="https://www.googletagmanager.com/gtag/js?id=G-2HP38G406P"></script>
<script>
  window.dataLayer = window.dataLayer || [];
  function gtag(){dataLayer.push(arguments);}
  gtag('js', new Date());

  gtag('config', 'G-2HP38G406P');
</script>



<script type="text/javascript">
function seleccionar_todo(){
   for (i=0;i<document.formulario_marcados.elements.length;i++)
      if(document.formulario_marcados.elements[i].type == "checkbox")
         document.formulario_marcados.elements[i].checked=1
}
function deseleccionar_todo(){
   for (i=0;i<document.formulario_marcados.elements.length;i++)
      if(document.formulario_marcados.elements[i].type == "checkbox")
         document.formulario_marcados.elements[i].checked=0
} 
</script>
	</head>
    <body style="min-width: 960px; background-color: #e2e2e2; height: 100%">	
		<a name="arriba"></a>
<div>
<div id="header">
    <table id="Tabla_01" height="105" border="0" cellpadding="0" cellspacing="0">
    <tr>
        <td>
            <img src="./publico/imagenes/cabecera/banner_new2_01.jpg" width="3" height="105" alt="">
        </td>
        <td>
            <img src="./publico/imagenes/cabecera/banner_new2_02.jpg" width="100" height="105" alt="">
        </td>
        <td width="100%" align="right"  style="min-width: 613px; background-repeat: no-repeat; background-position: top right; font-size: medium; font-weight: bold; color: white;">
            <table width="100%" height="105" border="0">
                <tr>
                    <td style=" vertical-align: top; width: 547px">
                        <!-- <a style="text-decoration: none" href="./"><img src="./publico/imagenes/cabecera/banner_new2_03.jpg"alt=""></a> -->
                        <a style="text-decoration: none" href="https://biblus.us.es/bibing/proyectos/"><img src="./publico/imagenes/cabecera/banner_new2_03.jpg"alt=""></a>
                     </td>
                     <td  background="./publico/imagenes/cabecera/banner_new2_04.jpg" style="background-repeat: no-repeat; background-color: #c6d4e8;"><br><br><br>
                         &nbsp;
                        </td>
                    </tr>
             </table>
        </td>
    </tr>
    <tr style="border-top: 1px solid #ffffff;">
        <td colspan="3" id="menu_sup" style="font-size: small; text-align: right; background-color: #9c1f2f; vertical-align: bottom">
                                     

                                <a href="https://biblus.us.es/bibing/proyectos/" class="enlace " style=" color: #ffffff">Inicio&nbsp;&nbsp;|</a>
                                <a href="https://biblus.us.es/bibing/proyectos/acerca/" style=" color: #ffffff" class="enlace ">&nbsp;&nbsp;Acerca de&nbsp;&nbsp;|</a>
                                
                                <a href="https://biblus.us.es/bibing/proyectos/contactar/" style=" color: #ffffff" class="enlace ">&nbsp;&nbsp;Contactar&nbsp;&nbsp;</a>
                         <!--       <a href="https://biblus.us.es/bibing/proyectos/ayuda/" style=" color: #ffffff" class="enlace  ">&nbsp;&nbsp;Ayuda&nbsp;&nbsp;|</a> -->
                                <!-- <form id="idioma" class=""  method="post" action="">
                                <div class="inline" style="color: white">
                                <input type="submit" class="tipo_letra " style=" color: #ffffff"name="idioma" value="Castellanossss" />&nbsp;&nbsp;
                                </div>	
                                </form> -->

                                </td>

    </tr>
    </table>
</div>
<br><br>
<div class="colmask rightmenu">
<div class="panel panel-default" style="margin:20px">
            <div class="panel-heading"><h4>AN�LISIS DE ACTUACIÓN DE LAS PROTECCIONES 21 Y 67N ANTE FALTAS SIMPLES, MULTIPLES E INTERCIRCUITO EN L�NEAS ELÉCTRICAS SOBRE LOS MISMOS APOYOS</h4><h5 style="color:grey"><span class="glyphicon glyphicon-user""></span> : <b>Ruiz Lozano, José Luis</b><br><span class="glyphicon glyphicon-book""></span>  :<b> Ingenier&iacute;a Industrial</b></h5></div><div class="panel-body"><div class="contenido_proyecto"><div class="titulo2_explorador">Contenido del proyecto:</div><div class="ruta"><a href="use/abreproy/4273">Directorio ra&iacute;z</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252F">z_Prog</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPC%252F">PC</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPC%252FPC162%252F">PC162</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPC%252FPC162%252FFC162.M%252F">FC162.M</a></div><div class="listing"></div></div></div></div>
</div>
</div>

        
	</body>
 
 
<!-- Inicio Google Analytics -->
<!--
	<script type="text/javascript">
		var gaJsHost = (("https:" == document.location.protocol) ? "https://ssl." : "http://www.");
		document.write(unescape("%3Cscript src='" + gaJsHost + "google-analytics.com/ga.js' type='text/javascript'%3E%3C/script%3E"));
	</script>
	
	<script type="text/javascript">
		try {
		var pageTracker = _gat._getTracker("UA-7789620-1");
		pageTracker._trackPageview();
	} 
		catch(err) {}
	</script>
-->
<!-- Fin Google Analytics -->

</html>