function[Im0,Im1,Im2,Vm0,Vm1,Vm2,Ip0,Ip1,Ip2,Vp0,Vp1,Vp2]=... maMDmp(Z0,Z1,Z2,E,falta,Rf,Rg,VbaseA,VbaseB,IbaseA,IbaseB); a=cos(2*pi/3)+sin(2*pi/3)*i; Mmp=zeros(12); switch(falta) case{1}, % Aa %Rf=input('Resistencia de Falta,por fase, Rf(Ohm) = '); Mmp(1,1)=1; Mmp(1,2)=-1; Mmp(2,1)=1; Mmp(2,3)=-1; Mmp(3,7)=1; Mmp(3,8)=-1; Mmp(4,7)=1; Mmp(4,9)=-1; Mmp(5,1)=IbaseA; Mmp(5,7)=IbaseB; Mmp(6,1)=-6*Rf*IbaseA; Mmp(6,4)=VbaseA; Mmp(6,5)=VbaseA; Mmp(6,6)=VbaseA; Mmp(6,10)=-VbaseB; Mmp(6,11)=-VbaseB; Mmp(6,12)=-VbaseB; case{2}, % Bc %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-a^2; Mmp(2,1)=1; Mmp(2,3)=-a; Mmp(3,7)=1; Mmp(3,8)=-a; Mmp(4,7)=1; Mmp(4,9)=-a^2; Mmp(5,1)=IbaseA; Mmp(5,7)=IbaseB; Mmp(6,1)=-6*Rf*IbaseA; Mmp(6,4)=VbaseA; Mmp(6,5)=a^2*VbaseA; Mmp(6,6)=a*VbaseA; Mmp(6,10)=-VbaseB; Mmp(6,11)=-a*VbaseB; Mmp(6,12)=-a^2*VbaseB; case{3}, % Cb %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-a; Mmp(2,1)=1; Mmp(2,3)=-a^2; Mmp(3,7)=1; Mmp(3,8)=-a^2; Mmp(4,7)=1; Mmp(4,9)=-a; Mmp(5,1)=IbaseA; Mmp(5,7)=IbaseB; Mmp(6,1)=-6*Rf*IbaseA; Mmp(6,4)=VbaseA; Mmp(6,5)=a*VbaseA; Mmp(6,6)=a^2*VbaseA; Mmp(6,10)=-VbaseB; Mmp(6,11)=-a^2*VbaseB; Mmp(6,12)=-a*VbaseB; case{4}, % Aag %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); %Rg=input('Resistencia de Falta,comun, Rg(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-1; Mmp(2,1)=1; Mmp(2,3)=-1; Mmp(3,7)=1; Mmp(3,8)=-1; Mmp(4,7)=1; Mmp(4,9)=-1; Mmp(5,1)=-3*(Rf+Rg)*IbaseA; Mmp(5,4)=VbaseA; Mmp(5,5)=VbaseA; Mmp(5,6)=VbaseA; Mmp(5,7)=-3*Rg*IbaseB; Mmp(6,1)=-3*Rg*IbaseA; Mmp(6,7)=-3*(Rf+Rg)*IbaseB; Mmp(6,10)=VbaseB; Mmp(6,11)=VbaseB; Mmp(6,12)=VbaseB; case{5}, % Bcg %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); %Rg=input('Resistencia de Falta,comun, Rg(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-a^2; Mmp(2,1)=1; Mmp(2,3)=-a; Mmp(3,7)=1; Mmp(3,8)=-a; Mmp(4,7)=1; Mmp(4,9)=-a^2; Mmp(5,1)=-3*(Rf+Rg)*IbaseA; Mmp(5,4)=VbaseA; Mmp(5,5)=a^2*VbaseA; Mmp(5,6)=a*VbaseA; Mmp(5,7)=-3*Rg*IbaseB; Mmp(6,1)=-3*Rg*IbaseA; Mmp(6,7)=-3*(Rf+Rg)*IbaseB; Mmp(6,10)=VbaseB; Mmp(6,11)=a*VbaseB; Mmp(6,12)=a^2*VbaseB; case{6}, % Cbg %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); %Rg=input('Resistencia de Falta,comun, Rg(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-a; Mmp(2,1)=1; Mmp(2,3)=-a^2; Mmp(3,7)=1; Mmp(3,8)=-a^2; Mmp(4,7)=1; Mmp(4,9)=-a; Mmp(5,1)=-3*(Rf+Rg)*IbaseA; Mmp(5,4)=VbaseA; Mmp(5,5)=a*VbaseA; Mmp(5,6)=a^2*VbaseA; Mmp(5,7)=-3*Rg*IbaseB; Mmp(6,1)=-3*Rg*IbaseA; Mmp(6,7)=-3*(Rf+Rg)*IbaseB; Mmp(6,10)=VbaseB; Mmp(6,11)=a^2*VbaseB; Mmp(6,12)=a*VbaseB; case{7}, % Aab %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-1; Mmp(2,1)=1; Mmp(2,3)=-1; Mmp(3,7)=1; Mmp(3,8)=a; Mmp(3,9)=a^2; Mmp(4,1)=IbaseA; Mmp(4,7)=IbaseB; Mmp(5,1)=-4*Rf*IbaseA; Mmp(5,4)=VbaseA; Mmp(5,5)=VbaseA; Mmp(5,6)=VbaseA; Mmp(5,8)=Rf*IbaseB; Mmp(5,9)=Rf*IbaseB; Mmp(5,10)=-VbaseB; Mmp(5,11)=-VbaseB; Mmp(5,12)=-VbaseB; Mmp(6,1)=-4*Rf*IbaseA; Mmp(6,4)=VbaseA; Mmp(6,5)=VbaseA; Mmp(6,6)=VbaseA; Mmp(6,8)=Rf*a^2*IbaseB; Mmp(6,9)=Rf*a*IbaseB; Mmp(6,10)=-VbaseB; Mmp(6,11)=-a^2*VbaseB; Mmp(6,12)=-a*VbaseB; case{8}, % Bac %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-a^2; Mmp(2,1)=1; Mmp(2,3)=-a; Mmp(3,7)=1; Mmp(3,8)=a^2; Mmp(3,9)=a; Mmp(4,1)=IbaseA; Mmp(4,7)=IbaseB; Mmp(5,1)=-4*Rf*IbaseA; Mmp(5,4)=VbaseA; Mmp(5,5)=a^2*VbaseA; Mmp(5,6)=a*VbaseA; Mmp(5,8)=Rf*IbaseB; Mmp(5,9)=Rf*IbaseB; Mmp(5,10)=-VbaseB; Mmp(5,11)=-VbaseB; Mmp(5,12)=-VbaseB; Mmp(6,1)=-4*Rf*IbaseA; Mmp(6,4)=VbaseA; Mmp(6,5)=a^2*VbaseA; Mmp(6,6)=a*VbaseA; Mmp(6,8)=Rf*a*IbaseB; Mmp(6,9)=Rf*a^2*IbaseB; Mmp(6,10)=-VbaseB; Mmp(6,11)=-a*VbaseB; Mmp(6,12)=-a^2*VbaseB; case{9}, % Cbc %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-a; Mmp(2,1)=1; Mmp(2,3)=-a^2; Mmp(3,7)=1; Mmp(3,8)=1; Mmp(3,9)=1; Mmp(4,1)=IbaseA; Mmp(4,7)=IbaseA; Mmp(5,1)=-4*Rf*IbaseA; Mmp(5,4)=VbaseA; Mmp(5,5)=a*VbaseA; Mmp(5,6)=a^2*VbaseA; Mmp(5,8)=Rf*a^2*IbaseB; Mmp(5,9)=Rf*a*IbaseB; Mmp(5,10)=-VbaseB; Mmp(5,11)=-a^2*VbaseB; Mmp(5,12)=-a*VbaseB; Mmp(6,1)=-4*Rf*IbaseA; Mmp(6,4)=VbaseA; Mmp(6,5)=a*VbaseA; Mmp(6,6)=a^2*VbaseA; Mmp(6,8)=Rf*a*IbaseB; Mmp(6,9)=Rf*a^2*IbaseB; Mmp(6,10)=-VbaseB; Mmp(6,11)=-a*VbaseB; Mmp(6,12)=-a^2*VbaseB; case{10}, % Aabg %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); %Rg=input('Resistencia de Falta,comun, Rg(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-1; Mmp(2,1)=1; Mmp(2,3)=-1; Mmp(3,7)=1; Mmp(3,8)=a; Mmp(3,9)=a^2; Mmp(4,1)=-3*(Rf+Rg)*IbaseA; Mmp(4,4)=VbaseA; Mmp(4,5)=VbaseA; Mmp(4,6)=VbaseA; Mmp(4,7)=-3*Rg*IbaseB; Mmp(5,1)=-3*Rf*IbaseA; Mmp(5,4)=VbaseA; Mmp(5,5)=VbaseA; Mmp(5,6)=VbaseA; Mmp(5,7)=Rf*IbaseB; Mmp(5,8)=Rf*IbaseB; Mmp(5,9)=Rf*IbaseB; Mmp(5,10)=-VbaseB; Mmp(5,11)=-VbaseB; Mmp(5,12)=-VbaseB; Mmp(6,1)=-3*Rf*IbaseA; Mmp(6,4)=VbaseA; Mmp(6,5)=VbaseA; Mmp(6,6)=VbaseA; Mmp(6,7)=Rf*(a+2)*IbaseB; Mmp(6,8)=-Rf*(a+2)*IbaseB; Mmp(6,10)=-VbaseB; Mmp(6,11)=-a^2*VbaseB; Mmp(6,12)=-a*VbaseB; case{11}, % Bacg %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); %Rg=input('Resistencia de Falta,comun, Rg(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-a^2; Mmp(2,1)=1; Mmp(2,3)=-a; Mmp(3,7)=1; Mmp(3,8)=a^2; Mmp(3,9)=a; Mmp(4,1)=-3*(Rf+Rg)*IbaseA; Mmp(4,4)=VbaseA; Mmp(4,5)=a^2*VbaseA; Mmp(4,6)=a*VbaseA; Mmp(4,7)=-3*Rg*IbaseB; Mmp(5,1)=-3*Rg*IbaseA; Mmp(5,7)=(Rf*(a^2-1)-3*Rg)*IbaseB; Mmp(5,8)=-Rf(1-a)*IbaseB; Mmp(5,10)=VbaseB; Mmp(5,11)=VbaseB; Mmp(5,12)=VbaseB; Mmp(6,1)=-3*Rg*IbaseA; Mmp(6,7)=(Rf*(a-1)-3*Rg)*IbaseB; Mmp(6,8)=Rf*(1-a)*IbaseB; Mmp(6,10)=VbaseB; Mmp(6,11)=a*VbaseB; Mmp(6,12)=a^2*VbaseB; case{12}, % Cbcg %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); %Rg=input('Resistencia de Falta,comun, Rg(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-a; Mmp(2,1)=1; Mmp(2,3)=-a^2; Mmp(3,7)=1; Mmp(3,8)=1; Mmp(3,9)=1; Mmp(4,1)=-3*(Rf+Rg)*IbaseA; Mmp(4,4)=VbaseA; Mmp(4,5)=a*VbaseA; Mmp(4,6)=a^2*VbaseA; Mmp(4,7)=-3*Rg*IbaseB; Mmp(5,1)=-3*Rg*IbaseA; Mmp(5,7)=(Rf*(a-1)-3*Rg)*IbaseB; Mmp(5,8)=-Rf*(a^2-a)*IbaseB; Mmp(5,10)=VbaseB; Mmp(5,11)=a^2*VbaseB; Mmp(5,12)=a*VbaseB; Mmp(6,1)=-3*Rg*IbaseA; Mmp(6,7)=(-Rf*(2+a)-3*Rg)*IbaseB; Mmp(6,8)=Rf*(a^2-a)*IbaseB; Mmp(6,10)=VbaseB; Mmp(6,11)=a*VbaseB; Mmp(6,12)=a^2*VbaseB; case{13}, % Aabc %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-1; Mmp(2,1)=1; Mmp(2,3)=-1; Mmp(3,1)=IbaseA; Mmp(3,7)=IbaseB; Mmp(4,1)=-4*Rf*IbaseA; Mmp(4,4)=VbaseA; Mmp(4,5)=VbaseA; Mmp(4,6)=VbaseA; Mmp(4,8)=Rf*IbaseB; Mmp(4,9)=Rf*IbaseB; Mmp(4,10)=-VbaseB; Mmp(4,11)=-VbaseB; Mmp(4,12)=-VbaseB; Mmp(5,1)=-4*Rf*IbaseA; Mmp(5,4)=VbaseA; Mmp(5,5)=VbaseA; Mmp(5,6)=VbaseA; Mmp(5,8)=Rf*a^2*IbaseB; Mmp(5,9)=Rf*a*IbaseB; Mmp(5,10)=-VbaseB; Mmp(5,11)=-a^2*VbaseB; Mmp(5,12)=-a*VbaseB; Mmp(6,1)=-4*Rf*IbaseA; Mmp(6,4)=VbaseA; Mmp(6,5)=VbaseA; Mmp(6,6)=VbaseA; Mmp(6,8)=Rf*a*IbaseB; Mmp(6,9)=Rf*a^2*IbaseB; Mmp(6,10)=-VbaseB; Mmp(6,11)=-a*VbaseB; Mmp(6,12)=-a^2*VbaseB; case{14}, % Aabcg %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); %Rg=input('Resistencia de Falta,comun, Rg(Ohm): '); Mmp(1,1)=1; Mmp(1,2)=-1; Mmp(2,1)=1; Mmp(2,3)=-1; Mmp(3,1)=-3*(Rf+Rg)*IbaseA; Mmp(3,4)=VbaseA; Mmp(3,5)=VbaseA; Mmp(3,6)=VbaseA; Mmp(3,7)=-3*Rg*IbaseB; Mmp(4,1)=-3*Rg*IbaseA; Mmp(4,7)=-(Rf+3*Rg)*IbaseB; Mmp(4,8)=-Rf*IbaseB; Mmp(4,9)=-Rf*IbaseB; Mmp(4,10)=VbaseB; Mmp(4,11)=VbaseB; Mmp(4,12)=VbaseB; Mmp(5,1)=-3*Rg*IbaseA; Mmp(5,7)=-(Rf+3*Rg); Mmp(5,8)=-a^2*Rf; Mmp(5,9)=-a*Rf; Mmp(5,10)=VbaseB; Mmp(5,11)=a^2*VbaseB; Mmp(5,12)=a*VbaseB; Mmp(6,1)=-3*Rg*IbaseA; Mmp(6,7)=-(Rf+3*Rg)*IbaseB; Mmp(6,8)=-a*Rf*IbaseB; Mmp(6,9)=-a^2*Rf*IbaseB; Mmp(6,10)=VbaseB; Mmp(6,11)=a*VbaseB; Mmp(6,12)=a^2*VbaseB; case{15}, % ABCabc %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); Mmp(1,1)=IbaseA; Mmp(1,7)=IbaseB; Mmp(2,2)=-Rf*IbaseA; Mmp(2,5)=VbaseA; Mmp(3,3)=-Rf*IbaseA; Mmp(3,6)=VbaseA; Mmp(4,1)=-Rf*IbaseA; Mmp(4,4)=VbaseA; Mmp(4,7)=Rf*IbaseB; Mmp(4,10)=-VbaseB; Mmp(5,8)=-Rf*IbaseB; Mmp(5,11)=VbaseB; Mmp(6,9)=-Rf*IbaseB; Mmp(6,12)=VbaseB; case{16}, % ABCabcg %Rf=input('Resistencia de Falta,por fase, Rf(Ohm): '); %Rg=input('Resistencia de Falta,comun, Rg(Ohm): '); Mmp(1,1)=-(Rf+3*Rg)*IbaseA; Mmp(1,4)=VbaseA; Mmp(1,7)=-3*Rg*IbaseB; Mmp(2,2)=-Rf*IbaseA; Mmp(2,5)=VbaseA; Mmp(3,3)=-Rf*IbaseA; Mmp(3,6)=VbaseA; Mmp(4,1)=-3*Rg*IbaseA; Mmp(4,7)=-(Rf+3*Rg)*IbaseB; Mmp(4,10)=VbaseB; Mmp(5,8)=-Rf*IbaseB; Mmp(5,11)=VbaseB; Mmp(6,9)=-Rf*IbaseB; Mmp(6,12)=VbaseB; end Mmp(7,1)=Z0(1,1);Mmp(7,4)=1;Mmp(7,7)=Z0(1,2); Mmp(8,1)=Z0(2,1);Mmp(8,7)=Z0(2,2);Mmp(8,10)=1; Mmp(9,2)=Z1(1,1);Mmp(9,5)=1;Mmp(9,8)=Z1(1,2); Mmp(10,2)=Z1(2,1);Mmp(10,8)=Z1(2,2);Mmp(10,11)=1; Mmp(11,3)=Z2(1,1);Mmp(11,6)=1;Mmp(11,9)=Z2(1,2); Mmp(12,3)=Z2(2,1);Mmp(12,9)=Z2(2,2);Mmp(12,12)=1; P=[0;0;0;0;0;0;0;0;E;E;0;0]; N=Mmp\P; % ------------------------------------------------------- % Im0=N(1);Im1=N(2);Im2=N(3); Vm0=N(4);Vm1=N(5);Vm2=N(6); Ip0=N(7);Ip1=N(8);Ip2=N(9); Vp0=N(10);Vp1=N(11);Vp2=N(12); 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  >  PGTf  >  PG128MTf  >  maMDmp.m