% -----------------------------------------------------------------------
clc
variasejecuciones=1;
while variasejecuciones==1
clear
% -----------------------------------------------------------------------
fprintf('PG124ni\n');
fprintf('Esquema Linea Doble SR\n');
fprintf('Unominal kV: 400,220,132,66\n');
fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n');
fprintf('Impedancia de Transferencia: Si,No\n');
fprintf('Compensacion: Si,No\n');
fprintf('Eleccion: kE,L,Sccs,Sccr,Kr0\n');
fprintf('*** Parametro: k0=zs0/zr0=0.1-1-10 ***\n');
fprintf('Falta Doble, m=p=0...1\n');
fprintf('Analisis Sm,Rm,Sp,Rp\n');
fprintf('zs0=k0*zr0');
fprintf('Salida Grafica: 1(1,2),2(1,2)\n');
fprintf('Salida Grafica 1: Sm,Rm(2 Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf(' : Sp,Rp(2 Subpl), Unidad Pref.(pu), k0=0.1-1-10\n');
fprintf('Salida Grafica 2: Sm,Rm(2 Subpl), Unidad Pref.(pua), k0=0.1-1-10\n');
fprintf(' : Sp,Rp(2 Subpl), Unidad Pref.(pua), k0=0.1-1-10\n');
fprintf('\n');
% -----------------------------------------------------------------------
%PG14ni
datos=input('datos: standard=1,input=2: ');
if datos==1
Un=220;
kE=1;
L=220;
Sccs=5000;
Sccr=5000;
%zsr1=100i;
%Ksr0=3;
imptransfer=2;
compens21=2;
Kr0=1;
% -----------------------------------------------------------------------
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): ');
imptransfer=input('imptransfer:si=1,no=2:');
if imptransfer==1
zsr1=input('zsr1=zsr2(Complejo r+xi en Ohms): ');
zsr2=zsr1;
Ksr0=input('Ksr0=zsr0/zsr1= ');
zsr0=Ksr0*zsr1;
end
compens21=input('Compensacion 21: Si=1, No=2 :');
Kr0=input('Kr0=zr0/zr1: ');
end
% -----------------------------------------------------------------------
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 Rgm=input('Resistencia de Falta en m. Rg(Ohm) = ');
elseif faltam==8
Rgm=input('Resistencia de Falta en m. Rg(Ohm) = ');
elseif faltam==9
Rgm=input('Resistencia de Falta en m. Rg(Ohm) = ');
elseif faltam==11
Rgm=input('Resistencia de Falta en m. Rg(Ohm) = ');
else
Rgm=0;
end
faltap=input('Falta en p: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :');
Rfp=input('Resistencia de Falta en p. Rf(Ohm) = ');
if faltap==7
Rgp=input('Resistencia de Falta en p. Rg(Ohm) = ');
elseif faltap==8
Rgp=input('Resistencia de Falta en p. Rg(Ohm) = ');
elseif faltap==9
Rgp=input('Resistencia de Falta en p. Rg(Ohm) = ');
elseif faltap==11
Rgp=input('Resistencia de Falta en p. Rg(Ohm) = ');
else
Rgp=0;
end
% -----------------------------------------------------------------------
pfaltam=1;
pfaltap=2;
% -----------------------------------------------------------------------
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;
%zsr2=zsr1;
%zsr0=Ksr0*zsr1;
Km=zlM0u/zl1u;
fila=1;
k0=0.1;
while k0<11
zs0=k0*zr0;
i=1;
m=0.01;
while m<1;
p=m;
switch(imptransfer)
case{1},
[Minc]=subMincT(1);
[z0]=subz0T(zs0,zr0,zsr0,zl0,zlM0,m,p);
[z1]=subz1T(zs1,zr1,zsr1,zl1,m,p);
[z2]=subz2T(zs2,zr2,zsr2,zl2,m,p);
case(2)
[Minc]=subMinc(1);
[z0]=subz0(zs0,zr0,zl0,zlM0,m,p);
[z1]=subz1(zs1,zr1,zl1,m,p);
[z2]=subz2(zs2,zr2,zl2,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]=...
mag4mp(Z0,Z1,Z2,E,faltam,faltap,Rfm,Rgm,Rfp,Rgp,pfaltam,pfaltap);
% IBUS(0,1,2)
I0=zeros(4,1);
I1=zeros(4,1);
I2=zeros(4,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);
% 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);
[Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2);
[Ispa,Ispb,Ispc]=abc(Isp0,Isp1,Isp2);
[Irpa,Irpb,Irpc]=abc(Irp0,Irp1,Irp2);
%
% CALCULO DE LAS MEDIDAS DE LAS PROTECCIONES DE DISTANCIA EN
% LOS EXTREMOS DE LINEA
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);
[Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=...
PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k);
[Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=...
PDist(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,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);
[Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0);
[Vspol,Ispop,AngDirSp]=PDire(Vs0,Isp0);
[Vrpol,Irpop,AngDirRp]=PDire(Vr0,Irp0);
if faltam==1
cfaltam='ag';
alcanceSm=abs(Sm_ZAG)/cos(angle(Sm_ZAG)-angle(zl1u));
alcanceRm=abs(Rm_ZAG)/cos(angle(Rm_ZAG)-angle(zl1u));
elementom='AG';
elseif faltam==2
cfaltam='bg';
alcanceSm=abs(Sm_ZBG)/cos(angle(Sm_ZBG)-angle(zl1u));
alcanceRm=abs(Rm_ZBG)/cos(angle(Rm_ZBG)-angle(zl1u));
elementom='BG';
elseif faltam==3
cfaltam='cg';
alcanceSm=abs(Sm_ZCG)/cos(angle(Sm_ZCG)-angle(zl1u));
alcanceRm=abs(Rm_ZCG)/cos(angle(Rm_ZCG)-angle(zl1u));
elementom='CG';
elseif faltam==4
cfaltam='ab';
alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
alcanceRm=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));
elementom='AB';
elseif faltam==5
cfaltam='bc';
alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
elementom='BC';
elseif faltam==6
cfaltam='ca';
alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
alcanceRm=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));
elementom='CA';
elseif faltam==7
cfaltam='abg';
alcanceSm=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u));
alcanceRm=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u));
elementom='AB';
elseif faltam==8
cfaltam='bcg';
alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
elementom='BC';
elseif faltam==9
cfaltam='cag';
alcanceSm=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u));
alcanceRm=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u));
elementom='CA';
elseif faltam==10
cfaltam='abc';
alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
elementom='BC';
elseif faltam==11
cfaltam='abcg';
alcanceSm=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u));
alcanceRm=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u));
elementom='BC';
end
alcanceSmpu=alcanceSm/abs(zl1u*m*L);
alcanceRmpu=alcanceRm/abs(zl1u*(1-m)*L);
alcanceSmpua=alcanceSm/abs(zl1u*L);
alcanceRmpua=alcanceRm/abs(zl1u*L);
if faltap==1
cfaltap='ag';
alcanceSp=abs(Sp_ZAG)/cos(angle(Sp_ZAG)-angle(zl1u));
alcanceRp=abs(Rp_ZAG)/cos(angle(Rp_ZAG)-angle(zl1u));
elementop='AG';
elseif faltap==2
cfaltap='bg';
alcanceSp=abs(Sp_ZBG)/cos(angle(Sp_ZBG)-angle(zl1u));
alcanceRp=abs(Rp_ZBG)/cos(angle(Rp_ZBG)-angle(zl1u));
elementop='BG';
elseif faltap==3
cfaltap='cg';
alcanceSp=abs(Sp_ZCG)/cos(angle(Sp_ZCG)-angle(zl1u));
alcanceRp=abs(Rp_ZCG)/cos(angle(Rp_ZCG)-angle(zl1u));
elementop='CG';
elseif faltap==4
cfaltap='ab';
alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
elementop='AB';
elseif faltap==5
cfaltap='bc';
alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
elementop='BC';
elseif faltap==6
cfaltap='ca';
alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
elementop='CA';
elseif faltap==7
cfaltap='abg';
alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u));
alcanceRp=abs(Rp_ZAB)/cos(angle(Rp_ZAB)-angle(zl1u));
elementop='AB';
elseif faltap==8
cfaltap='bcg';
alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
elementop='BC';
elseif faltap==9
cfaltap='cag';
alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u));
alcanceRp=abs(Rp_ZCA)/cos(angle(Rp_ZCA)-angle(zl1u));
elementop='CA';
elseif faltap==10
cfaltap='abc';
alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
elementop='BC';
elseif faltap==11
cfaltap='abcg';
alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u));
alcanceRp=abs(Rp_ZBC)/cos(angle(Rp_ZBC)-angle(zl1u));
elementop='BC';
end
alcanceSppu=alcanceSp/abs(zl1u*p*L);
alcanceRppu=alcanceRp/abs(zl1u*(1-p)*L);
alcanceSppua=alcanceSp/abs(zl1u*L);
alcanceRppua=alcanceRp/abs(zl1u*L);
Abscisa(fila,i)=m;
OrdenadaSm(fila,i)=alcanceSmpu;
OrdenadaRm(fila,i)=alcanceRmpu;
OrdenadaSp(fila,i)=alcanceSppu;
OrdenadaRp(fila,i)=alcanceRppu;
OrdenadaSma(fila,i)=alcanceSmpua;
OrdenadaRma(fila,i)=alcanceRmpua;
OrdenadaSpa(fila,i)=alcanceSppua;
OrdenadaRpa(fila,i)=alcanceRppua;
i=i+1;
m=m+0.01;
end
fila=fila+1;
k0=k0*10;
end
% -----------------------------------------------------------------------
if compens21==1
comp='21c';
elseif compens21==2
comp='21nc';
end
cadenam=sprintf('PG124ni,%s,%3.0fkV,kE=%-4.2f,%3.0fkm,Sccs=%-5.0f,Sccr=%-5.0f,Kr0=%-5.2f\n',...
comp,Un,kE,L,Sccs,Sccr,Kr0);
if imptransfer==1
cadenamT=sprintf('zsr1=%-6.2fi,Ksr0=%-4.2f,',imag(zsr1),Ksr0);
cadenam=sprintf('%s,%s',cadenam,cadenamT);
end
cadenam1=sprintf('m:%s,%s,Rfm=%-3.0f',cfaltam,elementom,Rfm);
cadenam=sprintf('%s %s',cadenam,cadenam1);
if faltam==7|faltam==8|faltam==9|faltam==11
cadenam2=sprintf('Rgm=%-3.0f',Rgm);
cadenam=sprintf('%s,%s',cadenam,cadenam2);
end
cadenap=sprintf('p: %s,%s,Rfp=%-3.0f',cfaltap,elementop,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);
fprintf('\n');
% -----------------------------------------------------------------------
variassalidas=1;
while variassalidas==1
salida=input('Salida: 1,2 :');
switch(salida)
case{1},
Abscisa1m=Abscisa(1:1,1:(i-1));
Ordenada1m=OrdenadaSm(1:1,1:(i-1));
Abscisa2m=Abscisa(2:2,1:(i-1));
Ordenada2m=OrdenadaSm(2:2,1:(i-1));
Abscisa3m=Abscisa(3:3,1:(i-1));
Ordenada3m=OrdenadaSm(3:3,1:(i-1));
Abscisa4m=Abscisa(1:1,1:(i-1));
Ordenada4m=OrdenadaRm(1:1,1:(i-1));
Abscisa5m=Abscisa(2:2,1:(i-1));
Ordenada5m=OrdenadaRm(2:2,1:(i-1));
Abscisa6m=Abscisa(3:3,1:(i-1));
Ordenada6m=OrdenadaRm(3:3,1:(i-1));
Abscisa1p=Abscisa(1:1,1:(i-1));
Ordenada1p=OrdenadaSp(1:1,1:(i-1));
Abscisa2p=Abscisa(2:2,1:(i-1));
Ordenada2p=OrdenadaSp(2:2,1:(i-1));
Abscisa3p=Abscisa(3:3,1:(i-1));
Ordenada3p=OrdenadaSp(3:3,1:(i-1));
Abscisa4p=Abscisa(1:1,1:(i-1));
Ordenada4p=OrdenadaRp(1:1,1:(i-1));
Abscisa5p=Abscisa(2:2,1:(i-1));
Ordenada5p=OrdenadaRp(2:2,1:(i-1));
Abscisa6p=Abscisa(3:3,1:(i-1));
Ordenada6p=OrdenadaRp(3:3,1:(i-1));
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=1.0;
Ordenadadi(cont)=0.8/mm;
Ordenadadd(cont)=0.8/(1-mm);
cont=cont+1;
mm=mm+0.01;
end
AnalisisLinea=input('Extremos Sm,Rm=1, Sp,Rp=2: ');
switch(AnalisisLinea)
case{1},
ejemasy1=input('Subplot 1: Eje +y: ');
ejemenosy1=input('Subplot 1: Eje -y (con signo): ');
ejemasy2=input('Subplot 2: Eje +y: ');
ejemenosy2=input('Subplot 2: Eje -y (con signo): ');
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'m',Abscisa1m,Ordenada1m,'r',Abscisa2m,Ordenada2m,'g',...
Abscisa3m,Ordenada3m,'b'),grid on
axis([0 1 ejemenosy1 ejemasy1]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSmpu')
%legend('Ref','Lim','k0=0.1','k0=1','k0=10')
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadd,'m',Abscisa4m,Ordenada4m,'r',Abscisa5m,Ordenada5m,'g',...
Abscisa6m,Ordenada6m,'b'),grid on
axis([0 1 ejemenosy2 ejemasy2]);
%title(cadenamp)
xlabel('m=p'),ylabel('alcanceRmpu')
legend('Ref','Lim','k0=0.1','k0=1','k0=10')
case{2},
ejemasy1=input('Subplot 1: Eje +y: ');
ejemenosy1=input('Subplot 1: Eje -y (con signo): ');
ejemasy2=input('Subplot 2: Eje +y: ');
ejemenosy2=input('Subplot 2: Eje -y (con signo): ');
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'m',Abscisa1p,Ordenada1p,'r',Abscisa2p,Ordenada2p,'g',...
Abscisa3p,Ordenada3p,'b'),grid on
axis([0 1 ejemenosy1 ejemasy1]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSppu')
%legend('Ref','Lim','k0=0.1','k0=1','k0=10')
subplot(2,1,2)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadd,'m',Abscisa4p,Ordenada4p,'r',Abscisa5p,Ordenada5p,'g',...
Abscisa6p,Ordenada6p,'b'),grid on
axis([0 1 ejemenosy2 ejemasy2]);
%title(cadenamp)
xlabel('m=p'),ylabel('alcanceRppu')
legend('Ref','Lim','k0=0.1','k0=1','k0=10')
end
case(2)
mm=0.01;
cont=1;
while mm<1;
Abscisa0(cont)=mm;
Ordenada0(cont)=mm;
Abscisa00(cont)=mm;
Ordenada00(cont)=1-mm;
Ordenadaz1(cont)=0.8;
cont=cont+1;
mm=mm+0.01;
end
Abscisa1ma=Abscisa(1:1,1:(i-1));
Ordenada1ma=OrdenadaSma(1:1,1:(i-1));
Abscisa2ma=Abscisa(2:2,1:(i-1));
Ordenada2ma=OrdenadaSma(2:2,1:(i-1));
Abscisa3ma=Abscisa(3:3,1:(i-1));
Ordenada3ma=OrdenadaSma(3:3,1:(i-1));
Abscisa4ma=Abscisa(1:1,1:(i-1));
Ordenada4ma=OrdenadaRma(1:1,1:(i-1));
Abscisa5ma=Abscisa(2:2,1:(i-1));
Ordenada5ma=OrdenadaRma(2:2,1:(i-1));
Abscisa6ma=Abscisa(3:3,1:(i-1));
Ordenada6ma=OrdenadaRma(3:3,1:(i-1));
Abscisa1pa=Abscisa(1:1,1:(i-1));
Ordenada1pa=OrdenadaSpa(1:1,1:(i-1));
Abscisa2pa=Abscisa(2:2,1:(i-1));
Ordenada2pa=OrdenadaSpa(2:2,1:(i-1));
Abscisa3pa=Abscisa(3:3,1:(i-1));
Ordenada3pa=OrdenadaSpa(3:3,1:(i-1));
Abscisa4pa=Abscisa(1:1,1:(i-1));
Ordenada4pa=OrdenadaRpa(1:1,1:(i-1));
Abscisa5pa=Abscisa(2:2,1:(i-1));
Ordenada5pa=OrdenadaRpa(2:2,1:(i-1));
Abscisa6pa=Abscisa(3:3,1:(i-1));
Ordenada6pa=OrdenadaRpa(3:3,1:(i-1));
AnalisisLineaa=input('Extremos Sm,Rm=1, Sp,Rp=2: ');
switch(AnalisisLineaa)
case{1},
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadaz1,'m',Abscisa1ma,Ordenada1ma,'r',Abscisa2ma,Ordenada2ma,'g',...
Abscisa3ma,Ordenada3ma,'b'),grid on
%axis([0 1 0 2]);
title(cadenamp);
xlabel('m=p'),ylabel('alcanceSmpua')
%legend('45','Lim','k0=0.1','k0=1','k0=10')
subplot(2,1,2)
plot(Abscisa00,Ordenada00,'c',Abscisa0,Ordenadaz1,'m',Abscisa4ma,Ordenada4ma,'r',Abscisa5ma,Ordenada5ma,'g',...
Abscisa6ma,Ordenada6ma,'b'),grid on
%axis([0 1 0 2]);
%title(cadenamp)
xlabel('m=p'),ylabel('alcanceRmpua')
legend('Ref','Lim','k0=0.1','k0=1','k0=10')
case{2},
subplot(2,1,1)
plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadaz1,'m',Abscisa1pa,Ordenada1pa,'r',Abscisa2pa,Ordenada2pa,'g',...
Abscisa3pa,Ordenada3pa,'b'),grid on
%axis([0 1 0 2]);
title(cadenamp)
xlabel('m=p'),ylabel('alcanceSppua')
%legend('45','Lim','k0=0.1','k0=1','k0=10')
subplot(2,1,2)
plot(Abscisa00,Ordenada00,'c',Abscisa0,Ordenadaz1,'m',Abscisa4pa,Ordenada4pa,'r',Abscisa5pa,Ordenada5pa,'g',...
Abscisa6pa,Ordenada6pa,'b'),grid on
%axis([0 1 0 2]);
%title(cadenamp)
xlabel('m=p'),ylabel('alcanceRppua')
legend('Ref','Lim','k0=0.1','k0=1','k0=10')
end
end
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.
