% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('\n');
fprintf('PG126MTf\n');
fprintf('Esquema Linea Doble SR\n');
fprintf('Unominal kV: 400,220,132,66\n');
fprintf('Eleccion: kE,L,Sccs,Sccr,Kr0,k0)\n');
fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf('*** Parametro zsr1=10i-100i-1000i ***\n');
fprintf('Elecion Ksr0=zsr0/zsr1');
fprintf('Compensacion 21, Si, No\n');
fprintf('Falta Intercircuito entre Lineas de igual Tension\n');
fprintf('m=p=0...1\n');
fprintf('Barra m=1, Barra p=2\n');
%fprintf('Salida Grafica 1: Analisis Sm, Elegir Unidad (pu)\n');
%fprintf(' : Analisis Sp, Elegir Unidad (pu)\n');
%fprintf('Salida Grafica 2: Analisis Sm, Elegir Unidad (pua)\n');
%fprintf(' : Analisis Sp, Elegir Unidad (pua)\n');
fprintf('\n');
% ----------------------------------------------------------------------------------
%PG126Mc
datos=input('datos: standard=1,input=2: ');
if datos==1
Un=220;
kE=1;
L=220;
Sccs=5000;
Sccr=5000;
Kr0=1;
zsr1=10i;
zsr2=zsr1;
Ksr0=3;
zsr0=Ksr0*zsr1;
imptransfer=1;
compens21=2;
else
% ----------------------------------------------------------------------------------
Un=input('U(Fase-Fase) kV:400,220,132,66: ');
kE=input('kE=E/Enominal: ');
L=input('Longitud Linea km: ');
Sccs=input('Sccs(MVA): ');
Sccr=input('Sccr(MVA): ');
Kr0=input('Kr0=zr0/zr1: ');
imptransfer=1;
zsr1=10i;
%zsr1=input('zsr1=zsr2(Complejo r+xi en Ohms): ');
zsr2=zsr1;
Ksr0=input('Ksr0=zsr0/zsr1= ');
zsr0=Ksr0*zsr1;
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
% ----------------------------------------------------------------------------------
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);
zl1=zl1u*L; zl2=zl2u*L; zl0=zl0u*L; zlM0=zlM0u*L;
zs1=(Un^2/Sccs)*i;
zs2=zs1;
zr1=(Un^2/Sccr)*i;
zr2=zr1;
zr0=Kr0*zr1;
%zs0=k0*zr0;
Km=zlM0u/zl1u;
% ----------------------------------------------------------------------------------
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
% ----------------------------------------------------------------------------------
fprintf('\n');
fila=1;
zs0=k0*zr0;
while abs(zsr1)<1001
i=1;
m=0.01;
while m<1;
p=m;
%
switch(imptransfer)
case{1},
[Minc]=subMincT(1);
[z0]=subz0T(zs0,zr0,zsr0,zl0u,zlM0u,L,m,p);
[z1]=subz1T(zs1,zr1,zsr1,zl1u,L,m,p);
[z2]=subz2T(zs2,zr2,zsr2,zl2u,L,m,p);
case(2)
[Minc]=subMinc(1);
[z0]=subz0(zs0,zr0,zl0u,zlM0u,L,m,p);
[z1]=subz1(zs1,zr1,zl1u,L,m,p);
[z2]=subz2(zs2,zr2,zl2u,L,m,p);
end
% MATRIZ DE ADMITANCIA DE BARRAS
[Y0]=subY0(Minc,z0);
[Y1]=subY1(Minc,z1);
[Y2]=subY2(Minc,z2);
% ZBUS(0,1,2)
Z0=inv(Y0);
Z1=inv(Y1);
Z2=inv(Y2);
% 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]=...
maMmp(Z0,Z1,Z2,E,falta,Rf,Rg);
% IBUS(0,1,2)
I0=[-Im0;-Ip0;0;0];
I1=[-Im1;-Ip1;0;0];
I2=[-Im2;-Ip2;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(Minc,V0,z0);
[Ie1]=subIe1(Minc,V1,z1);
[Ie2]=subIe2(Minc,V2,z2);
% V(0,1,2), I(0,1,2) EN LOS EXTREMOS DE LINEA
Vs0=V0(3); Vr0=V0(4);
Vs1=V1(3); Vr1=V1(4);
Vs2=V2(3); Vr2=V2(4);
Ism0=Ie0(3); Irm0=Ie0(5); Isp0=Ie0(4); Irp0=Ie0(6);
Ism1=Ie1(3); Irm1=Ie1(5); Isp1=Ie1(4); Irp1=Ie1(6);
Ism2=Ie2(3); Irm2=Ie2(5); Isp2=Ie2(4); Irp2=Ie2(6);
% 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
[Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2);
[Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2);
[Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2);
[Ispa,Ispb,Ispc]=abc(Isp0,Isp1,Isp2);
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2);
[Irpa,Irpb,Irpc]=abc(Irp0,Irp1,Irp2);
%
% PROTECCION DE DISTANCIA (21)
switch(compens21)
case{1},
% COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k,Km,Isp0);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
PDistComp(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k,Km,Ism0);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k,Km,Irp0);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
PDistComp(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k,Km,Irm0);
case{2},
% NO COMPENSADA
[Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=...
PDist(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
PDist(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k);
[Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=...
PDist(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k);
end
% CALCULO DE LAS MAGNITUDES DE ALIMETACION DE LA PROTECCION
% DIRECCIONAL DE TIERRA EN LOS EXTREMOS DE LINEA
[Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0); % kV, kA, grados
[Vspol,Ispop,AngDirSp]=PDire(Vs0,Isp0);
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0);
[Vrpol,Irpop,AngDirRp]=PDire(Vr0,Irp0);
% para Salida 1 y 2
alcSmAB=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
alcSmBC=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcSmCA=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
alcSmAG=abs(Sm_ZAG)/cos(angle(Sm_ZAG)-angle(zl1u));
alcSmBG=abs(Sm_ZBG)/cos(angle(Sm_ZBG)-angle(zl1u));
alcSmCG=abs(Sm_ZCG)/cos(angle(Sm_ZCG)-angle(zl1u));
alcSpAB=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
alcSpBC=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcSpCA=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
alcSpAG=abs(Sp_ZAG)/cos(angle(Sp_ZAG)-angle(zl1u));
alcSpBG=abs(Sp_ZBG)/cos(angle(Sp_ZBG)-angle(zl1u));
alcSpCG=abs(Sp_ZCG)/cos(angle(Sp_ZCG)-angle(zl1u));
alcRmAB=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));
alcRmBC=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
alcRmCA=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));
alcRmAG=abs(Rm_ZAG)/cos(angle(Rm_ZAG)-angle(zl1u));
alcRmBG=abs(Rm_ZBG)/cos(angle(Rm_ZBG)-angle(zl1u));
alcRmCG=abs(Rm_ZCG)/cos(angle(Rm_ZCG)-angle(zl1u));
alcRpAB=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
alcRpBC=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
alcRpCA=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
alcRpAG=abs(Rp_ZAG)/cos(angle(Rp_ZAG)-angle(zl1u));
alcRpBG=abs(Rp_ZBG)/cos(angle(Rp_ZBG)-angle(zl1u));
alcRpCG=abs(Rp_ZCG)/cos(angle(Rp_ZCG)-angle(zl1u));
% para Salida 1
alcSmABpu=alcSmAB/abs(zl1u*m*L);
alcSmBCpu=alcSmBC/abs(zl1u*m*L);
alcSmCApu=alcSmCA/abs(zl1u*m*L);
alcSmAGpu=alcSmAG/abs(zl1u*m*L);
alcSmBGpu=alcSmBG/abs(zl1u*m*L);
alcSmCGpu=alcSmCG/abs(zl1u*m*L);
alcSpABpu=alcSpAB/abs(zl1u*m*L);
alcSpBCpu=alcSpBC/abs(zl1u*m*L);
alcSpCApu=alcSpCA/abs(zl1u*m*L);
alcSpAGpu=alcSpAG/abs(zl1u*m*L);
alcSpBGpu=alcSpBG/abs(zl1u*m*L);
alcSpCGpu=alcSpCG/abs(zl1u*m*L);
alcRmABpu=alcRmAB/abs(zl1u*(1-m)*L);
alcRmBCpu=alcRmBC/abs(zl1u*(1-m)*L);
alcRmCApu=alcRmCA/abs(zl1u*(1-m)*L);
alcRmAGpu=alcRmAG/abs(zl1u*(1-m)*L);
alcRmBGpu=alcRmBG/abs(zl1u*(1-m)*L);
alcRmCGpu=alcRmCG/abs(zl1u*(1-m)*L);
alcRpABpu=alcRpAB/abs(zl1u*(1-m)*L);
alcRpBCpu=alcRpBC/abs(zl1u*(1-m)*L);
alcRpCApu=alcRpCA/abs(zl1u*(1-m)*L);
alcRpAGpu=alcRpAG/abs(zl1u*(1-m)*L);
alcRpBGpu=alcRpBG/abs(zl1u*(1-m)*L);
alcRpCGpu=alcRpCG/abs(zl1u*(1-m)*L);
cadylabel1='alcSmpu';
cadylabel2='alcSppu';
cadylabel3='alcRmpu';
cadylabel4='alcRppu';
% para Salida 2
alcSmABpua=alcSmAB/abs(zl1u*L);
alcSmBCpua=alcSmBC/abs(zl1u*L);
alcSmCApua=alcSmCA/abs(zl1u*L);
alcSmAGpua=alcSmAG/abs(zl1u*L);
alcSmBGpua=alcSmBG/abs(zl1u*L);
alcSmCGpua=alcSmCG/abs(zl1u*L);
alcSpABpua=alcSpAB/abs(zl1u*L);
alcSpBCpua=alcSpBC/abs(zl1u*L);
alcSpCApua=alcSpCA/abs(zl1u*L);
alcSpAGpua=alcSpAG/abs(zl1u*L);
alcSpBGpua=alcSpBG/abs(zl1u*L);
alcSpCGpua=alcSpCG/abs(zl1u*L);
alcRmABpua=alcRmAB/abs(zl1u*L);
alcRmBCpua=alcRmBC/abs(zl1u*L);
alcRmCApua=alcRmCA/abs(zl1u*L);
alcRmAGpua=alcRmAG/abs(zl1u*L);
alcRmBGpua=alcRmBG/abs(zl1u*L);
alcRmCGpua=alcRmCG/abs(zl1u*L);
alcRpABpua=alcRpAB/abs(zl1u*L);
alcRpBCpua=alcRpBC/abs(zl1u*L);
alcRpCApua=alcRpCA/abs(zl1u*L);
alcRpAGpua=alcRpAG/abs(zl1u*L);
alcRpBGpua=alcRpBG/abs(zl1u*L);
alcRpCGpua=alcRpCG/abs(zl1u*L);
cadylabel1a='alcSmpua';
cadylabel2a='alcSppua';
cadylabel3a='alcRmpua';
cadylabel4a='alcRppua';
Abscisa(fila,i)=m;
% para Salida 1
OrdenadaSmAB(fila,i)=alcSmABpu;
OrdenadaSmBC(fila,i)=alcSmBCpu;
OrdenadaSmCA(fila,i)=alcSmCApu;
OrdenadaSmAG(fila,i)=alcSmAGpu;
OrdenadaSmBG(fila,i)=alcSmBGpu;
OrdenadaSmCG(fila,i)=alcSmCGpu;
OrdenadaSpAB(fila,i)=alcSpABpu;
OrdenadaSpBC(fila,i)=alcSpBCpu;
OrdenadaSpCA(fila,i)=alcSpCApu;
OrdenadaSpAG(fila,i)=alcSpAGpu;
OrdenadaSpBG(fila,i)=alcSpBGpu;
OrdenadaSpCG(fila,i)=alcSpCGpu;
OrdenadaRmAB(fila,i)=alcRmABpu;
OrdenadaRmBC(fila,i)=alcRmBCpu;
OrdenadaRmCA(fila,i)=alcRmCApu;
OrdenadaRmAG(fila,i)=alcRmAGpu;
OrdenadaRmBG(fila,i)=alcRmBGpu;
OrdenadaRmCG(fila,i)=alcRmCGpu;
OrdenadaRpAB(fila,i)=alcRpABpu;
OrdenadaRpBC(fila,i)=alcRpBCpu;
OrdenadaRpCA(fila,i)=alcRpCApu;
OrdenadaRpAG(fila,i)=alcRpAGpu;
OrdenadaRpBG(fila,i)=alcRpBGpu;
OrdenadaRpCG(fila,i)=alcRpCGpu;
% para Salida 2
OrdenadaSmABa(fila,i)=alcSmABpua;
OrdenadaSmBCa(fila,i)=alcSmBCpua;
OrdenadaSmCAa(fila,i)=alcSmCApua;
OrdenadaSmAGa(fila,i)=alcSmAGpua;
OrdenadaSmBGa(fila,i)=alcSmBGpua;
OrdenadaSmCGa(fila,i)=alcSmCGpua;
OrdenadaSpABa(fila,i)=alcSpABpua;
OrdenadaSpBCa(fila,i)=alcSpBCpua;
OrdenadaSpCAa(fila,i)=alcSpCApua;
OrdenadaSpAGa(fila,i)=alcSpAGpua;
OrdenadaSpBGa(fila,i)=alcSpBGpua;
OrdenadaSpCGa(fila,i)=alcSpCGpua;
OrdenadaRmABa(fila,i)=alcRmABpua;
OrdenadaRmBCa(fila,i)=alcRmBCpua;
OrdenadaRmCAa(fila,i)=alcRmCApua;
OrdenadaRmAGa(fila,i)=alcRmAGpua;
OrdenadaRmBGa(fila,i)=alcRmBGpua;
OrdenadaRmCGa(fila,i)=alcRmCGpua;
OrdenadaRpABa(fila,i)=alcRpABpua;
OrdenadaRpBCa(fila,i)=alcRpBCpua;
OrdenadaRpCAa(fila,i)=alcRpCApua;
OrdenadaRpAGa(fila,i)=alcRpAGpua;
OrdenadaRpBGa(fila,i)=alcRpBGpua;
OrdenadaRpCGa(fila,i)=alcRpCGpua;
%***************************************
i=i+1;
m=m+0.01;
end %(while m)
fila=fila+1;
zsr1=zsr1*10;
zsr0=zsr1*Ksr0;
end % FIN DE BUCLES zsr1, m
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
%***************************************
% Titulo para las Salidas 1,2
%***************************************
cadenamc1=sprintf('PG126MTf,%s,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,Kr0=%-5.2f,k0=%-5.2f\n',...
comp,Un,kE,L,Sccs,Sccr,Kr0,k0);
cadenamT=sprintf('Ksr0=%-4.2f',Ksr0);
cadenamc1=sprintf('%s,%s',cadenamc1,cadenamT);
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('Salida Grafica 1: Analisis Sm, Elegir Unidad (pu)\n');
fprintf(' : Sp, " \n');
fprintf('Salida Grafica 2: Analisis Sm, Elegir Unidad (pua)\n');
fprintf(' : Sp, " \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: Sm=1, Sp=2, (Rm=3), (Rp=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 Sm, **unidad AB**
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmAB(1:1,1:(i-1));
Ordenada12=OrdenadaSmAB(2:2,1:(i-1));
Ordenada13=OrdenadaSmAB(3:3,1:(i-1));
cadenaylabel='alcanceSmABpu';
case{2},% salida 1 : alc...pu, extremo Sm, unidad BC
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmBC(1:1,1:(i-1));
Ordenada12=OrdenadaSmBC(2:2,1:(i-1));
Ordenada13=OrdenadaSmBC(3:3,1:(i-1));
cadenaylabel='alcanceSmBCpu';
case{3},% salida 1 : alc...pu, extremo Sm, unidad CA
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmCA(1:1,1:(i-1));
Ordenada12=OrdenadaSmCA(2:2,1:(i-1));
Ordenada13=OrdenadaSmCA(3:3,1:(i-1));
cadenaylabel='alcanceSmCApu';
case{4},% salida 1 : alc...pu, extremo Sm, unidad AG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmAG(1:1,1:(i-1));
Ordenada12=OrdenadaSmAG(2:2,1:(i-1));
Ordenada13=OrdenadaSmAG(3:3,1:(i-1));
cadenaylabel='alcanceSmAGpu';
case{5},% salida 1 : alc...pu, extremo Sm, unidad BG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmBG(1:1,1:(i-1));
Ordenada12=OrdenadaSmBG(2:2,1:(i-1));
Ordenada13=OrdenadaSmBG(3:3,1:(i-1));
cadenaylabel='alcanceSmBGpu';
case{6},% salida 1 : alc...pu, extremo Sm, unidad CG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmCG(1:1,1:(i-1));
Ordenada12=OrdenadaSmCG(2:2,1:(i-1));
Ordenada13=OrdenadaSmCG(3:3,1:(i-1));
cadenaylabel='alcanceSmCGpu';
end %end del switch(unidad)
case{2},% salida 1 : alc...pu, **extremo Sp**
unidad=input('unidad: AB=1, BC=2, CA=3, AG=4, BG=5, CG=6:');
switch(unidad)
case{1},% salida 1 : alc...pu, extremo Sp, **unidad AB**
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpAB(1:1,1:(i-1));
Ordenada12=OrdenadaSpAB(2:2,1:(i-1));
Ordenada13=OrdenadaSpAB(3:3,1:(i-1));
cadenaylabel='alcanceSpABpu';
case{2},% salida 1 : alc...pu, extremo Sp, unidad BC
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpBC(1:1,1:(i-1));
Ordenada12=OrdenadaSpBC(2:2,1:(i-1));
Ordenada13=OrdenadaSpBC(3:3,1:(i-1));
cadenaylabel='alcanceSpBCpu';
case{3},% salida 1 : alc...pu, extremo Sp, unidad CA
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpCA(1:1,1:(i-1));
Ordenada12=OrdenadaSpCA(2:2,1:(i-1));
Ordenada13=OrdenadaSpCA(3:3,1:(i-1));
cadenaylabel='alcanceSpCApu';
case{4},% salida 1 : alc...pu, extremo Sp, unidad AG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpAG(1:1,1:(i-1));
Ordenada12=OrdenadaSpAG(2:2,1:(i-1));
Ordenada13=OrdenadaSpAG(3:3,1:(i-1));
cadenaylabel='alcanceSpAGpu';
case{5},% salida 1 : alc...pu, extremo Sp, unidad BG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpBG(1:1,1:(i-1));
Ordenada12=OrdenadaSpBG(2:2,1:(i-1));
Ordenada13=OrdenadaSpBG(3:3,1:(i-1));
cadenaylabel='alcanceSpBGpu';
case{6},% salida 1 : alc...pu, extremo Sp, unidad CG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpCG(1:1,1:(i-1));
Ordenada12=OrdenadaSpCG(2:2,1:(i-1));
Ordenada13=OrdenadaSpCG(3:3,1:(i-1));
cadenaylabel='alcanceSpCGpu';
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','zsr1=10i','zsr1=100i','zsr1=1000i')
%******************************************************
case{2},% **salida 2** : alc...pua
extremo=input('extremo: Sm=1, Sp=2, (Rm=3), (Rp=4):');
switch(extremo)
case{1},% salida 2 : alc...pua, **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 Sm, **unidad AB**
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmABa(1:1,1:(i-1));
Ordenada12=OrdenadaSmABa(2:2,1:(i-1));
Ordenada13=OrdenadaSmABa(3:3,1:(i-1));
cadenaylabel='alcanceSmABpua';
case{2},% salida 2 : alc...pua, extremo Sm, unidad BC
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmBCa(1:1,1:(i-1));
Ordenada12=OrdenadaSmBCa(2:2,1:(i-1));
Ordenada13=OrdenadaSmBCa(3:3,1:(i-1));
cadenaylabel='alcanceSmBCpua';
case{3},% salida 2 : alc...pua, extremo Sm, unidad CA
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmCAa(1:1,1:(i-1));
Ordenada12=OrdenadaSmCAa(2:2,1:(i-1));
Ordenada13=OrdenadaSmCAa(3:3,1:(i-1));
cadenaylabel='alcanceSmCApua';
case{4},% salida 2 : alc...pua, extremo Sm, unidad AG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmAGa(1:1,1:(i-1));
Ordenada12=OrdenadaSmAGa(2:2,1:(i-1));
Ordenada13=OrdenadaSmAGa(3:3,1:(i-1));
cadenaylabel='alcanceSmAGpua';
case{5},% salida 2 : alc...pua, extremo Sm, unidad BG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmBGa(1:1,1:(i-1));
Ordenada12=OrdenadaSmBGa(2:2,1:(i-1));
Ordenada13=OrdenadaSmBGa(3:3,1:(i-1));
cadenaylabel='alcanceSmBGpua';
case{6},% salida 2 : alc...pua, extremo Sm, unidad CG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSmCGa(1:1,1:(i-1));
Ordenada12=OrdenadaSmCGa(2:2,1:(i-1));
Ordenada13=OrdenadaSmCGa(3:3,1:(i-1));
cadenaylabel='alcanceSmCGpua';
end %end del switch(unidad)
case{2},% salida 2 : alc...pua, **extremo Sp**
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 Sp, **unidad AB**
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpABa(1:1,1:(i-1));
Ordenada12=OrdenadaSpABa(2:2,1:(i-1));
Ordenada13=OrdenadaSpABa(3:3,1:(i-1));
cadenaylabel='alcanceSpABpua';
case{2},% salida 2 : alc...pua, extremo Sp, unidad BC
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpBCa(1:1,1:(i-1));
Ordenada12=OrdenadaSpBCa(2:2,1:(i-1));
Ordenada13=OrdenadaSpBCa(3:3,1:(i-1));
cadenaylabel='alcanceSpBCpua';
case{3},% salida 2 : alc...pua, extremo Sp, unidad CA
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpCAa(1:1,1:(i-1));
Ordenada12=OrdenadaSpCAa(2:2,1:(i-1));
Ordenada13=OrdenadaSpCAa(3:3,1:(i-1));
cadenaylabel='alcanceSpCApua';
case{4},% salida 2 : alc...pua, extremo Sp, unidad AG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpAGa(1:1,1:(i-1));
Ordenada12=OrdenadaSpAGa(2:2,1:(i-1));
Ordenada13=OrdenadaSpAGa(3:3,1:(i-1));
cadenaylabel='alcanceSpAGpua';
case{5},% salida 2 : alc...pua, extremo Sp, unidad BG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpBGa(1:1,1:(i-1));
Ordenada12=OrdenadaSpBGa(2:2,1:(i-1));
Ordenada13=OrdenadaSpBGa(3:3,1:(i-1));
cadenaylabel='alcanceSpBGpua';
case{6},% salida 2 : alc...pua, extremo Sp, unidad CG
Abscisa=Abscisa(1:1,1:(i-1));
Ordenada11=OrdenadaSpCGa(1:1,1:(i-1));
Ordenada12=OrdenadaSpCGa(2:2,1:(i-1));
Ordenada13=OrdenadaSpCGa(3:3,1:(i-1));
cadenaylabel='alcanceSpCGpua';
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','zsr1=10i','zsr1=100i','zsr1=1000i')
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.
