% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG124mlc\n');
fprintf('Esquema Linea Doble SR\n');
fprintf('Unominal kV: 400,220,132,66\n');
fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf('Impedancia de Transferencia: Si,No\n');
fprintf('Eleccion: (zsr1,Ksr0)\n');
fprintf('Eleccion: kE,L,Sccs,Sccr,Kr0\n');
fprintf('Eleccion: k0=zs0/zr0\n');
fprintf('Compensacion 21, Si, No\n');
fprintf('Falta Doble en la misma Linea, m=0...1, p=fijo a elegir\n');
%fprintf('Salida Grafica 1: Sm,Rm, Unidad Preferente (pu)\n');
%fprintf('Salida Grafica 2: Sm,Rm, Unidad Preferente (pua)\n');
%fprintf('Salida Grafica 3: Sm, Todas las Unidades (pua)\n');
%fprintf('                : Sr, Todas las Unidades (pua)\n');
%fprintf('                : Rm, Todas las Unidades (pua)\n');
%fprintf('                : Rs, Todas las Unidades (pua)\n');
%fprintf('Salida Grafica 4: S(Linea I), 67N\n');
%fprintf('                : R(Linea I), 67N\n');
%fprintf('                : S(Linea II), 67N\n');
%fprintf('                : R(Linea II), 67N\n');
fprintf('\n');
% -------------------------------------------------------------------
%PG124mlc
datos=input('datos: standard=1,input=2: ');
if datos==1 
Un=220;
kE=1;
L=220;
Sccs=5000;
Sccr=5000;
Kr0=1;
imptransfer=2;
compens21=2;
else
% -------------------------------------------------------------------
Un=input('U(Fase-Fase) kV:400,220,132,66: ');
kE=input('kE=E/Enominal: ');
L=input('Longitud Linea km:  ');
Sccs=input('Sccs(MVA): ');
Sccr=input('Sccr(MVA): ');
Kr0=input('Kr0=zr0/zr1: ');
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 :');
end
% -------------------------------------------------------------------
k0=input('k0=zs0/zr0: ');
p=input('p pu :  ');

faltam=input('Falta en m: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :');
Rfm=input('Resistencia de Falta en m. Rf(Ohm) = ');
if faltam==7|faltam==8|faltam==9|faltam==11
Rgm=input('Resistencia de Falta en m. Rg(Ohm) = ');
else 
   Rgm=0;
end

faltap=input('Falta en p: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :');
Rfp=input('Resistencia de Falta en p. Rf(Ohm) = ');
if faltap==7|faltap==8|faltap==9|faltap==11
Rgp=input('Resistencia de Falta en p. Rg(Ohm) = ');
else 
   Rgp=0;
end
% -------------------------------------------------------------------
pfaltam=3;
pfaltap=4;
% -------------------------------------------------------------------
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;
fprintf('\n');

fila=1;
i=1;
m=0.01;
while m<1
if m==p
m=m+0.01;
end

switch(imptransfer)
case{1},
if m<p
[Minc]=subMincT(1);
[z0]=subz0T(zs0,zr0,zsr0,zl0,zlM0,m,p);
[z1]=subz1T(zs1,zr1,zsr1,zl1,m,p);
[z2]=subz2T(zs2,zr2,zsr2,zl2,m,p);
else
[Minc]=subMincpmT(1);
[z0]=subz0pmT(zs0,zr0,zsr0,zl0,zlM0,m,p);
[z1]=subz1pmT(zs1,zr1,zsr1,zl1,m,p);
[z2]=subz2pmT(zs2,zr2,zsr2,zl2,m,p);
end
case(2)
if m<p
[Minc]=subMinc(1);
[z0]=subz0(zs0,zr0,zl0,zlM0,m,p);
[z1]=subz1(zs1,zr1,zl1,m,p);
[z2]=subz2(zs2,zr2,zl2,m,p);
else
[Minc]=subMincpm(1);
[z0]=subz0pm(zs0,zr0,zl0,zlM0,m,p);
[z1]=subz1pm(zs1,zr1,zl1,m,p);
[z2]=subz2pm(zs2,zr2,zl2,m,p);
end 
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);

% 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]=...
   mag4mp(Z0,Z1,Z2,E,faltam,faltap,Rfm,Rgm,Rfp,Rgp,pfaltam,pfaltap);

% IBUS(0,1,2)
I0=zeros(4,1);
I1=zeros(4,1);
I2=zeros(4,1);

