% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG73c\n');
fprintf('Esquema Linea TRiple SR\n');
fprintf('Unominal kV: 400,220,132,66\n');
fprintf('Ctes. de la Linea Triple, por Linea y por km, segun Unominal\n');
fprintf('Eleccion: kE,L,Sccs,Sccr,Kr0,k0=zs0/zr0\n');
fprintf('Impedancia de Transferencia Si,No\n');
fprintf('Eleccion:(zsr1,Ksr0)\n');
%fprintf('Compensacion 21, Si, No\n');
fprintf('Falta Triple, m=p=q=0...1\n');
fprintf('Barra m=1,p=2,q=3\n');
fprintf('Analisis: Extremos Sm,Sp,Sq\n');
%fprintf('Salida Grafica 1: Extremos Sm, Rm: Unidad Preferente (pu)\n');
%fprintf(' : Extremos Sp, Rp: Unidad Preferente (pu)\n');
%fprintf(' : Extremos Sq, Rq: Unidad Preferente (pu)\n');
%fprintf('Salida Grafica 2: Extremos Sm, Rm: Unidad Preferente (pua)\n');
%fprintf(' : Extremos Sp, Rp: Unidad Preferente (pua)\n');
%fprintf(' : Extremos Sq, Rq: Unidad Preferente (pua)\n');
%fprintf('Salida Grafica 3: Extremo Sm: Todas las Unidades (pua)\n');
%fprintf(' : Extremo Sp: Todas las Unidades (pua)\n');
%fprintf(' : Extremo Sq: Todas las Unidades (pua)\n');
%fprintf(' : Extremo Rm: Todas las Unidades (pua)\n');
%fprintf(' : Extremo Rp: Todas las Unidades (pua)\n');
%fprintf(' : Extremo Rq: Todas las Unidades (pua)\n');
%fprintf('Salida Grafica 4: Analisis Sm, 67N, AngDir, Vpol, Iop\n');
%fprintf(' : Analisis Rm, 67N, AngDir, Vpol, Iop\n');
%fprintf(' : Analisis Sp, 67N, AngDir, Vpol, Iop\n');
%fprintf(' : Analisis Rp, 67N, AngDir, Vpol, Iop\n');
%fprintf(' : Analisis Sq, 67N, AngDir, Vpol, Iop\n');
%fprintf(' : Analisis Rq, 67N, AngDir, Vpol, Iop\n');
%fprintf('Salida Grafica 5: 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(' : Analisis Sq, Todas las Unidades R-X\n');
%fprintf(' : Analisis Rq, Todas las Unidades R-X\n');
fprintf('\n');
% -------------------------------------------------------------------
datos=input('datos: standard=1,por input=2: ');
if datos==1
Un=220;
kE=1;
L=220;
Sccs=5000;
Sccr=5000;
Kr0=1;
imptransfer=2;
compens21=2;
% ---------------------------------------------------------------------
else
Un=input('Un(Fase-Fase) kV: ');
kE=input('kE: ');
L=input('Longitud Linea km: ');
Sccs=input('Sccs MVA: ');
Sccr=input('Sccr MVA: ');
Kr0=input('Fuente R, 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=2;
end
% ---------------------------------------------------------------------
k0=input('k0=zs0/zr0: ');
% ---------------------------------------------------------------------
% CALCULOS
% ---------------------------------------------------------------------
% PRELIMINARES
%Constantes de la Linea Doble, por Linea y por km
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
zr1=(Un^2/Sccr)*i; zr2=zr1; zr0=Kr0*zr1;
zs1=(Un^2/Sccs)*i; zs2=zs1; zs0=k0*zr0;
zl1=zl1u*L; zl2=zl1; zl0=zl0u*L; zlM0=zlM0u*L;
E=Un/sqrt(3)+0i;
k=(zl0u-zl1u)/(3*zl1u);
%Km=zlM0u/zl1u;
% -------------------------------------------------------------------
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
faltaq=input('Falta en q: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :');
Rfq=input('Resistencia de Falta en q. Rf(Ohm) = ');
if faltaq==7|faltaq==8|faltaq==9|faltaq==11
Rgq=input('Resistencia de Falta en q. Rg(Ohm) = ');
else
Rgq=0;
end
% -------------------------------------------------------------------
%faltam=1;
%faltap=2
%faltaq=3;
% -------------------------------------------------------------------
i=1;
m=0.01;
while m<1;
p=m;
q=m;
switch(imptransfer)
case{1},
[Minc]=subMincT(1);
[z0]=subz0T(zs0,zr0,zsr0,zl0,zlM0,m,p,q);
[z1]=subz1T(zs1,zr1,zsr1,zl1,m,p,q);
[z2]=subz2T(zs2,zr2,zsr2,zl2,m,p,q);
case(2)
[Minc]=subMinc(1);
[z0]=subz0(zs0,zr0,zl0,zlM0,m,p,q);
[z1]=subz1(zs1,zr1,zl1,m,p,q);
[z2]=subz2(zs2,zr2,zl2,m,p,q);
end
% YBUS(0,1,2)
[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,Iq0,Iq1,Iq2,Vq0,Vq1,Vq2]=...
ma2mp(Z0,Z1,Z2,E,faltam,faltap,faltaq,Rfm,Rgm,Rfp,Rgp,Rfq,Rgq);
% IBUS(0,1,2)
I0=[-Im0;-Ip0;-Iq0;0;0];
I1=[-Im1;-Ip1;-Iq1;0;0];
I2=[-Im2;-Ip2;-Iq2;0;0];
% VBUS(0,1,2)
V0=Z0*I0;
V1=Z1*I1+E;
V2=Z2*I2;
% Ie(0,1,2) DE LOS ELEMENTOS
[Ie0]=subIe0(Minc,V0,z0);
[Ie1]=subIe1(Minc,V1,z1);
[Ie2]=subIe2(Minc,V2,z2);
% ASIGNACIONES V(0,1,2)
Vs0=V0(4); Vr0=V0(5);
Vs1=V1(4); Vr1=V1(5);
Vs2=V2(4); Vr2=V2(5);
% ASIGNACIONES I(0,1,2)
Ism0=Ie0(3); Irm0=Ie0(6); Isp0=Ie0(4); Irp0=Ie0(7); Isq0=Ie0(5); Irq0=Ie0(8);
Ism1=Ie1(3); Irm1=Ie1(6); Isp1=Ie1(4); Irp1=Ie1(7); Isq1=Ie1(5); Irq1=Ie1(8);
Ism2=Ie2(3); Irm2=Ie2(6); Isp2=Ie2(4); Irp2=Ie2(7); Isq2=Ie2(5); Irq2=Ie2(8);
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS PUNTOS m,p,q 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);
[Vqa,Vqb,Vqc]=abc(Vq0,Vq1,Vq2);
[Iqa,Iqb,Iqc]=abc(Iq0,Iq1,Iq2);
%
% 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);
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2);
[Ispa,Ispb,Ispc]=abc(Isp0,Isp1,Isp2);
[Irpa,Irpb,Irpc]=abc(Irp0,Irp1,Irp2);
[Isqa,Isqb,Isqc]=abc(Isq0,Isq1,Isq2);
[Irqa,Irqb,Irqc]=abc(Irq0,Irq1,Irq2);
%
% CALCULO DE LAS MEDIDAS DE LAS PROTECCIONES DE DISTANCIA EN
% LOS EXTREMOS DE LINEA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
PDist(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
PDist(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k);
[Sq_ZAB,Sq_ZBC,Sq_ZCA,Sq_ZAG,Sq_ZBG,Sq_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isqa,Isqb,Isqc,Isq0,k);
[Rq_ZAB,Rq_ZBC,Rq_ZCA,Rq_ZAG,Rq_ZBG,Rq_ZCG]=...
PDist(Vra,Vrb,Vrc,Irqa,Irqb,Irqc,Irq0,k);
%
% CALCULO DE LAS MAGNITUDES DE ALIMETACION DE LA PROTECCION
% DIRECCIONAL DE TIERRA EN LOS EXTREMOS DE LINEA
[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0);
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0);
[Vspol,Ispop,AngDirSp]=PDire(Vs0,Isp0);
[Vrpol,Irpop,AngDirRp]=PDire(Vr0,Irp0);
[Vspol,Isqop,AngDirSq]=PDire(Vs0,Isq0);
[Vrpol,Irqop,AngDirRq]=PDire(Vr0,Irq0);
% PREPARACION SALIDAS
% -------------------------------------------
% SALIDAS
% -------------------------------------------
% para Salida 1,2
if faltam==1
cfaltam='ag';
alcanceSm=abs(Sm_ZAG)/cos(angle(Sm_ZAG)-angle(zl1u));
alcanceRm=abs(Rm_ZAG)/cos(angle(Rm_ZAG)-angle(zl1u));
elementom='AG';
cadenalabelSm='alcanceSmAGpu';
cadenalabelRm='alcanceRmAGpu';
cadenalabelSma='alcanceSmAGpua';
cadenalabelRma='alcanceRmAGpua';
elseif faltam==2
cfaltam='bg';
alcanceSm=abs(Sm_ZBG)/cos(angle(Sm_ZBG)-angle(zl1u));
alcanceRm=abs(Rm_ZBG)/cos(angle(Rm_ZBG)-angle(zl1u));
elementom='BG';
cadenalabelSm='alcanceSmBGpu';
cadenalabelRm='alcanceRmBGpu';
cadenalabelSma='alcanceSmBGpua';
cadenalabelRma='alcanceRmBGpua';
elseif faltam==3
cfaltam='cg';
alcanceSm=abs(Sm_ZCG)/cos(angle(Sm_ZCG)-angle(zl1u));
alcanceRm=abs(Rm_ZCG)/cos(angle(Rm_ZCG)-angle(zl1u));
elementom='CG';
cadenalabelSm='alcanceSmCGpu';
cadenalabelRm='alcanceRmCGpu';
cadenalabelSma='alcanceSmCGpua';
cadenalabelRma='alcanceRmCGpua';
elseif faltam==4
cfaltam='ab';
alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
alcanceRm=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));
elementom='AB';
cadenalabelSm='alcanceSmABpu';
cadenalabelRm='alcanceRmABpu';
cadenalabelSma='alcanceSmABpua';
cadenalabelRma='alcanceRmABpua';
elseif faltam==5
cfaltam='bc';
alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
elementom='BC';
cadenalabelSm='alcanceSmBCpu';
cadenalabelRm='alcanceRmBCpu';
cadenalabelSma='alcanceSmBCpua';
cadenalabelRma='alcanceRmBCpua';
elseif faltam==6
cfaltam='ca';
alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
alcanceRm=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));
elementom='CA';
cadenalabelSm='alcanceSmCApu';
cadenalabelRm='alcanceRmCApu';
cadenalabelSma='alcanceSmCApua';
cadenalabelRma='alcanceRmCApua';
elseif faltam==7
cfaltam='abg';
alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
alcanceRm=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));
elementom='AB';
cadenalabelSm='alcanceSmABpu';
cadenalabelRm='alcanceRmABpu';
cadenalabelSma='alcanceSmABpua';
cadenalabelRma='alcanceRmABpua';
elseif faltam==8
cfaltam='bcg';
alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
elementom='BC';
cadenalabelSm='alcanceSmBCpu';
cadenalabelRm='alcanceRmBCpu';
cadenalabelSma='alcanceSmBCpua';
cadenalabelRma='alcanceRmBCpua';
elseif faltam==9
cfaltam='cag';
alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
alcanceRm=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));
elementom='CA';
cadenalabelSm='alcanceSmCApu';
cadenalabelRm='alcanceRmCApu';
cadenalabelSma='alcanceSmCApua';
cadenalabelRma='alcanceRmCApua';
elseif faltam==10
cfaltam='abc';
alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
elementom='BC';
cadenalabelSm='alcanceSmBCpu';
cadenalabelRm='alcanceRmBCpu';
cadenalabelSma='alcanceSmBCpua';
cadenalabelRma='alcanceRmBCpua';
elseif faltam==11
cfaltam='abcg';
alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
elementom='BC';
cadenalabelSm='alcanceSmBCpu';
cadenalabelRm='alcanceRmBCpu';
cadenalabelSma='alcanceSmBCpua';
cadenalabelRma='alcanceRmBCpua';
end
if faltap==1
cfaltap='ag';
alcanceSp=abs(Sp_ZAG)/cos(angle(Sp_ZAG)-angle(zl1u));
alcanceRp=abs(Rp_ZAG)/cos(angle(Rp_ZAG)-angle(zl1u));
elementop='AG';
cadenalabelSp='alcanceSpAGpu';
cadenalabelRp='alcanceRpAGpu';
cadenalabelSpa='alcanceSpAGpua';
cadenalabelRpa='alcanceRpAGpua';
elseif faltap==2
cfaltap='bg';
alcanceSp=abs(Sp_ZBG)/cos(angle(Sp_ZBG)-angle(zl1u));
alcanceRp=abs(Rp_ZBG)/cos(angle(Rp_ZBG)-angle(zl1u));
elementop='BG';
cadenalabelSp='alcanceSpBGpu';
cadenalabelRp='alcanceRpBGpu';
cadenalabelSpa='alcanceSpBGpua';
cadenalabelRpa='alcanceRpBGpua';
elseif faltap==3
cfaltap='cg';
alcanceSp=abs(Sp_ZCG)/cos(angle(Sp_ZCG)-angle(zl1u));
alcanceRp=abs(Rp_ZCG)/cos(angle(Rp_ZCG)-angle(zl1u));
elementop='CG';
cadenalabelSp='alcanceSpCGpu';
cadenalabelRp='alcanceRpCGpu';
cadenalabelSpa='alcanceSpCGpua';
cadenalabelRpa='alcanceRpCGpua';
elseif faltap==4
cfaltap='ab';
alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
elementop='AB';
cadenalabelSp='alcanceSpABpu';
cadenalabelRp='alcanceRpABpu';
cadenalabelSpa='alcanceSpABpua';
cadenalabelRpa='alcanceRpABpua';
elseif faltap==5
cfaltap='bc';
alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
elementop='BC';
cadenalabelSp='alcanceSpBCpu';
cadenalabelRp='alcanceRpBCpu';
cadenalabelSpa='alcanceSpBCpua';
cadenalabelRpa='alcanceRpBCpua';
elseif faltap==6
cfaltap='ca';
alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
elementop='CA';
cadenalabelSp='alcanceSpCApu';
cadenalabelRp='alcanceRpCApu';
cadenalabelSpa='alcanceSpCApua';
cadenalabelRpa='alcanceRpCApua';
elseif faltap==7
cfaltap='abg';
alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
elementop='AB';
cadenalabelSp='alcanceSpABpu';
cadenalabelRp='alcanceRpABpu';
cadenalabelSpa='alcanceSpABpua';
cadenalabelRpa='alcanceRpABpua';
elseif faltap==8
cfaltap='bcg';
alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
elementop='BC';
cadenalabelSp='alcanceSpBCpu';
cadenalabelRp='alcanceRpBCpu';
cadenalabelSpa='alcanceSpBCpua';
cadenalabelRpa='alcanceRpBCpua';
elseif faltap==9
cfaltap='cag';
alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
elementop='CA';
cadenalabelSp='alcanceSpCApu';
cadenalabelRp='alcanceRpCApu';
cadenalabelSpa='alcanceSpCApua';
cadenalabelRpa='alcanceRpCApua';
elseif faltap==10
cfaltap='abc';
alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
elementop='BC';
cadenalabelSp='alcanceSpBCpu';
cadenalabelRp='alcanceRpBCpu';
cadenalabelSpa='alcanceSpBCpua';
cadenalabelRpa='alcanceRpBCpua';
elseif faltap==11
cfaltap='abcg';
alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
elementop='BC';
cadenalabelSp='alcanceSpBCpu';
cadenalabelRp='alcanceRpBCpu';
cadenalabelSpa='alcanceSpBCpua';
cadenalabelRpa='alcanceRpBCpua';
end
if faltaq==1
cfaltaq='ag';
alcanceSq=abs(Sq_ZAG)/cos(angle(Sq_ZAG)-angle(zl1u));
alcanceRq=abs(Rq_ZAG)/cos(angle(Rq_ZAG)-angle(zl1u));
elementoq='AG';
cadenalabelSq='alcanceSqAGpu';
cadenalabelRq='alcanceRqAGpu';
cadenalabelSqa='alcanceSqAGpua';
cadenalabelRqa='alcanceRqAGpua';
elseif faltaq==2
cfaltaq='bg';
alcanceSq=abs(Sq_ZBG)/cos(angle(Sq_ZBG)-angle(zl1u));
alcanceRq=abs(Rq_ZBG)/cos(angle(Rq_ZBG)-angle(zl1u));
elementoq='BG';
cadenalabelSq='alcanceSqBGpu';
cadenalabelRq='alcanceRqBGpu';
cadenalabelSqa='alcanceSqBGpua';
cadenalabelRqa='alcanceRqBGpua';
elseif faltaq==3
cfaltaq='cg';
alcanceSq=abs(Sq_ZCG)/cos(angle(Sq_ZCG)-angle(zl1u));
alcanceRq=abs(Rq_ZCG)/cos(angle(Rq_ZCG)-angle(zl1u));
elementoq='CG';
cadenalabelSq='alcanceSqCGpu';
cadenalabelRq='alcanceRqCGpu';
cadenalabelSqa='alcanceSqCGpua';
cadenalabelRqa='alcanceRqCGpua';
elseif faltaq==4
cfaltaq='ab';
alcanceSq=abs(Sq_ZAB)/cos(angle(Sq_ZAB)-angle(zl1u));
alcanceRq=abs(Rq_ZAB)/cos(angle(Rq_ZAB)-angle(zl1u));
elementoq='AB';
cadenalabelSq='alcanceSqABpu';
cadenalabelRq='alcanceRqABpu';
cadenalabelSqa='alcanceSqABpua';
cadenalabelRqa='alcanceRqABpua';
elseif faltaq==5
cfaltaq='bc';
alcanceSq=abs(Sq_ZBC)/cos(angle(Sq_ZBC)-angle(zl1u));
alcanceRq=abs(Rq_ZBC)/cos(angle(Rq_ZBC)-angle(zl1u));
elementoq='BC';
cadenalabelSq='alcanceSqBCpu';
cadenalabelRq='alcanceRqBCpu';
cadenalabelSqa='alcanceSqBCpua';
cadenalabelRqa='alcanceRqBCpua';
elseif faltaq==6
cfaltap='ca';
alcanceSq=abs(Sq_ZCA)/cos(angle(Sq_ZCA)-angle(zl1u));
alcanceRq=abs(Rq_ZCA)/cos(angle(Rq_ZCA)-angle(zl1u));
elementoq='CA';
cadenalabelSq='alcanceSqCApu';
cadenalabelRq='alcanceRqCApu';
cadenalabelSqa='alcanceSqCApua';
cadenalabelRqa='alcanceRqCApua';
elseif faltaq==7
cfaltaq='abg';
alcanceSq=abs(Sq_ZAB)/cos(angle(Sq_ZAB)-angle(zl1u));
alcanceRq=abs(Rq_ZAB)/cos(angle(Rq_ZAB)-angle(zl1u));
elementoq='AB';
cadenalabelSq='alcanceSqABpu';
cadenalabelRq='alcanceRqABpu';
cadenalabelSqa='alcanceSqABpua';
cadenalabelRqa='alcanceRqABpua';
elseif faltaq==8
cfaltaq='bcg';
alcanceSq=abs(Sq_ZBC)/cos(angle(Sq_ZBC)-angle(zl1u));
alcanceRq=abs(Rq_ZBC)/cos(angle(Rq_ZBC)-angle(zl1u));
elementoq='BC';
cadenalabelSq='alcanceSqBCpu';
cadenalabelRq='alcanceRqBCpu';
cadenalabelSqa='alcanceSqBCpua';
cadenalabelRqa='alcanceRqBCpua';
elseif faltaq==9
cfaltaq='cag';
alcanceSq=abs(Sq_ZCA)/cos(angle(Sq_ZCA)-angle(zl1u));
alcanceRq=abs(Rq_ZCA)/cos(angle(Rq_ZCA)-angle(zl1u));
elementoq='CA';
cadenalabelSq='alcanceSqCApu';
cadenalabelRq='alcanceRqCApu';
cadenalabelSqa='alcanceSqCApua';
cadenalabelRqa='alcanceRqCApua';
elseif faltaq==10
cfaltaq='abc';
alcanceSq=abs(Sq_ZBC)/cos(angle(Sq_ZBC)-angle(zl1u));
alcanceRq=abs(Rq_ZBC)/cos(angle(Rq_ZBC)-angle(zl1u));
elementoq='BC';
cadenalabelSq='alcanceSqBCpu';
cadenalabelRq='alcanceRqBCpu';
cadenalabelSqa='alcanceSqBCpua';
cadenalabelRqa='alcanceRqBCpua';
elseif faltaq==11
cfaltaq='abcg';
alcanceSq=abs(Sq_ZBC)/cos(angle(Sq_ZBC)-angle(zl1u));
alcanceRq=abs(Rq_ZBC)/cos(angle(Rq_ZBC)-angle(zl1u));
elementoq='BC';
cadenalabelSq='alcanceSqBCpu';
cadenalabelRq='alcanceRqBCpu';
cadenalabelSqa='alcanceSqBCpua';
cadenalabelRqa='alcanceRqBCpua';
end
alcanceSmpu=alcanceSm/abs(zl1u*m*L);
alcanceRmpu=alcanceRm/abs(zl1u*(1-m)*L);
alcanceSppu=alcanceSp/abs(zl1u*p*L);
alcanceRppu=alcanceRp/abs(zl1u*(1-p)*L);
alcanceSqpu=alcanceSq/abs(zl1u*q*L);
alcanceRqpu=alcanceRq/abs(zl1u*(1-q)*L);
alcanceSmpua=alcanceSm/abs(zl1u*L);
alcanceRmpua=alcanceRm/abs(zl1u*L);
alcanceSppua=alcanceSp/abs(zl1u*L);
alcanceRppua=alcanceRp/abs(zl1u*L);
alcanceSqpua=alcanceSq/abs(zl1u*L);
alcanceRqpua=alcanceRq/abs(zl1u*L);
Abscisa(1,i)=m;
OrdenadaSm(1,i)=alcanceSmpu;
OrdenadaRm(1,i)=alcanceRmpu;
OrdenadaSp(1,i)=alcanceSppu;
OrdenadaRp(1,i)=alcanceRppu;
OrdenadaSq(1,i)=alcanceSqpu;
OrdenadaRq(1,i)=alcanceRqpu;
OrdenadaSma(1,i)=alcanceSmpua;
OrdenadaRma(1,i)=alcanceRmpua;
OrdenadaSpa(1,i)=alcanceSppua;
OrdenadaRpa(1,i)=alcanceRppua;
OrdenadaSqa(1,i)=alcanceSqpua;
OrdenadaRqa(1,i)=alcanceRqpua;
% Para salida 3
alcanceSm_ZAB=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
alcanceSm_ZBC=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcanceSm_ZCA=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
alcanceSm_ZAG=abs(Sm_ZAG)/cos(angle(Sm_ZAG)-angle(zl1u));
alcanceSm_ZBG=abs(Sm_ZBG)/cos(angle(Sm_ZBG)-angle(zl1u));
alcanceSm_ZCG=abs(Sm_ZCG)/cos(angle(Sm_ZCG)-angle(zl1u));
alcanceSp_ZAB=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
alcanceSp_ZBC=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcanceSp_ZCA=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
alcanceSp_ZAG=abs(Sp_ZAG)/cos(angle(Sp_ZAG)-angle(zl1u));
alcanceSp_ZBG=abs(Sp_ZBG)/cos(angle(Sp_ZBG)-angle(zl1u));
alcanceSp_ZCG=abs(Sp_ZCG)/cos(angle(Sp_ZCG)-angle(zl1u));
alcanceSq_ZAB=abs(Sq_ZAB)/cos(angle(Sq_ZAB)-angle(zl1u));
alcanceSq_ZBC=abs(Sq_ZBC)/cos(angle(Sq_ZBC)-angle(zl1u));
alcanceSq_ZCA=abs(Sq_ZCA)/cos(angle(Sq_ZCA)-angle(zl1u));
alcanceSq_ZAG=abs(Sq_ZAG)/cos(angle(Sq_ZAG)-angle(zl1u));
alcanceSq_ZBG=abs(Sq_ZBG)/cos(angle(Sq_ZBG)-angle(zl1u));
alcanceSq_ZCG=abs(Sq_ZCG)/cos(angle(Sq_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);
alcanceSp_ZABpua=alcanceSp_ZAB/abs(zl1u*L);
alcanceSp_ZBCpua=alcanceSp_ZBC/abs(zl1u*L);
alcanceSp_ZCApua=alcanceSp_ZCA/abs(zl1u*L);
alcanceSp_ZAGpua=alcanceSp_ZAG/abs(zl1u*L);
alcanceSp_ZBGpua=alcanceSp_ZBG/abs(zl1u*L);
alcanceSp_ZCGpua=alcanceSp_ZCG/abs(zl1u*L);
alcanceSq_ZABpua=alcanceSq_ZAB/abs(zl1u*L);
alcanceSq_ZBCpua=alcanceSq_ZBC/abs(zl1u*L);
alcanceSq_ZCApua=alcanceSq_ZCA/abs(zl1u*L);
alcanceSq_ZAGpua=alcanceSq_ZAG/abs(zl1u*L);
alcanceSq_ZBGpua=alcanceSq_ZBG/abs(zl1u*L);
alcanceSq_ZCGpua=alcanceSq_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;
OrdenadaSp_ZABpua(1,i)=alcanceSp_ZABpua;
OrdenadaSp_ZBCpua(1,i)=alcanceSp_ZBCpua;
OrdenadaSp_ZCApua(1,i)=alcanceSp_ZCApua;
OrdenadaSp_ZAGpua(1,i)=alcanceSp_ZAGpua;
OrdenadaSp_ZBGpua(1,i)=alcanceSp_ZBGpua;
OrdenadaSp_ZCGpua(1,i)=alcanceSp_ZCGpua;
OrdenadaSq_ZABpua(1,i)=alcanceSq_ZABpua;
OrdenadaSq_ZBCpua(1,i)=alcanceSq_ZBCpua;
OrdenadaSq_ZCApua(1,i)=alcanceSq_ZCApua;
OrdenadaSq_ZAGpua(1,i)=alcanceSq_ZAGpua;
OrdenadaSq_ZBGpua(1,i)=alcanceSq_ZBGpua;
OrdenadaSq_ZCGpua(1,i)=alcanceSq_ZCGpua;
%***************************************
alcanceRm_ZAB=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));
alcanceRm_ZBC=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
alcanceRm_ZCA=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));
alcanceRm_ZAG=abs(Rm_ZAG)/cos(angle(Rm_ZAG)-angle(zl1u));
alcanceRm_ZBG=abs(Rm_ZBG)/cos(angle(Rm_ZBG)-angle(zl1u));
alcanceRm_ZCG=abs(Rm_ZCG)/cos(angle(Rm_ZCG)-angle(zl1u));
alcanceRp_ZAB=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
alcanceRp_ZBC=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
alcanceRp_ZCA=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
alcanceRp_ZAG=abs(Rp_ZAG)/cos(angle(Rp_ZAG)-angle(zl1u));
alcanceRp_ZBG=abs(Rp_ZBG)/cos(angle(Rp_ZBG)-angle(zl1u));
alcanceRp_ZCG=abs(Rp_ZCG)/cos(angle(Rp_ZCG)-angle(zl1u));
alcanceRq_ZAB=abs(Rq_ZAB)/cos(angle(Rq_ZAB)-angle(zl1u));
alcanceRq_ZBC=abs(Rq_ZBC)/cos(angle(Rq_ZBC)-angle(zl1u));
alcanceRq_ZCA=abs(Rq_ZCA)/cos(angle(Rq_ZCA)-angle(zl1u));
alcanceRq_ZAG=abs(Rq_ZAG)/cos(angle(Rq_ZAG)-angle(zl1u));
alcanceRq_ZBG=abs(Rq_ZBG)/cos(angle(Rq_ZBG)-angle(zl1u));
alcanceRq_ZCG=abs(Rq_ZCG)/cos(angle(Rq_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);
alcanceRp_ZABpua=alcanceRp_ZAB/abs(zl1u*L);
alcanceRp_ZBCpua=alcanceRp_ZBC/abs(zl1u*L);
alcanceRp_ZCApua=alcanceRp_ZCA/abs(zl1u*L);
alcanceRp_ZAGpua=alcanceRp_ZAG/abs(zl1u*L);
alcanceRp_ZBGpua=alcanceRp_ZBG/abs(zl1u*L);
alcanceRp_ZCGpua=alcanceRp_ZCG/abs(zl1u*L);
alcanceRq_ZABpua=alcanceRq_ZAB/abs(zl1u*L);
alcanceRq_ZBCpua=alcanceRq_ZBC/abs(zl1u*L);
alcanceRq_ZCApua=alcanceRq_ZCA/abs(zl1u*L);
alcanceRq_ZAGpua=alcanceRq_ZAG/abs(zl1u*L);
alcanceRq_ZBGpua=alcanceRq_ZBG/abs(zl1u*L);
alcanceRq_ZCGpua=alcanceRq_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;
OrdenadaRp_ZABpua(1,i)=alcanceRp_ZABpua;
OrdenadaRp_ZBCpua(1,i)=alcanceRp_ZBCpua;
OrdenadaRp_ZCApua(1,i)=alcanceRp_ZCApua;
OrdenadaRp_ZAGpua(1,i)=alcanceRp_ZAGpua;
OrdenadaRp_ZBGpua(1,i)=alcanceRp_ZBGpua;
OrdenadaRp_ZCGpua(1,i)=alcanceRp_ZCGpua;
OrdenadaRq_ZABpua(1,i)=alcanceRq_ZABpua;
OrdenadaRq_ZBCpua(1,i)=alcanceRq_ZBCpua;
OrdenadaRq_ZCApua(1,i)=alcanceRq_ZCApua;
OrdenadaRq_ZAGpua(1,i)=alcanceRq_ZAGpua;
OrdenadaRq_ZBGpua(1,i)=alcanceRq_ZBGpua;
OrdenadaRq_ZCGpua(1,i)=alcanceRq_ZCGpua;
%***************************************
%Para salida 4, 67N
%***************************************
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;
vectorVspol(1,i)=Vspol;
vectorIsqop(1,i)=Isqop;
vectorAngDirSq(1,i)=AngDirSq;
vectorVrpol(1,i)=Vrpol;
vectorIrqop(1,i)=Irqop;
vectorAngDirRq(1,i)=AngDirRq;
%***************************************
% Para Salida 4 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;
vectorSq_ZAB(1,i)=Sq_ZAB;
vectorSq_ZBC(1,i)=Sq_ZBC;
vectorSq_ZCA(1,i)=Sq_ZCA;
vectorSq_ZAG(1,i)=Sq_ZAG;
vectorSq_ZBG(1,i)=Sq_ZBG;
vectorSq_ZCG(1,i)=Sq_ZCG;
vectorRq_ZAB(1,i)=Rq_ZAB;
vectorRq_ZBC(1,i)=Rq_ZBC;
vectorRq_ZCA(1,i)=Rq_ZCA;
vectorRq_ZAG(1,i)=Rq_ZAG;
vectorRq_ZBG(1,i)=Rq_ZBG;
vectorRq_ZCG(1,i)=Rq_ZCG;
%***************************************
i=i+1;
m=m+0.01;
end
% ------------------------------------------------------
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
% ------------------------------------------------------
cadenam1=sprintf('PG173ic,%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);
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: %s,%s,Rfp=%-3.0f',cfaltap,elementop,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
cadenap1=sprintf('Rgp=%-3.0f',Rgp);
cadenap=sprintf('%s,%s',cadenap,cadenap1);
end
cadenamp=sprintf('%s, %s',cadenam,cadenap);
cadenaq=sprintf('q: %s,%s,Rfq=%-3.0f',cfaltaq,elementoq,Rfq);
if faltaq==7|faltaq==8|faltaq==9|faltaq==11
cadenaq1=sprintf('Rgq=%-3.0f',Rgq);
cadenaq=sprintf('%s,%s',cadenaq,cadenaq1);
end
cadenampq=sprintf('%s, %s, %s',cadenam,cadenap,cadenaq);
% ------------------------------------------------------
cadenamc1=sprintf('PG173ic,%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: %s,Rfp=%-3.0f',cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
cadenapc1=sprintf('Rgp=%-3.0f',Rgp);
cadenapc=sprintf('%s,%s',cadenapc,cadenapc1);
end
cadenampc=sprintf('%s, %s',cadenamc,cadenapc);
cadenaqc=sprintf('q: %s,Rfq=%-3.0f',cfaltaq,Rfq);
if faltaq==7|faltaq==8|faltaq==9|faltaq==11
cadenaqc1=sprintf('Rgq=%-3.0f',Rgq);
cadenaqc=sprintf('%s,%s',cadenaqc,cadenaqc1);
end
cadenampqc=sprintf('%s, %s',cadenampc,cadenaqc);
%***************************************
% Titulo para las gráficas 67N
%***************************************
dircadenamc1=sprintf('PG173ic,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: %s,Rfp=%-3.0f',cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
dircadenapc1=sprintf('Rgp=%-3.0f',Rgp);
dircadenapc=sprintf('%s,%s',dircadenapc,dircadenapc1);
end
dircadenampc=sprintf('%s, %s',dircadenamc,dircadenapc);
dircadenaqc=sprintf('q: %s,Rfq=%-3.0f',cfaltaq,Rfq);
if faltaq==7|faltaq==8|faltaq==9|faltaq==11
dircadenaqc1=sprintf('Rgq=%-3.0f',Rgq);
dircadenaqc=sprintf('%s,%s',dircadenaqc,dircadenaqc1);
end
dircadenampqc=sprintf('%s, %s',dircadenampc,dircadenaqc);
%*************************************
variassalidas=1;
while variassalidas==1
fprintf('Salida Grafica 1: Extremos Sm, Rm: Unidad Preferente (pu)\n');
fprintf(' : Extremos Sp, Rp: Unidad Preferente (pu)\n');
fprintf(' : Extremos Sq, Rq: Unidad Preferente (pu)\n');
fprintf('Salida Grafica 2: Extremos Sm, Rm: Unidad Preferente (pua)\n');
fprintf(' : Extremos Sp, Rp: Unidad Preferente (pua)\n');
fprintf(' : Extremos Sq, Rq: Unidad Preferente (pua)\n');
fprintf('Salida Grafica 3: Extremo Sm: Todas las Unidades (pua)\n');
fprintf(' : Extremo Sp: Todas las Unidades (pua)\n');
fprintf(' : Extremo Sq: Todas las Unidades (pua)\n');
fprintf(' : Extremo Rm: Todas las Unidades (pua)\n');
fprintf(' : Extremo Rp: Todas las Unidades (pua)\n');
fprintf(' : Extremo Rq: Todas las Unidades (pua)\n');
fprintf('Salida Grafica 4: Analisis Sm, 67N\n');
fprintf(' : Analisis Rm, 67N\n');
fprintf(' : Analisis Sp, 67N\n');
fprintf(' : Analisis Rp, 67N\n');
fprintf(' : Analisis Sq, 67N\n');
fprintf(' : Analisis Rq, 67N\n');
fprintf('Salida Grafica 5: 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(' : Analisis Sq, Todas las Unidades R-X\n');
fprintf(' : Analisis Rq, Todas las Unidades R-X\n');
fprintf('\n');
salida=input('Salida: 1,2,3,4,5 :');
switch(salida)
case{1}, % ******** Unidad Preferente (pu) ********
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=1.0;
Ordenadadi(cont)=0.8/mm;
Ordenadadd(cont)=0.8/(1-mm);
cont=cont+1;
mm=mm+0.01;
end
extremos=input('Extremos: m=1, p=2, q=3: ');
switch(extremos)
case{1},
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'m',Abscisa,OrdenadaSm,'r'),grid on
axis([0 1 0 3]);
title(cadenampq)
xlabel('m=p=q'),ylabel(cadenalabelSm)
legend('Ref','Lim',cadenalabelSm)
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadd,'m',Abscisa,OrdenadaRm,'r'),grid on
axis([0 1 0 3]);
%title(cadenamp)
xlabel('m=p=q'),ylabel(cadenalabelRm)
legend('Ref','Lim',cadenalabelRm)
case{2},
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'m',Abscisa,OrdenadaSp,'r'),grid on
axis([0 1 0 3]);
title(cadenampq)
xlabel('m=p=q'),ylabel(cadenalabelSp)
legend('Ref','Lim',cadenalabelSp)
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadd,'m',Abscisa,OrdenadaRp,'r'),grid on
axis([0 1 0 3]);
%title(cadenamp)
xlabel('m=p=q'),ylabel(cadenalabelRp)
legend('Ref','Lim',cadenalabelRp)
case{3},
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'m',Abscisa,OrdenadaSq,'r'),grid on
axis([0 1 0 3]);
title(cadenampq)
xlabel('m=p=q'),ylabel(cadenalabelSq)
legend('Ref','Lim',cadenalabelSq)
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadd,'m',Abscisa,OrdenadaRq,'r'),grid on
axis([0 1 0 3]);
%title(cadenamp)
xlabel('m=p=q'),ylabel(cadenalabelRq)
legend('Ref','Lim',cadenalabelRq)
end % extremos
case(2) % ******** Unidad Preferente (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
extremos=input('Extremos: m=1, p=2, q=3: ');
switch(extremos)
case{1},
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadaz1,'g',Abscisa,OrdenadaSma,'r'),grid on
title(cadenampq)
xlabel('m=p=q'),ylabel(cadenalabelSma)
legend('Ref','z1',cadenalabelSma)
subplot(2,1,2)
plot(Abscisa00,Ordenada00,'c',Abscisa0,Ordenadaz1,'g',Abscisa,OrdenadaRma,'r'),grid on
%title(cadenamp)
xlabel('m=p=q'),ylabel(cadenalabelRma)
legend('Ref','z1',cadenalabelRma)
case{2},
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadaz1,'g',Abscisa,OrdenadaSpa,'r'),grid on
title(cadenampq)
xlabel('m=p=q'),ylabel(cadenalabelSpa)
legend('Ref','z1',cadenalabelSpa)
subplot(2,1,2)
plot(Abscisa00,Ordenada00,'c',Abscisa0,Ordenadaz1,'g',Abscisa,OrdenadaRpa,'r'),grid on
%title(cadenamp)
xlabel('m=p=q'),ylabel(cadenalabelRpa)
legend('Ref','z1',cadenalabelRpa)
case{3},
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadaz1,'g',Abscisa,OrdenadaSqa,'r'),grid on
title(cadenampq)
xlabel('m=p=q'),ylabel(cadenalabelSqa)
legend('Ref','z1',cadenalabelSqa)
subplot(2,1,2)
plot(Abscisa00,Ordenada00,'c',Abscisa0,Ordenadaz1,'g',Abscisa,OrdenadaRqa,'r'),grid on
%title(cadenamp)
xlabel('m=p=q'),ylabel(cadenalabelRqa)
legend('Ref','z1',cadenalabelRqa)
end
case(3) % ******** Todas las Unidades (pua) ********
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenada0a(cont)=mm;
Ordenada0b(cont)=1-mm;
cont=cont+1;
mm=mm+0.01;
end
extremos=input('Extremo Sm=1, Sp=2, Sq=3, Rm=4, Rp=5, Rq=6 :');
switch(extremos)
case{1}, % salida 3, Extremo Sm
cadylabel1a='alcSmpua';
% Colores
b=plot(Abscisa0,Ordenada0a);
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,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 -2 2]);
title(cadenampqc)
xlabel('m=p=q'),ylabel(cadylabel1a)
legend('45º','AB','BC','CA','AG','BG','CG','z1')
case{2}, % salida 3, Extremo Sp
cadylabel1a='alcSppua';
% Colores
b=plot(Abscisa0,Ordenada0a);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaSp_ZABpua,'r',Abscisa,OrdenadaSp_ZBCpua,'g',Abscisa,OrdenadaSp_ZCApua,'b',...
Abscisa,OrdenadaSp_ZAGpua,'m',Abscisa,OrdenadaSp_ZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSp_ZCGpua);
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(cadenampqc)
xlabel('m=p=q'),ylabel(cadylabel1a)
legend('45º','AB','BC','CA','AG','BG','CG','z1')
case{3}, % salida 3, Extremo Sq
cadylabel1a='alcSqpua';
% Colores
b=plot(Abscisa0,Ordenada0a);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaSq_ZABpua,'r',Abscisa,OrdenadaSq_ZBCpua,'g',Abscisa,OrdenadaSq_ZCApua,'b',...
Abscisa,OrdenadaSq_ZAGpua,'m',Abscisa,OrdenadaSq_ZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSq_ZCGpua);
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(cadenampqc)
xlabel('m=p=q'),ylabel(cadylabel1a)
legend('45º','AB','BC','CA','AG','BG','CG','z1')
case{4}, % salida 3, Extremo Rm
cadylabel1a='alcRmpua';
% Colores
b=plot(Abscisa0,Ordenada0b);
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,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 -2 2]);
title(cadenampqc)
xlabel('m=p=q'),ylabel(cadylabel1a)
legend('Ref','AB','BC','CA','AG','BG','CG','z1')
case{5}, % salida 3, Extremo Rp
cadylabel1a='alcRppua';
% Colores
b=plot(Abscisa0,Ordenada0b);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaRp_ZABpua,'r',Abscisa,OrdenadaRp_ZBCpua,'g',Abscisa,OrdenadaRp_ZCApua,'b',...
Abscisa,OrdenadaRp_ZAGpua,'m',Abscisa,OrdenadaRp_ZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRp_ZCGpua);
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(cadenampqc)
xlabel('m=p=q'),ylabel(cadylabel1a)
legend('Ref','AB','BC','CA','AG','BG','CG','z1')
case{6}, % salida 3, Extremo Rq
cadylabel1a='alcRqpua';
% Colores
b=plot(Abscisa0,Ordenada0b);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaRq_ZABpua,'r',Abscisa,OrdenadaRq_ZBCpua,'g',Abscisa,OrdenadaRq_ZCApua,'b',...
Abscisa,OrdenadaRq_ZAGpua,'m',Abscisa,OrdenadaRq_ZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaRq_ZCGpua);
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(cadenampqc)
xlabel('m=p=q'),ylabel(cadylabel1a)
legend('Ref','AB','BC','CA','AG','BG','CG','z1')
end
%******************************************
case{4} % ******** Salida 67N ********
extr=input('Extremo: Sm=1,Rm=2,Sp=3,Rp=4,Sq=5,Rq=6: ');
switch(extr)
case{1}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSm,'b');
axis([0 1 -180 180]);
xlabel('m=p=q'),ylabel('AngDirSm(grados)')
title(dircadenampqc)
rayar(1);
legend('AngDirSm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenampqc)
xlabel('m=p=q'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIsmop,'b');
grid
xlabel('m=p=q'),ylabel('modulo Ismop(kA)')
case{2}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRm,'b');
axis([0 1 -180 180]);
xlabel('m=p=q'),ylabel('AngDirRm(grados)')
title(dircadenampqc)
rayar(1);
legend('AngDirRm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenampqc)
xlabel('m=p=q'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIrmop,'b');
grid
xlabel('m=p=q'),ylabel('modulo Irmop(kA)')
case{3}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSp,'b');
axis([0 1 -180 180]);
xlabel('m=p=q'),ylabel('AngDirSp(grados)')
title(dircadenampqc)
rayar(1);
legend('AngDirSp','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenampqc)
xlabel('m=p=q'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIspop,'b');
grid
xlabel('m=p=q'),ylabel('modulo Ispop(kA)')
case{4}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRp,'b');
axis([0 1 -180 180]);
xlabel('m=p=q'),ylabel('AngDirRp(grados)')
title(dircadenampqc)
rayar(1);
legend('AngDirRp','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenampqc)
xlabel('m=p=q'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIrpop,'b');
grid
xlabel('m=p=q'),ylabel('modulo Irpop(kA)')
case{5}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSq,'b');
axis([0 1 -180 180]);
xlabel('m=p=q'),ylabel('AngDirSq(grados)')
title(dircadenampqc)
rayar(1);
legend('AngDirSq','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenampqc)
xlabel('m=p=q'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIsqop,'b');
grid
xlabel('m=p=q'),ylabel('modulo Isqop(kA)')
case{6}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRq,'b');
axis([0 1 -180 180]);
xlabel('m=p=q'),ylabel('AngDirRq(grados)')
title(dircadenampqc)
rayar(1);
legend('AngDirRq','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenampqc)
xlabel('m=p=q'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIrqop,'b');
grid
xlabel('m=p=q'),ylabel('modulo Irqop(kA)')
end % end del switch extr
%****************************************
case{5}, % ************ Salida R-X ***************
extr=input('Extremo: Sm=1, Rm=2, Sp=3, Rp=4, Sq=5, Rq=6: ');
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
RXcadenam=sprintf('Proteccion: %s',proteccion);
RXcadenampqc=sprintf('%s,%s',cadenampqc,RXcadenam);
title(RXcadenampqc)
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
RXcadenam=sprintf('Proteccion: %s',proteccion);
RXcadenampqc=sprintf('%s,%s',cadenampqc,RXcadenam);
title(RXcadenampqc)
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
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenampqc=sprintf('%s,%s',cadenampqc,RXcadenap);
title(RXcadenampqc)
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
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenampqc=sprintf('%s,%s',cadenampqc,RXcadenap);
title(RXcadenampqc)
legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG')
case{5}, % Proteccion Sq
proteccion='Sq';
pzona=0.8;
hold on
plot19(L,zl1u,pzona,proteccion)
hold off
hold on
% Colores
a=plot(vectorSq_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorSq_ZBC,'m')
plot(vectorSq_ZCA,'c')
plot(vectorSq_ZAG,'r')
plot(vectorSq_ZBG,'k')
% Colores
c=plot(vectorSq_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
[subAbscisa,subvectorSq_ZAB]=subconjunto21(Abscisa,vectorSq_ZAB);
subAbscisa(:,1)=[];
subvectorSq_ZAB(:,1)=[];
hold on
b=plot(subvectorSq_ZAB,'*');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,subvectorSq_ZAB,'naranja')
hold off
[subAbscisa,subvectorSq_ZBC]=subconjunto21(Abscisa,vectorSq_ZBC);
subAbscisa(:,1)=[];
subvectorSq_ZBC(:,1)=[];
hold on
plot(subvectorSq_ZBC,'m*');
plot20(L,zl1u,subvectorSq_ZBC,'m')
hold off
[subAbscisa,subvectorSq_ZCA]=subconjunto21(Abscisa,vectorSq_ZCA);
subAbscisa(:,1)=[];
subvectorSq_ZCA(:,1)=[];
hold on
plot(subvectorSq_ZCA,'c*');
plot20(L,zl1u,subvectorSq_ZCA,'c')
hold off
[subAbscisa,subvectorSq_ZAG]=subconjunto21(Abscisa,vectorSq_ZAG);
subAbscisa(:,1)=[];
subvectorSq_ZAG(:,1)=[];
hold on
plot(subvectorSq_ZAG,'r*');
plot20(L,zl1u,subvectorSq_ZAG,'r')
hold off
[subAbscisa,subvectorSq_ZBG]=subconjunto21(Abscisa,vectorSq_ZBG);
subAbscisa(:,1)=[];
subvectorSq_ZBG(:,1)=[];
hold on
plot(subvectorSq_ZBG,'k*');
plot20(L,zl1u,subvectorSq_ZBG,'k')
hold off
[subAbscisa,subvectorSq_ZCG]=subconjunto21(Abscisa,vectorSq_ZCG);
subAbscisa(:,1)=[];
subvectorSq_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorSq_ZCG,'*');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,subvectorSq_ZCG,'marron')
hold off
RXcadenaq=sprintf('Proteccion: %s',proteccion);
RXcadenampqc=sprintf('%s,%s',cadenampqc,RXcadenaq);
title(RXcadenampqc)
legend('Mho','Linea','m=0..1','AB','BC','CA','AG','BG','CG')
case{6} % Proteccion Rq
proteccion='Rq';
pzona=0.8;
hold on
plot19(L,zl1u,pzona,proteccion)
hold off
hold on
% Colores
a=plot(vectorRq_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorRq_ZBC,'m')
plot(vectorRq_ZCA,'c')
plot(vectorRq_ZAG,'r')
plot(vectorRq_ZBG,'k')
% Colores
c=plot(vectorRq_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
[subAbscisa,subvectorRq_ZAB]=subconjunto21(Abscisa,vectorRq_ZAB);
subAbscisa(:,11)=[];
subvectorRq_ZAB(:,11)=[];
invsubvectorRq_ZAB=subvectorRq_ZAB(1:1,10:-1:1);
hold on
b=plot(invsubvectorRq_ZAB,'*');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,invsubvectorRq_ZAB,'naranja')
hold off
[subAbscisa,subvectorRq_ZBC]=subconjunto21(Abscisa,vectorRq_ZBC);
subAbscisa(:,11)=[];
subvectorRq_ZBC(:,11)=[];
invsubvectorRq_ZBC=subvectorRq_ZBC(1:1,10:-1:1);
hold on
plot(subvectorRq_ZBC,'m*');
plot20(L,zl1u,invsubvectorRq_ZBC,'m')
hold off
[subAbscisa,subvectorRq_ZCA]=subconjunto21(Abscisa,vectorRq_ZCA);
subAbscisa(:,11)=[];
subvectorRq_ZCA(:,11)=[];
invsubvectorRq_ZCA=subvectorRq_ZCA(1:1,10:-1:1);
hold on
plot(subvectorRq_ZCA,'c*');
plot20(L,zl1u,invsubvectorRq_ZCA,'c')
hold off
[subAbscisa,subvectorRq_ZAG]=subconjunto21(Abscisa,vectorRq_ZAG);
subAbscisa(:,11)=[];
subvectorRq_ZAG(:,11)=[];
invsubvectorRq_ZAG=subvectorRq_ZAG(1:1,10:-1:1);
hold on
plot(subvectorRq_ZAG,'r*');
plot20(L,zl1u,invsubvectorRq_ZAG,'r')
hold off
[subAbscisa,subvectorRq_ZBG]=subconjunto21(Abscisa,vectorRq_ZBG);
subAbscisa(:,11)=[];
subvectorRq_ZBG(:,11)=[];
invsubvectorRq_ZBG=subvectorRq_ZBG(1:1,10:-1:1);
hold on
plot(subvectorRq_ZBG,'k*');
plot20(L,zl1u,invsubvectorRq_ZBG,'k')
hold off
[subAbscisa,subvectorRq_ZCG]=subconjunto21(Abscisa,vectorRq_ZCG);
subAbscisa(:,11)=[];
subvectorRq_ZCG(:,11)=[];
invsubvectorRq_ZCG=subvectorRq_ZCG(1:1,10:-1:1);
hold on
% Colores
c=plot(invsubvectorRq_ZCG,'*');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,invsubvectorRq_ZCG,'marron')
hold off
RXcadenaq=sprintf('Proteccion: %s',proteccion);
RXcadenampqc=sprintf('%s,%s',cadenampqc,RXcadenaq);
title(RXcadenampqc)
legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG')
end % end del switch extr de Salida 5
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.
