Location via proxy:   [ UP ]  
[Report a bug]   [Manage cookies]                

Codigo en Matlab de Lineas de Transmisio

Descargar como docx, pdf o txt
Descargar como docx, pdf o txt
Está en la página 1de 2

function [A,B,C,D,V_S,I_S,RT,Ps,Qs] = LineasT(z,y,l,s,V,tn,fp)

%Linea corta
if l<80
Z=z*l;
A=1;
B=Z;
C=0;
D=1;
end

%Linea media
if 80<l<240
Y=y*l;
Z=z*l;
A=1+(Z*Y*0.5);
B=Z;
C=Y*(1+Z*Y*0.25);
D=1+(Z*Y*0.5);
end

%Linea larga
if l>240
gamma=sqrt(z*y);
Zc=sqrt(z/y);
Z_prima=Zc*sinh(gamma*l);
Y_prima=(1/Zc)*tanh(gamma*l*0.5);

A=1+(Z_prima*Y_prima);
B=Z_prima;
C=Y_prima*(1+Z_prima*Y_prima*0.5);
D=1+(Z_prima*(Y_prima));
end

V_r=tn*V;
V_R=V_r/sqrt(3);

%Atraso
if fp<0
theta=acos(abs(fp));
I_r=s/(sqrt(3)*V_r*abs(fp));
I_R=(I_r*cos(-theta))+(i*I_r*sin(-theta));
V_S=A*V_R+B*I_R;
I_S=C*V_R+D*I_R;
end

%Adelanto
if fp>0
theta=acos(fp);
I_r=s/(sqrt(3)*V_r*fp);
I_R=(I_r*cos(theta))+(i*I_r*sin(theta));
V_S=A*V_R+B*I_R;
I_S=C*V_R+D*I_R;
end

%Voltaje en vacio o sin carga


V_0=abs(V_S)/abs(A);
RT=(abs(V_0-V_R)/V_R)*100;

%Potencia activa y reactiva, en el suministro


Ps=sqrt(3)*(sqrt(3)*abs(V_S))*abs(I_S)*cos(angle(V_S)-angle(I_S));
Qs=sqrt(3)*(sqrt(3)*abs(V_S))*abs(I_S)*sin(angle(V_S)-angle(I_S));
end

También podría gustarte