I0(pfaltam,1)=-Im0; 
I1(pfaltam,1)=-Im1; 
I2(pfaltam,1)=-Im2; 

I0(pfaltap,1)=-Ip0;
I1(pfaltap,1)=-Ip1;
I2(pfaltap,1)=-Ip2;

% VBUS(0,1,2)
V0=Z0*I0;
V1=Z1*I1+E;
V2=Z2*I2;

% I(0,1,2) DE LOS ELEMENTOS
[Ie0]=subIe0(Minc,V0,z0);
[Ie1]=subIe1(Minc,V1,z1);
[Ie2]=subIe2(Minc,V2,z2);

% V(0,1,2), I(0,1,2) EN LOS EXTREMOS DE LINEA
Vs0=V0(1);  Vr0=V0(2);
Vs1=V1(1);  Vr1=V1(2);
Vs2=V2(1);  Vr2=V2(2);

IsI0=Ie0(3); IrI0=Ie0(5);  IsII0=Ie0(4); IrII0=-Ie0(4);
IsI1=Ie1(3); IrI1=Ie1(5);  IsII1=Ie1(4); IrII1=-Ie1(4);  
IsI2=Ie2(3); IrI2=Ie2(5);  IsII2=Ie2(4); IrII2=-Ie2(4);
%
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS PUNTOS m,p DE FALTA
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2);
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);
[Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2);
%
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS EXTREMOS DE LINEA
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[IsIa,IsIb,IsIc]=abc(IsI0,IsI1,IsI2);
[IrIa,IrIb,IrIc]=abc(IrI0,IrI1,IrI2);
[IsIIa,IsIIb,IsIIc]=abc(IsII0,IsII1,IsII2);
[IrIIa,IrIIb,IrIIc]=abc(IrII0,IrII1,IrII2);
%
% CALCULO DE LAS MEDIDAS DE LAS PROTECCIONES DE DISTANCIA EN
% LOS EXTREMOS DE LINEA

% PROTECCION DE DISTANCIA (21)
switch(compens21)
case{1},
% COMPENSADA
[SI_ZAB,SI_ZBC,SI_ZCA,SI_ZAG,SI_ZBG,SI_ZCG]=...
   PDistComp(Vsa,Vsb,Vsc,IsIa,IsIb,IsIc,IsI0,k,Km,IsII0);
[SII_ZAB,SII_ZBC,SII_ZCA,SII_ZAG,SII_ZBG,SII_ZCG]=...
  PDistComp(Vsa,Vsb,Vsc,IsIIa,IsIIb,IsIIc,IsII0,k,Km,IsI0);
[RI_ZAB,RI_ZBC,RI_ZCA,RI_ZAG,RI_ZBG,RI_ZCG]=...
   PDistComp(Vra,Vrb,Vrc,IrIa,IrIb,IrIc,IrI0,k,Km,IrII0);
[RII_ZAB,RII_ZBC,RII_ZCA,RII_ZAG,RII_ZBG,RII_ZCG]=...
   PDistComp(Vra,Vrb,Vrc,IrIIa,IrIIb,IrIIc,IrII0,k,Km,IrI0);
case{2},
% NO COMPENSADA
[SI_ZAB,SI_ZBC,SI_ZCA,SI_ZAG,SI_ZBG,SI_ZCG]=...
   PDist(Vsa,Vsb,Vsc,IsIa,IsIb,IsIc,IsI0,k);
[SII_ZAB,SII_ZBC,SII_ZCA,SII_ZAG,SII_ZBG,SII_ZCG]=...
  PDist(Vsa,Vsb,Vsc,IsIIa,IsIIb,IsIIc,IsII0,k);
[RI_ZAB,RI_ZBC,RI_ZCA,RI_ZAG,RI_ZBG,RI_ZCG]=...
   PDist(Vra,Vrb,Vrc,IrIa,IrIb,IrIc,IrI0,k);
[RII_ZAB,RII_ZBC,RII_ZCA,RII_ZAG,RII_ZBG,RII_ZCG]=...
   PDist(Vra,Vrb,Vrc,IrIIa,IrIIb,IrIIc,IrII0,k);
end

%
% CALCULO DE LAS MAGNITUDES DE ALIMETACION DE LA PROTECCION
% DIRECCIONAL DE TIERRA EN LOS EXTREMOS DE LINEA
[Vspol,IsIop,AngDirSI]=PDire(Vs0,IsI0);
[Vrpol,IrIop,AngDirRI]=PDire(Vr0,IrI0);
[Vspol,IsIIop,AngDirSII]=PDire(Vs0,IsII0);
[Vrpol,IrIIop,AngDirRII]=PDire(Vr0,IrII0);

