function[Im0,Im1,Im2,Vm0,Vm1,Vm2,faltam,Rfm,Rgm]=ma2m(Z0,Z1,Z2,E,zbase,pfaltam) % PFCc24
Rgm=0;
a=cos(2*pi/3)+sin(2*pi/3)*i;
Mm=zeros(6);
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}, % ag
Rfm=input('Resistencia de Falta en m. Rf(Ohm) = '); Rfm=Rfm/zbase;
Mm(1,1)=1;Mm(1,2)=-1;
Mm(2,1)=1;Mm(2,3)=-1;
Mm(3,1)=-3*Rfm;Mm(3,4)=1;Mm(3,5)=1;Mm(3,6)=1;
case{2}, % bg
Rfm=input('Resistencia de Falta en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mm(1,1)=1;Mm(1,2)=-a^2;
Mm(2,1)=1;Mm(2,3)=-a;
Mm(3,1)=-Rfm;Mm(3,2)=-a^2*Rfm;Mm(3,3)=-a*Rfm;Mm(3,4)=1;Mm(3,5)=a^2;Mm(3,6)=a;
case{3}, % cg
Rfm=input('Resistencia de Falta en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mm(1,1)=1;Mm(1,2)=-a;
Mm(2,1)=1;Mm(2,3)=-a^2;
Mm(3,1)=-Rfm;Mm(3,2)=-a*Rfm;Mm(3,3)=-a^2*Rfm;Mm(3,4)=1;Mm(3,5)=a;Mm(3,6)=a^2;
case{4}, % ab
Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mm(1,1)=1;
Mm(2,2)=1;Mm(2,3)=a;
Mm(3,2)=-Rfm;Mm(3,3)=a*Rfm;Mm(3,5)=1;Mm(3,6)=-a;
case{5}, % bc
Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mm(1,1)=1;
Mm(2,2)=1;Mm(2,3)=1;
Mm(3,2)=-Rfm;Mm(3,3)=Rfm;Mm(3,5)=1;Mm(3,6)=-1;
case{6}, % ca
Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm): ');Rfm=Rfm/zbase;
Mm(1,1)=1;
Mm(2,2)=1;Mm(2,3)=a^2;
Mm(3,2)=-Rfm;Mm(3,3)=a^2*Rfm;Mm(3,5)=1;Mm(3,6)=-a^2;
case{7}, % abg
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;
Mm(1,1)=1;Mm(1,2)=a;Mm(1,3)=a^2;
Mm(2,2)=-a^2*Rfm;Mm(2,3)=Rfm;Mm(2,5)=a^2;Mm(2,6)=-1;
Mm(3,1)=-(Rfm+3*Rgm);Mm(3,2)=-Rfm;Mm(3,3)=-Rfm;
Mm(3,4)=1;Mm(3,5)=1;Mm(3,6)=1;
case{8}, % bcg
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;
Mm(1,1)=1;Mm(1,2)=1;Mm(1,3)=1;
Mm(2,2)=-Rfm;Mm(2,3)=Rfm;Mm(2,5)=1;Mm(2,6)=-1;
Mm(3,1)=-(Rfm+3*Rgm);Mm(3,2)=-a^2*Rfm;Mm(3,3)=-a*Rfm;
Mm(3,4)=1;Mm(3,5)=a^2;Mm(3,6)=a;
case{9}, % cag
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;
Mm(1,1)=1;Mm(1,2)=a^2;Mm(1,3)=a;
Mm(2,2)=-Rfm;Mm(2,3)=a^2*Rfm;Mm(2,5)=1;Mm(2,6)=-a^2;
Mm(3,1)=-(Rfm+3*Rgm);Mm(3,2)=-a*Rfm;Mm(3,3)=-a^2*Rfm;
Mm(3,4)=1;Mm(3,5)=a;Mm(3,6)=a^2;
case{10}, % abc
Rfm=input('Resistencia de Falta,por fase, en m. Rf(Ohm) = ');Rfm=Rfm/zbase;
Mm(1,1)=1;
Mm(2,2)=-Rfm; Mm(2,5)=1;
Mm(3,3)=-Rfm; Mm(3,6)=1;
case{11}, % abcg
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;
Mm(1,1)=-(Rfm+3*Rgm); Mm(1,4)=1;
Mm(2,2)=-Rfm; Mm(2,5)=1;
Mm(3,3)=-Rfm; Mm(3,6)=1;
end
Mm(4,1)=Z0(pfaltam,pfaltam); Mm(4,4)=1;
Mm(5,2)=Z1(pfaltam,pfaltam); Mm(5,5)=1;
Mm(6,3)=Z2(pfaltam,pfaltam); Mm(6,6)=1;
P=[0;0;0;0;E;0];
N=Mm\P;
% -------------------------------------------------------
%
Im0=N(1); Im1=N(2); Im2=N(3);
Vm0=N(4); Vm1=N(5); Vm2=N(6);
e-REdING. Biblioteca de la Escuela Superior de Ingenieros de Sevilla.
