clc
variasejecuciones=1;
while variasejecuciones==1
clear
% ---------------------------------------------------
fprintf('\n');
fprintf('PG54ic\n');
fprintf('Esquema Linea Doble SR, Bifurcando a X,Y\n');
fprintf('5 Tramos\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,L1,L2,L3,L4,L5,Sccs,Sccx,Sccr,Sccy\n');
fprintf(' Kx0,Kr0,Ky0,k0=zs0/zr0\n');
fprintf('Compensacion 21: Si,No\n');
fprintf('Falta Doble m=p\n');
fprintf('Eleccion del Tramo en Falta\n');
%fprintf('Salida Grafica 1: alcance pu (Extremos Faltas m,p, Todas las Unidades)\n');
%fprintf('Salida Grafica 2: alcance pua, Todas las Unidades: Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys\n');
%fprintf('Salida Grafica 3: 67N, AngDir,Vpol,Iop: Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys\n');
%fprintf('Salida Grafica 4: R-X, Todas las Unidades: Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys\n');
fprintf('\n');
fprintf('Para cada Tramo: Barra m=1, Barra p=2\n');
% ---------------------------------------------------
datos=input('datos: standard=1,input=2: ');
if datos==1
Un=220; % kV
kE=1;
L1=100; % km
L2=100;
L3=100;
L4=100;
L5=100;
compens21=2;
Sccs=5000; % MVA
Sccx=5000;
Sccr=5000;
Sccy=5000;
Kx0=1;
Kr0=1;
Ky0=1;
% ---------------------------------------------------
%Un=220; % kV
%kE=1.1;
%L1=100; % km
%L2=100;
%L3=10;
%L4=100;
%L5=100;
%compens21=2;
%Sccs=10000; % MVA
%Sccx=10000;
%Sccr=10000;
%Sccy=10000;
%Kx0=0.1;
%Kr0=0.1;
%Ky0=0.1;
% ---------------------------------------------------
else
Un=input('Unominal kV:400,220,132,66: ');
kE=input('kE=E/Enominal: ');
L1=input('L1(km)= ');
L2=input('L2(km)= ');
L3=input('L3(km)= ');
L4=input('L4(km)= ');
L5=input('L5(km)= ');
compens21=input('Compensacion 21:Si=1,No=2: ');
Sccs=input('Sccs(MVA)= ');
Sccx=input('Sccx(MVA)= ');
Sccr=input('Sccr(MVA)= ');
Sccy=input('Sccy(MVA)= ');
Kx0=input('Kx0= ');
Kr0=input('Kr0= ');
Ky0=input('Ky0= ');
end
% ---------------------------------------------------
k0=input('k0=zs0/zr0= ');
fprintf('Tramo en Falta: SE1:1, E1E2:2, E2R:3, E1X=4, E2Y:5 :\n');
tramo=input('Tramo: ');
switch(tramo)
case{1},
pfaltam=11;
pfaltap=20;
case{2},
pfaltam=14;
pfaltap=19;
case{3},
pfaltam=15;
pfaltap=16;
case{4},
pfaltam=12;
pfaltap=13;
case{5},
pfaltam=18;
pfaltap=17;
end % end del switch tramo
%-----------------------------------------------------------------------
if pfaltam==11
ppfaltam='SE1(SX)';
elseif pfaltam==12
ppfaltam='E1X(SX)';
elseif pfaltam==14
ppfaltam='E1E2(XR)';
elseif pfaltam==15
ppfaltam='E2R(XR)';
elseif pfaltam==18
ppfaltam='E2Y(SY)';
end
if pfaltap==13
ppfaltap='E1X(XR)';
elseif pfaltap==16
ppfaltap='E2R(YR)';
elseif pfaltap==17
ppfaltap='E2Y(RY)';
elseif pfaltap==19
ppfaltap='E1E2(SY)';
elseif pfaltap==20
ppfaltap='SE1(SY)';
end
%----------------------------------------------------------------------
fprintf('Falta en m, en %s: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11\n',ppfaltam);
faltam=input('Falta en m: ');
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
if faltam==1 cfaltam='ag';elseif faltam==2 cfaltam='bg';elseif faltam==3 cfaltam='cg';
elseif faltam==4 cfaltam='ab';elseif faltam==5 cfaltam='bc';elseif faltam==6 cfaltam='ca';
elseif faltam==7 cfaltam='abg';elseif faltam==8 cfaltam='bcg';elseif faltam==9 cfaltam='cag';
elseif faltam==10 cfaltam='abc';elseif faltam==11 cfaltam='abcg';
end
fprintf('\n');
fprintf('Falta en p, en %s: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11\n',ppfaltap);
faltap=input('Falta en p: ');
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
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
% ------------------------------------------------------------------
% DATOS
% ------------------------------------------------------------------
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
% ----------------------------------------------------------------------
a=0.5;b=0.5;c=0.5;d=0.5;e=0.5;f=0.5;g=0.5;h=0.5;q=0.5;r=0.5;
% ----------------------------------------------------------------------
%CALCULOS PRELIMINARES
% ----------------------------------------------------------------------
U=kE*Un; % kV
E=U/sqrt(3)+0i; % kV
k=(zl0u-zl1u)/(3*zl1u);
Km=zlM0u/zl1u;
zs1=(Un^2/Sccs)*i; zx1=(Un^2/Sccx)*i; zr1=(Un^2/Sccr)*i; zy1=(Un^2/Sccy)*i; % Ohm
zs2=zs1; zx2=zx1; zr2=zr1; zy2=zy1;
zx0=Kx0*zx1; zr0=Kr0*zr1; zy0=Ky0*zy1;
zs0=k0*zr0;
% ----------------------------------------------------------------------
% ----------------------------------------------------------------------
% CALCULOS
% ----------------------------------------------------------------------
[Minc]=subMinc(1);
i=1;
m=0.01;
while m<1;
p=m;
if pfaltam==11
a=m;
r=p;
elseif pfaltam==12
b=m;
c=p;
elseif pfaltam==14
d=m;
q=p;
elseif pfaltam==15
e=m;
f=p;
elseif pfaltam==18
h=m;
g=p;
end
[z0]=subz0(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r); % Ohm
[z1]=subz1(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r);
[z2]=subz2(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r);
% MATRIZ DE ADMITANCIA DE BARRAS
[Y0]=subY0(Minc,z0);
[Y1]=subY1(Minc,z1);
[Y2]=subY2(Minc,z2);
% ZBUS(0,1,2)
Z0=inv(Y0);
Z1=inv(Y1);
Z2=inv(Y2);
% CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p DE FALTA
[Im0,Im1,Im2,Vm0,Vm1,Vm2,Ip0,Ip1,Ip2,Vp0,Vp1,Vp2]=...
mag4mp(Z0,Z1,Z2,E,faltam,faltap,Rfm,Rgm,Rfp,Rgp,pfaltam,pfaltap); % kA, kV
% IBUS(0,1,2)
I0=zeros(20,1); I1=zeros(20,1); I2=zeros(20,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;
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,I
Vs0=V0(1);Vs1=V1(1);Vs2=V2(1);
Vx0=V0(2);Vx1=V1(2);Vx2=V2(2);
Vr0=V0(3);Vr1=V1(3);Vr2=V2(3);
Vy0=V0(4);Vy1=V1(4);Vy2=V2(4);
Is0=Ie0(1);Is1=Ie1(1);Is2=Ie2(1);
Ix0=Ie0(2);Ix1=Ie1(2);Ix2=Ie2(2);
Ir0=Ie0(3);Ir1=Ie1(3);Ir2=Ie2(3);
Iy0=Ie0(4);Iy1=Ie1(4);Iy2=Ie2(4);
Isx0=Ie0(5);Isx1=Ie1(5);Isx2=Ie2(5);
Isy0=Ie0(6);Isy1=Ie1(6);Isy2=Ie2(6);
Ixs0=Ie0(7);Ixs1=Ie1(7);Ixs2=Ie2(7);
Ixr0=Ie0(8);Ixr1=Ie1(8);Ixr2=Ie2(8);
Irx0=Ie0(9);Irx1=Ie1(9);Irx2=Ie2(9);
Iry0=Ie0(10);Iry1=Ie1(10);Iry2=Ie2(10);
Iys0=Ie0(11);Iys1=Ie1(11);Iys2=Ie2(11);
Iyr0=Ie0(12);Iyr1=Ie1(12);Iyr2=Ie2(12);
% CALCULO DE LAS COMPONENTES DE FASE DE I,V EN LOS PUNTOS m,p DE FALTA
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);
[Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2);
[Isa,Isb,Isc]=abc(Is0,Is1,Is2);
[Ixa,Ixb,Ixc]=abc(Ix0,Ix1,Ix2);
[Ira,Irb,Irc]=abc(Ir0,Ir1,Ir2);
[Iya,Iyb,Iyc]=abc(Iy0,Iy1,Iy2);
[Isxa,Isxb,Isxc]=abc(Isx0,Isx1,Isx2);
[Isya,Isyb,Isyc]=abc(Isy0,Isy1,Isy2);
[Ixsa,Ixsb,Ixsc]=abc(Ixs0,Ixs1,Ixs2);
[Ixra,Ixrb,Ixrc]=abc(Ixr0,Ixr1,Ixr2);
[Irxa,Irxb,Irxc]=abc(Irx0,Irx1,Irx2);
[Irya,Iryb,Iryc]=abc(Iry0,Iry1,Iry2);
[Iyra,Iyrb,Iyrc]=abc(Iyr0,Iyr1,Iyr2);
[Iysa,Iysb,Iysc]=abc(Iys0,Iys1,Iys2);
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2);
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vxa,Vxb,Vxc]=abc(Vx0,Vx1,Vx2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[Vya,Vyb,Vyc]=abc(Vy0,Vy1,Vy2);
% PROTECCION DE DISTANCIA (21)
switch(compens21)
case{1},
% COMPENSADA
[Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isxa,Isxb,Isxc,Isx0,k,Km,Isy0);
[Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isya,Isyb,Isyc,Isy0,k,Km,Isx0);
[Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG]=...
PDistComp(Vxa,Vxb,Vxc,Ixsa,Ixsb,Ixsc,Ixs0,k,Km,Ixr0);
[Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG]=...
PDistComp(Vxa,Vxb,Vxc,Ixra,Ixrb,Ixrc,Ixr0,k,Km,Ixs0);
[Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irxa,Irxb,Irxc,Irx0,k,Km,Iry0);
[Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irya,Iryb,Iryc,Iry0,k,Km,Irx0);
[Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG]=...
PDistComp(Vya,Vyb,Vyc,Iyra,Iyrb,Iyrc,Iyr0,k,Km,Iys0);
[Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG]=...
PDistComp(Vya,Vyb,Vyc,Iysa,Iysb,Iysc,Iys0,k,Km,Iyr0);
case{2},
% NO COMPENSADA
[Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isxa,Isxb,Isxc,Isx0,k);
[Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isya,Isyb,Isyc,Isy0,k);
[Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG]=...
PDist(Vxa,Vxb,Vxc,Ixsa,Ixsb,Ixsc,Ixs0,k);
[Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG]=...
PDist(Vxa,Vxb,Vxc,Ixra,Ixrb,Ixrc,Ixr0,k);
[Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG]=...
PDist(Vra,Vrb,Vrc,Irxa,Irxb,Irxc,Irx0,k);
[Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG]=...
PDist(Vra,Vrb,Vrc,Irya,Iryb,Iryc,Iry0,k);
[Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG]=...
PDist(Vya,Vyb,Vyc,Iyra,Iyrb,Iyrc,Iyr0,k);
[Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG]=...
PDist(Vya,Vyb,Vyc,Iysa,Iysb,Iysc,Iys0,k);
end
% 67N
[Vspol,Isxop,AngDirSx]=PDire(Vs0,Isx0);
[Vspol,Isyop,AngDirSy]=PDire(Vs0,Isy0);
[Vxpol,Ixsop,AngDirXs]=PDire(Vx0,Ixs0);
[Vxpol,Ixrop,AngDirXr]=PDire(Vx0,Ixr0);
[Vrpol,Irxop,AngDirRx]=PDire(Vr0,Irx0);
[Vrpol,Iryop,AngDirRy]=PDire(Vr0,Iry0);
[Vypol,Iyrop,AngDirYr]=PDire(Vy0,Iyr0);
[Vypol,Iysop,AngDirYs]=PDire(Vy0,Iys0);
%***************************************
switch(pfaltam)
case{11},
em1='Extremo S';
em2='Extremo X';
L=m*L1;
La=L1+L4;
[alcSxZABpu,alcSxZBCpu,alcSxZCApu,alcSxZAGpu,alcSxZBGpu,alcSxZCGpu]=...
alcpu(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,L);
L=L4+(1-m)*L1;
[alcXsZABpu,alcXsZBCpu,alcXsZCApu,alcXsZAGpu,alcXsZBGpu,alcXsZCGpu]=...
alcpu(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,L);
cadylabel1='alcSmpu';
cadylabel2='alcXmpu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=0;
bp=L1/La;
as=(L4+L1)/La;
bs=-L1/La;
case{12},
em1='Extremo S';
em2='Extremo X';
L=L1+m*L4;
La=L1+L4;
[alcSxZABpu,alcSxZBCpu,alcSxZCApu,alcSxZAGpu,alcSxZBGpu,alcSxZCGpu]=...
alcpu(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,L);
L=(1-m)*L4;
[alcXsZABpu,alcXsZBCpu,alcXsZCApu,alcXsZAGpu,alcXsZBGpu,alcXsZCGpu]=...
alcpu(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,L);
cadylabel1='alcSmpu';
cadylabel2='alcXmpu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L1/La;
bp=L4/La;
as=L4/La;
bs=-L4/La;
case{14},
em1='Extremo X';
em2='Extremo R';
L=L4+m*L2;
La=L4+L2+L3;
[alcXrZABpu,alcXrZBCpu,alcXrZCApu,alcXrZAGpu,alcXrZBGpu,alcXrZCGpu]=...
alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L);
L=L3+(1-m)*L2;
[alcRxZABpu,alcRxZBCpu,alcRxZCApu,alcRxZAGpu,alcRxZBGpu,alcRxZCGpu]=...
alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L);
cadylabel1='alcXmpu';
cadylabel2='alcRmpu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L4/La;
bp=L2/La;
as=(L3+L2)/La;
bs=-L2/La;
case{15},
em1='Extremo X';
em2='Extremo R';
L=L4+L2+m*L3;
La=L4+L2+L3;
[alcXrZABpu,alcXrZBCpu,alcXrZCApu,alcXrZAGpu,alcXrZBGpu,alcXrZCGpu]=...
alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L);
L=(1-m)*L3;
[alcRxZABpu,alcRxZBCpu,alcRxZCApu,alcRxZAGpu,alcRxZBGpu,alcRxZCGpu]=...
alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L);
cadylabel1='alcXmpu';
cadylabel2='alcRmpu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=(L4+L2)/La;
bp=L3/La;
as=L3/La;
bs=-L3/La;
case{18},
em1='Extremo S';
em2='Extremo Y';
L=L1+L2+m*L5;
La=L5+L2+L1;
[alcSyZABpu,alcSyZBCpu,alcSyZCApu,alcSyZAGpu,alcSyZBGpu,alcSyZCGpu]=...
alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L);
L=(1-m)*L5;
[alcYsZABpu,alcYsZBCpu,alcYsZCApu,alcYsZAGpu,alcYsZBGpu,alcYsZCGpu]=...
alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L);
cadylabel1='alcSmpu';
cadylabel2='alcYmpu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L5/La;
bp=-L5/La;
as=(L1+L2)/La;
bs=L5/La;
end
switch(pfaltap)
case{13},
em1='Extremo R';
em2='Extremo X';
L=L3+L2+p*L4;
La=L4+L2+L3;
[alcRxZABpu,alcRxZBCpu,alcRxZCApu,alcRxZAGpu,alcRxZBGpu,alcRxZCGpu]=...
alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L);
L=(1-p)*L4;
[alcXrZABpu,alcXrZBCpu,alcXrZCApu,alcXrZAGpu,alcXrZBGpu,alcXrZCGpu]=...
alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L);
cadylabel1='alcRppu';
cadylabel2='alcXppu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L4/La;
bp=-L4/La;
as=(L3+L2)/La;
bs=L4/La;
case{16},
em1='Extremo Y';
em2='Extremo R';
L=L5+p*L3;
La=L3+L5;
[alcYrZABpu,alcYrZBCpu,alcYrZCApu,alcYrZAGpu,alcYrZBGpu,alcYrZCGpu]=...
alcpu(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,L);
L=(1-p)*L3;
[alcRyZABpu,alcRyZBCpu,alcRyZCApu,alcRyZAGpu,alcRyZBGpu,alcRyZCGpu]=...
alcpu(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,L);
cadylabel1='alcYppu';
cadylabel2='alcRppu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L3/La;
bp=-L3/La;
as=L5/La;
bs=L3/La;
case{17},
em1='Extremo R';
em2='Extremo Y';
L=L3+p*L5;
La=L3+L5;
[alcRyZABpu,alcRyZBCpu,alcRyZCApu,alcRyZAGpu,alcRyZBGpu,alcRyZCGpu]=...
alcpu(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,L);
L=(1-p)*L5;
[alcYrZABpu,alcYrZBCpu,alcYrZCApu,alcYrZAGpu,alcYrZBGpu,alcYrZCGpu]=...
alcpu(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,L);
cadylabel1='alcRppu';
cadylabel2='alcYppu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=L3/La;
bp=L5/La;
as=L5/La;
bs=-L5/La;
case{19},
em1='Extremo S';
em2='Extremo Y';
L=L1+p*L2;
La=L5+L2+L1;
[alcSyZABpu,alcSyZBCpu,alcSyZCApu,alcSyZAGpu,alcSyZBGpu,alcSyZCGpu]=...
alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L);
L=L5+(1-p)*L2;
[alcYsZABpu,alcYsZBCpu,alcYsZCApu,alcYsZAGpu,alcYsZBGpu,alcYsZCGpu]=...
alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L);
cadylabel1='alcSppu';
cadylabel2='alcYppu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=(L5+L2)/La;
bp=-L2/La;
as=L1/La;
bs=L2/La;
case{20},
em1='Extremo S';
em2='Extremo Y';
L=p*L1;
La=L5+L2+L1;
[alcSyZABpu,alcSyZBCpu,alcSyZCApu,alcSyZAGpu,alcSyZBGpu,alcSyZCGpu]=...
alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L);
L=L5+L2+(1-p)*L1;
[alcYsZABpu,alcYsZBCpu,alcYsZCApu,alcYsZAGpu,alcYsZBGpu,alcYsZCGpu]=...
alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L);
cadylabel1='alcSppu';
cadylabel2='alcYppu';
% Cobertura Extremo Primero: cp=ap+bp*m;
% Cobertura Extremo Segundo: cs=as+bs*m;
ap=(L5+L2+L1)/La;
bp=-L1/La;
as=0;
bs=L1/La;
end
%***************************************
[alcSxZABpua,alcSxZBCpua,alcSxZCApua,alcSxZAGpua,alcSxZBGpua,alcSxZCGpua]=...
alcpua(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,L1+L4);
[alcSyZABpua,alcSyZBCpua,alcSyZCApua,alcSyZAGpua,alcSyZBGpua,alcSyZCGpua]=...
alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L1+L2+L5);
[alcXsZABpua,alcXsZBCpua,alcXsZCApua,alcXsZAGpua,alcXsZBGpua,alcXsZCGpua]=...
alcpua(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,L1+L4);
[alcXrZABpua,alcXrZBCpua,alcXrZCApua,alcXrZAGpua,alcXrZBGpua,alcXrZCGpua]=...
alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L4+L2+L3);
[alcRxZABpua,alcRxZBCpua,alcRxZCApua,alcRxZAGpua,alcRxZBGpua,alcRxZCGpua]=...
alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L4+L2+L3);
[alcRyZABpua,alcRyZBCpua,alcRyZCApua,alcRyZAGpua,alcRyZBGpua,alcRyZCGpua]=...
alcpua(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,L3+L5);
[alcYrZABpua,alcYrZBCpua,alcYrZCApua,alcYrZAGpua,alcYrZBGpua,alcYrZCGpua]=...
alcpua(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,L3+L5);
[alcYsZABpua,alcYsZBCpua,alcYsZCApua,alcYsZAGpua,alcYsZBGpua,alcYsZCGpua]=...
alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L1+L2+L5);
%***************************************
%Para salida 1, alcance pu
%***************************************
switch(tramo)
case {1},
vectalcSxZABpu(1,i)=alcSxZABpu;
vectalcSxZBCpu(1,i)=alcSxZBCpu;
vectalcSxZCApu(1,i)=alcSxZCApu;
vectalcSxZAGpu(1,i)=alcSxZAGpu;
vectalcSxZBGpu(1,i)=alcSxZBGpu;
vectalcSxZCGpu(1,i)=alcSxZCGpu;
vectalcXsZABpu(1,i)=alcXsZABpu;
vectalcXsZBCpu(1,i)=alcXsZBCpu;
vectalcXsZCApu(1,i)=alcXsZCApu;
vectalcXsZAGpu(1,i)=alcXsZAGpu;
vectalcXsZBGpu(1,i)=alcXsZBGpu;
vectalcXsZCGpu(1,i)=alcXsZCGpu;
vectalcSyZABpu(1,i)=alcSyZABpu;
vectalcSyZBCpu(1,i)=alcSyZBCpu;
vectalcSyZCApu(1,i)=alcSyZCApu;
vectalcSyZAGpu(1,i)=alcSyZAGpu;
vectalcSyZBGpu(1,i)=alcSyZBGpu;
vectalcSyZCGpu(1,i)=alcSyZCGpu;
vectalcYsZABpu(1,i)=alcYsZABpu;
vectalcYsZBCpu(1,i)=alcYsZBCpu;
vectalcYsZCApu(1,i)=alcYsZCApu;
vectalcYsZAGpu(1,i)=alcYsZAGpu;
vectalcYsZBGpu(1,i)=alcYsZBGpu;
vectalcYsZCGpu(1,i)=alcYsZCGpu;
case {2},
vectalcXrZABpu(1,i)=alcXrZABpu;
vectalcXrZBCpu(1,i)=alcXrZBCpu;
vectalcXrZCApu(1,i)=alcXrZCApu;
vectalcXrZAGpu(1,i)=alcXrZAGpu;
vectalcXrZBGpu(1,i)=alcXrZBGpu;
vectalcXrZCGpu(1,i)=alcXrZCGpu;
vectalcRxZABpu(1,i)=alcRxZABpu;
vectalcRxZBCpu(1,i)=alcRxZBCpu;
vectalcRxZCApu(1,i)=alcRxZCApu;
vectalcRxZAGpu(1,i)=alcRxZAGpu;
vectalcRxZBGpu(1,i)=alcRxZBGpu;
vectalcRxZCGpu(1,i)=alcRxZCGpu;
vectalcSyZABpu(1,i)=alcSyZABpu;
vectalcSyZBCpu(1,i)=alcSyZBCpu;
vectalcSyZCApu(1,i)=alcSyZCApu;
vectalcSyZAGpu(1,i)=alcSyZAGpu;
vectalcSyZBGpu(1,i)=alcSyZBGpu;
vectalcSyZCGpu(1,i)=alcSyZCGpu;
vectalcYsZABpu(1,i)=alcYsZABpu;
vectalcYsZBCpu(1,i)=alcYsZBCpu;
vectalcYsZCApu(1,i)=alcYsZCApu;
vectalcYsZAGpu(1,i)=alcYsZAGpu;
vectalcYsZBGpu(1,i)=alcYsZBGpu;
vectalcYsZCGpu(1,i)=alcYsZCGpu;
case {3},
vectalcXrZABpu(1,i)=alcXrZABpu;
vectalcXrZBCpu(1,i)=alcXrZBCpu;
vectalcXrZCApu(1,i)=alcXrZCApu;
vectalcXrZAGpu(1,i)=alcXrZAGpu;
vectalcXrZBGpu(1,i)=alcXrZBGpu;
vectalcXrZCGpu(1,i)=alcXrZCGpu;
vectalcRxZABpu(1,i)=alcRxZABpu;
vectalcRxZBCpu(1,i)=alcRxZBCpu;
vectalcRxZCApu(1,i)=alcRxZCApu;
vectalcRxZAGpu(1,i)=alcRxZAGpu;
vectalcRxZBGpu(1,i)=alcRxZBGpu;
vectalcRxZCGpu(1,i)=alcRxZCGpu;
vectalcRyZABpu(1,i)=alcRyZABpu;
vectalcRyZBCpu(1,i)=alcRyZBCpu;
vectalcRyZCApu(1,i)=alcRyZCApu;
vectalcRyZAGpu(1,i)=alcRyZAGpu;
vectalcRyZBGpu(1,i)=alcRyZBGpu;
vectalcRyZCGpu(1,i)=alcRyZCGpu;
vectalcYrZABpu(1,i)=alcYrZABpu;
vectalcYrZBCpu(1,i)=alcYrZBCpu;
vectalcYrZCApu(1,i)=alcYrZCApu;
vectalcYrZAGpu(1,i)=alcYrZAGpu;
vectalcYrZBGpu(1,i)=alcYrZBGpu;
vectalcYrZCGpu(1,i)=alcYrZCGpu;
case {4},
vectalcSxZABpu(1,i)=alcSxZABpu;
vectalcSxZBCpu(1,i)=alcSxZBCpu;
vectalcSxZCApu(1,i)=alcSxZCApu;
vectalcSxZAGpu(1,i)=alcSxZAGpu;
vectalcSxZBGpu(1,i)=alcSxZBGpu;
vectalcSxZCGpu(1,i)=alcSxZCGpu;
vectalcXsZABpu(1,i)=alcXsZABpu;
vectalcXsZBCpu(1,i)=alcXsZBCpu;
vectalcXsZCApu(1,i)=alcXsZCApu;
vectalcXsZAGpu(1,i)=alcXsZAGpu;
vectalcXsZBGpu(1,i)=alcXsZBGpu;
vectalcXsZCGpu(1,i)=alcXsZCGpu;
vectalcXrZABpu(1,i)=alcXrZABpu;
vectalcXrZBCpu(1,i)=alcXrZBCpu;
vectalcXrZCApu(1,i)=alcXrZCApu;
vectalcXrZAGpu(1,i)=alcXrZAGpu;
vectalcXrZBGpu(1,i)=alcXrZBGpu;
vectalcXrZCGpu(1,i)=alcXrZCGpu;
vectalcRxZABpu(1,i)=alcRxZABpu;
vectalcRxZBCpu(1,i)=alcRxZBCpu;
vectalcRxZCApu(1,i)=alcRxZCApu;
vectalcRxZAGpu(1,i)=alcRxZAGpu;
vectalcRxZBGpu(1,i)=alcRxZBGpu;
vectalcRxZCGpu(1,i)=alcRxZCGpu;
case {5},
vectalcSyZABpu(1,i)=alcSyZABpu;
vectalcSyZBCpu(1,i)=alcSyZBCpu;
vectalcSyZCApu(1,i)=alcSyZCApu;
vectalcSyZAGpu(1,i)=alcSyZAGpu;
vectalcSyZBGpu(1,i)=alcSyZBGpu;
vectalcSyZCGpu(1,i)=alcSyZCGpu;
vectalcYsZABpu(1,i)=alcYsZABpu;
vectalcYsZBCpu(1,i)=alcYsZBCpu;
vectalcYsZCApu(1,i)=alcYsZCApu;
vectalcYsZAGpu(1,i)=alcYsZAGpu;
vectalcYsZBGpu(1,i)=alcYsZBGpu;
vectalcYsZCGpu(1,i)=alcYsZCGpu;
vectalcRyZABpu(1,i)=alcRyZABpu;
vectalcRyZBCpu(1,i)=alcRyZBCpu;
vectalcRyZCApu(1,i)=alcRyZCApu;
vectalcRyZAGpu(1,i)=alcRyZAGpu;
vectalcRyZBGpu(1,i)=alcRyZBGpu;
vectalcRyZCGpu(1,i)=alcRyZCGpu;
vectalcYrZABpu(1,i)=alcYrZABpu;
vectalcYrZBCpu(1,i)=alcYrZBCpu;
vectalcYrZCApu(1,i)=alcYrZCApu;
vectalcYrZAGpu(1,i)=alcYrZAGpu;
vectalcYrZBGpu(1,i)=alcYrZBGpu;
vectalcYrZCGpu(1,i)=alcYrZCGpu;
end % end del switch tramo
%***************************************
%***************************************
%Para salida 2, alcance pua
%***************************************
vectalcSxZABpua(1,i)=alcSxZABpua;
vectalcSxZBCpua(1,i)=alcSxZBCpua;
vectalcSxZCApua(1,i)=alcSxZCApua;
vectalcSxZAGpua(1,i)=alcSxZAGpua;
vectalcSxZBGpua(1,i)=alcSxZBGpua;
vectalcSxZCGpua(1,i)=alcSxZCGpua;
vectalcSyZABpua(1,i)=alcSyZABpua;
vectalcSyZBCpua(1,i)=alcSyZBCpua;
vectalcSyZCApua(1,i)=alcSyZCApua;
vectalcSyZAGpua(1,i)=alcSyZAGpua;
vectalcSyZBGpua(1,i)=alcSyZBGpua;
vectalcSyZCGpua(1,i)=alcSyZCGpua;
vectalcXsZABpua(1,i)=alcXsZABpua;
vectalcXsZBCpua(1,i)=alcXsZBCpua;
vectalcXsZCApua(1,i)=alcXsZCApua;
vectalcXsZAGpua(1,i)=alcXsZAGpua;
vectalcXsZBGpua(1,i)=alcXsZBGpua;
vectalcXsZCGpua(1,i)=alcXsZCGpua;
vectalcXrZABpua(1,i)=alcXrZABpua;
vectalcXrZBCpua(1,i)=alcXrZBCpua;
vectalcXrZCApua(1,i)=alcXrZCApua;
vectalcXrZAGpua(1,i)=alcXrZAGpua;
vectalcXrZBGpua(1,i)=alcXrZBGpua;
vectalcXrZCGpua(1,i)=alcXrZCGpua;
vectalcRxZABpua(1,i)=alcRxZABpua;
vectalcRxZBCpua(1,i)=alcRxZBCpua;
vectalcRxZCApua(1,i)=alcRxZCApua;
vectalcRxZAGpua(1,i)=alcRxZAGpua;
vectalcRxZBGpua(1,i)=alcRxZBGpua;
vectalcRxZCGpua(1,i)=alcRxZCGpua;
vectalcRyZABpua(1,i)=alcRyZABpua;
vectalcRyZBCpua(1,i)=alcRyZBCpua;
vectalcRyZCApua(1,i)=alcRyZCApua;
vectalcRyZAGpua(1,i)=alcRyZAGpua;
vectalcRyZBGpua(1,i)=alcRyZBGpua;
vectalcRyZCGpua(1,i)=alcRyZCGpua;
vectalcYrZABpua(1,i)=alcYrZABpua;
vectalcYrZBCpua(1,i)=alcYrZBCpua;
vectalcYrZCApua(1,i)=alcYrZCApua;
vectalcYrZAGpua(1,i)=alcYrZAGpua;
vectalcYrZBGpua(1,i)=alcYrZBGpua;
vectalcYrZCGpua(1,i)=alcYrZCGpua;
vectalcYsZABpua(1,i)=alcYsZABpua;
vectalcYsZBCpua(1,i)=alcYsZBCpua;
vectalcYsZCApua(1,i)=alcYsZCApua;
vectalcYsZAGpua(1,i)=alcYsZAGpua;
vectalcYsZBGpua(1,i)=alcYsZBGpua;
vectalcYsZCGpua(1,i)=alcYsZCGpua;
%***************************************
%Para salida 3, 67N
%***************************************
Abscisa(1,i)=m;
vectorVspol(1,i)=Vspol;
vectorIsxop(1,i)=Isxop;
vectorAngDirSx(1,i)=AngDirSx;
vectorVspol(1,i)=Vspol;
vectorIsyop(1,i)=Isyop;
vectorAngDirSy(1,i)=AngDirSy;
vectorVxpol(1,i)=Vxpol;
vectorIxsop(1,i)=Ixsop;
vectorAngDirXs(1,i)=AngDirXs;
vectorVxpol(1,i)=Vxpol;
vectorIxrop(1,i)=Ixrop;
vectorAngDirXr(1,i)=AngDirXr;
vectorVrpol(1,i)=Vrpol;
vectorIrxop(1,i)=Irxop;
vectorAngDirRx(1,i)=AngDirRx;
vectorVrpol(1,i)=Vrpol;
vectorIryop(1,i)=Iryop;
vectorAngDirRy(1,i)=AngDirRy;
vectorVypol(1,i)=Vypol;
vectorIyrop(1,i)=Iyrop;
vectorAngDirYr(1,i)=AngDirYr;
vectorVypol(1,i)=Vypol;
vectorIysop(1,i)=Iysop;
vectorAngDirYs(1,i)=AngDirYs;
%***************************************
% Para Salida 4, R-X
%***************************************
vectorSx_ZAB(1,i)=Sx_ZAB;
vectorSx_ZBC(1,i)=Sx_ZBC;
vectorSx_ZCA(1,i)=Sx_ZCA;
vectorSx_ZAG(1,i)=Sx_ZAG;
vectorSx_ZBG(1,i)=Sx_ZBG;
vectorSx_ZCG(1,i)=Sx_ZCG;
vectorSy_ZAB(1,i)=Sy_ZAB;
vectorSy_ZBC(1,i)=Sy_ZBC;
vectorSy_ZCA(1,i)=Sy_ZCA;
vectorSy_ZAG(1,i)=Sy_ZAG;
vectorSy_ZBG(1,i)=Sy_ZBG;
vectorSy_ZCG(1,i)=Sy_ZCG;
vectorXs_ZAB(1,i)=Xs_ZAB;
vectorXs_ZBC(1,i)=Xs_ZBC;
vectorXs_ZCA(1,i)=Xs_ZCA;
vectorXs_ZAG(1,i)=Xs_ZAG;
vectorXs_ZBG(1,i)=Xs_ZBG;
vectorXs_ZCG(1,i)=Xs_ZCG;
vectorXr_ZAB(1,i)=Xr_ZAB;
vectorXr_ZBC(1,i)=Xr_ZBC;
vectorXr_ZCA(1,i)=Xr_ZCA;
vectorXr_ZAG(1,i)=Xr_ZAG;
vectorXr_ZBG(1,i)=Xr_ZBG;
vectorXr_ZCG(1,i)=Xr_ZCG;
vectorRx_ZAB(1,i)=Rx_ZAB;
vectorRx_ZBC(1,i)=Rx_ZBC;
vectorRx_ZCA(1,i)=Rx_ZCA;
vectorRx_ZAG(1,i)=Rx_ZAG;
vectorRx_ZBG(1,i)=Rx_ZBG;
vectorRx_ZCG(1,i)=Rx_ZCG;
vectorRy_ZAB(1,i)=Ry_ZAB;
vectorRy_ZBC(1,i)=Ry_ZBC;
vectorRy_ZCA(1,i)=Ry_ZCA;
vectorRy_ZAG(1,i)=Ry_ZAG;
vectorRy_ZBG(1,i)=Ry_ZBG;
vectorRy_ZCG(1,i)=Ry_ZCG;
vectorYr_ZAB(1,i)=Yr_ZAB;
vectorYr_ZBC(1,i)=Yr_ZBC;
vectorYr_ZCA(1,i)=Yr_ZCA;
vectorYr_ZAG(1,i)=Yr_ZAG;
vectorYr_ZBG(1,i)=Yr_ZBG;
vectorYr_ZCG(1,i)=Yr_ZCG;
vectorYs_ZAB(1,i)=Ys_ZAB;
vectorYs_ZBC(1,i)=Ys_ZBC;
vectorYs_ZCA(1,i)=Ys_ZCA;
vectorYs_ZAG(1,i)=Ys_ZAG;
vectorYs_ZBG(1,i)=Ys_ZBG;
vectorYs_ZCG(1,i)=Ys_ZCG;
%***************************************
i=i+1;
m=m+0.01;
end % end del while m
% --------------------------------------------------------
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
% --------------------------------------------------------
%***************************************
% Titulo para las Salidas 1,2,3
%***************************************
cadenam1=sprintf('PG54ic,m=p,%s,%3.0fkV,kE=%-4.2f,L1=%3.0fkm,L2=%3.0fkm,L3=%3.0fkm,L4=%3.0fkm,L5=%3.0fkm',...
comp,Un,kE,L1,L2,L3,L4,L5);
cadenam2=sprintf('Sccs=%-5.0f,Sccx=%-5.0f,Sccr=%-5.0f\nSccy=%-5.0f',...
Sccs,Sccx,Sccr,Sccy);
cadenam3=sprintf('Kx0=%-6.2f,Kr0=%-6.2f,Ky0=%-6.2f,k0=%-6.2f',...
Kx0,Kr0,Ky0,k0);
cadenam4=sprintf('m: %s,%s,Rfm=%-3.0f',ppfaltam,cfaltam,Rfm);
if faltam==7|faltam==8|faltam==9|faltam==11
cadenam5=sprintf('Rgm=%-3.0f',Rgm);
cadenam4=sprintf('%s,%s',cadenam4,cadenam5);
end
cadenam=sprintf('%s,%s,%s,%s',cadenam1,cadenam2,cadenam3,cadenam4);
cadenap=sprintf('p: %s,%s,Rfp=%-3.0f',ppfaltap,cfaltap,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);
%***************************************
% Titulo para la Salida 2, gráficas 67N
%***************************************
dircadenam1=sprintf('PG54ic,m=p,67N,%3.0fkV,kE=%-4.2f,L1=%3.0fkm,L2=%3.0fkm,L3=%3.0fkm,L4=%3.0fkm,L5=%3.0fkm',...
Un,kE,L1,L2,L3,L4,L5);
dircadenam2=sprintf('Sccs=%-5.0f,Sccx=%-5.0f,Sccr=%-5.0f\nSccy=%-5.0f',...
Sccs,Sccx,Sccr,Sccy);
dircadenam3=sprintf('Kx0=%-6.2f,Kr0=%-6.2f,Ky0=%-6.2f,k0=%-6.2f',...
Kx0,Kr0,Ky0,k0);
dircadenam4=sprintf('m:%s,%s, Rfm=%-3.0f',ppfaltam,cfaltam,Rfm);
if faltam==7|faltam==8|faltam==9|faltam==11
dircadenam5=sprintf('Rgm=%-3.0f',Rgm);
dircadenam4=sprintf('%s,%s',dircadenam4,dircadenam5);
end
dircadenam=sprintf('%s,%s,%s,%s',dircadenam1,dircadenam2,dircadenam3,dircadenam4);
dircadenap=sprintf('p: %s,%s,Rfp=%-3.0f',ppfaltap,cfaltap,Rfp);
if faltap==7|faltap==8|faltap==9|faltap==11
dircadenap1=sprintf('Rgp=%-3.0f',Rgp);
dircadenap=sprintf('%s,%s',dircadenap,dircadenap1);
end
dircadenamp=sprintf('%s, %s',dircadenam,dircadenap);
%*************************************
% ------------------------------------------------------
% SALIDA
% ------------------------------------------------------
variassalidas=1;
while variassalidas==1
fprintf('\n');
fprintf('Salida Grafica 1: alcance pu (Extremos Faltas m,p)\n');
fprintf('Salida Grafica 2: alcance pua, Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys\n');
fprintf('Salida Grafica 3: 67N, Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys\n');
fprintf('Salida Grafica 4: R-X, Sx,Sy,Xs,Xr,Rx,Ry,Yr,Ys\n');
fprintf('\n');
salida=input('Salida: 1,2,3,4: ');
switch(salida)
case{1}, % ********** salida 1, alcance pu **********
switch(tramo)
case {1}, % Sal1, tramo 1
extr=input('Extremo: Sx=1, Sy=2, Xs=3, Ys=4: ');
switch(extr)
case {1}, % Sal1, tramo 1, Extremo Sx
ap=0;
bp=L1/(L1+L4);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimSx(cont)=(0.8*(L1+L4))/(mm*L1);
Ordenadacp(cont)=ap+bp*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimSx);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcSxZABpu,'r',Abscisa,vectalcSxZBCpu,'g',Abscisa,vectalcSxZCApu,'b',...
Abscisa,vectalcSxZAGpu,'m',Abscisa,vectalcSxZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcSxZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSxpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {2}, % Sal1, tramo 1, Extremo Sy
ap=0;
bp=L1/(L1+L2+L5);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimSy(cont)=(0.8*(L1+L2+L5))/(mm*L1);
Ordenadacp(cont)=ap+bp*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimSy);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcSyZABpu,'r',Abscisa,vectalcSyZBCpu,'g',Abscisa,vectalcSyZCApu,'b',...
Abscisa,vectalcSyZAGpu,'m',Abscisa,vectalcSyZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcSyZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSypu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {3}, % Sal1, tramo 1, Extremo Xs
ap=1;
bp=-L1/(L4+L1);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimXs(cont)=(0.8*(L1+L4))/(L4+(1-mm)*L1);
Ordenadacp(cont)=ap+bp*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimXs);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcXsZABpu,'r',Abscisa,vectalcXsZBCpu,'g',Abscisa,vectalcXsZCApu,'b',...
Abscisa,vectalcXsZAGpu,'m',Abscisa,vectalcXsZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcXsZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceXspu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {4}, % Sal1, tramo 1, Extremo Ys
at=1;
bt=-L1/(L5+L2+L1);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimYs(cont)=(0.8*(L1+L2+L5))/(L5+L2+(1-mm)*L1);
Ordenadact(cont)=at+bt*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadact);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimYs);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcYsZABpu,'r',Abscisa,vectalcYsZBCpu,'g',Abscisa,vectalcYsZCApu,'b',...
Abscisa,vectalcYsZAGpu,'m',Abscisa,vectalcYsZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcYsZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceYspu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
end % Sal1, end del switch extr del caso 1
case {2}, % Sal1, tramo 2
extr=input('Extremo: Xr=1, Rx=2, Sy=3, Ys=4: ');
switch(extr)
case {1}, % Sal1, tramo 2, Extremo Xr
as=L4/(L4+L2+L3);
bs=L2/(L4+L2+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimXr(cont)=(0.8*(L4+L2+L3))/(L4+mm*L2);
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimXr);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcXrZABpu,'r',Abscisa,vectalcXrZBCpu,'g',Abscisa,vectalcXrZCApu,'b',...
Abscisa,vectalcXrZAGpu,'m',Abscisa,vectalcXrZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcXrZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceXrpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {2}, % Sal1, tramo 2, Extremo Rx
as=(L3+L2)/(L4+L2+L3);
bs=-L2/(L4+L2+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimRx(cont)=(0.8*(L4+L2+L3))/(L3+(1-mm)*L2);
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimRx);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcRxZABpu,'r',Abscisa,vectalcRxZBCpu,'g',Abscisa,vectalcRxZCApu,'b',...
Abscisa,vectalcRxZAGpu,'m',Abscisa,vectalcRxZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcRxZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceRxpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {3}, % Sal1, tramo 2, Extremo Sy
as=L1/(L1+L2+L5);
bs=L2/(L1+L2+L5);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimSy(cont)=(0.8*(L1+L2+L5))/(L1+mm*L2);
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimSy);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcSyZABpu,'r',Abscisa,vectalcSyZBCpu,'g',Abscisa,vectalcSyZCApu,'b',...
Abscisa,vectalcSyZAGpu,'m',Abscisa,vectalcSyZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcSyZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSypu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {4}, % Sal1, tramo 2, Extremo Ys
as=(L5+L2)/(L1+L2+L5);
bs=-L2/(L1+L2+L5);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimYs(cont)=(0.8*(L1+L2+L5))/(L5+(1-mm)*L2);
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimYs);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcYsZABpu,'r',Abscisa,vectalcYsZBCpu,'g',Abscisa,vectalcYsZCApu,'b',...
Abscisa,vectalcYsZAGpu,'m',Abscisa,vectalcYsZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcYsZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceYspu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
end % end del switch extr
case {3}, % Sal1, tramo 3
extr=input('Extremo: Xr=1, Rx=2, Yr=3, Ry=4: ');
switch(extr)
case {1}, % Sal1, tramo 3, Extremo Xr
at=(L4+L2)/(L4+L2+L3);
bs=L3/(L4+L2+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimXr(cont)=(0.8*(L4+L2+L3))/(L4+L2+mm*L3);
Ordenadact(cont)=at+bt*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimXr);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcXrZABpu,'r',Abscisa,vectalcXrZBCpu,'g',Abscisa,vectalcXrZCApu,'b',...
Abscisa,vectalcXrZAGpu,'m',Abscisa,vectalcXrZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcXrZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceXrpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {2}, % Sal1, tramo 3, Extremo Rx
ap=L3/(L4+L2+L3);
bp=-L3/(L4+L2+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimRx(cont)=(0.8*(L4+L2+L3))/((1-mm)*L3);
Ordenadacp(cont)=ap+bp*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimRx);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcRxZABpu,'r',Abscisa,vectalcRxZBCpu,'g',Abscisa,vectalcRxZCApu,'b',...
Abscisa,vectalcRxZAGpu,'m',Abscisa,vectalcRxZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcRxZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceRxpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {3}, % Sal1, tramo 3, Extremo Yr
as=L5/(L5+L3);
bs=L3/(L5+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimYr(cont)=(0.8*(L5+L3))/(L5+mm*L3);
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimYr);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcYrZABpu,'r',Abscisa,vectalcYrZBCpu,'g',Abscisa,vectalcYrZCApu,'b',...
Abscisa,vectalcYrZAGpu,'m',Abscisa,vectalcYrZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcYrZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceYrpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {4}, % Sal1, tramo 3, Extremo Ry
ap=L3/(L5+L3);
bp=-L3/(L5+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimRy(cont)=(0.8*(L5+L3))/((1-mm)*L3);
Ordenadacp(cont)=ap+bp*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimRy);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcRyZABpu,'r',Abscisa,vectalcRyZBCpu,'g',Abscisa,vectalcRyZCApu,'b',...
Abscisa,vectalcRyZAGpu,'m',Abscisa,vectalcRyZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcRyZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceRypu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
end % end del switch extr
case {4}, % Sal1, tramo 4
extr=input('Extremo: Sx=1, Xs=2, Xr=3, Rx=4: ');
switch(extr)
case {1}, % Sal1, tramo 4, Extremo Sx
as=L1/(L1+L4);
bs=L4/(L1+L4);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimSx(cont)=(0.8*(L1+L4))/(L1+mm*L4);
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimSx);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcSxZABpu,'r',Abscisa,vectalcSxZBCpu,'g',Abscisa,vectalcSxZCApu,'b',...
Abscisa,vectalcSxZAGpu,'m',Abscisa,vectalcSxZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcSxZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSxpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {2}, % Sal1, tramo 4, Extremo Xs
ap=1;
bp=-L1/(L1+L4);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimXs(cont)=(0.8*(L1+L4))/((1-mm)*L4);
Ordenadacp(cont)=ap+bp*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimXs);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcXsZABpu,'r',Abscisa,vectalcXsZBCpu,'g',Abscisa,vectalcXsZCApu,'b',...
Abscisa,vectalcXsZAGpu,'m',Abscisa,vectalcXsZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcXsZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceXspu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {3}, % Sal1, tramo 4, Extremo Xr
ap=L4/(L4+L2+L3);
bp=-L4/(L4+L2+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimXr(cont)=(0.8*(L4+L2+L3))/((1-mm)*L4);
Ordenadacp(cont)=ap+bp*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimXr);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcXrZABpu,'r',Abscisa,vectalcXrZBCpu,'g',Abscisa,vectalcXrZCApu,'b',...
Abscisa,vectalcXrZAGpu,'m',Abscisa,vectalcXrZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcXrZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceXrpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {4}, % Sal1, tramo 4, Extremo Rx
at=(L3+L2)/(L4+L2+L3);
bt=L4/(L4+L2+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimRx(cont)=(0.8*(L4+L2+L3))/(L3+L2+mm*L4);
Ordenadact(cont)=at+bt*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadact);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimRx);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcRxZABpu,'r',Abscisa,vectalcRxZBCpu,'g',Abscisa,vectalcRxZCApu,'b',...
Abscisa,vectalcRxZAGpu,'m',Abscisa,vectalcRxZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcRxZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceRxpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
end % end del switch extr
case {5}, % Sal1, tramo 5
extr=input('Extremo: Sy=1, Ys=2, Ry=3, Yr=4: ');
switch(extr)
case {1}, % Sal1, tramo 5, Extremo Sy
at=(L1+L2)/(L1+L2+L5);
bt=L5/(L1+L2+L5);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimSy(cont)=(0.8*(L1+L2+L5))/(L1+L2+mm*L5);
Ordenadact(cont)=at+bt*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadact);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimSy);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcSyZABpu,'r',Abscisa,vectalcSyZBCpu,'g',Abscisa,vectalcSyZCApu,'b',...
Abscisa,vectalcSyZAGpu,'m',Abscisa,vectalcSyZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcSyZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSypu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {2}, % Sal1, tramo 5, Extremo Ys
ap=L5/(L1+L2+L5);
bp=-L5/(L1+L2+L5);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimYs(cont)=(0.8*(L1+L2+L5))/((1-mm)*L5);
Ordenadacp(cont)=ap+bp*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimYs);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcYsZABpu,'r',Abscisa,vectalcYsZBCpu,'g',Abscisa,vectalcYsZCApu,'b',...
Abscisa,vectalcYsZAGpu,'m',Abscisa,vectalcYsZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcYsZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceYspu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {3}, % Sal1, tramo 5, Extremo Ry
as=L3/(L3+L5);
bs=L5/(L3+L5);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimRy(cont)=(0.8*(L3+L5))/(L3+mm*L5);
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimRy);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcRyZABpu,'r',Abscisa,vectalcRyZBCpu,'g',Abscisa,vectalcRyZCApu,'b',...
Abscisa,vectalcRyZAGpu,'m',Abscisa,vectalcRyZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcRyZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceRypu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
case {4}, % Sal1, tramo 5, Extremo Yr
ap=L5/(L3+L5);
bp=-L5/(L3+L5);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
LimYr(cont)=(0.8*(L3+L5))/((1-mm)*L5);
Ordenadacp(cont)=ap+bp*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color',[1 0.8431 0])
% Fin Colores
hold off
hold on
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
h=plot(Abscisa0,LimYr);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
hold on
plot(Abscisa,vectalcYrZABpu,'r',Abscisa,vectalcYrZBCpu,'g',Abscisa,vectalcYrZCApu,'b',...
Abscisa,vectalcYrZAGpu,'m',Abscisa,vectalcYrZBGpu,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcYrZCGpu);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceYrpu')
legend('Cob','Ref','Lim','AB','BC','CA','AG','BG','CG')
end % end del switch extr
end % end del switch tramo
case{2}, % ********** salida 2, alcance pua **************
extr=input('Extremo: Sx=1,Sy=2,Xs=3,Xr=4,Rx=5,Ry=6,Yr=7,Ys=8: ');
switch(extr)
case{1}, % Sx
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
cont=cont+1;
mm=mm+0.01;
end
ap=0;
bp=L1/(L1+L4);
as=L1/(L1+L4);
bs=L4/(L1+L4);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
if pfaltam==11
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color','c')
% Fin Colores
elseif pfaltam==12
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color','c')
% Fin Colores
end
hold off
hold on
plot(Abscisa,vectalcSxZABpua,'r',Abscisa,vectalcSxZBCpua,'g',Abscisa,vectalcSxZCApua,'b',...
Abscisa,vectalcSxZAGpua,'m',Abscisa,vectalcSxZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcSxZCGpua);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcSxpua')
if (pfaltam==11) | (pfaltam==12)
legend('Cob','AB','BC','CA','AG','BG','CG','z1')
else
legend('AB','BC','CA','AG','BG','CG','z1')
end
case{2}, % Sy
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
%Ordenadacp(cont)=ap+bp*mm;
%Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ap=0;
bp=L1/(L1+L2+L5);
as=L1/(L1+L2+L5);
bs=L2/(L1+L2+L5);
at=(L1+L2)/(L1+L2+L5);
bt=L5/(L1+L2+L5);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
Ordenadact(cont)=at+bt*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
if pfaltap==20
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color','c')
% Fin Colores
elseif pfaltap==19
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color','c')
% Fin Colores
elseif pfaltam==18
% Colores
b=plot(Abscisa0,Ordenadact);
set(b,'Color','c')
% Fin Colores
end
hold off
hold on
plot(Abscisa,vectalcSyZABpua,'r',Abscisa,vectalcSyZBCpua,'g',Abscisa,vectalcSyZCApua,'b',...
Abscisa,vectalcSyZAGpua,'m',Abscisa,vectalcSyZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcSyZCGpua);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcSypua')
if (pfaltap==20) | (pfaltap==19) | (pfaltam==18)
legend('Cob','AB','BC','CA','AG','BG','CG','z1')
else
legend('AB','BC','CA','AG','BG','CG','z1')
end
case{3}, % Xs
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
%Ordenadacp(cont)=ap+bp*mm;
%Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ap=1;
bp=-L1/(L4+L1);
as=L4/(L4+L1);
bs=-L4/(L4+L1);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
if pfaltam==11
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color','c')
% Fin Colores
elseif pfaltam==12
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color','c')
% Fin Colores
end
hold off
hold on
plot(Abscisa,vectalcXsZABpua,'r',Abscisa,vectalcXsZBCpua,'g',Abscisa,vectalcXsZCApua,'b',...
Abscisa,vectalcXsZAGpua,'m',Abscisa,vectalcXsZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcXsZCGpua);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcXspua')
if (pfaltam==11) | (pfaltam==12)
legend('Cob','AB','BC','CA','AG','BG','CG','z1')
else
legend('AB','BC','CA','AG','BG','CG','z1')
end
case{4}, % Xr
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
%Ordenadacp(cont)=ap+bp*mm;
%Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ap=L4/(L4+L2+L3);
bp=-L4/(L4+L2+L3);
as=L4/(L4+L2+L3);
bs=L2/(L4+L2+L3);
at=(L4+L2)/(L4+L2+L3);
bt=L3/(L4+L2+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
Ordenadact(cont)=at+bt*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
if pfaltap==13
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color','c')
% Fin Colores
elseif pfaltam==14
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color','c')
% Fin Colores
elseif pfaltam==15
% Colores
b=plot(Abscisa0,Ordenadact);
set(b,'Color','c')
% Fin Colores
end
hold off
hold on
plot(Abscisa,vectalcXrZABpua,'r',Abscisa,vectalcXrZBCpua,'g',Abscisa,vectalcXrZCApua,'b',...
Abscisa,vectalcXrZAGpua,'m',Abscisa,vectalcXrZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcXrZCGpua);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcXrpua')
if (pfaltap==13) | (pfaltam==14) | (pfaltam==15)
legend('Cob','AB','BC','CA','AG','BG','CG','z1')
else
legend('AB','BC','CA','AG','BG','CG','z1')
end
case{5}, % Rx
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
%Ordenadacp(cont)=ap+bp*mm;
%Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ap=L3/(L4+L2+L3);
bp=-L3/(L4+L2+L3);
as=(L3+L2)/(L3+L2+L4);
bs=-L2/(L3+L2+L4);
at=(L3+L2)/(L3+L2+L4);
bt=L4/(L3+L2+L4);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
Ordenadact(cont)=at+bt*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
if pfaltam==15
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color','c')
% Fin Colores
elseif pfaltam==14
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color','c')
% Fin Colores
elseif pfaltap==13
% Colores
b=plot(Abscisa0,Ordenadact);
set(b,'Color','c')
% Fin Colores
end
hold off
hold on
plot(Abscisa,vectalcRxZABpua,'r',Abscisa,vectalcRxZBCpua,'g',Abscisa,vectalcRxZCApua,'b',...
Abscisa,vectalcRxZAGpua,'m',Abscisa,vectalcRxZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcRxZCGpua);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcRxpua')
if (pfaltap==13) | (pfaltam==14) | (pfaltam==15)
legend('Cob','AB','BC','CA','AG','BG','CG','z1')
else
legend('AB','BC','CA','AG','BG','CG','z1')
end
case{6} % Ry
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
%Ordenadacp(cont)=ap+bp*mm;
%Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ap=L3/(L3+L5);
bp=-L3/(L3+L5);
as=L3/(L3+L5);
bs=L5/(L3+L5);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
if pfaltap==16
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color','c')
% Fin Colores
elseif pfaltap==17
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color','c')
% Fin Colores
end
hold off
hold on
plot(Abscisa,vectalcRyZABpua,'r',Abscisa,vectalcRyZBCpua,'g',Abscisa,vectalcRyZCApua,'b',...
Abscisa,vectalcRyZAGpua,'m',Abscisa,vectalcRyZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcRyZCGpua);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcRypua')
if (pfaltap==16) | (pfaltap==17)
legend('Cob','AB','BC','CA','AG','BG','CG','z1')
else
legend('AB','BC','CA','AG','BG','CG','z1')
end
case{7} % Yr
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
%Ordenadacp(cont)=ap+bp*mm;
%Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ap=L5/(L5+L3);
bp=-L5/(L5+L3);
as=L5/(L5+L3);
bs=L3/(L5+L3);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
if pfaltap==17
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color','c')
% Fin Colores
elseif pfaltap==16
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color','c')
% Fin Colores
end
hold off
hold on
plot(Abscisa,vectalcYrZABpua,'r',Abscisa,vectalcYrZBCpua,'g',Abscisa,vectalcYrZCApua,'b',...
Abscisa,vectalcYrZAGpua,'m',Abscisa,vectalcYrZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcYrZCGpua);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcYrpua')
if (pfaltap==16) | (pfaltap==17)
legend('Cob','AB','BC','CA','AG','BG','CG','z1')
else
legend('AB','BC','CA','AG','BG','CG','z1')
end
case{8} % Ys
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
%Ordenadacp(cont)=ap+bp*mm;
%Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
ap=L5/(L5+L2+L1);
bp=-L5/(L5+L2+L1);
as=(L5+L2)/(L5+L2+L1);
bs=-L2/(L5+L2+L1);
at=1;
bt=-L1/(L5+L2+L1);
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
Ordenadact(cont)=at+bt*mm;
cont=cont+1;
mm=mm+0.01;
end
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
hold on
if pfaltam==18
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color','c')
% Fin Colores
elseif pfaltap==19
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color','c')
% Fin Colores
elseif pfaltap==20
% Colores
b=plot(Abscisa0,Ordenadact);
set(b,'Color','c')
% Fin Colores
end
hold off
hold on
plot(Abscisa,vectalcYsZABpua,'r',Abscisa,vectalcYsZBCpua,'g',Abscisa,vectalcYsZCApua,'b',...
Abscisa,vectalcYsZAGpua,'m',Abscisa,vectalcYsZBGpua,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,vectalcYsZCGpua);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 ejemenosy ejemasy]);
title(cadenamp)
xlabel('m=p'),ylabel('alcYspua')
if (pfaltam==18) | (pfaltap==19) | (pfaltap==20)
legend('Cob','AB','BC','CA','AG','BG','CG','z1')
else
legend('AB','BC','CA','AG','BG','CG','z1')
end
end % end del switch extremo de salida 2
case{3}, % salida 3, 67N
extr=input('Extremo: Sx=1,Sy=2,Xs=3,Xr=4,Rx=5,Ry=6,Yr=7,Ys=8: ');
switch(extr)
case{1}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSx,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirSx(º)')
title(dircadenamp)
rayar(1);
legend('AngDirSx','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenam)
xlabel('m=p'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIsxop,'b');
grid
xlabel('m=p'),ylabel('modulo Isxop(kA)')
case{2}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSy,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirSy(º)')
title(dircadenamp)
rayar(1);
legend('AngDirSy','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenam)
xlabel('m=p'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIsyop,'b');
grid
xlabel('m=p'),ylabel('modulo Isyop(kA)')
case{3}
subplot(3,1,1)
plot(Abscisa,vectorAngDirXs,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirXs(º)')
title(dircadenamp)
rayar(1);
legend('AngDirXs','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVxpol,'b');
grid
%title(dircadenam)
xlabel('m=p'),ylabel('modulo Vxpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIxsop,'b');
grid
xlabel('m=p'),ylabel('modulo Ixsop(kA)')
case{4}
subplot(3,1,1)
plot(Abscisa,vectorAngDirXr,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirXr(º)')
title(dircadenamp)
rayar(1);
legend('AngDirXr','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVxpol,'b');
grid
%title(dircadenam)
xlabel('m=p'),ylabel('modulo Vxpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIxrop,'b');
grid
xlabel('m=p'),ylabel('modulo Ixrop(kA)')
case{5}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRx,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirRx(º)')
title(dircadenamp)
rayar(1);
legend('AngDirRx','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenam)
xlabel('m=p'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIrxop,'b');
grid
xlabel('m=p'),ylabel('modulo Irxop(kA)')
case{6}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRy,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirRy(º)')
title(dircadenamp)
rayar(1);
legend('AngDirRy','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenam)
xlabel('m=p'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIryop,'b');
grid
xlabel('m=p'),ylabel('modulo Iryop(kA)')
case{7}
subplot(3,1,1)
plot(Abscisa,vectorAngDirYr,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirYr(º)')
title(dircadenamp)
rayar(1);
legend('AngDirYr','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVypol,'b');
grid
%title(dircadenam)
xlabel('m=p'),ylabel('modulo Vypol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIyrop,'b');
grid
xlabel('m=p'),ylabel('modulo Iyrop(kA)')
case{8}
subplot(3,1,1)
plot(Abscisa,vectorAngDirYs,'b');
axis([0 1 -180 180]);
xlabel('m=p'),ylabel('AngDirYs(º)')
title(dircadenamp)
rayar(1);
legend('AngDirYs','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVypol,'b');
grid
%title(dircadenam)
xlabel('m=p'),ylabel('modulo Vypol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIysop,'b');
grid
xlabel('m=p'),ylabel('modulo Iysop(kA)')
end % end del switch extr
case{4}, % salida 4, R-X
extr=input('Extremo: Sx=1,Sy=2,Xs=3,Xr=4,Rx=5,Ry=6,Yr=7,Ys=8: ');
switch(extr)
case{1} % Extremo Sx
extremo='Sx';
pzona=0.8;
L=L1+L4;
hold on
[eje_x,eje_menos_y,eje_mas_y]=plot19d(L,zl1u,pzona);
hold off
hold on
% Colores
a=plot(vectorSx_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorSx_ZBC,'m')
plot(vectorSx_ZCA,'c')
plot(vectorSx_ZAG,'r')
plot(vectorSx_ZBG,'k')
% Colores
c=plot(vectorSx_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
plot22Sx(L,zl1u,L1,L2,L3,L4,L5,pfaltam,pfaltap)
[subAbscisa,subvectorSx_ZAB]=subconjunto21(Abscisa,vectorSx_ZAB);
subAbscisa(:,1)=[];
subvectorSx_ZAB(:,1)=[];
hold on
b=plot(subvectorSx_ZAB,'+');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,subvectorSx_ZAB,'naranja',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSx_ZBC]=subconjunto21(Abscisa,vectorSx_ZBC);
subAbscisa(:,1)=[];
subvectorSx_ZBC(:,1)=[];
hold on
plot(subvectorSx_ZBC,'m+');
plot20(L,zl1u,subvectorSx_ZBC,'m',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSx_ZCA]=subconjunto21(Abscisa,vectorSx_ZCA);
subAbscisa(:,1)=[];
subvectorSx_ZCA(:,1)=[];
hold on
plot(subvectorSx_ZCA,'c+');
plot20(L,zl1u,subvectorSx_ZCA,'c',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSx_ZAG]=subconjunto21(Abscisa,vectorSx_ZAG);
subAbscisa(:,1)=[];
subvectorSx_ZAG(:,1)=[];
hold on
plot(subvectorSx_ZAG,'r+');
plot20(L,zl1u,subvectorSx_ZAG,'r',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSx_ZBG]=subconjunto21(Abscisa,vectorSx_ZBG);
subAbscisa(:,1)=[];
subvectorSx_ZBG(:,1)=[];
hold on
plot(subvectorSx_ZBG,'k+');
plot20(L,zl1u,subvectorSx_ZBG,'k',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSx_ZCG]=subconjunto21(Abscisa,vectorSx_ZCG);
subAbscisa(:,1)=[];
subvectorSx_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorSx_ZCG,'+');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,subvectorSx_ZCG,'marron',eje_x,eje_menos_y,eje_mas_y)
hold off
RXcadenap=sprintf('Extr: %s',extremo);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
title(RXcadenamp)
legend('Mho','Lin','AB','BC','CA','AG','BG','CG')
case{2} % Extremo Sy
extremo='Sy';
pzona=0.8;
L=L1+L2+L5;
hold on
[eje_x,eje_menos_y,eje_mas_y]=plot19d(L,zl1u,pzona);
hold off
hold on
% Colores
a=plot(vectorSy_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorSy_ZBC,'m')
plot(vectorSy_ZCA,'c')
plot(vectorSy_ZAG,'r')
plot(vectorSy_ZBG,'k')
% Colores
c=plot(vectorSy_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
plot22Sy(L,zl1u,L1,L2,L3,L4,L5,pfaltam,pfaltap)
[subAbscisa,subvectorSy_ZAB]=subconjunto21(Abscisa,vectorSy_ZAB);
subAbscisa(:,1)=[];
subvectorSy_ZAB(:,1)=[];
hold on
b=plot(subvectorSy_ZAB,'+');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,subvectorSy_ZAB,'naranja',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSy_ZBC]=subconjunto21(Abscisa,vectorSy_ZBC);
subAbscisa(:,1)=[];
subvectorSy_ZBC(:,1)=[];
hold on
plot(subvectorSy_ZBC,'m+');
plot20(L,zl1u,subvectorSy_ZBC,'m',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSy_ZCA]=subconjunto21(Abscisa,vectorSy_ZCA);
subAbscisa(:,1)=[];
subvectorSy_ZCA(:,1)=[];
hold on
plot(subvectorSy_ZCA,'c+');
plot20(L,zl1u,subvectorSy_ZCA,'c',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSy_ZAG]=subconjunto21(Abscisa,vectorSy_ZAG);
subAbscisa(:,1)=[];
subvectorSy_ZAG(:,1)=[];
hold on
plot(subvectorSy_ZAG,'r+');
plot20(L,zl1u,subvectorSy_ZAG,'r',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSy_ZBG]=subconjunto21(Abscisa,vectorSy_ZBG);
subAbscisa(:,1)=[];
subvectorSy_ZBG(:,1)=[];
hold on
plot(subvectorSy_ZBG,'k+');
plot20(L,zl1u,subvectorSy_ZBG,'k',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorSy_ZCG]=subconjunto21(Abscisa,vectorSy_ZCG);
subAbscisa(:,1)=[];
subvectorSy_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorSy_ZCG,'+');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,subvectorSy_ZCG,'marron',eje_x,eje_menos_y,eje_mas_y)
hold off
RXcadenap=sprintf('Extr: %s',extremo);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
title(RXcadenamp)
legend('Mho','Lin','AB','BC','CA','AG','BG','CG')
case{3} % Extremo Xs
extremo='Xs';
pzona=0.8;
L=L1+L4;
hold on
[eje_x,eje_menos_y,eje_mas_y]=plot19d(L,zl1u,pzona);
hold off
hold on
% Colores
a=plot(vectorXs_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorXs_ZBC,'m')
plot(vectorXs_ZCA,'c')
plot(vectorXs_ZAG,'r')
plot(vectorXs_ZBG,'k')
% Colores
c=plot(vectorXs_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
plot22Xs(L,zl1u,L1,L2,L3,L4,L5,pfaltam,pfaltap)
[subAbscisa,subvectorXs_ZAB]=subconjunto21(Abscisa,vectorXs_ZAB);
subAbscisa(:,1)=[];
subvectorXs_ZAB(:,1)=[];
hold on
b=plot(subvectorXs_ZAB,'+');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,subvectorXs_ZAB,'naranja',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXs_ZBC]=subconjunto21(Abscisa,vectorXs_ZBC);
subAbscisa(:,1)=[];
subvectorXs_ZBC(:,1)=[];
hold on
plot(subvectorXs_ZBC,'m+');
plot20(L,zl1u,subvectorXs_ZBC,'m',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXs_ZCA]=subconjunto21(Abscisa,vectorXs_ZCA);
subAbscisa(:,1)=[];
subvectorXs_ZCA(:,1)=[];
hold on
plot(subvectorXs_ZCA,'c+');
plot20(L,zl1u,subvectorXs_ZCA,'c',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXs_ZAG]=subconjunto21(Abscisa,vectorXs_ZAG);
subAbscisa(:,1)=[];
subvectorXs_ZAG(:,1)=[];
hold on
plot(subvectorXs_ZAG,'r+');
plot20(L,zl1u,subvectorXs_ZAG,'r',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXs_ZBG]=subconjunto21(Abscisa,vectorXs_ZBG);
subAbscisa(:,1)=[];
subvectorXs_ZBG(:,1)=[];
hold on
plot(subvectorXs_ZBG,'k+');
plot20(L,zl1u,subvectorXs_ZBG,'k',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXs_ZCG]=subconjunto21(Abscisa,vectorXs_ZCG);
subAbscisa(:,1)=[];
subvectorXs_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorXs_ZCG,'+');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,subvectorXs_ZCG,'marron',eje_x,eje_menos_y,eje_mas_y)
hold off
RXcadenap=sprintf('Extr: %s',extremo);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
title(RXcadenamp)
legend('Mho','Lin','AB','BC','CA','AG','BG','CG')
case{4}, %Extremo Xr
extremo='Xr';
pzona=0.8;
L=L4+L2+L3;
hold on
[eje_x,eje_menos_y,eje_mas_y]=plot19d(L,zl1u,pzona);
hold off
hold on
% Colores
a=plot(vectorXr_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorXr_ZBC,'m')
plot(vectorXr_ZCA,'c')
plot(vectorXr_ZAG,'r')
plot(vectorXr_ZBG,'k')
% Colores
c=plot(vectorXr_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
plot22Xr(L,zl1u,L1,L2,L3,L4,L5,pfaltam,pfaltap)
[subAbscisa,subvectorXr_ZAB]=subconjunto21(Abscisa,vectorXr_ZAB);
subAbscisa(:,1)=[];
subvectorXr_ZAB(:,1)=[];
hold on
b=plot(subvectorXr_ZAB,'+');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,subvectorXr_ZAB,'naranja',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXr_ZBC]=subconjunto21(Abscisa,vectorXr_ZBC);
subAbscisa(:,1)=[];
subvectorXr_ZBC(:,1)=[];
hold on
plot(subvectorXr_ZBC,'m+');
plot20(L,zl1u,subvectorXr_ZBC,'m',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXr_ZCA]=subconjunto21(Abscisa,vectorXr_ZCA);
subAbscisa(:,1)=[];
subvectorXr_ZCA(:,1)=[];
hold on
plot(subvectorXr_ZCA,'c+');
plot20(L,zl1u,subvectorXr_ZCA,'c',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXr_ZAG]=subconjunto21(Abscisa,vectorXr_ZAG);
subAbscisa(:,1)=[];
subvectorXr_ZAG(:,1)=[];
hold on
plot(subvectorXr_ZAG,'r+');
plot20(L,zl1u,subvectorXr_ZAG,'r',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXr_ZBG]=subconjunto21(Abscisa,vectorXr_ZBG);
subAbscisa(:,1)=[];
subvectorXr_ZBG(:,1)=[];
hold on
plot(subvectorXr_ZBG,'k+');
plot20(L,zl1u,subvectorXr_ZBG,'k',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorXr_ZCG]=subconjunto21(Abscisa,vectorXr_ZCG);
subAbscisa(:,1)=[];
subvectorXr_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorXr_ZCG,'+');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,subvectorXr_ZCG,'marron',eje_x,eje_menos_y,eje_mas_y)
hold off
RXcadenap=sprintf('Extr: %s',extremo);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
title(RXcadenamp)
legend('Mho','Lin','AB','BC','CA','AG','BG','CG')
case{5}, %Extremo Rx
extremo='Rx';
pzona=0.8;
L=L4+L2+L3;
hold on
[eje_x,eje_menos_y,eje_mas_y]=plot19d(L,zl1u,pzona);
hold off
hold on
% Colores
a=plot(vectorRx_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorRx_ZBC,'m')
plot(vectorRx_ZCA,'c')
plot(vectorRx_ZAG,'r')
plot(vectorRx_ZBG,'k')
% Colores
c=plot(vectorRx_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
plot22Rx(L,zl1u,L1,L2,L3,L4,L5,pfaltam,pfaltap)
[subAbscisa,subvectorRx_ZAB]=subconjunto21(Abscisa,vectorRx_ZAB);
subAbscisa(:,1)=[];
subvectorRx_ZAB(:,1)=[];
hold on
b=plot(subvectorRx_ZAB,'+');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,subvectorRx_ZAB,'naranja',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRx_ZBC]=subconjunto21(Abscisa,vectorRx_ZBC);
subAbscisa(:,1)=[];
subvectorRx_ZBC(:,1)=[];
hold on
plot(subvectorRx_ZBC,'m+');
plot20(L,zl1u,subvectorRx_ZBC,'m',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRx_ZCA]=subconjunto21(Abscisa,vectorRx_ZCA);
subAbscisa(:,1)=[];
subvectorRx_ZCA(:,1)=[];
hold on
plot(subvectorRx_ZCA,'c+');
plot20(L,zl1u,subvectorRx_ZCA,'c',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRx_ZAG]=subconjunto21(Abscisa,vectorRx_ZAG);
subAbscisa(:,1)=[];
subvectorRx_ZAG(:,1)=[];
hold on
plot(subvectorRx_ZAG,'r+');
plot20(L,zl1u,subvectorRx_ZAG,'r',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRx_ZBG]=subconjunto21(Abscisa,vectorRx_ZBG);
subAbscisa(:,1)=[];
subvectorRx_ZBG(:,1)=[];
hold on
plot(subvectorRx_ZBG,'k+');
plot20(L,zl1u,subvectorRx_ZBG,'k',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRx_ZCG]=subconjunto21(Abscisa,vectorRx_ZCG);
subAbscisa(:,1)=[];
subvectorRx_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorRx_ZCG,'+');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,subvectorRx_ZCG,'marron',eje_x,eje_menos_y,eje_mas_y)
hold off
RXcadenap=sprintf('Extr: %s',extremo);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
title(RXcadenamp)
legend('Mho','Lin','AB','BC','CA','AG','BG','CG')
case{6}, %Extremo Ry
extremo='Ry';
pzona=0.8;
L=L3+L5;
hold on
[eje_x,eje_menos_y,eje_mas_y]=plot19d(L,zl1u,pzona);
hold off
hold on
% Colores
a=plot(vectorRy_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorRy_ZBC,'m')
plot(vectorRy_ZCA,'c')
plot(vectorRy_ZAG,'r')
plot(vectorRy_ZBG,'k')
% Colores
c=plot(vectorRy_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
plot22Ry(L,zl1u,L1,L2,L3,L4,L5,pfaltam,pfaltap)
[subAbscisa,subvectorRy_ZAB]=subconjunto21(Abscisa,vectorRy_ZAB);
subAbscisa(:,1)=[];
subvectorRy_ZAB(:,1)=[];
hold on
b=plot(subvectorRy_ZAB,'+');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,subvectorRy_ZAB,'naranja',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRy_ZBC]=subconjunto21(Abscisa,vectorRy_ZBC);
subAbscisa(:,1)=[];
subvectorRy_ZBC(:,1)=[];
hold on
plot(subvectorRy_ZBC,'m+');
plot20(L,zl1u,subvectorRy_ZBC,'m',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRy_ZCA]=subconjunto21(Abscisa,vectorRy_ZCA);
subAbscisa(:,1)=[];
subvectorRy_ZCA(:,1)=[];
hold on
plot(subvectorRy_ZCA,'c+');
plot20(L,zl1u,subvectorRy_ZCA,'c',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRy_ZAG]=subconjunto21(Abscisa,vectorRy_ZAG);
subAbscisa(:,1)=[];
subvectorRy_ZAG(:,1)=[];
hold on
plot(subvectorRy_ZAG,'r+');
plot20(L,zl1u,subvectorRy_ZAG,'r',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRy_ZBG]=subconjunto21(Abscisa,vectorRy_ZBG);
subAbscisa(:,1)=[];
subvectorRy_ZBG(:,1)=[];
hold on
plot(subvectorRy_ZBG,'k+');
plot20(L,zl1u,subvectorRy_ZBG,'k',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorRy_ZCG]=subconjunto21(Abscisa,vectorRy_ZCG);
subAbscisa(:,1)=[];
subvectorRy_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorRy_ZCG,'+');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,subvectorRy_ZCG,'marron',eje_x,eje_menos_y,eje_mas_y)
hold off
RXcadenap=sprintf('Extr: %s',extremo);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
title(RXcadenamp)
legend('Mho','Lin','AB','BC','CA','AG','BG','CG')
case{7}, %Extremo Yr
extremo='Yr';
pzona=0.8;
L=L3+L5;
hold on
[eje_x,eje_menos_y,eje_mas_y]=plot19d(L,zl1u,pzona);
hold off
hold on
% Colores
a=plot(vectorYr_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorYr_ZBC,'m')
plot(vectorYr_ZCA,'c')
plot(vectorYr_ZAG,'r')
plot(vectorYr_ZBG,'k')
% Colores
c=plot(vectorYr_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
plot22Yr(L,zl1u,L1,L2,L3,L4,L5,pfaltam,pfaltap)
[subAbscisa,subvectorYr_ZAB]=subconjunto21(Abscisa,vectorYr_ZAB);
subAbscisa(:,1)=[];
subvectorYr_ZAB(:,1)=[];
hold on
b=plot(subvectorYr_ZAB,'+');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,subvectorYr_ZAB,'naranja',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYr_ZBC]=subconjunto21(Abscisa,vectorYr_ZBC);
subAbscisa(:,1)=[];
subvectorYr_ZBC(:,1)=[];
hold on
plot(subvectorYr_ZBC,'m+');
plot20(L,zl1u,subvectorYr_ZBC,'m',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYr_ZCA]=subconjunto21(Abscisa,vectorYr_ZCA);
subAbscisa(:,1)=[];
subvectorYr_ZCA(:,1)=[];
hold on
plot(subvectorYr_ZCA,'c+');
plot20(L,zl1u,subvectorYr_ZCA,'c',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYr_ZAG]=subconjunto21(Abscisa,vectorYr_ZAG);
subAbscisa(:,1)=[];
subvectorYr_ZAG(:,1)=[];
hold on
plot(subvectorYr_ZAG,'r+');
plot20(L,zl1u,subvectorYr_ZAG,'r',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYr_ZBG]=subconjunto21(Abscisa,vectorYr_ZBG);
subAbscisa(:,1)=[];
subvectorYr_ZBG(:,1)=[];
hold on
plot(subvectorYr_ZBG,'k+');
plot20(L,zl1u,subvectorYr_ZBG,'k',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYr_ZCG]=subconjunto21(Abscisa,vectorYr_ZCG);
subAbscisa(:,1)=[];
subvectorYr_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorYr_ZCG,'+');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,subvectorYr_ZCG,'marron',eje_x,eje_menos_y,eje_mas_y)
hold off
RXcadenap=sprintf('Extr: %s',extremo);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
title(RXcadenamp)
legend('Mho','Lin','AB','BC','CA','AG','BG','CG')
case{8}, %Extremo Ys
extremo='Ys';
pzona=0.8;
L=L5+L2+L1;
hold on
[eje_x,eje_menos_y,eje_mas_y]=plot19d(L,zl1u,pzona);
hold off
hold on
% Colores
a=plot(vectorYs_ZAB);
set(a,'Color',[1 0.60 0])
% Fin Colores
plot(vectorYs_ZBC,'m')
plot(vectorYs_ZCA,'c')
plot(vectorYs_ZAG,'r')
plot(vectorYs_ZBG,'k')
% Colores
c=plot(vectorYs_ZCG);
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
hold off
plot22Ys(L,zl1u,L1,L2,L3,L4,L5,pfaltam,pfaltap)
[subAbscisa,subvectorYs_ZAB]=subconjunto21(Abscisa,vectorYs_ZAB);
subAbscisa(:,1)=[];
subvectorYs_ZAB(:,1)=[];
hold on
b=plot(subvectorYs_ZAB,'+');
set(b,'Color',[1 0.60 0])
plot20(L,zl1u,subvectorYs_ZAB,'naranja',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYs_ZBC]=subconjunto21(Abscisa,vectorYs_ZBC);
subAbscisa(:,1)=[];
subvectorYs_ZBC(:,1)=[];
hold on
plot(subvectorYs_ZBC,'m+');
plot20(L,zl1u,subvectorYs_ZBC,'m',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYs_ZCA]=subconjunto21(Abscisa,vectorYs_ZCA);
subAbscisa(:,1)=[];
subvectorYs_ZCA(:,1)=[];
hold on
plot(subvectorYs_ZCA,'c+');
plot20(L,zl1u,subvectorYs_ZCA,'c',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYs_ZAG]=subconjunto21(Abscisa,vectorYs_ZAG);
subAbscisa(:,1)=[];
subvectorYs_ZAG(:,1)=[];
hold on
plot(subvectorYs_ZAG,'r+');
plot20(L,zl1u,subvectorYs_ZAG,'r',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYs_ZBG]=subconjunto21(Abscisa,vectorYs_ZBG);
subAbscisa(:,1)=[];
subvectorYs_ZBG(:,1)=[];
hold on
plot(subvectorYs_ZBG,'k+');
plot20(L,zl1u,subvectorYs_ZBG,'k',eje_x,eje_menos_y,eje_mas_y)
hold off
[subAbscisa,subvectorYs_ZCG]=subconjunto21(Abscisa,vectorYs_ZCG);
subAbscisa(:,1)=[];
subvectorYs_ZCG(:,1)=[];
hold on
% Colores
c=plot(subvectorYs_ZCG,'+');
set(c,'Color',[0.7216 0.4510 0.20])
% Fin Colores
plot20(L,zl1u,subvectorYs_ZCG,'marron',eje_x,eje_menos_y,eje_mas_y)
hold off
RXcadenap=sprintf('Extr: %s',extremo);
RXcadenamp=sprintf('%s,%s',cadenamp,RXcadenap);
title(RXcadenamp)
legend('Mho','Lin','AB','BC','CA','AG','BG','CG')
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.