% Para Salidas 1 y 2

if faltam==1
   cfaltam='ag';
   alcanceSm=abs(SI_ZAG)/cos(angle(SI_ZAG)-angle(zl1u));
   alcanceRm=abs(RI_ZAG)/cos(angle(RI_ZAG)-angle(zl1u));
   elementom='AG';

elseif faltam==2
   cfaltam='bg';
   alcanceSm=abs(SI_ZBG)/cos(angle(SI_ZBG)-angle(zl1u));
   alcanceRm=abs(RI_ZBG)/cos(angle(RI_ZBG)-angle(zl1u));
   elementom='BG';

elseif faltam==3
   cfaltam='cg';
   alcanceSm=abs(SI_ZCG)/cos(angle(SI_ZCG)-angle(zl1u));
   alcanceRm=abs(RI_ZCG)/cos(angle(RI_ZCG)-angle(zl1u));
   elementom='CG';

elseif faltam==4
   cfaltam='ab';
   alcanceSm=abs(SI_ZAB)/cos(angle(SI_ZAB)-angle(zl1u));
   alcanceRm=abs(RI_ZAB)/cos(angle(RI_ZAB)-angle(zl1u));
   elementom='AB';

elseif faltam==5
   cfaltam='bc';
   alcanceSm=abs(SI_ZBC)/cos(angle(SI_ZBC)-angle(zl1u));
   alcanceRm=abs(RI_ZBC)/cos(angle(RI_ZBC)-angle(zl1u));
   elementom='BC';

elseif faltam==6
   cfaltam='ca';
   alcanceSm=abs(SI_ZCA)/cos(angle(SI_ZCA)-angle(zl1u));
   alcanceRm=abs(RI_ZCA)/cos(angle(RI_ZCA)-angle(zl1u));
   elementom='CA';

elseif faltam==7
   cfaltam='abg';
   alcanceSm=abs(SI_ZAB)/cos(angle(SI_ZAB)-angle(zl1u));
   alcanceRm=abs(RI_ZAB)/cos(angle(RI_ZAB)-angle(zl1u));
   elementom='AB';

elseif faltam==8
   cfaltam='bcg';
   alcanceSm=abs(SI_ZBC)/cos(angle(SI_ZBC)-angle(zl1u));
   alcanceRm=abs(RI_ZBC)/cos(angle(RI_ZBC)-angle(zl1u));
   elementom='BC';

elseif faltam==9
   cfaltam='cag';
   alcanceSm=abs(SI_ZCA)/cos(angle(SI_ZCA)-angle(zl1u));
   alcanceRm=abs(RI_ZCA)/cos(angle(RI_ZCA)-angle(zl1u));
   elementom='CA';
  
elseif faltam==10
   cfaltam='abc';
   alcanceSm=abs(SI_ZBC)/cos(angle(SI_ZBC)-angle(zl1u));
   alcanceRm=abs(RI_ZBC)/cos(angle(RI_ZBC)-angle(zl1u));
   elementom='BC';
   
elseif faltam==11
   cfaltam='abcg';
   alcanceSm=abs(SI_ZBC)/cos(angle(SI_ZBC)-angle(zl1u));
   alcanceRm=abs(RI_ZBC)/cos(angle(RI_ZBC)-angle(zl1u));
   elementom='BC';
   
end

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

alcanceSmpu=alcanceSm/abs(zl1u*m*L);
alcanceRmpu=alcanceRm/abs(zl1u*(1-m)*L);
   
alcanceSmpua=alcanceSm/abs(zl1u*L);
alcanceRmpua=alcanceRm/abs(zl1u*L);

Abscisa(fila,i)=m;
OrdenadaSm(fila,i)=alcanceSmpu;
OrdenadaRm(fila,i)=alcanceRmpu;
OrdenadaSma(fila,i)=alcanceSmpua;
OrdenadaRma(fila,i)=alcanceRmpua;

% para Salida 3

