clc
variasejecuciones=1;
while variasejecuciones==1
clear
% --------------------------------------------------
fprintf('PG44ic\n');
fprintf('Esquema Linea Doble SE\n');
fprintf('Lineas Separadas ER y ET\n');
fprintf('Unominal kV: 400,220,132,66\n');
fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf('Compensacion 21: Si,No\n');
fprintf('Eleccion: kE,L1,L2,L3,Sccs,Sccr,Scct,Kr0,Ks0\n');
fprintf('Eleccion: k0=zt0/zs0\n');
fprintf('Falta Doble en SE, m=p=0...1\n');
fprintf('Analisis Sm,Rm\n');
%fprintf('Salida Grafica 1: Sm,Rm(2 Subpl),Unidad Pref.(pu)\n');
%fprintf(' : Sp,Tp(2 Subpl),Unidad Pref.(pu)\n');
%fprintf('Salida Grafica 2: Sm,Rm(2 Subpl),Unidad Pref.(pua)\n');
%fprintf(' : Sp,Tp(2 Subpl),Unidad Pref.(pua)\n');
%fprintf('Salida Grafica 3: Analisis Sm, Todas las Unidades (pua)\n');
%fprintf(' : Analisis Rm, Todas las Unidades (pua)\n');
%fprintf(' : Analisis Sp, Todas las Unidades (pua)\n');
%fprintf(' : Analisis Tp, 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 Tp, 67N\n');
%fprintf('Salida Grafica 5: Analisis Sm, R-X\n');
%fprintf(' : Analisis Rm, R-X\n');
%fprintf(' : Analisis Sp, R-X\n');
%fprintf(' : Analisis Tp, R-X\n');
fprintf('\n');
% --------------------------------------------------
%PG44ic
datos=input('datos: standard=1,input=2: ');
if datos==1
Un=220;
kE=1;
L1=200;
L2=100;
%L3=10;
L3=50;
compens21=2;
Sccs=1000;
Sccr=1000;
Scct=5000;
Ks0=1;
Kr0=1;
% --------------------------------------------------
else
Un=input('Un kV: 400,220,132,66 :');
kE=input('kE=E/Enominal: ');
L1=input('L1 km: ');
L2=input('L2 km: ');
L3=input('L3 km: ');
compens21=input('Compensacion 21: Si=1, No=2 :');
Sccs=input('Sccs(MVA): ');
Sccr=input('Sccr(MVA): ');
Scct=input('Scct(MVA): ');
Ks0=input('Ks0: ');
Kr0=input('Kr0: ');
end
% -----------------------------------------------------------------
a=0.5; b=0.5; c=0.5; d=0.5;
% -----------------------------------------------------------------
pfaltam=6;
pfaltap=8;
% -----------------------------------------------------------------
k0=input('k0=zt0/zs0: ');
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
% -----------------------------------------------------------------
% DATOS
% -----------------------------------------------------------------
Ubase=Un; % kV
Vbase=Ubase/(sqrt(3)); % kV
MVAbase=100;
Ibase=MVAbase/(3*Vbase); % kA
zbase=Vbase/Ibase; % Ohm
if Un==400
zl1u=(0.0346+0.3208i)/zbase; zl2u=zl1u; zl0u=(0.1828+1.1482i)/zbase;
zlM0u=(0.1482+0.7383i)/zbase;
elseif Un==220
zl1u=(0.0829+0.4135i)/zbase; zl2u=zl1u; zl0u=(0.2311+1.3085i)/zbase;
zlM0u=(0.1482+0.8031i)/zbase;
elseif Un==132
zl1u=(0.1342+0.4042i)/zbase; zl2u=zl1u; zl0u=(0.2824+1.3718i)/zbase;
zlM0u=(0.1482+0.8947i)/zbase;
elseif Un==66
zl1u=(0.2156+0.3952i)/zbase; zl2u=zl1u; zl0u=(0.3636+1.4278i)/zbase;
zlM0u=(0.1482+0.9288i)/zbase;
end
% ---------------------------------------------------
% ---------------------------------------------------
%CALCULOS PRELIMINARES
% ---------------------------------------------------
U=kE*Un; % kV
E=U/sqrt(3)+0i; % kV
E=E/Vbase; % E pu
k=(zl0u-zl1u)/(3*zl1u);
Km=zlM0u/zl1u;
zs1=(MVAbase/Sccs)*i; % pu
zr1=(MVAbase/Sccr)*i; % pu
zt1=(MVAbase/Scct)*i; % pu
zs2=zs1; % pu
zr2=zr1; % pu
zt2=zt1; % pu
zs0=Ks0*zs1; % pu
zr0=Kr0*zr1; % pu
zt0=k0*zs0;
% ----------------------------------------------------
% CALCULOS
% -----------------------------------------------------
i=1;
m=0.01;
while m<1;
p=m;
a=m;
c=p;
[Minc]=subMinc(1);
[Minca]=subMinca(1);
[z0]=subz0(zs0,zr0,zt0,zl0u,zlM0u,L1,L2,L3,a,b,c,d); % pu
[z1]=subz1(zs1,zr1,zt1,zl1u,L1,L2,L3,a,b,c,d); % pu
[z2]=subz2(zs2,zr2,zt2,zl2u,L1,L2,L3,a,b,c,d); % pu
% MATRIZ DE ADMITANCIA DE BARRAS
[Y0]=subY0(Minc,z0); % pu
[Y1]=subY1(Minc,z1); % pu
[Y2]=subY2(Minc,z2); % pu
% ZBUS(0,1,2)
Z0=inv(Y0); % pu
Z1=inv(Y1); % pu
Z2=inv(Y2); % pu
% CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p DE FALTA
[Im0,Im1,Im2,Vm0,Vm1,Vm2,Ip0,Ip1,Ip2,Vp0,Vp1,Vp2]=...
mag4mp(Z0,Z1,Z2,E,faltam,faltap,Rfm,Rgm,Rfp,Rgp,pfaltam,pfaltap);
if faltap==1 cfaltap='ag';elseif faltap==2 cfaltap='bg';elseif faltap==3 cfaltap='cg';
elseif faltap==4 cfaltap='ab';elseif faltap==5 cfaltap='bc';elseif faltap==6 cfaltap='ca';
elseif faltap==7 cfaltap='abg';elseif faltap==8 cfaltap='bcg';elseif faltap==9 cfaltap='cag';
elseif faltap==10 cfaltap='abc';elseif faltap==11 cfaltap='abcg';
end
% IBUS(0,1,2)
I0=zeros(9,1);
I1=zeros(9,1);
I2=zeros(9,1);
I0(pfaltam,1)=-Im0;
I1(pfaltam,1)=-Im1;
I2(pfaltam,1)=-Im2;
I0(pfaltap,1)=-Ip0;
I1(pfaltap,1)=-Ip1;
I2(pfaltap,1)=-Ip2;
% VBUS(0,1,2)
V0=Z0*I0;
V1=Z1*I1+E; V1a=[V1;E];
V2=Z2*I2;
% I(0,1,2) DE LOS ELEMENTOS
[Ie0]=subIe0(Minc,V0,z0);
[Ie1]=subIe1(Minca,V1a,z1);
[Ie2]=subIe2(Minc,V2,z2);
% ASIGNACIONES V,I
Vs0=V0(1);Vs1=V1(1);Vs2=V2(1);
Vr0=V0(2);Vr1=V1(2);Vr2=V2(2);
Vt0=V0(3);Vt1=V1(3);Vt2=V2(3);
Is0=Ie0(1);Is1=Ie1(1);Is2=Ie2(1);
Ir0=Ie0(2);Ir1=Ie1(2);Ir2=Ie2(2);
It0=Ie0(3);It1=Ie1(3);It2=Ie2(3);
Isr0=Ie0(4);Isr1=Ie1(4);Isr2=Ie2(4);
Ist0=Ie0(5);Ist1=Ie1(5);Ist2=Ie2(5);
Irs0=Ie0(6);Irs1=Ie1(6);Irs2=Ie2(6);
Its0=Ie0(7);Its1=Ie1(7);Its2=Ie2(7);
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);
[Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2);
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2);
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[Vta,Vtb,Vtc]=abc(Vt0,Vt1,Vt2);
[Isa,Isb,Isc]=abc(Is0,Is1,Is2);
[Ira,Irb,Irc]=abc(Ir0,Ir1,Ir2);
[Ita,Itb,Itc]=abc(It0,It1,It2);
[Isra,Isrb,Isrc]=abc(Isr0,Isr1,Isr2);
[Ista,Istb,Istc]=abc(Ist0,Ist1,Ist2);
[Irsa,Irsb,Irsc]=abc(Irs0,Irs1,Irs2);
[Itsa,Itsb,Itsc]=abc(Its0,Its1,Its2);
% PROTECCION DE DISTANCIA (21)
switch(compens21)
case{1},
% COMPENSADA
[Sr_ZAB,Sr_ZBC,Sr_ZCA,Sr_ZAG,Sr_ZBG,Sr_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isra,Isrb,Isrc,Isr0,k,Km,Ist0);
[St_ZAB,St_ZBC,St_ZCA,St_ZAG,St_ZBG,St_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Ista,Istb,Istc,Ist0,k,Km,Isr0);
%[Rs_ZAB,Rs_ZBC,Rs_ZCA,Rs_ZAG,Rs_ZBG,Rs_ZCG]=...
% PDistComp(Vra,Vrb,Vrc,Irsa,Irsb,Irsc,Irs0,k,Km,Its0);
%[Ts_ZAB,Ts_ZBC,Ts_ZCA,Ts_ZAG,Ts_ZBG,Ts_ZCG]=...
% PDistComp(Vta,Vtb,Vtc,Itsa,Itsb,Itsc,Its0,k,Km,Irs0);
[Rs_ZAB,Rs_ZBC,Rs_ZCA,Rs_ZAG,Rs_ZBG,Rs_ZCG]=...
PDist(Vra,Vrb,Vrc,Irsa,Irsb,Irsc,Irs0,k);
[Ts_ZAB,Ts_ZBC,Ts_ZCA,Ts_ZAG,Ts_ZBG,Ts_ZCG]=...
PDist(Vta,Vtb,Vtc,Itsa,Itsb,Itsc,Its0,k);
case{2},
% NO COMPENSADA
[Sr_ZAB,Sr_ZBC,Sr_ZCA,Sr_ZAG,Sr_ZBG,Sr_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isra,Isrb,Isrc,Isr0,k);
[St_ZAB,St_ZBC,St_ZCA,St_ZAG,St_ZBG,St_ZCG]=...
PDist(Vsa,Vsb,Vsc,Ista,Istb,Istc,Ist0,k);
[Rs_ZAB,Rs_ZBC,Rs_ZCA,Rs_ZAG,Rs_ZBG,Rs_ZCG]=...
PDist(Vra,Vrb,Vrc,Irsa,Irsb,Irsc,Irs0,k);
[Ts_ZAB,Ts_ZBC,Ts_ZCA,Ts_ZAG,Ts_ZBG,Ts_ZCG]=...
PDist(Vta,Vtb,Vtc,Itsa,Itsb,Itsc,Its0,k);
end
% 67N
[Vspol,Isrop,AngDirSr]=PDire(Vs0,Isr0); VVspol=Vspol*Vbase;IIsrop=Isrop*Ibase;
[Vspol,Istop,AngDirSt]=PDire(Vs0,Ist0); IIstop=Istop*Ibase;
[Vrpol,Irsop,AngDirRs]=PDire(Vr0,Irs0); VVrpol=Vrpol*Vbase;IIrsop=Irsop*Ibase;
[Vtpol,Itsop,AngDirTs]=PDire(Vt0,Its0); VVtpol=Vtpol*Vbase;IItsop=Itsop*Ibase;
% ******************************************************************
if faltam==1
cfaltam='ag';
alcanceSm=abs(Sr_ZAG)/cos(angle(Sr_ZAG)-angle(zl1u));
alcanceRm=abs(Rs_ZAG)/cos(angle(Rs_ZAG)-angle(zl1u));
elementom='AG';
elseif faltam==2
cfaltam='bg';
alcanceSm=abs(Sr_ZBG)/cos(angle(Sr_ZBG)-angle(zl1u));
alcanceRm=abs(Rs_ZBG)/cos(angle(Rs_ZBG)-angle(zl1u));
elementom='BG';
elseif faltam==3
cfaltam='cg';
alcanceSm=abs(Sr_ZCG)/cos(angle(Sr_ZCG)-angle(zl1u));
alcanceRm=abs(Rs_ZCG)/cos(angle(Rs_ZCG)-angle(zl1u));
elementom='CG';
elseif faltam==4
cfaltam='ab';
alcanceSm=abs(Sr_ZAB)/cos(angle(Sr_ZAB)-angle(zl1u));
alcanceRm=abs(Rs_ZAB)/cos(angle(Rs_ZAB)-angle(zl1u));
elementom='AB';
elseif faltam==5
cfaltam='bc';
alcanceSm=abs(Sr_ZBC)/cos(angle(Sr_ZBC)-angle(zl1u));
alcanceRm=abs(Rs_ZBC)/cos(angle(Rs_ZBC)-angle(zl1u));
elementom='BC';
elseif faltam==6
cfaltam='ca';
alcanceSm=abs(Sr_ZCA)/cos(angle(Sr_ZCA)-angle(zl1u));
alcanceRm=abs(Rs_ZCA)/cos(angle(Rs_ZCA)-angle(zl1u));
elementom='CA';
elseif faltam==7
cfaltam='abg';
alcanceSm=abs(Sr_ZAB)/cos(angle(Sr_ZAB)-angle(zl1u));
alcanceRm=abs(Rs_ZAB)/cos(angle(Rs_ZAB)-angle(zl1u));
elementom='AB';
elseif faltam==8
cfaltam='bcg';
alcanceSm=abs(Sr_ZBC)/cos(angle(Sr_ZBC)-angle(zl1u));
alcanceRm=abs(Rs_ZBC)/cos(angle(Rs_ZBC)-angle(zl1u));
elementom='BC';
elseif faltam==9
cfaltam='cag';
alcanceSm=abs(Sr_ZCA)/cos(angle(Sr_ZCA)-angle(zl1u));
alcanceRm=abs(Rs_ZCA)/cos(angle(Rs_ZCA)-angle(zl1u));
elementom='CA';
elseif faltam==10
cfaltam='abc';
alcanceSm=abs(Sr_ZBC)/cos(angle(Sr_ZBC)-angle(zl1u));
alcanceRm=abs(Rs_ZBC)/cos(angle(Rs_ZBC)-angle(zl1u));
elementom='BC';
elseif faltam==11
cfaltam='abcg';
alcanceSm=abs(Sr_ZBC)/cos(angle(Sr_ZBC)-angle(zl1u));
alcanceRm=abs(Rs_ZBC)/cos(angle(Rs_ZBC)-angle(zl1u));
elementom='BC';
end
if faltap==1
cfaltap='ag';
alcanceSp=abs(St_ZAG)/cos(angle(St_ZAG)-angle(zl1u));
alcanceTp=abs(Ts_ZAG)/cos(angle(Ts_ZAG)-angle(zl1u));
elementop='AG';
elseif faltap==2
cfaltap='bg';
alcanceSp=abs(St_ZBG)/cos(angle(St_ZBG)-angle(zl1u));
alcanceTp=abs(Ts_ZBG)/cos(angle(Ts_ZBG)-angle(zl1u));
elementop='BG';
elseif faltap==3
cfaltap='cg';
alcanceSp=abs(St_ZCG)/cos(angle(St_ZCG)-angle(zl1u));
alcanceTp=abs(Ts_ZCG)/cos(angle(Ts_ZCG)-angle(zl1u));
elementop='CG';
elseif faltap==4
cfaltap='ab';
alcanceSp=abs(St_ZAB)/cos(angle(St_ZAB)-angle(zl1u));
alcanceTp=abs(Ts_ZAB)/cos(angle(Ts_ZAB)-angle(zl1u));
elementop='AB';
elseif faltap==5
cfaltap='bc';
alcanceSp=abs(St_ZBC)/cos(angle(St_ZBC)-angle(zl1u));
alcanceTp=abs(Ts_ZBC)/cos(angle(Ts_ZBC)-angle(zl1u));
elementop='BC';
elseif faltap==6
cfaltap='ca';
alcanceSp=abs(St_ZCA)/cos(angle(St_ZCA)-angle(zl1u));
alcanceTp=abs(Ts_ZCA)/cos(angle(Ts_ZCA)-angle(zl1u));
elementop='CA';
elseif faltap==7
cfaltap='abg';
alcanceSp=abs(St_ZAB)/cos(angle(St_ZAB)-angle(zl1u));
alcanceTp=abs(Ts_ZAB)/cos(angle(Ts_ZAB)-angle(zl1u));
elementop='AB';
elseif faltap==8
cfaltap='bcg';
alcanceSp=abs(St_ZBC)/cos(angle(St_ZBC)-angle(zl1u));
alcanceTp=abs(Ts_ZBC)/cos(angle(Ts_ZBC)-angle(zl1u));
elementop='BC';
elseif faltap==9
cfaltap='cag';
alcanceSp=abs(St_ZCA)/cos(angle(St_ZCA)-angle(zl1u));
alcanceTp=abs(Ts_ZCA)/cos(angle(Ts_ZCA)-angle(zl1u));
elementop='CA';
elseif faltap==10
cfaltap='abc';
alcanceSp=abs(St_ZBC)/cos(angle(St_ZBC)-angle(zl1u));
alcanceTp=abs(Ts_ZBC)/cos(angle(Ts_ZBC)-angle(zl1u));
elementop='BC';
elseif faltap==11
cfaltap='abcg';
alcanceSp=abs(St_ZBC)/cos(angle(St_ZBC)-angle(zl1u));
alcanceTp=abs(Ts_ZBC)/cos(angle(Ts_ZBC)-angle(zl1u));
elementop='BC';
end
alcanceSmpu=alcanceSm/abs(zl1u*m*L1);
alcanceRmpu=alcanceRm/abs(zl1u*(L2+(1-m)*L1));
alcanceSmpua=alcanceSm/abs(zl1u*(L1+L2));
alcanceRmpua=alcanceRm/abs(zl1u*(L1+L2));
alcanceSppu=alcanceSp/abs(zl1u*p*L1);
alcanceTppu=alcanceTp/abs(zl1u*(L3+(1-p)*L1));
alcanceSppua=alcanceSp/abs(zl1u*(L1+L3));
alcanceTppua=alcanceTp/abs(zl1u*(L1+L3));
Abscisa(1,i)=m;
OrdenadaSm(1,i)=alcanceSmpu;
OrdenadaRm(1,i)=alcanceRmpu;
OrdenadaSma(1,i)=alcanceSmpua;
OrdenadaRma(1,i)=alcanceRmpua;
OrdenadaSp(1,i)=alcanceSppu;
OrdenadaTp(1,i)=alcanceTppu;
OrdenadaSpa(1,i)=alcanceSppua;
OrdenadaTpa(1,i)=alcanceTppua;
alcanceSm_ZAB=abs(Sr_ZAB)/cos(angle(Sr_ZAB)-angle(zl1u));
alcanceSm_ZBC=abs(Sr_ZBC)/cos(angle(Sr_ZBC)-angle(zl1u));
alcanceSm_ZCA=abs(Sr_ZCA)/cos(angle(Sr_ZCA)-angle(zl1u));
alcanceSm_ZAG=abs(Sr_ZAG)/cos(angle(Sr_ZAG)-angle(zl1u));
alcanceSm_ZBG=abs(Sr_ZBG)/cos(angle(Sr_ZBG)-angle(zl1u));
alcanceSm_ZCG=abs(Sr_ZCG)/cos(angle(Sr_ZCG)-angle(zl1u));
alcanceSp_ZAB=abs(St_ZAB)/cos(angle(St_ZAB)-angle(zl1u));
alcanceSp_ZBC=abs(St_ZBC)/cos(angle(St_ZBC)-angle(zl1u));
alcanceSp_ZCA=abs(St_ZCA)/cos(angle(St_ZCA)-angle(zl1u));
alcanceSp_ZAG=abs(St_ZAG)/cos(angle(St_ZAG)-angle(zl1u));
alcanceSp_ZBG=abs(St_ZBG)/cos(angle(St_ZBG)-angle(zl1u));
alcanceSp_ZCG=abs(St_ZCG)/cos(angle(St_ZCG)-angle(zl1u));
alcanceRm_ZAB=abs(Rs_ZAB)/cos(angle(Rs_ZAB)-angle(zl1u));
alcanceRm_ZBC=abs(Rs_ZBC)/cos(angle(Rs_ZBC)-angle(zl1u));
alcanceRm_ZCA=abs(Rs_ZCA)/cos(angle(Rs_ZCA)-angle(zl1u));
alcanceRm_ZAG=abs(Rs_ZAG)/cos(angle(Rs_ZAG)-angle(zl1u));
alcanceRm_ZBG=abs(Rs_ZBG)/cos(angle(Rs_ZBG)-angle(zl1u));
alcanceRm_ZCG=abs(Rs_ZCG)/cos(angle(Rs_ZCG)-angle(zl1u));
alcanceTp_ZAB=abs(Ts_ZAB)/cos(angle(Ts_ZAB)-angle(zl1u));
alcanceTp_ZBC=abs(Ts_ZBC)/cos(angle(Ts_ZBC)-angle(zl1u));
alcanceTp_ZCA=abs(Ts_ZCA)/cos(angle(Ts_ZCA)-angle(zl1u));
alcanceTp_ZAG=abs(Ts_ZAG)/cos(angle(Ts_ZAG)-angle(zl1u));
alcanceTp_ZBG=abs(Ts_ZBG)/cos(angle(Ts_ZBG)-angle(zl1u));
alcanceTp_ZCG=abs(Ts_ZCG)/cos(angle(Ts_ZCG)-angle(zl1u));
alcanceSm_ZABpua=alcanceSm_ZAB/abs(zl1u*(L1+L2));
alcanceSm_ZBCpua=alcanceSm_ZBC/abs(zl1u*(L1+L2));
alcanceSm_ZCApua=alcanceSm_ZCA/abs(zl1u*(L1+L2));
alcanceSm_ZAGpua=alcanceSm_ZAG/abs(zl1u*(L1+L2));
alcanceSm_ZBGpua=alcanceSm_ZBG/abs(zl1u*(L1+L2));
alcanceSm_ZCGpua=alcanceSm_ZCG/abs(zl1u*(L1+L2));
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;
alcanceRm_ZABpua=alcanceRm_ZAB/abs(zl1u*(L1+L2));
alcanceRm_ZBCpua=alcanceRm_ZBC/abs(zl1u*(L1+L2));
alcanceRm_ZCApua=alcanceRm_ZCA/abs(zl1u*(L1+L2));
alcanceRm_ZAGpua=alcanceRm_ZAG/abs(zl1u*(L1+L2));
alcanceRm_ZBGpua=alcanceRm_ZBG/abs(zl1u*(L1+L2));
alcanceRm_ZCGpua=alcanceRm_ZCG/abs(zl1u*(L1+L2));
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;
alcanceSp_ZABpua=alcanceSp_ZAB/abs(zl1u*(L1+L3));
alcanceSp_ZBCpua=alcanceSp_ZBC/abs(zl1u*(L1+L3));
alcanceSp_ZCApua=alcanceSp_ZCA/abs(zl1u*(L1+L3));
alcanceSp_ZAGpua=alcanceSp_ZAG/abs(zl1u*(L1+L3));
alcanceSp_ZBGpua=alcanceSp_ZBG/abs(zl1u*(L1+L3));
alcanceSp_ZCGpua=alcanceSp_ZCG/abs(zl1u*(L1+L3));
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;
alcanceTp_ZABpua=alcanceTp_ZAB/abs(zl1u*(L1+L3));
alcanceTp_ZBCpua=alcanceTp_ZBC/abs(zl1u*(L1+L3));
alcanceTp_ZCApua=alcanceTp_ZCA/abs(zl1u*(L1+L3));
alcanceTp_ZAGpua=alcanceTp_ZAG/abs(zl1u*(L1+L3));
alcanceTp_ZBGpua=alcanceTp_ZBG/abs(zl1u*(L1+L3));
alcanceTp_ZCGpua=alcanceTp_ZCG/abs(zl1u*(L1+L3));
OrdenadaTp_ZABpua(1,i)=alcanceTp_ZABpua;
OrdenadaTp_ZBCpua(1,i)=alcanceTp_ZBCpua;
OrdenadaTp_ZCApua(1,i)=alcanceTp_ZCApua;
OrdenadaTp_ZAGpua(1,i)=alcanceTp_ZAGpua;
OrdenadaTp_ZBGpua(1,i)=alcanceTp_ZBGpua;
OrdenadaTp_ZCGpua(1,i)=alcanceTp_ZCGpua;
%***************************************
%Para salida 4, 67N
%***************************************
vectorVspol(1,i)=VVspol;
vectorIsmop(1,i)=IIsrop;
vectorAngDirSm(1,i)=AngDirSr;
vectorVrpol(1,i)=VVrpol;
vectorIrmop(1,i)=IIrsop;
vectorAngDirRm(1,i)=AngDirRs;
vectorVspol(1,i)=VVspol;
vectorIspop(1,i)=IIstop;
vectorAngDirSp(1,i)=AngDirSt;
vectorVtpol(1,i)=VVtpol;
vectorItpop(1,i)=IItsop;
vectorAngDirTp(1,i)=AngDirTs;
%***************************************
%***************************************
% Para Salida 5 R-X
%***************************************
vectorSm_ZAB(1,i)=Sr_ZAB*zbase;
vectorSm_ZBC(1,i)=Sr_ZBC*zbase;
vectorSm_ZCA(1,i)=Sr_ZCA*zbase;
vectorSm_ZAG(1,i)=Sr_ZAG*zbase;
vectorSm_ZBG(1,i)=Sr_ZBG*zbase;
vectorSm_ZCG(1,i)=Sr_ZCG*zbase;
vectorRm_ZAB(1,i)=Rs_ZAB*zbase;
vectorRm_ZBC(1,i)=Rs_ZBC*zbase;
vectorRm_ZCA(1,i)=Rs_ZCA*zbase;
vectorRm_ZAG(1,i)=Rs_ZAG*zbase;
vectorRm_ZBG(1,i)=Rs_ZBG*zbase;
vectorRm_ZCG(1,i)=Rs_ZCG*zbase;
vectorSp_ZAB(1,i)=St_ZAB*zbase;
vectorSp_ZBC(1,i)=St_ZBC*zbase;
vectorSp_ZCA(1,i)=St_ZCA*zbase;
vectorSp_ZAG(1,i)=St_ZAG*zbase;
vectorSp_ZBG(1,i)=St_ZBG*zbase;
vectorSp_ZCG(1,i)=St_ZCG*zbase;
vectorTp_ZAB(1,i)=Ts_ZAB*zbase;
vectorTp_ZBC(1,i)=Ts_ZBC*zbase;
vectorTp_ZCA(1,i)=Ts_ZCA*zbase;
vectorTp_ZAG(1,i)=Ts_ZAG*zbase;
vectorTp_ZBG(1,i)=Ts_ZBG*zbase;
vectorTp_ZCG(1,i)=Ts_ZCG*zbase;
%***************************************
i=i+1;
m=m+0.01;
end
% -----------------------------------------------------------------------
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
% -----------------------------------------------------------------------
% Para Salidas 1 y 2
cadena0=sprintf('PG44ic,%s,%3.0fkV,kE=%-4.2f,L1=%3.0fkm,L2=%3.0fkm,L3=%3.0fkm',...
comp,Un,kE,L1,L2,L3);
cadena1=sprintf('Sccs=%-5.0f,Sccr=%-5.0f,Scct=%-5.0f\n',Sccs,Sccr,Scct);
cadena2=sprintf('Ks0=%-4.2f,Kr0=%-4.2f,k0=%-4.1f',Ks0,Kr0,k0);
cadena3=sprintf('SE(SR),m:%s,%s,Rfm=%-3.0f',cfaltam,elementom,Rfm);
cadenam=sprintf('%s,%s,%s,%s',cadena0,cadena1,cadena2,cadena3);
if faltam==7|faltam==8|faltam==9|faltam==11
cadena4=sprintf('Rgm=%-3.0f',Rgm);
cadenam=sprintf('%s,%s',cadenam,cadena4);
end
cadenap=sprintf('SE(ST),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);
% Para Salida 3 y 5
cadena0c=sprintf('PG44ic,%s,%3.0fkV,kE=%-4.2f,L1=%3.0fkm,L2=%3.0fkm,L3=%3.0fkm',...
comp,Un,kE,L1,L2,L3);
cadena1c=sprintf('Sccs=%-5.0f,Sccr=%-5.0f,Scct=%-5.0f\n',Sccs,Sccr,Scct);
cadena2c=sprintf('Ks0=%-4.2f,Kr0=%-4.2f,k0=%-4.1f',Ks0,Kr0,k0);
cadena3c=sprintf('SE(SR),m:%s,Rfm=%-3.0f',cfaltam,Rfm);
cadenamc=sprintf('%s,%s,%s,%s',cadena0c,cadena1c,cadena2c,cadena3c);
if faltam==7|faltam==8|faltam==9|faltam==11
cadena4c=sprintf('Rgm=%-3.0f',Rgm);
cadenamc=sprintf('%s,%s',cadenamc,cadena4c);
end
cadenapc=sprintf('SE(ST),p:%s,Rfp=%-3.0f',cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
cadenap1c=sprintf('Rgp=%-3.0f',Rgp);
cadenapc=sprintf('%s,%s',cadenapc,cadenap1c);
end
cadenampc=sprintf('%s,%s',cadenamc,cadenapc);
% Para Salida 4
dircadena0c=sprintf('PG44ic,67N,%3.0fkV,kE=%-4.2f,L1=%3.0fkm,L2=%3.0fkm,L3=%3.0fkm',...
Un,kE,L1,L2,L3);
dircadena1c=sprintf('Sccs=%-5.0f,Sccr=%-5.0f,Scct=%-5.0f\n',Sccs,Sccr,Scct);
dircadena2c=sprintf('Ks0=%-4.2f,Kr0=%-4.2f,k0=%-4.1f',Ks0,Kr0,k0);
dircadena3c=sprintf('SE(SR),m:%s,Rfm=%-3.0f',cfaltam,Rfm);
dircadenamc=sprintf('%s,%s,%s,%s',dircadena0c,dircadena1c,dircadena2c,dircadena3c);
if faltam==7|faltam==8|faltam==9|faltam==11
dircadena4c=sprintf('Rgm=%-3.0f',Rgm);
dircadenamc=sprintf('%s,%s',dircadenamc,dircadena4c);
end
dircadenapc=sprintf('SE(ST),p:%s,Rfp=%-3.0f',cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
dircadenap1c=sprintf('Rgp=%-3.0f',Rgp);
dircadenapc=sprintf('%s,%s',dircadenapc,dircadenap1c);
end
dircadenampc=sprintf('%s,%s',dircadenamc,dircadenapc);
% -----------------------------------------------------------------------
variassalidas=1;
while variassalidas==1
fprintf('\n');
fprintf('Salida Grafica 1: Sm,Rm(2 Subpl),Unidad Pref.(pu)\n');
fprintf(' : Sp,Tp(2 Subpl),Unidad Pref.(pu)\n');
fprintf('Salida Grafica 2: Sm,Rm(2 Subpl),Unidad Pref.(pua)\n');
fprintf(' : Sp,Tp(2 Subpl),Unidad Pref.(pua)\n');
fprintf('Salida Grafica 3: Analisis Sm, Todas las Unidades (pua)\n');
fprintf(' : Analisis Rm, Todas las Unidades (pua)\n');
fprintf(' : Analisis Sp, Todas las Unidades (pua)\n');
fprintf(' : Analisis Tp, 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 Tp, 67N\n');
fprintf('Salida Grafica 5: Analisis Sm, R-X\n');
fprintf(' : Analisis Rm, R-X\n');
fprintf(' : Analisis Sp, R-X\n');
fprintf(' : Analisis Tp, R-X\n');
fprintf('\n');
salida=input('Salida: 1,2,3,4,5 :');
switch(salida)
case{1}, % Salida 1
extremos=input('Extremos Sm,Rm=1, Sp,Tp=2: ');
switch(extremos)
case{1}, %Extremos Sm,Rm
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimSr(cont)=(0.8*(L1+L2))/(mm*L1);
LimRs(cont)=(0.8*(L1+L2))/(L2+(1-mm)*L1);
cobSm(cont)=(mm*L1)/(L1+L2);
cobRm(cont)=(L2+(1-mm)*L1)/(L1+L2);
cont=cont+1;
mm=mm+0.01;
end
Abscisa1m=Abscisa(1:1,1:(i-1));
Ordenada1m=OrdenadaSm(1:1,1:(i-1));
Abscisa4m=Abscisa(1:1,1:(i-1));
Ordenada4m=OrdenadaRm(1:1,1:(i-1));
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,LimSr,'g',Abscisa0,cobSm,'m',Abscisa1m,Ordenada1m,'r'),grid on
axis([0 1 0 3]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSmpu')
legend('Ref','Lim','Cob')
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,LimRs,'g',Abscisa0,cobRm,'m',Abscisa4m,Ordenada4m,'r'),grid on
axis([0 1 0 3]);
%title(cadenamp)
xlabel('m=p'),ylabel('alcanceRmpu')
legend('Ref','Lim','Cob')
case{2}, % Extremos Sp,Tp
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimSt(cont)=(0.8*(L1+L3))/(mm*L1);
LimTs(cont)=(0.8*(L1+L3))/(L3+(1-mm)*L1);
cobSp(cont)=(mm*L1)/(L1+L3);
cobTp(cont)=(L3+(1-mm)*L1)/(L1+L3);
cont=cont+1;
mm=mm+0.01;
end
Abscisa1p=Abscisa(1:1,1:(i-1));
Ordenada1p=OrdenadaSp(1:1,1:(i-1));
Abscisa4p=Abscisa(1:1,1:(i-1));
Ordenada4p=OrdenadaTp(1:1,1:(i-1));
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,LimSt,'g',Abscisa0,cobSp,'m',Abscisa1p,Ordenada1p,'r'),grid on
axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSppu')
legend('Ref','Lim','Cob')
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,LimTs,'g',Abscisa0,cobTp,'m',Abscisa4p,Ordenada4p,'r'),grid on
axis([0 1 0 2]);
%title(cadenamp)
xlabel('m=p'),ylabel('alcanceTppu')
legend('Ref','Lim','Cob')
end % fin switch extremos m
case(2) % Salida 2
extremos=input('Extremos Sm,Rm=1, Sp,Tp=2: ');
switch(extremos)
case{1}, % Extremos Sm,Rm
mm=0.01;
cont=1;
while mm<1;
Abscisam0(cont)=mm;
Ordenadam0(cont)=mm*(L1/(L1+L2));
Ordenadamz1(cont)=0.8;
cont=cont+1;
mm=mm+0.01;
end
mm=0.01;
cont=1;
while mm<1;
Abscisam00(cont)=mm;
Ordenadam00(cont)=1-mm*(L1/(L1+L2));
cont=cont+1;
mm=mm+0.01;
end
Abscisa1ma=Abscisa(1:1,1:(i-1));
Ordenada1ma=OrdenadaSma(1:1,1:(i-1));
Abscisa4ma=Abscisa(1:1,1:(i-1));
Ordenada4ma=OrdenadaRma(1:1,1:(i-1));
subplot(2,1,1)
plot(Abscisam0,Ordenadam0,'c',Abscisam0,Ordenadamz1,'m',Abscisa1ma,Ordenada1ma,'r'),grid on
axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSmpua')
legend('45','z1')
subplot(2,1,2)
plot(Abscisam00,Ordenadam00,'c',Abscisam0,Ordenadamz1,'m',Abscisa4ma,Ordenada4ma,'r'),grid on
axis([0 1 0 2]);
%title(cadenamp)
xlabel('m=p'),ylabel('alcanceRmpua')
legend('-45','z1')
case{2}, % Extremos Sp,Tp
mm=0.01;
cont=1;
while mm<1;
Abscisap0(cont)=mm;
Ordenadap0(cont)=mm*(L1/(L1+L3));
Ordenadapz1(cont)=0.8;
cont=cont+1;
mm=mm+0.01;
end
mm=0.01;
cont=1;
while mm<1;
Abscisap00(cont)=mm;
Ordenadap00(cont)=1-mm*(L1/(L1+L3));
cont=cont+1;
mm=mm+0.01;
end
Abscisa1pa=Abscisa(1:1,1:(i-1));
Ordenada1pa=OrdenadaSpa(1:1,1:(i-1));
Abscisa4pa=Abscisa(1:1,1:(i-1));
Ordenada4pa=OrdenadaTpa(1:1,1:(i-1));
subplot(2,1,1)
plot(Abscisap0,Ordenadap0,'c',Abscisap0,Ordenadapz1,'m',Abscisa1pa,Ordenada1pa,'r'),grid on
axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSppua')
legend('45','z1')
subplot(2,1,2)
plot(Abscisap00,Ordenadap00,'c',Abscisap0,Ordenadapz1,'m',Abscisa4pa,Ordenada4pa,'r'),grid on
axis([0 1 0 2]);
%title(cadenamp)
xlabel('m=p'),ylabel('alcanceTppua')
legend('-45','z1')
end % fin switch extremos p
case{3} % Salida 3
extremos=input('Extremos: m=1, p=2 :');
switch(extremos)
case{1}, % Extremos m
mm=0.02;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenada0a(cont)=mm*(L1/(L1+L2));
Ordenada00a(cont)=1-mm*(L1/(L1+L2));
cont=cont+1;
mm=mm+0.01;
end
extr=input('Extremo: Sm=1, Rm=2 :');
switch(extr)
case{1}, % 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])
% Fin Colores
axis([0 1 -2 2]);
title(cadenampc)
xlabel('m=p'),ylabel(cadylabel1a)
legend('45','AB','BC','CA','AG','BG','CG','z1')
case{2}, % Extremo Rm
cadylabel1a='alcRmpua';
% Colores
b=plot(Abscisa0,Ordenada00a);
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 0 2]);
title(cadenampc)
xlabel('m=p'),ylabel(cadylabel1a)
legend('-45','AB','BC','CA','AG','BG','CG','z1')
end
case{2}, % Extremos p
mm=0.02;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenada0a(cont)=mm*(L1/(L1+L3));
Ordenada00a(cont)=1-mm*(L1/(L1+L3));
cont=cont+1;
mm=mm+0.01;
end
extr=input('Extremo: Sp=1, Rp=2 :');
switch(extr)
case{1}, % 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(cadenampc)
xlabel('m=p'),ylabel(cadylabel1a)
legend('45','AB','BC','CA','AG','BG','CG','z1')
case{2}, % Extremo Tp
cadylabel1a='alcTppua';
% Colores
b=plot(Abscisa0,Ordenada00a);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaTp_ZABpua,'r',Abscisa,OrdenadaTp_ZBCpua,'g',Abscisa,OrdenadaTp_ZCApua,'b',...
Abscisa,OrdenadaTp_ZAGpua,'m',Abscisa,OrdenadaTp_ZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaTp_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 0 2]);
title(cadenampc)
xlabel('m=p'),ylabel(cadylabel1a)
legend('-45','AB','BC','CA','AG','BG','CG','z1')
end
end % end extremos
case{4} % salida 4 (Direccional)
extr=input('Extremo: Sm=1, Rm=2, Sp=3, Tp=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(dircadenampc)
rayar(1);
legend('AngDirSm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenampc)
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(dircadenampc)
rayar(1);
legend('AngDirRm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenampc)
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(dircadenampc)
rayar(1);
legend('AngDirSp','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenampc)
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,vectorAngDirTp,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirTp(grados)')
title(dircadenampc)
rayar(1);
legend('AngDirTp','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVtpol,'b');
grid
%title(dircadenampc)
xlabel('m=p'),ylabel('modulo Vtpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorItpop,'b');
grid
xlabel('m=p'),ylabel('modulo Itpop(kA)')
end % end del switch extr
case{5}, % ********* Salida 5 ************
extr=input('Extremo: Sm=1, Rm=2, Sp=3, Rp=4: ');
switch(extr)
case{1}
proteccion='Sm';
L=L1+L2;
pzona=0.8;
hold on
plot19b(L1,L2,L3,zl1u*zbase,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*zbase,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*zbase,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*zbase,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*zbase,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*zbase,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*zbase,subvectorSm_ZCG,'marron')
hold off
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenampc=sprintf('%s,%s',cadenampc,RXcadenap);
title(RXcadenampc)
legend('Mho','Linea','m=0..1','AB','BC','CA','AG','BG','CG')
case{2}
L=L1+L2;
proteccion='Rm';
pzona=0.8;
hold on
plot19b(L1,L2,L3,zl1u*zbase,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*zbase,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*zbase,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*zbase,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*zbase,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*zbase,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*zbase,invsubvectorRm_ZCG,'marron')
hold off
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenampc=sprintf('%s,%s',cadenampc,RXcadenap);
title(RXcadenampc)
legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG')
case{3}
L=L1+L3;
proteccion='Sp';
pzona=0.8;
hold on
plot19b(L1,L2,L3,zl1u*zbase,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*zbase,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*zbase,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*zbase,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*zbase,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*zbase,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*zbase,subvectorSp_ZCG,'marron')
hold off
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenampc=sprintf('%s,%s',cadenampc,RXcadenap);
title(RXcadenampc)
legend('Mho','Linea','m=0..1','AB','BC','CA','AG','BG','CG')
case{4}
L=L1+L3;
proteccion='Tp';
pzona=0.8;
hold on
plot19b(L1,L2,L3,zl1u*zbase,pzona,proteccion)
hold off
hold on
% Colores
a=plot(vectorTp_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorTp_ZBC,'m')
plot(vectorTp_ZCA,'c')
plot(vectorTp_ZAG,'r')
plot(vectorTp_ZBG,'k')
% Colores
c=plot(vectorTp_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
[subAbscisa,subvectorTp_ZAB]=subconjunto21(Abscisa,vectorTp_ZAB);
subAbscisa(:,11)=[];
subvectorTp_ZAB(:,11)=[];
invsubvectorTp_ZAB=subvectorTp_ZAB(1:1,10:-1:1);
hold on
b=plot(invsubvectorTp_ZAB,'*');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u*zbase,invsubvectorTp_ZAB,'naranja')
hold off
[subAbscisa,subvectorTp_ZBC]=subconjunto21(Abscisa,vectorTp_ZBC);
subAbscisa(:,11)=[];
subvectorTp_ZBC(:,11)=[];
invsubvectorTp_ZBC=subvectorTp_ZBC(1:1,10:-1:1);
hold on
plot(subvectorTp_ZBC,'m*');
plot20(L,zl1u*zbase,invsubvectorTp_ZBC,'m')
hold off
[subAbscisa,subvectorTp_ZCA]=subconjunto21(Abscisa,vectorTp_ZCA);
subAbscisa(:,11)=[];
subvectorTp_ZCA(:,11)=[];
invsubvectorTp_ZCA=subvectorTp_ZCA(1:1,10:-1:1);
hold on
plot(subvectorTp_ZCA,'c*');
plot20(L,zl1u*zbase,invsubvectorTp_ZCA,'c')
hold off
[subAbscisa,subvectorTp_ZAG]=subconjunto21(Abscisa,vectorTp_ZAG);
subAbscisa(:,11)=[];
subvectorTp_ZAG(:,11)=[];
invsubvectorTp_ZAG=subvectorTp_ZAG(1:1,10:-1:1);
hold on
plot(subvectorTp_ZAG,'r*');
plot20(L,zl1u*zbase,invsubvectorTp_ZAG,'r')
hold off
[subAbscisa,subvectorTp_ZBG]=subconjunto21(Abscisa,vectorTp_ZBG);
subAbscisa(:,11)=[];
subvectorTp_ZBG(:,11)=[];
invsubvectorTp_ZBG=subvectorTp_ZBG(1:1,10:-1:1);
hold on
plot(subvectorTp_ZBG,'k*');
plot20(L,zl1u*zbase,invsubvectorTp_ZBG,'k')
hold off
[subAbscisa,subvectorTp_ZCG]=subconjunto21(Abscisa,vectorTp_ZCG);
subAbscisa(:,11)=[];
subvectorTp_ZCG(:,11)=[];
invsubvectorTp_ZCG=subvectorTp_ZCG(1:1,10:-1:1);
hold on
% Colores
c=plot(invsubvectorTp_ZCG,'*');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u*zbase,invsubvectorTp_ZCG,'marron')
hold off
RXcadenap=sprintf('Proteccion: %s',proteccion);
RXcadenampc=sprintf('%s,%s',cadenampc,RXcadenap);
title(RXcadenampc)
legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG')
end % end del switch extr de Salida 5
end % end 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.
