% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG128Mca\n');
fprintf('Esquema Linea Doble SR\n');
fprintf('Unominal AT kV: 220\n');
fprintf('Unominal BT kV: 66\n');
fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf('Eleccion Ctes. del Tf\n');
fprintf('Tf AT/BT/T: Neutro a Tierra en BT\n');
fprintf('Eleccion: kE,LA=LB=L,Sccs,Sccr,Kr0,k0\n');
fprintf('Impedancia de Transferencia Si, No\n');
fprintf('Eleccion: (zsr1,Ksr0)\n');
fprintf('Compensacion 21, Si, No\n');
fprintf('Falta Intercircuito entre Lineas de Diferente Tension\n');
fprintf('m=p=0...1\n');
fprintf('Barra m=1, Barra p=2\n');
%fprintf('Salida Grafica 1: Analisis SAm, Todas las Unidades (pu)\n');
%fprintf(' : Analisis SBp, Todas las Unidades (pu)\n');
%fprintf(' : Analisis RAm, Todas las Unidades (pu)\n');
%fprintf(' : Analisis RBp, Todas las Unidades (pu)\n');
%fprintf('Salida Grafica 2: Analisis SAm, Todas las Unidades (pua)\n');
%fprintf(' : Analisis SBp, Todas las Unidades (pua)\n');
%fprintf(' : Analisis RAm, Todas las Unidades (pua)\n');
%fprintf(' : Analisis RBp, Todas las Unidades (pua)\n');
%fprintf('Salida Grafica 3: Analisis SAm, 67N\n');
%fprintf(' : Analisis RAm, 67N\n');
%fprintf(' : Analisis SBp, 67N\n');
%fprintf(' : Analisis RBp, 67N\n');
%fprintf('Salida Grafica 4: Analisis SAm, Todas las Unidades R-X\n');
%fprintf(' : Analisis RAm, Todas las Unidades R-X\n');
%fprintf(' : Analisis SBp, Todas las Unidades R-X\n');
%fprintf(' : Analisis RBp, Todas las Unidades R-X\n');
fprintf('\n');
% -----------------------------------------------------------------------
%PG128Mca
datos=input('datos: standard=1,input=2: ');
if datos==1
UA=220;
UB=66;
L=66;
kE=1;
Sp=120; Ss=120; St=40; Up=220; Us=66; Ut=33;
tccps=12; tccpt=5; tccst=2;
Sccs=5000;
Sccr=5000;
imptransfer=2;
compens21=2;
Kr0=1;
else
% -----------------------------------------------------------------------
UA=input('UA(Fase-Fase) kV, poner 220: ');
UB=input('UB(Fase-Fase) kV, poner 66: ');
Ut=input('Ut(Fase-Fase) kV, poner 33: ');
kE=input('kE: ');
L=input('Longitud L=LA=LB km: ');
Sccs=input('Sccs MVA: ');
Sccr=input('Sccr MVA: ');
Kr0=input('Kr0=zr0/zr1: ');
Sp=input('Sp MVA, poner 120: ');
Ss=input('Ss MVA, poner 120: ');
St=input('St MVA, poner 40: ');
tccps=input('tccps %, Ref a Sp, poner 12: ');
tccpt=input('tccpt %, Ref a St, poner 5: ');
tccst=input('tccst %, Ref a St, poner 2: ');
imptransfer=input('imptransfer:si=1,no=2: ');
compens21=input('Compensacion 21: Si=1, No=2 :');
end
% -----------------------------------------------------------------------
k0=input('k0=zs0/zr0 (Tipicos=0.1-1-10): ');
fprintf('Tipo de Falta: \n');
fprintf('Aa=1,Bc=2,Cb=3,Aag=4,Bcg=5,Cbg=6,Aab=7,Bac=8,Cbc=9\n');
fprintf('Aabg=10,Bacg=11,Cbcg=12,Aabc=13,Aabcg=14,ABCabc=15,ABCabcg=16\n');
falta=input('Elegir Tipo de Falta: ');
Rf=input('Resistencia de Falta. Rf(Ohm) = ');
if falta==5|falta==6|falta==10|falta==11|falta==12|falta==14|falta==16
Rg=input('Resistencia comun de Falta. Rg(Ohm) = ');
else
Rg=0;
end
% -----------------------------------------------------------------------
if UA==400
zlA1u=0.0346+0.3208i; zlA2u=zlA1u; zlA0u=0.1828+1.1482i; zlM0u=0.1482+0.7383i;
elseif UA==220
zlA1u=0.0829+0.4135i; zlA2u=zlA1u; zlA0u=0.2311+1.3085i; zlM0u=0.1482+0.8031i;
elseif UA==132
zlA1u=0.1342+0.4042i; zlA2u=zl1u; zlA0u=0.2824+1.3718i; zlM0u=0.1482+0.8947i;
elseif UA==66
zlA1u=0.2156+0.3952i; zlA2u=zl1u; zlA0u=0.3636+1.4278i; zlM0u=0.1482+0.9288i;
end
if UB==400
zlB1u=0.0346+0.3208i; zlB2u=zlB1u; zlB0u=0.1828+1.1482i;
elseif UB==220
zlB1u=0.0829+0.4135i; zlB2u=zlB1u; zlB0u=0.2311+1.3085i;
elseif UB==132
zlB1u=0.1342+0.4042i; zlB2u=zlB1u; zlB0u=0.2824+1.3718i;
elseif UB==66
zlB1u=0.2156+0.3952i; zlB2u=zlB1u; zlB0u=0.3636+1.4278i;
end
% -----------------------------------------------------------------------
LA=L;
LB=L;
%EAv=UA*1000/sqrt(3)+0i; EA=1.0;
%EBv=UB*1000/sqrt(3)+0i; EB=1.0;
E=1.0;
E=kE*E;
kA=(zlA0u-zlA1u)/(3*zlA1u);
kB=(zlB0u-zlB1u)/(3*zlB1u);
Sbase=100;
UbaseA=UA;
UbaseB=UB;
VbaseA=UbaseA/sqrt(3);
VbaseB=UbaseB/sqrt(3);
zbaseA=UA^2/Sbase;
zbaseB=UB^2/Sbase;
IbaseA=(UA/sqrt(3))/zbaseA;
IbaseB=(UB/sqrt(3))/zbaseB;
KmA=zlM0u/zlA1u; KmA=KmA*(IbaseB/IbaseA);
KmB=zlM0u/zlB1u; KmB=KmB*(IbaseA/IbaseB);
fprintf('\n');
if imptransfer==1
zsr1=input('zsr1=zsr2(Complejo r+xi en Ohms): ');
zsr1=zsr1/zbaseA;
zsr2=zsr1;
Ksr0=input('Ksr0=zsr0/zsr1= ');
zsr0=Ksr0*zsr1;
end
if falta==2 cfalta='Bc';elseif falta==3 cfalta='Cb';elseif falta==5 cfalta='Bcg';
elseif falta==6 cfalta='Cbg';elseif falta==7 cfalta='Aab';elseif falta==8 cfalta='Bac';
elseif falta==9 cfalta='Cbc';elseif falta==10 cfalta='Aabg';elseif falta==11 cfalta='Bacg';
elseif falta==12 cfalta='Cbcg';elseif falta==13 cfalta='Aabc';elseif falta==14 cfalta='Aabcg';
elseif falta==15 cfalta='ABCabc';elseif falta==16 cfalta='ABCabcg';
end
i=1;
m=0.01;
while m<1;
p=m;
% REDES DE SECUENCIA 0,1,2 pu
[zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1]=pusec1(UA,Sccs,Sccr,zlA1u,m,LA,zlB1u,p,LB,zbaseA,zbaseB);
[zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2]=pusec2(UA,Sccs,Sccr,zlA2u,m,LA,zlB2u,p,LB,zbaseA,zbaseB);
[zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0]=...
pusec0(k0,Kr0,zr1,zlA0u,m,LA,zlB0u,p,LB,zlM0u,zbaseA,zbaseB,UbaseA,UbaseB,Sbase);
zpTS=1/2*((tccps/100)*UA^2/Sp+(tccpt/100)*UA^2/St-(tccst/100)*UA^2/St); zpTS=zpTS/zbaseA;
zsTS=1/2*((tccps/100)*UA^2/Sp-(tccpt/100)*UA^2/St+(tccst/100)*UA^2/St); zsTS=zsTS/zbaseA;
ztTS=1/2*(-(tccps/100)*UA^2/Sp+(tccpt/100)*UA^2/St+(tccst/100)*UA^2/St); ztTS=ztTS/zbaseA;
zpTR=zpTS; zsTR=zsTS; ztTR=ztTS;
% CALCULO DE LA REDUCCION DE REDES 0,1,2 A LOS PUNTOS m,p DE FALTA
switch(imptransfer)
case{1},
[Minc0]=subMinc0T(1);
[Minc1]=subMinc1T(1);
[Minc2]=subMinc2T(1);
[z1]=subz1T(zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1,zpTS,zsTS,zpTR,zsTR,zsr1);
[z2]=subz2T(zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2,zpTS,zsTS,zpTR,zsTR,zsr2);
[z0]=subz0T(zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0,zsTS,ztTS,zsTR,ztTR,zsr0);
case{2},
[Minc0]=subMinc0(1);
[Minc1]=subMinc1(1);
[Minc2]=subMinc2(1);
[z1]=subz1(zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1,zpTS,zsTS,zpTR,zsTR);
[z2]=subz2(zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2,zpTS,zsTS,zpTR,zsTR);
[z0]=subz0(zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0,zsTS,ztTS,zsTR,ztTR);
end
% MATRIZ DE ADMITANCIA DE BARRAS
[Y0]=subY0(Minc0,z0);
[Y1]=subY1(Minc1,z1);
[Y2]=subY2(Minc2,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]=...
maMDmp(Z0,Z1,Z2,E,falta,Rf,Rg,VbaseA,VbaseB,IbaseA,IbaseB);
% IBUS(0,1,2)
I0=[-Im0;-Ip0;0;0;0;0];
I1=[-Im1;-Ip1;0;0;0;0];
I2=[-Im2;-Ip2;0;0;0;0];
% VBUS(0,1,2)
V0=Z0*I0;
V1=Z1*I1+E;
V2=Z2*I2;
% I(0,1,2) DE LOS ELEMENTOS
[Ie0]=subIe0(Minc0,V0,z0);
[Ie1]=subIe1(Minc1,V1,z1);
[Ie2]=subIe2(Minc2,V2,z2);
VSA0=V0(3); VSA1=V1(3); VSA2=V2(3);
VRA0=V0(4); VRA1=V1(4); VRA2=V2(4);
VSB0=V0(5); VSB1=V1(5); VSB2=V2(5);
VRB0=V0(6); VRB1=V1(6); VRB2=V2(6);
ISAm0=Ie0(3); ISAm1=Ie1(3); ISAm2=Ie2(3);
IRAm0=Ie0(4); IRAm1=Ie1(4); IRAm2=Ie2(4);
ISBp0=Ie0(7); ISBp1=Ie1(7); ISBp2=Ie2(7);
IRBp0=Ie0(8); IRBp1=Ie1(8); IRBp2=Ie2(8);
% 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
[VSAa,VSAb,VSAc]=abc(VSA0,VSA1,VSA2);
[VRAa,VRAb,VRAc]=abc(VRA0,VRA1,VRA2);
[VSBa,VSBb,VSBc]=abc(VSB0,VSB1,VSB2);
[VRBa,VRBb,VRBc]=abc(VRB0,VRB1,VRB2);
[ISAma,ISAmb,ISAmc]=abc(ISAm0,ISAm1,ISAm2);
[ISBpa,ISBpb,ISBpc]=abc(ISBp0,ISBp1,ISBp2);
[IRAma,IRAmb,IRAmc]=abc(IRAm0,IRAm1,IRAm2);
[IRBpa,IRBpb,IRBpc]=abc(IRBp0,IRBp1,IRBp2);
%
% PROTECCION DE DISTANCIA (21)
switch(compens21)
case{1},
% COMPENSADA
[SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG]=...
PDistComp(VSAa,VSAb,VSAc,ISAma,ISAmb,ISAmc,ISAm0,kA,KmA,ISBp0);
[RAm_ZAB,RAm_ZBC,RAm_ZCA,RAm_ZAG,RAm_ZBG,RAm_ZCG]=...
PDistComp(VRAa,VRAb,VRAc,IRAma,IRAmb,IRAmc,IRAm0,kA,KmA,IRBp0);
[SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG]=...
PDistComp(VSBa,VSBb,VSBc,ISBpa,ISBpb,ISBpc,ISBp0,kB,KmB,ISAm0);
[RBp_ZAB,RBp_ZBC,RBp_ZCA,RBp_ZAG,RBp_ZBG,RBp_ZCG]=...
PDistComp(VRBa,VRBb,VRBc,IRBpa,IRBpb,IRBpc,IRBp0,kB,KmB,IRAm0);
case{2},
% NO COMPENSADA
[SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG]=...
PDist(VSAa,VSAb,VSAc,ISAma,ISAmb,ISAmc,ISAm0,kA);
[RAm_ZAB,RAm_ZBC,RAm_ZCA,RAm_ZAG,RAm_ZBG,RAm_ZCG]=...
PDist(VRAa,VRAb,VRAc,IRAma,IRAmb,IRAmc,IRAm0,kA);
[SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG]=...
PDist(VSBa,VSBb,VSBc,ISBpa,ISBpb,ISBpc,ISBp0,kB);
[RBp_ZAB,RBp_ZBC,RBp_ZCA,RBp_ZAG,RBp_ZBG,RBp_ZCG]=...
PDist(VRBa,VRBb,VRBc,IRBpa,IRBpb,IRBpc,IRBp0,kB);
end
%
% CALCULO DE LAS MAGNITUDES DE ALIMETACION DE LA PROTECCION
% DIRECCIONAL DE TIERRA EN LOS EXTREMOS DE LINEA (pu)
%
[VSApol,ISAmop,AngDirSAm]=PDire(VSA0,ISAm0);
[VSBpol,ISBpop,AngDirSBp]=PDire(VSB0,ISBp0);
[VRApol,IRAmop,AngDirRAm]=PDire(VRA0,IRAm0);
[VRBpol,IRBpop,AngDirRBp]=PDire(VRB0,IRBp0);
%
% Conversion I(0,1,2)pu a I(0,1,2) kA
[IIm0,IIm1,IIm2]=conv3(Im0,Im1,Im2,IbaseA);
[IIp0,IIp1,IIp2]=conv3(Ip0,Ip1,Ip2,IbaseB);
[IISAm0,IISAm1,IISAm2,IIRAm0,IIRAm1,IIRAm2]=...
conv6(ISAm0,ISAm1,ISAm2,IRAm0,IRAm1,IRAm2,IbaseA);
[IISBp0,IISBp1,IISBp2,IIRBp0,IIRBp1,IIRBp2]=...
conv6(ISBp0,ISBp1,ISBp2,IRBp0,IRBp1,IRBp2,IbaseB);
% Conversion V(0,1,2)pu a V(0,1,2) kV
[VVm0,VVm1,VVm2,VVSA0,VVSA1,VVSA2,VVRA0,VVRA1,VVRA2]=...
conv9(Vm0,Vm1,Vm2,VSA0,VSA1,VSA2,VRA0,VRA1,VRA2,VbaseA);
[VVp0,VVp1,VVp2,VVSB0,VVSB1,VVSB2,VVRB0,VVRB1,VVRB2]=...
conv9(Vp0,Vp1,Vp2,VSB0,VSB1,VSB2,VRB0,VRB1,VRB2,VbaseB);
% Conversion I(a,b,c)pu a I(a,b,c) kA
[IIma,IImb,IImc]=conv3(Ima,Imb,Imc,IbaseA);
[IIpa,IIpb,IIpc]=conv3(Ipa,Ipb,Ipc,IbaseB);
[IISAma,IISAmb,IISAmc,IIRAma,IIRAmb,IIRAmc]=...
conv6(ISAma,ISAmb,ISAmc,IRAma,IRAmb,IRAmc,IbaseA);
[IISBpa,IISBpb,IISBpc,IIRBpa,IIRBpb,IIRBpc]=...
conv6(ISBpa,ISBpb,ISBpc,IRBpa,IRBpb,IRBpc,IbaseB);
% Conversion V(a,b,c)pu a V(a,b,c) kV
[VVma,VVmb,VVmc]=conv3(Vma,Vmb,Vmc,VbaseA);
[VVpa,VVpb,VVpc]=conv3(Vpa,Vpb,Vpc,VbaseB);
[VVSAa,VVSAb,VVSAc,VVRAa,VVRAb,VVRAc]=...
conv6(VSAa,VSAb,VSAc,VRAa,VRAb,VRAc,VbaseA);
[VVSBa,VVSBb,VVSBc,VVRBa,VVRBb,VVRBc]=...
conv6(VSBa,VSBb,VSBc,VRBa,VRBb,VRBc,VbaseB);
% 21: Conversion Z(pu) a Z(Ohm)
[SAm_ZZAB,SAm_ZZBC,SAm_ZZCA,SAm_ZZAG,SAm_ZZBG,SAm_ZZCG]=...
conv6(SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG,zbaseA);
[RAm_ZZAB,RAm_ZZBC,RAm_ZZCA,RAm_ZZAG,RAm_ZZBG,RAm_ZZCG]=...
conv6(RAm_ZAB,RAm_ZBC,RAm_ZCA,RAm_ZAG,RAm_ZBG,RAm_ZCG,zbaseA);
[SBp_ZZAB,SBp_ZZBC,SBp_ZZCA,SBp_ZZAG,SBp_ZZBG,SBp_ZZCG]=...
conv6(SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG,zbaseB);
[RBp_ZZAB,RBp_ZZBC,RBp_ZZCA,RBp_ZZAG,RBp_ZZBG,RBp_ZZCG]=...
conv6(RBp_ZAB,RBp_ZBC,RBp_ZCA,RBp_ZAG,RBp_ZBG,RBp_ZCG,zbaseB);
% 67N: Conversion pu a kV
[VVSApol,VVRApol]=conv2(VSApol,VRApol,VbaseA);
[VVSBpol,VVRBpol]=conv2(VSBpol,VRBpol,VbaseB);
% 67N: Conversion pu a kA
[IISAmop,IIRAmop]=conv2(ISAmop,IRAmop,IbaseA);
[IISBpop,IIRBpop]=conv2(ISBpop,IRBpop,IbaseB);
alcSAmAB=abs(SAm_ZZAB)/cos(angle(SAm_ZZAB)-angle(zlA1u));
alcSAmBC=abs(SAm_ZZBC)/cos(angle(SAm_ZZBC)-angle(zlA1u));
alcSAmCA=abs(SAm_ZZCA)/cos(angle(SAm_ZZCA)-angle(zlA1u));
alcSAmAG=abs(SAm_ZZAG)/cos(angle(SAm_ZZAG)-angle(zlA1u));
alcSAmBG=abs(SAm_ZZBG)/cos(angle(SAm_ZZBG)-angle(zlA1u));
alcSAmCG=abs(SAm_ZZCG)/cos(angle(SAm_ZZCG)-angle(zlA1u));
alcSBpAB=abs(SBp_ZZAB)/cos(angle(SBp_ZZAB)-angle(zlB1u));
alcSBpBC=abs(SBp_ZZBC)/cos(angle(SBp_ZZBC)-angle(zlB1u));
alcSBpCA=abs(SBp_ZZCA)/cos(angle(SBp_ZZCA)-angle(zlB1u));
alcSBpAG=abs(SBp_ZZAG)/cos(angle(SBp_ZZAG)-angle(zlB1u));
alcSBpBG=abs(SBp_ZZBG)/cos(angle(SBp_ZZBG)-angle(zlB1u));
alcSBpCG=abs(SBp_ZZCG)/cos(angle(SBp_ZZCG)-angle(zlB1u));
alcRAmAB=abs(RAm_ZZAB)/cos(angle(RAm_ZZAB)-angle(zlA1u));
alcRAmBC=abs(RAm_ZZBC)/cos(angle(RAm_ZZBC)-angle(zlA1u));
alcRAmCA=abs(RAm_ZZCA)/cos(angle(RAm_ZZCA)-angle(zlA1u));
alcRAmAG=abs(RAm_ZZAG)/cos(angle(RAm_ZZAG)-angle(zlA1u));
alcRAmBG=abs(RAm_ZZBG)/cos(angle(RAm_ZZBG)-angle(zlA1u));
alcRAmCG=abs(RAm_ZZCG)/cos(angle(RAm_ZZCG)-angle(zlA1u));
alcRBpAB=abs(RBp_ZZAB)/cos(angle(RBp_ZZAB)-angle(zlB1u));
alcRBpBC=abs(RBp_ZZBC)/cos(angle(RBp_ZZBC)-angle(zlB1u));
alcRBpCA=abs(RBp_ZZCA)/cos(angle(RBp_ZZCA)-angle(zlB1u));
alcRBpAG=abs(RBp_ZZAG)/cos(angle(RBp_ZZAG)-angle(zlB1u));
alcRBpBG=abs(RBp_ZZBG)/cos(angle(RBp_ZZBG)-angle(zlB1u));
alcRBpCG=abs(RBp_ZZCG)/cos(angle(RBp_ZZCG)-angle(zlB1u));
% para Salida 1
alcSAmABpu=alcSAmAB/abs(zlA1u*m*L);
alcSAmBCpu=alcSAmBC/abs(zlA1u*m*L);
alcSAmCApu=alcSAmCA/abs(zlA1u*m*L);
alcSAmAGpu=alcSAmAG/abs(zlA1u*m*L);
alcSAmBGpu=alcSAmBG/abs(zlA1u*m*L);
alcSAmCGpu=alcSAmCG/abs(zlA1u*m*L);
alcSBpABpu=alcSBpAB/abs(zlB1u*m*L);
alcSBpBCpu=alcSBpBC/abs(zlB1u*m*L);
alcSBpCApu=alcSBpCA/abs(zlB1u*m*L);
alcSBpAGpu=alcSBpAG/abs(zlB1u*m*L);
alcSBpBGpu=alcSBpBG/abs(zlB1u*m*L);
alcSBpCGpu=alcSBpCG/abs(zlB1u*m*L);
alcRAmABpu=alcRAmAB/abs(zlA1u*(1-m)*L);
alcRAmBCpu=alcRAmBC/abs(zlA1u*(1-m)*L);
alcRAmCApu=alcRAmCA/abs(zlA1u*(1-m)*L);
alcRAmAGpu=alcRAmAG/abs(zlA1u*(1-m)*L);
alcRAmBGpu=alcRAmBG/abs(zlA1u*(1-m)*L);
alcRAmCGpu=alcRAmCG/abs(zlA1u*(1-m)*L);
alcRBpABpu=alcRBpAB/abs(zlB1u*(1-m)*L);
alcRBpBCpu=alcRBpBC/abs(zlB1u*(1-m)*L);
alcRBpCApu=alcRBpCA/abs(zlB1u*(1-m)*L);
alcRBpAGpu=alcRBpAG/abs(zlB1u*(1-m)*L);
alcRBpBGpu=alcRBpBG/abs(zlB1u*(1-m)*L);
alcRBpCGpu=alcRBpCG/abs(zlB1u*(1-m)*L);
cadylabel1='alcSAmpu';
cadylabel2='alcSBppu';
cadylabel3='alcRAmpu';
cadylabel4='alcRBppu';
% para Salida 2
alcSAmABpua=alcSAmAB/abs(zlA1u*L);
alcSAmBCpua=alcSAmBC/abs(zlA1u*L);
alcSAmCApua=alcSAmCA/abs(zlA1u*L);
alcSAmAGpua=alcSAmAG/abs(zlA1u*L);
alcSAmBGpua=alcSAmBG/abs(zlA1u*L);
alcSAmCGpua=alcSAmCG/abs(zlA1u*L);
alcSBpABpua=alcSBpAB/abs(zlB1u*L);
alcSBpBCpua=alcSBpBC/abs(zlB1u*L);
alcSBpCApua=alcSBpCA/abs(zlB1u*L);
alcSBpAGpua=alcSBpAG/abs(zlB1u*L);
alcSBpBGpua=alcSBpBG/abs(zlB1u*L);
alcSBpCGpua=alcSBpCG/abs(zlB1u*L);
alcRAmABpua=alcRAmAB/abs(zlA1u*L);
alcRAmBCpua=alcRAmBC/abs(zlA1u*L);
alcRAmCApua=alcRAmCA/abs(zlA1u*L);
alcRAmAGpua=alcRAmAG/abs(zlA1u*L);
alcRAmBGpua=alcRAmBG/abs(zlA1u*L);
alcRAmCGpua=alcRAmCG/abs(zlA1u*L);
alcRBpABpua=alcRBpAB/abs(zlB1u*L);
alcRBpBCpua=alcRBpBC/abs(zlB1u*L);
alcRBpCApua=alcRBpCA/abs(zlB1u*L);
alcRBpAGpua=alcRBpAG/abs(zlB1u*L);
alcRBpBGpua=alcRBpBG/abs(zlB1u*L);
alcRBpCGpua=alcRBpCG/abs(zlB1u*L);
cadylabel1a='alcSAmpua';
cadylabel2a='alcSBppua';
cadylabel3a='alcRAmpua';
cadylabel4a='alcRBppua';
Abscisa(i)=m;
OrdenadaSAmAB(1,i)=alcSAmABpu;
OrdenadaSAmBC(1,i)=alcSAmBCpu;
OrdenadaSAmCA(1,i)=alcSAmCApu;
OrdenadaSAmAG(1,i)=alcSAmAGpu;
OrdenadaSAmBG(1,i)=alcSAmBGpu;
OrdenadaSAmCG(1,i)=alcSAmCGpu;
OrdenadaSBpAB(1,i)=alcSBpABpu;
OrdenadaSBpBC(1,i)=alcSBpBCpu;
OrdenadaSBpCA(1,i)=alcSBpCApu;
OrdenadaSBpAG(1,i)=alcSBpAGpu;
OrdenadaSBpBG(1,i)=alcSBpBGpu;
OrdenadaSBpCG(1,i)=alcSBpCGpu;
OrdenadaRAmAB(1,i)=alcRAmABpu;
OrdenadaRAmBC(1,i)=alcRAmBCpu;
OrdenadaRAmCA(1,i)=alcRAmCApu;
OrdenadaRAmAG(1,i)=alcRAmAGpu;
OrdenadaRAmBG(1,i)=alcRAmBGpu;
OrdenadaRAmCG(1,i)=alcRAmCGpu;
OrdenadaRBpAB(1,i)=alcRBpABpu;
OrdenadaRBpBC(1,i)=alcRBpBCpu;
OrdenadaRBpCA(1,i)=alcRBpCApu;
OrdenadaRBpAG(1,i)=alcRBpAGpu;
OrdenadaRBpBG(1,i)=alcRBpBGpu;
OrdenadaRBpCG(1,i)=alcRBpCGpu;
OrdenadaSAmABa(1,i)=alcSAmABpua;
OrdenadaSAmBCa(1,i)=alcSAmBCpua;
OrdenadaSAmCAa(1,i)=alcSAmCApua;
OrdenadaSAmAGa(1,i)=alcSAmAGpua;
OrdenadaSAmBGa(1,i)=alcSAmBGpua;
OrdenadaSAmCGa(1,i)=alcSAmCGpua;
OrdenadaSBpABa(1,i)=alcSBpABpua;
OrdenadaSBpBCa(1,i)=alcSBpBCpua;
OrdenadaSBpCAa(1,i)=alcSBpCApua;
OrdenadaSBpAGa(1,i)=alcSBpAGpua;
OrdenadaSBpBGa(1,i)=alcSBpBGpua;
OrdenadaSBpCGa(1,i)=alcSBpCGpua;
OrdenadaRAmABa(1,i)=alcRAmABpua;
OrdenadaRAmBCa(1,i)=alcRAmBCpua;
OrdenadaRAmCAa(1,i)=alcRAmCApua;
OrdenadaRAmAGa(1,i)=alcRAmAGpua;
OrdenadaRAmBGa(1,i)=alcRAmBGpua;
OrdenadaRAmCGa(1,i)=alcRAmCGpua;
OrdenadaRBpABa(1,i)=alcRBpABpua;
OrdenadaRBpBCa(1,i)=alcRBpBCpua;
OrdenadaRBpCAa(1,i)=alcRBpCApua;
OrdenadaRBpAGa(1,i)=alcRBpAGpua;
OrdenadaRBpBGa(1,i)=alcRBpBGpua;
OrdenadaRBpCGa(1,i)=alcRBpCGpua;
%***************************************
%Para salida 3
%***************************************
vectorVVSApol(1,i)=VVSApol;
vectorIISAmop(1,i)=IISAmop;
vectorAngDirSAm(1,i)=AngDirSAm;
vectorVVRApol(1,i)=VVRApol;
vectorIIRAmop(1,i)=IIRAmop;
vectorAngDirRAm(1,i)=AngDirRAm;
vectorVVSBpol(1,i)=VVSBpol;
vectorIISBpop(1,i)=IISBpop;
vectorAngDirSBp(1,i)=AngDirSBp;
vectorVVRBpol(1,i)=VVRBpol;
vectorIIRBpop(1,i)=IIRBpop;
vectorAngDirRBp(1,i)=AngDirRBp;
%***************************************
% Para Salida 4 R-X
%***************************************
vectorSm_ZAB(1,i)=SAm_ZZAB;
vectorSm_ZBC(1,i)=SAm_ZZBC;
vectorSm_ZCA(1,i)=SAm_ZZCA;
vectorSm_ZAG(1,i)=SAm_ZZAG;
vectorSm_ZBG(1,i)=SAm_ZZBG;
vectorSm_ZCG(1,i)=SAm_ZZCG;
vectorRm_ZAB(1,i)=RAm_ZZAB;
vectorRm_ZBC(1,i)=RAm_ZZBC;
vectorRm_ZCA(1,i)=RAm_ZZCA;
vectorRm_ZAG(1,i)=RAm_ZZAG;
vectorRm_ZBG(1,i)=RAm_ZZBG;
vectorRm_ZCG(1,i)=RAm_ZZCG;
vectorSp_ZAB(1,i)=SBp_ZZAB;
vectorSp_ZBC(1,i)=SBp_ZZBC;
vectorSp_ZCA(1,i)=SBp_ZZCA;
vectorSp_ZAG(1,i)=SBp_ZZAG;
vectorSp_ZBG(1,i)=SBp_ZZBG;
vectorSp_ZCG(1,i)=SBp_ZZCG;
vectorRp_ZAB(1,i)=RBp_ZZAB;
vectorRp_ZBC(1,i)=RBp_ZZBC;
vectorRp_ZCA(1,i)=RBp_ZZCA;
vectorRp_ZAG(1,i)=RBp_ZZAG;
vectorRp_ZBG(1,i)=RBp_ZZBG;
vectorRp_ZCG(1,i)=RBp_ZZCG;
%***************************************
i=i+1;
m=m+0.01;
end %(while m)
% ----------------------------------------------------------
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
% ----------------------------------------------------------
%***************************************
% Titulo para las Salidas 1,2,4
%***************************************
cadenamc1=sprintf('PG128Mca,%s,%3.0fkV,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,Kr0=%-5.2f,k0=%-4.1f\n',...
comp,UA,UB,kE,L,Sccs,Sccr,Kr0,k0);
if imptransfer==1
cadenamT=sprintf('zsr1=%-6.2fi,Ksr0=%-4.2f',imag(zsr1)*zbaseA,Ksr0);
cadenamc1=sprintf('%s,%s',cadenamc1,cadenamT);
end
cadenamc2=sprintf('%s,Rf=%-3.0f',cfalta,Rf);
if falta==1|falta==2|falta==3|falta==4|falta==7|falta==8|falta==9|falta==13|falta==15
cadenamc=sprintf('%s,%s',cadenamc1,cadenamc2);
elseif falta==5|falta==6|falta==10|falta==11|falta==12|falta==14|falta==16
cadenamc3=sprintf('Rg=%-3.0f',Rg);
cadenamc=sprintf('%s,%s,%s',cadenamc1,cadenamc2,cadenamc3);
end
%***************************************
% Titulo para las gráficas 67N
%***************************************
dircadenamc1=sprintf('PG128Mca,67N,%3.0fkV,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,Kr0=%-5.2f,k0=%-4.1f\n',...
UA,UB,kE,L,Sccs,Sccr,Kr0,k0);
if imptransfer==1
dircadenamT=sprintf('zsr1=%-6.2fi,Ksr0=%-4.2f',imag(zsr1)*zbaseA,Ksr0);
dircadenamc1=sprintf('%s,%s',dircadenamc1,dircadenamT);
end
dircadenamc2=sprintf('%s,Rf=%-3.0f',cfalta,Rf);
if falta==1|falta==2|falta==3|falta==4|falta==7|falta==8|falta==9|falta==13|falta==15
dircadenamc=sprintf('%s,%s',dircadenamc1,dircadenamc2);
elseif falta==5|falta==6|falta==10|falta==11|falta==12|falta==14|falta==16
dircadenamc3=sprintf('Rg=%-3.0f',Rg);
dircadenamc=sprintf('%s,%s,%s',dircadenamc1,dircadenamc2,dircadenamc3);
end
%*************************************
% ----------------------------------------------------------------------------------
% SALIDA GRAFICA
% ----------------------------------------------------------------------------------
variassalidas=1;
while variassalidas==1
fprintf('Salida Grafica 1 : Analisis SAm, Todas las Unidades (pu)\n');
fprintf(' Analisis RBp, Todas las Unidades (pu)\n');
fprintf(' Analisis RAm, Todas las Unidades (pu)\n');
fprintf(' Analisis RBp, Todas las Unidades (pu)\n');
fprintf('Salida Grafica 2 : Analisis SAm, Todas las Unidades (pua)\n');
fprintf(' : Analisis SBp, Todas las Unidades (pua)\n');
fprintf(' Analisis RAm, Todas las Unidades (pua)\n');
fprintf(' Analisis RBp, Todas las Unidades (pua)\n');
fprintf('Salida Grafica 3 : Analisis SAm, 67N\n');
fprintf(' : Analisis RAm, 67N\n');
fprintf(' : Analisis SBp, 67N\n');
fprintf(' : Analisis RBp, 67N\n');
fprintf('Salida Grafica 4: Analisis SAm, Todas las Unidades R-X\n');
fprintf(' : Analisis RAm, Todas las Unidades R-X\n');
fprintf(' : Analisis SBp, Todas las Unidades R-X\n');
fprintf(' : Analisis RBp, Todas las Unidades R-X\n');
salida=input('Salida: 1,2,3,4:');
switch(salida)
case{1}, % **************** Salida 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
extremo=input('SAm=1,SBp=2,RAm=3,RBp=4 : ');
switch(extremo)
case{1},
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
hold on
h=plot(Abscisa0,Ordenadadi);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,OrdenadaSAmAB,'r',Abscisa,OrdenadaSAmBC,'g',Abscisa,OrdenadaSAmCA,'b',...
Abscisa,OrdenadaSAmAG,'m',Abscisa,OrdenadaSAmBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSAmCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadylabel1)
legend('Ref','Lim','AB','BC','CA','AG','BG','CG')
case{2},
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
hold on
h=plot(Abscisa0,Ordenadadi);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,OrdenadaSBpAB,'r',Abscisa,OrdenadaSBpBC,'g',Abscisa,OrdenadaSBpCA,'b',...
Abscisa,OrdenadaSBpAG,'m',Abscisa,OrdenadaSBpBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSBpCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadylabel2)
legend('Ref','Lim','AB','BC','CA','AG','BG','CG')
case{3},
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
hold on
h=plot(Abscisa0,Ordenadadd);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,OrdenadaRAmAB,'r',Abscisa,OrdenadaRAmBC,'g',Abscisa,OrdenadaRAmCA,'b',...
Abscisa,OrdenadaRAmAG,'m',Abscisa,OrdenadaRAmBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRAmCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadylabel3)
legend('Ref','Lim','AB','BC','CA','AG','BG','CG')
case{4},
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
hold on
h=plot(Abscisa0,Ordenadadd);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,OrdenadaRBpAB,'r',Abscisa,OrdenadaRBpBC,'g',Abscisa,OrdenadaRBpCA,'b',...
Abscisa,OrdenadaRBpAG,'m',Abscisa,OrdenadaRBpBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRBpCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadylabel4)
legend('Ref','Lim','AB','BC','CA','AG','BG','CG')
end
case(2) % pua, Todas las Unidades
extremo=input('SAm=1,SBp=2,RAm=3,RBp=4 : ');
switch(extremo)
case{1},
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenada45(cont)=mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaSAmABa,'r',Abscisa,OrdenadaSAmBCa,'g',Abscisa,OrdenadaSAmCAa,'b',...
Abscisa,OrdenadaSAmAGa,'m',Abscisa,OrdenadaSAmBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSAmCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadylabel1a)
legend('Ref','AB','BC','CA','AG','BG','CG','z1')
case{2},
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenada45(cont)=mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaSBpABa,'r',Abscisa,OrdenadaSBpBCa,'g',Abscisa,OrdenadaSBpCAa,'b',...
Abscisa,OrdenadaSBpAGa,'m',Abscisa,OrdenadaSBpBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSBpCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadylabel2a)
legend('Ref','AB','BC','CA','AG','BG','CG','z1')
case{3},
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenada45(cont)=1-mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaRAmABa,'r',Abscisa,OrdenadaRAmBCa,'g',Abscisa,OrdenadaRAmCAa,'b',...
Abscisa,OrdenadaRAmAGa,'m',Abscisa,OrdenadaRAmBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRAmCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadylabel3a)
legend('Ref','AB','BC','CA','AG','BG','CG','z1')
case{4},
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenada45(cont)=1-mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaRBpABa,'r',Abscisa,OrdenadaRBpBCa,'g',Abscisa,OrdenadaRBpCAa,'b',...
Abscisa,OrdenadaRBpAGa,'m',Abscisa,OrdenadaRBpBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRBpCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadylabel4a)
legend('Ref','AB','BC','CA','AG','BG','CG','z1')
end
%**************************************
case{3} % Direccionales
extr=input('Extremo: SAm=1, RAm=2, SBp=3, RBp=4: ');
switch(extr)
case{1}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSAm,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirSAm(grados)')
title(dircadenamc)
rayar(1);
legend('AngDirSAm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVVSApol,'b');
grid
%title(dircadenamc)
xlabel('m=p'),ylabel('modulo VSApol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIISAmop,'b');
grid
xlabel('m=p'),ylabel('modulo ISAmop(kA)')
case{2}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRAm,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirRAm(grados)')
title(dircadenamc)
rayar(1);
legend('AngDirRAm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVVRApol,'b');
grid
%title(dircadenamc)
xlabel('m=p'),ylabel('modulo VRApol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIIRAmop,'b');
grid
xlabel('m=p'),ylabel('modulo IRAmop(kA)')
case{3}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSBp,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirSBp(grados)')
title(dircadenamc)
rayar(1);
legend('AngDirSBp','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVVSBpol,'b');
grid
%title(dircadenamc)
xlabel('m=p'),ylabel('modulo VSBpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIISBpop,'b');
grid
xlabel('m=p'),ylabel('modulo ISBpop(kA)')
case{4}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRBp,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirRBp(grados)')
title(dircadenamc)
rayar(1);
legend('AngDirRBp','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVVRBpol,'b');
grid
%title(dircadenamc)
xlabel('m=p'),ylabel('modulo VRBpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIIRBpop,'b');
grid
xlabel('m=p'),ylabel('modulo IRBpop(kA)')
end % end del switch extr
%*************************************
%*************************************
case{4}, % ********* Salida 4 ************
extr=input('Extremo: SAm=1, RAm=2, SBp=3, RBp=4: ');
switch(extr)
case{1}
proteccion='Sm';
pzona=0.8;
hold on
plot19(L,zlA1u,pzona,proteccion)
hold off
hold on
% Colores
a=plot(vectorSm_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorSm_ZBC,'m')
plot(vectorSm_ZCA,'c')
plot(vectorSm_ZAG,'r')
plot(vectorSm_ZBG,'k')
% Colores
c=plot(vectorSm_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
[subAbscisa,subvectorSm_ZAB]=subconjunto21(Abscisa,vectorSm_ZAB);
subAbscisa(:,1)=[];
subvectorSm_ZAB(:,1)=[];
hold on
b=plot(subvectorSm_ZAB,'*');
set(b,'Color',[1 0.60 0])
plot20(L,zlA1u,subvectorSm_ZAB,'naranja')
hold off
[subAbscisa,subvectorSm_ZBC]=subconjunto21(Abscisa,vectorSm_ZBC);
subAbscisa(:,1)=[];
subvectorSm_ZBC(:,1)=[];
hold on
plot(subvectorSm_ZBC,'m*');
plot20(L,zlA1u,subvectorSm_ZBC,'m')
hold off
[subAbscisa,subvectorSm_ZCA]=subconjunto21(Abscisa,vectorSm_ZCA);
subAbscisa(:,1)=[];
subvectorSm_ZCA(:,1)=[];
hold on
plot(subvectorSm_ZCA,'c*');
plot20(L,zlA1u,subvectorSm_ZCA,'c')
hold off
[subAbscisa,subvectorSm_ZAG]=subconjunto21(Abscisa,vectorSm_ZAG);
subAbscisa(:,1)=[];
subvectorSm_ZAG(:,1)=[];
hold on
plot(subvectorSm_ZAG,'r*');
plot20(L,zlA1u,subvectorSm_ZAG,'r')
hold off
[subAbscisa,subvectorSm_ZBG]=subconjunto21(Abscisa,vectorSm_ZBG);
subAbscisa(:,1)=[];
subvectorSm_ZBG(:,1)=[];
hold on
plot(subvectorSm_ZBG,'k*');
plot20(L,zlA1u,subvectorSm_ZBG,'k')
hold off
[subAbscisa,subvectorSm_ZCG]=subconjunto21(Abscisa,vectorSm_ZCG);
subAbscisa(:,1)=[];
subvectorSm_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorSm_ZCG,'*');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zlA1u,subvectorSm_ZCG,'marron')
hold off
proteccionc='SAm';
RXcadenamc=sprintf('%s,Proteccion:%s',cadenamc,proteccionc);
title(RXcadenamc)
legend('Mho','Linea','m=0..1','AB','BC','CA','AG','BG','CG')
case{2}
proteccion='Rm';
pzona=0.8;
hold on
plot19(L,zlA1u,pzona,proteccion)
hold off
hold on
% Colores
a=plot(vectorRm_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorRm_ZBC,'m')
plot(vectorRm_ZCA,'c')
plot(vectorRm_ZAG,'r')
plot(vectorRm_ZBG,'k')
% Colores
c=plot(vectorRm_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
[subAbscisa,subvectorRm_ZAB]=subconjunto21(Abscisa,vectorRm_ZAB);
subAbscisa(:,11)=[];
subvectorRm_ZAB(:,11)=[];
invsubvectorRm_ZAB=subvectorRm_ZAB(1:1,10:-1:1);
hold on
b=plot(invsubvectorRm_ZAB,'*');
set(b,'Color',[1 0.60 0])
plot20(L,zlA1u,invsubvectorRm_ZAB,'naranja')
hold off
[subAbscisa,subvectorRm_ZBC]=subconjunto21(Abscisa,vectorRm_ZBC);
subAbscisa(:,11)=[];
subvectorRm_ZBC(:,11)=[];
invsubvectorRm_ZBC=subvectorRm_ZBC(1:1,10:-1:1);
hold on
plot(subvectorRm_ZBC,'m*');
plot20(L,zlA1u,invsubvectorRm_ZBC,'m')
hold off
[subAbscisa,subvectorRm_ZCA]=subconjunto21(Abscisa,vectorRm_ZCA);
subAbscisa(:,11)=[];
subvectorRm_ZCA(:,11)=[];
invsubvectorRm_ZCA=subvectorRm_ZCA(1:1,10:-1:1);
hold on
plot(subvectorRm_ZCA,'c*');
plot20(L,zlA1u,invsubvectorRm_ZCA,'c')
hold off
[subAbscisa,subvectorRm_ZAG]=subconjunto21(Abscisa,vectorRm_ZAG);
subAbscisa(:,11)=[];
subvectorRm_ZAG(:,11)=[];
invsubvectorRm_ZAG=subvectorRm_ZAG(1:1,10:-1:1);
hold on
plot(subvectorRm_ZAG,'r*');
plot20(L,zlA1u,invsubvectorRm_ZAG,'r')
hold off
[subAbscisa,subvectorRm_ZBG]=subconjunto21(Abscisa,vectorRm_ZBG);
subAbscisa(:,11)=[];
subvectorRm_ZBG(:,11)=[];
invsubvectorRm_ZBG=subvectorRm_ZBG(1:1,10:-1:1);
hold on
plot(subvectorRm_ZBG,'k*');
plot20(L,zlA1u,invsubvectorRm_ZBG,'k')
hold off
[subAbscisa,subvectorRm_ZCG]=subconjunto21(Abscisa,vectorRm_ZCG);
subAbscisa(:,11)=[];
subvectorRm_ZCG(:,11)=[];
invsubvectorRm_ZCG=subvectorRm_ZCG(1:1,10:-1:1);
hold on
% Colores
c=plot(invsubvectorRm_ZCG,'*');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zlA1u,invsubvectorRm_ZCG,'marron')
hold off
proteccionc='RAm';
RXcadenamc=sprintf('%s,Proteccion:%s',cadenamc,proteccionc);
title(RXcadenamc)
legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG')
case{3}
proteccion='Sp';
pzona=0.8;
hold on
plot19(L,zlB1u,pzona,proteccion)
hold off
hold on
% Colores
a=plot(vectorSp_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorSp_ZBC,'m')
plot(vectorSp_ZCA,'c')
plot(vectorSp_ZAG,'r')
plot(vectorSp_ZBG,'k')
% Colores
c=plot(vectorSp_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
[subAbscisa,subvectorSp_ZAB]=subconjunto21(Abscisa,vectorSp_ZAB);
subAbscisa(:,1)=[];
subvectorSp_ZAB(:,1)=[];
hold on
b=plot(subvectorSp_ZAB,'*');
set(b,'Color',[1 0.60 0])
plot20(L,zlB1u,subvectorSp_ZAB,'naranja')
hold off
[subAbscisa,subvectorSp_ZBC]=subconjunto21(Abscisa,vectorSp_ZBC);
subAbscisa(:,1)=[];
subvectorSp_ZBC(:,1)=[];
hold on
plot(subvectorSp_ZBC,'m*');
plot20(L,zlB1u,subvectorSp_ZBC,'m')
hold off
[subAbscisa,subvectorSp_ZCA]=subconjunto21(Abscisa,vectorSp_ZCA);
subAbscisa(:,1)=[];
subvectorSp_ZCA(:,1)=[];
hold on
plot(subvectorSp_ZCA,'c*');
plot20(L,zlB1u,subvectorSp_ZCA,'c')
hold off
[subAbscisa,subvectorSp_ZAG]=subconjunto21(Abscisa,vectorSp_ZAG);
subAbscisa(:,1)=[];
subvectorSp_ZAG(:,1)=[];
hold on
plot(subvectorSp_ZAG,'r*');
plot20(L,zlB1u,subvectorSp_ZAG,'r')
hold off
[subAbscisa,subvectorSp_ZBG]=subconjunto21(Abscisa,vectorSp_ZBG);
subAbscisa(:,1)=[];
subvectorSp_ZBG(:,1)=[];
hold on
plot(subvectorSp_ZBG,'k*');
plot20(L,zlB1u,subvectorSp_ZBG,'k')
hold off
[subAbscisa,subvectorSp_ZCG]=subconjunto21(Abscisa,vectorSp_ZCG);
subAbscisa(:,1)=[];
subvectorSp_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorSp_ZCG,'*');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zlB1u,subvectorSp_ZCG,'marron')
hold off
proteccionc='SBp';
RXcadenamc=sprintf('%s,Proteccion:%s',cadenamc,proteccionc);
title(RXcadenamc)
legend('Mho','Linea','m=0..1','AB','BC','CA','AG','BG','CG')
case{4}
proteccion='Rp';
pzona=0.8;
hold on
plot19(L,zlB1u,pzona,proteccion)
hold off
hold on
% Colores
a=plot(vectorRp_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorRp_ZBC,'m')
plot(vectorRp_ZCA,'c')
plot(vectorRp_ZAG,'r')
plot(vectorRp_ZBG,'k')
% Colores
c=plot(vectorRp_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
[subAbscisa,subvectorRp_ZAB]=subconjunto21(Abscisa,vectorRp_ZAB);
subAbscisa(:,11)=[];
subvectorRp_ZAB(:,11)=[];
invsubvectorRp_ZAB=subvectorRp_ZAB(1:1,10:-1:1);
hold on
b=plot(invsubvectorRp_ZAB,'*');
set(b,'Color',[1 0.60 0])
plot20(L,zlB1u,invsubvectorRp_ZAB,'naranja')
hold off
[subAbscisa,subvectorRp_ZBC]=subconjunto21(Abscisa,vectorRp_ZBC);
subAbscisa(:,11)=[];
subvectorRp_ZBC(:,11)=[];
invsubvectorRp_ZBC=subvectorRp_ZBC(1:1,10:-1:1);
hold on
plot(subvectorRp_ZBC,'m*');
plot20(L,zlB1u,invsubvectorRp_ZBC,'m')
hold off
[subAbscisa,subvectorRp_ZCA]=subconjunto21(Abscisa,vectorRp_ZCA);
subAbscisa(:,11)=[];
subvectorRp_ZCA(:,11)=[];
invsubvectorRp_ZCA=subvectorRp_ZCA(1:1,10:-1:1);
hold on
plot(subvectorRp_ZCA,'c*');
plot20(L,zlB1u,invsubvectorRp_ZCA,'c')
hold off
[subAbscisa,subvectorRp_ZAG]=subconjunto21(Abscisa,vectorRp_ZAG);
subAbscisa(:,11)=[];
subvectorRp_ZAG(:,11)=[];
invsubvectorRp_ZAG=subvectorRp_ZAG(1:1,10:-1:1);
hold on
plot(subvectorRp_ZAG,'r*');
plot20(L,zlB1u,invsubvectorRp_ZAG,'r')
hold off
[subAbscisa,subvectorRp_ZBG]=subconjunto21(Abscisa,vectorRp_ZBG);
subAbscisa(:,11)=[];
subvectorRp_ZBG(:,11)=[];
invsubvectorRp_ZBG=subvectorRp_ZBG(1:1,10:-1:1);
hold on
plot(subvectorRp_ZBG,'k*');
plot20(L,zlB1u,invsubvectorRp_ZBG,'k')
hold off
[subAbscisa,subvectorRp_ZCG]=subconjunto21(Abscisa,vectorRp_ZCG);
subAbscisa(:,11)=[];
subvectorRp_ZCG(:,11)=[];
invsubvectorRp_ZCG=subvectorRp_ZCG(1:1,10:-1:1);
hold on
% Colores
c=plot(invsubvectorRp_ZCG,'*');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zlB1u,invsubvectorRp_ZCG,'marron')
hold off
proteccionc='RBp';
RXcadenamc=sprintf('%s,Proteccion:%s',cadenamc,proteccionc);
title(RXcadenamc)
legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG')
end % end del switch extr de Salida 4
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');
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