alcanceSm_ZAB=abs(SI_ZAB)/cos(angle(SI_ZAB)-angle(zl1u));
alcanceSm_ZBC=abs(SI_ZBC)/cos(angle(SI_ZBC)-angle(zl1u));
alcanceSm_ZCA=abs(SI_ZCA)/cos(angle(SI_ZCA)-angle(zl1u));
alcanceSm_ZAG=abs(SI_ZAG)/cos(angle(SI_ZAG)-angle(zl1u));
alcanceSm_ZBG=abs(SI_ZBG)/cos(angle(SI_ZBG)-angle(zl1u));
alcanceSm_ZCG=abs(SI_ZCG)/cos(angle(SI_ZCG)-angle(zl1u));

alcanceSm_ZABpua=alcanceSm_ZAB/abs(zl1u*L); 
alcanceSm_ZBCpua=alcanceSm_ZBC/abs(zl1u*L);
alcanceSm_ZCApua=alcanceSm_ZCA/abs(zl1u*L);
alcanceSm_ZAGpua=alcanceSm_ZAG/abs(zl1u*L);
alcanceSm_ZBGpua=alcanceSm_ZBG/abs(zl1u*L);
alcanceSm_ZCGpua=alcanceSm_ZCG/abs(zl1u*L);

OrdenadaSm_ZABpua(1,i)=alcanceSm_ZABpua;
OrdenadaSm_ZBCpua(1,i)=alcanceSm_ZBCpua;
OrdenadaSm_ZCApua(1,i)=alcanceSm_ZCApua;
OrdenadaSm_ZAGpua(1,i)=alcanceSm_ZAGpua;
OrdenadaSm_ZBGpua(1,i)=alcanceSm_ZBGpua;
OrdenadaSm_ZCGpua(1,i)=alcanceSm_ZCGpua;

% ---------------------------------------------------------

alcanceSr_ZAB=abs(SII_ZAB)/cos(angle(SII_ZAB)-angle(zl1u));
alcanceSr_ZBC=abs(SII_ZBC)/cos(angle(SII_ZBC)-angle(zl1u));
alcanceSr_ZCA=abs(SII_ZCA)/cos(angle(SII_ZCA)-angle(zl1u));
alcanceSr_ZAG=abs(SII_ZAG)/cos(angle(SII_ZAG)-angle(zl1u));
alcanceSr_ZBG=abs(SII_ZBG)/cos(angle(SII_ZBG)-angle(zl1u));
alcanceSr_ZCG=abs(SII_ZCG)/cos(angle(SII_ZCG)-angle(zl1u));

alcanceSr_ZABpua=alcanceSr_ZAB/abs(zl1u*L); 
alcanceSr_ZBCpua=alcanceSr_ZBC/abs(zl1u*L);
alcanceSr_ZCApua=alcanceSr_ZCA/abs(zl1u*L);
alcanceSr_ZAGpua=alcanceSr_ZAG/abs(zl1u*L);
alcanceSr_ZBGpua=alcanceSr_ZBG/abs(zl1u*L);
alcanceSr_ZCGpua=alcanceSr_ZCG/abs(zl1u*L);

OrdenadaSr_ZABpua(1,i)=alcanceSr_ZABpua;
OrdenadaSr_ZBCpua(1,i)=alcanceSr_ZBCpua;
OrdenadaSr_ZCApua(1,i)=alcanceSr_ZCApua;
OrdenadaSr_ZAGpua(1,i)=alcanceSr_ZAGpua;
OrdenadaSr_ZBGpua(1,i)=alcanceSr_ZBGpua;
OrdenadaSr_ZCGpua(1,i)=alcanceSr_ZCGpua;

% ------------------------------------------------------------

alcanceRm_ZAB=abs(RI_ZAB)/cos(angle(RI_ZAB)-angle(zl1u));
alcanceRm_ZBC=abs(RI_ZBC)/cos(angle(RI_ZBC)-angle(zl1u));
alcanceRm_ZCA=abs(RI_ZCA)/cos(angle(RI_ZCA)-angle(zl1u));
alcanceRm_ZAG=abs(RI_ZAG)/cos(angle(RI_ZAG)-angle(zl1u));
alcanceRm_ZBG=abs(RI_ZBG)/cos(angle(RI_ZBG)-angle(zl1u));
alcanceRm_ZCG=abs(RI_ZCG)/cos(angle(RI_ZCG)-angle(zl1u));

alcanceRm_ZABpua=alcanceRm_ZAB/abs(zl1u*L); 
alcanceRm_ZBCpua=alcanceRm_ZBC/abs(zl1u*L);
alcanceRm_ZCApua=alcanceRm_ZCA/abs(zl1u*L);
alcanceRm_ZAGpua=alcanceRm_ZAG/abs(zl1u*L);
alcanceRm_ZBGpua=alcanceRm_ZBG/abs(zl1u*L);
alcanceRm_ZCGpua=alcanceRm_ZCG/abs(zl1u*L);

