% ----------------------------------------------------------------------- clc variasejecuciones=1; while variasejecuciones==1 clear % ----------------------------------------------------------------------- fprintf('\n'); fprintf('PG52c\n'); fprintf('Esquema Linea Doble SR, Bifurcando a X,Y\n'); fprintf('5 Tramos\n'); fprintf('Unominal kV: 400,220,132,66\n'); fprintf('Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n'); fprintf('Eleccion: kE,L1,L2,L3,L4,L5,Sccs,Sccx,Sccr,Sccy\n'); fprintf(' Kx0,Kr0,Ky0,k0=zs0/zr0\n'); fprintf('Compensacion 21: Si,No\n'); fprintf('Falta Simple. Para cada Tramo: m=0...1\n'); fprintf('Barra m a elegir entre a,b,c,d,e,f,g,h,q,r\n'); fprintf('Analisis: Extremos de la Linea en Falta \n'); fprintf('Salida Grafica 1: Extremos m, Todas las Unidades (pu)\n'); fprintf('Salida Grafica 2: Extremos m, Todas las Unidades (pua)\n'); fprintf('Salida Grafica 3: 67N, Todas las Posiciones\n'); fprintf('\n'); % ---------------------------------------------- %PG52c datos=input('datos: standard=1,input=2: '); if datos==1 Un=220; kE=1; L1=100; L2=100; L3=100; L4=100; L5=100; Sccs=5000; Sccx=5000; Sccr=5000; Sccy=5000; Kx0=1; Kr0=1; Ky0=1; compens21=2; else % ---------------------------------------------- Un=input('U(Fase-Fase) kV:400,220,132,66: '); kE=input('kE=E/Enominal: '); L1=input('Longitud Tramo 1 km: '); L2=input('Longitud Tramo 2 km: '); L3=input('Longitud Tramo 3 km: '); L4=input('Longitud Tramo 4 km: '); L5=input('Longitud Tramo 5 km: '); Sccs=input('Sccs(MVA): '); Sccx=input('Sccx(MVA): '); Sccr=input('Sccr(MVA): '); Sccy=input('Sccy(MVA): '); Kx0=input('Kx0: '); Kr0=input('Kr0: '); Ky0=input('Ky0: '); compens21=input('Compensacion 21: Si=1, No=2 :'); end %---------------------------------------------------------------- k0=input('k0=zs0/zr0: '); fprintf('Punto en Falta: T1sup=11,T4izq=12,T4der=13,T2sup=14,T3sup=15\n'); fprintf(' : T3inf=16,T5der=17,T5izq=18,T2inf=19,T1inf=20\n'); pfaltam=input('Punto en Falta: '); faltam=input('Falta en m: ag=1,bg=2,cg=3,ab=4,bc=5,ca=6,abg=7,bcg=8,cag=9,abc=10,abcg=11 :'); Rfm=input('Resistencia de Falta en m. Rf(Ohm) = '); if faltam==7|faltam==8|faltam==9|faltam==11 Rgm=input('Resistencia de Falta en m. Rg(Ohm) = '); else Rgm=0; end %---------------------------------------------------------------- En=Un/sqrt(3); U=kE*Un; E=kE*En; if Un==400 zl1u=0.0346+0.3208i; zl2u=zl1u; zl0u=0.1828+1.1482i; zlM0u=0.1482+0.7383i; elseif Un==220 zl1u=0.0829+0.4135i; zl2u=zl1u; zl0u=0.2311+1.3085i; zlM0u=0.1482+0.8031i; elseif Un==132 zl1u=0.1342+0.4042i; zl2u=zl1u; zl0u=0.2824+1.3718i; zlM0u=0.1482+0.8947i; elseif Un==66 zl1u=0.2156+0.3952i; zl2u=zl1u; zl0u=0.3636+1.4278i; zlM0u=0.1482+0.9288i; end k=(zl0u-zl1u)/(3*zl1u); zs1=(Un^2/Sccs)*i; zs2=zs1; zx1=(Un^2/Sccx)*i; zx2=zx1; zr1=(Un^2/Sccr)*i; zr2=zr1; zy1=(Un^2/Sccy)*i; zy2=zy1; zx0=Kx0*zx1; zr0=Kr0*zr1; zy0=Ky0*zy1; zs0=k0*zr0; Km=zlM0u/zl1u; fprintf('\n'); % --------------------------------------------------------------- a=0.5;b=0.5;c=0.5;d=0.5;e=0.5;f=0.5;g=0.5;h=0.5;q=0.5;r=0.5; % --------------------------------------------------------------- [Minc]=subMinc(1); i=1; m=0.01; while m<1; if pfaltam==11 a=m;elseif pfaltam==12 b=m;elseif pfaltam==13 c=m;elseif pfaltam==14 d=m; elseif pfaltam==15 e=m;elseif pfaltam==16 f=m;elseif pfaltam==17 g=m; elseif pfaltam==18 h=m;elseif pfaltam==19 q=m;elseif pfaltam==20 r=m; end [z0]=subz0(zs0,zx0,zr0,zy0,zl0u,zlM0u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r); [z1]=subz1(zs1,zx1,zr1,zy1,zl1u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r); [z2]=subz2(zs2,zx2,zr2,zy2,zl2u,L1,L2,L3,L4,L5,a,b,c,d,e,f,g,h,q,r); % MATRIZ DE ADMITANCIA DE BARRAS [Y0]=subY0(Minc,z0); [Y1]=subY1(Minc,z1); [Y2]=subY2(Minc,z2); % ZBUS(0,1,2) Z0=inv(Y0); Z0(11,11); Z1=inv(Y1); Z1(11,11); Z2=inv(Y2); Z2(11,11); % CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p DE FALTA [Im0,Im1,Im2,Vm0,Vm1,Vm2]=mag2m(Z0,Z1,Z2,E,faltam,Rfm,Rgm,pfaltam); % IBUS(0,1,2) I0=zeros(20,1); I1=zeros(20,1); I2=zeros(20,1); I0(pfaltam,1)=-Im0; I1(pfaltam,1)=-Im1; I2(pfaltam,1)=-Im2; % VBUS(0,1,2) V0=Z0*I0; V1=Z1*I1+E; V2=Z2*I2; % I(0,1,2) DE LOS ELEMENTOS [Ie0]=subIe0(Minc,V0,z0); [Ie1]=subIe1(Minc,V1,z1); [Ie2]=subIe2(Minc,V2,z2); % CALCULO DE LAS COMPONENTES DE FASE DE V,I EN EL PUNTO m DE FALTA [Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2); [Ima,Imb,Imc]=abc(Im0,Im1,Im2); % % CALCULO DE LAS COMPONENTES DE FASE DE V,I EN LOS EXTREMOS DE LINEA Vs0=V0(1);Vs1=V1(1);Vs2=V2(1); Vx0=V0(2);Vx1=V1(2);Vx2=V2(2); Vr0=V0(3);Vr1=V1(3);Vr2=V2(3); Vy0=V0(4);Vy1=V1(4);Vy2=V2(4); [Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2); [Vxa,Vxb,Vxc]=abc(Vx0,Vx1,Vx2); [Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2); [Vya,Vyb,Vyc]=abc(Vy0,Vy1,Vy2); Isx0=Ie0(5);Isx1=Ie1(5);Isx2=Ie2(5); Isy0=Ie0(6);Isy1=Ie1(6);Isy2=Ie2(6); Ixs0=Ie0(7);Ixs1=Ie1(7);Ixs2=Ie2(7); Ixr0=Ie0(8);Ixr1=Ie1(8);Ixr2=Ie2(8); Irx0=Ie0(9);Irx1=Ie1(9);Irx2=Ie2(9); Iry0=Ie0(10);Iry1=Ie1(10);Iry2=Ie2(10); Iys0=Ie0(11);Iys1=Ie1(11);Iys2=Ie2(11); Iyr0=Ie0(12);Iyr1=Ie1(12);Iyr2=Ie2(12); [Isxa,Isxb,Isxc]=abc(Isx0,Isx1,Isx2); [Isya,Isyb,Isyc]=abc(Isy0,Isy1,Isy2); [Ixsa,Ixsb,Ixsc]=abc(Ixs0,Ixs1,Ixs2); [Ixra,Ixrb,Ixrc]=abc(Ixr0,Ixr1,Ixr2); [Irxa,Irxb,Irxc]=abc(Irx0,Irx1,Irx2); [Irya,Iryb,Iryc]=abc(Iry0,Iry1,Iry2); [Iyra,Iyrb,Iyrc]=abc(Iyr0,Iyr1,Iyr2); [Iysa,Iysb,Iysc]=abc(Iys0,Iys1,Iys2); % PROTECCION DE DISTANCIA (21) switch(compens21) case{1}, % COMPENSADA [Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG]=... PDistComp(Vsa,Vsb,Vsc,Isxa,Isxb,Isxc,Isx0,k,Km,Isy0); [Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG]=... PDistComp(Vsa,Vsb,Vsc,Isya,Isyb,Isyc,Isy0,k,Km,Isx0); [Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG]=... PDistComp(Vxa,Vxb,Vxc,Ixsa,Ixsb,Ixsc,Ixs0,k,Km,Ixr0); [Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG]=... PDistComp(Vxa,Vxb,Vxc,Ixra,Ixrb,Ixrc,Ixr0,k,Km,Ixs0); [Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG]=... PDistComp(Vra,Vrb,Vrc,Irxa,Irxb,Irxc,Irx0,k,Km,Iry0); [Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG]=... PDistComp(Vra,Vrb,Vrc,Irya,Iryb,Iryc,Iry0,k,Km,Irx0); [Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG]=... PDistComp(Vya,Vyb,Vyc,Iyra,Iyrb,Iyrc,Iyr0,k,Km,Iys0); [Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG]=... PDistComp(Vya,Vyb,Vyc,Iysa,Iysb,Iysc,Iys0,k,Km,Iyr0); case{2}, % NO COMPENSADA [Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG]=... PDist(Vsa,Vsb,Vsc,Isxa,Isxb,Isxc,Isx0,k); [Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG]=... PDist(Vsa,Vsb,Vsc,Isya,Isyb,Isyc,Isy0,k); [Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG]=... PDist(Vxa,Vxb,Vxc,Ixsa,Ixsb,Ixsc,Ixs0,k); [Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG]=... PDist(Vxa,Vxb,Vxc,Ixra,Ixrb,Ixrc,Ixr0,k); [Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG]=... PDist(Vra,Vrb,Vrc,Irxa,Irxb,Irxc,Irx0,k); [Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG]=... PDist(Vra,Vrb,Vrc,Irya,Iryb,Iryc,Iry0,k); [Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG]=... PDist(Vya,Vyb,Vyc,Iyra,Iyrb,Iyrc,Iyr0,k); [Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG]=... PDist(Vya,Vyb,Vyc,Iysa,Iysb,Iysc,Iys0,k); end % 67N [Vspol,Isxop,AngDirSx]=PDire(Vs0,Isx0); [Vspol,Isyop,AngDirSy]=PDire(Vs0,Isy0); [Vxpol,Ixsop,AngDirXs]=PDire(Vx0,Ixs0); [Vxpol,Ixrop,AngDirXr]=PDire(Vx0,Ixr0); [Vrpol,Irxop,AngDirRx]=PDire(Vr0,Irx0); [Vrpol,Iryop,AngDirRy]=PDire(Vr0,Iry0); [Vypol,Iyrop,AngDirYr]=PDire(Vy0,Iyr0); [Vypol,Iysop,AngDirYs]=PDire(Vy0,Iys0); % --------------------------------------------------- switch(pfaltam) case{11}, em1='Extremo S'; em2='Extremo X'; L=m*L1; La=L1+L4; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,La); L=L4+(1-m)*L1; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,La); cadylabel1='alcSmpu'; cadylabel2='alcXmpu'; cadylabel1a='alcSmpua'; cadylabel2a='alcXmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=0; bp=L1/La; as=(L4+L1)/La; bs=-L1/La; % Salida pu Lim1=0.8*(L1+L4)/(m*L1); Lim2=0.8*(L1+L4)/(L4+(1-m)*L1); case{12}, em1='Extremo S'; em2='Extremo X'; L=L1+m*L4; La=L1+L4; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Sx_ZAB,Sx_ZBC,Sx_ZCA,Sx_ZAG,Sx_ZBG,Sx_ZCG,zl1u,La); L=(1-m)*L4; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Xs_ZAB,Xs_ZBC,Xs_ZCA,Xs_ZAG,Xs_ZBG,Xs_ZCG,zl1u,La); cadylabel1='alcSmpu'; cadylabel2='alcXmpu'; cadylabel1a='alcSmpua'; cadylabel2a='alcXmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=L1/La; bp=L4/La; as=L4/La; bs=-L4/La; % Salida pu Lim1=0.8*(L1+L4)/(L1+m*L4); Lim2=0.8*(L1+L4)/((1-m)*L4); case{13}, em1='Extremo X'; em2='Extremo R'; L=(1-m)*L4; La=L4+L2+L3; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,La); L=L3+L2+m*L4; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,La); cadylabel1='alcXmpu'; cadylabel2='alcRmpu'; cadylabel1a='alcXmpua'; cadylabel2a='alcRmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=L4/La; bp=-L4/La; as=(L3+L2)/La; bs=L4/La; % Salida pu Lim1=0.8*(L4+L2+L3)/((1-m)*L4); Lim2=0.8*(L4+L2+L3)/(L3+L2+m*L4); case{14}, em1='Extremo X'; em2='Extremo R'; L=L4+m*L2; La=L4+L2+L3; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,La); L=L3+(1-m)*L2; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,La); cadylabel1='alcXmpu'; cadylabel2='alcRmpu'; cadylabel1a='alcXmpua'; cadylabel2a='alcRmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=L4/La; bp=L2/La; as=(L3+L2)/La; bs=-L2/La; Lim1=0.8*(L4+L2+L3)/(L4+m*L2); Lim2=0.8*(L4+L2+L3)/(L3+(1-m)*L2); case{15}, em1='Extremo X'; em2='Extremo R'; L=L4+L2+m*L3; La=L4+L2+L3; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Xr_ZAB,Xr_ZBC,Xr_ZCA,Xr_ZAG,Xr_ZBG,Xr_ZCG,zl1u,La); L=(1-m)*L3; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Rx_ZAB,Rx_ZBC,Rx_ZCA,Rx_ZAG,Rx_ZBG,Rx_ZCG,zl1u,La); cadylabel1='alcXmpu'; cadylabel2='alcRmpu'; cadylabel1a='alcXmpua'; cadylabel2a='alcRmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=(L4+L2)/La; bp=L3/La; as=L3/La; bs=-L3/La; Lim1=0.8*(L4+L2+L3)/(L4+L2+m*L3); Lim2=0.8*(L4+L2+L3)/((1-m)*L3); case{16}, em1='Extremo R'; em2='Extremo Y'; L=(1-m)*L3; La=L3+L5; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,La); L=L5+m*L3; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,La); cadylabel1='alcRmpu'; cadylabel2='alcYmpu'; cadylabel1a='alcRmpua'; cadylabel2a='alcYmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=L3/La; bp=-L3/La; as=L5/La; bs=L3/La; Lim1=0.8*(L3+L5)/((1-m)*L3); Lim2=0.8*(L3+L5)/(L5+m*L3); case{17}, em1='Extremo R'; em2='Extremo Y'; L=L3+m*L5; La=L3+L5; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Ry_ZAB,Ry_ZBC,Ry_ZCA,Ry_ZAG,Ry_ZBG,Ry_ZCG,zl1u,La); L=(1-m)*L5; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Yr_ZAB,Yr_ZBC,Yr_ZCA,Yr_ZAG,Yr_ZBG,Yr_ZCG,zl1u,La); cadylabel1='alcRmpu'; cadylabel2='alcYmpu'; cadylabel1a='alcRmpua'; cadylabel2a='alcYmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=L3/La; bp=L5/La; as=L5/La; bs=-L5/La; Lim1=0.8*(L3+L5)/(L3+m*L5); Lim2=0.8*(L3+L5)/((1-m)*L5); case{18}, em1='Extremo Y'; em2='Extremo S'; L=(1-m)*L5; La=L5+L2+L1; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,La); L=L1+L2+m*L5; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,La); cadylabel1='alcYmpu'; cadylabel2='alcSmpu'; cadylabel1a='alcYmpua'; cadylabel2a='alcSmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=L5/La; bp=-L5/La; as=(L1+L2)/La; bs=L5/La; Lim1=0.8*(L1+L2+L5)/((1-m)*L5); Lim2=0.8*(L1+L2+L5)/(L1+L2+m*L5); case{19}, em1='Extremo Y'; em2='Extremo S'; L=L5+(1-m)*L2; La=L5+L2+L1; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,La); L=L1+m*L2; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,La); cadylabel1='alcYmpu'; cadylabel2='alcSmpu'; cadylabel1a='alcYmpua'; cadylabel2a='alcSmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=(L5+L2)/La; bp=-L2/La; as=L1/La; bs=L2/La; Lim1=0.8*(L1+L2+L5)/(L5+(1-m)*L2); Lim2=0.8*(L1+L2+L5)/(L1+m*L2); case{20}, em1='Extremo Y'; em2='Extremo S'; L=L5+L2+(1-m)*L1; La=L5+L2+L1; [alcSmABpu,alcSmBCpu,alcSmCApu,alcSmAGpu,alcSmBGpu,alcSmCGpu]=... alcpu(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,L); [alcSmABpua,alcSmBCpua,alcSmCApua,alcSmAGpua,alcSmBGpua,alcSmCGpua]=... alcpua(Ys_ZAB,Ys_ZBC,Ys_ZCA,Ys_ZAG,Ys_ZBG,Ys_ZCG,zl1u,La); L=m*L1; [alcXmABpu,alcXmBCpu,alcXmCApu,alcXmAGpu,alcXmBGpu,alcXmCGpu]=... alcpu(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,L); [alcXmABpua,alcXmBCpua,alcXmCApua,alcXmAGpua,alcXmBGpua,alcXmCGpua]=... alcpua(Sy_ZAB,Sy_ZBC,Sy_ZCA,Sy_ZAG,Sy_ZBG,Sy_ZCG,zl1u,La); cadylabel1='alcYmpu'; cadylabel2='alcSmpu'; cadylabel1a='alcYmpua'; cadylabel2a='alcSmpua'; % Cobertura Extremo Primero: cp=ap+bp*m; % Cobertura Extremo Segundo: cs=as+bs*m; ap=(L5+L2+L1)/La; bp=-L1/La; as=0; bs=L1/La; Lim1=0.8*(L1+L2+L5)/(L5+L2+(1-m)*L1); Lim2=0.8*(L1+L2+L5)/(m*L1); end Abscisa(i)=m; OrdenadaSmAB(1,i)=alcSmABpu; OrdenadaSmBC(1,i)=alcSmBCpu; OrdenadaSmCA(1,i)=alcSmCApu; OrdenadaSmAG(1,i)=alcSmAGpu; OrdenadaSmBG(1,i)=alcSmBGpu; OrdenadaSmCG(1,i)=alcSmCGpu; OrdenadaXmAB(1,i)=alcXmABpu; OrdenadaXmBC(1,i)=alcXmBCpu; OrdenadaXmCA(1,i)=alcXmCApu; OrdenadaXmAG(1,i)=alcXmAGpu; OrdenadaXmBG(1,i)=alcXmBGpu; OrdenadaXmCG(1,i)=alcXmCGpu; OrdenadaSmABa(1,i)=alcSmABpua; OrdenadaSmBCa(1,i)=alcSmBCpua; OrdenadaSmCAa(1,i)=alcSmCApua; OrdenadaSmAGa(1,i)=alcSmAGpua; OrdenadaSmBGa(1,i)=alcSmBGpua; OrdenadaSmCGa(1,i)=alcSmCGpua; OrdenadaXmABa(1,i)=alcXmABpua; OrdenadaXmBCa(1,i)=alcXmBCpua; OrdenadaXmCAa(1,i)=alcXmCApua; OrdenadaXmAGa(1,i)=alcXmAGpua; OrdenadaXmBGa(1,i)=alcXmBGpua; OrdenadaXmCGa(1,i)=alcXmCGpua; %*************************************** %Para salida 2, 67N %*************************************** vectorVspol(1,i)=Vspol; vectorIsxop(1,i)=Isxop; vectorAngDirSx(1,i)=AngDirSx; vectorIsyop(1,i)=Isyop; vectorAngDirSy(1,i)=AngDirSy; vectorVxpol(1,i)=Vxpol; vectorIxsop(1,i)=Ixsop; vectorAngDirXs(1,i)=AngDirXs; vectorIxrop(1,i)=Ixrop; vectorAngDirXr(1,i)=AngDirXr; vectorVrpol(1,i)=Vrpol; vectorIrxop(1,i)=Irxop; vectorAngDirRx(1,i)=AngDirRx; vectorIryop(1,i)=Iryop; vectorAngDirRy(1,i)=AngDirRy; vectorVypol(1,i)=Vypol; vectorIyrop(1,i)=Iyrop; vectorAngDirYr(1,i)=AngDirYr; vectorIysop(1,i)=Iysop; vectorAngDirYs(1,i)=AngDirYs; %***************************** % Formar Vectores para LimIzq y LimDer vectorLim1(1,i)=Lim1; vectorLim2(1,i)=Lim2; i=i+1; m=m+0.01; end %(while m) if faltam==1 cfaltam='ag';elseif faltam==2 cfaltam='bg';elseif faltam==3 cfaltam='cg'; elseif faltam==4 cfaltam='ab';elseif faltam==5 cfaltam='bc';elseif faltam==6 cfaltam='ca'; elseif faltam==7 cfaltam='abg';elseif faltam==8 cfaltam='bcg';elseif faltam==9 cfaltam='cag'; elseif faltam==10 cfaltam='abc';elseif faltam==11 cfaltam='abcg'; end % ---------------------------------------------------------------------------------- % SALIDA GRAFICA % ---------------------------------------------------------------------------------- if compens21==1 comp='21c'; elseif compens21==2 comp='21nc'; end cadenam1=sprintf('PG52c,%s,%3.0fkV,kE=%-4.2f,L1=%3.0fkm,L2=%3.0fkm,L3=%3.0fkm,L4=%3.0fkm,L5=%3.0fkm,',... comp,Un,kE,L1,L2,L3,L4,L5); cadenam2=sprintf('Sccs=%-5.0f \n Sccx=%-5.0f,Sccr=%-5.0f,Sccy=%-5.0f,',Sccs,Sccx,Sccr,Sccy); cadenam3=sprintf('Kx0=%-4.1f,Kr0=%-4.1f,Ky0=%-4.1f,k0=%-4.1f,T=%3.0f,m:%s,Rfm=%-3.0f',... Kx0,Kr0,Ky0,k0,pfaltam,cfaltam,Rfm); if faltam==1|faltam==2|faltam==3|faltam==4|faltam==5|faltam==6|faltam==10 cadenam=sprintf('%s %s %s',cadenam1,cadenam2,cadenam3); elseif faltam==7|faltam==8|faltam==9|faltam==11 cadenam4=sprintf('Rgm=%-3.0f',Rgm); cadenam=sprintf('%s %s %s %s',cadenam1,cadenam2,cadenam3,cadenam4); end %*************************************** % Titulo para las gráficas 67N %*************************************** dircadenam1=sprintf('PG52c,67N,%3.0fkV,kE=%-4.2f,L1=%3.0fkm,L2=%3.0fkm,L3=%3.0fkm,L4=%3.0fkm,L5=%3.0fkm,',... Un,kE,L1,L2,L3,L4,L5); dircadenam2=sprintf('Sccs=%-5.0f \n Sccx=%-5.0f,Sccr=%-5.0f,Sccy=%-5.0f,',Sccs,Sccx,Sccr,Sccy); dircadenam3=sprintf('Kx0=%-4.1f,Kr0=%-4.1f,Ky0=%-4.1f,k0=%-4.1f,T=%3.0f,m:%s,Rfm=%-3.0f',... Kx0,Kr0,Ky0,k0,pfaltam,cfaltam,Rfm); if faltam==1|faltam==2|faltam==3|faltam==4|faltam==5|faltam==6|faltam==10 dircadenam=sprintf('%s %s %s',dircadenam1,dircadenam2,dircadenam3); elseif faltam==7|faltam==8|faltam==9|faltam==11 dircadenam4=sprintf('Rgm=%-3.0f',Rgm); dircadenam=sprintf('%s %s %s %s',dircadenam1,dircadenam2,dircadenam3,dircadenam4); end %************************************* variassalidas=1; while variassalidas==1 fprintf('Salida Grafica 1: Extremos m, Todas las Unidades (pu)\n'); fprintf('Salida Grafica 2: Extremos m, Todas las Unidades (pua)\n'); fprintf('Salida Grafica 3: 67N, Todas las Posiciones\n'); fprintf('\n'); salida=input('Salida: 1,2,3:'); switch(salida) case{1}, % ********* Salida 1 ************ mm=0.01; cont=1; while mm<1 Abscisa0(cont)=mm; Ordenada0(cont)=1; Ordenadacp(cont)=ap+bp*mm; Ordenadacs(cont)=as+bs*mm; cont=cont+1; mm=mm+0.01; end fprintf('%s =1,%s =2\n',em1,em2); extr=input('Extremo:'); %extr=input('Extremo prim.=1, seg.=2: '); switch(extr) case{1}, % Colores b=plot(Abscisa0,Ordenada0); set(b,'Color','c') hold on h=plot(Abscisa,vectorLim1); set(h,'Color',[0.7216 0.4510 0.20]) h=plot(Abscisa0,Ordenadacp); set(h,'Color',[1 0.8431 0]) hold off % Fin Colores hold on plot(Abscisa,OrdenadaSmAB,'r',Abscisa,OrdenadaSmBC,'g',Abscisa,OrdenadaSmCA,'b',... Abscisa,OrdenadaSmAG,'m',Abscisa,OrdenadaSmBG,'k'),grid on hold off % Colores hold on a=plot(Abscisa,OrdenadaSmCG); set(a,'Color',[1 0.5490 0]) hold off % Fin Colores axis([0 1 0 30]); title(cadenam) xlabel('m'),ylabel(cadylabel1) legend('Ref','Lim','Cob','AB','BC','CA','AG','BG','CG') case{2}, % Colores b=plot(Abscisa0,Ordenada0); set(b,'Color','c') hold on h=plot(Abscisa,vectorLim2); set(h,'Color',[0.7216 0.4510 0.20]) h=plot(Abscisa0,Ordenadacs); set(h,'Color',[1 0.8431 0]) hold off % Fin Colores hold on plot(Abscisa,OrdenadaXmAB,'r',Abscisa,OrdenadaXmBC,'g',Abscisa,OrdenadaXmCA,'b',... Abscisa,OrdenadaXmAG,'m',Abscisa,OrdenadaXmBG,'k'),grid on hold off % Colores hold on a=plot(Abscisa,OrdenadaXmCG); set(a,'Color',[1 0.5490 0]) hold off % Fin Colores axis([0 1 0 3]); title(cadenam) xlabel('m'),ylabel(cadylabel2) legend('Ref','Lim','Cob','AB','BC','CA','AG','BG','CG') end case(2) % ********* Salida 2 ************ mm=0.01; cont=1; while mm<1 Abscisa0(cont)=mm; Ordenada0(cont)=0.8; Ordenadacp(cont)=ap+bp*mm; Ordenadacs(cont)=as+bs*mm; cont=cont+1; mm=mm+0.01; end fprintf('%s =1,%s =2\n',em1,em2); extr=input('Extremo:'); %extr=input('Extremo prim.=1, seg.=2: '); switch(extr) case{1}, % Colores b=plot(Abscisa0,Ordenadacp); set(b,'Color','c') % Fin Colores hold on plot(Abscisa,OrdenadaSmABa,'r',Abscisa,OrdenadaSmBCa,'g',Abscisa,OrdenadaSmCAa,'b',... Abscisa,OrdenadaSmAGa,'m',Abscisa,OrdenadaSmBGa,'k'),grid on hold off % Colores hold on a=plot(Abscisa,OrdenadaSmCGa); set(a,'Color',[1 0.5490 0]) h=plot(Abscisa0,Ordenada0); set(h,'Color',[0.7216 0.4510 0.20]) hold off % Fin Colores axis([0 1 -2 2]); title(cadenam) xlabel('m'),ylabel(cadylabel1a) legend('Ref,Cob','AB','BC','CA','AG','BG','CG','z1') case{2}, % Colores b=plot(Abscisa0,Ordenadacs); set(b,'Color','c') % Fin Colores hold on plot(Abscisa,OrdenadaXmABa,'r',Abscisa,OrdenadaXmBCa,'g',Abscisa,OrdenadaXmCAa,'b',... Abscisa,OrdenadaXmAGa,'m',Abscisa,OrdenadaXmBGa,'k'),grid on hold off % Colores hold on a=plot(Abscisa,OrdenadaXmCGa); set(a,'Color',[1 0.5490 0]) h=plot(Abscisa0,Ordenada0); set(h,'Color',[0.7216 0.4510 0.20]) hold off % Fin Colores axis([0 1 -2 2]); title(cadenam) xlabel('m'),ylabel(cadylabel2a) legend('Ref,Cob','AB','BC','CA','AG','BG','CG','z1') end % end del switch extremo %************************************************ case{3} % ********* Salida 3, 67N ************ extr=input('Extremo: Sx=1,Sy=2,Xs=3,Xr=4,Rx=5,Ry=6,Yr=7,Ys=8: '); switch(extr) case{1} subplot(3,1,1) plot(Abscisa,vectorAngDirSx,'b'); axis([0 1 -180 180]); xlabel('m'),ylabel('AngDirSx(º)') title(dircadenam) rayar(1); legend('AngDirSx','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVspol,'b'); grid %title(dircadenam) xlabel('m'),ylabel('modulo Vspol(kV)') subplot(3,1,3) plot(Abscisa,vectorIsxop,'b'); grid xlabel('m'),ylabel('modulo Isxop(kA)') case{2} subplot(3,1,1) plot(Abscisa,vectorAngDirSy,'b'); axis([0 1 -180 180]); xlabel('m'),ylabel('AngDirSy(º)') title(dircadenam) rayar(1); legend('AngDirSy','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVspol,'b'); grid %title(dircadenam) xlabel('m'),ylabel('modulo Vspol(kV)') subplot(3,1,3) plot(Abscisa,vectorIsyop,'b'); grid xlabel('m'),ylabel('modulo Isyop(kA)') case{3} subplot(3,1,1) plot(Abscisa,vectorAngDirXs,'b'); axis([0 1 -180 180]); xlabel('m'),ylabel('AngDirXs(º)') title(dircadenam) rayar(1); legend('AngDirXs','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVxpol,'b'); grid %title(dircadenam) xlabel('m'),ylabel('modulo Vxpol(kV)') subplot(3,1,3) plot(Abscisa,vectorIxsop,'b'); grid xlabel('m'),ylabel('modulo Ixsop(kA)') case{4} subplot(3,1,1) plot(Abscisa,vectorAngDirXr,'b'); axis([0 1 -180 180]); xlabel('m'),ylabel('AngDirXr(º)') title(dircadenam) rayar(1); legend('AngDirXr','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVxpol,'b'); grid %title(dircadenam) xlabel('m'),ylabel('modulo Vxpol(kV)') subplot(3,1,3) plot(Abscisa,vectorIxrop,'b'); grid xlabel('m'),ylabel('modulo Ixrop(kA)') case{5} subplot(3,1,1) plot(Abscisa,vectorAngDirRx,'b'); axis([0 1 -180 180]); xlabel('m'),ylabel('AngDirRx(º)') title(dircadenam) rayar(1); legend('AngDirRx','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVrpol,'b'); grid %title(dircadenam) xlabel('m'),ylabel('modulo Vrpol(kV)') subplot(3,1,3) plot(Abscisa,vectorIrxop,'b'); grid xlabel('m'),ylabel('modulo Irxop(kA)') case{6} subplot(3,1,1) plot(Abscisa,vectorAngDirRy,'b'); axis([0 1 -180 180]); xlabel('m'),ylabel('AngDirRy(º)') title(dircadenam) rayar(1); legend('AngDirRy','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVrpol,'b'); grid %title(dircadenam) xlabel('m'),ylabel('modulo Vrpol(kV)') subplot(3,1,3) plot(Abscisa,vectorIryop,'b'); grid xlabel('m'),ylabel('modulo Iryop(kA)') case{7} subplot(3,1,1) plot(Abscisa,vectorAngDirYr,'b'); axis([0 1 -180 180]); xlabel('m'),ylabel('AngDirYr(º)') title(dircadenam) rayar(1); legend('AngDirYr','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVypol,'b'); grid %title(dircadenam) xlabel('m'),ylabel('modulo Vypol(kV)') subplot(3,1,3) plot(Abscisa,vectorIyrop,'b'); grid xlabel('m'),ylabel('modulo Iyrop(kA)') case{8} subplot(3,1,1) plot(Abscisa,vectorAngDirYs,'b'); axis([0 1 -180 180]); xlabel('m'),ylabel('AngDirYs(º)') title(dircadenam) rayar(1); legend('AngDirYs','AngMaxPar') subplot(3,1,2) plot(Abscisa,vectorVypol,'b'); grid %title(dircadenam) xlabel('m'),ylabel('modulo Vypol(kV)') subplot(3,1,3) plot(Abscisa,vectorIysop,'b'); grid xlabel('m'),ylabel('modulo Iysop(kA)') end % end del switch extr %************************************************ end % end del switch salida variassalidas=input('Mas Salidas=1,Acabar=2 :'); end % end del while variassalidas variasejecuciones=input('Mas Ejecuciones=1,Acabar=2 :'); end % end del while variasejecuciones fprintf(1,'Pulse Enter\n');