function[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);
Rgm=0;Rgp=0;
a=cos(2*pi/3)+sin(2*pi/3)*i;
Mmp=zeros(12);
%faltam=input('Tipo de 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 :');
switch(faltam)
case{1},
%Rfm=input('Resistencia de Falta en m. Rf(Ohm) = '); Rfm=Rfm/zbase;
Mmp(1,1)=1;Mmp(1,2)=-1;
Mmp(2,1)=1;Mmp(2,3)=-1;
Mmp(3,1)=-3*Rfm;Mmp(3,4)=1;Mmp(3,5)=1;Mmp(3,6)=1;
%fprintf('\t\tFalta:ag');
%fprintf('\tRfm:%-10.2f\n',Rfm);
case{2},
%Rfm=input('Resistencia de Falta en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mmp(1,1)=1;Mmp(1,2)=-a^2;
Mmp(2,1)=1;Mmp(2,3)=-a;
Mmp(3,1)=-Rfm;Mmp(3,2)=-a^2*Rfm;Mmp(3,3)=-a*Rfm;Mmp(3,4)=1;Mmp(3,5)=a^2;Mmp(3,6)=a;
%fprintf('\t\tFalta:bg');
%fprintf('\tRfm:%-10.2f\n',Rfm);
case{3},
%Rfm=input('Resistencia de Falta en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mmp(1,1)=1;Mmp(1,2)=-a;
Mmp(2,1)=1;Mmp(2,3)=-a^2;
Mmp(3,1)=-Rfm;Mmp(3,2)=-a*Rfm;Mmp(3,3)=-a^2*Rfm;Mmp(3,4)=1;Mmp(3,5)=a;Mmp(3,6)=a^2;
%fprintf('\t\tFalta:cg');
%fprintf('\tRfm:%-10.2f\n',Rfm);
case{4},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mmp(1,1)=1;
Mmp(2,2)=1;Mmp(2,3)=a;
Mmp(3,2)=-Rfm;Mmp(3,3)=a*Rfm;Mmp(3,5)=1;Mmp(3,6)=-a;
%fprintf('\t\tFalta:ab');
%fprintf('\tRfm:%-10.2f\n',Rfm);
case{5},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mmp(1,1)=1;
Mmp(2,2)=1;Mmp(2,3)=1;
Mmp(3,2)=-Rfm;Mmp(3,3)=Rfm;Mmp(3,5)=1;Mmp(3,6)=-1;
%fprintf('\t\tFalta:bc');
%fprintf('\tRfm:%-10.2f\n',Rfm);
case{6},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mmp(1,1)=1;
Mmp(2,2)=1;Mmp(2,3)=a^2;
Mmp(3,2)=-Rfm;Mmp(3,3)=a^2*Rfm;Mmp(3,5)=1;Mmp(3,6)=-a^2;
%fprintf('\t\tFalta:ca');
%fprintf('\tRfm:%-10.2f\n',Rfm);
case{7},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');Rfm=Rfm/zbase;
%Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): ');Rgm=Rgm/zbase;
Mmp(1,1)=1;Mmp(1,2)=a;Mmp(1,3)=a^2;
Mmp(2,2)=-a^2*Rfm;Mmp(2,3)=Rfm;Mmp(2,5)=a^2;Mmp(2,6)=-1;
Mmp(3,1)=-(Rfm+3*Rgm);Mmp(3,2)=-Rfm;Mmp(3,3)=-Rfm;
Mmp(3,4)=1;Mmp(3,5)=1;Mmp(3,6)=1;
%fprintf('\t\tFalta:abg');
%fprintf('\tRfm:%-10.2f',Rfm);
%fprintf('\tRgm:%-10.2f\n',Rgm);
case{8},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');Rfm=Rfm/zbase;
%Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): ');Rgm=Rgm/zbase;
Mmp(1,1)=1;Mmp(1,2)=1;Mmp(1,3)=1;
Mmp(2,2)=-Rfm;Mmp(2,3)=Rfm;Mmp(2,5)=1;Mmp(2,6)=-1;
Mmp(3,1)=-(Rfm+3*Rgm);Mmp(3,2)=-a^2*Rfm;Mmp(3,3)=-a*Rfm;
Mmp(3,4)=1;Mmp(3,5)=a^2;Mmp(3,6)=a;
%fprintf('\t\tFalta:bcg');
%fprintf('\tRfm:%-10.2f',Rfm);
%fprintf('\tRgm:%-10.2f\n',Rgm);
case{9},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');Rfm=Rfm/zbase;
%Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): ');Rgm=Rgm/zbase;
Mmp(1,1)=1;Mmp(1,2)=a^2;Mmp(1,3)=a;
Mmp(2,2)=-Rfm;Mmp(2,3)=a^2*Rfm;Mmp(2,5)=1;Mmp(2,6)=-a^2;
Mmp(3,1)=-(Rfm+3*Rgm);Mmp(3,2)=-a*Rfm;Mmp(3,3)=-a^2*Rfm;
Mmp(3,4)=1;Mmp(3,5)=a;Mmp(3,6)=a^2;
%fprintf('\t\tFalta:cag');
%fprintf('\tRfm:%-10.2f',Rfm);
%fprintf('\tRgm:%-10.2f\n',Rgm);
case{10},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm) = ');Rfm=Rfm/zbase;
Mmp(1,1)=1;
Mmp(2,2)=Rfm*a^2;Mmp(2,3)=-Rfm;Mmp(2,5)=-a^2;Mmp(2,6)=1;
Mmp(3,2)=-Rfm;Mmp(3,3)=Rfm*a^2;Mmp(3,5)=1;Mmp(3,6)=-a^2;
%fprintf('\t\tFalta:abc');
%fprintf('\tRfm:%-10.2f\n',Rfm);
case{11},
%Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm) = ');Rfm=Rfm/zbase;
%Rgm=input('Resistencia de Falta,comun en m. Rg(Ohm): ');Rgm=Rgm/zbase;
Mmp(1,1)=-(Rfm+3*Rgm);Mmp(1,4)=1;
Mmp(2,2)=-Rfm;Mmp(2,3)=-Rfm;Mmp(2,5)=1;Mmp(2,6)=1;
Mmp(3,2)=-a*Rfm;Mmp(3,3)=-Rfm;Mmp(3,5)=a;Mmp(3,6)=1;
%fprintf('\t\tFalta:abcg');
%fprintf('\tRfm:%-10.2f',Rfm);
%fprintf('\tRgm:%-10.2f\n',Rgm);
end
%faltap=input('Tipo de 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 :');
switch(faltap)
case{1},
%Rfp=input('Resistencia de Falta en p. Rf(Ohm): ');Rfp=Rfp/zbase;
Mmp(4,7)=1;Mmp(4,8)=-1;
Mmp(5,7)=1;Mmp(5,9)=-1;
Mmp(6,7)=-3*Rfp;Mmp(6,10)=1;Mmp(6,11)=1;Mmp(6,12)=1;
%fprintf('\t\tFalta:ag');
%fprintf('\tRfp:%-10.2f\n',Rfm);
case{2},
%Rfp=input('Resistencia de Falta en p. Rf(Ohm): ');Rfp=Rfp/zbase;
Mmp(4,7)=1;Mmp(4,8)=-a^2;
Mmp(5,7)=1;Mmp(5,9)=-a;
Mmp(6,7)=-Rfp;Mmp(6,8)=-a^2*Rfp;Mmp(6,9)=-a*Rfp;Mmp(6,10)=1;Mmp(6,11)=a^2;Mmp(6,12)=a;
%fprintf('\t\tFalta:bg');
%fprintf('\tRfp:%-10.2f\n',Rfm);
case{3},
%Rfp=input('Resistencia de Falta en p. Rf(Ohm): ');Rfp=Rfp/zbase;
Mmp(4,7)=1;Mmp(4,8)=-a;
Mmp(5,7)=1;Mmp(5,9)=-a^2;
Mmp(6,7)=-Rfp;Mmp(6,8)=-a*Rfp;Mmp(6,9)=-a^2*Rfp;Mmp(6,10)=1;Mmp(6,11)=a;Mmp(6,12)=a^2;
%fprintf('\t\tFalta:cg');
%fprintf('\tRfp:%-10.2f\n',Rfm);
case{4},
%Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): ');Rfp=Rfp/zbase;
Mmp(4,7)=1;
Mmp(5,8)=1;Mmp(5,9)=a;
Mmp(6,8)=-Rfp;Mmp(6,9)=a*Rfp;Mmp(6,11)=1;Mmp(6,12)=-a;
%fprintf('\t\tFalta:ab');
%fprintf('\tRfp:%-10.2f\n',Rfm);
case{5},
%Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): ');Rfp=Rfp/zbase;
Mmp(4,7)=1;
Mmp(5,8)=1;Mmp(5,9)=1;
Mmp(6,8)=-Rfp;Mmp(6,9)=Rfp;Mmp(6,11)=1;Mmp(6,12)=-1;
%fprintf('\t\tFalta:bc');
%fprintf('\tRfp:%-10.2f\n',Rfm);
case{6},
%Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): ');Rfp=Rfp/zbase;
Mmp(4,7)=1;
Mmp(5,8)=1;Mmp(5,9)=a^2;
Mmp(6,8)=-Rfp;Mmp(6,9)=a^2*Rfp;Mmp(6,11)=1;Mmp(6,12)=-a^2;
%fprintf('\t\tFalta:ca');
%fprintf('\tRfp:%-10.2f\n',Rfm);
case{7},
%Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): ');Rfp=Rfp/zbase;
%Rgp=input('Resistencia de Falta,comun en p. Rg(Ohm): ');Rgp=Rgp/zbase;
Mmp(4,7)=1;Mmp(4,8)=a;Mmp(4,9)=a^2;
Mmp(5,8)=-a^2*Rfp;Mmp(5,9)=Rfp;Mmp(5,11)=a^2;Mmp(5,12)=-1;
Mmp(6,7)=-(Rfp+3*Rgp);Mmp(6,8)=-Rfp;Mmp(6,9)=-Rfp;
Mmp(6,10)=1;Mmp(6,11)=1;Mmp(6,12)=1;
%fprintf('\t\tFalta:abg');
%fprintf('\tRfp:%-10.2f',Rfm);
%fprintf('\tRgp:%-10.2f\n',Rgm);
case{8},
%Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): ');Rfp=Rfp/zbase;
%Rgp=input('Resistencia de Falta,comun en p. Rg(Ohm): ');Rgp=Rgp/zbase;
Mmp(4,7)=1;Mmp(4,8)=1;Mmp(4,9)=1;
Mmp(5,8)=-Rfp;Mmp(5,9)=Rfp;Mmp(5,11)=1;Mmp(5,12)=-1;
Mmp(6,7)=-(Rfp+3*Rgp);Mmp(6,8)=-a^2*Rfp;Mmp(6,9)=-a*Rfp;
Mmp(6,10)=1;Mmp(6,11)=a^2;Mmp(6,12)=a;
%fprintf('\t\tFalta:bcg');
%fprintf('\tRfp:%-10.2f',Rfm);
%fprintf('\tRgp:%-10.2f\n',Rgm);
case{9},
%Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm): ');Rfp=Rfp/zbase;
%Rgp=input('Resistencia de Falta,comun en p. Rg(Ohm): ');Rgp=Rgp/zbase;
Mmp(4,7)=1;Mmp(4,8)=a^2;Mmp(4,9)=a;
Mmp(5,8)=-Rfp;Mmp(5,9)=a^2*Rfp;Mmp(5,11)=1;Mmp(5,12)=-a^2;
Mmp(6,7)=-(Rfp+3*Rgp);Mmp(6,8)=-a*Rfp;Mmp(6,9)=-a^2*Rfp;
Mmp(6,10)=1;Mmp(6,11)=a;Mmp(6,12)=a^2;
%fprintf('\t\tFalta:cag');
%fprintf('\tRfp:%-10.2f',Rfm);
%fprintf('\tRgp:%-10.2f\n',Rgm);
case{10},
%Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm) = ');Rfp=Rfp/zbase;
Mmp(4,7)=1;
Mmp(5,8)=Rfp*a^2;Mmp(5,9)=-Rfp;Mmp(5,11)=-a^2;Mmp(5,12)=1;
Mmp(6,8)=-Rfp;Mmp(6,9)=Rfp*a^2;Mmp(6,11)=1;Mmp(6,12)=-a^2;
%fprintf('\t\tFalta:abc');
%fprintf('\tRfp:%-10.2f',Rfm);
case{11},
%Rfp=input('Resistencia de Falta,por fase, en p. Rf(Ohm) = ');Rfp=Rfp/zbase;
%Rgp=input('Resistencia de Falta,comun en p. Rg(Ohm): ');Rgp=Rgp/zbase;
Mmp(4,7)=-(Rfp+3*Rgp);Mmp(4,10)=1;
Mmp(5,8)=-Rfp;Mmp(5,9)=-Rfp;Mmp(5,11)=1;Mmp(5,12)=1;
Mmp(6,8)=-a*Rfp;Mmp(6,9)=-Rfp;Mmp(6,11)=a;Mmp(6,12)=1;
%fprintf('\t\tFalta:abcg');
%fprintf('\tRfp:%-10.2f',Rfm);
%fprintf('\tRgp:%-10.2f\n',Rgm);
end
Mmp(7,1)=Z0(pfaltam,pfaltam); Mmp(7,4)=1; Mmp(7,7)=Z0(pfaltam,pfaltap);
Mmp(8,1)=Z0(pfaltap,pfaltam); Mmp(8,7)=Z0(pfaltap,pfaltap); Mmp(8,10)=1;
Mmp(9,2)=Z1(pfaltam,pfaltam); Mmp(9,5)=1; Mmp(9,8)=Z1(pfaltam,pfaltap);
Mmp(10,2)=Z1(pfaltap,pfaltam); Mmp(10,8)=Z1(pfaltap,pfaltap); Mmp(10,11)=1;
Mmp(11,3)=Z2(pfaltam,pfaltam); Mmp(11,6)=1; Mmp(11,9)=Z2(pfaltam,pfaltap);
Mmp(12,3)=Z2(pfaltap,pfaltam); Mmp(12,9)=Z2(pfaltap,pfaltap); 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.