OrdenadaRm_ZABpua(1,i)=alcanceRm_ZABpua;
OrdenadaRm_ZBCpua(1,i)=alcanceRm_ZBCpua;
OrdenadaRm_ZCApua(1,i)=alcanceRm_ZCApua;
OrdenadaRm_ZAGpua(1,i)=alcanceRm_ZAGpua;
OrdenadaRm_ZBGpua(1,i)=alcanceRm_ZBGpua;
OrdenadaRm_ZCGpua(1,i)=alcanceRm_ZCGpua;

% ------------------------------------------------------------

alcanceRs_ZAB=abs(RII_ZAB)/cos(angle(RII_ZAB)-angle(zl1u));
alcanceRs_ZBC=abs(RII_ZBC)/cos(angle(RII_ZBC)-angle(zl1u));
alcanceRs_ZCA=abs(RII_ZCA)/cos(angle(RII_ZCA)-angle(zl1u));
alcanceRs_ZAG=abs(RII_ZAG)/cos(angle(RII_ZAG)-angle(zl1u));
alcanceRs_ZBG=abs(RII_ZBG)/cos(angle(RII_ZBG)-angle(zl1u));
alcanceRs_ZCG=abs(RII_ZCG)/cos(angle(RII_ZCG)-angle(zl1u));

alcanceRs_ZABpua=alcanceRs_ZAB/abs(zl1u*L); 
alcanceRs_ZBCpua=alcanceRs_ZBC/abs(zl1u*L);
alcanceRs_ZCApua=alcanceRs_ZCA/abs(zl1u*L);
alcanceRs_ZAGpua=alcanceRs_ZAG/abs(zl1u*L);
alcanceRs_ZBGpua=alcanceRs_ZBG/abs(zl1u*L);
alcanceRs_ZCGpua=alcanceRs_ZCG/abs(zl1u*L);

OrdenadaRs_ZABpua(1,i)=alcanceRs_ZABpua;
OrdenadaRs_ZBCpua(1,i)=alcanceRs_ZBCpua;
OrdenadaRs_ZCApua(1,i)=alcanceRs_ZCApua;
OrdenadaRs_ZAGpua(1,i)=alcanceRs_ZAGpua;
OrdenadaRs_ZBGpua(1,i)=alcanceRs_ZBGpua;
OrdenadaRs_ZCGpua(1,i)=alcanceRs_ZCGpua;

%***************************************
%Para salida 4
%***************************************
vectorVspol(1,i)=Vspol;
vectorIsIop(1,i)=IsIop;
vectorAngDirSI(1,i)=AngDirSI;

vectorVrpol(1,i)=Vrpol;
vectorIrIop(1,i)=IrIop;
vectorAngDirRI(1,i)=AngDirRI;

vectorVspol(1,i)=Vspol;
vectorIsIIop(1,i)=IsIIop;
vectorAngDirSII(1,i)=AngDirSII;

vectorVrpol(1,i)=Vrpol;
vectorIrIIop(1,i)=IrIIop;
vectorAngDirRII(1,i)=AngDirRII;

%***************************************

i=i+1;
m=m+0.01;
end

% -----------------------------------------------------------------------
if compens21==1
   comp='21c';
elseif compens21==2
   comp='21nc';
end

cadenam1=sprintf('PG124mlc,%s,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,Kr0=%-5.2f',...
   comp,Un,kE,L,Sccs,Sccr,Kr0);
cadenamk=sprintf('k0=%-4.1f',k0);
cadenam1=sprintf('%s,%s\n',cadenam1,cadenamk);

