clc variasejecuciones=1; while variasejecuciones==1 clear % -------------------------------------------------- salida=input('Pantalla=1, Archivo=2 :'); if salida==2 cadena=input('Poner NombreArchivo.txt:','s'); cadena1=strcat('Salida\',cadena); end if salida==1 fid=1; elseif salida==2 fid=fopen(cadena1,'w+'); end % -------------------------------------------------- fprintf(fid,'PC44\n'); fprintf(fid,'Esquema Linea Doble SE\n'); fprintf(fid,'Lineas Separadas ER y ET\n'); fprintf(fid,'Unominal kV: 400,220,132,66\n'); fprintf(fid,'Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n'); fprintf(fid,'Compensacion 21: Si,No\n'); fprintf(fid,'Eleccion: kE,L1,L2,L3,Sccs,Sccr,Scct,Kr0,Ks0,k0=zt0/zs0\n'); fprintf(fid,'Falta Simple m, o Doble m,p\n'); fprintf(fid,'No situar m y p en el mismo Tramo\n'); fprintf(fid,'Barras Ficticias Iniciales m,p\n'); fprintf(fid,'\n'); % -------------------------------------------------- %PC44 datos=input('datos: standard=1,input=2: '); if datos==1 Un=220; kE=1; L1=200; L2=100; L3=10; Compens21=2; Sccs=1000; Sccr=1000; Scct=5000; Ks0=1; Kr0=1; % -------------------------------------------------- else Un=input('Un kV: 400,220,132,66 :'); kE=input('kE=E/Enominal: '); L1=input('L1 km: '); L2=input('L2 km: '); L3=input('L3 km: '); Compens21=input('Compensacion 21: Si=1, No=2 :'); Sccs=input('Sccs(MVA): '); Sccr=input('Sccr(MVA): '); Scct=input('Scct(MVA): '); Ks0=input('Ks0: '); Kr0=input('Kr0: '); end k0=input('k0=zt0/zs0: '); % ----------------------------------------------------------------- % DATOS % ----------------------------------------------------------------- Ubase=Un; % kV Vbase=Ubase/(sqrt(3)); % kV MVAbase=100; Ibase=MVAbase/(3*Vbase); % kA zbase=Vbase/Ibase; % Ohm if Un==400 zl1u=(0.0346+0.3208i)/zbase; zl2u=zl1u; zl0u=(0.1828+1.1482i)/zbase; zlM0u=(0.1482+0.7383i)/zbase; elseif Un==220 zl1u=(0.0829+0.4135i)/zbase; zl2u=zl1u; zl0u=(0.2311+1.3085i)/zbase; zlM0u=(0.1482+0.8031i)/zbase; elseif Un==132 zl1u=(0.1342+0.4042i)/zbase; zl2u=zl1u; zl0u=(0.2824+1.3718i)/zbase; zlM0u=(0.1482+0.8947i)/zbase; elseif Un==66 zl1u=(0.2156+0.3952i)/zbase; zl2u=zl1u; zl0u=(0.3636+1.4278i)/zbase; zlM0u=(0.1482+0.9288i)/zbase; end % --------------------------------------------------- a=0.5; b=0.5; c=0.5; d=0.5; % --------------------------------------------------- %CALCULOS PRELIMINARES % --------------------------------------------------- U=kE*Un; % kV E=U/sqrt(3)+0i; % kV E=E/Vbase; % E pu k=(zl0u-zl1u)/(3*zl1u); Km=zlM0u/zl1u; zs1=(MVAbase/Sccs)*i; % pu zr1=(MVAbase/Sccr)*i; % pu zt1=(MVAbase/Scct)*i; % pu zs2=zs1; % pu zr2=zr1; % pu zt2=zt1; % pu zs0=Ks0*zs1; % pu zr0=Kr0*zr1; % pu zt0=k0*zs0; % pu % ---------------------------------------------------- % CALCULOS % ----------------------------------------------------- TipoFalta=input('Falta Simple m=1, Doble m,p=2 :'); if fid==1 else fprintf(fid,'Tramo con Falta m: T11=6,T12=7,T21=8,T22=9\n'); end pfaltam=input('Tramo con Falta m: T11=6,T12=7,T21=8,T22=9:'); if pfaltam==6 fprintf(fid,'Falta m en Linea SR, entre S y E\n'); a=input('m =(0...1)='); m=a; if fid==1 else fprintf(fid,'m=%4.2f\n',m); end Sm=zl1u*m*L1; Rm=zl1u*(L2+(1-m)*L1); % pu Sm=Sm*zbase; Rm=Rm*zbase; % Ohm fprintf(fid,'Sm =%-5.2f+i(%-5.2f)Ohm\n',real(Sm),imag(Sm)); fprintf(fid,'Rm =%-5.2f+i(%-5.2f)Ohm\n',real(Rm),imag(Rm)); fprintf(fid,'\n'); elseif pfaltam==7 fprintf(fid,'Falta m en Linea SR, entre E y R\n'); b=input('m=(0...1)='); m=b; if fid==1 else fprintf(fid,'m=%4.2f\n',m); end Sm=zl1u*(L1+m*L2); Rm=zl1u*(1-m)*L2; Sm=Sm*zbase; Rm=Rm*zbase; fprintf(fid,'Sm =%-5.2f+i(%-5.2f)Ohm\n',real(Sm),imag(Sm)); fprintf(fid,'Rm =%-5.2f+i(%-5.2f)Ohm\n',real(Rm),imag(Rm)); fprintf(fid,'\n'); elseif pfaltam==8 fprintf(fid,'Falta m en Linea ST, entre S y E\n'); c=input('m=(0...1)='); m=c; if fid==1 else fprintf(fid,'m=%4.2f\n',m); end Sm=zl1u*m*L1; Tm=zl1u*(L3+(1-m)*L1); Sm=Sm*zbase; Tm=Tm*zbase; fprintf(fid,'Sm =%-5.2f+i(%-5.2f)Ohm\n',real(Sm),imag(Sm)); fprintf(fid,'Tm =%-5.2f+i(%-5.2f)Ohm\n',real(Tm),imag(Tm)); fprintf(fid,'\n'); elseif pfaltam==9 fprintf(fid,'Falta m en Linea ST, entre E y T\n'); d=input('m=(0...1)='); m=d; if fid==1 else fprintf(fid,'m=%4.2f\n',m); end Sm=zl1u*(L1+m*L3); Tm=zl1u*(1-m)*L3; Sm=Sm*zbase; Tm=Tm*zbase; fprintf(fid,'Sm =%-5.2f+i(%-5.2f)Ohm\n',real(Sm),imag(Sm)); fprintf(fid,'Tm =%-5.2f+i(%-5.2f)Ohm\n',real(Tm),imag(Tm)); fprintf(fid,'\n'); end if TipoFalta==2 if fid==1 else fprintf(fid,'Tramo con Falta p: T11=6,T12=7,T21=8,T22=9\n'); end pfaltap=input('Tramo con Falta p: T11=6,T12=7,T21=8,T22=9: '); if pfaltap==6 fprintf(fid,'Falta p en Linea SR, entre S y E\n'); a=input('p =(0...1)='); p=a; if fid==1 else fprintf(fid,'p=%4.2f\n',p); end Sp=zl1u*p*L1; Rp=zl1u*(L2+(1-p)*L1); % pu Sp=Sp*zbase; Rp=Rp*zbase; % Ohm fprintf(fid,'Sp =%-5.2f+i(%-5.2f)Ohm\n',real(Sp),imag(Sp)); fprintf(fid,'Rp =%-5.2f+i(%-5.2f)Ohm\n',real(Rp),imag(Rp)); fprintf(fid,'\n'); elseif pfaltap==7 fprintf(fid,'Falta p en Linea SR, entre E y R\n'); b=input('p=(0...1)='); p=b; if fid==1 else fprintf(fid,'p=%4.2f\n',m); end Sp=zl1u*(L1+p*L2); Rp=zl1u*(1-p)*L2; Sp=Sp*zbase; Rp=Rp*zbase; fprintf(fid,'Sp =%-5.2f+i(%-5.2f)Ohm\n',real(Sp),imag(Sp)); fprintf(fid,'Rp =%-5.2f+i(%-5.2f)Ohm\n',real(Rp),imag(Rp)); fprintf(fid,'\n'); elseif pfaltap==8 fprintf(fid,'Falta p en Linea ST, entre S y E\n'); c=input('p=(0...1)='); p=c; if fid==1 else fprintf(fid,'p=%4.2f\n',m); end Sp=zl1u*p*L1; Tp=zl1u*(L3+(1-p)*L1); Sp=Sp*zbase; Tp=Tp*zbase; fprintf(fid,'Sp =%-5.2f+i(%-5.2f)Ohm\n',real(Sp),imag(Sp)); fprintf(fid,'Tp =%-5.2f+i(%-5.2f)Ohm\n',real(Tp),imag(Tp)); fprintf(fid,'\n'); elseif pfaltap==9 fprintf(fid,'Falta p en Linea ST, entre E y T\n'); d=input('p=(0...1)='); p=d; if fid==1 else fprintf(fid,'p=%4.2f\n',m); end Sp=zl1u*(L1+p*L3); Tp=zl1u*(1-p)*L3; Sp=Sp*zbase; Tp=Tp*zbase; fprintf(fid,'Sp =%-5.2f+i(%-5.2f)Ohm\n',real(Sp),imag(Sp)); fprintf(fid,'Tp =%-5.2f+i(%-5.2f)Ohm\n',real(Tp),imag(Tp)); fprintf(fid,'\n'); end end [Minc]=subMinc(1); [Minca]=subMinca(1); [z0]=subz0(zs0,zr0,zt0,zl0u,zlM0u,L1,L2,L3,a,b,c,d); % pu [z1]=subz1(zs1,zr1,zt1,zl1u,L1,L2,L3,a,b,c,d); % pu [z2]=subz2(zs2,zr2,zt2,zl2u,L1,L2,L3,a,b,c,d); % pu % MATRIZ DE ADMITANCIA DE BARRAS [Y0]=subY0(Minc,z0); % pu [Y1]=subY1(Minc,z1); % pu [Y2]=subY2(Minc,z2); % pu % ZBUS(0,1,2) Z0=inv(Y0); % pu Z1=inv(Y1); % pu Z2=inv(Y2); % pu if TipoFalta==1 % CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN EL PUNTO m DE FALTA [Im0,Im1,Im2,Vm0,Vm1,Vm2,faltam,Rfm,Rgm]=ma2m(Z0,Z1,Z2,E,zbase,pfaltam); %pu elseif TipoFalta==2 % CALCULO DE LAS COMPONENTES SIMETRICAS DE V,I EN LOS PUNTOS m,p DE FALTA [Im0,Im1,Im2,Vm0,Vm1,Vm2,Ip0,Ip1,Ip2,Vp0,Vp1,Vp2,faltam,faltap,Rfm,Rgm,Rfp,Rgp]=... ma4mp(Z0,Z1,Z2,E,zbase,pfaltam,pfaltap); % pu end if faltam==1 cfaltam='ag';elseif faltam==2 cfaltam='bg';elseif faltam==3 cfaltam='cg'; elseif faltam==4 cfaltam='ab';elseif faltam==5 cfaltam='bc';elseif faltam==6 cfaltam='ca'; elseif faltam==7 cfaltam='abg';elseif faltam==8 cfaltam='bcg';elseif faltam==9 cfaltam='cag'; elseif faltam==10 cfaltam='abc';elseif faltam==11 cfaltam='abcg'; end if TipoFalta==2 if faltap==1 cfaltap='ag';elseif faltap==2 cfaltap='bg';elseif faltap==3 cfaltap='cg'; elseif faltap==4 cfaltap='ab';elseif faltap==5 cfaltap='bc';elseif faltap==6 cfaltap='ca'; elseif faltap==7 cfaltap='abg';elseif faltap==8 cfaltap='bcg';elseif faltap==9 cfaltap='cag'; elseif faltap==10 cfaltap='abc';elseif faltap==11 cfaltap='abcg'; end end fprintf(fid,'\n'); % IBUS(0,1,2) I0=zeros(9,1); I1=zeros(9,1); I2=zeros(9,1); I0(pfaltam,1)=-Im0; I1(pfaltam,1)=-Im1; I2(pfaltam,1)=-Im2; if TipoFalta==2 I0(pfaltap,1)=-Ip0; I1(pfaltap,1)=-Ip1; I2(pfaltap,1)=-Ip2; end % VBUS(0,1,2) V0=Z0*I0; V1=Z1*I1+E; V1a=[V1;E]; V2=Z2*I2; % I(0,1,2) DE LOS ELEMENTOS [Ie0]=subIe0(Minc,V0,z0); [Ie1]=subIe1(Minca,V1a,z1); [Ie2]=subIe2(Minc,V2,z2); % ASIGNACIONES V,I Vbase=Vbase*1000; % V Ibase=Ibase*1000; % A VVm0=Vm0*Vbase;VVm1=Vm1*Vbase;VVm2=Vm2*Vbase; if TipoFalta==2 VVp0=Vp0*Vbase;VVp1=Vp1*Vbase;VVp2=Vp2*Vbase; end Vs0=V0(1);Vs1=V1(1);Vs2=V2(1); VVs0=Vs0*Vbase;VVs1=Vs1*Vbase;VVs2=Vs2*Vbase; Vr0=V0(2);Vr1=V1(2);Vr2=V2(2); VVr0=Vr0*Vbase;VVr1=Vr1*Vbase;VVr2=Vr2*Vbase; Vt0=V0(3);Vt1=V1(3);Vt2=V2(3); VVt0=Vt0*Vbase;VVt1=Vt1*Vbase;VVt2=Vt2*Vbase; IIm0=Im0*Ibase;IIm1=Im1*Ibase;IIm2=Im2*Ibase; if TipoFalta==2 IIp0=Ip0*Ibase;IIp1=Ip1*Ibase;IIp2=Ip2*Ibase; end Is0=Ie0(1);Is1=Ie1(1);Is2=Ie2(1); IIs0=Is0*Ibase;IIs1=Is1*Ibase;IIs2=Is2*Ibase; Ir0=Ie0(2);Ir1=Ie1(2);Ir2=Ie2(2); IIr0=Ir0*Ibase;IIr1=Ir1*Ibase;IIr2=Ir2*Ibase; It0=Ie0(3);It1=Ie1(3);It2=Ie2(3); IIt0=It0*Ibase;IIt1=It1*Ibase;IIt2=It2*Ibase; Isr0=Ie0(4);Isr1=Ie1(4);Isr2=Ie2(4); IIsr0=Isr0*Ibase;IIsr1=Isr1*Ibase;IIsr2=Isr2*Ibase; Ist0=Ie0(5);Ist1=Ie1(5);Ist2=Ie2(5); IIst0=Ist0*Ibase;IIst1=Ist1*Ibase;IIst2=Ist2*Ibase; Irs0=Ie0(6);Irs1=Ie1(6);Irs2=Ie2(6); IIrs0=Irs0*Ibase;IIrs1=Irs1*Ibase;IIrs2=Irs2*Ibase; Its0=Ie0(7);Its1=Ie1(7);Its2=Ie2(7); IIts0=Its0*Ibase;IIts1=Its1*Ibase;IIts2=Its2*Ibase; [Ima,Imb,Imc]=abc(Im0,Im1,Im2); IIma=Ima*Ibase;IImb=Imb*Ibase;IImc=Imc*Ibase; if TipoFalta==2 [Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2); IIpa=Ipa*Ibase;IIpb=Ipb*Ibase;IIpc=Ipc*Ibase; end [Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2); VVma=Vma*Vbase;VVmb=Vmb*Vbase;VVmc=Vmc*Vbase; if TipoFalta==2 [Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2); VVpa=Vpa*Vbase;VVpb=Vpb*Vbase;VVpc=Vpc*Vbase; end [Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2); VVsa=Vsa*Vbase;VVsb=Vsb*Vbase;VVsc=Vsc*Vbase; [Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2); VVra=Vra*Vbase;VVrb=Vrb*Vbase;VVrc=Vrc*Vbase; [Vta,Vtb,Vtc]=abc(Vt0,Vt1,Vt2); VVta=Vta*Vbase;VVtb=Vtb*Vbase;VVtc=Vtc*Vbase; [Isa,Isb,Isc]=abc(Is0,Is1,Is2); IIsa=Isa*Ibase;IIsb=Isb*Ibase;IIsc=Isc*Ibase; [Ira,Irb,Irc]=abc(Ir0,Ir1,Ir2); IIra=Ira*Ibase;IIrb=Irb*Ibase;IIrc=Irc*Ibase; [Ita,Itb,Itc]=abc(It0,It1,It2); IIta=Ita*Ibase;IItb=Itb*Ibase;IItc=Itc*Ibase; [Isra,Isrb,Isrc]=abc(Isr0,Isr1,Isr2); IIsra=Isra*Ibase;IIsrb=Isrb*Ibase;IIsrc=Isrc*Ibase; [Ista,Istb,Istc]=abc(Ist0,Ist1,Ist2); IIsta=Ista*Ibase;IIstb=Istb*Ibase;IIstc=Istc*Ibase; [Irsa,Irsb,Irsc]=abc(Irs0,Irs1,Irs2); IIrsa=Irsa*Ibase;IIrsb=Irsb*Ibase;IIrsc=Irsc*Ibase; [Itsa,Itsb,Itsc]=abc(Its0,Its1,Its2); IItsa=Itsa*Ibase;IItsb=Itsb*Ibase;IItsc=Itsc*Ibase; % PROTECCION DE DISTANCIA (21) switch(Compens21) case{1}, % COMPENSADA [Sr_ZAB,Sr_ZBC,Sr_ZCA,Sr_ZAG,Sr_ZBG,Sr_ZCG]=... PDistComp(Vsa,Vsb,Vsc,Isra,Isrb,Isrc,Isr0,k,Km,Ist0); [St_ZAB,St_ZBC,St_ZCA,St_ZAG,St_ZBG,St_ZCG]=... PDistComp(Vsa,Vsb,Vsc,Ista,Istb,Istc,Ist0,k,Km,Isr0); %[Rs_ZAB,Rs_ZBC,Rs_ZCA,Rs_ZAG,Rs_ZBG,Rs_ZCG]=... % PDistComp(Vra,Vrb,Vrc,Irsa,Irsb,Irsc,Irs0,k,Km,Its0); %[Ts_ZAB,Ts_ZBC,Ts_ZCA,Ts_ZAG,Ts_ZBG,Ts_ZCG]=... % PDistComp(Vta,Vtb,Vtc,Itsa,Itsb,Itsc,Its0,k,Km,Irs0); [Rs_ZAB,Rs_ZBC,Rs_ZCA,Rs_ZAG,Rs_ZBG,Rs_ZCG]=... PDist(Vra,Vrb,Vrc,Irsa,Irsb,Irsc,Irs0,k); [Ts_ZAB,Ts_ZBC,Ts_ZCA,Ts_ZAG,Ts_ZBG,Ts_ZCG]=... PDist(Vta,Vtb,Vtc,Itsa,Itsb,Itsc,Its0,k); case{2}, % NO COMPENSADA [Sr_ZAB,Sr_ZBC,Sr_ZCA,Sr_ZAG,Sr_ZBG,Sr_ZCG]=... PDist(Vsa,Vsb,Vsc,Isra,Isrb,Isrc,Isr0,k); [St_ZAB,St_ZBC,St_ZCA,St_ZAG,St_ZBG,St_ZCG]=... PDist(Vsa,Vsb,Vsc,Ista,Istb,Istc,Ist0,k); [Rs_ZAB,Rs_ZBC,Rs_ZCA,Rs_ZAG,Rs_ZBG,Rs_ZCG]=... PDist(Vra,Vrb,Vrc,Irsa,Irsb,Irsc,Irs0,k); [Ts_ZAB,Ts_ZBC,Ts_ZCA,Ts_ZAG,Ts_ZBG,Ts_ZCG]=... PDist(Vta,Vtb,Vtc,Itsa,Itsb,Itsc,Its0,k); end % 67N [Vspol,Isrop,AngDirSr]=PDire(Vs0,Isr0); VVspol=Vspol*Vbase;IIsrop=Isrop*Ibase; [Vspol,Istop,AngDirSt]=PDire(Vs0,Ist0); IIstop=Istop*Ibase; [Vrpol,Irsop,AngDirRs]=PDire(Vr0,Irs0); VVrpol=Vrpol*Vbase;IIrsop=Irsop*Ibase; [Vtpol,Itsop,AngDirTs]=PDire(Vt0,Its0); VVtpol=Vtpol*Vbase;IItsop=Itsop*Ibase; % ****************************************************************** % *************************** SALIDA ****************************** % ****************************************************************** fprintf(fid,'\n'); fprintf(fid,'***************************************************\n'); fprintf(fid,'Un=%-4.0f(kV)\n',Un); fprintf(fid,'kE=E/En=%-4.0f\n',kE); fprintf(fid,'L1%-5.1f(km)\n',L1); fprintf(fid,'L2=%-5.1f(km)\n',L2); fprintf(fid,'L3=%-5.1f(km)\n',L3); fprintf(fid,'Ibase=%-8.2f(kA)\n',Ibase); fprintf(fid,'Tipo de Falta, Simple(1), Doble(2) =%-1.0f\n',TipoFalta); fprintf(fid,'Compensacion 21, Si(1), No(2) =%-1.0f\n',Compens21); fprintf(fid,'Sccs =%-8.0f(MVA)\n',Sccs); fprintf(fid,'Sccr =%-8.0f(MVA)\n',Sccr); fprintf(fid,'Scct =%-8.0f(MVA)\n',Scct); fprintf(fid,'Ks0 =%-4.1f\n',Ks0); fprintf(fid,'Kr0 =%-4.1f\n',Kr0); fprintf(fid,'k0=zt0/zs0=%-4.1f\n',k0); fprintf(fid,'\n'); fprintf(fid,'m =%-3.2f\t',m); fprintf(fid,'Falta:%s\t',cfaltam); fprintf(fid,'Rfm(Ohm):%-10.2f',Rfm*zbase); if (faltam==7)|(faltam==8)|(faltam==9)|(faltam==11) fprintf(fid,'\tRgm(Ohm):%-10.2f',Rgm*zbase); end if TipoFalta==2 fprintf(fid,'\n'); fprintf(fid,'p =%-3.2f\t',p); fprintf(fid,'Falta:%s\t',cfaltap); fprintf(fid,'Rfp(Ohm):%-10.2f',Rfp*zbase); if (faltap==7)|(faltap==8)|(faltap==9)|(faltap==11) fprintf(fid,'\tRgp(Ohm):%-10.2f',Rgp*zbase); end end fprintf(fid,'\n'); fprintf(fid,'***************************************************\n'); variassalidas=1; while variassalidas==1 fprintf(fid,'\n'); output=input('Salida: pu(cart)=1, pu(polar)=2, Ohm-Volt-Amp(cart)=3, Ohm-Volt-Amp(polar)=4 :'); if salida==2 fprintf(fid,'Salida: pu(cart)=1, pu(polar)=2, Ohm-Volt-Amp(cart)=3, Ohm-Volt-Amp(polar)=4 : %1.0f',output); fprintf(fid,'\n'); end %output=menu('Elegir Salida:','pu(cartesiana)','pu(polar)','Ohm-Volt-Amp(cartesiana)','Ohm-Volt-Amp(polar)'); switch(output) % OPCIONES DE SALIDA % -------------------------------- % 21 67N V(0,1,2) V(a,b,c) % I(0,1,2) I(a,b,c) % 1: pu(cartesiana) *** *** % 2: pu(polar) *** *** *** % 3: OVA(cartesiana) *** *** *** % 4: OVA(polar) *** *** *** *** % ********************************************************************** case{1}, % pu(cartesiana): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c) % ----------------- ----------------- % ********************************************************************** % ----------------------------------------------------------------------- % ** C1 ************** SALIDA V(0,1,2),I(0,1,2) pu cartesiana *********** % ----------------------------------------------------------------------- % Vm,Vp_2,Vs,Vr % Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr % ----------------------------------------------------------------------- fprintf(fid,'\n'); fprintf(fid,'SALIDA V(0,1,2),I(0,1,2) pu cartesiana\n\n'); fprintf(fid,'\n'); fprintf(fid,'Vm0 =%8.4f+i(%-8.4f)pu\t',real(Vm0),imag(Vm0)); fprintf(fid,'Vm1 =%8.4f+i(%-8.4f)pu\t',real(Vm1),imag(Vm1)); fprintf(fid,'Vm2 =%8.4f+i(%-8.4f)pu\t',real(Vm2),imag(Vm2)); fprintf(fid,'\n'); if TipoFalta==2 fprintf(fid,'Vp0 =%8.4f+i(%-8.4f)pu\t',real(Vp0),imag(Vp0)); fprintf(fid,'Vp1 =%8.4f+i(%-8.4f)pu\t',real(Vp1),imag(Vp1)); fprintf(fid,'Vp2 =%8.4f+i(%-8.4f)pu\t',real(Vp2),imag(Vp2)); fprintf(fid,'\n'); end fprintf(fid,'Vs0 =%8.4f+i(%-8.4f)pu\t',real(Vs0),imag(Vs0)); fprintf(fid,'Vs1 =%8.4f+i(%-8.4f)pu\t',real(Vs1),imag(Vs1)); fprintf(fid,'Vs2 =%8.4f+i(%-8.4f)pu\t',real(Vs2),imag(Vs2)); fprintf(fid,'\n'); fprintf(fid,'Vr0 =%8.4f+i(%-8.4f)pu\t',real(Vr0),imag(Vr0)); fprintf(fid,'Vr1 =%8.4f+i(%-8.4f)pu\t',real(Vr1),imag(Vr1)); fprintf(fid,'Vr2 =%8.4f+i(%-8.4f)pu\t',real(Vr2),imag(Vr2)); fprintf(fid,'\n'); fprintf(fid,'Vt0 =%8.4f+i(%-8.4f)pu\t',real(Vt0),imag(Vt0)); fprintf(fid,'Vt1 =%8.4f+i(%-8.4f)pu\t',real(Vt1),imag(Vt1)); fprintf(fid,'Vt2 =%8.4f+i(%-8.4f)pu\t',real(Vt2),imag(Vt2)); fprintf(fid,'\n\n'); % ---------------------------------------------------------- fprintf(fid,'Im0 =%8.4f+i(%-8.4f)pu\t',real(Im0),imag(Im0)); fprintf(fid,'Im1 =%8.4f+i(%-8.4f)pu\t',real(Im1),imag(Im1)); fprintf(fid,'Im2 =%8.4f+i(%-8.4f)pu\t',real(Im2),imag(Im2)); fprintf(fid,'\n'); if TipoFalta==2 fprintf(fid,'Ip0 =%8.4f+i(%-8.4f)pu\t',real(Ip0),imag(Ip0)); fprintf(fid,'Ip1 =%8.4f+i(%-8.4f)pu\t',real(Ip1),imag(Ip1)); fprintf(fid,'Ip2 =%8.4f+i(%-8.4f)pu\t',real(Ip2),imag(Ip2)); fprintf(fid,'\n'); end fprintf(fid,'Is0 =%8.4f+i(%-8.4f)pu\t',real(Is0),imag(Is0)); fprintf(fid,'Is1 =%8.4f+i(%-8.4f)pu\t',real(Is1),imag(Is1)); fprintf(fid,'Is2 =%8.4f+i(%-8.4f)pu\t',real(Is2),imag(Is2)); fprintf(fid,'\n'); fprintf(fid,'Ir0 =%8.4f+i(%-8.4f)pu\t',real(Ir0),imag(Ir0)); fprintf(fid,'Ir1 =%8.4f+i(%-8.4f)pu\t',real(Ir1),imag(Ir1)); fprintf(fid,'Ir2 =%8.4f+i(%-8.4f)pu\t',real(Ir2),imag(Ir2)); fprintf(fid,'\n'); fprintf(fid,'It0 =%8.4f+i(%-8.4f)pu\t',real(It0),imag(It0)); fprintf(fid,'It1 =%8.4f+i(%-8.4f)pu\t',real(It1),imag(It1)); fprintf(fid,'It2 =%8.4f+i(%-8.4f)pu\t',real(It2),imag(It2)); fprintf(fid,'\n'); fprintf(fid,'Isr0 =%8.4f+i(%-8.4f)pu\t',real(Isr0),imag(Isr0)); fprintf(fid,'Isr1 =%8.4f+i(%-8.4f)pu\t',real(Isr1),imag(Isr1)); fprintf(fid,'Isr2 =%8.4f+i(%-8.4f)pu\t',real(Isr2),imag(Isr2)); fprintf(fid,'\n'); fprintf(fid,'Ist0 =%8.4f+i(%-8.4f)pu\t',real(Ist0),imag(Ist0)); fprintf(fid,'Ist1 =%8.4f+i(%-8.4f)pu\t',real(Ist1),imag(Ist1)); fprintf(fid,'Ist2 =%8.4f+i(%-8.4f)pu\t',real(Ist2),imag(Ist2)); fprintf(fid,'\n'); fprintf(fid,'Irs0 =%8.4f+i(%-8.4f)pu\t',real(Irs0),imag(Irs0)); fprintf(fid,'Irs1 =%8.4f+i(%-8.4f)pu\t',real(Irs1),imag(Irs1)); fprintf(fid,'Irs2 =%8.4f+i(%-8.4f)pu\t',real(Irs2),imag(Irs2)); fprintf(fid,'\n'); fprintf(fid,'Its0 =%8.4f+i(%-8.4f)pu\t',real(Its0),imag(Its0)); fprintf(fid,'Its1 =%8.4f+i(%-8.4f)pu\t',real(Its1),imag(Its1)); fprintf(fid,'Its2 =%8.4f+i(%-8.4f)pu\t',real(Its2),imag(Its2)); fprintf(fid,'\n'); % ------------------------------------------------------------------ % ** C1 *********** SALIDA V(a,b,c),I(a,b,c) pu cartesiana ********* % ------------------------------------------------------------------ % Vm,Vp_2,Vs,Vr % Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr % ------------------------------------------------------------------ fprintf(fid,'\n'); fprintf(fid,'SALIDA V(a,b,c),I(a,b,c) pu cartesiana\n\n'); fprintf(fid,'\n'); fprintf(fid,'Vma =%8.4f+i(%-8.4f)pu\t',real(Vma),imag(Vma)); fprintf(fid,'Vmb =%8.4f+i(%-8.4f)pu\t',real(Vmb),imag(Vmb)); fprintf(fid,'Vmc =%8.4f+i(%-8.4f)pu\t',real(Vmc),imag(Vmc)); if TipoFalta==2 fprintf(fid,'\n'); fprintf(fid,'Vpa =%8.4f+i(%-8.4f)pu\t',real(Vpa),imag(Vpa)); fprintf(fid,'Vpb =%8.4f+i(%-8.4f)pu\t',real(Vpb),imag(Vpb)); fprintf(fid,'Vpc =%8.4f+i(%-8.4f)pu\t',real(Vpc),imag(Vpc)); end fprintf(fid,'\n'); fprintf(fid,'Vsa =%8.4f+i(%-8.4f)pu\t',real(Vsa),imag(Vsa)); fprintf(fid,'Vsb =%8.4f+i(%-8.4f)pu\t',real(Vsb),imag(Vsb)); fprintf(fid,'Vsc =%8.4f+i(%-8.4f)pu\t',real(Vsc),imag(Vsc)); fprintf(fid,'\n'); fprintf(fid,'Vra =%8.4f+i(%-8.4f)pu\t',real(Vra),imag(Vra)); fprintf(fid,'Vrb =%8.4f+i(%-8.4f)pu\t',real(Vrb),imag(Vrb)); fprintf(fid,'Vrc =%8.4f+i(%-8.4f)pu\t',real(Vrc),imag(Vrc)); fprintf(fid,'\n'); fprintf(fid,'Vta =%8.4f+i(%-8.4f)pu\t',real(Vta),imag(Vta)); fprintf(fid,'Vtb =%8.4f+i(%-8.4f)pu\t',real(Vtb),imag(Vtb)); fprintf(fid,'Vtc =%8.4f+i(%-8.4f)pu\t',real(Vtc),imag(Vtc)); fprintf(fid,'\n\n'); fprintf(fid,'Ima =%8.4f+i(%-8.4f)pu\t',real(Ima),imag(Ima)); fprintf(fid,'Imb =%8.4f+i(%-8.4f)pu\t',real(Imb),imag(Imb)); fprintf(fid,'Imc =%8.4f+i(%-8.4f)pu\t',real(Imc),imag(Imc)); if TipoFalta==2 fprintf(fid,'\n'); fprintf(fid,'Ipa =%8.4f+i(%-8.4f)pu\t',real(Ipa),imag(Ipa)); fprintf(fid,'Ipb =%8.4f+i(%-8.4f)pu\t',real(Ipb),imag(Ipb)); fprintf(fid,'Ipc =%8.4f+i(%-8.4f)pu\t',real(Ipc),imag(Ipc)); end fprintf(fid,'\n'); fprintf(fid,'Isa =%8.4f+i(%-8.4f)pu\t',real(Isa),imag(Isa)); fprintf(fid,'Isb =%8.4f+i(%-8.4f)pu\t',real(Isb),imag(Isb)); fprintf(fid,'Isc =%8.4f+i(%-8.4f)pu\t',real(Isc),imag(Isc)); fprintf(fid,'\n'); fprintf(fid,'Ira =%8.4f+i(%-8.4f)pu\t',real(Ira),imag(Ira)); fprintf(fid,'Irb =%8.4f+i(%-8.4f)pu\t',real(Irb),imag(Irb)); fprintf(fid,'Irc =%8.4f+i(%-8.4f)pu\t',real(Irc),imag(Irc)); fprintf(fid,'\n'); fprintf(fid,'Ita =%8.4f+i(%-8.4f)pu\t',real(Ita),imag(Ita)); fprintf(fid,'Itb =%8.4f+i(%-8.4f)pu\t',real(Itb),imag(Itb)); fprintf(fid,'Itc =%8.4f+i(%-8.4f)pu\t',real(Itc),imag(Itc)); fprintf(fid,'\n'); fprintf(fid,'Isra =%8.4f+i(%-8.4f)pu\t',real(Isra),imag(Isra)); fprintf(fid,'Isrb =%8.4f+i(%-8.4f)pu\t',real(Isrb),imag(Isrb)); fprintf(fid,'Isrc =%8.4f+i(%-8.4f)pu\t',real(Isrc),imag(Isrc)); fprintf(fid,'\n'); fprintf(fid,'Ista =%8.4f+i(%-8.4f)pu\t',real(Ista),imag(Ista)); fprintf(fid,'Istb =%8.4f+i(%-8.4f)pu\t',real(Istb),imag(Istb)); fprintf(fid,'Istc =%8.4f+i(%-8.4f)pu\t',real(Istc),imag(Istc)); fprintf(fid,'\n'); fprintf(fid,'Irsa =%8.4f+i(%-8.4f)pu\t',real(Irsa),imag(Irsa)); fprintf(fid,'Irsb =%8.4f+i(%-8.4f)pu\t',real(Irsb),imag(Irsb)); fprintf(fid,'Irsc =%8.4f+i(%-8.4f)pu\t',real(Irsc),imag(Irsc)); fprintf(fid,'\n'); fprintf(fid,'Itsa =%8.4f+i(%-8.4f)pu\t',real(Itsa),imag(Itsa)); fprintf(fid,'Itsb =%8.4f+i(%-8.4f)pu\t',real(Itsb),imag(Itsb)); fprintf(fid,'Itsc =%8.4f+i(%-8.4f)pu\t',real(Itsc),imag(Itsc)); fprintf(fid,'\n'); % ****************************************************************** case{2}, % pu(polar): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c) % --- ----------------- ---------------- % ****************************************************************** % ------------------------------------------------------------------ % ** C2 ******************* SALIDA 67N pu polar ******************** % ------------------------------------------------------------------ % Sm(e1,e3),Sp_2,Sr_1,Rm(e1,e2),Rp_2,Rs_1 % ------------------------------------------------------------------ fprintf(fid,'SALIDA 67N, pu, polar\n'); fprintf(fid,'\n'); fprintf(fid,'Sr_67N:\tVspol=%8.4f\tIsrop=%8.4f\tAngDirSr=%5.1f\n',Vspol,Isrop,AngDirSr); fprintf(fid,'St_67N:\tVspol=%8.4f\tIstop=%8.4f\tAngDirSt=%5.1f\n',Vspol,Istop,AngDirSt); fprintf(fid,'Rs_67N:\tVrpol=%8.4f\tIrsop=%8.4f\tAngDirRs=%5.1f\n',Vrpol,Irsop,AngDirRs); fprintf(fid,'Ts_67N:\tVtpol=%8.4f\tItsop=%8.4f\tAngDirTs=%5.1f\n',Vtpol,Itsop,AngDirTs); fprintf(fid,'\n'); % ------------------------------------------------------------------ % ** C2 ************* SALIDA V(0,1,2), I(0,1,2) pu polar *********** % ------------------------------------------------------------------ % Vm,Vp_2,Vs,Vr % Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr % ------------------------------------------------------------------ fprintf(fid,'SALIDA V(0,1,2), I(0,1,2) pu polar\n'); fprintf(fid,'\n'); [mod,arg]=mod_arg(Vm0);fprintf(fid,'Vm0 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vm1);fprintf(fid,'Vm1 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vm2);fprintf(fid,'Vm2 =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); if TipoFalta==2 [mod,arg]=mod_arg(Vp0);fprintf(fid,'Vp0 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vp1);fprintf(fid,'Vp1 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vp2);fprintf(fid,'Vp2 =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); end [mod,arg]=mod_arg(Vs0);fprintf(fid,'Vs0 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vs1);fprintf(fid,'Vs1 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vs2);fprintf(fid,'Vs2 =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Vr0);fprintf(fid,'Vr0 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vr1);fprintf(fid,'Vr1 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vr2);fprintf(fid,'Vr2 =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Vt0);fprintf(fid,'Vt0 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vt1);fprintf(fid,'Vt1 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vt2);fprintf(fid,'Vt2 =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); %------------------------------------------------------------------------------- fprintf(fid,'\n'); %------------------------------------------------------------------------------- [mod,arg]=mod_arg(Im0);fprintf(fid,'Im0 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Im1);fprintf(fid,'Im1 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Im2);fprintf(fid,'Im2 =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); if TipoFalta==2 [mod,arg]=mod_arg(Ip0);fprintf(fid,'Ip0 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Ip1);fprintf(fid,'Ip1 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Ip2);fprintf(fid,'Ip2 =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); end [mod,arg]=mod_arg(Is0);fprintf(fid,'Is0 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Is1);fprintf(fid,'Is1 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Is2);fprintf(fid,'Is2 =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Ir0);fprintf(fid,'Ir0 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Ir1);fprintf(fid,'Ir1 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Ir2);fprintf(fid,'Ir2 =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(It0);fprintf(fid,'It0 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(It1);fprintf(fid,'It1 =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(It2);fprintf(fid,'It2 =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Isr0);fprintf(fid,'Isr0 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Isr1);fprintf(fid,'Isr1 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Isr2);fprintf(fid,'Isr2 =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Ist0);fprintf(fid,'Ist0 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Ist1);fprintf(fid,'Ist1 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Ist2);fprintf(fid,'Ist2 =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Irs0);fprintf(fid,'Irs0 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Irs1);fprintf(fid,'Irs1 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Irs2);fprintf(fid,'Irs2 =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Its0);fprintf(fid,'Its0 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Its1);fprintf(fid,'Its1 =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Its2);fprintf(fid,'Its2 =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); fprintf(fid,'\n'); % ------------------------------------------------------------------ % ** C2 ************ SALIDA V(a,b,c), I(a,b,c) pu Polar ************ % ------------------------------------------------------------------ % Vm,Vp_2,Vs,Vr % Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr % ------------------------------------------------------------------ fprintf(fid,'SALIDA V(a,b,c), I(a,b,c) pu Polar\n'); fprintf(fid,'\n'); [mod,arg]=mod_arg(Vma);fprintf(fid,'Vma =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vmb);fprintf(fid,'Vmb =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vmc);fprintf(fid,'Vmc =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); if TipoFalta==2 [mod,arg]=mod_arg(Vpa);fprintf(fid,'Vpa =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vpb);fprintf(fid,'Vpb =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vpc);fprintf(fid,'Vpc =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); end [mod,arg]=mod_arg(Vsa);fprintf(fid,'Vsa =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vsb);fprintf(fid,'Vsb =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vsc);fprintf(fid,'Vsc =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Vra);fprintf(fid,'Vra =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vrb);fprintf(fid,'Vrb =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vrc);fprintf(fid,'Vrc =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Vta);fprintf(fid,'Vta =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vtb);fprintf(fid,'Vtb =%8.4f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(Vtc);fprintf(fid,'Vtc =%8.4f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); %------------------------------------------------------------------------------- fprintf(fid,'\n'); %------------------------------------------------------------------------------- [mod,arg]=mod_arg(Ima);fprintf(fid,'Ima =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Imb);fprintf(fid,'Imb =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Imc);fprintf(fid,'Imc =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); if TipoFalta==2 [mod,arg]=mod_arg(Ipa);fprintf(fid,'Ipa =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Ipb);fprintf(fid,'Ipb =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Ipc);fprintf(fid,'Ipc =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); end [mod,arg]=mod_arg(Isa);fprintf(fid,'Isa =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Isb);fprintf(fid,'Isb =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Isc);fprintf(fid,'Isc =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Ira);fprintf(fid,'Ira =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Irb);fprintf(fid,'Irb =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Irc);fprintf(fid,'Irc =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Ita);fprintf(fid,'Ita =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Itb);fprintf(fid,'Itb =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Itc);fprintf(fid,'Itc =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n\n'); [mod,arg]=mod_arg(Isra);fprintf(fid,'Isra =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Isrb);fprintf(fid,'Isrb =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Isrc);fprintf(fid,'Isrc =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Ista);fprintf(fid,'Ista =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Istb);fprintf(fid,'Istb =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Istc);fprintf(fid,'Istc =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Irsa);fprintf(fid,'Irsa =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Irsb);fprintf(fid,'Irsb =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Irsc);fprintf(fid,'Irsc =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(Itsa);fprintf(fid,'Itsa =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Itsb);fprintf(fid,'Itsb =%8.4f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(Itsc);fprintf(fid,'Itsc =%8.4f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); % *********************************************************************** case{3}, % OVA cartesiana: 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c) % --- ----------------- ----------------- % *********************************************************************** % ----------------------------------------------------------------------- % C3.1 SALIDA 21, Ohm, cartesiana % ----------------------------------------------------------------------- % Sm(e1,e3),Sp_2,Sr_1,Rm(e1,e2),Rp_2,Rs_1 % ----------------------------------------------------------------------- fprintf(fid,'\n'); fprintf(fid,'SALIDA 21,Ohm, cartesiana \n'); fprintf(fid,'\n'); fprintf(fid,'Nota: Componentes de Z >=10000 Ohm se representan por Inf\n'); fprintf(fid,'Nota: Componentes de Z <=-10000 Ohm se representan por -Inf\n'); fprintf(fid,'\n'); if pfaltam==6 | pfaltam==7 [mod,arg]=mod_arg(Sm);fprintf(fid,'Sm=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_arg(Rm);fprintf(fid,'Rm=%7.2f(%-5.1f)(Ohm)\n',mod,arg); elseif pfaltam==8 | pfaltam==9 [mod,arg]=mod_arg(Sm);fprintf(fid,'Sm=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_arg(Tm);fprintf(fid,'Tm=%7.2f(%-5.1f)(Ohm)\n',mod,arg); end if TipoFalta==2 if pfaltap==6 | pfaltap==7 [mod,arg]=mod_arg(Sp);fprintf(fid,'Sp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_arg(Rp);fprintf(fid,'Rp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); elseif pfaltam==8 | pfaltam==9 [mod,arg]=mod_arg(Sp);fprintf(fid,'Sp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_arg(Tp);fprintf(fid,'Tp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); end end fprintf(fid,'\n'); BarraInt=input('Supone una Barra Intermedia A, en Linea 1, un poco antes de m?: Si=1,No=2 : '); if BarraInt==1 SApu=input('SA/SR = '); if salida ==2 fprintf(fid,'Suponer Barra Intermedia A, Linea 1: SA/SR =%4.2f\n',SApu); end SA=SApu*zl1*zbase; fprintf(fid,'SA =%-7.2f+i(%-7.2f)Ohm\n',real(SA),imag(SA)); SAzona1=0.85*SA; fprintf(fid,'SAL1(zona1) =%-7.2f+i(%-7.2f)Ohm\n',real(SAzona1),imag(SAzona1)); fprintf(fid,'\n'); end [Sr_ZABre,Sr_ZABim]=zmed(Sr_ZAB,zbase); fprintf(fid,'Sr_ZAB=%7.2f+i(%-7.2f)Ohm\n',Sr_ZABre,Sr_ZABim); [Sr_ZBCre,Sr_ZBCim]=zmed(Sr_ZBC,zbase); fprintf(fid,'Sr_ZBC=%7.2f+i(%-7.2f)Ohm\n',Sr_ZBCre,Sr_ZBCim); [Sr_ZCAre,Sr_ZCAim]=zmed(Sr_ZCA,zbase); fprintf(fid,'Sr_ZCA=%7.2f+i(%-7.2f)Ohm\n',Sr_ZCAre,Sr_ZCAim); [Sr_ZAGre,Sr_ZAGim]=zmed(Sr_ZAG,zbase); fprintf(fid,'Sr_ZAG=%7.2f+i(%-7.2f)Ohm\n',Sr_ZAGre,Sr_ZAGim); [Sr_ZBGre,Sr_ZBGim]=zmed(Sr_ZBG,zbase); fprintf(fid,'Sr_ZBG=%7.2f+i(%-7.2f)Ohm\n',Sr_ZBGre,Sr_ZBGim); [Sr_ZCGre,Sr_ZCGim]=zmed(Sr_ZCG,zbase); fprintf(fid,'Sr_ZCG=%7.2f+i(%-7.2f)Ohm\n',Sr_ZCGre,Sr_ZCGim); fprintf(fid,'\n'); [St_ZABre,St_ZABim]=zmed(St_ZAB,zbase);fprintf(fid,'St_ZAB=%7.2f+i(%-7.2f)Ohm\n',St_ZABre,St_ZABim); [St_ZBCre,St_ZBCim]=zmed(St_ZBC,zbase);fprintf(fid,'St_ZBC=%7.2f+i(%-7.2f)Ohm\n',St_ZBCre,St_ZBCim); [St_ZCAre,St_ZCAim]=zmed(St_ZCA,zbase);fprintf(fid,'St_ZCA=%7.2f+i(%-7.2f)Ohm\n',St_ZCAre,St_ZCAim); [St_ZAGre,St_ZAGim]=zmed(St_ZAG,zbase);fprintf(fid,'St_ZAG=%7.2f+i(%-7.2f)Ohm\n',St_ZAGre,St_ZAGim); [St_ZBGre,St_ZBGim]=zmed(St_ZBG,zbase);fprintf(fid,'St_ZBG=%7.2f+i(%-7.2f)Ohm\n',St_ZBGre,St_ZBGim); [St_ZCGre,St_ZCGim]=zmed(St_ZCG,zbase);fprintf(fid,'St_ZCG=%7.2f+i(%-7.2f)Ohm\n',St_ZCGre,St_ZCGim); fprintf(fid,'\n'); [Rs_ZABre,Rs_ZABim]=zmed(Rs_ZAB,zbase);fprintf(fid,'Rs_ZAB=%7.2f+i(%-7.2f)Ohm\n',Rs_ZABre,Rs_ZABim); [Rs_ZBCre,Rs_ZBCim]=zmed(Rs_ZBC,zbase);fprintf(fid,'Rs_ZBC=%7.2f+i(%-7.2f)Ohm\n',Rs_ZBCre,Rs_ZBCim); [Rs_ZCAre,Rs_ZCAim]=zmed(Rs_ZCA,zbase);fprintf(fid,'Rs_ZCA=%7.2f+i(%-7.2f)Ohm\n',Rs_ZCAre,Rs_ZCAim); [Rs_ZAGre,Rs_ZAGim]=zmed(Rs_ZAG,zbase);fprintf(fid,'Rs_ZAG=%7.2f+i(%-7.2f)Ohm\n',Rs_ZAGre,Rs_ZAGim); [Rs_ZBGre,Rs_ZBGim]=zmed(Rs_ZBG,zbase);fprintf(fid,'Rs_ZBG=%7.2f+i(%-7.2f)Ohm\n',Rs_ZBGre,Rs_ZBGim); [Rs_ZCGre,Rs_ZCGim]=zmed(Rs_ZCG,zbase);fprintf(fid,'Rs_ZCG=%7.2f+i(%-7.2f)Ohm\n',Rs_ZCGre,Rs_ZCGim); fprintf(fid,'\n'); [Ts_ZABre,Ts_ZABim]=zmed(Ts_ZAB,zbase);fprintf(fid,'Ts_ZAB=%7.2f+i(%-7.2f)Ohm\n',Ts_ZABre,Ts_ZABim); [Ts_ZBCre,Ts_ZBCim]=zmed(Ts_ZBC,zbase);fprintf(fid,'Ts_ZBC=%7.2f+i(%-7.2f)Ohm\n',Ts_ZBCre,Ts_ZBCim); [Ts_ZCAre,Ts_ZCAim]=zmed(Ts_ZCA,zbase);fprintf(fid,'Ts_ZCA=%7.2f+i(%-7.2f)Ohm\n',Ts_ZCAre,Ts_ZCAim); [Ts_ZAGre,Ts_ZAGim]=zmed(Ts_ZAG,zbase);fprintf(fid,'Ts_ZAG=%7.2f+i(%-7.2f)Ohm\n',Ts_ZAGre,Ts_ZAGim); [Ts_ZBGre,Ts_ZBGim]=zmed(Ts_ZBG,zbase);fprintf(fid,'Ts_ZBG=%7.2f+i(%-7.2f)Ohm\n',Ts_ZBGre,Ts_ZBGim); [Ts_ZCGre,Ts_ZCGim]=zmed(Ts_ZCG,zbase);fprintf(fid,'Ts_ZCG=%7.2f+i(%-7.2f)Ohm\n',Ts_ZCGre,Ts_ZCGim); fprintf(fid,'\n\n'); % ------------------------------------------------------------------ % C3.3 SALIDA V(0,1,2),I(0,1,2) Volt-Amp cartesiana % ------------------------------------------------------------------ % Vm,Vp_2,Vs,Vr % Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr % ------------------------------------------------------------------ fprintf(fid,'SALIDA V(0,1,2),I(0,1,2) Volt, Amp, cartesiana \n'); fprintf(fid,'\n'); fprintf(fid,'Vm0 =%8.0f+i(%-8.0f)V\t',real(VVm0),imag(VVm0)); fprintf(fid,'Vm1 =%8.0f+i(%-8.0f)V\t',real(VVm1),imag(VVm1)); fprintf(fid,'Vm2 =%8.0f+i(%-8.0f)V\t',real(VVm2),imag(VVm2)); if TipoFalta==2 fprintf(fid,'\n'); fprintf(fid,'Vp0 =%8.0f+i(%-8.0f)V\t',real(VVp0),imag(VVp0)); fprintf(fid,'Vp1 =%8.0f+i(%-8.0f)V\t',real(VVp1),imag(VVp1)); fprintf(fid,'Vp2 =%8.0f+i(%-8.0f)V\t',real(VVp2),imag(VVp2)); end fprintf(fid,'\n'); fprintf(fid,'Vs0 =%8.0f+i(%-8.0f)V\t',real(VVs0),imag(VVs0)); fprintf(fid,'Vs1 =%8.0f+i(%-8.0f)V\t',real(VVs1),imag(VVs1)); fprintf(fid,'Vs2 =%8.0f+i(%-8.0f)V\t',real(VVs2),imag(VVs2)); fprintf(fid,'\n'); fprintf(fid,'Vr0 =%8.0f+i(%-8.0f)V\t',real(VVr0),imag(VVr0)); fprintf(fid,'Vr1 =%8.0f+i(%-8.0f)V\t',real(VVr1),imag(VVr1)); fprintf(fid,'Vr2 =%8.0f+i(%-8.0f)V\t',real(VVr2),imag(VVr2)); fprintf(fid,'\n'); fprintf(fid,'Vt0 =%8.0f+i(%-8.0f)V\t',real(VVt0),imag(VVt0)); fprintf(fid,'Vt1 =%8.0f+i(%-8.0f)V\t',real(VVt1),imag(VVt1)); fprintf(fid,'Vt2 =%8.0f+i(%-8.0f)V\t',real(VVt2),imag(VVt2)); fprintf(fid,'\n\n'); % --------------------------------------------------------------- fprintf(fid,'Im0 =%8.0f+i(%-8.0f)A\t',real(IIm0),imag(IIm0)); fprintf(fid,'Im1 =%8.0f+i(%-8.0f)A\t',real(IIm1),imag(IIm1)); fprintf(fid,'Im2 =%8.0f+i(%-8.0f)A\t',real(IIm2),imag(IIm2)); if TipoFalta==2 fprintf(fid,'\n'); fprintf(fid,'Ip0 =%8.0f+i(%-8.0f)A\t',real(IIp0),imag(IIp0)); fprintf(fid,'Ip1 =%8.0f+i(%-8.0f)A\t',real(IIp1),imag(IIp1)); fprintf(fid,'Ip2 =%8.0f+i(%-8.0f)A\t',real(IIp2),imag(IIp2)); end fprintf(fid,'\n'); fprintf(fid,'Is0 =%8.0f+i(%-8.0f)A\t',real(IIs0),imag(IIs0)); fprintf(fid,'Is1 =%8.0f+i(%-8.0f)A\t',real(IIs1),imag(IIs1)); fprintf(fid,'Is2 =%8.0f+i(%-8.0f)A\t',real(IIs2),imag(IIs2)); fprintf(fid,'\n'); fprintf(fid,'Ir0 =%8.0f+i(%-8.0f)A\t',real(IIr0),imag(IIr0)); fprintf(fid,'Ir1 =%8.0f+i(%-8.0f)A\t',real(IIr1),imag(IIr1)); fprintf(fid,'Ir2 =%8.0f+i(%-8.0f)A\t',real(IIr2),imag(IIr2)); fprintf(fid,'\n'); fprintf(fid,'It0 =%8.0f+i(%-8.0f)A\t',real(IIt0),imag(IIt0)); fprintf(fid,'It1 =%8.0f+i(%-8.0f)A\t',real(IIt1),imag(IIt1)); fprintf(fid,'It2 =%8.0f+i(%-8.0f)A\t',real(IIt2),imag(IIt2)); fprintf(fid,'\n'); fprintf(fid,'Isr0 =%8.0f+i(%-8.0f)A\t',real(IIsr0),imag(IIsr0)); fprintf(fid,'Isr1 =%8.0f+i(%-8.0f)A\t',real(IIsr1),imag(IIsr1)); fprintf(fid,'Isr2 =%8.0f+i(%-8.0f)A\t',real(IIsr2),imag(IIsr2)); fprintf(fid,'\n'); fprintf(fid,'Ist0 =%8.0f+i(%-8.0f)A\t',real(IIst0),imag(IIst0)); fprintf(fid,'Ist1 =%8.0f+i(%-8.0f)A\t',real(IIst1),imag(IIst1)); fprintf(fid,'Ist2 =%8.0f+i(%-8.0f)A\t',real(IIst2),imag(IIst2)); fprintf(fid,'\n'); fprintf(fid,'Irs0 =%8.0f+i(%-8.0f)A\t',real(IIrs0),imag(IIrs0)); fprintf(fid,'Irs1 =%8.0f+i(%-8.0f)A\t',real(IIrs1),imag(IIrs1)); fprintf(fid,'Irs2 =%8.0f+i(%-8.0f)A\t',real(IIrs2),imag(IIrs2)); fprintf(fid,'\n'); fprintf(fid,'Its0 =%8.0f+i(%-8.0f)A\t',real(IIts0),imag(IIts0)); fprintf(fid,'Its1 =%8.0f+i(%-8.0f)A\t',real(IIts1),imag(IIts1)); fprintf(fid,'Its2 =%8.0f+i(%-8.0f)A\t',real(IIts2),imag(IIts2)); fprintf(fid,'\n'); % ------------------------------------------------------------------ % C3.4 SALIDA V(a,b,c),I(a,b,c) Volt,Amp cartesiana % ------------------------------------------------------------------ % Vm,Vp_2,Vs,Vr % Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr % ------------------------------------------------------------------ fprintf(fid,'\n'); fprintf(fid,'SALIDA V(a,b,c),I(a,b,c) Volt,Amp cartesiana\n'); fprintf(fid,'\n'); fprintf(fid,'Vma =%8.0f+i(%-8.0f)V\t',real(VVma),imag(VVma)); fprintf(fid,'Vmb =%8.0f+i(%-8.0f)V\t',real(VVmb),imag(VVmb)); fprintf(fid,'Vmc =%8.0f+i(%-8.0f)V\t',real(VVmc),imag(VVmc)); if TipoFalta==2 fprintf(fid,'\n'); fprintf(fid,'Vpa =%8.0f+i(%-8.0f)V\t',real(VVpa),imag(VVpa)); fprintf(fid,'Vpb =%8.0f+i(%-8.0f)V\t',real(VVpb),imag(VVpb)); fprintf(fid,'Vpc =%8.0f+i(%-8.0f)V\t',real(VVpc),imag(VVpc)); end fprintf(fid,'\n'); fprintf(fid,'Vsa =%8.0f+i(%-8.0f)V\t',real(VVsa),imag(VVsa)); fprintf(fid,'Vsb =%8.0f+i(%-8.0f)V\t',real(VVsb),imag(VVsb)); fprintf(fid,'Vsc =%8.0f+i(%-8.0f)V\t',real(VVsc),imag(VVsc)); fprintf(fid,'\n'); fprintf(fid,'Vra =%8.0f+i(%-8.0f)V\t',real(VVra),imag(VVra)); fprintf(fid,'Vrb =%8.0f+i(%-8.0f)V\t',real(VVrb),imag(VVrb)); fprintf(fid,'Vrc =%8.0f+i(%-8.0f)V\t',real(VVrc),imag(VVrc)); fprintf(fid,'\n'); fprintf(fid,'Vta =%8.0f+i(%-8.0f)V\t',real(VVta),imag(VVta)); fprintf(fid,'Vtb =%8.0f+i(%-8.0f)V\t',real(VVtb),imag(VVtb)); fprintf(fid,'Vtc =%8.0f+i(%-8.0f)V\t',real(VVtc),imag(VVtc)); fprintf(fid,'\n\n'); fprintf(fid,'Ima =%8.0f+i(%-8.0f)A\t',real(IIma),imag(IIma)); fprintf(fid,'Imb =%8.0f+i(%-8.0f)A\t',real(IImb),imag(IImb)); fprintf(fid,'Imc =%8.0f+i(%-8.0f)A\t',real(IImc),imag(IImc)); if TipoFalta==2 fprintf(fid,'\n'); fprintf(fid,'Ipa =%8.0f+i(%-8.0f)A\t',real(IIpa),imag(IIpa)); fprintf(fid,'Ipb =%8.0f+i(%-8.0f)A\t',real(IIpb),imag(IIpb)); fprintf(fid,'Ipc =%8.0f+i(%-8.0f)A\t',real(IIpc),imag(IIpc)); end fprintf(fid,'\n'); fprintf(fid,'Isa =%8.0f+i(%-8.0f)A\t',real(IIsa),imag(IIsa)); fprintf(fid,'Isb =%8.0f+i(%-8.0f)A\t',real(IIsb),imag(IIsb)); fprintf(fid,'Isc =%8.0f+i(%-8.0f)A\t',real(IIsc),imag(IIsc)); fprintf(fid,'\n'); fprintf(fid,'Ira =%8.0f+i(%-8.0f)A\t',real(IIra),imag(IIra)); fprintf(fid,'Irb =%8.0f+i(%-8.0f)A\t',real(IIrb),imag(IIrb)); fprintf(fid,'Irc =%8.0f+i(%-8.0f)A\t',real(IIrc),imag(IIrc)); fprintf(fid,'\n'); fprintf(fid,'Ita =%8.0f+i(%-8.0f)A\t',real(IIta),imag(IIta)); fprintf(fid,'Itb =%8.0f+i(%-8.0f)A\t',real(IItb),imag(IItb)); fprintf(fid,'Itc =%8.0f+i(%-8.0f)A\t',real(IItc),imag(IItc)); fprintf(fid,'\n'); fprintf(fid,'Isra =%8.0f+i(%-8.0f)A\t',real(IIsra),imag(IIsra)); fprintf(fid,'Isrb =%8.0f+i(%-8.0f)A\t',real(IIsrb),imag(IIsrb)); fprintf(fid,'Isrc =%8.0f+i(%-8.0f)A\t',real(IIsrc),imag(IIsrc)); fprintf(fid,'\n'); fprintf(fid,'Ista =%8.0f+i(%-8.0f)A\t',real(IIsta),imag(IIsta)); fprintf(fid,'Istb =%8.0f+i(%-8.0f)A\t',real(IIstb),imag(IIstb)); fprintf(fid,'Istc =%8.0f+i(%-8.0f)A\t',real(IIstc),imag(IIstc)); fprintf(fid,'\n'); fprintf(fid,'Irsa =%8.0f+i(%-8.0f)A\t',real(IIrsa),imag(IIrsa)); fprintf(fid,'Irsb =%8.0f+i(%-8.0f)A\t',real(IIrsb),imag(IIrsb)); fprintf(fid,'Irsc =%8.0f+i(%-8.0f)A\t',real(IIrsc),imag(IIrsc)); fprintf(fid,'\n'); fprintf(fid,'Itsa =%8.0f+i(%-8.0f)A\t',real(IItsa),imag(IItsa)); fprintf(fid,'Itsb =%8.0f+i(%-8.0f)A\t',real(IItsb),imag(IItsb)); fprintf(fid,'Itsc =%8.0f+i(%-8.0f)A\t',real(IItsc),imag(IItsc)); fprintf(fid,'\n'); % ****************************************************************** case{4}, % OVA polar: 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c) % --- --- % ****************************************************************** % ------------------------------------------------------------------ % ** C4 ****************** SALIDA 21, Ohm , polar ****************** % ------------------------------------------------------------------ % Sm(e1,e3),Sp_2,Sr_1,Rm(e1,e2),Rp_2,Rs_1 % ------------------------------------------------------------------ fprintf(fid,'\n'); fprintf(fid,'SALIDA 21, Ohm , polar \n'); fprintf(fid,'\n'); fprintf(fid,'Nota: Componentes de Z >=10000 Ohm se representan por Inf\n'); fprintf(fid,'Nota: Componentes de Z <=-10000 Ohm se representan por -Inf\n'); fprintf(fid,'\n'); if pfaltam==6 | pfaltam==7 [mod,arg]=mod_arg(Sm);fprintf(fid,'Sm=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_arg(Rm);fprintf(fid,'Rm=%7.2f(%-5.1f)(Ohm)\n',mod,arg); elseif pfaltam==8 | pfaltam==9 [mod,arg]=mod_arg(Sm);fprintf(fid,'Sm=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_arg(Tm);fprintf(fid,'Tm=%7.2f(%-5.1f)(Ohm)\n',mod,arg); end if TipoFalta==2 if pfaltap==6 | pfaltap==7 [mod,arg]=mod_arg(Sp);fprintf(fid,'Sp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_arg(Rp);fprintf(fid,'Rp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); elseif pfaltap==8 | pfaltap==9 [mod,arg]=mod_arg(Sp);fprintf(fid,'Sp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_arg(Tp);fprintf(fid,'Tp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); end end fprintf(fid,'\n'); [mod,arg]=mod_argz(Sr_ZAB*zbase);fprintf(fid,'Sr_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sr_ZBC*zbase);fprintf(fid,'Sr_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sr_ZCA*zbase);fprintf(fid,'Sr_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sr_ZAG*zbase);fprintf(fid,'Sr_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sr_ZBG*zbase);fprintf(fid,'Sr_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sr_ZCG*zbase);fprintf(fid,'Sr_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_argz(St_ZAB*zbase);fprintf(fid,'St_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(St_ZBC*zbase);fprintf(fid,'St_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(St_ZCA*zbase);fprintf(fid,'St_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(St_ZAG*zbase);fprintf(fid,'St_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(St_ZBG*zbase);fprintf(fid,'St_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(St_ZCG*zbase);fprintf(fid,'St_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_argz(Rs_ZAB*zbase);fprintf(fid,'Rs_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rs_ZBC*zbase);fprintf(fid,'Rs_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rs_ZCA*zbase);fprintf(fid,'Rs_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rs_ZAG*zbase);fprintf(fid,'Rs_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rs_ZBG*zbase);fprintf(fid,'Rs_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rs_ZCG*zbase);fprintf(fid,'Rs_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_argz(Ts_ZAB*zbase);fprintf(fid,'Ts_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Ts_ZBC*zbase);fprintf(fid,'Ts_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Ts_ZCA*zbase);fprintf(fid,'Ts_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Ts_ZAG*zbase);fprintf(fid,'Ts_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Ts_ZBG*zbase);fprintf(fid,'Ts_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Ts_ZCG*zbase);fprintf(fid,'Ts_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); fprintf(fid,'\n'); % ----------------------------------------------------------------------- % ** C4 ****************** SALIDA 67N, Volt, Amp, Polar ***************** % ----------------------------------------------------------------------- % Sm(e1,e3),Sp_2,Sr_1,Rm(e1,e2),Rp_2,Rs_1 % ----------------------------------------------------------------------- fprintf(fid,'\n'); fprintf(fid,'SALIDA 67N, Volt, Amp, polar\n'); fprintf(fid,'\n'); fprintf(fid,'Sr_67N:\tVspol=%-8.0f(V)\tIsrop=%-8.0f(A)\tAngDirSr=%-5.1f\n',VVspol,IIsrop,AngDirSr); fprintf(fid,'St_67N:\tVspol=%-8.0f(V)\tIstop=%-8.0f(A)\tAngDirSt=%-5.1f\n',VVspol,IIstop,AngDirSt); fprintf(fid,'Rs_67N:\tVrpol=%-8.0f(V)\tIrsop=%-8.0f(A)\tAngDirRs=%-5.1f\n',VVrpol,IIrsop,AngDirRs); fprintf(fid,'Ts_67N:\tVtpol=%-8.0f(V)\tItsop=%-8.0f(A)\tAngDirTs=%-5.1f\n',VVtpol,IItsop,AngDirTs); fprintf(fid,'\n'); % ------------------------------------------------------------------ % ** C4 ********* SALIDA V(0,1,2), I(0,1,2) Volt,Amp polar ******** % ------------------------------------------------------------------ % Vm,Vp_2,Vs,Vr % Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr % ------------------------------------------------------------------ fprintf(fid,'\n'); fprintf(fid,'SALIDA V(0,1,2), I(0,1,2) Volt-Amp polar\n'); fprintf(fid,'\n'); [mod,arg]=mod_arg(VVm0);fprintf(fid,'Vm0 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVm1);fprintf(fid,'Vm1 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVm2);fprintf(fid,'Vm2 =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); if TipoFalta==2 [mod,arg]=mod_arg(VVp0);fprintf(fid,'Vp0 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVp1);fprintf(fid,'Vp1 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVp2);fprintf(fid,'Vp2 =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); end [mod,arg]=mod_arg(VVs0);fprintf(fid,'Vs0 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVs1);fprintf(fid,'Vs1 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVs2);fprintf(fid,'Vs2 =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(VVr0);fprintf(fid,'Vr0 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVr1);fprintf(fid,'Vr1 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVr2);fprintf(fid,'Vr2 =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(VVt0);fprintf(fid,'Vt0 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVt1);fprintf(fid,'Vt1 =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVt2);fprintf(fid,'Vt2 =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); %------------------------------------------------------------------------------- fprintf(fid,'\n'); %------------------------------------------------------------------------------- [mod,arg]=mod_arg(IIm0);fprintf(fid,'Im0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIm1);fprintf(fid,'Im1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIm2);fprintf(fid,'Im2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); if TipoFalta==2 [mod,arg]=mod_arg(IIp0);fprintf(fid,'Ip0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIp1);fprintf(fid,'Ip1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIp2);fprintf(fid,'Ip2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); end [mod,arg]=mod_arg(IIs0);fprintf(fid,'Is0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIs1);fprintf(fid,'Is1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIs2);fprintf(fid,'Is2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIr0);fprintf(fid,'Ir0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIr1);fprintf(fid,'Ir1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIr2);fprintf(fid,'Ir2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIt0);fprintf(fid,'It0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIt1);fprintf(fid,'It1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIt2);fprintf(fid,'It2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIsr0);fprintf(fid,'Isr0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsr1);fprintf(fid,'Isr1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsr2);fprintf(fid,'Isr2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIst0);fprintf(fid,'Ist0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIst1);fprintf(fid,'Ist1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIst2);fprintf(fid,'Ist2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIrs0);fprintf(fid,'Irs0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrs1);fprintf(fid,'Irs1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrs2);fprintf(fid,'Irs2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIts0);fprintf(fid,'Its0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIts1);fprintf(fid,'Its1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIts2);fprintf(fid,'Its2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); % ------------------------------------------------------------------ % ** C4 ******** SALIDA V(a,b,c), I(a,b,c) Volt,Amp polar ********* % ------------------------------------------------------------------ % Vm,Vp_2,Vs,Vr % Im,Ip_2,Is,Ir,Ism(e1,e3),Isp_2,Isr_1,Irm(e1,e2),Irp_2,Irs_1,Izsr % ------------------------------------------------------------------ fprintf(fid,'\n'); fprintf(fid,'SALIDA V(a,b,c), I(a,b,c) Volt,Amp Polar\n'); fprintf(fid,'\n'); [mod,arg]=mod_arg(VVma);fprintf(fid,'Vma =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVmb);fprintf(fid,'Vmb =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVmc);fprintf(fid,'Vmc =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); if TipoFalta==2 [mod,arg]=mod_arg(VVpa);fprintf(fid,'Vpa =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVpb);fprintf(fid,'Vpb =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVpc);fprintf(fid,'Vpc =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); end [mod,arg]=mod_arg(VVsa);fprintf(fid,'Vsa =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVsb);fprintf(fid,'Vsb =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVsc);fprintf(fid,'Vsc =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(VVra);fprintf(fid,'Vra =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVrb);fprintf(fid,'Vrb =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVrc);fprintf(fid,'Vrc =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(VVta);fprintf(fid,'Vta =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVtb);fprintf(fid,'Vtb =%8.0f(%-5.1f)(V)\t',mod,arg); [mod,arg]=mod_arg(VVtc);fprintf(fid,'Vtc =%8.0f(%-5.1f)(V)',mod,arg); fprintf(fid,'\n'); %------------------------------------------------------------------------------- fprintf(fid,'\n'); %------------------------------------------------------------------------------- [mod,arg]=mod_arg(IIma);fprintf(fid,'Ima =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IImb);fprintf(fid,'Imb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IImc);fprintf(fid,'Imc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); if TipoFalta==2 [mod,arg]=mod_arg(IIpa);fprintf(fid,'Ipa =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIpb);fprintf(fid,'Ipb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIpc);fprintf(fid,'Ipc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); end [mod,arg]=mod_arg(IIsa);fprintf(fid,'Isa =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsb);fprintf(fid,'Isb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsc);fprintf(fid,'Isc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIra);fprintf(fid,'Ira =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrb);fprintf(fid,'Irb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrc);fprintf(fid,'Irc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIta);fprintf(fid,'Ita =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IItb);fprintf(fid,'Itb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IItc);fprintf(fid,'Itc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n\n'); [mod,arg]=mod_arg(IIsra);fprintf(fid,'Isra =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsrb);fprintf(fid,'Isrb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsrc);fprintf(fid,'Isrc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIsta);fprintf(fid,'Ista =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIstb);fprintf(fid,'Istb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIstc);fprintf(fid,'Istc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIrsa);fprintf(fid,'Irsa =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrsb);fprintf(fid,'Irsb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrsc);fprintf(fid,'Irsc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IItsa);fprintf(fid,'Itsa =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IItsb);fprintf(fid,'Itsb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IItsc);fprintf(fid,'Itsc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); end % end del switch(output) fprintf(fid,'\n'); variassalidas=input('Mas Salidas=1,Acabar=2 :'); fprintf(fid,'\n'); end % end del while variassalidas if salida==2 fclose(fid); end fprintf(1,'\n'); variasejecuciones=input('Mas Ejecuciones=1,Acabar=2 :'); end % end del while variasejecuciones fprintf(1,'Pulse Enter\n'); 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  >  PC  >  PC44  >  FC44.m