% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG128Mn\n');
fprintf('Esquema Linea Doble SR\n');
fprintf('Unominal AT kV: 220\n');
fprintf('Unominal BT kV: 66\n');
fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf('Eleccion Ctes. del Tf\n');
fprintf('Tf AT/BT/T: Neutro a Tierra en AT y BT\n');
fprintf('Eleccion: kE,LA=LB=L,Sccs,Sccr,Kr0,k0\n');
fprintf('*** Parametro k0=zs0/zr0=0,1-1-10 ***\n');
fprintf('Impedancia de Transferencia Si, No\n');
fprintf('Eleccion: (zsr1,Ksr0)\n');
fprintf('Compensacion 21, Si, No\n');
fprintf('Falta Intercircuito entre Lineas de Diferente Tension\n');
fprintf('m=p=0...1\n');
fprintf('Barra m=1, Barra p=2\n');
%fprintf('Salida Grafica 1: Analisis SAm, Elegir Unidad (pu)\n');
%fprintf(' : Analisis SBp, Elegir Unidad (pu)\n');
%fprintf('Salida Grafica 2: Analisis SAm, Elegir Unidad (pua)\n');
%fprintf(' : Analisis SBp, Elegir Unidad (pua)\n');
fprintf('\n');
% -----------------------------------------------------------------------
%PG128Mc
datos=input('datos: standard=1,input=2: ');
% -----------------------------------------------------------------------
if datos==1
UA=220;
UB=66;
kE=1;
L=66;
Sp=120; Ss=120; St=40; Up=220; Us=66; Ut=33;
tccps=12; tccpt=5; tccst=2;
Sccs=5000;
Sccr=5000;
Kr0=1;
imptransfer=2;
compens21=2;
% ------------------------------------------------------------------------
else
UA=input('UA(Fase-Fase) kV, poner 220: ');
UB=input('UB(Fase-Fase) kV, poner 66: ');
Ut=input('Ut(Fase-Fase) kV, poner 33: ');
kE=input('kE= ');
L=input('Longitud L=LA=LB km: ');
Sccs=input('Sccs MVA: ');
Sccr=input('Sccr MVA: ');
Kr0=input('Fuente R, zr0/zr1: ');
Sp=input('Sp MVA, poner 120: ');
Ss=input('Ss MVA, poner 120: ');
St=input('St MVA, poner 40: ');
tccps=input('tccps %, Ref a Sp, poner 12: ');
tccpt=input('tccpt %, Ref a St, poner 5: ');
tccst=input('tccst %, Ref a St, poner 2: ');
imptransfer=input('imptransfer:si=1,no=2:');
compens21=input('Compensacion 21: Si=1, No=2 :');
end
% ------------------------------------------------------------------------
%k0=input('k0=zs0/zr0 (Tipicos=0.1-1-10): ');
fprintf('Tipo de Falta: \n');
fprintf('Aa=1,Bc=2,Cb=3,Aag=4,Bcg=5,Cbg=6,Aab=7,Bac=8,Cbc=9\n');
fprintf('Aabg=10,Bacg=11,Cbcg=12,Aabc=13,Aabcg=14,ABCabc=15,ABCabcg=16\n');
falta=input('Elegir Tipo de Falta: ');
Rf=input('Resistencia de Falta. Rf(Ohm) = ');
if falta==5|falta==6|falta==10|falta==11|falta==12|falta==14|falta==16
Rg=input('Resistencia comun de Falta. Rg(Ohm) = ');
else
Rg=0;
end
% ------------------------------------------------------------------------
if UA==400
zlA1u=0.0346+0.3208i; zlA2u=zlA1u; zlA0u=0.1828+1.1482i; zlM0u=0.1482+0.7383i;
elseif UA==220
zlA1u=0.0829+0.4135i; zlA2u=zlA1u; zlA0u=0.2311+1.3085i; zlM0u=0.1482+0.8031i;
elseif UA==132
zlA1u=0.1342+0.4042i; zlA2u=zl1u; zlA0u=0.2824+1.3718i; zlM0u=0.1482+0.8947i;
elseif UA==66
zlA1u=0.2156+0.3952i; zlA2u=zl1u; zlA0u=0.3636+1.4278i; zlM0u=0.1482+0.9288i;
end
if UB==400
zlB1u=0.0346+0.3208i; zlB2u=zlB1u; zlB0u=0.1828+1.1482i;
elseif UB==220
zlB1u=0.0829+0.4135i; zlB2u=zlB1u; zlB0u=0.2311+1.3085i;
elseif UB==132
zlB1u=0.1342+0.4042i; zlB2u=zlB1u; zlB0u=0.2824+1.3718i;
elseif UB==66
zlB1u=0.2156+0.3952i; zlB2u=zlB1u; zlB0u=0.3636+1.4278i;
end
% ------------------------------------------------------------------------
LA=L;
LB=L;
%EAv=UA*1000/sqrt(3)+0i; EA=1.0;
%EBv=UB*1000/sqrt(3)+0i; EB=1.0;
E=1.0;
E=kE*E;
kA=(zlA0u-zlA1u)/(3*zlA1u);
kB=(zlB0u-zlB1u)/(3*zlB1u);
Sbase=100; % MVA
UbaseA=UA; % kV
UbaseB=UB; % kV
VbaseA=UbaseA/sqrt(3); % kV
VbaseB=UbaseB/sqrt(3); % kV
zbaseA=UA^2/Sbase; % Ohm
zbaseB=UB^2/Sbase; % Ohm
IbaseA=(UA/sqrt(3))/zbaseA; % kA
IbaseB=(UB/sqrt(3))/zbaseB; %kA
KmA=zlM0u/zlA1u; KmA=KmA*(IbaseB/IbaseA);
KmB=zlM0u/zlB1u; KmB=KmB*(IbaseA/IbaseB);
fprintf('\n');
if imptransfer==1
zsr1=input('zsr1=zsr2(Complejo r+xi en Ohms): ');
zsr1=zsr1/zbaseA;
zsr2=zsr1;
Ksr0=input('Ksr0=zsr0/zsr1= ');
zsr0=Ksr0*zsr1;
end
if falta==2 cfalta='Bc';elseif falta==3 cfalta='Cb';elseif falta==5 cfalta='Bcg';
elseif falta==6 cfalta='Cbg';elseif falta==7 cfalta='Aab';elseif falta==8 cfalta='Bac';
elseif falta==9 cfalta='Cbc';elseif falta==10 cfalta='Aabg';elseif falta==11 cfalta='Bacg';
elseif falta==12 cfalta='Cbcg';elseif falta==13 cfalta='Aabc';elseif falta==14 cfalta='Aabcg';
elseif falta==15 cfalta='ABCabc';elseif falta==16 cfalta='ABCabcg';
end
% ------------------------------------------------------------------------
fila=1;
k0=0.1;
while k0<11
i=1;
m=0.01;
while m<1;
p=m;
% REDES DE SECUENCIA 0,1,2 pu
[zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1]=pusec1(UA,Sccs,Sccr,zlA1u,m,LA,zlB1u,p,LB,zbaseA,zbaseB);
[zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2]=pusec2(UA,Sccs,Sccr,zlA2u,m,LA,zlB2u,p,LB,zbaseA,zbaseB);
[zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0]=...
pusec0(k0,Kr0,zr1,zlA0u,m,LA,zlB0u,p,LB,zlM0u,zbaseA,zbaseB,UbaseA,UbaseB,Sbase);
zpTS=1/2*((tccps/100)*UA^2/Sp+(tccpt/100)*UA^2/St-(tccst/100)*UA^2/St); zpTS=zpTS/zbaseA;
zsTS=1/2*((tccps/100)*UA^2/Sp-(tccpt/100)*UA^2/St+(tccst/100)*UA^2/St); zsTS=zsTS/zbaseA;
ztTS=1/2*(-(tccps/100)*UA^2/Sp+(tccpt/100)*UA^2/St+(tccst/100)*UA^2/St); ztTS=ztTS/zbaseA;
zpTR=zpTS; zsTR=zsTS; ztTR=ztTS;
% CALCULO DE LA REDUCCION DE REDES 0,1,2 A LOS PUNTOS m,p DE FALTA
switch(imptransfer)
case{1},
[Minc0]=subMinc0T(1);
[Minc1]=subMinc1T(1);
[Minc2]=subMinc2T(1);
[z1]=subz1T(zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1,zpTS,zsTS,zpTR,zsTR,zsr1);
[z2]=subz2T(zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2,zpTS,zsTS,zpTR,zsTR,zsr2);
[z0]=subz0T(zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0,zpTS,zsTS,ztTS,zpTR,zsTR,ztTR,zsr0);
case(2)
[Minc0]=subMinc0(1);
[Minc1]=subMinc1(1);
[Minc2]=subMinc2(1);
[z1]=subz1(zs1,zr1,zSAm1,zRAm1,zSBp1,zRBp1,zpTS,zsTS,zpTR,zsTR);
[z2]=subz2(zs2,zr2,zSAm2,zRAm2,zSBp2,zRBp2,zpTS,zsTS,zpTR,zsTR);
[z0]=subz0(zs0,zr0,zSAm0,zRAm0,zSBp0,zRBp0,zMSm0,zMRm0,zpTS,zsTS,ztTS,zpTR,zsTR,ztTR);
end
% MATRIZ DE ADMITANCIA DE BARRAS
[Y0]=subY0(Minc0,z0);
[Y1]=subY1(Minc1,z1);
[Y2]=subY2(Minc2,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]=...
maMDmp(Z0,Z1,Z2,E,falta,Rf,Rg,VbaseA,VbaseB,IbaseA,IbaseB);
% IBUS(0,1,2)
I0=[-Im0;-Ip0;0;0;0;0;0;0];
I1=[-Im1;-Ip1;0;0;0;0];
I2=[-Im2;-Ip2;0;0;0;0];
% VBUS(0,1,2)
V0=Z0*I0;
V1=Z1*I1+E;
V2=Z2*I2;
% I(0,1,2) DE LOS ELEMENTOS
[Ie0]=subIe0(Minc0,V0,z0);
[Ie1]=subIe1(Minc1,V1,z1);
[Ie2]=subIe2(Minc2,V2,z2);
VSA0=V0(3); VSA1=V1(3); VSA2=V2(3);
VRA0=V0(4); VRA1=V1(4); VRA2=V2(4);
VSB0=V0(5); VSB1=V1(5); VSB2=V2(5);
VRB0=V0(6); VRB1=V1(6); VRB2=V2(6);
ISAm0=Ie0(3); ISAm1=Ie1(3); ISAm2=Ie2(3);
IRAm0=Ie0(4); IRAm1=Ie1(4); IRAm2=Ie2(4);
ISBp0=Ie0(7); ISBp1=Ie1(7); ISBp2=Ie2(7);
IRBp0=Ie0(8); IRBp1=Ie1(8); IRBp2=Ie2(8);
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS PUNTOS m,p DE FALTA
[Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2);
[Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2);
[Ima,Imb,Imc]=abc(Im0,Im1,Im2);
[Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2);
% CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS EXTREMOS DE LINEA
[VSAa,VSAb,VSAc]=abc(VSA0,VSA1,VSA2);
[VRAa,VRAb,VRAc]=abc(VRA0,VRA1,VRA2);
[VSBa,VSBb,VSBc]=abc(VSB0,VSB1,VSB2);
[VRBa,VRBb,VRBc]=abc(VRB0,VRB1,VRB2);
[ISAma,ISAmb,ISAmc]=abc(ISAm0,ISAm1,ISAm2);
[ISBpa,ISBpb,ISBpc]=abc(ISBp0,ISBp1,ISBp2);
[IRAma,IRAmb,IRAmc]=abc(IRAm0,IRAm1,IRAm2);
[IRBpa,IRBpb,IRBpc]=abc(IRBp0,IRBp1,IRBp2);
%
% PROTECCION DE DISTANCIA (21)
switch(compens21)
case{1},
% COMPENSADA
[SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG]=...
PDistComp(VSAa,VSAb,VSAc,ISAma,ISAmb,ISAmc,ISAm0,kA,KmA,ISBp0);
[RAm_ZAB,RAm_ZBC,RAm_ZCA,RAm_ZAG,RAm_ZBG,RAm_ZCG]=...
PDistComp(VRAa,VRAb,VRAc,IRAma,IRAmb,IRAmc,IRAm0,kA,KmA,IRBp0);
[SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG]=...
PDistComp(VSBa,VSBb,VSBc,ISBpa,ISBpb,ISBpc,ISBp0,kB,KmB,ISAm0);
[RBp_ZAB,RBp_ZBC,RBp_ZCA,RBp_ZAG,RBp_ZBG,RBp_ZCG]=...
PDistComp(VRBa,VRBb,VRBc,IRBpa,IRBpb,IRBpc,IRBp0,kB,KmB,IRAm0);
case{2},
% NO COMPENSADA
[SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG]=...
PDist(VSAa,VSAb,VSAc,ISAma,ISAmb,ISAmc,ISAm0,kA);
[RAm_ZAB,RAm_ZBC,RAm_ZCA,RAm_ZAG,RAm_ZBG,RAm_ZCG]=...
PDist(VRAa,VRAb,VRAc,IRAma,IRAmb,IRAmc,IRAm0,kA);
[SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG]=...
PDist(VSBa,VSBb,VSBc,ISBpa,ISBpb,ISBpc,ISBp0,kB);
[RBp_ZAB,RBp_ZBC,RBp_ZCA,RBp_ZAG,RBp_ZBG,RBp_ZCG]=...
PDist(VRBa,VRBb,VRBc,IRBpa,IRBpb,IRBpc,IRBp0,kB);
end
% CALCULO DE LAS MAGNITUDES DE ALIMETACION DE LA PROTECCION
% DIRECCIONAL DE TIERRA EN LOS EXTREMOS DE LINEA
%
[VSApol,ISAmop,AngDirSAm]=PDire(VSA0,ISAm0); % pu
[VSBpol,ISBpop,AngDirSBp]=PDire(VSB0,ISBp0);
[VRApol,IRAmop,AngDirRAm]=PDire(VRA0,IRAm0);
[VRBpol,IRBpop,AngDirRBp]=PDire(VRB0,IRBp0);
%
% Conversion I(0,1,2)pu a I(0,1,2)kA
[IIm0,IIm1,IIm2]=conv3(Im0,Im1,Im2,IbaseA);
[IIp0,IIp1,IIp2]=conv3(Ip0,Ip1,Ip2,IbaseB);
[IISAm0,IISAm1,IISAm2,IIRAm0,IIRAm1,IIRAm2]=...
conv6(ISAm0,ISAm1,ISAm2,IRAm0,IRAm1,IRAm2,IbaseA);
[IISBp0,IISBp1,IISBp2,IIRBp0,IIRBp1,IIRBp2]=...
conv6(ISBp0,ISBp1,ISBp2,IRBp0,IRBp1,IRBp2,IbaseB);
% Conversion V(0,1,2)pu a V(0,1,2)kV
[VVm0,VVm1,VVm2,VVSA0,VVSA1,VVSA2,VVRA0,VVRA1,VVRA2]=...
conv9(Vm0,Vm1,Vm2,VSA0,VSA1,VSA2,VRA0,VRA1,VRA2,VbaseA);
[VVp0,VVp1,VVp2,VVSB0,VVSB1,VVSB2,VVRB0,VVRB1,VVRB2]=...
conv9(Vp0,Vp1,Vp2,VSB0,VSB1,VSB2,VRB0,VRB1,VRB2,VbaseB);
% Conversion I(a,b,c)pu a I(a,b,c)kA
[IIma,IImb,IImc]=conv3(Ima,Imb,Imc,IbaseA);
[IIpa,IIpb,IIpc]=conv3(Ipa,Ipb,Ipc,IbaseB);
[IISAma,IISAmb,IISAmc,IIRAma,IIRAmb,IIRAmc]=...
conv6(ISAma,ISAmb,ISAmc,IRAma,IRAmb,IRAmc,IbaseA);
[IISBpa,IISBpb,IISBpc,IIRBpa,IIRBpb,IIRBpc]=...
conv6(ISBpa,ISBpb,ISBpc,IRBpa,IRBpb,IRBpc,IbaseB);
% Conversion V(a,b,c)pu a V(a,b,c)Volt
[VVma,VVmb,VVmc]=conv3(Vma,Vmb,Vmc,VbaseA);
[VVpa,VVpb,VVpc]=conv3(Vpa,Vpb,Vpc,VbaseB);
[VVSAa,VVSAb,VVSAc,VVRAa,VVRAb,VVRAc]=...
conv6(VSAa,VSAb,VSAc,VRAa,VRAb,VRAc,VbaseA);
[VVSBa,VVSBb,VVSBc,VVRBa,VVRBb,VVRBc]=...
conv6(VSBa,VSBb,VSBc,VRBa,VRBb,VRBc,VbaseB);
% 21: Conversion Z(pu) a Z(Ohm)
[SAm_ZZAB,SAm_ZZBC,SAm_ZZCA,SAm_ZZAG,SAm_ZZBG,SAm_ZZCG]=...
conv6(SAm_ZAB,SAm_ZBC,SAm_ZCA,SAm_ZAG,SAm_ZBG,SAm_ZCG,zbaseA);
[RAm_ZZAB,RAm_ZZBC,RAm_ZZCA,RAm_ZZAG,RAm_ZZBG,RAm_ZZCG]=...
conv6(RAm_ZAB,RAm_ZBC,RAm_ZCA,RAm_ZAG,RAm_ZBG,RAm_ZCG,zbaseA);
[SBp_ZZAB,SBp_ZZBC,SBp_ZZCA,SBp_ZZAG,SBp_ZZBG,SBp_ZZCG]=...
conv6(SBp_ZAB,SBp_ZBC,SBp_ZCA,SBp_ZAG,SBp_ZBG,SBp_ZCG,zbaseB);
[RBp_ZZAB,RBp_ZZBC,RBp_ZZCA,RBp_ZZAG,RBp_ZZBG,RBp_ZZCG]=...
conv6(RBp_ZAB,RBp_ZBC,RBp_ZCA,RBp_ZAG,RBp_ZBG,RBp_ZCG,zbaseB);
% 67N: Conversion pu a kV
[VVSApol,VVRApol]=conv2(VSApol,VRApol,VbaseA);
[VVSBpol,VVRBpol]=conv2(VSBpol,VRBpol,VbaseB);
% 67N: Conversion pu a kA
[IISAmop,IIRAmop]=conv2(ISAmop,IRAmop,IbaseA); % kA
[IISBpop,IIRBpop]=conv2(ISBpop,IRBpop,IbaseB);
alcSAmAB=abs(SAm_ZZAB)/cos(angle(SAm_ZZAB)-angle(zlA1u));
alcSAmBC=abs(SAm_ZZBC)/cos(angle(SAm_ZZBC)-angle(zlA1u));
alcSAmCA=abs(SAm_ZZCA)/cos(angle(SAm_ZZCA)-angle(zlA1u));
alcSAmAG=abs(SAm_ZZAG)/cos(angle(SAm_ZZAG)-angle(zlA1u));
alcSAmBG=abs(SAm_ZZBG)/cos(angle(SAm_ZZBG)-angle(zlA1u));
alcSAmCG=abs(SAm_ZZCG)/cos(angle(SAm_ZZCG)-angle(zlA1u));
alcSBpAB=abs(SBp_ZZAB)/cos(angle(SBp_ZZAB)-angle(zlB1u));
alcSBpBC=abs(SBp_ZZBC)/cos(angle(SBp_ZZBC)-angle(zlB1u));
alcSBpCA=abs(SBp_ZZCA)/cos(angle(SBp_ZZCA)-angle(zlB1u));
alcSBpAG=abs(SBp_ZZAG)/cos(angle(SBp_ZZAG)-angle(zlB1u));
alcSBpBG=abs(SBp_ZZBG)/cos(angle(SBp_ZZBG)-angle(zlB1u));
alcSBpCG=abs(SBp_ZZCG)/cos(angle(SBp_ZZCG)-angle(zlB1u));
alcRAmAB=abs(RAm_ZZAB)/cos(angle(RAm_ZZAB)-angle(zlA1u));
alcRAmBC=abs(RAm_ZZBC)/cos(angle(RAm_ZZBC)-angle(zlA1u));
alcRAmCA=abs(RAm_ZZCA)/cos(angle(RAm_ZZCA)-angle(zlA1u));
alcRAmAG=abs(RAm_ZZAG)/cos(angle(RAm_ZZAG)-angle(zlA1u));
alcRAmBG=abs(RAm_ZZBG)/cos(angle(RAm_ZZBG)-angle(zlA1u));
alcRAmCG=abs(RAm_ZZCG)/cos(angle(RAm_ZZCG)-angle(zlA1u));
alcRBpAB=abs(RBp_ZZAB)/cos(angle(RBp_ZZAB)-angle(zlB1u));
alcRBpBC=abs(RBp_ZZBC)/cos(angle(RBp_ZZBC)-angle(zlB1u));
alcRBpCA=abs(RBp_ZZCA)/cos(angle(RBp_ZZCA)-angle(zlB1u));
alcRBpAG=abs(RBp_ZZAG)/cos(angle(RBp_ZZAG)-angle(zlB1u));
alcRBpBG=abs(RBp_ZZBG)/cos(angle(RBp_ZZBG)-angle(zlB1u));
alcRBpCG=abs(RBp_ZZCG)/cos(angle(RBp_ZZCG)-angle(zlB1u));
% para Salida 1
alcSAmABpu=alcSAmAB/abs(zlA1u*m*L);
alcSAmBCpu=alcSAmBC/abs(zlA1u*m*L);
alcSAmCApu=alcSAmCA/abs(zlA1u*m*L);
alcSAmAGpu=alcSAmAG/abs(zlA1u*m*L);
alcSAmBGpu=alcSAmBG/abs(zlA1u*m*L);
alcSAmCGpu=alcSAmCG/abs(zlA1u*m*L);
alcSBpABpu=alcSBpAB/abs(zlB1u*m*L);
alcSBpBCpu=alcSBpBC/abs(zlB1u*m*L);
alcSBpCApu=alcSBpCA/abs(zlB1u*m*L);
alcSBpAGpu=alcSBpAG/abs(zlB1u*m*L);
alcSBpBGpu=alcSBpBG/abs(zlB1u*m*L);
alcSBpCGpu=alcSBpCG/abs(zlB1u*m*L);
alcRAmABpu=alcRAmAB/abs(zlA1u*(1-m)*L);
alcRAmBCpu=alcRAmBC/abs(zlA1u*(1-m)*L);
alcRAmCApu=alcRAmCA/abs(zlA1u*(1-m)*L);
alcRAmAGpu=alcRAmAG/abs(zlA1u*(1-m)*L);
alcRAmBGpu=alcRAmBG/abs(zlA1u*(1-m)*L);
alcRAmCGpu=alcRAmCG/abs(zlA1u*(1-m)*L);
alcRBpABpu=alcRBpAB/abs(zlB1u*(1-m)*L);
alcRBpBCpu=alcRBpBC/abs(zlB1u*(1-m)*L);
alcRBpCApu=alcRBpCA/abs(zlB1u*(1-m)*L);
alcRBpAGpu=alcRBpAG/abs(zlB1u*(1-m)*L);
alcRBpBGpu=alcRBpBG/abs(zlB1u*(1-m)*L);
alcRBpCGpu=alcRBpCG/abs(zlB1u*(1-m)*L);
cadylabel1='alcSAmpu';
cadylabel2='alcSBppu';
cadylabel3='alcRAmpu';
cadylabel4='alcRBppu';
% para Salida 2
alcSAmABpua=alcSAmAB/abs(zlA1u*L);
alcSAmBCpua=alcSAmBC/abs(zlA1u*L);
alcSAmCApua=alcSAmCA/abs(zlA1u*L);
alcSAmAGpua=alcSAmAG/abs(zlA1u*L);
alcSAmBGpua=alcSAmBG/abs(zlA1u*L);
alcSAmCGpua=alcSAmCG/abs(zlA1u*L);
alcSBpABpua=alcSBpAB/abs(zlB1u*L);
alcSBpBCpua=alcSBpBC/abs(zlB1u*L);
alcSBpCApua=alcSBpCA/abs(zlB1u*L);
alcSBpAGpua=alcSBpAG/abs(zlB1u*L);
alcSBpBGpua=alcSBpBG/abs(zlB1u*L);
alcSBpCGpua=alcSBpCG/abs(zlB1u*L);
alcRAmABpua=alcRAmAB/abs(zlA1u*L);
alcRAmBCpua=alcRAmBC/abs(zlA1u*L);
alcRAmCApua=alcRAmCA/abs(zlA1u*L);
alcRAmAGpua=alcRAmAG/abs(zlA1u*L);
alcRAmBGpua=alcRAmBG/abs(zlA1u*L);
alcRAmCGpua=alcRAmCG/abs(zlA1u*L);
alcRBpABpua=alcRBpAB/abs(zlB1u*L);
alcRBpBCpua=alcRBpBC/abs(zlB1u*L);
alcRBpCApua=alcRBpCA/abs(zlB1u*L);
alcRBpAGpua=alcRBpAG/abs(zlB1u*L);
alcRBpBGpua=alcRBpBG/abs(zlB1u*L);
alcRBpCGpua=alcRBpCG/abs(zlB1u*L);
cadylabel1a='alcSAmpua';
cadylabel2a='alcSBppua';
cadylabel3a='alcRAmpua';
cadylabel4a='alcRBppua';
Abscisa(fila,i)=m;
OrdenadaSAmAB(fila,i)=alcSAmABpu;
OrdenadaSAmBC(fila,i)=alcSAmBCpu;
OrdenadaSAmCA(fila,i)=alcSAmCApu;
OrdenadaSAmAG(fila,i)=alcSAmAGpu;
OrdenadaSAmBG(fila,i)=alcSAmBGpu;
OrdenadaSAmCG(fila,i)=alcSAmCGpu;
OrdenadaSBpAB(fila,i)=alcSBpABpu;
OrdenadaSBpBC(fila,i)=alcSBpBCpu;
OrdenadaSBpCA(fila,i)=alcSBpCApu;
OrdenadaSBpAG(fila,i)=alcSBpAGpu;
OrdenadaSBpBG(fila,i)=alcSBpBGpu;
OrdenadaSBpCG(fila,i)=alcSBpCGpu;
OrdenadaRAmAB(fila,i)=alcRAmABpu;
OrdenadaRAmBC(fila,i)=alcRAmBCpu;
OrdenadaRAmCA(fila,i)=alcRAmCApu;
OrdenadaRAmAG(fila,i)=alcRAmAGpu;
OrdenadaRAmBG(fila,i)=alcRAmBGpu;
OrdenadaRAmCG(fila,i)=alcRAmCGpu;
OrdenadaRBpAB(fila,i)=alcRBpABpu;
OrdenadaRBpBC(fila,i)=alcRBpBCpu;
OrdenadaRBpCA(fila,i)=alcRBpCApu;
OrdenadaRBpAG(fila,i)=alcRBpAGpu;
OrdenadaRBpBG(fila,i)=alcRBpBGpu;
OrdenadaRBpCG(fila,i)=alcRBpCGpu;
OrdenadaSAmABa(fila,i)=alcSAmABpua;
OrdenadaSAmBCa(fila,i)=alcSAmBCpua;
OrdenadaSAmCAa(fila,i)=alcSAmCApua;
OrdenadaSAmAGa(fila,i)=alcSAmAGpua;
OrdenadaSAmBGa(fila,i)=alcSAmBGpua;
OrdenadaSAmCGa(fila,i)=alcSAmCGpua;
OrdenadaSBpABa(fila,i)=alcSBpABpua;
OrdenadaSBpBCa(fila,i)=alcSBpBCpua;
OrdenadaSBpCAa(fila,i)=alcSBpCApua;
OrdenadaSBpAGa(fila,i)=alcSBpAGpua;
OrdenadaSBpBGa(fila,i)=alcSBpBGpua;
OrdenadaSBpCGa(fila,i)=alcSBpCGpua;
OrdenadaRAmABa(fila,i)=alcRAmABpua;
OrdenadaRAmBCa(fila,i)=alcRAmBCpua;
OrdenadaRAmCAa(fila,i)=alcRAmCApua;
OrdenadaRAmAGa(fila,i)=alcRAmAGpua;
OrdenadaRAmBGa(fila,i)=alcRAmBGpua;
OrdenadaRAmCGa(fila,i)=alcRAmCGpua;
OrdenadaRBpABa(fila,i)=alcRBpABpua;
OrdenadaRBpBCa(fila,i)=alcRBpBCpua;
OrdenadaRBpCAa(fila,i)=alcRBpCApua;
OrdenadaRBpAGa(fila,i)=alcRBpAGpua;
OrdenadaRBpBGa(fila,i)=alcRBpBGpua;
OrdenadaRBpCGa(fila,i)=alcRBpCGpua;
i=i+1;
m=m+0.01;
end %(while m)
fila=fila+1;
k0=k0*10;
end % FIN DE BUCLES k0, m
%***************************************
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
%***************************************
% Titulo para las Salidas 1,2
%***************************************
cadenamc1=sprintf('PG128Mn,%s,%3.0fkV,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,Kr0=%-5.2f\n',...
comp,UA,UB,kE,L,Sccs,Sccr,Kr0);
if imptransfer==1
cadenamT=sprintf('zsr1=%-6.2fi,Ksr0=%-4.2f',imag(zsr1)*zbaseA,Ksr0);
cadenamc1=sprintf('%s,%s',cadenamc1,cadenamT);
end
cadenamc2=sprintf('%s,Rf=%-3.0f',cfalta,Rf);
if falta==1|falta==2|falta==3|falta==4|falta==7|falta==8|falta==9|falta==13|falta==15
cadenamc=sprintf('%s,%s',cadenamc1,cadenamc2);
elseif falta==5|falta==6|falta==10|falta==11|falta==12|falta==14|falta==16
cadenamc3=sprintf('Rg=%-3.0f',Rg);
cadenamc=sprintf('%s,%s,%s',cadenamc1,cadenamc2,cadenamc3);
end
% ----------------------------------------------------------------------------------
% SALIDA GRAFICA
% ----------------------------------------------------------------------------------
variassalidas=1;
while variassalidas==1
fprintf('\n');
fprintf('Salida Grafica 1: Analisis SAm, Elegir Unidad (pu)\n');
fprintf(' : SBp, " \n');
fprintf('Salida Grafica 2: Analisis SAm, Elegir Unidad (pua)\n');
fprintf(' : SBp, " \n');
fprintf('\n');
salida=input('salida: 1, 2 :');
switch(salida)
case{1},% **salida 1** : alc...pu
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=1;
Ordenadadi(cont)=0.8/mm;
Ordenadadd(cont)=0.8/(1-mm);
cont=cont+1;
mm=mm+0.01;
end
extremo=input('extremo: SAm=1, SBp=2, (RAm=3), (RBp=4):');
switch(extremo)
case{1},% salida 1 : alc...pu, **extremo Sm**
unidad=input('unidad: AB=1, BC=2, CA=3, AG=4, BG=5, CG=6:');
switch(unidad)
case{1},% salida 1 : alc...pu, extremo SAm, **unidad AB**
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmAB(1:1,1:(i-1));
Ordenada12=OrdenadaSAmAB(2:2,1:(i-1));
Ordenada13=OrdenadaSAmAB(3:3,1:(i-1));
cadenaylabel='alcanceSAmABpu';
case{2},% salida 1 : alc...pu, extremo SAm, unidad BC
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmBC(1:1,1:(i-1));
Ordenada12=OrdenadaSAmBC(2:2,1:(i-1));
Ordenada13=OrdenadaSAmBC(3:3,1:(i-1));
cadenaylabel='alcanceSAmBCpu';
case{3},% salida 1 : alc...pu, extremo SAm, unidad CA
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmCA(1:1,1:(i-1));
Ordenada12=OrdenadaSAmCA(2:2,1:(i-1));
Ordenada13=OrdenadaSAmCA(3:3,1:(i-1));
cadenaylabel='alcanceSAmCApu';
case{4},% salida 1 : alc...pu, extremo SAm, unidad AG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmAG(1:1,1:(i-1));
Ordenada12=OrdenadaSAmAG(2:2,1:(i-1));
Ordenada13=OrdenadaSAmAG(3:3,1:(i-1));
cadenaylabel='alcanceSAmAGpu';
case{5},% salida 1 : alc...pu, extremo SAm, unidad BG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmBG(1:1,1:(i-1));
Ordenada12=OrdenadaSAmBG(2:2,1:(i-1));
Ordenada13=OrdenadaSAmBG(3:3,1:(i-1));
cadenaylabel='alcanceSAmBGpu';
case{6},% salida 1 : alc...pu, extremo SAm, unidad CG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmCG(1:1,1:(i-1));
Ordenada12=OrdenadaSAmCG(2:2,1:(i-1));
Ordenada13=OrdenadaSAmCG(3:3,1:(i-1));
cadenaylabel='alcanceSAmCGpu';
end %end del switch(unidad)
case{2},% salida 1 : alc...pu, **extremo SBp**
unidad=input('unidad: AB=1, BC=2, CA=3, AG=4, BG=5, CG=6:');
switch(unidad)
case{1},% salida 1 : alc...pu, extremo SBp, **unidad AB**
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpAB(1:1,1:(i-1));
Ordenada12=OrdenadaSBpAB(2:2,1:(i-1));
Ordenada13=OrdenadaSBpAB(3:3,1:(i-1));
cadenaylabel='alcanceSBpABpu';
case{2},% salida 1 : alc...pu, extremo SBp, unidad BC
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpBC(1:1,1:(i-1));
Ordenada12=OrdenadaSBpBC(2:2,1:(i-1));
Ordenada13=OrdenadaSBpBC(3:3,1:(i-1));
cadenaylabel='alcanceSBpBCpu';
case{3},% salida 1 : alc...pu, extremo SBp, unidad CA
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpCA(1:1,1:(i-1));
Ordenada12=OrdenadaSBpCA(2:2,1:(i-1));
Ordenada13=OrdenadaSBpCA(3:3,1:(i-1));
cadenaylabel='alcanceSBpCApu';
case{4},% salida 1 : alc...pu, extremo SBp, unidad AG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpAG(1:1,1:(i-1));
Ordenada12=OrdenadaSBpAG(2:2,1:(i-1));
Ordenada13=OrdenadaSBpAG(3:3,1:(i-1));
cadenaylabel='alcanceSBpAGpu';
case{5},% salida 1 : alc...pu, extremo SBp, unidad BG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpBG(1:1,1:(i-1));
Ordenada12=OrdenadaSBpBG(2:2,1:(i-1));
Ordenada13=OrdenadaSBpBG(3:3,1:(i-1));
cadenaylabel='alcanceSBpBGpu';
case{6},% salida 1 : alc...pu, extremo SBp, unidad CG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpCG(1:1,1:(i-1));
Ordenada12=OrdenadaSBpCG(2:2,1:(i-1));
Ordenada13=OrdenadaSBpCG(3:3,1:(i-1));
cadenaylabel='alcanceSBpCGpu';
end %end del switch(unidad)
end %end del switch(extremo)
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'m',Abscisa,Ordenada11,'r',Abscisa,Ordenada12,'g',...
Abscisa,Ordenada13,'b'),grid on
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadenaylabel)
legend('Ref','Lim','k0=0.1','k0=1','k0=10')
% **********************************************************
case{2},% **salida 2** : alc...pua
extremo=input('extremo: SAm=1, SBp=2, (RAm=3), (RBp=4):');
switch(extremo)
case{1},% salida 1 : alc...pu, **extremo Sm**
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenada45(cont)=mm;
cont=cont+1;
mm=mm+0.01;
end
unidad=input('unidad: AB=1, BC=2, CA=3, AG=4, BG=5, CG=6:');
switch(unidad)
case{1},% salida 2 : alc...pua, extremo SAm, **unidad AB**
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmABa(1:1,1:(i-1));
Ordenada12=OrdenadaSAmABa(2:2,1:(i-1));
Ordenada13=OrdenadaSAmABa(3:3,1:(i-1));
cadenaylabel='alcanceSAmABpua';
case{2},% salida 2 : alc...pua, extremo SAm, unidad BC
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmBCa(1:1,1:(i-1));
Ordenada12=OrdenadaSAmBCa(2:2,1:(i-1));
Ordenada13=OrdenadaSAmBCa(3:3,1:(i-1));
cadenaylabel='alcanceSAmBCpua';
case{3},% salida 2 : alc...pua, extremo SAm, unidad CA
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmCAa(1:1,1:(i-1));
Ordenada12=OrdenadaSAmCAa(2:2,1:(i-1));
Ordenada13=OrdenadaSAmCAa(3:3,1:(i-1));
cadenaylabel='alcanceSAmCApua';
case{4},% salida 2 : alc...pua, extremo SAm, unidad AG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmAGa(1:1,1:(i-1));
Ordenada12=OrdenadaSAmAGa(2:2,1:(i-1));
Ordenada13=OrdenadaSAmAGa(3:3,1:(i-1));
cadenaylabel='alcanceSAmAGpua';
case{5},% salida 2 : alc...pua, extremo SAm, unidad BG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmBGa(1:1,1:(i-1));
Ordenada12=OrdenadaSAmBGa(2:2,1:(i-1));
Ordenada13=OrdenadaSAmBGa(3:3,1:(i-1));
cadenaylabel='alcanceSAmBGpua';
case{6},% salida 2 : alc...pua, extremo SAm, unidad CG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSAmCGa(1:1,1:(i-1));
Ordenada12=OrdenadaSAmCGa(2:2,1:(i-1));
Ordenada13=OrdenadaSAmCGa(3:3,1:(i-1));
cadenaylabel='alcanceSAmCGpua';
end %end del switch(unidad)
case{2},% salida 2 : alc...pua, **extremo SBp**
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=0.8;
Ordenada45(cont)=mm;
cont=cont+1;
mm=mm+0.01;
end
unidad=input('unidad: AB=1, BC=2, CA=3, AG=4, BG=5, CG=6:');
switch(unidad)
case{1},% salida 2 : alc...pua, extremo SBp, **unidad AB**
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpABa(1:1,1:(i-1));
Ordenada12=OrdenadaSBpABa(2:2,1:(i-1));
Ordenada13=OrdenadaSBpABa(3:3,1:(i-1));
cadenaylabel='alcanceSBpABpua';
case{2},% salida 2 : alc...pua, extremo SBp, unidad BC
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpBCa(1:1,1:(i-1));
Ordenada12=OrdenadaSBpBCa(2:2,1:(i-1));
Ordenada13=OrdenadaSBpBCa(3:3,1:(i-1));
cadenaylabel='alcanceSBpBCpua';
case{3},% salida 2 : alc...pua, extremo SBp, unidad CA
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpCAa(1:1,1:(i-1));
Ordenada12=OrdenadaSBpCAa(2:2,1:(i-1));
Ordenada13=OrdenadaSBpCAa(3:3,1:(i-1));
cadenaylabel='alcanceSBpCApua';
case{4},% salida 2 : alc...pua, extremo SBp, unidad AG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpAGa(1:1,1:(i-1));
Ordenada12=OrdenadaSBpAGa(2:2,1:(i-1));
Ordenada13=OrdenadaSBpAGa(3:3,1:(i-1));
cadenaylabel='alcanceSBpAGpua';
case{5},% salida 2 : alc...pua, extremo SBp, unidad BG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpBGa(1:1,1:(i-1));
Ordenada12=OrdenadaSBpBGa(2:2,1:(i-1));
Ordenada13=OrdenadaSBpBGa(3:3,1:(i-1));
cadenaylabel='alcanceSBpBGpua';
case{6},% salida 2 : alc...pua, extremo SBp, unidad CG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSBpCGa(1:1,1:(i-1));
Ordenada12=OrdenadaSBpCGa(2:2,1:(i-1));
Ordenada13=OrdenadaSBpCGa(3:3,1:(i-1));
cadenaylabel='alcanceSBpCGpua';
end %end del switch(unidad)
end %end del switch(extremo)
ejemasy=input('Eje +y: ');
ejemenosy=input('Eje -y (con signo): ');
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenada45,'m',Abscisa,Ordenada11,'r',Abscisa,Ordenada12,'g',...
Abscisa,Ordenada13,'b'),grid on
axis([0 1 ejemenosy ejemasy]);
title(cadenamc)
xlabel('m=p'),ylabel(cadenaylabel)
legend('Cob','Ref','k0=0.1','k0=1','k0=10')
end %end del switch(salida)
variassalidas=input('Mas Salidas=1,Acabar=2 :');
end % end del while variassalidas
variasejecuciones=input('Mas Ejecuciones=1,Acabar=2 :');
fprintf('\n');
end % end del while variasejecuciones
fprintf('Pulse Enter\n');
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