if imptransfer==1
cadenamT=sprintf('zsr1=%-6.2fi,Ksr0=%-4.2f',imag(zsr1),Ksr0);
cadenam1=sprintf('%s,%s',cadenam1,cadenamT);   
end
cadenam2=sprintf('m:%s,%s,Rfm=%-3.0f',cfaltam,elementom,Rfm);
if faltam==1|faltam==2|faltam==3|faltam==4|faltam==5|faltam==6|faltam==10
cadenam=sprintf('%s,%s',cadenam1,cadenam2);
elseif faltam==7|faltam==8|faltam==9|faltam==11
cadenam3=sprintf('Rgm=%-3.0f',Rgm);
cadenam=sprintf('%s,%s,%s',cadenam1,cadenam2,cadenam3);
end
cadenap=sprintf('p=%-4.2f, %s, Rfp=%-3.0f',p,cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
cadenap=sprintf('p=%-4.2f, %s, Rfp=%-3.0f, Rgp=%-3.0f',p,cfaltap,Rfp,Rgp);
end
cadenamp=sprintf('%s, %s',cadenam,cadenap);
% -----------------------------------------------------------------------
cadenamc1=sprintf('PG124mlc,%s,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,Kr0=%-5.2f,k0=%-4.1f\n',...
   comp,Un,kE,L,Sccs,Sccr,Kr0,k0);
if imptransfer==1
cadenamcT=sprintf('zsr1=%-6.2f,Ksr0=%-4.2f',imag(zsr1),Ksr0);
cadenamc1=sprintf('%s,%s',cadenamc1,cadenamcT);   
end
cadenamc2=sprintf('m:%s,Rfm=%-3.0f',cfaltam,Rfm);
if faltam==1|faltam==2|faltam==3|faltam==4|faltam==5|faltam==6|faltam==10
cadenamc=sprintf('%s,%s',cadenamc1,cadenamc2);
   
elseif faltam==7|faltam==8|faltam==9|faltam==11
cadenamc3=sprintf('Rgm=%-3.0f',Rgm);
cadenamc=sprintf('%s,%s,%s',cadenamc1,cadenamc2,cadenamc3);
end
cadenapc=sprintf('p=%-4.2f, %s, Rfp=%-3.0f',p,cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
cadenapc=sprintf('p=%-4.2f, %s, Rfp=%-3.0f, Rgp=%-3.0f',p,cfaltap,Rfp,Rgp);
end
cadenampc=sprintf('%s, %s',cadenamc,cadenapc);

% -----------------------------------------------------------------------

%***************************************
% Titulo para las gráficas 67N
%***************************************
dircadenamc1=sprintf('PG124mlc,67N,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,Kr0=%-5.2f,k0=%-4.1f\n',...
   Un,kE,L,Sccs,Sccr,Kr0,k0);
if imptransfer==1
dircadenamcT=sprintf('zsr1=%-6.2f,Ksr0=%-4.2f',imag(zsr1),Ksr0);
dircadenamc1=sprintf('%s,%s',dircadenamc1,dircadenamcT);   
end
dircadenamc2=sprintf('m:%s,Rfm=%-3.0f',cfaltam,Rfm);
if faltam==1|faltam==2|faltam==3|faltam==4|faltam==5|faltam==6|faltam==10
dircadenamc=sprintf('%s,%s',dircadenamc1,dircadenamc2);
   
elseif faltam==7|faltam==8|faltam==9|faltam==11
dircadenamc3=sprintf('Rgm=%-3.0f',Rgm);
dircadenamc=sprintf('%s,%s,%s',dircadenamc1,dircadenamc2,dircadenamc3);
end
dircadenapc=sprintf('p=%-4.2f, %s, Rfp=%-3.0f',p,cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
dircadenapc=sprintf('p=%-4.2f, %s, Rfp=%-3.0f, Rgp=%-3.0f',p,cfaltap,Rfp,Rgp);
end
dircadenampc=sprintf('%s, %s',dircadenamc,dircadenapc);

%*************************************

variassalidas=1;
while variassalidas==1
   
fprintf('Salida Grafica 1: Sm,Rm, Unidad Preferente (pu)\n');
fprintf('Salida Grafica 2: Sm,Rm, Unidad Preferente (pua)\n');
fprintf('Salida Grafica 3: Sm, Todas las Unidades (pua)\n');
fprintf('                : Sr, Todas las Unidades (pua)\n');
fprintf('                : Rm, Todas las Unidades (pua)\n');
fprintf('                : Rs, Todas las Unidades (pua)\n');
fprintf('Salida Grafica 4: S(Linea I), 67N\n');
fprintf('                : R(Linea I), 67N\n');
fprintf('                : S(Linea II), 67N\n');
fprintf('                : R(Linea II), 67N\n');
   
salida=input('Salida: 1,2,3,4: ');
switch(salida)
       case{1}, % *******  Salida Grafica 1 : pu **************

mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
Ordenadadi(cont)=0.8/mm;
Ordenadadd(cont)=0.8/(1-mm);
cont=cont+1;
mm=mm+0.01;
end

ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');

subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'g',Abscisa,OrdenadaSm,'r'),grid on
title(cadenamp)
axis([0 1 ejemenosy ejemasy]);
xlabel('m'),ylabel('alcanceSmpu')
legend('Ref','Lim','alcanceSmpu')

subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadd,'g',Abscisa,OrdenadaRm,'r'),grid on
%title(cadenamp)
axis([0 1 ejemenosy ejemasy]);
xlabel('m'),ylabel('alcanceRmpu')
legend('Ref','Lim','alcanceRmpu')

       case(2)  % *******  Salida Grafica 2 : pua **************
   
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=mm;
Abscisa00(cont)=mm;

Ordenada00(cont)=1-mm;
Ordenadaz1(cont)=0.8;
cont=cont+1;
mm=mm+0.01;
end

subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadaz1,'g',Abscisa,OrdenadaSma,'r'),grid on
title(cadenamp)
%axis([0 1 0 2]);
xlabel('m'),ylabel('alcanceSmpua')
legend('Ref','z1','alcanceSmpua')

