% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG52c\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 Simple. Para cada Tramo: m=0...1\n');
fprintf('Barra m a elegir entre a,b,c,d,e,f,g,h,q,r\n');
fprintf('Analisis: Extremos de la Linea en Falta \n');
fprintf('Salida Grafica 1: Extremos m, Todas las Unidades (pu)\n');
fprintf('Salida Grafica 2: Extremos m, Todas las Unidades (pua)\n');
fprintf('Salida Grafica 3: 67N, Todas las Posiciones\n');
fprintf('\n');
% ----------------------------------------------
%PG52c
datos=input('datos: standard=1,input=2: ');
if datos==1
Un=220;
kE=1;
L1=100;
L2=100;
L3=100;
L4=100;
L5=100;
Sccs=5000;
Sccx=5000;
Sccr=5000;
Sccy=5000;
Kx0=1;
Kr0=1;
Ky0=1;
compens21=2;
else
% ----------------------------------------------
Un=input('U(Fase-Fase) kV:400,220,132,66: ');
kE=input('kE=E/Enominal: ');
L1=input('Longitud Tramo 1 km: ');
L2=input('Longitud Tramo 2 km: ');
L3=input('Longitud Tramo 3 km: ');
L4=input('Longitud Tramo 4 km: ');
L5=input('Longitud Tramo 5 km: ');
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: ');
compens21=input('Compensacion 21: Si=1, No=2 :');
end
%----------------------------------------------------------------
k0=input('k0=zs0/zr0: ');
fprintf('Punto en Falta: T1sup=11,T4izq=12,T4der=13,T2sup=14,T3sup=15\n');
fprintf(' : T3inf=16,T5der=17,T5izq=18,T2inf=19,T1inf=20\n');
pfaltam=input('Punto en Falta: ');
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
%----------------------------------------------------------------
En=Un/sqrt(3);
U=kE*Un;
E=kE*En;
if Un==400
zl1u=0.0346+0.3208i; zl2u=zl1u; zl0u=0.1828+1.1482i; zlM0u=0.1482+0.7383i;
elseif Un==220
zl1u=0.0829+0.4135i; zl2u=zl1u; zl0u=0.2311+1.3085i; zlM0u=0.1482+0.8031i;
elseif Un==132
zl1u=0.1342+0.4042i; zl2u=zl1u; zl0u=0.2824+1.3718i; zlM0u=0.1482+0.8947i;
elseif Un==66
zl1u=0.2156+0.3952i; zl2u=zl1u; zl0u=0.3636+1.4278i; zlM0u=0.1482+0.9288i;
end
k=(zl0u-zl1u)/(3*zl1u);
zs1=(Un^2/Sccs)*i;
zs2=zs1;
zx1=(Un^2/Sccx)*i;
zx2=zx1;
zr1=(Un^2/Sccr)*i;
zr2=zr1;
zy1=(Un^2/Sccy)*i;
zy2=zy1;
zx0=Kx0*zx1;
zr0=Kr0*zr1;
zy0=Ky0*zy1;
zs0=k0*zr0;
Km=zlM0u/zl1u;
fprintf('\n');
% ---------------------------------------------------------------
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;
% ---------------------------------------------------------------
[Minc]=subMinc(1);
i=1;
m=0.01;
while m<1;
if pfaltam==11 a=m;elseif pfaltam==12 b=m;elseif pfaltam==13 c=m;elseif pfaltam==14 d=m;
elseif pfaltam==15 e=m;elseif pfaltam==16 f=m;elseif pfaltam==17 g=m;
elseif pfaltam==18 h=m;elseif pfaltam==19 q=m;elseif pfaltam==20 r=m;
end
[z0]=subz0(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r);
[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); Z0(11,11);
Z1=inv(Y1); Z1(11,11);
Z2=inv(Y2); Z2(11,11);
% CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p DE FALTA
[Im0,Im1,Im2,Vm0,Vm1,Vm2]=mag2m(Z0,Z1,Z2,E,faltam,Rfm,Rgm,pfaltam);
% 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;
% 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);
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN EL PUNTO m DE FALTA
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);
%
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS EXTREMOS DE LINEA
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);
[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);
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);
[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);
% 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;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,La);
L=L4+(1-m)*L1;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,La);
cadylabel1='alcSmpu';
cadylabel2='alcXmpu';
cadylabel1a='alcSmpua';
cadylabel2a='alcXmpua';
% 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;
% Salida pu
Lim1=0.8*(L1+L4)/(m*L1);
Lim2=0.8*(L1+L4)/(L4+(1-m)*L1);
case{12},
em1='Extremo S';
em2='Extremo X';
L=L1+m*L4;
La=L1+L4;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,La);
L=(1-m)*L4;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,La);
cadylabel1='alcSmpu';
cadylabel2='alcXmpu';
cadylabel1a='alcSmpua';
cadylabel2a='alcXmpua';
% 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;
% Salida pu
Lim1=0.8*(L1+L4)/(L1+m*L4);
Lim2=0.8*(L1+L4)/((1-m)*L4);
case{13},
em1='Extremo X';
em2='Extremo R';
L=(1-m)*L4;
La=L4+L2+L3;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,La);
L=L3+L2+m*L4;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,La);
cadylabel1='alcXmpu';
cadylabel2='alcRmpu';
cadylabel1a='alcXmpua';
cadylabel2a='alcRmpua';
% 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;
% Salida pu
Lim1=0.8*(L4+L2+L3)/((1-m)*L4);
Lim2=0.8*(L4+L2+L3)/(L3+L2+m*L4);
case{14},
em1='Extremo X';
em2='Extremo R';
L=L4+m*L2;
La=L4+L2+L3;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,La);
L=L3+(1-m)*L2;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,La);
cadylabel1='alcXmpu';
cadylabel2='alcRmpu';
cadylabel1a='alcXmpua';
cadylabel2a='alcRmpua';
% 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;
Lim1=0.8*(L4+L2+L3)/(L4+m*L2);
Lim2=0.8*(L4+L2+L3)/(L3+(1-m)*L2);
case{15},
em1='Extremo X';
em2='Extremo R';
L=L4+L2+m*L3;
La=L4+L2+L3;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,La);
L=(1-m)*L3;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,La);
cadylabel1='alcXmpu';
cadylabel2='alcRmpu';
cadylabel1a='alcXmpua';
cadylabel2a='alcRmpua';
% 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;
Lim1=0.8*(L4+L2+L3)/(L4+L2+m*L3);
Lim2=0.8*(L4+L2+L3)/((1-m)*L3);
case{16},
em1='Extremo R';
em2='Extremo Y';
L=(1-m)*L3;
La=L3+L5;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,La);
L=L5+m*L3;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,La);
cadylabel1='alcRmpu';
cadylabel2='alcYmpu';
cadylabel1a='alcRmpua';
cadylabel2a='alcYmpua';
% 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;
Lim1=0.8*(L3+L5)/((1-m)*L3);
Lim2=0.8*(L3+L5)/(L5+m*L3);
case{17},
em1='Extremo R';
em2='Extremo Y';
L=L3+m*L5;
La=L3+L5;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,La);
L=(1-m)*L5;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,La);
cadylabel1='alcRmpu';
cadylabel2='alcYmpu';
cadylabel1a='alcRmpua';
cadylabel2a='alcYmpua';
% 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;
Lim1=0.8*(L3+L5)/(L3+m*L5);
Lim2=0.8*(L3+L5)/((1-m)*L5);
case{18},
em1='Extremo Y';
em2='Extremo S';
L=(1-m)*L5;
La=L5+L2+L1;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,La);
L=L1+L2+m*L5;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,La);
cadylabel1='alcYmpu';
cadylabel2='alcSmpu';
cadylabel1a='alcYmpua';
cadylabel2a='alcSmpua';
% 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;
Lim1=0.8*(L1+L2+L5)/((1-m)*L5);
Lim2=0.8*(L1+L2+L5)/(L1+L2+m*L5);
case{19},
em1='Extremo Y';
em2='Extremo S';
L=L5+(1-m)*L2;
La=L5+L2+L1;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,La);
L=L1+m*L2;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,La);
cadylabel1='alcYmpu';
cadylabel2='alcSmpu';
cadylabel1a='alcYmpua';
cadylabel2a='alcSmpua';
% 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;
Lim1=0.8*(L1+L2+L5)/(L5+(1-m)*L2);
Lim2=0.8*(L1+L2+L5)/(L1+m*L2);
case{20},
em1='Extremo Y';
em2='Extremo S';
L=L5+L2+(1-m)*L1;
La=L5+L2+L1;
[alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=...
alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L);
[alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=...
alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,La);
L=m*L1;
[alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=...
alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L);
[alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=...
alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,La);
cadylabel1='alcYmpu';
cadylabel2='alcSmpu';
cadylabel1a='alcYmpua';
cadylabel2a='alcSmpua';
% 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;
Lim1=0.8*(L1+L2+L5)/(L5+L2+(1-m)*L1);
Lim2=0.8*(L1+L2+L5)/(m*L1);
end
Abscisa(i)=m;
OrdenadaSmAB(1,i)=alcSmABpu;
OrdenadaSmBC(1,i)=alcSmBCpu;
OrdenadaSmCA(1,i)=alcSmCApu;
OrdenadaSmAG(1,i)=alcSmAGpu;
OrdenadaSmBG(1,i)=alcSmBGpu;
OrdenadaSmCG(1,i)=alcSmCGpu;
OrdenadaXmAB(1,i)=alcXmABpu;
OrdenadaXmBC(1,i)=alcXmBCpu;
OrdenadaXmCA(1,i)=alcXmCApu;
OrdenadaXmAG(1,i)=alcXmAGpu;
OrdenadaXmBG(1,i)=alcXmBGpu;
OrdenadaXmCG(1,i)=alcXmCGpu;
OrdenadaSmABa(1,i)=alcSmABpua;
OrdenadaSmBCa(1,i)=alcSmBCpua;
OrdenadaSmCAa(1,i)=alcSmCApua;
OrdenadaSmAGa(1,i)=alcSmAGpua;
OrdenadaSmBGa(1,i)=alcSmBGpua;
OrdenadaSmCGa(1,i)=alcSmCGpua;
OrdenadaXmABa(1,i)=alcXmABpua;
OrdenadaXmBCa(1,i)=alcXmBCpua;
OrdenadaXmCAa(1,i)=alcXmCApua;
OrdenadaXmAGa(1,i)=alcXmAGpua;
OrdenadaXmBGa(1,i)=alcXmBGpua;
OrdenadaXmCGa(1,i)=alcXmCGpua;
%***************************************
%Para salida 2, 67N
%***************************************
vectorVspol(1,i)=Vspol;
vectorIsxop(1,i)=Isxop;
vectorAngDirSx(1,i)=AngDirSx;
vectorIsyop(1,i)=Isyop;
vectorAngDirSy(1,i)=AngDirSy;
vectorVxpol(1,i)=Vxpol;
vectorIxsop(1,i)=Ixsop;
vectorAngDirXs(1,i)=AngDirXs;
vectorIxrop(1,i)=Ixrop;
vectorAngDirXr(1,i)=AngDirXr;
vectorVrpol(1,i)=Vrpol;
vectorIrxop(1,i)=Irxop;
vectorAngDirRx(1,i)=AngDirRx;
vectorIryop(1,i)=Iryop;
vectorAngDirRy(1,i)=AngDirRy;
vectorVypol(1,i)=Vypol;
vectorIyrop(1,i)=Iyrop;
vectorAngDirYr(1,i)=AngDirYr;
vectorIysop(1,i)=Iysop;
vectorAngDirYs(1,i)=AngDirYs;
%*****************************
% Formar Vectores para LimIzq y LimDer
vectorLim1(1,i)=Lim1;
vectorLim2(1,i)=Lim2;
i=i+1;
m=m+0.01;
end %(while m)
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
% ----------------------------------------------------------------------------------
% SALIDA GRAFICA
% ----------------------------------------------------------------------------------
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
cadenam1=sprintf('PG52c,%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 \n Sccx=%-5.0f,Sccr=%-5.0f,Sccy=%-5.0f,',Sccs,Sccx,Sccr,Sccy);
cadenam3=sprintf('Kx0=%-4.1f,Kr0=%-4.1f,Ky0=%-4.1f,k0=%-4.1f,T=%3.0f,m:%s,Rfm=%-3.0f',...
Kx0,Kr0,Ky0,k0,pfaltam,cfaltam,Rfm);
if faltam==1|faltam==2|faltam==3|faltam==4|faltam==5|faltam==6|faltam==10
cadenam=sprintf('%s %s %s',cadenam1,cadenam2,cadenam3);
elseif faltam==7|faltam==8|faltam==9|faltam==11
cadenam4=sprintf('Rgm=%-3.0f',Rgm);
cadenam=sprintf('%s %s %s %s',cadenam1,cadenam2,cadenam3,cadenam4);
end
%***************************************
% Titulo para las gráficas 67N
%***************************************
dircadenam1=sprintf('PG52c,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 \n Sccx=%-5.0f,Sccr=%-5.0f,Sccy=%-5.0f,',Sccs,Sccx,Sccr,Sccy);
dircadenam3=sprintf('Kx0=%-4.1f,Kr0=%-4.1f,Ky0=%-4.1f,k0=%-4.1f,T=%3.0f,m:%s,Rfm=%-3.0f',...
Kx0,Kr0,Ky0,k0,pfaltam,cfaltam,Rfm);
if faltam==1|faltam==2|faltam==3|faltam==4|faltam==5|faltam==6|faltam==10
dircadenam=sprintf('%s %s %s',dircadenam1,dircadenam2,dircadenam3);
elseif faltam==7|faltam==8|faltam==9|faltam==11
dircadenam4=sprintf('Rgm=%-3.0f',Rgm);
dircadenam=sprintf('%s %s %s %s',dircadenam1,dircadenam2,dircadenam3,dircadenam4);
end
%*************************************
variassalidas=1;
while variassalidas==1
fprintf('Salida Grafica 1: Extremos m, Todas las Unidades (pu)\n');
fprintf('Salida Grafica 2: Extremos m, Todas las Unidades (pua)\n');
fprintf('Salida Grafica 3: 67N, Todas las Posiciones\n');
fprintf('\n');
salida=input('Salida: 1,2,3:');
switch(salida)
case{1}, % ********* Salida 1 ************
mm=0.01;
cont=1;
while mm<1
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
Ordenadacp(cont)=ap+bp*mm;
Ordenadacs(cont)=as+bs*mm;
cont=cont+1;
mm=mm+0.01;
end
fprintf('%s =1,%s =2\n',em1,em2);
extr=input('Extremo:');
%extr=input('Extremo prim.=1, seg.=2: ');
switch(extr)
case{1},
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
hold on
h=plot(Abscisa,vectorLim1);
set(h,'Color',[0.7216 0.4510 0.20])
h=plot(Abscisa0,Ordenadacp);
set(h,'Color',[1 0.8431 0])
hold off
% Fin Colores
hold on
plot(Abscisa,OrdenadaSmAB,'r',Abscisa,OrdenadaSmBC,'g',Abscisa,OrdenadaSmCA,'b',...
Abscisa,OrdenadaSmAG,'m',Abscisa,OrdenadaSmBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSmCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 0 30]);
title(cadenam)
xlabel('m'),ylabel(cadylabel1)
legend('Ref','Lim','Cob','AB','BC','CA','AG','BG','CG')
case{2},
% Colores
b=plot(Abscisa0,Ordenada0);
set(b,'Color','c')
hold on
h=plot(Abscisa,vectorLim2);
set(h,'Color',[0.7216 0.4510 0.20])
h=plot(Abscisa0,Ordenadacs);
set(h,'Color',[1 0.8431 0])
hold off
% Fin Colores
hold on
plot(Abscisa,OrdenadaXmAB,'r',Abscisa,OrdenadaXmBC,'g',Abscisa,OrdenadaXmCA,'b',...
Abscisa,OrdenadaXmAG,'m',Abscisa,OrdenadaXmBG,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaXmCG);
set(a,'Color',[1 0.5490 0])
hold off
% Fin Colores
axis([0 1 0 3]);
title(cadenam)
xlabel('m'),ylabel(cadylabel2)
legend('Ref','Lim','Cob','AB','BC','CA','AG','BG','CG')
end
case(2) % ********* Salida 2 ************
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
fprintf('%s =1,%s =2\n',em1,em2);
extr=input('Extremo:');
%extr=input('Extremo prim.=1, seg.=2: ');
switch(extr)
case{1},
% Colores
b=plot(Abscisa0,Ordenadacp);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaSmABa,'r',Abscisa,OrdenadaSmBCa,'g',Abscisa,OrdenadaSmCAa,'b',...
Abscisa,OrdenadaSmAGa,'m',Abscisa,OrdenadaSmBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaSmCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 -2 2]);
title(cadenam)
xlabel('m'),ylabel(cadylabel1a)
legend('Ref,Cob','AB','BC','CA','AG','BG','CG','z1')
case{2},
% Colores
b=plot(Abscisa0,Ordenadacs);
set(b,'Color','c')
% Fin Colores
hold on
plot(Abscisa,OrdenadaXmABa,'r',Abscisa,OrdenadaXmBCa,'g',Abscisa,OrdenadaXmCAa,'b',...
Abscisa,OrdenadaXmAGa,'m',Abscisa,OrdenadaXmBGa,'k'),grid on
hold off
% Colores
hold on
a=plot(Abscisa,OrdenadaXmCGa);
set(a,'Color',[1 0.5490 0])
h=plot(Abscisa0,Ordenada0);
set(h,'Color',[0.7216 0.4510 0.20])
hold off
% Fin Colores
axis([0 1 -2 2]);
title(cadenam)
xlabel('m'),ylabel(cadylabel2a)
legend('Ref,Cob','AB','BC','CA','AG','BG','CG','z1')
end % end del switch extremo
%************************************************
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'),ylabel('AngDirSx(º)')
title(dircadenam)
rayar(1);
legend('AngDirSx','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIsxop,'b');
grid
xlabel('m'),ylabel('modulo Isxop(kA)')
case{2}
subplot(3,1,1)
plot(Abscisa,vectorAngDirSy,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirSy(º)')
title(dircadenam)
rayar(1);
legend('AngDirSy','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVspol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vspol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIsyop,'b');
grid
xlabel('m'),ylabel('modulo Isyop(kA)')
case{3}
subplot(3,1,1)
plot(Abscisa,vectorAngDirXs,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirXs(º)')
title(dircadenam)
rayar(1);
legend('AngDirXs','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVxpol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vxpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIxsop,'b');
grid
xlabel('m'),ylabel('modulo Ixsop(kA)')
case{4}
subplot(3,1,1)
plot(Abscisa,vectorAngDirXr,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirXr(º)')
title(dircadenam)
rayar(1);
legend('AngDirXr','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVxpol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vxpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIxrop,'b');
grid
xlabel('m'),ylabel('modulo Ixrop(kA)')
case{5}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRx,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirRx(º)')
title(dircadenam)
rayar(1);
legend('AngDirRx','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIrxop,'b');
grid
xlabel('m'),ylabel('modulo Irxop(kA)')
case{6}
subplot(3,1,1)
plot(Abscisa,vectorAngDirRy,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirRy(º)')
title(dircadenam)
rayar(1);
legend('AngDirRy','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVrpol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vrpol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIryop,'b');
grid
xlabel('m'),ylabel('modulo Iryop(kA)')
case{7}
subplot(3,1,1)
plot(Abscisa,vectorAngDirYr,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirYr(º)')
title(dircadenam)
rayar(1);
legend('AngDirYr','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVypol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vypol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIyrop,'b');
grid
xlabel('m'),ylabel('modulo Iyrop(kA)')
case{8}
subplot(3,1,1)
plot(Abscisa,vectorAngDirYs,'b');
axis([0 1 -180 180]);
xlabel('m'),ylabel('AngDirYs(º)')
title(dircadenam)
rayar(1);
legend('AngDirYs','AngMaxPar')
subplot(3,1,2)
plot(Abscisa,vectorVypol,'b');
grid
%title(dircadenam)
xlabel('m'),ylabel('modulo Vypol(kV)')
subplot(3,1,3)
plot(Abscisa,vectorIysop,'b');
grid
xlabel('m'),ylabel('modulo Iysop(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.
