clc variasejecuciones=1; while variasejecuciones==1 clear % -------------------------------------------------- fprintf('PG34ic\n'); fprintf('Esquema Linea Doble: SR,ST\n'); fprintf('Barras separadas(R,T) en el Extremo remoto\n'); fprintf('Unominal kV: 400,220,132,66\n'); fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n'); fprintf('Compensacion 21: Si,No\n'); fprintf('Eleccion: kE,L,Sccs,Sccr,Scct,Kr0,Ks0\n'); fprintf('Eleccion: k0=zt0/zs0\n'); fprintf('Falta Doble, m=p=0...1\n'); fprintf('Compensacion 21, Si,No\n'); fprintf('Analisis Sm,Rm\n'); %fprintf('Salida Grafica 1: Sm,Rm(2 Subpl),Unidad Pref.(pu)\n'); %fprintf(' Sp,Tp(2 Subpl),Unidad Pref.(pu)\n'); %fprintf('Salida Grafica 2: Sm,Rm(2 Subpl),Unidad Pref.(pua)\n'); %fprintf(' Sp,Tp(2 Subpl),Unidad Pref.(pua)\n'); %fprintf('Salida Grafica 3: Analisis Sm, Todas las Unidades (pua)\n'); %fprintf(' : Analisis Rm, Todas las Unidades (pua)\n'); %fprintf(' : Analisis Sp, Todas las Unidades (pua)\n'); %fprintf(' : Analisis Tp, Todas las Unidades (pua)\n'); %fprintf('Salida Grafica 4: Analisis Sm, 67N\n'); %fprintf(' : Analisis Rm, 67N\n'); %fprintf(' : Analisis Sp, 67N\n'); %fprintf(' : Analisis Tp, 67N\n'); %fprintf('Salida Grafica 5: Analisis Sm, R-X\n'); %fprintf(' : Analisis Rm, R-X\n'); %fprintf(' : Analisis Sp, R-X\n'); %fprintf(' : Analisis Tp, R-X\n'); fprintf('\n'); % -------------------------------------------------- %PG34ic datos=input('datos: standard=1,input=2: '); if datos==1 Un=220; kE=1; L=220; compens21=2; Sccs=1000; Sccr=1000; Scct=5000; Ks0=1; Kr0=1; % -------------------------------------------------- else Un=input('Un kV: 400,220,132,66 :'); kE=input('kE=E/Enominal: '); L=input('L km: '); compens21=input('Compensacion 21: Si=1, No=2 :'); Sccs=input('Sccs(MVA): '); Sccr=input('Sccr(MVA): '); Scct=input('Scct(MVA): '); Ks0=input('Ks0: '); Kr0=input('Kr0: '); end % -------------------------------------------------- k0=input('k0=zt0/zs0: '); 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 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|faltap==8|faltap==9|faltap==11 Rgp=input('Resistencia de Falta en p. Rg(Ohm) = '); else Rgp=0; end fprintf('\n'); % ----------------------------------------------------------------- pfaltam=4; pfaltap=5; % ----------------------------------------------------------------- % DATOS % ----------------------------------------------------------------- Ubase=Un; % kV Vbase=Ubase/(sqrt(3)); % kV MVAbase=100; Ibase=MVAbase/(3*Vbase); % kA zbase=Vbase/Ibase; % Ohm pzona=0.8; if Un==400 zl1u=(0.0346+0.3208i)/zbase; zl2u=zl1u; zl0u=(0.1828+1.1482i)/zbase; zlM0u=(0.1482+0.7383i)/zbase; elseif Un==220 zl1u=(0.0829+0.4135i)/zbase; zl2u=zl1u; zl0u=(0.2311+1.3085i)/zbase; zlM0u=(0.1482+0.8031i)/zbase; elseif Un==132 zl1u=(0.1342+0.4042i)/zbase; zl2u=zl1u; zl0u=(0.2824+1.3718i)/zbase; zlM0u=(0.1482+0.8947i)/zbase; elseif Un==66 zl1u=(0.2156+0.3952i)/zbase; zl2u=zl1u; zl0u=(0.3636+1.4278i)/zbase; zlM0u=(0.1482+0.9288i)/zbase; end % --------------------------------------------------- %CALCULOS PRELIMINARES % --------------------------------------------------- U=kE*Un; % kV E=U/sqrt(3)+0i; % kV E=E/Vbase; % E pu k=(zl0u-zl1u)/(3*zl1u); Km=zlM0u/zl1u; zs1=(MVAbase/Sccs)*i; % pu zr1=(MVAbase/Sccr)*i; % pu zt1=(MVAbase/Scct)*i; % pu zs2=zs1; % pu zr2=zr1; % pu zt2=zt1; % pu zs0=Ks0*zs1; % pu zr0=Kr0*zr1; % pu zt0=k0*zs0; % pu zl0=zl0u*L; zl1=zl1u*L; zl2=zl2u*L; zlM0=zlM0u*L; % pu % ---------------------------------------------------- % CALCULOS % ----------------------------------------------------- i=1; m=0.01; while m<1; p=m; [Minc]=subMinc(1); [Minca]=subMinca(1); [z0]=subz0(zs0,zr0,zt0,zl0,zlM0,m,p); [z1]=subz1(zs1,zr1,zt1,zl1,m,p); [z2]=subz2(zs2,zr2,zt2,zl2,m,p); % MATRIZ DE ADMITANCIA DE BARRAS [Y0]=subY0(Minc,z0); % pu [Y1]=subY1(Minc,z1); % pu [Y2]=subY2(Minc,z2); % pu % ZBUS(0,1,2) Z0=inv(Y0); % pu Z1=inv(Y1); % pu Z2=inv(Y2); % pu % 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); 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 % IBUS(0,1,2) I0=zeros(5,1); I1=zeros(5,1); I2=zeros(5,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; V1a=[V1;E]; V2=Z2*I2; % Ie(0,1,2) DE LOS ELEMENTOS [Ie0]=subIe0(Minc,V0,z0); [Ie1]=subIe1(Minca,V1a,z1); [Ie2]=subIe2(Minc,V2,z2); % ASIGNACIONES V,I Vs0=V0(1);Vs1=V1(1);Vs2=V2(1); Vr0=V0(2);Vr1=V1(2);Vr2=V2(2); Vt0=V0(3);Vt1=V1(3);Vt2=V2(3); Is0=Ie0(1);Is1=Ie1(1);Is2=Ie2(1); Ir0=Ie0(2);Ir1=Ie1(2);Ir2=Ie2(2); It0=Ie0(3);It1=Ie1(3);It2=Ie2(3); Ism0=Ie0(4);Ism1=Ie1(4);Ism2=Ie2(4); Isp0=Ie0(5);Isp1=Ie1(5);Isp2=Ie2(5); Irm0=Ie0(6);Irm1=Ie1(6);Irm2=Ie2(6); Itp0=Ie0(7);Itp1=Ie1(7);Itp2=Ie2(7); [Ima,Imb,Imc]=abc(Im0,Im1,Im2); [Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2); [Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2); [Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2); [Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2); [Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2); [Vta,Vtb,Vtc]=abc(Vt0,Vt1,Vt2); [Isa,Isb,Isc]=abc(Is0,Is1,Is2); [Ira,Irb,Irc]=abc(Ir0,Ir1,Ir2); [Ita,Itb,Itc]=abc(It0,It1,It2); [Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2); [Ispa,Ispb,Ispc]=abc(Isp0,Isp1,Isp2); [Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2); [Itpa,Itpb,Itpc]=abc(Itp0,Itp1,Itp2); % 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,Itp0); [Tp_ZAB,Tp_ZBC,Tp_ZCA,Tp_ZAG,Tp_ZBG,Tp_ZCG]=... PDistComp(Vta,Vtb,Vtc,Itpa,Itpb,Itpc,Itp0,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); [Tp_ZAB,Tp_ZBC,Tp_ZCA,Tp_ZAG,Tp_ZBG,Tp_ZCG]=... PDist(Vta,Vtb,Vtc,Itpa,Itpb,Itpc,Itp0,k); end % 67N [Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0); VVspol=Vspol*Vbase;IIsmop=Ismop*Ibase; [Vspol,Ispop,AngDirSp]=PDire(Vs0,Isp0); IIspop=Ispop*Ibase; [Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0); VVrpol=Vrpol*Vbase;IIrmop=Irmop*Ibase; [Vtpol,Itpop,AngDirTp]=PDire(Vt0,Itp0); VVtpol=Vtpol*Vbase;IItpop=Itpop*Ibase; alcanceSm_ZAB=abs(Sm_ZAB)/cos(angle(Sm_ZAB)-angle(zl1u)); alcanceSm_ZBC=abs(Sm_ZBC)/cos(angle(Sm_ZBC)-angle(zl1u)); alcanceSm_ZCA=abs(Sm_ZCA)/cos(angle(Sm_ZCA)-angle(zl1u)); alcanceSm_ZAG=abs(Sm_ZAG)/cos(angle(Sm_ZAG)-angle(zl1u)); alcanceSm_ZBG=abs(Sm_ZBG)/cos(angle(Sm_ZBG)-angle(zl1u)); alcanceSm_ZCG=abs(Sm_ZCG)/cos(angle(Sm_ZCG)-angle(zl1u)); alcanceSp_ZAB=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u)); alcanceSp_ZBC=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u)); alcanceSp_ZCA=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u)); alcanceSp_ZAG=abs(Sp_ZAG)/cos(angle(Sp_ZAG)-angle(zl1u)); alcanceSp_ZBG=abs(Sp_ZBG)/cos(angle(Sp_ZBG)-angle(zl1u)); alcanceSp_ZCG=abs(Sp_ZCG)/cos(angle(Sp_ZCG)-angle(zl1u)); alcanceRm_ZAB=abs(Rm_ZAB)/cos(angle(Rm_ZAB)-angle(zl1u)); alcanceRm_ZBC=abs(Rm_ZBC)/cos(angle(Rm_ZBC)-angle(zl1u)); alcanceRm_ZCA=abs(Rm_ZCA)/cos(angle(Rm_ZCA)-angle(zl1u)); alcanceRm_ZAG=abs(Rm_ZAG)/cos(angle(Rm_ZAG)-angle(zl1u)); alcanceRm_ZBG=abs(Rm_ZBG)/cos(angle(Rm_ZBG)-angle(zl1u)); alcanceRm_ZCG=abs(Rm_ZCG)/cos(angle(Rm_ZCG)-angle(zl1u)); alcanceTp_ZAB=abs(Tp_ZAB)/cos(angle(Tp_ZAB)-angle(zl1u)); alcanceTp_ZBC=abs(Tp_ZBC)/cos(angle(Tp_ZBC)-angle(zl1u)); alcanceTp_ZCA=abs(Tp_ZCA)/cos(angle(Tp_ZCA)-angle(zl1u)); alcanceTp_ZAG=abs(Tp_ZAG)/cos(angle(Tp_ZAG)-angle(zl1u)); alcanceTp_ZBG=abs(Tp_ZBG)/cos(angle(Tp_ZBG)-angle(zl1u)); alcanceTp_ZCG=abs(Tp_ZCG)/cos(angle(Tp_ZCG)-angle(zl1u)); if faltam==1 cfaltam='ag'; medidaSm21=Sm_ZAG; medidaRm21=Rm_ZAG; 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 if faltap==1 cfaltap='ag'; alcanceSp=abs(Sp_ZAG)/cos(angle(Sp_ZAG)-angle(zl1u)); alcanceTp=abs(Tp_ZAG)/cos(angle(Tp_ZAG)-angle(zl1u)); elementop='AG'; elseif faltap==2 cfaltap='bg'; alcanceSp=abs(Sp_ZBG)/cos(angle(Sp_ZBG)-angle(zl1u)); alcanceTp=abs(Tp_ZBG)/cos(angle(Tp_ZBG)-angle(zl1u)); elementop='BG'; elseif faltap==3 cfaltap='cg'; alcanceSp=abs(Sp_ZCG)/cos(angle(Sp_ZCG)-angle(zl1u)); alcanceTp=abs(Tp_ZCG)/cos(angle(Tp_ZCG)-angle(zl1u)); elementop='CG'; elseif faltap==4 cfaltap='ab'; alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u)); alcanceTp=abs(Tp_ZAB)/cos(angle(Tp_ZAB)-angle(zl1u)); elementop='AB'; elseif faltap==5 cfaltap='bc'; alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u)); alcanceTp=abs(Tp_ZBC)/cos(angle(Tp_ZBC)-angle(zl1u)); elementop='BC'; elseif faltap==6 cfaltap='ca'; alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u)); alcanceTp=abs(Tp_ZCA)/cos(angle(Tp_ZCA)-angle(zl1u)); elementop='CA'; elseif faltap==7 cfaltap='abg'; alcanceSp=abs(Sp_ZAB)/cos(angle(Sp_ZAB)-angle(zl1u)); alcanceTp=abs(Tp_ZAB)/cos(angle(Tp_ZAB)-angle(zl1u)); elementop='AB'; elseif faltap==8 cfaltap='bcg'; alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u)); alcanceTp=abs(Tp_ZBC)/cos(angle(Tp_ZBC)-angle(zl1u)); elementop='BC'; elseif faltap==9 cfaltap='cag'; alcanceSp=abs(Sp_ZCA)/cos(angle(Sp_ZCA)-angle(zl1u)); alcanceTp=abs(Tp_ZCA)/cos(angle(Tp_ZCA)-angle(zl1u)); elementop='CA'; elseif faltap==10 cfaltap='abc'; alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u)); alcanceTp=abs(Tp_ZBC)/cos(angle(Tp_ZBC)-angle(zl1u)); elementop='BC'; elseif faltap==11 cfaltap='abcg'; alcanceSp=abs(Sp_ZBC)/cos(angle(Sp_ZBC)-angle(zl1u)); alcanceTp=abs(Tp_ZBC)/cos(angle(Tp_ZBC)-angle(zl1u)); elementop='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); alcanceSppu=alcanceSp/abs(zl1u*p*L); alcanceTppu=alcanceTp/abs(zl1u*(1-p)*L); alcanceSppua=alcanceSp/abs(zl1u*L); alcanceTppua=alcanceTp/abs(zl1u*L); Abscisa(1,i)=m; OrdenadaSm(1,i)=alcanceSmpu; OrdenadaRm(1,i)=alcanceRmpu; OrdenadaSma(1,i)=alcanceSmpua; OrdenadaRma(1,i)=alcanceRmpua; OrdenadaSp(1,i)=alcanceSppu; OrdenadaTp(1,i)=alcanceTppu; OrdenadaSpa(1,i)=alcanceSppua; OrdenadaTpa(1,i)=alcanceTppua; Sm21(1,i)=medidaSm21; Rm21(1,i)=medidaRm21; alcanceSm_ZABpua=alcanceSm_ZAB/abs(zl1u*L); alcanceSm_ZBCpua=alcanceSm_ZBC/abs(zl1u*L); alcanceSm_ZCApua=alcanceSm_ZCA/abs(zl1u*L); alcanceSm_ZAGpua=alcanceSm_ZAG/abs(zl1u*L); alcanceSm_ZBGpua=alcanceSm_ZBG/abs(zl1u*L); alcanceSm_ZCGpua=alcanceSm_ZCG/abs(zl1u*L); OrdenadaSm_ZABpua(1,i)=alcanceSm_ZABpua; OrdenadaSm_ZBCpua(1,i)=alcanceSm_ZBCpua; OrdenadaSm_ZCApua(1,i)=alcanceSm_ZCApua; OrdenadaSm_ZAGpua(1,i)=alcanceSm_ZAGpua; OrdenadaSm_ZBGpua(1,i)=alcanceSm_ZBGpua; OrdenadaSm_ZCGpua(1,i)=alcanceSm_ZCGpua; alcanceRm_ZABpua=alcanceRm_ZAB/abs(zl1u*L); alcanceRm_ZBCpua=alcanceRm_ZBC/abs(zl1u*L); alcanceRm_ZCApua=alcanceRm_ZCA/abs(zl1u*L); alcanceRm_ZAGpua=alcanceRm_ZAG/abs(zl1u*L); alcanceRm_ZBGpua=alcanceRm_ZBG/abs(zl1u*L); alcanceRm_ZCGpua=alcanceRm_ZCG/abs(zl1u*L); OrdenadaRm_ZABpua(1,i)=alcanceRm_ZABpua; OrdenadaRm_ZBCpua(1,i)=alcanceRm_ZBCpua; OrdenadaRm_ZCApua(1,i)=alcanceRm_ZCApua; OrdenadaRm_ZAGpua(1,i)=alcanceRm_ZAGpua; OrdenadaRm_ZBGpua(1,i)=alcanceRm_ZBGpua; OrdenadaRm_ZCGpua(1,i)=alcanceRm_ZCGpua; alcanceSp_ZABpua=alcanceSp_ZAB/abs(zl1u*L); alcanceSp_ZBCpua=alcanceSp_ZBC/abs(zl1u*L); alcanceSp_ZCApua=alcanceSp_ZCA/abs(zl1u*L); alcanceSp_ZAGpua=alcanceSp_ZAG/abs(zl1u*L); alcanceSp_ZBGpua=alcanceSp_ZBG/abs(zl1u*L); alcanceSp_ZCGpua=alcanceSp_ZCG/abs(zl1u*L); OrdenadaSp_ZABpua(1,i)=alcanceSp_ZABpua; OrdenadaSp_ZBCpua(1,i)=alcanceSp_ZBCpua; OrdenadaSp_ZCApua(1,i)=alcanceSp_ZCApua; OrdenadaSp_ZAGpua(1,i)=alcanceSp_ZAGpua; OrdenadaSp_ZBGpua(1,i)=alcanceSp_ZBGpua; OrdenadaSp_ZCGpua(1,i)=alcanceSp_ZCGpua; alcanceTp_ZABpua=alcanceTp_ZAB/abs(zl1u*L); alcanceTp_ZBCpua=alcanceTp_ZBC/abs(zl1u*L); alcanceTp_ZCApua=alcanceTp_ZCA/abs(zl1u*L); alcanceTp_ZAGpua=alcanceTp_ZAG/abs(zl1u*L); alcanceTp_ZBGpua=alcanceTp_ZBG/abs(zl1u*L); alcanceTp_ZCGpua=alcanceTp_ZCG/abs(zl1u*L); OrdenadaTp_ZABpua(1,i)=alcanceTp_ZABpua; OrdenadaTp_ZBCpua(1,i)=alcanceTp_ZBCpua; OrdenadaTp_ZCApua(1,i)=alcanceTp_ZCApua; OrdenadaTp_ZAGpua(1,i)=alcanceTp_ZAGpua; OrdenadaTp_ZBGpua(1,i)=alcanceTp_ZBGpua; OrdenadaTp_ZCGpua(1,i)=alcanceTp_ZCGpua; %*************************************** %Para salida 4 %*************************************** vectorVspol(1,i)=VVspol; vectorIsmop(1,i)=IIsmop; vectorAngDirSm(1,i)=AngDirSm; vectorVrpol(1,i)=VVrpol; vectorIrmop(1,i)=IIrmop; vectorAngDirRm(1,i)=AngDirRm; vectorVspol(1,i)=VVspol; vectorIspop(1,i)=IIspop; vectorAngDirSp(1,i)=AngDirSp; vectorVtpol(1,i)=VVtpol; vectorItpop(1,i)=IItpop; vectorAngDirTp(1,i)=AngDirTp; %*************************************** % Para Salida 5 R-X %*************************************** vectorSm_ZAB(1,i)=Sm_ZAB*zbase; vectorSm_ZBC(1,i)=Sm_ZBC*zbase; vectorSm_ZCA(1,i)=Sm_ZCA*zbase; vectorSm_ZAG(1,i)=Sm_ZAG*zbase; vectorSm_ZBG(1,i)=Sm_ZBG*zbase; vectorSm_ZCG(1,i)=Sm_ZCG*zbase; vectorRm_ZAB(1,i)=Rm_ZAB*zbase; vectorRm_ZBC(1,i)=Rm_ZBC*zbase; vectorRm_ZCA(1,i)=Rm_ZCA*zbase; vectorRm_ZAG(1,i)=Rm_ZAG*zbase; vectorRm_ZBG(1,i)=Rm_ZBG*zbase; vectorRm_ZCG(1,i)=Rm_ZCG*zbase; vectorSp_ZAB(1,i)=Sp_ZAB*zbase; vectorSp_ZBC(1,i)=Sp_ZBC*zbase; vectorSp_ZCA(1,i)=Sp_ZCA*zbase; vectorSp_ZAG(1,i)=Sp_ZAG*zbase; vectorSp_ZBG(1,i)=Sp_ZBG*zbase; vectorSp_ZCG(1,i)=Sp_ZCG*zbase; vectorTp_ZAB(1,i)=Tp_ZAB*zbase; vectorTp_ZBC(1,i)=Tp_ZBC*zbase; vectorTp_ZCA(1,i)=Tp_ZCA*zbase; vectorTp_ZAG(1,i)=Tp_ZAG*zbase; vectorTp_ZBG(1,i)=Tp_ZBG*zbase; vectorTp_ZCG(1,i)=Tp_ZCG*zbase; %*************************************** i=i+1; m=m+0.01; end % ----------------------------------------------------------------------- if compens21==1 comp='21c'; elseif compens21==2 comp='21nc'; end % ----------------------------------------------------------------------- % Para Salidas 1 y 2 cadena0=sprintf('PG34ic,%s,%4.0f(KV),kE=%-4.2f,%4.0f(km)',comp,Un,kE,L); cadena1=sprintf('Sccs=%-5.0f,Sccr=%-5.0f,Scct=%-5.0f\n',Sccs,Sccr,Scct); cadena2=sprintf('Ks0=%-4.2f,Kr0=%-4.2f,k0=%-4.1f',Ks0,Kr0,k0); cadena3=sprintf('SR,m:%s,%s,Rfm=%-3.0f',cfaltam,elementom,Rfm); cadenam=sprintf('%s,%s,%s,%s',cadena0,cadena1,cadena2,cadena3); if faltam==7|faltam==8|faltam==9|faltam==11 cadena4=sprintf('Rgm=%-3.0f',Rgm); cadenam=sprintf('%s,%s',cadenam,cadena4); end cadenap=sprintf('ST,p:%s,%s,Rfp=%-3.0f',cfaltap,elementop,Rfp); if faltap==7|faltap==8|faltap==9|faltap==11 cadenap=sprintf('ST,p:%s,%s,Rfp=%-3.0f,Rgp=%-3.0f',cfaltap,elementop,Rfp,Rgp); end cadenamp=sprintf('%s,%s',cadenam,cadenap); % Para Salida 3 y 5 cadena0c=sprintf('PG34ic,m=p,%s,%4.0f(KV),kE=%-4.2f,%4.0f(km)',comp,Un,kE,L); cadena1c=sprintf('Sccs=%-5.0f,Sccr=%-5.0f,Scct=%-5.0f\n',Sccs,Sccr,Scct); cadena2c=sprintf('Ks0=%-4.2f,Kr0=%-4.2f,k0=%-4.1f',Ks0,Kr0,k0); cadena3c=sprintf('SR,m:%s,Rfm=%-3.0f',cfaltam,Rfm); cadenamc=sprintf('%s,%s,%s,%s',cadena0c,cadena1c,cadena2c,cadena3c); if faltam==7|faltam==8|faltam==9|faltam==11 cadena4c=sprintf('Rgm=%-3.0f',Rgm); cadenamc=sprintf('%s,%s',cadenamc,cadena4c); end cadenapc=sprintf('ST,p:%s,Rfp=%-3.0f',cfaltap,Rfp); if faltap==7|faltap==8|faltap==9|faltap==11 cadenapc=sprintf('ST,p:%s,Rfp=%-3.0f, Rgp=%-3.0f',cfaltap,Rfp,Rgp); end cadenampc=sprintf('%s,%s',cadenamc,cadenapc); % Para Salida 4 dircadena0c=sprintf('PG34ic,67N,%4.0f(KV),kE=%-4.2f,%4.0f(km)',Un,kE,L); dircadena1c=sprintf('Sccs=%-5.0f,Sccr=%-5.0f,Scct=%-5.0f\n',Sccs,Sccr,Scct); dircadena2c=sprintf('Ks0=%-4.2f,Kr0=%-4.2f,k0=%-4.1f',Ks0,Kr0,k0); dircadena3c=sprintf('SR,m:%s,Rfm=%-3.0f',cfaltam,Rfm); dircadenamc=sprintf('%s,%s,%s,%s',dircadena0c,dircadena1c,dircadena2c,dircadena3c); if faltam==7|faltam==8|faltam==9|faltam==11 dircadena4c=sprintf('Rgm=%-3.0f',Rgm); dircadenamc=sprintf('%s,%s',dircadenamc,dircadena4c); end dircadenapc=sprintf('ST,p:%s,Rfp=%-3.0f',cfaltap,Rfp); if faltap==7|faltap==8|faltap==9|faltap==11 dircadenapc=sprintf('ST,p:%s,Rfp=%-3.0f, Rgp=%-3.0f',cfaltap,Rfp,Rgp); end dircadenampc=sprintf('%s,%s',dircadenamc,dircadenapc); % ----------------------------------------------------------------------- variassalidas=1; while variassalidas==1 fprintf('\n'); fprintf('Salida Grafica 1: Sm,Rm(2 Subpl),Unidad Pref.(pu)\n'); fprintf(' Sp,Rp(2 Subpl),Unidad Pref.(pu)\n'); fprintf('Salida Grafica 1: Sm,Rm(2 Subpl),Unidad Pref.(pua)\n'); fprintf(' Sp,Rp(2 Subpl),Unidad Pref.(pua)\n'); fprintf('Salida Grafica 3: Analisis Sm, Todas las Unidades (pua)\n'); fprintf(' : Analisis Rm, Todas las Unidades (pua)\n'); fprintf(' : Analisis Sp, Todas las Unidades (pua)\n'); fprintf(' : Analisis Tp, Todas las Unidades (pua)\n'); fprintf('Salida Grafica 4: Analisis Sm, 67N\n'); fprintf(' : Analisis Rm, 67N\n'); fprintf(' : Analisis Sp, 67N\n'); fprintf(' : Analisis Tp, 67N\n'); fprintf('Salida Grafica 5: Analisis Sm, R-X\n'); fprintf(' : Analisis Rm, R-X\n'); fprintf(' : Analisis Sp, R-X\n'); fprintf(' : Analisis Tp, R-X\n'); fprintf('\n'); salida=input('Salida: 1,2,3,4,5 :'); switch(salida) case{1}, % Salida Grafica 1 (pu) Abscisa1m=Abscisa(1:1,1:(i-1)); Ordenada1m=OrdenadaSm(1:1,1:(i-1)); Abscisa4m=Abscisa(1:1,1:(i-1)); Ordenada4m=OrdenadaRm(1:1,1:(i-1)); Abscisa1p=Abscisa(1:1,1:(i-1)); Ordenada1p=OrdenadaSp(1:1,1:(i-1)); Abscisa4p=Abscisa(1:1,1:(i-1)); Ordenada4p=OrdenadaTp(1:1,1:(i-1)); 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 AnalisisLinea=input('Extremos m=1, p=2: '); switch(AnalisisLinea) case{1}, % salida 1, Extremos m subplot(2,1,1) plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'g',Abscisa1m,Ordenada1m,'r'),grid on axis([0 1 0 3]); title(cadenamp) xlabel('m=p'),ylabel('alcanceSmpu') legend('Ref','Lim') subplot(2,1,2) plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadd,'g',Abscisa4m,Ordenada4m,'r'),grid on axis([0 1 0 3]); %title(cadenamp) xlabel('m=p'),ylabel('alcanceRmpu') legend('Ref','Lim') case{2}, % salida 1, Extremos p subplot(2,1,1) plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadi,'g',Abscisa1p,Ordenada1p,'r'),grid on axis([0 1 0 3]); title(cadenamp) xlabel('m=p'),ylabel('alcanceSppu') legend('Ref','Lim') subplot(2,1,2) plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadadd,'g',Abscisa4p,Ordenada4p,'r'),grid on axis([0 1 0 3]); %title(cadenamp) xlabel('m=p'),ylabel('alcanceTppu') legend('Ref','Lim') end case(2) % Salida Grafica 2 (pua) 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)); Abscisa4ma=Abscisa(1:1,1:(i-1)); Ordenada4ma=OrdenadaRma(1:1,1:(i-1)); Abscisa1pa=Abscisa(1:1,1:(i-1)); Ordenada1pa=OrdenadaSpa(1:1,1:(i-1)); Abscisa4pa=Abscisa(1:1,1:(i-1)); Ordenada4pa=OrdenadaTpa(1:1,1:(i-1)); AnalisisLinea=input('Extremos m=1, p=2: '); switch(AnalisisLinea) case{1}, % salida 2, Extremos m subplot(2,1,1) plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadaz1,'m',Abscisa1ma,Ordenada1ma,'r'),grid on axis([0 1 0 2]); title(cadenamp) xlabel('m=p'),ylabel('alcanceSmpua') legend('45','z1') subplot(2,1,2) plot(Abscisa00,Ordenada00,'c',Abscisa0,Ordenadaz1,'m',Abscisa4ma,Ordenada4ma,'r'),grid on axis([0 1 0 2]); %title(cadenamp) xlabel('m=p'),ylabel('alcanceRmpua') legend('-45','z1') case{2}, % salida 2, Extremos p subplot(2,1,1) plot(Abscisa0,Ordenada0,'c',Abscisa0,Ordenadaz1,'m',Abscisa1pa,Ordenada1pa,'r'),grid on axis([0 1 0 2]); title(cadenamp) xlabel('m=p'),ylabel('alcanceSppua') legend('45','z1') subplot(2,1,2) plot(Abscisa00,Ordenada00,'c',Abscisa0,Ordenadaz1,'m',Abscisa4pa,Ordenada4pa,'r'),grid on axis([0 1 0 2]); %title(cadenamp) xlabel('m=p'),ylabel('alcanceTppua') legend('-45','z1') end case{3} % Salida Grafica 3 (Todas las Unidades pua) mm=0.01; cont=1; while mm<1; Abscisa0(cont)=mm; Ordenada0(cont)=0.8; Ordenada0a(cont)=mm; Ordenada00a(cont)=1-mm; cont=cont+1; mm=mm+0.01; end extremos=input('Extremos: m=1, p=2 :'); switch(extremos) case{1}, extr=input('Extremo: Sm=1, Rm=2 :'); switch(extr) case{1}, cadylabel1a='alcSmpua'; % Colores b=plot(Abscisa0,Ordenada0a); set(b,'Color','c') % Fin Colores hold on plot(Abscisa,OrdenadaSm_ZABpua,'r',Abscisa,OrdenadaSm_ZBCpua,'g',Abscisa,OrdenadaSm_ZCApua,'b',... Abscisa,OrdenadaSm_ZAGpua,'m',Abscisa,OrdenadaSm_ZBGpua,'k'),grid on hold off % Colores hold on a=plot(Abscisa,OrdenadaSm_ZCGpua); set(a,'Color',[1 0.5490 0]) h=plot(Abscisa0,Ordenada0); set(h,'Color',[0.7216 0.4510 0.20]) % Fin Colores axis([0 1 -2 2]); title(cadenampc) xlabel('m=p'),ylabel(cadylabel1a) legend('45','AB','BC','CA','AG','BG','CG','z1') case{2}, cadylabel1a='alcRmpua'; % Colores b=plot(Abscisa0,Ordenada00a); set(b,'Color','c') % Fin Colores hold on plot(Abscisa,OrdenadaRm_ZABpua,'r',Abscisa,OrdenadaRm_ZBCpua,'g',Abscisa,OrdenadaRm_ZCApua,'b',... Abscisa,OrdenadaRm_ZAGpua,'m',Abscisa,OrdenadaRm_ZBGpua,'k'),grid on hold off % Colores hold on a=plot(Abscisa,OrdenadaRm_ZCGpua); 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(cadenampc) xlabel('m=p'),ylabel(cadylabel1a) legend('-45','AB','BC','CA','AG','BG','CG','z1') end case{2}, extr=input('Extremo: Sp=1, Tp=2 :'); switch(extr) case{1}, cadylabel1a='alcSppua'; % Colores b=plot(Abscisa0,Ordenada0a); set(b,'Color','c') % Fin Colores hold on plot(Abscisa,OrdenadaSp_ZABpua,'r',Abscisa,OrdenadaSp_ZBCpua,'g',Abscisa,OrdenadaSp_ZCApua,'b',... Abscisa,OrdenadaSp_ZAGpua,'m',Abscisa,OrdenadaSp_ZBGpua,'k'),grid on hold off % Colores hold on a=plot(Abscisa,OrdenadaSp_ZCGpua); 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(cadenampc) xlabel('m=p'),ylabel(cadylabel1a) legend('45','AB','BC','CA','AG','BG','CG','z1') case{2}, cadylabel1a='alcTppua'; % Colores b=plot(Abscisa0,Ordenada00a); set(b,'Color','c') % Fin Colores hold on plot(Abscisa,OrdenadaTp_ZABpua,'r',Abscisa,OrdenadaTp_ZBCpua,'g',Abscisa,OrdenadaTp_ZCApua,'b',... Abscisa,OrdenadaTp_ZAGpua,'m',Abscisa,OrdenadaTp_ZBGpua,'k'),grid on hold off % Colores hold on a=plot(Abscisa,OrdenadaTp_ZCGpua); 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(cadenampc) xlabel('m=p'),ylabel(cadylabel1a) legend('-45','AB','BC','CA','AG','BG','CG','z1') end end %extremo case{4} % salida 4 (Direccional) extr=input('Extremo: Sm=1, Rm=2, Sp=3, Tp=4: '); switch(extr) case{1} subplot(3,1,1) plot(Abscisa,vectorAngDirSm,'b'); axis([0 1 -180 180]); xlabel('m=p'),ylabel('AngDirSm(grados)') title(dircadenampc) rayar(1); legend('AngDirSm','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVspol,'b'); grid %title(dircadenampc) xlabel('m=p'),ylabel('modulo Vspol(kV)') subplot(3,1,3) plot(Abscisa,vectorIsmop,'b'); grid xlabel('m=p'),ylabel('modulo Ismop(kA)') case{2} subplot(3,1,1) plot(Abscisa,vectorAngDirRm,'b'); axis([0 1 -180 180]); xlabel('m=p'),ylabel('AngDirRm(grados)') title(dircadenampc) rayar(1); legend('AngDirRm','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVrpol,'b'); grid %title(dircadenampc) xlabel('m=p'),ylabel('modulo Vrpol(kV)') subplot(3,1,3) plot(Abscisa,vectorIrmop,'b'); grid xlabel('m=p'),ylabel('modulo Irmop(kA)') case{3} subplot(3,1,1) plot(Abscisa,vectorAngDirSp,'b'); axis([0 1 -180 180]); xlabel('m=p'),ylabel('AngDirSp(grados)') title(dircadenampc) rayar(1); legend('AngDirSp','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVspol,'b'); grid %title(dircadenampc) xlabel('m=p'),ylabel('modulo Vspol(kV)') subplot(3,1,3) plot(Abscisa,vectorIspop,'b'); grid xlabel('m=p'),ylabel('modulo Ispop(kA)') case{4} subplot(3,1,1) plot(Abscisa,vectorAngDirTp,'b'); axis([0 1 -180 180]); xlabel('m=p'),ylabel('AngDirTp(grados)') title(dircadenampc) rayar(1); legend('AngDirTp','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVtpol,'b'); grid %title(dircadenampc) xlabel('m=p'),ylabel('modulo Vtpol(kV)') subplot(3,1,3) plot(Abscisa,vectorItpop,'b'); grid xlabel('m=p'),ylabel('modulo Itpop(kA)') end % end del switch extr case{5}, % ********* Salida 5 ************ extr=input('Extremo: Sm=1, Rm=2, Sp=3, Tp=4: '); switch(extr) case{1} proteccion='Sm'; pzona=0.8; hold on plot19a(L,zl1u*zbase,pzona,proteccion) hold off hold on % Colores a=plot(vectorSm_ZAB); set(a,'Color',[1 0.60 0]) % Fin Colores plot(vectorSm_ZBC,'m') plot(vectorSm_ZCA,'c') plot(vectorSm_ZAG,'r') plot(vectorSm_ZBG,'k') % Colores c=plot(vectorSm_ZCG); set(c,'Color',[0.7216 0.4510 0.20]) % Fin Colores hold off [subAbscisa,subvectorSm_ZAB]=subconjunto21a(Abscisa,vectorSm_ZAB); subAbscisa(:,1)=[]; subvectorSm_ZAB(:,1)=[]; hold on b=plot(subvectorSm_ZAB,'*'); set(b,'Color',[1 0.60 0]) plot20(L,zl1u*zbase,subvectorSm_ZAB,'naranja') hold off [subAbscisa,subvectorSm_ZBC]=subconjunto21a(Abscisa,vectorSm_ZBC); subAbscisa(:,1)=[]; subvectorSm_ZBC(:,1)=[]; hold on plot(subvectorSm_ZBC,'m*'); plot20(L,zl1u*zbase,subvectorSm_ZBC,'m') hold off [subAbscisa,subvectorSm_ZCA]=subconjunto21a(Abscisa,vectorSm_ZCA); subAbscisa(:,1)=[]; subvectorSm_ZCA(:,1)=[]; hold on plot(subvectorSm_ZCA,'c*'); plot20(L,zl1u*zbase,subvectorSm_ZCA,'c') hold off [subAbscisa,subvectorSm_ZAG]=subconjunto21a(Abscisa,vectorSm_ZAG); subAbscisa(:,1)=[]; subvectorSm_ZAG(:,1)=[]; hold on plot(subvectorSm_ZAG,'r*'); plot20(L,zl1u*zbase,subvectorSm_ZAG,'r') hold off [subAbscisa,subvectorSm_ZBG]=subconjunto21a(Abscisa,vectorSm_ZBG); subAbscisa(:,1)=[]; subvectorSm_ZBG(:,1)=[]; hold on plot(subvectorSm_ZBG,'k*'); plot20(L,zl1u*zbase,subvectorSm_ZBG,'k') hold off [subAbscisa,subvectorSm_ZCG]=subconjunto21a(Abscisa,vectorSm_ZCG); subAbscisa(:,1)=[]; subvectorSm_ZCG(:,1)=[]; hold on % Colores c=plot(subvectorSm_ZCG,'*'); set(c,'Color',[0.7216 0.4510 0.20]) % Fin Colores plot20(L,zl1u*zbase,subvectorSm_ZCG,'marron') hold off RXcadenap=sprintf('Proteccion: %s',proteccion); RXcadenampc=sprintf('%s,%s',cadenampc,RXcadenap); title(RXcadenampc) legend('Mho','Linea','m=0..1','AB','BC','CA','AG','BG','CG') case{2} proteccion='Rm'; pzona=0.8; hold on plot19a(L,zl1u*zbase,pzona,proteccion) hold off hold on % Colores a=plot(vectorRm_ZAB); set(a,'Color',[1 0.60 0]) % Fin Colores plot(vectorRm_ZBC,'m') plot(vectorRm_ZCA,'c') plot(vectorRm_ZAG,'r') plot(vectorRm_ZBG,'k') % Colores c=plot(vectorRm_ZCG); set(c,'Color',[0.7216 0.4510 0.20]) % Fin Colores hold off [subAbscisa,subvectorRm_ZAB]=subconjunto21a(Abscisa,vectorRm_ZAB); subAbscisa(:,11)=[]; subvectorRm_ZAB(:,11)=[]; invsubvectorRm_ZAB=subvectorRm_ZAB(1:1,10:-1:1); hold on b=plot(invsubvectorRm_ZAB,'*'); set(b,'Color',[1 0.60 0]) plot20(L,zl1u*zbase,invsubvectorRm_ZAB,'naranja') hold off [subAbscisa,subvectorRm_ZBC]=subconjunto21a(Abscisa,vectorRm_ZBC); subAbscisa(:,11)=[]; subvectorRm_ZBC(:,11)=[]; invsubvectorRm_ZBC=subvectorRm_ZBC(1:1,10:-1:1); hold on plot(subvectorRm_ZBC,'m*'); plot20(L,zl1u*zbase,invsubvectorRm_ZBC,'m') hold off [subAbscisa,subvectorRm_ZCA]=subconjunto21a(Abscisa,vectorRm_ZCA); subAbscisa(:,11)=[]; subvectorRm_ZCA(:,11)=[]; invsubvectorRm_ZCA=subvectorRm_ZCA(1:1,10:-1:1); hold on plot(subvectorRm_ZCA,'c*'); plot20(L,zl1u*zbase,invsubvectorRm_ZCA,'c') hold off [subAbscisa,subvectorRm_ZAG]=subconjunto21a(Abscisa,vectorRm_ZAG); subAbscisa(:,11)=[]; subvectorRm_ZAG(:,11)=[]; invsubvectorRm_ZAG=subvectorRm_ZAG(1:1,10:-1:1); hold on plot(subvectorRm_ZAG,'r*'); plot20(L,zl1u*zbase,invsubvectorRm_ZAG,'r') hold off [subAbscisa,subvectorRm_ZBG]=subconjunto21a(Abscisa,vectorRm_ZBG); subAbscisa(:,11)=[]; subvectorRm_ZBG(:,11)=[]; invsubvectorRm_ZBG=subvectorRm_ZBG(1:1,10:-1:1); hold on plot(subvectorRm_ZBG,'k*'); plot20(L,zl1u*zbase,invsubvectorRm_ZBG,'k') hold off [subAbscisa,subvectorRm_ZCG]=subconjunto21a(Abscisa,vectorRm_ZCG); subAbscisa(:,11)=[]; subvectorRm_ZCG(:,11)=[]; invsubvectorRm_ZCG=subvectorRm_ZCG(1:1,10:-1:1); hold on % Colores c=plot(invsubvectorRm_ZCG,'*'); set(c,'Color',[0.7216 0.4510 0.20]) % Fin Colores plot20(L,zl1u*zbase,invsubvectorRm_ZCG,'marron') hold off RXcadenap=sprintf('Proteccion: %s',proteccion); RXcadenampc=sprintf('%s,%s',cadenampc,RXcadenap); title(RXcadenampc) legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG') case{3} proteccion='Sp'; pzona=0.8; hold on plot19a(L,zl1u*zbase,pzona,proteccion) hold off hold on % Colores a=plot(vectorSp_ZAB); set(a,'Color',[1 0.60 0]) % Fin Colores plot(vectorSp_ZBC,'m') plot(vectorSp_ZCA,'c') plot(vectorSp_ZAG,'r') plot(vectorSp_ZBG,'k') % Colores c=plot(vectorSp_ZCG); set(c,'Color',[0.7216 0.4510 0.20]) % Fin Colores hold off [subAbscisa,subvectorSp_ZAB]=subconjunto21a(Abscisa,vectorSp_ZAB); subAbscisa(:,1)=[]; subvectorSp_ZAB(:,1)=[]; hold on b=plot(subvectorSp_ZAB,'*'); set(b,'Color',[1 0.60 0]) plot20(L,zl1u*zbase,subvectorSp_ZAB,'naranja') hold off [subAbscisa,subvectorSp_ZBC]=subconjunto21a(Abscisa,vectorSp_ZBC); subAbscisa(:,1)=[]; subvectorSp_ZBC(:,1)=[]; hold on plot(subvectorSp_ZBC,'m*'); plot20(L,zl1u*zbase,subvectorSp_ZBC,'m') hold off [subAbscisa,subvectorSp_ZCA]=subconjunto21a(Abscisa,vectorSp_ZCA); subAbscisa(:,1)=[]; subvectorSp_ZCA(:,1)=[]; hold on plot(subvectorSp_ZCA,'c*'); plot20(L,zl1u*zbase,subvectorSp_ZCA,'c') hold off [subAbscisa,subvectorSp_ZAG]=subconjunto21a(Abscisa,vectorSp_ZAG); subAbscisa(:,1)=[]; subvectorSp_ZAG(:,1)=[]; hold on plot(subvectorSp_ZAG,'r*'); plot20(L,zl1u*zbase,subvectorSp_ZAG,'r') hold off [subAbscisa,subvectorSp_ZBG]=subconjunto21a(Abscisa,vectorSp_ZBG); subAbscisa(:,1)=[]; subvectorSp_ZBG(:,1)=[]; hold on plot(subvectorSp_ZBG,'k*'); plot20(L,zl1u*zbase,subvectorSp_ZBG,'k') hold off [subAbscisa,subvectorSp_ZCG]=subconjunto21a(Abscisa,vectorSp_ZCG); subAbscisa(:,1)=[]; subvectorSp_ZCG(:,1)=[]; hold on % Colores c=plot(subvectorSp_ZCG,'*'); set(c,'Color',[0.7216 0.4510 0.20]) % Fin Colores plot20(L,zl1u*zbase,subvectorSp_ZCG,'marron') hold off RXcadenap=sprintf('Proteccion: %s',proteccion); RXcadenampc=sprintf('%s,%s',cadenampc,RXcadenap); title(RXcadenampc) legend('Mho','Linea','m=0..1','AB','BC','CA','AG','BG','CG') case{4} proteccion='Tp'; pzona=0.8; hold on plot19a(L,zl1u*zbase,pzona,proteccion) hold off hold on % Colores a=plot(vectorTp_ZAB); set(a,'Color',[1 0.60 0]) % Fin Colores plot(vectorTp_ZBC,'m') plot(vectorTp_ZCA,'c') plot(vectorTp_ZAG,'r') plot(vectorTp_ZBG,'k') % Colores c=plot(vectorTp_ZCG); set(c,'Color',[0.7216 0.4510 0.20]) % Fin Colores hold off [subAbscisa,subvectorTp_ZAB]=subconjunto21a(Abscisa,vectorTp_ZAB); subAbscisa(:,11)=[]; subvectorTp_ZAB(:,11)=[]; invsubvectorTp_ZAB=subvectorTp_ZAB(1:1,10:-1:1); hold on b=plot(invsubvectorTp_ZAB,'*'); set(b,'Color',[1 0.60 0]) plot20(L,zl1u*zbase,invsubvectorTp_ZAB,'naranja') hold off [subAbscisa,subvectorTp_ZBC]=subconjunto21a(Abscisa,vectorTp_ZBC); subAbscisa(:,11)=[]; subvectorTp_ZBC(:,11)=[]; invsubvectorTp_ZBC=subvectorTp_ZBC(1:1,10:-1:1); hold on plot(subvectorTp_ZBC,'m*'); plot20(L,zl1u*zbase,invsubvectorTp_ZBC,'m') hold off [subAbscisa,subvectorTp_ZCA]=subconjunto21a(Abscisa,vectorTp_ZCA); subAbscisa(:,11)=[]; subvectorTp_ZCA(:,11)=[]; invsubvectorTp_ZCA=subvectorTp_ZCA(1:1,10:-1:1); hold on plot(subvectorTp_ZCA,'c*'); plot20(L,zl1u*zbase,invsubvectorTp_ZCA,'c') hold off [subAbscisa,subvectorTp_ZAG]=subconjunto21a(Abscisa,vectorTp_ZAG); subAbscisa(:,11)=[]; subvectorTp_ZAG(:,11)=[]; invsubvectorTp_ZAG=subvectorTp_ZAG(1:1,10:-1:1); hold on plot(subvectorTp_ZAG,'r*'); plot20(L,zl1u*zbase,invsubvectorTp_ZAG,'r') hold off [subAbscisa,subvectorTp_ZBG]=subconjunto21a(Abscisa,vectorTp_ZBG); subAbscisa(:,11)=[]; subvectorTp_ZBG(:,11)=[]; invsubvectorTp_ZBG=subvectorTp_ZBG(1:1,10:-1:1); hold on plot(subvectorTp_ZBG,'k*'); plot20(L,zl1u*zbase,invsubvectorTp_ZBG,'k') hold off [subAbscisa,subvectorTp_ZCG]=subconjunto21a(Abscisa,vectorTp_ZCG); subAbscisa(:,11)=[]; subvectorTp_ZCG(:,11)=[]; invsubvectorTp_ZCG=subvectorTp_ZCG(1:1,10:-1:1); hold on % Colores c=plot(invsubvectorTp_ZCG,'*'); set(c,'Color',[0.7216 0.4510 0.20]) % Fin Colores plot20(L,zl1u*zbase,invsubvectorTp_ZCG,'marron') hold off RXcadenap=sprintf('Proteccion: %s',proteccion); RXcadenampc=sprintf('%s,%s',cadenampc,RXcadenap); title(RXcadenampc) legend('Mho','Linea','1-m=0..1','AB','BC','CA','AG','BG','CG') end % end del switch extr de la salida 5 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.


ANÁLISIS DE ACTUACIÓN DE LAS PROTECCIONES 21 Y 67N ANTE FALTAS SIMPLES, MULTIPLES E INTERCIRCUITO EN LÍNEAS ELÉCTRICAS SOBRE LOS MISMOS APOYOS

: Ruiz Lozano, José Luis
: Ingeniería Industrial
Contenido del proyecto:
Directorio raíz  >  z_Prog  >  PGc  >  PG34ic  >  FG34ic.m