% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG126Mc\n');
fprintf('Esquema Linea Doble SR\n');
fprintf('Unominal kV: 400,220,132,66\n');
fprintf('Eleccion: kE,L,Sccs,Sccr,Kr0,k0)\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('Compensacion 21, Si, No\n');
fprintf('Falta Intercircuito entre Lineas de igual Tension\n');
fprintf('m=p=0...1\n');
fprintf('Barra m=1, Barra p=2\n');
%fprintf('Salida Grafica 1: Analisis Sm, Todas las Unidades (pu)\n');
%fprintf(' : Analisis Sp, Todas las Unidades (pu)\n');
%fprintf(' : Analisis Rm, Todas las Unidades (pu)\n');
%fprintf(' : Analisis Rp, Todas las Unidades (pu)\n');
%fprintf('Salida Grafica 2: Analisis Sm, Todas las Unidades (pua)\n');
%fprintf(' : Analisis Sp, Todas las Unidades (pua)\n');
%fprintf(' Analisis Rm, Todas las Unidades (pua)\n');
%fprintf(' Analisis Rp, Todas las Unidades (pua)\n');
%fprintf('Salida Grafica 3: Analisis Sm, 67N\n');
%fprintf(' : Analisis Rm, 67N\n');
%fprintf(' : Analisis Sp, 67N\n');
%fprintf(' : Analisis Rp, 67N\n');
%fprintf('Salida Grafica 4: Analisis Sm, Todas las Unidades R-X\n');
%fprintf(' : Analisis Rm, Todas las Unidades R-X\n');
%fprintf(' : Analisis Sp, Todas las Unidades R-X\n');
%fprintf(' : Analisis Rp, Todas las Unidades R-X\n');
fprintf('\n');
% ----------------------------------------------------------------------------------
%PG126Mc
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=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
% ----------------------------------------------------------------------------------
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;
% ----------------------------------------------------------------------------------
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
% ----------------------------------------------------------------------------------
fprintf('\n');
i=1;
m=0.01;
while m<1;
p=m;
%
switch(imptransfer)
case{1},
[Minc]=subMincT(1);
[z0]=subz0T(zs0,zr0,zsr0,zl0u,zlM0u,L,m,p);
[z1]=subz1T(zs1,zr1,zsr1,zl1u,L,m,p);
[z2]=subz2T(zs2,zr2,zsr2,zl2u,L,m,p);
case(2)
[Minc]=subMinc(1);
[z0]=subz0(zs0,zr0,zl0u,zlM0u,L,m,p);
[z1]=subz1(zs1,zr1,zl1u,L,m,p);
[z2]=subz2(zs2,zr2,zl2u,L,m,p);
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]=...
maMmp(Z0,Z1,Z2,E,falta,Rf,Rg);
% IBUS(0,1,2)
I0=[-Im0;-Ip0;0;0];
I1=[-Im1;-Ip1;0;0];
I2=[-Im2;-Ip2;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(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(3); Vr0=V0(4);
Vs1=V1(3); Vr1=V1(4);
Vs2=V2(3); Vr2=V2(4);
Ism0=Ie0(3); Irm0=Ie0(5); Isp0=Ie0(4); Irp0=Ie0(6);
Ism1=Ie1(3); Irm1=Ie1(5); Isp1=Ie1(4); Irp1=Ie1(6);
Ism2=Ie2(3); Irm2=Ie2(5); Isp2=Ie2(4); Irp2=Ie2(6);
% 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);
[Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2);
[Ispa,Ispb,Ispc]=abc(Isp0,Isp1,Isp2);
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2);
[Irpa,Irpb,Irpc]=abc(Irp0,Irp1,Irp2);
%
% PROTECCION DE DISTANCIA (21)
switch(compens21)
case{1},
% COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k,Km,Isp0);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k,Km,Ism0);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k,Km,Irp0);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k,Km,Irm0);
case{2},
% NO COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
PDist(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
PDist(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k);
end
% CALCULO DE LAS MAGNITUDES DE ALIMETACION DE LA PROTECCION
% DIRECCIONAL DE TIERRA EN LOS EXTREMOS DE LINEA
[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0); % kV, kA, grados
[Vspol,Ispop,AngDirSp]=PDire(Vs0,Isp0);
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0);
[Vrpol,Irpop,AngDirRp]=PDire(Vr0,Irp0);
% para Salida 1 y 2
alcSmAB=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
alcSmBC=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcSmCA=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
alcSmAG=abs(Sm_ZAG)/cos(angle(Sm_ZAG)-angle(zl1u));
alcSmBG=abs(Sm_ZBG)/cos(angle(Sm_ZBG)-angle(zl1u));
alcSmCG=abs(Sm_ZCG)/cos(angle(Sm_ZCG)-angle(zl1u));
alcSpAB=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
alcSpBC=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcSpCA=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
alcSpAG=abs(Sp_ZAG)/cos(angle(Sp_ZAG)-angle(zl1u));
alcSpBG=abs(Sp_ZBG)/cos(angle(Sp_ZBG)-angle(zl1u));
alcSpCG=abs(Sp_ZCG)/cos(angle(Sp_ZCG)-angle(zl1u));
alcRmAB=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));
alcRmBC=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
alcRmCA=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));
alcRmAG=abs(Rm_ZAG)/cos(angle(Rm_ZAG)-angle(zl1u));
alcRmBG=abs(Rm_ZBG)/cos(angle(Rm_ZBG)-angle(zl1u));
alcRmCG=abs(Rm_ZCG)/cos(angle(Rm_ZCG)-angle(zl1u));
alcRpAB=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
alcRpBC=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
alcRpCA=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
alcRpAG=abs(Rp_ZAG)/cos(angle(Rp_ZAG)-angle(zl1u));
alcRpBG=abs(Rp_ZBG)/cos(angle(Rp_ZBG)-angle(zl1u));
alcRpCG=abs(Rp_ZCG)/cos(angle(Rp_ZCG)-angle(zl1u));
% para Salida 1
alcSmABpu=alcSmAB/abs(zl1u*m*L);
alcSmBCpu=alcSmBC/abs(zl1u*m*L);
alcSmCApu=alcSmCA/abs(zl1u*m*L);
alcSmAGpu=alcSmAG/abs(zl1u*m*L);
alcSmBGpu=alcSmBG/abs(zl1u*m*L);
alcSmCGpu=alcSmCG/abs(zl1u*m*L);
alcSpABpu=alcSpAB/abs(zl1u*m*L);
alcSpBCpu=alcSpBC/abs(zl1u*m*L);
alcSpCApu=alcSpCA/abs(zl1u*m*L);
alcSpAGpu=alcSpAG/abs(zl1u*m*L);
alcSpBGpu=alcSpBG/abs(zl1u*m*L);
alcSpCGpu=alcSpCG/abs(zl1u*m*L);
alcRmABpu=alcRmAB/abs(zl1u*(1-m)*L);
alcRmBCpu=alcRmBC/abs(zl1u*(1-m)*L);
alcRmCApu=alcRmCA/abs(zl1u*(1-m)*L);
alcRmAGpu=alcRmAG/abs(zl1u*(1-m)*L);
alcRmBGpu=alcRmBG/abs(zl1u*(1-m)*L);
alcRmCGpu=alcRmCG/abs(zl1u*(1-m)*L);
alcRpABpu=alcRpAB/abs(zl1u*(1-m)*L);
alcRpBCpu=alcRpBC/abs(zl1u*(1-m)*L);
alcRpCApu=alcRpCA/abs(zl1u*(1-m)*L);
alcRpAGpu=alcRpAG/abs(zl1u*(1-m)*L);
alcRpBGpu=alcRpBG/abs(zl1u*(1-m)*L);
alcRpCGpu=alcRpCG/abs(zl1u*(1-m)*L);
cadylabel1='alcSmpu';
cadylabel2='alcSppu';
cadylabel3='alcRmpu';
cadylabel4='alcRppu';
% para Salida 2
alcSmABpua=alcSmAB/abs(zl1u*L);
alcSmBCpua=alcSmBC/abs(zl1u*L);
alcSmCApua=alcSmCA/abs(zl1u*L);
alcSmAGpua=alcSmAG/abs(zl1u*L);
alcSmBGpua=alcSmBG/abs(zl1u*L);
alcSmCGpua=alcSmCG/abs(zl1u*L);
alcSpABpua=alcSpAB/abs(zl1u*L);
alcSpBCpua=alcSpBC/abs(zl1u*L);
alcSpCApua=alcSpCA/abs(zl1u*L);
alcSpAGpua=alcSpAG/abs(zl1u*L);
alcSpBGpua=alcSpBG/abs(zl1u*L);
alcSpCGpua=alcSpCG/abs(zl1u*L);
alcRmABpua=alcRmAB/abs(zl1u*L);
alcRmBCpua=alcRmBC/abs(zl1u*L);
alcRmCApua=alcRmCA/abs(zl1u*L);
alcRmAGpua=alcRmAG/abs(zl1u*L);
alcRmBGpua=alcRmBG/abs(zl1u*L);
alcRmCGpua=alcRmCG/abs(zl1u*L);
alcRpABpua=alcRpAB/abs(zl1u*L);
alcRpBCpua=alcRpBC/abs(zl1u*L);
alcRpCApua=alcRpCA/abs(zl1u*L);
alcRpAGpua=alcRpAG/abs(zl1u*L);
alcRpBGpua=alcRpBG/abs(zl1u*L);
alcRpCGpua=alcRpCG/abs(zl1u*L);
cadylabel1a='alcSmpua';
cadylabel2a='alcSppua';
cadylabel3a='alcRmpua';
cadylabel4a='alcRppua';
Abscisa(i)=m;
% para Salida 1
OrdenadaSmAB(1,i)=alcSmABpu;
OrdenadaSmBC(1,i)=alcSmBCpu;
OrdenadaSmCA(1,i)=alcSmCApu;
OrdenadaSmAG(1,i)=alcSmAGpu;
OrdenadaSmBG(1,i)=alcSmBGpu;
OrdenadaSmCG(1,i)=alcSmCGpu;
OrdenadaSpAB(1,i)=alcSpABpu;
OrdenadaSpBC(1,i)=alcSpBCpu;
OrdenadaSpCA(1,i)=alcSpCApu;
OrdenadaSpAG(1,i)=alcSpAGpu;
OrdenadaSpBG(1,i)=alcSpBGpu;
OrdenadaSpCG(1,i)=alcSpCGpu;
OrdenadaRmAB(1,i)=alcRmABpu;
OrdenadaRmBC(1,i)=alcRmBCpu;
OrdenadaRmCA(1,i)=alcRmCApu;
OrdenadaRmAG(1,i)=alcRmAGpu;
OrdenadaRmBG(1,i)=alcRmBGpu;
OrdenadaRmCG(1,i)=alcRmCGpu;
OrdenadaRpAB(1,i)=alcRpABpu;
OrdenadaRpBC(1,i)=alcRpBCpu;
OrdenadaRpCA(1,i)=alcRpCApu;
OrdenadaRpAG(1,i)=alcRpAGpu;
OrdenadaRpBG(1,i)=alcRpBGpu;
OrdenadaRpCG(1,i)=alcRpCGpu;
% para Salida 2
OrdenadaSmABa(1,i)=alcSmABpua;
OrdenadaSmBCa(1,i)=alcSmBCpua;
OrdenadaSmCAa(1,i)=alcSmCApua;
OrdenadaSmAGa(1,i)=alcSmAGpua;
OrdenadaSmBGa(1,i)=alcSmBGpua;
OrdenadaSmCGa(1,i)=alcSmCGpua;
OrdenadaSpABa(1,i)=alcSpABpua;
OrdenadaSpBCa(1,i)=alcSpBCpua;
OrdenadaSpCAa(1,i)=alcSpCApua;
OrdenadaSpAGa(1,i)=alcSpAGpua;
OrdenadaSpBGa(1,i)=alcSpBGpua;
OrdenadaSpCGa(1,i)=alcSpCGpua;
OrdenadaRmABa(1,i)=alcRmABpua;
OrdenadaRmBCa(1,i)=alcRmBCpua;
OrdenadaRmCAa(1,i)=alcRmCApua;
OrdenadaRmAGa(1,i)=alcRmAGpua;
OrdenadaRmBGa(1,i)=alcRmBGpua;
OrdenadaRmCGa(1,i)=alcRmCGpua;
OrdenadaRpABa(1,i)=alcRpABpua;
OrdenadaRpBCa(1,i)=alcRpBCpua;
OrdenadaRpCAa(1,i)=alcRpCApua;
OrdenadaRpAGa(1,i)=alcRpAGpua;
OrdenadaRpBGa(1,i)=alcRpBGpua;
OrdenadaRpCGa(1,i)=alcRpCGpua;
%***************************************
%Para salida 4
%***************************************
vectorVspol(1,i)=Vspol;
vectorIsmop(1,i)=Ismop;
vectorAngDirSm(1,i)=AngDirSm;
vectorVrpol(1,i)=Vrpol;
vectorIrmop(1,i)=Irmop;
vectorAngDirRm(1,i)=AngDirRm;
vectorVspol(1,i)=Vspol;
vectorIspop(1,i)=Ispop;
vectorAngDirSp(1,i)=AngDirSp;
vectorVrpol(1,i)=Vrpol;
vectorIrpop(1,i)=Irpop;
vectorAngDirRp(1,i)=AngDirRp;
%***************************************
% Para Salida 5 R-X
%***************************************
vectorSm_ZAB(1,i)=Sm_ZAB;
vectorSm_ZBC(1,i)=Sm_ZBC;
vectorSm_ZCA(1,i)=Sm_ZCA;
vectorSm_ZAG(1,i)=Sm_ZAG;
vectorSm_ZBG(1,i)=Sm_ZBG;
vectorSm_ZCG(1,i)=Sm_ZCG;
vectorRm_ZAB(1,i)=Rm_ZAB;
vectorRm_ZBC(1,i)=Rm_ZBC;
vectorRm_ZCA(1,i)=Rm_ZCA;
vectorRm_ZAG(1,i)=Rm_ZAG;
vectorRm_ZBG(1,i)=Rm_ZBG;
vectorRm_ZCG(1,i)=Rm_ZCG;
vectorSp_ZAB(1,i)=Sp_ZAB;
vectorSp_ZBC(1,i)=Sp_ZBC;
vectorSp_ZCA(1,i)=Sp_ZCA;
vectorSp_ZAG(1,i)=Sp_ZAG;
vectorSp_ZBG(1,i)=Sp_ZBG;
vectorSp_ZCG(1,i)=Sp_ZCG;
vectorRp_ZAB(1,i)=Rp_ZAB;
vectorRp_ZBC(1,i)=Rp_ZBC;
vectorRp_ZCA(1,i)=Rp_ZCA;
vectorRp_ZAG(1,i)=Rp_ZAG;
vectorRp_ZBG(1,i)=Rp_ZBG;
vectorRp_ZCG(1,i)=Rp_ZCG;
%***************************************
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('PG126Mc,%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
cadenamT=sprintf('zsr1=%-6.2fi,Ksr0=%-4.2f',imag(zsr1),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 Salida 3, gráficas 67N
%***************************************
dircadenamc1=sprintf('PG126Mc,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
dircadenamT=sprintf('zsr1=%-6.2fi,Ksr0=%-4.2f',imag(zsr1),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 Sm, Todas las Unidades (pu)\n');
fprintf(' : Analisis Sp, Todas las Unidades (pu)\n');
fprintf(' : Analisis Rm, Todas las Unidades (pu)\n');
fprintf(' : Analisis Rp, Todas las Unidades (pu)\n');
fprintf('Salida Grafica 2: Analisis Sm, Todas las Unidades (pua)\n');
fprintf(' : Analisis Sp, Todas las Unidades (pua)\n');
fprintf(' Analisis Rm, Todas las Unidades (pua)\n');
fprintf(' Analisis Rp, Todas las Unidades (pua)\n');
fprintf('Salida Grafica 3: Analisis Sm, 67N\n');
fprintf(' : Analisis Rm, 67N\n');
fprintf(' : Analisis Sp, 67N\n');
fprintf(' : Analisis Rp, 67N\n');
fprintf('Salida Grafica 4: Analisis Sm, Todas las Unidades R-X\n');
fprintf(' : Analisis Rm, Todas las Unidades R-X\n');
fprintf(' : Analisis Sp, Todas las Unidades R-X\n');
fprintf(' : Analisis Rp, Todas las Unidades R-X\n');
fprintf('\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
extremo=input('Sm=1,Sp=2,Rm=3,Rs=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,OrdenadaSmAB,'r',Abscisa,OrdenadaSmBC,'g',Abscisa,OrdenadaSmCA,'b',...
Abscisa,OrdenadaSmAG,'m',Abscisa,OrdenadaSmBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSmCG);
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,OrdenadaSpAB,'r',Abscisa,OrdenadaSpBC,'g',Abscisa,OrdenadaSpCA,'b',...
Abscisa,OrdenadaSpAG,'m',Abscisa,OrdenadaSpBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSpCG);
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,OrdenadaRmAB,'r',Abscisa,OrdenadaRmBC,'g',Abscisa,OrdenadaRmCA,'b',...
Abscisa,OrdenadaRmAG,'m',Abscisa,OrdenadaRmBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRmCG);
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,OrdenadaRpAB,'r',Abscisa,OrdenadaRpBC,'g',Abscisa,OrdenadaRpCA,'b',...
Abscisa,OrdenadaRpAG,'m',Abscisa,OrdenadaRpBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRpCG);
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) % *********** Salida Grafica 2 (pua) ****************
extremo=input('Sm=1,Sp=2,Rm=3,Rp=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
% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaSmABa,'r',Abscisa,OrdenadaSmBCa,'g',Abscisa,OrdenadaSmCAa,'b',...
Abscisa,OrdenadaSmAGa,'m',Abscisa,OrdenadaSmBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSmCGa);
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 -2 2]);
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
% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaSpABa,'r',Abscisa,OrdenadaSpBCa,'g',Abscisa,OrdenadaSpCAa,'b',...
Abscisa,OrdenadaSpAGa,'m',Abscisa,OrdenadaSpBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSpCGa);
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 -2 2]);
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
% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaRmABa,'r',Abscisa,OrdenadaRmBCa,'g',Abscisa,OrdenadaRmCAa,'b',...
Abscisa,OrdenadaRmAGa,'m',Abscisa,OrdenadaRmBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRmCGa);
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 -2 2]);
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
% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaRpABa,'r',Abscisa,OrdenadaRpBCa,'g',Abscisa,OrdenadaRpCAa,'b',...
Abscisa,OrdenadaRpAGa,'m',Abscisa,OrdenadaRpBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRpCGa);
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 -2 2]);
title(cadenamc)
xlabel('m=p'),ylabel(cadylabel4a)
legend('Ref','AB','BC','CA','AG','BG','CG','z1')
end
%**************************************
case{3} % Salida 67N
extr=input('Extremo: Sm=1, Rm=2, Sp=3, Rp=4: ');
switch(extr)
case{1}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSm,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirSm(grados)')
title(dircadenamc)
rayar(1);
legend('AngDirSm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenamc)
xlabel('m=p'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIsmop,'b');
grid
xlabel('m=p'),ylabel('modulo Ismop(kA)')
case{2}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRm,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirRm(grados)')
title(dircadenamc)
rayar(1);
legend('AngDirRm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenamc)
xlabel('m=p'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIrmop,'b');
grid
xlabel('m=p'),ylabel('modulo Irmop(kA)')
case{3}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSp,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirSp(grados)')
title(dircadenamc)
rayar(1);
legend('AngDirSp','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenamc)
xlabel('m=p'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIspop,'b');
grid
xlabel('m=p'),ylabel('modulo Ispop(kA)')
case{4}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRp,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirRp(grados)')
title(dircadenamc)
rayar(1);
legend('AngDirRp','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenamc)
xlabel('m=p'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIrpop,'b');
grid
xlabel('m=p'),ylabel('modulo Irpop(kA)')
end % end del switch extr
case{4}, % ********* Salida 4 ************
extr=input('Extremo: Sm=1, Rm=2, Sp=3, Rp=4: ');
switch(extr)
case{1}
proteccion='Sm';
pzona=0.8;
hold on
plot19(L,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,subvectorSm_ZCG,'marron')
hold off
RXcadenamc=sprintf('%s,Proteccion:%s',cadenamc,proteccion);
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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,invsubvectorRm_ZCG,'marron')
hold off
RXcadenamc=sprintf('%s,Proteccion:%s',cadenamc,proteccion);
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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,subvectorSp_ZCG,'marron')
hold off
RXcadenamc=sprintf('%s,Proteccion:%s',cadenamc,proteccion);
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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,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,zl1u,invsubvectorRp_ZCG,'marron')
hold off
RXcadenamc=sprintf('%s,Proteccion:%s',cadenamc,proteccion);
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.