subplot(2,1,2)
plot(Abscisa00,Ordenada00,'c',Abscisa0,Ordenadaz1,'g',Abscisa,OrdenadaRma,'r'),grid on
axis([0 1 -2 2]);
%title(cadenamp)
xlabel('m'),ylabel('alcanceRmpua')
legend('Ref','z1','alcanceRmpua')

       case(3)  % *******  Salida Grafica 3 : Todas Unidades pua **************
          
       extremo=input('Extremo: Sm=1, Sr=2, Rm=3, Rs=4: ');

       switch(extremo)
       case{1},  % Sm

mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenadaz1(cont)=0.8;
Ordenada45(cont)=mm;
cont=cont+1;
mm=mm+0.01;
end

cadylabel1a='alcSmpua';

% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores

hold on
plot(Abscisa,OrdenadaSm_ZABpua,'r',Abscisa,OrdenadaSm_ZBCpua,'g',Abscisa,OrdenadaSm_ZCApua,'b',...
   Abscisa,OrdenadaSm_ZAGpua,'m',Abscisa,OrdenadaSm_ZBGpua,'k'),grid on
hold off

% Colores
hold on
a=plot(Abscisa,OrdenadaSm_ZCGpua);
set(a,'Color',[1 0.5490 0])

h=plot(Abscisa0,Ordenadaz1);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores


axis([0 1 -5 5]);
title(cadenampc)
xlabel('m'),ylabel(cadylabel1a)
legend('45º','AB','BC','CA','AG','BG','CG','z1')

       case(2) %  Sr
          
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenadaz1(cont)=0.8;
Ordenada45(cont)=mm;
cont=cont+1;
mm=mm+0.01;
end

cadylabel1a='alcSrpua';

plot(Abscisa,OrdenadaSr_ZABpua,'r',Abscisa,OrdenadaSr_ZBCpua,'g',Abscisa,OrdenadaSr_ZCApua,'b',...
   Abscisa,OrdenadaSr_ZAGpua,'m',Abscisa,OrdenadaSr_ZBGpua,'k'),grid on

% Colores
hold on
a=plot(Abscisa,OrdenadaSr_ZCGpua);
set(a,'Color',[1 0.5490 0])

h=plot(Abscisa0,Ordenadaz1);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores


axis([0 1 -20 20]);
%axis([0.8 1 0.8 1.2]);
title(cadenampc)
xlabel('m'),ylabel(cadylabel1a)
legend('AB','BC','CA','AG','BG','CG','z1')

       case{3}, % Rm

mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenadaz1(cont)=0.8;
Ordenada45(cont)=1-mm;
cont=cont+1;
mm=mm+0.01;
end

cadylabel1a='alcRmpua';

% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores

hold on
plot(Abscisa,OrdenadaRm_ZABpua,'r',Abscisa,OrdenadaRm_ZBCpua,'g',Abscisa,OrdenadaRm_ZCApua,'b',...
   Abscisa,OrdenadaRm_ZAGpua,'m',Abscisa,OrdenadaRm_ZBGpua,'k'),grid on
hold off

% Colores
hold on
a=plot(Abscisa,OrdenadaRm_ZCGpua);
set(a,'Color',[1 0.5490 0])

h=plot(Abscisa0,Ordenadaz1);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores


axis([0 1 -5 5]);
title(cadenampc)
xlabel('m'),ylabel(cadylabel1a)
legend('-45º','AB','BC','CA','AG','BG','CG','z1')

       case{4},  % Rs

mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenadaz1(cont)=0.8;
Ordenada45(cont)=1-mm;
cont=cont+1;
mm=mm+0.01;
end

cadylabel1a='alcRspua';

% Colores
%b=plot(Abscisa0,Ordenada45);
%set(b,'Color','c')
% Fin Colores

%hold on
plot(Abscisa,OrdenadaRs_ZABpua,'r',Abscisa,OrdenadaRs_ZBCpua,'g',Abscisa,OrdenadaRs_ZCApua,'b',...
   Abscisa,OrdenadaRs_ZAGpua,'m',Abscisa,OrdenadaRs_ZBGpua,'k'),grid on
%hold off

% Colores
hold on
a=plot(Abscisa,OrdenadaRs_ZCGpua);
set(a,'Color',[1 0.5490 0])

h=plot(Abscisa0,Ordenadaz1);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores


axis([0 1 -5 5]);
title(cadenampc)
xlabel('m'),ylabel(cadylabel1a)
legend('AB','BC','CA','AG','BG','CG','z1')

end  % end del switch extremo de la salida 3

       case{4},   % Proteccion Direccional
          
         extr=input('Extremo: S(Linea I)=1, R(Linea I)=2,S(Linea II)=3, R(Linea II)=4: ');
   
         switch(extr)
         case{1}
            
   subplot(3,1,1)
   plot(Abscisa,vectorAngDirSI,'b');
   axis([0 1 -180 180]);
   xlabel('m'),ylabel('AngDirSI(grados)')
   title(dircadenampc)
   rayar(1);
   legend('AngDirSI','AngMaxPar')
   
   subplot(3,1,2)
   plot(Abscisa,vectorVspol,'b');
   grid   
   %title(dircadenampc)
   xlabel('m'),ylabel('modulo Vspol(kV)')

   subplot(3,1,3)
   plot(Abscisa,vectorIsIop,'b');
   grid   
   xlabel('m'),ylabel('modulo IsIop(kA)')
   
case{2}
   
   subplot(3,1,1)
   plot(Abscisa,vectorAngDirRI,'b');
   axis([0 1 -180 180]);
   xlabel('m'),ylabel('AngDirRI(grados)')
   title(dircadenampc)
   rayar(1);
   legend('AngDirRI','AngMaxPar')
   
   subplot(3,1,2)
   plot(Abscisa,vectorVrpol,'b');
   grid   
   %title(dircadenampc)
   xlabel('m'),ylabel('modulo Vrpol(kV)')

   subplot(3,1,3)
   plot(Abscisa,vectorIrIop,'b');
   grid   
   xlabel('m'),ylabel('modulo IrIop(kA)')
   
            case{3}
               
   subplot(3,1,1)
   plot(Abscisa,vectorAngDirSII,'b');
   axis([0 1 -180 180]);
   xlabel('m'),ylabel('AngDirSII(grados)')
   title(dircadenampc)
   rayar(1);
   legend('AngDirSII','AngMaxPar')
   
   subplot(3,1,2)
   plot(Abscisa,vectorVspol,'b');
   grid   
   %title(dircadenampc)
   xlabel('m'),ylabel('modulo Vspol(kV)')

   subplot(3,1,3)
   plot(Abscisa,vectorIsIIop,'b');
   grid   
   xlabel('m'),ylabel('modulo IsIIop(kA)')
   
case{4}
   
 subplot(3,1,1)
      plot(Abscisa,vectorAngDirRII,'b');
   axis([0 1 -180 180]);
   xlabel('m'),ylabel('AngDirRII(grados)')
   title(dircadenampc)
   rayar(1);
   legend('AngDirRII','AngMaxPar')
   
   subplot(3,1,2)
   plot(Abscisa,vectorVrpol,'b');
   grid   
   %title(dircadenampc)
   xlabel('m'),ylabel('modulo Vrpol(kV)')

   subplot(3,1,3)
   plot(Abscisa,vectorIrIIop,'b');
   grid   
   xlabel('m'),ylabel('modulo IrIIop(kA)')

end %end del switch extr
          
          
end  % end del switch salida

variassalidas=input('Mas Salidas=1,Acabar=2 :');
end  % end del while variassalidas

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%252FPGc%252F">PGc</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPGc%252FPG124mlc%252F">PG124mlc</a>&nbsp;&nbsp;>&nbsp;&nbsp;<a href="use/abreproy/4273/direccion/z_Prog%252FPGc%252FPG124mlc%252FFG124mlc.m%252F">FG124mlc.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>