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,'\n'); fprintf(fid,'PC126M\n'); fprintf(fid,'Esquema Linea Doble SR\n'); fprintf(fid,'Unominal kV: 400,220,132,66\n'); fprintf(fid,'Eleccion: kE,L,Sccs,Sccr,Kr0,k0)\n'); fprintf(fid,'Ctes. de la Linea Doble, por Linea y por km, segun Unominal\n'); fprintf(fid,'Impedancia de Transferencia Si, No\n'); fprintf(fid,'Eleccion: (zsr1,Ksr0)\n'); fprintf(fid,'Compensacion 21, Si, No\n'); fprintf(fid,'Falta Intercircuito entre Lineas de igual Tension\n'); fprintf(fid,'m=p=0...1\n'); fprintf(fid,'Barra m=1, Barra p=2\n'); fprintf(fid,'\n'); % ------------------------------------------------------ %PC26 datos=input('datos: standard=1,input=2: '); if datos==1 Unominal=220; kE=1; L=220; imptransfer=2; Compens21=2; Sccs=5000; Sccr=5000; Kr0=1; % ------------------------------------------------------ else Unominal=input('Unominal kV: 400,220,132,66: '); kE=input('kE=E/Enominal: '); L=input('L(km): '); imptransfer=input('zsr Si=1, No=2: '); if imptransfer==1 zsr1=input('zsr1=zsr2(Complejo r+xi en Ohms): '); zsr2=zsr1; Ksr0=input('Ksr0: '); zsr0=Ksr0*zsr1; end Compens21=input('Compensacion 21:Si=1,No=2: '); Sccs=input('Sccs(MVA): '); Sccr=input('Sccr(MVA): '); Kr0=input('Kr0: '); end k0=input('k0=Ks0/Kr0: '); m=input('m=p pu : '); p=m; % ------------------------------------------------------- Ubase=Unominal; Vbase=Ubase/(sqrt(3)); MVAbase=100; Ibase=MVAbase/(3*Vbase); zbase=Vbase/Ibase; if Unominal==400 zl1u=(0.0346+0.3208i)/zbase; zl2u=zl1u; zl0u=(0.1828+1.1482i)/zbase; zlM0u=(0.1482+0.7383i)/zbase; elseif Unominal==220 zl1u=(0.0829+0.4135i)/zbase; zl2u=zl1u; zl0u=(0.2311+1.3085i)/zbase; zlM0u=(0.1482+0.8031i)/zbase; elseif Unominal==132 zl1u=(0.1342+0.4042i)/zbase; zl2u=zl1u; zl0u=(0.2824+1.3718i)/zbase; zlM0u=(0.1482+0.8947i)/zbase; elseif Unominal==66 zl1u=(0.2156+0.3952i)/zbase; zl2u=zl1u; zl0u=(0.3636+1.4278i)/zbase; zlM0u=(0.1482+0.9288i)/zbase; end % ------------------------------------------------------- % CALCULOS % ------------------------------------------------------- U=kE*Unominal; % kV E=U/sqrt(3)+0i; E=E/Vbase; % E pu k=(zl0u-zl1u)/(3*zl1u); Km=zlM0u/zl1u; zs1=(MVAbase/Sccs)*i; % pu zr1=(MVAbase/Sccr)*i; % pu zs2=zs1; % pu zr2=zr1; % pu zr0=Kr0*zr1; % pu zs0=k0*zr0; % pu zl0=zl0u*L; zl1=zl1u*L; zl2=zl2u*L; zlM0=zlM0u*L; % pu Sm=m*zl1; Rm=(1-m)*zl1; % pu Sp=p*zl1; Rp=(1-p)*zl1; % pu SSm=Sm*zbase; RRm=Rm*zbase; % Ohm SSp=Sp*zbase; RRp=Rp*zbase; % Ohm if imptransfer==1 zsr1=zsr1/zbase; zsr2=zsr1; zsr0=zsr0/zbase; end if imptransfer==1 [Minc]=subMincsr(1); [Minca]=subMincasr(1); [z0]=subz0sr(zs0,zsr0,zr0,zl0,zlM0,m,p); [z1]=subz1sr(zs1,zsr1,zr1,zl1,m,p); [z2]=subz2sr(zs2,zsr2,zr2,zl2,m,p); elseif imptransfer==2 [Minc]=subMinc(1); [Minca]=subMinca(1); [z0]=subz0(zs0,zr0,zl0,zlM0,m,p); [z1]=subz1(zs1,zr1,zl1,m,p); [z2]=subz2(zs2,zr2,zl2,m,p); end % MATRIZ DE ADMITANCIA DE BARRAS [Y0]=subY0(Minc,z0); [Y1]=subY1(Minc,z1); [Y2]=subY2(Minc,z2); % ZBUS(0,1,2) Z0=inv(Y0); Z1=inv(Y1); Z2=inv(Y2); % 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,falta,Rf,Rg]=... maMmp(Z0,Z1,Z2,E,zbase,fid); % pu % IBUS(0,1,2) I0=[0;0;-Im0;-Ip0]; I1=[0;0;-Im1;-Ip1]; I2=[0;0;-Im2;-Ip2]; % 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 (0,1,2) Vbase=Vbase*1000; Ibase=Ibase*1000; VVm0=Vm0*Vbase; VVm1=Vm1*Vbase; VVm2=Vm2*Vbase; VVp0=Vp0*Vbase; VVp1=Vp1*Vbase; VVp2=Vp2*Vbase; IIm0=Im0*Ibase; IIm1=Im1*Ibase; IIm2=Im2*Ibase; IIp0=Ip0*Ibase; IIp1=Ip1*Ibase; IIp2=Ip2*Ibase; 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; Is0=Ie0(1); Is1=Ie1(1); Is2=Ie2(1); Ir0=Ie0(2); Ir1=Ie1(2); Ir2=Ie2(2); Ism0=Ie0(3); Ism1=Ie1(3); Ism2=Ie2(3); Isp0=Ie0(4); Isp1=Ie1(4); Isp2=Ie2(4); Irm0=Ie0(5); Irm1=Ie1(5); Irm2=Ie2(5); Irp0=Ie0(6); Irp1=Ie1(6); Irp2=Ie2(6); if imptransfer==1 Izsr0=Ie0(7); Izsr1=Ie1(7); Izsr2=Ie2(7); end IIs0=Is0*Ibase; IIs1=Is1*Ibase; IIs2=Is2*Ibase; IIr0=Ir0*Ibase; IIr1=Ir1*Ibase; IIr2=Ir2*Ibase; IIsm0=Ism0*Ibase; IIsm1=Ism1*Ibase; IIsm2=Ism2*Ibase; IIsp0=Isp0*Ibase; IIsp1=Isp1*Ibase; IIsp2=Isp2*Ibase; IIrm0=Irm0*Ibase; IIrm1=Irm1*Ibase; IIrm2=Irm2*Ibase; IIrp0=Irp0*Ibase; IIrp1=Irp1*Ibase; IIrp2=Irp2*Ibase; if imptransfer==1 IIzsr0=Izsr0*Ibase; IIzsr1=Izsr1*Ibase; IIzsr2=Izsr2*Ibase; end % ASIGNACIONES V,I (a,b,c) [Vma,Vmb,Vmc]=abc(Vm0,Vm1,Vm2); [Vpa,Vpb,Vpc]=abc(Vp0,Vp1,Vp2); [Vsa,Vsb,Vsc]=abc(Vs0,Vs1,Vs2); [Vra,Vrb,Vrc]=abc(Vr0,Vr1,Vr2); [Ima,Imb,Imc]=abc(Im0,Im1,Im2); [Ipa,Ipb,Ipc]=abc(Ip0,Ip1,Ip2); [Isa,Isb,Isc]=abc(Is0,Is1,Is2); [Ira,Irb,Irc]=abc(Ir0,Ir1,Ir2); [Isma,Ismb,Ismc]=abc(Ism0,Ism1,Ism2); [Ispa,Ispb,Ispc]=abc(Isp0,Isp1,Isp2); [Irma,Irmb,Irmc]=abc(Irm0,Irm1,Irm2); [Irpa,Irpb,Irpc]=abc(Irp0,Irp1,Irp2); if imptransfer==1 [Izsra,Izsrb,Izsrc]=abc(Izsr0,Izsr1,Izsr2); end [VVma,VVmb,VVmc]=abc(VVm0,VVm1,VVm2); [VVpa,VVpb,VVpc]=abc(VVp0,VVp1,VVp2); [VVsa,VVsb,VVsc]=abc(VVs0,VVs1,VVs2); [VVra,VVrb,VVrc]=abc(VVr0,VVr1,VVr2); [IIma,IImb,IImc]=abc(IIm0,IIm1,IIm2); [IIpa,IIpb,IIpc]=abc(IIp0,IIp1,IIp2); [IIsa,IIsb,IIsc]=abc(IIs0,IIs1,IIs2); [IIra,IIrb,IIrc]=abc(IIr0,IIr1,IIr2); [IIsma,IIsmb,IIsmc]=abc(IIsm0,IIsm1,IIsm2); [IIspa,IIspb,IIspc]=abc(IIsp0,IIsp1,IIsp2); [IIrma,IIrmb,IIrmc]=abc(IIrm0,IIrm1,IIrm2); [IIrpa,IIrpb,IIrpc]=abc(IIrp0,IIrp1,IIrp2); if imptransfer==1 [IIzsra,IIzsrb,IIzsrc]=abc(IIzsr0,IIzsr1,IIzsr2); end % PROTECCION DE DISTANCIA (21) switch(Compens21) case{1}, % COMPENSADA [Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=... PDistComp(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k,Km,Isp0); [Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=... PDistComp(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k,Km,Ism0); [Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=... PDistComp(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k,Km,Irp0); [Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=... PDistComp(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k,Km,Irm0); case{2}, % NO COMPENSADA [Sm_ZAB,Sm_ZBC,Sm_ZCA,Sm_ZAG,Sm_ZBG,Sm_ZCG]=... PDist(Vsa,Vsb,Vsc,Isma,Ismb,Ismc,Ism0,k); [Sp_ZAB,Sp_ZBC,Sp_ZCA,Sp_ZAG,Sp_ZBG,Sp_ZCG]=... PDist(Vsa,Vsb,Vsc,Ispa,Ispb,Ispc,Isp0,k); [Rm_ZAB,Rm_ZBC,Rm_ZCA,Rm_ZAG,Rm_ZBG,Rm_ZCG]=... PDist(Vra,Vrb,Vrc,Irma,Irmb,Irmc,Irm0,k); [Rp_ZAB,Rp_ZBC,Rp_ZCA,Rp_ZAG,Rp_ZBG,Rp_ZCG]=... PDist(Vra,Vrb,Vrc,Irpa,Irpb,Irpc,Irp0,k); end % 67N [Vspol,Ismop,AngDirSm]=PDire(Vs0,Ism0); [Vspol,Ispop,AngDirSp]=PDire(Vs0,Isp0); [Vrpol,Irmop,AngDirRm]=PDire(Vr0,Irm0); [Vrpol,Irpop,AngDirRp]=PDire(Vr0,Irp0); VVspol=Vspol*Vbase; IIsmop=Ismop*Ibase; IIspop=Ispop*Ibase; VVrpol=Vrpol*Vbase; IIrmop=Irmop*Ibase; IIrpop=Irpop*Ibase; fprintf(fid,'\n'); % ----------------------------------------------------------------------- % SALIDAS % ----------------------------------------------------------------------- fprintf(fid,'********************************************************\n'); fprintf(fid,'Unominal =%-4.0f(kV)\n',Unominal); fprintf(fid,'kE=E/Enominal =%-4.0f\n',kE); fprintf(fid,'Impedancia de Transferencia,Si(1),No(2)=%-1.0f\n',imptransfer); fprintf(fid,'Compensacion 21,Si(1),No(2)=%-1.0f\n',Compens21); fprintf(fid,'L=%-5.1f(km)\n',L); fprintf(fid,'Sccs=%-5.0f(MVA)\n',Sccs); fprintf(fid,'Sccr=%-5.0f(MVA)\n',Sccr); fprintf(fid,'Kr0=%-6.2f\n',Kr0); fprintf(fid,'k0=%-6.2f\n',k0); fprintf(fid,'Ibase=%-8.2f(A)\n',Ibase*1000); if imptransfer==1 fprintf(fid,'zsr1=zsr2=%-7.2f+i(%-7.2f)Ohm\n',real(zsr1),imag(zsr1)); fprintf(fid,'zsr0=%-7.2f+i(%-7.2f)Ohm\n',real(zsr0),imag(zsr0)); end if falta==1 cfalta='Aa';elseif falta==2 cfalta='Bc';elseif falta==3 cfalta='Cb';elseif falta==4 cfalta='Aag'; elseif falta==5 cfalta='Bcg';elseif falta==6 cfalta='Cbg';elseif falta==7 cfalta='Aab'; elseif falta==8 cfalta='Bac';elseif falta==9 cfalta='Cbc';elseif falta==10 cfalta='Aabg'; elseif falta==11 cfalta='Bacg';elseif falta==12 cfalta='Cbcg';elseif falta==13 cfalta='Aabc'; elseif falta==14 cfalta='Aabcg';elseif falta==15 cfalta='ABCabc';elseif falta==16 cfalta='ABCabcg'; end fprintf(fid,'\n'); fprintf(fid,'Falta:%s\n',cfalta); fprintf(fid,'Rf:%-10.2f(Ohm)\n',Rf*zbase); if (falta==4)|(falta==5)|(falta==6)|(falta==10)|(falta==11)|(falta==12)|(falta==14)|(falta==16) fprintf(fid,'Rg:%-10.2f(Ohm)\n',Rg*zbase); end fprintf(fid,'\n'); fprintf(fid,'********************************************************\n'); variassalidas=1; while variassalidas==1 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 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.1 SALIDA 21 pu cartesiana % ----------------------------------------------------------------------- % Sm,Sp,Rm,Rp % ----------------------------------------------------------------------- fprintf(fid,'\n'); fprintf(fid,'Sm =%-8.4f+i(%-8.4f)\n',real(Sm),imag(Sm)); fprintf(fid,'Rm =%-8.4f+i(%-8.4f)\n',real(Rm),imag(Rm)); fprintf(fid,'Sp = %-8.4f+i(%-8.4f)\n',real(Sp),imag(Sp)); fprintf(fid,'Rp = %8.4f+i(%-8.4f)\n',real(Rp),imag(Rp)); fprintf(fid,'\n'); fprintf(fid,'Sm_ZAB:%8.4f+i(%-8.4f)\n',real(Sm_ZAB),imag(Sm_ZAB)); fprintf(fid,'Sm_ZBC:%8.4f+i(%-8.4f)\n',real(Sm_ZBC),imag(Sm_ZBC)); fprintf(fid,'Sm_ZCA:%8.4f+i(%-8.4f)\n',real(Sm_ZCA),imag(Sm_ZCA)); fprintf(fid,'Sm_ZAG:%8.4f+i(%-8.4f)\n',real(Sm_ZAG),imag(Sm_ZAG)); fprintf(fid,'Sm_ZBG:%8.4f+i(%-8.4f)\n',real(Sm_ZBG),imag(Sm_ZBG)); fprintf(fid,'Sm_ZCG:%8.4f+i(%-8.4f)\n',real(Sm_ZCG),imag(Sm_ZCG)); fprintf(fid,'\n'); fprintf(fid,'Rm_ZAB:%8.4f+i(%-8.4f)\n',real(Rm_ZAB),imag(Rm_ZAB)); fprintf(fid,'Rm_ZBC:%8.4f+i(%-8.4f)\n',real(Rm_ZBC),imag(Rm_ZBC)); fprintf(fid,'Rm_ZCA:%8.4f+i(%-8.4f)\n',real(Rm_ZCA),imag(Rm_ZCA)); fprintf(fid,'Rm_ZAG:%8.4f+i(%-8.4f)\n',real(Rm_ZAG),imag(Rm_ZAG)); fprintf(fid,'Rm_ZBG:%8.4f+i(%-8.4f)\n',real(Rm_ZBG),imag(Rm_ZBG)); fprintf(fid,'Rm_ZCG:%8.4f+i(%-8.4f)\n',real(Rm_ZCG),imag(Rm_ZCG)); fprintf(fid,'\n'); fprintf(fid,'Sp_ZAB:%8.4f+i(%-8.4f)\n',real(Sp_ZAB),imag(Sp_ZAB)); fprintf(fid,'Sp_ZBC:%8.4f+i(%-8.4f)\n',real(Sp_ZBC),imag(Sp_ZBC)); fprintf(fid,'Sp_ZCA:%8.4f+i(%-8.4f)\n',real(Sp_ZCA),imag(Sp_ZCA)); fprintf(fid,'Sp_ZAG:%8.4f+i(%-8.4f)\n',real(Sp_ZAG),imag(Sp_ZAG)); fprintf(fid,'Sp_ZBG:%8.4f+i(%-8.4f)\n',real(Sp_ZBG),imag(Sp_ZBG)); fprintf(fid,'Sp_ZCG:%8.4f+i(%-8.4f)\n',real(Sp_ZCG),imag(Sp_ZCG)); fprintf(fid,'\n'); fprintf(fid,'Rp_ZAB:%8.4f+i(%-8.4f)\n',real(Rp_ZAB),imag(Rp_ZAB)); fprintf(fid,'Rp_ZBC:%8.4f+i(%-8.4f)\n',real(Rp_ZBC),imag(Rp_ZBC)); fprintf(fid,'Rp_ZCA:%8.4f+i(%-8.4f)\n',real(Rp_ZCA),imag(Rp_ZCA)); fprintf(fid,'Rp_ZAG:%8.4f+i(%-8.4f)\n',real(Rp_ZAG),imag(Rp_ZAG)); fprintf(fid,'Rp_ZBG:%8.4f+i(%-8.4f)\n',real(Rp_ZBG),imag(Rp_ZBG)); fprintf(fid,'Rp_ZCG:%8.4f+i(%-8.4f)\n',real(Rp_ZCG),imag(Rp_ZCG)); fprintf(fid,'\n'); % ----------------------------------------------------------------------- % C1.2 SALIDA 67N pu cartesiana % ----------------------------------------------------------------------- % Sm,Sp,Rm,Rp % ----------------------------------------------------------------------- % ----------------------------------------------------------------------- % C1.3 SALIDA V(0,1,2),I(0,1,2) pu cartesiana % ----------------------------------------------------------------------- % Vm,Vp,Vs,Vr % Im,Ip,Is,Ir,Ism,Isp,Irm,Irp,(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)',real(Vm0),imag(Vm0)); fprintf(fid,'\tVm1:%8.4f+i(%-8.4f)',real(Vm1),imag(Vm1)); fprintf(fid,'\tVm2:%8.4f+i(%-8.4f)\n',real(Vm2),imag(Vm2)); fprintf(fid,'Vp0:%8.4f+i(%-8.4f)',real(Vp0),imag(Vp0)); fprintf(fid,'\tVp1:%8.4f+i(%-8.4f)',real(Vp1),imag(Vp1)); fprintf(fid,'\tVp2:%8.4f+i(%-8.4f)\n',real(Vp2),imag(Vp2)); fprintf(fid,'Vs0:%8.4f+i(%-8.4f)',real(Vs0),imag(Vs0)); fprintf(fid,'\tVs1:%8.4f+i(%-8.4f)',real(Vs1),imag(Vs1)); fprintf(fid,'\tVs2:%8.4f+i(%-8.4f)\n',real(Vs2),imag(Vs2)); fprintf(fid,'Vr0:%8.4f+i(%-8.4f)',real(Vr0),imag(Vr0)); fprintf(fid,'\tVr1:%8.4f+i(%-8.4f)',real(Vr1),imag(Vr1)); fprintf(fid,'\tVr2:%8.4f+i(%-8.4f)',real(Vr2),imag(Vr2)); fprintf(fid,'\n\n'); fprintf(fid,'Im0:%8.4f+i(%-8.4f)',real(Im0),imag(Im0)); fprintf(fid,'\tIm1:%8.4f+i(%-8.4f)',real(Im1),imag(Im1)); fprintf(fid,'\tIm2:%8.4f+i(%-8.4f)\n',real(Im2),imag(Im2)); fprintf(fid,'Ip0:%8.4f+i(%-8.4f)',real(Ip0),imag(Ip0)); fprintf(fid,'\tIp1:%8.4f+i(%-8.4f)',real(Ip1),imag(Ip1)); fprintf(fid,'\tIp2:%8.4f+i(%-8.4f)\n',real(Ip2),imag(Ip2)); fprintf(fid,'Is0:%8.4f+i(%-8.4f)',real(Is0),imag(Is0)); fprintf(fid,'\tIs1:%8.4f+i(%-8.4f)',real(Is1),imag(Is1)); fprintf(fid,'\tIs2:%8.4f+i(%-8.4f)\n',real(Is2),imag(Is2)); fprintf(fid,'Ir0:%8.4f+i(%-8.4f)',real(Ir0),imag(Ir0)); fprintf(fid,'\tIr1:%8.4f+i(%-8.4f)',real(Ir1),imag(Ir1)); fprintf(fid,'\tIr2:%8.4f+i(%-8.4f)\n',real(Ir2),imag(Ir2)); fprintf(fid,'Ism0:%8.4f+i(%-8.4f)',real(Ism0),imag(Ism0)); fprintf(fid,'\tIsm1:%8.4f+i(%-8.4f)',real(Ism1),imag(Ism1)); fprintf(fid,'\tIsm2:%8.4f+i(%-8.4f)\n',real(Ism2),imag(Ism2)); fprintf(fid,'Isp0:%8.4f+i(%-8.4f)',real(Isp0),imag(Isp0)); fprintf(fid,'\tIsp1:%8.4f+i(%-8.4f)',real(Isp1),imag(Isp1)); fprintf(fid,'\tIsp2:%8.4f+i(%-8.4f)\n',real(Isp2),imag(Isp2)); fprintf(fid,'Irm0:%8.4f+i(%-8.4f)',real(Irm0),imag(Irm0)); fprintf(fid,'\tIrm1:%8.4f+i(%-8.4f)',real(Irm1),imag(Irm1)); fprintf(fid,'\tIrm2:%8.4f+i(%-8.4f)\n',real(Irm2),imag(Irm2)); fprintf(fid,'Irp0:%8.4f+i(%-8.4f)',real(Irp0),imag(Irp0)); fprintf(fid,'\tIrp1:%8.4f+i(%-8.4f)',real(Irp1),imag(Irp1)); fprintf(fid,'\tIrp2:%8.4f+i(%-8.4f)',real(Irp2),imag(Irp2)); fprintf(fid,'\n'); if imptransfer==1 fprintf(fid,'Izsr0:%8.4f+i(%-8.4f)A',real(Izsr0),imag(Izsr0)); fprintf(fid,'\tIzsr1:%8.4f+i(%-8.4f)A',real(Izsr1),imag(Izsr1)); fprintf(fid,'\tIzsr2:%8.4f+i(%-8.4f)A',real(Izsr2),imag(Izsr2)); fprintf(fid,'\n'); end % ----------------------------------------------------------------------- % C1.4 SALIDA V(a,b,c),I(a,b,c) pu cartesiana % ----------------------------------------------------------------------- % Vm,Vp,Vs,Vr % Im,Ip,Is,Ir,Ism,Isp,Irm,Irp %,Izsr % ----------------------------------------------------------------------- fprintf(fid,'\n'); fprintf(fid,'Vma:%8.4f+i(%-8.4f)',real(Vma),imag(Vma)); fprintf(fid,'\tVmb:%8.4f+i(%-8.4f)',real(Vmb),imag(Vmb)); fprintf(fid,'\tVmc:%8.4f+i(%-8.4f)\n',real(Vmc),imag(Vmc)); fprintf(fid,'Vpa:%8.4f+i(%-8.4f)',real(Vpa),imag(Vpa)); fprintf(fid,'\tVpb:%8.4f+i(%-8.4f)',real(Vpb),imag(Vpb)); fprintf(fid,'\tVpc:%8.4f+i(%-8.4f)\n',real(Vpc),imag(Vpc)); fprintf(fid,'Vsa:%8.4f+i(%-8.4f)',real(Vsa),imag(Vsa)); fprintf(fid,'\tVsb:%8.4f+i(%-8.4f)',real(Vsb),imag(Vsb)); fprintf(fid,'\tVsc:%8.4f+i(%-8.4f)\n',real(Vsc),imag(Vsc)); fprintf(fid,'Vra:%8.4f+i(%-8.4f)',real(Vra),imag(Vra)); fprintf(fid,'\tVrb:%8.4f+i(%-8.4f)',real(Vrb),imag(Vrb)); fprintf(fid,'\tVrc:%8.4f+i(%-8.4f)\n',real(Vrc),imag(Vrc)); fprintf(fid,'\n'); fprintf(fid,'Ima:%8.4f+i(%-8.4f)',real(Ima),imag(Ima)); fprintf(fid,'\tImb:%8.4f+i(%-8.4f)',real(Imb),imag(Imb)); fprintf(fid,'\tImc:%8.4f+i(%-8.4f)\n',real(Imc),imag(Imc)); fprintf(fid,'Ipa:%8.4f+i(%-8.4f)',real(Ipa),imag(Ipa)); fprintf(fid,'\tIpb:%8.4f+i(%-8.4f)',real(Ipb),imag(Ipb)); fprintf(fid,'\tIpc:%8.4f+i(%-8.4f)\n',real(Ipc),imag(Ipc)); fprintf(fid,'Isa:%8.4f+i(%-8.4f)',real(Isa),imag(Isa)); fprintf(fid,'\tIsb:%8.4f+i(%-8.4f)',real(Isb),imag(Isb)); fprintf(fid,'\tIsc:%8.4f+i(%-8.4f)\n',real(Isc),imag(Isc)); fprintf(fid,'Ira:%8.4f+i(%-8.4f)',real(Ira),imag(Ira)); fprintf(fid,'\tIrb:%8.4f+i(%-8.4f)',real(Irb),imag(Irb)); fprintf(fid,'\tIrc:%8.4f+i(%-8.4f)\n',real(Irc),imag(Irc)); fprintf(fid,'Isma:%8.4f+i(%-8.4f)',real(Isma),imag(Isma)); fprintf(fid,'\tIsmb:%8.4f+i(%-8.4f)',real(Ismb),imag(Ismb)); fprintf(fid,'\tIsmc:%8.4f+i(%-8.4f)\n',real(Ismc),imag(Ismc)); fprintf(fid,'Ispa:%8.4f+i(%-8.4f)',real(Ispa),imag(Ispa)); fprintf(fid,'\tIspb:%8.4f+i(%-8.4f)',real(Ispb),imag(Ispb)); fprintf(fid,'\tIspc:%8.4f+i(%-8.4f)\n',real(Ispc),imag(Ispc)); fprintf(fid,'Irma:%8.4f+i(%-8.4f)',real(Irma),imag(Irma)); fprintf(fid,'\tIrmb:%8.4f+i(%-8.4f)',real(Irmb),imag(Irmb)); fprintf(fid,'\tIrmc:%8.4f+i(%-8.4f)\n',real(Irmc),imag(Irmc)); fprintf(fid,'Irpa:%8.4f+i(%-8.4f)',real(Irpa),imag(Irpa)); fprintf(fid,'\tIrpb:%8.4f+i(%-8.4f)',real(Irpb),imag(Irpb)); fprintf(fid,'\tIrpc:%8.4f+i(%-8.4f)\n',real(Irpc),imag(Irpc)); fprintf(fid,'\n'); if imptransfer==1 fprintf(fid,'Izsra:%8.4f+i(%-8.4f)A',real(Izsra),imag(Izsra)); fprintf(fid,'\tIzsrb:%8.4f+i(%-8.4f)A',real(Izsrb),imag(Izsrb)); fprintf(fid,'\tIzsrc:%8.4f+i(%-8.4f)A',real(Izsrc),imag(Izsrc)); fprintf(fid,'\n'); end % ********************************************************************** case{2}, % pu(polar): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c) % % ********************************************************************** % ********************************************************************** 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,Sp,Rm,Rp % ----------------------------------------------------------------------- fprintf(fid,'\n'); fprintf(fid,'Sm=Sp=%-5.2f+i(%-5.2f)\n',real(SSm),imag(SSm)); fprintf(fid,'Rm=Rp=%-5.2f+i(%-5.2f)\n',real(RRm),imag(RRm)); fprintf(fid,'\n'); [Sm_ZABre,Sm_ZABim]=zmed(Sm_ZAB,zbase); fprintf(fid,'Sm_ZAB:%7.2f+i(%-7.2f)Ohm\n',Sm_ZABre,Sm_ZABim); [Sm_ZBCre,Sm_ZBCim]=zmed(Sm_ZBC,zbase); fprintf(fid,'Sm_ZBC:%7.2f+i(%-7.2f)Ohm\n',Sm_ZBCre,Sm_ZBCim); [Sm_ZCAre,Sm_ZCAim]=zmed(Sm_ZCA,zbase); fprintf(fid,'Sm_ZCA:%7.2f+i(%-7.2f)Ohm\n',Sm_ZCAre,Sm_ZCAim); [Sm_ZAGre,Sm_ZAGim]=zmed(Sm_ZAG,zbase); fprintf(fid,'Sm_ZAG:%7.2f+i(%-7.2f)Ohm\n',Sm_ZAGre,Sm_ZAGim); [Sm_ZBGre,Sm_ZBGim]=zmed(Sm_ZBG,zbase); fprintf(fid,'Sm_ZBG:%7.2f+i(%-7.2f)Ohm\n',Sm_ZBGre,Sm_ZBGim); [Sm_ZCGre,Sm_ZCGim]=zmed(Sm_ZCG,zbase); fprintf(fid,'Sm_ZCG:%7.2f+i(%-7.2f)Ohm\n',Sm_ZCGre,Sm_ZCGim); fprintf(fid,'\n'); [Sp_ZABre,Sp_ZABim]=zmed(Sp_ZAB,zbase);fprintf(fid,'Sp_ZAB:%7.2f+i(%-7.2f)Ohm\n',Sp_ZABre,Sp_ZABim); [Sp_ZBCre,Sp_ZBCim]=zmed(Sp_ZBC,zbase);fprintf(fid,'Sp_ZBC:%7.2f+i(%-7.2f)Ohm\n',Sp_ZBCre,Sp_ZBCim); [Sp_ZCAre,Sp_ZCAim]=zmed(Sp_ZCA,zbase);fprintf(fid,'Sp_ZCA:%7.2f+i(%-7.2f)Ohm\n',Sp_ZCAre,Sp_ZCAim); [Sp_ZAGre,Sp_ZAGim]=zmed(Sp_ZAG,zbase);fprintf(fid,'Sp_ZAG:%7.2f+i(%-7.2f)Ohm\n',Sp_ZAGre,Sp_ZAGim); [Sp_ZBGre,Sp_ZBGim]=zmed(Sp_ZBG,zbase);fprintf(fid,'Sp_ZBG:%7.2f+i(%-7.2f)Ohm\n',Sp_ZBGre,Sp_ZBGim); [Sp_ZCGre,Sp_ZCGim]=zmed(Sp_ZCG,zbase);fprintf(fid,'Sp_ZCG:%7.2f+i(%-7.2f)Ohm\n',Sp_ZCGre,Sp_ZCGim); fprintf(fid,'\n'); [Rm_ZABre,Rm_ZABim]=zmed(Rm_ZAB,zbase);fprintf(fid,'Rm_ZAB:%7.2f+i(%-7.2f)Ohm\n',Rm_ZABre,Rm_ZABim); [Rm_ZBCre,Rm_ZBCim]=zmed(Rm_ZBC,zbase);fprintf(fid,'Rm_ZBC:%7.2f+i(%-7.2f)Ohm\n',Rm_ZBCre,Rm_ZBCim); [Rm_ZCAre,Rm_ZCAim]=zmed(Rm_ZCA,zbase);fprintf(fid,'Rm_ZCA:%7.2f+i(%-7.2f)Ohm\n',Rm_ZCAre,Rm_ZCAim); [Rm_ZAGre,Rm_ZAGim]=zmed(Rm_ZAG,zbase);fprintf(fid,'Rm_ZAG:%7.2f+i(%-7.2f)Ohm\n',Rm_ZAGre,Rm_ZAGim); [Rm_ZBGre,Rm_ZBGim]=zmed(Rm_ZBG,zbase);fprintf(fid,'Rm_ZBG:%7.2f+i(%-7.2f)Ohm\n',Rm_ZBGre,Rm_ZBGim); [Rm_ZCGre,Rm_ZCGim]=zmed(Rm_ZCG,zbase);fprintf(fid,'Rm_ZCG:%7.2f+i(%-7.2f)Ohm\n',Rm_ZCGre,Rm_ZCGim); fprintf(fid,'\n'); [Rp_ZABre,Rp_ZABim]=zmed(Rp_ZAB,zbase);fprintf(fid,'Rp_ZAB:%7.2f+i(%-7.2f)Ohm\n',Rp_ZABre,Rp_ZABim); [Rp_ZBCre,Rp_ZBCim]=zmed(Rp_ZBC,zbase);fprintf(fid,'Rp_ZBC:%7.2f+i(%-7.2f)Ohm\n',Rp_ZBCre,Rp_ZBCim); [Rp_ZCAre,Rp_ZCAim]=zmed(Rp_ZCA,zbase);fprintf(fid,'Rp_ZCA:%7.2f+i(%-7.2f)Ohm\n',Rp_ZCAre,Rp_ZCAim); [Rp_ZAGre,Rp_ZAGim]=zmed(Rp_ZAG,zbase);fprintf(fid,'Rp_ZAG:%7.2f+i(%-7.2f)Ohm\n',Rp_ZAGre,Rp_ZAGim); [Rp_ZBGre,Rp_ZBGim]=zmed(Rp_ZBG,zbase);fprintf(fid,'Rp_ZBG:%7.2f+i(%-7.2f)Ohm\n',Rp_ZBGre,Rp_ZBGim); [Rp_ZCGre,Rp_ZCGim]=zmed(Rp_ZCG,zbase);fprintf(fid,'Rp_ZCG:%7.2f+i(%-7.2f)Ohm\n',Rp_ZCGre,Rp_ZCGim); fprintf(fid,'\n'); % ------------------------------------------------------------------ % C3.3 SALIDA V(0,1,2),I(0,1,2) Volt-Amp cartesiana % ------------------------------------------------------------------ % Vm,Vp,Vs,Vr % Im,Ip,Is,Ir,Ism,Isp,Irm,Irp %,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',real(VVm0),imag(VVm0)); fprintf(fid,'\tVm1 :%8.0f+i(%-8.0f)V',real(VVm1),imag(VVm1)); fprintf(fid,'\tVm2 :%8.0f+i(%-8.0f)V',real(VVm2),imag(VVm2)); fprintf(fid,'\n'); fprintf(fid,'Vp0 :%8.0f+i(%-8.0f)V',real(VVp0),imag(VVp0)); fprintf(fid,'\tVp1 :%8.0f+i(%-8.0f)V',real(VVp1),imag(VVp1)); fprintf(fid,'\tVp2 :%8.0f+i(%-8.0f)V',real(VVp2),imag(VVp2)); fprintf(fid,'\n'); fprintf(fid,'Vs0 :%8.0f+i(%-8.0f)V',real(VVs0),imag(VVs0)); fprintf(fid,'\tVs1 :%8.0f+i(%-8.0f)V',real(VVs1),imag(VVs1)); fprintf(fid,'\tVs2 :%8.0f+i(%-8.0f)V',real(VVs2),imag(VVs2)); fprintf(fid,'\n'); fprintf(fid,'Vr0 :%8.0f+i(%-8.0f)V',real(VVr0),imag(VVr0)); fprintf(fid,'\tVr1 :%8.0f+i(%-8.0f)V',real(VVr1),imag(VVr1)); fprintf(fid,'\tVr2 :%8.0f+i(%-8.0f)V',real(VVr2),imag(VVr2)); fprintf(fid,'\n\n'); % --------------------------------------------------------------- fprintf(fid,'Im0 :%8.0f+i(%-8.0f)A',real(IIm0),imag(IIm0)); fprintf(fid,'\tIm1 :%8.0f+i(%-8.0f)A',real(IIm1),imag(IIm1)); fprintf(fid,'\tIm2 :%8.0f+i(%-8.0f)A',real(IIm2),imag(IIm2)); fprintf(fid,'\n'); fprintf(fid,'Ip0 :%8.0f+i(%-8.0f)A',real(IIp0),imag(IIp0)); fprintf(fid,'\tIp1 :%8.0f+i(%-8.0f)A',real(IIp1),imag(IIp1)); fprintf(fid,'\tIp2 :%8.0f+i(%-8.0f)A',real(IIp2),imag(IIp2)); fprintf(fid,'\n'); fprintf(fid,'Is0 :%8.0f+i(%-8.0f)A',real(IIs0),imag(IIs0)); fprintf(fid,'\tIs1 :%8.0f+i(%-8.0f)A',real(IIs1),imag(IIs1)); fprintf(fid,'\tIs2 :%8.0f+i(%-8.0f)A',real(IIs2),imag(IIs2)); fprintf(fid,'\n'); fprintf(fid,'Ir0 :%8.0f+i(%-8.0f)A',real(IIr0),imag(IIr0)); fprintf(fid,'\tIr1 :%8.0f+i(%-8.0f)A',real(IIr1),imag(IIr1)); fprintf(fid,'\tIr2 :%8.0f+i(%-8.0f)A',real(IIr2),imag(IIr2)); fprintf(fid,'\n'); fprintf(fid,'Ism0 :%8.0f+i(%-8.0f)A',real(IIsm0),imag(IIsm0)); fprintf(fid,'\tIsm1 :%8.0f+i(%-8.0f)A',real(IIsm1),imag(IIsm1)); fprintf(fid,'\tIsm2 :%8.0f+i(%-8.0f)A',real(IIsm2),imag(IIsm2)); fprintf(fid,'\n'); fprintf(fid,'Isp0 :%8.0f+i(%-8.0f)A',real(IIsp0),imag(IIsp0)); fprintf(fid,'\tIsp1 :%8.0f+i(%-8.0f)A',real(IIsp1),imag(IIsp1)); fprintf(fid,'\tIsp2 :%8.0f+i(%-8.0f)A',real(IIsp2),imag(IIsp2)); fprintf(fid,'\n'); fprintf(fid,'Irm0 :%8.0f+i(%-8.0f)A',real(IIrm0),imag(IIrm0)); fprintf(fid,'\tIrm1 :%8.0f+i(%-8.0f)A',real(IIrm1),imag(IIrm1)); fprintf(fid,'\tIrm2 :%8.0f+i(%-8.0f)A',real(IIrm2),imag(IIrm2)); fprintf(fid,'\n'); fprintf(fid,'Irp0 :%8.0f+i(%-8.0f)A',real(IIrp0),imag(IIrp0)); fprintf(fid,'\tIrp1 :%8.0f+i(%-8.0f)A',real(IIrp1),imag(IIrp1)); fprintf(fid,'\tIrp2 :%8.0f+i(%-8.0f)A',real(IIrp2),imag(IIrp2)); fprintf(fid,'\n'); if imptransfer==1 fprintf(fid,'Izsr0:%8.0f+i(%-8.0f)A',real(IIzsr0),imag(IIzsr0)); fprintf(fid,'\tIzsr1:%8.0f+i(%-8.0f)A',real(IIzsr1),imag(IIzsr1)); fprintf(fid,'\tIzsr2:%8.0f+i(%-8.0f)A',real(IIzsr2),imag(IIzsr2)); fprintf(fid,'\n'); end % ------------------------------------------------------------------ % C3.4 SALIDA V(a,b,c),I(a,b,c) Volt,Amp cartesiana % ------------------------------------------------------------------ % Vm,Vp,Vs,Vr % Im,Ip,Is,Ir,Ism,Isp,Irm,Irp %,Izsr % ------------------------------------------------------------------ fprintf(fid,'\n'); fprintf(fid,'SALIDA V(a,b,c),I(a,b,c) Volt,Amp\n'); fprintf(fid,'\n'); fprintf(fid,'Vma :%8.0f+i(%-8.0f)V',real(VVma),imag(VVma)); fprintf(fid,'\tVmb :%8.0f+i(%-8.0f)V',real(VVmb),imag(VVmb)); fprintf(fid,'\tVmc :%8.0f+i(%-8.0f)V',real(VVmc),imag(VVmc)); fprintf(fid,'\n'); fprintf(fid,'Vpa :%8.0f+i(%-8.0f)V',real(VVpa),imag(VVpa)); fprintf(fid,'\tVpb :%8.0f+i(%-8.0f)V',real(VVpb),imag(VVpb)); fprintf(fid,'\tVpc :%8.0f+i(%-8.0f)V',real(VVpc),imag(VVpc)); fprintf(fid,'\n'); fprintf(fid,'Vsa :%8.0f+i(%-8.0f)V',real(VVsa),imag(VVsa)); fprintf(fid,'\tVsb :%8.0f+i(%-8.0f)V',real(VVsb),imag(VVsb)); fprintf(fid,'\tVsc :%8.0f+i(%-8.0f)V',real(VVsc),imag(VVsc)); fprintf(fid,'\n'); fprintf(fid,'Vra :%8.0f+i(%-8.0f)V',real(VVra),imag(VVra)); fprintf(fid,'\tVrb :%8.0f+i(%-8.0f)V',real(VVrb),imag(VVrb)); fprintf(fid,'\tVrc :%8.0f+i(%-8.0f)V',real(VVrc),imag(VVrc)); fprintf(fid,'\n\n'); fprintf(fid,'Ima :%8.0f+i(%-8.0f)A',real(IIma),imag(IIma)); fprintf(fid,'\tImb :%8.0f+i(%-8.0f)A',real(IImb),imag(IImb)); fprintf(fid,'\tImc :%8.0f+i(%-8.0f)A',real(IImc),imag(IImc)); fprintf(fid,'\n'); fprintf(fid,'Ipa :%8.0f+i(%-8.0f)A',real(IIpa),imag(IIpa)); fprintf(fid,'\tIpb :%8.0f+i(%-8.0f)A',real(IIpb),imag(IIpb)); fprintf(fid,'\tIpc :%8.0f+i(%-8.0f)A',real(IIpc),imag(IIpc)); fprintf(fid,'\n'); fprintf(fid,'Isa :%8.0f+i(%-8.0f)A',real(IIsa),imag(IIsa)); fprintf(fid,'\tIsb :%8.0f+i(%-8.0f)A',real(IIsb),imag(IIsb)); fprintf(fid,'\tIsc :%8.0f+i(%-8.0f)A',real(IIsc),imag(IIsc)); fprintf(fid,'\n'); fprintf(fid,'Ira :%8.0f+i(%-8.0f)A',real(IIra),imag(IIra)); fprintf(fid,'\tIrb :%8.0f+i(%-8.0f)A',real(IIrb),imag(IIrb)); fprintf(fid,'\tIrc :%8.0f+i(%-8.0f)A',real(IIrc),imag(IIrc)); fprintf(fid,'\n'); fprintf(fid,'Isma :%8.0f+i(%-8.0f)A',real(IIsma),imag(IIsma)); fprintf(fid,'\tIsmb :%8.0f+i(%-8.0f)A',real(IIsmb),imag(IIsmb)); fprintf(fid,'\tIsmc :%8.0f+i(%-8.0f)A',real(IIsmc),imag(IIsmc)); fprintf(fid,'\n'); fprintf(fid,'Ispa :%8.0f+i(%-8.0f)A',real(IIspa),imag(IIspa)); fprintf(fid,'\tIspb :%8.0f+i(%-8.0f)A',real(IIspb),imag(IIspb)); fprintf(fid,'\tIspc :%8.0f+i(%-8.0f)A',real(IIspc),imag(IIspc)); fprintf(fid,'\n'); fprintf(fid,'Irma :%8.0f+i(%-8.0f)A',real(IIrma),imag(IIrma)); fprintf(fid,'\tIrmb :%8.0f+i(%-8.0f)A',real(IIrmb),imag(IIrmb)); fprintf(fid,'\tIrmc :%8.0f+i(%-8.0f)A',real(IIrmc),imag(IIrmc)); fprintf(fid,'\n'); fprintf(fid,'Irpa :%8.0f+i(%-8.0f)A',real(IIrpa),imag(IIrpa)); fprintf(fid,'\tIrpb :%8.0f+i(%-8.0f)A',real(IIrpb),imag(IIrpb)); fprintf(fid,'\tIrpc :%8.0f+i(%-8.0f)A',real(IIrpc),imag(IIrpc)); fprintf(fid,'\n'); if imptransfer==1 fprintf(fid,'Izsra:%8.0f+i(%-8.0f)A',real(IIzsra),imag(IIzsra)); fprintf(fid,'\tIzsrb:%8.0f+i(%-8.0f)A',real(IIzsrb),imag(IIzsrb)); fprintf(fid,'\tIzsrc:%8.0f+i(%-8.0f)A',real(IIzsrc),imag(IIzsrc)); fprintf(fid,'\n'); end % ********************************************************************** case{4}, % OVA (polar): 21, 67N, V(0,1,2),I(0,1,2), V(a,b,c),I(a,b,c) % --- ----------------- ----------------- % ********************************************************************** % ------------------------------------------------------------------ % C4.1 SALIDA 21, Ohm , polar % ------------------------------------------------------------------ % Sm,Sp,Rm,Rp % ------------------------------------------------------------------ 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'); [mod,arg]=mod_arg(SSm);fprintf(fid,'Sm=Sp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_arg(RRm);fprintf(fid,'Rm=Rp=%7.2f(%-5.1f)(Ohm)\n',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_argz(Sm_ZAB*zbase);fprintf(fid,'Sm_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sm_ZBC*zbase);fprintf(fid,'Sm_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sm_ZCA*zbase);fprintf(fid,'Sm_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sm_ZAG*zbase);fprintf(fid,'Sm_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sm_ZBG*zbase);fprintf(fid,'Sm_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sm_ZCG*zbase);fprintf(fid,'Sm_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_argz(Sp_ZAB*zbase);fprintf(fid,'Sp_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sp_ZBC*zbase);fprintf(fid,'Sp_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sp_ZCA*zbase);fprintf(fid,'Sp_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sp_ZAG*zbase);fprintf(fid,'Sp_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sp_ZBG*zbase);fprintf(fid,'Sp_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Sp_ZCG*zbase);fprintf(fid,'Sp_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_argz(Rm_ZAB*zbase);fprintf(fid,'Rm_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rm_ZBC*zbase);fprintf(fid,'Rm_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rm_ZCA*zbase);fprintf(fid,'Rm_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rm_ZAG*zbase);fprintf(fid,'Rm_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rm_ZBG*zbase);fprintf(fid,'Rm_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rm_ZCG*zbase);fprintf(fid,'Rm_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_argz(Rp_ZAB*zbase);fprintf(fid,'Rp_ZAB=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rp_ZBC*zbase);fprintf(fid,'Rp_ZBC=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rp_ZCA*zbase);fprintf(fid,'Rp_ZCA=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rp_ZAG*zbase);fprintf(fid,'Rp_ZAG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rp_ZBG*zbase);fprintf(fid,'Rp_ZBG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); [mod,arg]=mod_argz(Rp_ZCG*zbase);fprintf(fid,'Rp_ZCG=%7.2f(%-5.1f)(Ohm)\n',mod,arg); fprintf(fid,'\n'); % ----------------------------------------------------------------------- % C4.2 SALIDA 67N Volt,Amp polar % ----------------------------------------------------------------------- % Sm,Sp,Rm,Rp % ----------------------------------------------------------------------- fprintf(fid,'\n'); fprintf(fid,'SALIDA 67N, Volt, Amp, polar\n'); fprintf(fid,'\n'); fprintf(fid,'Sm_67N:\tVspol=%-8.0f(V)\tIsmop=%-8.0f(A)\tAngDirSm=%-5.1f\n',VVspol,IIsmop,AngDirSm); fprintf(fid,'Sp_67N:\tVspol=%-8.0f(V)\tIspop=%-8.0f(A)\tAngDirSp=%-5.1f\n',VVspol,IIspop,AngDirSp); fprintf(fid,'Rm_67N:\tVrpol=%-8.0f(V)\tIrmop=%-8.0f(A)\tAngDirRm=%-5.1f\n',VVrpol,IIrmop,AngDirRm); fprintf(fid,'Rp_67N:\tVrpol=%-8.0f(V)\tIrpop=%-8.0f(A)\tAngDirRp=%-5.1f\n',VVrpol,IIrpop,AngDirRp); % ------------------------------------------------------------------ % C4.3 SALIDA V(0,1,2), I(0,1,2) Volt,Amp polar % ------------------------------------------------------------------ % Vm,Vp,Vs,Vr % Im,Ip,Is,Ir,Ism,Isp,Irm,Irp %,Izsr % ------------------------------------------------------------------ fprintf(fid,'\n'); fprintf(fid,'SALIDA V(0,1,2), I(0,1,2) pu 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'); [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'); [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'); %------------------------------------------------------------------------------- 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'); [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'); [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(IIsm0);fprintf(fid,'Ism0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsm1);fprintf(fid,'Ism1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsm2);fprintf(fid,'Ism2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIsp0);fprintf(fid,'Isp0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsp1);fprintf(fid,'Isp1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsp2);fprintf(fid,'Isp2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIrm0);fprintf(fid,'Irm0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrm1);fprintf(fid,'Irm1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrm2);fprintf(fid,'Irm2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIrp0);fprintf(fid,'Irp0 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrp1);fprintf(fid,'Irp1 =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrp2);fprintf(fid,'Irp2 =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); if imptransfer==1 [mod,arg]=mod_arg(IIzsr0);fprintf(fid,'Izsr0=%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIzsr1);fprintf(fid,'Izsr1=%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIzsr2);fprintf(fid,'Izsr2=%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); end % ------------------------------------------------------------------ % C4.4 SALIDA V(a,b,c), I(a,b,c) Volt,Amp polar % ------------------------------------------------------------------ % Vm,Vp,Vs,Vr % Im,Ip,Is,Ir,Ism,Isp,Irm,Irp %,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'); [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'); [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'); %------------------------------------------------------------------------------- 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'); [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'); [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\n'); [mod,arg]=mod_arg(IIsma);fprintf(fid,'Isma =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsmb);fprintf(fid,'Ismb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIsmc);fprintf(fid,'Ismc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIspa);fprintf(fid,'Ispa =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIspb);fprintf(fid,'Ispb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIspc);fprintf(fid,'Ispc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIrma);fprintf(fid,'Irma =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrmb);fprintf(fid,'Irmb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrmc);fprintf(fid,'Irmc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); [mod,arg]=mod_arg(IIrpa);fprintf(fid,'Irpa =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrpb);fprintf(fid,'Irpb =%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIrpc);fprintf(fid,'Irpc =%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); if imptransfer==1 [mod,arg]=mod_arg(IIzsra);fprintf(fid,'Izsra=%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIzsrb);fprintf(fid,'Izsrb=%8.0f(%-5.1f)(A)\t',mod,arg); [mod,arg]=mod_arg(IIzsrc);fprintf(fid,'Izsrc=%8.0f(%-5.1f)(A)',mod,arg); fprintf(fid,'\n'); end 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  >  PC126M  >  FC126.m