% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG122c\n');
fprintf('Esquema Linea Doble SR\n');
fprintf('Unominal kV: 400,220,132,66\n');
fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf('Eleccion: kE,L,Sccs,Sccr,Kr0\n');
fprintf('Eleccion: k0=zs0/zr0\n');
fprintf('Impedancia de Transferencia Si,No\n');
fprintf('Eleccion:(zsr1,Ksr0)\n');
fprintf('Compensacion 21, Si, No\n');
fprintf('Falta Simple, m=0...1\n');
%fprintf('Salida Grafica 1: Analisis Sm, Rm, Unidad Preferente (pu)\n');
%fprintf('Salida Grafica 2: Analisis Sm, Rm, Unidad Preferente (pua)\n');
%fprintf('Salida Grafica 3: Analisis Sm, Todas las Unidades (pua)\n');
%fprintf(' Analisis Rm, Todas las Unidades (pua)\n');
%fprintf('Salida Grafica 4: Analisis Sm, 67N\n');
%fprintf(' Analisis Rm, 67N\n');
%fprintf(' Analisis Sr, 67N\n');
%fprintf(' Analisis Rs, 67N\n');
fprintf('\n');
% -------------------------------------------------------------------
%PG122c
datos=input('datos: standard=1,input=2: ');
% -------------------------------------------------------------------
if datos==1
Un=220;
kE=1;
L=220;
imptransfer=2;
compens21=2;
Sccs=5000;
Sccr=5000;
Kr0=1;
% -------------------------------------------------------------------
else
Un=input('U(Fase-Fase) kV:400,220,132,66: ');
kE=input('kE=E/Enominal: ');
L=input('Longitud Linea km: ');
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 :');
Sccs=input('Sccs(MVA): ');
Sccr=input('Sccr(MVA): ');
Kr0=input('Kr0=zr0/zr1: ');
end
% -------------------------------------------------------------------
k0=input('k0=zs0/zr0= ');
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
% -------------------------------------------------------------------
pfaltam=1;
% -------------------------------------------------------------------
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
% -------------------------------------------------------------------
En=Un/sqrt(3);
U=kE*Un;
E=kE*En;
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;
fprintf('\n');
Km=zlM0u/zl1u;
fprintf('\n');
% -------------------------------------------------------------------
i=1;
m=0.01;
while m<1;
switch(imptransfer)
case{1},
[Minc]=subMincT(1);
[z0]=subz0T(zs0,zr0,zsr0,zl0,zlM0,m);
[z1]=subz1T(zs1,zr1,zsr1,zl1,m);
[z2]=subz2T(zs2,zr2,zsr2,zl2,m);
case(2)
[Minc]=subMinc(1);
[z0]=subz0(zs0,zr0,zl0,zlM0,m);
[z1]=subz1(zs1,zr1,zl1,m);
[z2]=subz2(zs2,zr2,zl2,m);
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);
% V(0,1,2), I(0,1,2) EN EL PUNTO m DE LA FALTA
[Im0,Im1,Im2,Vm0,Vm1,Vm2]=mag2m(Z0,Z1,Z2,E,faltam,Rfm,Rgm,pfaltam);
% IBUS(0,1,2)
I0=zeros(3,1);
I1=zeros(3,1);
I2=zeros(3,1);
I0(pfaltam,1)=-Im0;
I1(pfaltam,1)=-Im1;
I2(pfaltam,1)=-Im2;
% 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);
% ASIGNACIONES V(0,1,2)
Vs0=V0(2); Vr0=V0(3);
Vs1=V1(2); Vr1=V1(3);
Vs2=V2(2); Vr2=V2(3);
% ASIGNACIONES I(0,1,2)
Ism0=Ie0(3); Irm0=Ie0(5); Isr0=Ie0(4); Irs0=-Isr0;
Ism1=Ie1(3); Irm1=Ie1(5); Isr1=Ie1(4); Irs1=-Isr1;
Ism2=Ie2(3); Irm2=Ie2(5); Isr2=Ie2(4); Irs2=-Isr2;
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2);
[Isra,Isrb,Isrc]=abc(Isr0,Isr1,Isr2);
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2);
[Irsa,Irsb,Irsc]=abc(Irs0,Irs1,Irs2);
% 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,Isr0);
[Sr_ZAB,Sr_ZBC,Sr_ZCA,Sr_ZAG,Sr_ZBG,Sr_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isra,Isrb,Isrc,Isr0,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,Irs0);
[Rs_ZAB,Rs_ZBC,Rs_ZCA,Rs_ZAG,Rs_ZBG,Rs_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irsa,Irsb,Irsc,Irs0,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);
[Sr_ZAB,Sr_ZBC,Sr_ZCA,Sr_ZAG,Sr_ZBG,Sr_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isra,Isrb,Isrc,Isr0,k);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k);
[Rs_ZAB,Rs_ZBC,Rs_ZCA,Rs_ZAG,Rs_ZBG,Rs_ZCG]=...
PDist(Vra,Vrb,Vrc,Irsa,Irsb,Irsc,Irs0,k);
end
[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0);
[Vspol,Isrop,AngDirSr]=PDire(Vs0,Isr0);
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0);
[Vrpol,Irsop,AngDirRs]=PDire(Vr0,Irs0);
%
% 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));
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));
% 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
alcanceSmpu=alcanceSm/abs(zl1u*m*L);
alcanceRmpu=alcanceRm/abs(zl1u*(1-m)*L);
alcanceSmpua=alcanceSm/abs(zl1u*L);
alcanceRmpua=alcanceRm/abs(zl1u*L);
Abscisa(1,i)=m;
OrdenadaSm(1,i)=alcanceSmpu;
OrdenadaRm(1,i)=alcanceRmpu;
OrdenadaSma(1,i)=alcanceSmpua;
OrdenadaRma(1,i)=alcanceRmpua;
% para Salida 3
alcanceSm_ZABpua=alcanceSm_ZAB/abs(zl1u*L);
alcanceSm_ZBCpua=alcanceSm_ZBC/abs(zl1u*L);
alcanceSm_ZCApua=alcanceSm_ZCA/abs(zl1u*L);
alcanceSm_ZAGpua=alcanceSm_ZAG/abs(zl1u*L);
alcanceSm_ZBGpua=alcanceSm_ZBG/abs(zl1u*L);
alcanceSm_ZCGpua=alcanceSm_ZCG/abs(zl1u*L);
OrdenadaSm_ZABpua(1,i)=alcanceSm_ZABpua;
OrdenadaSm_ZBCpua(1,i)=alcanceSm_ZBCpua;
OrdenadaSm_ZCApua(1,i)=alcanceSm_ZCApua;
OrdenadaSm_ZAGpua(1,i)=alcanceSm_ZAGpua;
OrdenadaSm_ZBGpua(1,i)=alcanceSm_ZBGpua;
OrdenadaSm_ZCGpua(1,i)=alcanceSm_ZCGpua;
alcanceRm_ZABpua=alcanceRm_ZAB/abs(zl1u*L);
alcanceRm_ZBCpua=alcanceRm_ZBC/abs(zl1u*L);
alcanceRm_ZCApua=alcanceRm_ZCA/abs(zl1u*L);
alcanceRm_ZAGpua=alcanceRm_ZAG/abs(zl1u*L);
alcanceRm_ZBGpua=alcanceRm_ZBG/abs(zl1u*L);
alcanceRm_ZCGpua=alcanceRm_ZCG/abs(zl1u*L);
OrdenadaRm_ZABpua(1,i)=alcanceRm_ZABpua;
OrdenadaRm_ZBCpua(1,i)=alcanceRm_ZBCpua;
OrdenadaRm_ZCApua(1,i)=alcanceRm_ZCApua;
OrdenadaRm_ZAGpua(1,i)=alcanceRm_ZAGpua;
OrdenadaRm_ZBGpua(1,i)=alcanceRm_ZBGpua;
OrdenadaRm_ZCGpua(1,i)=alcanceRm_ZCGpua;
%***************************************
%Para salida 4
%***************************************
vectorVspol(1,i)=Vspol;
vectorIsmop(1,i)=Ismop;
vectorAngDirSm(1,i)=AngDirSm;
%vectorcte1(1,i)=30;
%vectorcte2(1,i)=-150;
%vectorcte3(1,i)=-60;
vectorVrpol(1,i)=Vrpol;
vectorIrmop(1,i)=Irmop;
vectorAngDirRm(1,i)=AngDirRm;
vectorVspol(1,i)=Vspol;
vectorIsrop(1,i)=Isrop;
vectorAngDirSr(1,i)=AngDirSr;
vectorVrpol(1,i)=Vrpol;
vectorIrsop(1,i)=Irsop;
vectorAngDirRs(1,i)=AngDirRs;
%***************************************
i=i+1;
m=m+0.01;
end
% -------------------------------------------------------------------
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
% -------------------------------------------------------------------
cadenam1=sprintf('PG122c,%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
% -------------------------------------------------------------------
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',cadenam1,cadenamc2);
elseif faltam==7|faltam==8|faltam==9|faltam==11
cadenamc3=sprintf('Rgm=%-3.0f',Rgm);
cadenamc=sprintf('%s,%s,%s',cadenam1,cadenamc2,cadenamc3);
end
%***************************************
% Titulo para las gráficas 67N
%***************************************
dircadenam1=sprintf('PG122c,67N,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,\nKr0=%-5.2f,k0=%-4.1f',...
Un,kE,L,Sccs,Sccr,Kr0,k0);
if imptransfer==1
dircadenamT=sprintf('zsr1=%-6.2fi,Ksr0=%-4.2f',imag(zsr1),Ksr0);
dircadenam1=sprintf('%s,%s',dircadenam1,dircadenamT);
end
dircadenam2=sprintf('m:%s,Rfm=%-3.0f',cfaltam,Rfm);
if faltam==1|faltam==2|faltam==3|faltam==4|faltam==5|faltam==6|faltam==10
dircadenam=sprintf('%s,%s',dircadenam1,dircadenam2);
elseif faltam==7|faltam==8|faltam==9|faltam==11
dircadenam3=sprintf('Rgm=%-3.0f',Rgm);
dircadenam=sprintf('%s,%s,%s',dircadenam1,dircadenam2,dircadenam3);
end
%*************************************
variassalidas=1;
while variassalidas==1
fprintf('\n');
fprintf('Salida Grafica 1: Analisis Sm, Rm, Unidad Preferente (pu)\n');
fprintf('Salida Grafica 2: Analisis Sm, Rm, Unidad Preferente (pua)\n');
fprintf('Salida Grafica 3: Analisis Sm, Rm, Todas las Unidades (pua)\n');
fprintf('Salida Grafica 4: Analisis Sm, 67N\n');
fprintf(' : Analisis Rm, 67N\n');
fprintf(' : Analisis Sr, 67N\n');
fprintf(' : Analisis Rs, 67N\n');
fprintf('\n');
salida=input('Salida: 1,2,3,4: ');
switch(salida)
case{1}, % Salida 1: Sm, Rm, Unidad Preferente (pu)\n');
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
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'g',Abscisa,OrdenadaSm,'r'),grid on
axis([0 1 ejemenosy ejemasy]);
title(cadenam)
xlabel('m'),ylabel(cadenalabelSm)
legend('Ref','Lim',cadenalabelSm)
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadd,'g',Abscisa,OrdenadaRm,'r'),grid on
axis([0 1 ejemenosy ejemasy]);
%title(cadenam)
xlabel('m'),ylabel(cadenalabelRm)
legend('Ref','Lim',cadenalabelRm)
case(2)
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=mm;
Abscisa00(cont)=mm;
Ordenada00(cont)=1-mm;
Ordenadaz1(cont)=0.8;
cont=cont+1;
mm=mm+0.01;
end
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadaz1,'m',Abscisa,OrdenadaSma,'r'),grid on
title(cadenam)
xlabel('m'),ylabel(cadenalabelSma)
legend('Ref','z1',cadenalabelSma)
subplot(2,1,2)
plot(Abscisa00,Ordenada00,'c',Abscisa0,Ordenadaz1,'m',Abscisa,OrdenadaRma,'r'),grid on
%title(cadenam)
xlabel('m'),ylabel(cadenalabelRma)
legend('Ref','z1',cadenalabelRma)
case(3)
extr=input('Extremo: Sm=1, Rm=2: ');
switch(extr)
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
cadylabel1a='alcSmpua';
% Colores
b=plot(Abscisa0,Ordenada45);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaSm_ZABpua,'r',Abscisa,OrdenadaSm_ZBCpua,'g',Abscisa,OrdenadaSm_ZCApua,'b',...
Abscisa,OrdenadaSm_ZAGpua,'m',Abscisa,OrdenadaSm_ZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSm_ZCGpua);
set(a,'Color',[1 0.60 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 -5 5]);
%axis([0.8 1 0.8 1.2]);
title(cadenamc)
xlabel('m'),ylabel(cadylabel1a)
%legend('AB','BC','CA','AG','BG','CG','z1','ref')
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;
Ordenadamenos45(cont)=1-mm;
cont=cont+1;
mm=mm+0.01;
end
cadylabel1a='alcRmpua';
% Colores
b=plot(Abscisa0,Ordenadamenos45);
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.60 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 -5 5]);
%axis([0.8 1 0.8 1.2]);
title(cadenamc)
xlabel('m'),ylabel(cadylabel1a)
%legend('AB','BC','CA','AG','BG','CG','z1','ref')
legend('Ref','AB','BC','CA','AG','BG','CG','z1')
end
case{4}
extr=input('Extremo: Sm=1, Rm=2, Sr=3, Rs=4: ');
switch(extr)
case{1}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSm,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirSm(grados)')
title(dircadenam)
rayar(1);
legend('AngDirSm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIsmop,'b');
grid
xlabel('m'),ylabel('modulo Ismop(kA)')
case{2}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRm,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirRm(grados)')
title(dircadenam)
rayar(1);
legend('AngDirRm','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIrmop,'b');
grid
xlabel('m'),ylabel('modulo Irmop(kA)')
case{3}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSr,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirSr(grados)')
title(dircadenam)
rayar(1);
legend('AngDirSr','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIsrop,'b');
grid
xlabel('m'),ylabel('modulo Isrop(kA)')
case{4}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRs,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirRs(grados)')
title(dircadenam)
rayar(1);
legend('AngDirRs','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIrsop,'b');
grid
xlabel('m'),ylabel('modulo Irsop(kA)')
end % end del switch extr
end % end del switch salida
variassalidas=input('Mas Salidas=1,Acabar=2 :');
end % end del while variassalidas
variasejecuciones=input('Mas Ejecuciones=1,Acabar=2 :');
end % end del while variasejecuciones
fprintf(1,'Pulse Enter\n');
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
