Subprogramas de PFCc24 ********************************* function [Aa,Ab,Ac]=abc(A0,A1,A2) a=cos(2*pi/3)+sin(2*pi/3)*i; T=zeros(3); T(1,1)=1;T(1,2)=1;T(1,3)=1; T(2,1)=1;T(2,2)=a^2;T(2,3)=a; T(3,1)=1;T(3,2)=a;T(3,3)=a^2; A=[A0;A1;A2]; B=T*A; Aa=B(1);Ab=B(2);Ac=B(3); ********************************* function [mod,arg]=mod_arg(complejocart) mod=abs(complejocart); argrad=angle(complejocart); arg=argrad*(180/pi); ********************************* function [mod,arg]=mod_argz(complejocart) mod=abs(complejocart); if mod>=10000 mod=exp(1000); end argrad=angle(complejocart); arg=argrad*(180/pi); ********************************* function[Vpol,Iop,AngDir]=PDire(V0,I0) Vpol=-3*V0; Iop=3*I0; AngDir=angle(Vpol/Iop); Vpol=abs(Vpol); Iop=abs(Iop); AngDir=-(AngDir/pi)*180; ********************************* %[Xm_ZAB,Xm_ZBC,Xm_ZCA,Xm_ZAG,Xm_ZBG,Xm_ZCG]=... % PDist(Vxa,Vxb,Vxc,Ixma,Ixmb,Ixmc,Ixm0,k); function[ZAB,ZBC,ZCA,ZAG,ZBG,ZCG]=PDist(Va,Vb,Vc,Ia,Ib,Ic,I0,k); ZAB=(Va-Vb)/(Ia-Ib); ZBC=(Vb-Vc)/(Ib-Ic); ZCA=(Vc-Va)/(Ic-Ia); ZAG=Va/(Ia+k*3*I0); ZBG=Vb/(Ib+k*3*I0); ZCG=Vc/(Ic+k*3*I0); ********************************* function[ZAB,ZBC,ZCA,ZAG,ZBG,ZCG]=PDistComp(Va,Vb,Vc,Ia,Ib,Ic,I0,k,Km,I0p) ZAB=(Va-Vb)/(Ia-Ib); ZBC=(Vb-Vc)/(Ib-Ic); ZCA=(Vc-Va)/(Ic-Ia); ZAG=Va/(Ia+k*3*I0+Km*I0p); ZBG=Vb/(Ib+k*3*I0+Km*I0p); ZCG=Vc/(Ic+k*3*I0+Km*I0p); ********************************* function [Ie0]=subIe0(Minc,V0,z0) Ve0=Minc*V0; Ie0=z0\Ve0; ********************************* function[Ie1]=subIe1(Minca,V1a,z1) Ve1=Minca*V1a; Ie1=z1\Ve1; ********************************* function [Ie2]=subIe2(Minc,V2,z2); Ve2=Minc*V2; Ie2=z2\Ve2; ********************************* function [Minc]=subMinc1(u) % Red Normal Minc=zeros(6,4); Minc(1,1)=-1; Minc(2,2)=-1; Minc(3,3)=-1; Minc(3,1)=1; Minc(4,4)=-1; Minc(4,1)=1; Minc(5,3)=-1; Minc(5,2)=1; Minc(6,4)=-1; Minc(6,2)=1; ********************************* function [Minc]=subMinc1sr(u) % Red Normal. Con zsr Minc=zeros(7,4); Minc(1,1)=-1; Minc(2,2)=-1; Minc(3,1)=1; Minc(3,3)=-1; Minc(4,1)=1; Minc(4,4)=-1; Minc(5,2)=1; Minc(5,3)=-1; Minc(6,2)=1; Minc(6,4)=-1; Minc(7,1)=1; Minc(7,2)=-1; % zsr ********************************* function [Minc]=subMinc2(u) % Sm fuera Minc=zeros(7,4); Minc(1,1)=-1; Minc(2,2)=-1; Minc(3,1)=1; Minc(3,3)=-1; Minc(4,1)=1; Minc(4,4)=-1; Minc(5,2)=1; Minc(5,3)=-1; Minc(6,2)=1; Minc(6,4)=-1; Minc(7,1)=1; Minc(7,3)=-1; % anula Sm ********************************* function [Minc]=subMinc2sr(u) % Sm fuera. Con zsr Minc=zeros(8,4); Minc(1,1)=-1; Minc(2,2)=-1; Minc(3,1)=1; Minc(3,3)=-1; Minc(4,1)=1; Minc(4,4)=-1; Minc(5,2)=1; Minc(5,3)=-1; Minc(6,2)=1; Minc(6,4)=-1; Minc(7,1)=1; Minc(7,2)=-1; % zsr Minc(8,1)=1; Minc(8,3)=-1; % anula Sm ********************************* function [Minc]=subMinc3(u) % Rm fuera Minc=zeros(7,4); Minc(1,1)=-1; Minc(2,2)=-1; Minc(3,1)=1; Minc(3,3)=-1; Minc(4,1)=1; Minc(4,4)=-1; Minc(5,2)=1; Minc(5,3)=-1; Minc(6,2)=1; Minc(6,4)=-1; Minc(7,2)=1; Minc(7,3)=-1; % Anula Rm ********************************* function [Minc]=subMinc3sr(u) % Rm fuera. Con zsr Minc=zeros(7,4); Minc(1,1)=-1; Minc(2,2)=-1; Minc(3,1)=1; Minc(3,3)=-1; Minc(4,1)=1; Minc(4,4)=-1; Minc(5,2)=1; Minc(5,3)=-1; Minc(6,2)=1; Minc(6,4)=-1; Minc(7,1)=1; Minc(7,2)=-1; % zsr Minc(8,2)=1; Minc(8,3)=-1; % anula Rm ********************************* function [Minca]=subMinca1(u); % Red Normal Minca=zeros(6,5); Minca(1,5)=1; Minca(1,1)=-1; Minca(2,5)=1; Minca(2,2)=-1; Minca(3,1)=1; Minca(3,3)=-1; Minca(4,1)=1; Minca(4,4)=-1; Minca(5,2)=1; Minca(5,3)=-1; Minca(6,2)=1; Minca(6,4)=-1; ********************************* function [Minca]=subMinca1sr(u) % Red Normal. Con zsr Minca=zeros(7,5); Minca(1,5)=1; Minca(1,1)=-1; Minca(2,5)=1; Minca(2,2)=-1; Minca(3,1)=1; Minca(3,3)=-1; Minca(4,1)=1; Minca(4,4)=-1; Minca(5,2)=1; Minca(5,3)=-1; Minca(6,2)=1; Minca(6,4)=-1; Minca(7,1)=1; Minca(7,2)=-1; % zsr ********************************* function [Minca]=subMinca2(u) % Sm fuera Minca=zeros(7,5); Minca(1,5)=1; Minca(1,1)=-1; Minca(2,5)=1; Minca(2,2)=-1; Minca(3,1)=1; Minca(3,3)=-1; Minca(4,1)=1; Minca(4,4)=-1; Minca(5,2)=1; Minca(5,3)=-1; Minca(6,2)=1; Minca(6,4)=-1; Minca(7,1)=1; Minc(7,3)=-1; % Anula Sm ********************************* function [Minca]=subMinca2sr(u) % Sm fuera . Con zsr Minca=zeros(8,5); Minca(1,5)=1; Minca(1,1)=-1; Minca(2,5)=1; Minca(2,2)=-1; Minca(3,1)=1; Minca(3,3)=-1; Minca(4,1)=1; Minca(4,4)=-1; Minca(5,2)=1; Minca(5,3)=-1; Minca(6,2)=1; Minca(6,4)=-1; Minca(7,1)=1; Minca(7,2)=-1; % zsr Minca(8,1)=1; Minca(8,3)=-1; % Anula Sm ********************************* function [Minca]=subMinca3(u) % Rm fuera Minca=zeros(7,5); Minca(1,5)=1; Minca(1,1)=-1; Minca(2,5)=1; Minca(2,2)=-1; Minca(3,1)=1; Minca(3,3)=-1; Minca(4,1)=1; Minca(4,4)=-1; Minca(5,2)=1; Minca(5,3)=-1; Minca(6,2)=1; Minca(6,4)=-1; Minc(7,2)=1; Minc(7,1)=-1; % Anula Rm ********************************* function [Minca]=subMinca3sr(u) % Rm fuera . Con zsr Minca=zeros(8,5); Minca(1,5)=1; Minca(1,1)=-1; Minca(2,5)=1; Minca(2,2)=-1; Minca(3,1)=1; Minca(3,3)=-1; Minca(4,1)=1; Minca(4,4)=-1; Minca(5,2)=1; Minca(5,3)=-1; Minca(6,2)=1; Minca(6,4)=-1; Minca(7,1)=1; Minca(7,2)=-1; % zsr Minca(8,2)=1; Minca(8,1)=-1; % Anula Rm ********************************* function [V0]=subV0(Y0,Im0,Ip0) I0=[-Im0;-Ip0;0;0;0;0]; V0=Y0\I0; ********************************* function [V1]=subV1(Y1,Im1,Ip1,E) I1=[-Im1;-Ip1;0;0;0;0]; V1=Y1\I1; V1=V1+E; ********************************* function [V2]=subV2(Y2,Im2,Ip2) I2=[-Im2;-Ip2;0;0;0;0]; V2=Y2\I2; ********************************* function [Y0]=subY0(Minc,z0); y0=inv(z0); Y0=Minc.'*y0*Minc; ********************************* function[Y1]=subY1(Minc,z1) y1=inv(z1); Y1=Minc.'*y1*Minc; ********************************* function[Y2]=subY2(Minc,z2) y2=inv(z2); Y2=Minc.'*y2*Minc; ********************************* function[z0]=subz0_1(zs0,zr0,zl0,zlM0,m,p) z0=zeros(6); z0(1,1)=zs0; z0(2,2)=zr0; z0(3,3)=zl0*m; z0(4,4)=zl0*p; z0(5,5)=zl0*(1-m); z0(6,6)=zl0*(1-p); if m<=p z0(3,4)=zlM0*m; z0(4,5)=-zlM0*(p-m); z0(5,6)=zlM0*(1-p); z0(4,3)=z0(3,4); z0(5,4)=z0(4,5); z0(6,5)=z0(5,6); else z0(3,4)=zlM0*p; z0(3,6)=-zlM0*(m-p); z0(5,6)=zlM0*(1-m); z0(4,3)=z0(3,4); z0(6,3)=z0(3,6); z0(6,5)=z0(5,6); end ********************************* function[z0]=subz0_1sr(zs0,zsr0,zr0,zl0,zlM0,m,p) z0=zeros(7); z0(1,1)=zs0; z0(2,2)=zr0; z0(3,3)=m*zl0; z0(4,4)=zl0*p; z0(5,5)=zl0*(1-m); z0(6,6)=zl0*(1-p); z0(7,7)=zsr0; if m<=p z0(3,4)=zlM0*m; z0(4,5)=-zlM0*(p-m); z0(5,6)=zlM0*(1-p); z0(4,3)=z0(3,4); z0(5,4)=z0(4,5); z0(6,5)=z0(5,6); else z0(3,4)=zlM0*p; z0(3,6)=-zlM0*(m-p); z0(5,6)=zlM0*(1-m); z0(4,3)=z0(3,4); z0(6,3)=z0(3,6); z0(6,5)=z0(5,6); end ********************************* function[z0]=subz0_2(zs0,zr0,zl0,zlM0,m,p) z0=zeros(7); z0(1,1)=zs0; z0(2,2)=zr0; z0(3,3)=m*zl0; z0(4,4)=zl0*p; z0(5,5)=zl0*(1-m); z0(6,6)=zl0*(1-p); z0(7,7)=-z0(3,3); if m<=p z0(3,4)=zlM0*m; z0(4,5)=-zlM0*(p-m); z0(5,6)=zlM0*(1-p); z0(7,4)=z0(3,4); z0(4,3)=z0(3,4); z0(5,4)=z0(4,5); z0(6,5)=z0(5,6); z0(4,7)=z0(7,4); else z0(3,4)=zlM0*p; z0(3,6)=-zlM0*(m-p); z0(5,6)=zlM0*(1-m); z0(7,4)=z0(3,4); z0(7,6)=z0(3,6); z0(4,3)=z0(3,4); z0(6,3)=z0(3,6); z0(6,5)=z0(5,6); z0(4,7)=z0(7,4); z0(6,7)=z0(7,6); end ********************************* function[z0]=subz0_2(zs0,zsr0,zr0,zl0,zlM0,m,p) z0=zeros(8); z0(1,1)=zs0; z0(2,2)=zr0; z0(3,3)=m*zl0; z0(4,4)=zl0*p; z0(5,5)=zl0*(1-m); z0(6,6)=zl0*(1-p); z0(7,7)=-z0(3,3); z0(8,8)=zsr0; if m<=p z0(3,4)=zlM0*m; z0(4,5)=-zlM0*(p-m); z0(5,6)=zlM0*(1-p); z0(7,4)=z0(3,4); z0(4,3)=z0(3,4); z0(5,4)=z0(4,5); z0(6,5)=z0(5,6); z0(4,7)=z0(7,4); else z0(3,4)=zlM0*p; z0(3,6)=-zlM0*(m-p); z0(5,6)=zlM0*(1-m); z0(7,4)=z0(3,4); z0(7,6)=z0(3,6); z0(4,3)=z0(3,4); z0(6,3)=z0(3,6); z0(6,5)=z0(5,6); z0(4,7)=z0(7,4); z0(6,7)=z0(7,6); end ********************************* function[z0]=subz0_3(zs0,zr0,zl0,zlM0,m,p) z0=zeros(7); z0(1,1)=zs0; z0(2,2)=zr0; z0(3,3)=m*zl0; z0(4,4)=zl0*p; z0(5,5)=zl0*(1-m); z0(6,6)=zl0*(1-p); z0(7,7)=-z0(5,5); if m<=p z0(3,4)=zlM0*m; z0(4,5)=-zlM0*(p-m); z0(5,6)=zlM0*(1-p); z0(4,7)=z0(4,5); z0(7,6)=z0(5,6); z0(4,3)=z0(3,4); z0(5,4)=z0(4,5); z0(6,5)=z0(5,6); z0(7,4)=z0(4,7); z0(6,7)=z0(7,6); else z0(3,4)=zlM0*p; z0(3,6)=-zlM0*(m-p); z0(5,6)=zlM0*(1-m); z0(7,6)=z0(5,6); z0(4,3)=z0(3,4); z0(6,3)=z0(3,6); z0(6,5)=z0(5,6); z0(6,7)=z0(7,6); end ********************************* function[z0]=subz0_3sr(zs0,zsr0,zr0,zl0,zlM0,m,p) z0=zeros(8); z0(1,1)=zs0; z0(2,2)=zr0; z0(3,3)=m*zl0; z0(4,4)=p*zl0; z0(5,5)=(1-m)*zl0; z0(6,6)=zl0*(1-p); z0(7,7)=-z0(5,5); z0(8,8)=zsr0; if m<=p z0(3,4)=zlM0*m; z0(4,5)=-zlM0*(p-m); z0(5,6)=zlM0*(1-p); z0(4,7)=z0(4,5); z0(7,6)=z0(5,6); z0(4,3)=z0(3,4); z0(5,4)=z0(4,5); z0(6,5)=z0(5,6); z0(7,4)=z0(4,7); z0(6,7)=z0(7,6); else z0(3,4)=zlM0*p; z0(3,6)=-zlM0*(m-p); z0(5,6)=zlM0*(1-m); z0(7,6)=z0(5,6); z0(4,3)=z0(3,4); z0(6,3)=z0(3,6); z0(6,5)=z0(5,6); z0(6,7)=z0(7,6); end ********************************* function[z1]=subz1_1(zs1,zr1,zl1,m,p) z1=zeros(6); z1(1,1)=zs1; z1(2,2)=zr1; z1(3,3)=zl1*m; z1(4,4)=zl1*p; z1(5,5)=zl1*(1-m); z1(6,6)=zl1*(1-p); ********************************* function[z1]=subz1_1sr(zs1,zsr1,zr1,zl1,m,p); z1=zeros(7); z1(1,1)=zs1; z1(2,2)=zr1; z1(3,3)=zl1*m; z1(4,4)=zl1*p; z1(5,5)=zl1*(1-m); z1(6,6)=zl1*(1-p); z1(7,7)=zsr1; ********************************* function[z1]=subz1_2(zs1,zr1,zl1,m,p); z1=zeros(7); z1(1,1)=zs1; z1(2,2)=zr1; z1(3,3)=zl1*m; z1(4,4)=zl1*p; z1(5,5)=zl1*(1-m); z1(6,6)=zl1*(1-p); z1(7,7)=-z1(3,3); ********************************* function[z1]=subz1_2sr(zs1,zsr1,zr1,zl1,m,p); z1=zeros(8); z1(1,1)=zs1; z1(2,2)=zr1; z1(3,3)=zl1*m; z1(4,4)=zl1*p; z1(5,5)=zl1*(1-m); z1(6,6)=zl1*(1-p); z1(7,7)=-z1(3,3); z1(8,8)=zsr1; ********************************* function[z1]=subz1_3(zs1,zr1,zl1,m,p); z1=zeros(7); z1(1,1)=zs1; z1(2,2)=zr1; z1(3,3)=zl1*m; z1(4,4)=zl1*p; z1(5,5)=zl1*(1-m); z1(6,6)=zl1*(1-p); z1(7,7)=-z1(5,5); ********************************* function[z1]=subz1_3sr(zs1,zsr1,zr1,zl1,m,p); z1=zeros(8); z1(1,1)=zs1; z1(2,2)=zr1; z1(3,3)=zl1*m; z1(4,4)=zl1*p; z1(5,5)=zl1*(1-m); z1(6,6)=zl1*(1-p); z1(7,7)=-z1(5,5); z1(8,8)=zsr1; ********************************* function[z2]=subz2_1(zs2,zr2,zl2,m,p) z2=zeros(6); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=zl2*m; z2(4,4)=zl2*p; z2(5,5)=zl2*(1-m); z2(6,6)=zl2*(1-p); ********************************* function[z2]=subz2_1sr(zs2,zsr2,zr2,zl2,m,p) z2=zeros(7); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=zl2*m; z2(4,4)=zl2*p; z2(5,5)=zl2*(1-m); z2(6,6)=zl2*(1-p); z2(7,7)=zsr2; ********************************* function[z2]=subz2_2(zs2,zr2,zl2,m,p) z2=zeros(7); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=zl2*m; z2(4,4)=zl2*p; z2(5,5)=zl2*(1-m); z2(6,6)=zl2*(1-p); z2(7,7)=-z2(3,3); ********************************* function[z2]=subz2_2sr(zs2,zsr2,zr2,zl2,m,p) z2=zeros(8); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=zl2*m; z2(4,4)=zl2*p; z2(5,5)=zl2*(1-m); z2(6,6)=zl2*(1-p); z2(7,7)=-z2(3,3); z2(8,8)=zsr2; ********************************* function[z2]=subz2_3(zs2,zr2,zl2,m,p) z2=zeros(7); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=zl2*m; z2(4,4)=zl2*p; z2(5,5)=zl2*(1-m); z2(6,6)=zl2*(1-p); z2(7,7)=-z2(5,5); ********************************* function[z2]=subz2_3sr(zs2,zsr2,zr2,zl2,m,p) z2=zeros(8); z2(1,1)=zs2; z2(2,2)=zr2; z2(3,3)=zl2*m; z2(4,4)=zl2*p; z2(5,5)=zl2*(1-m); z2(6,6)=zl2*(1-p); z2(7,7)=-z2(5,5); z2(8,8)=zsr2; ********************************* function[Zmedreal,Zmedimag]=zmed(Zmedida,zbase); Zmedreal=real(Zmedida)*zbase; if Zmedreal>=10000 Zmedreal=exp(1000); elseif Zmedreal<=-10000 Zmedreal=log(0); end Zmedimag=imag(Zmedida)*zbase; if Zmedimag>=10000 Zmedimag=exp(1000); elseif Zmedimag<=-10000 Zmedimag=log(0); end ********************************* 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  >  PGn  >  PG34ni  >  Subprogr.m