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