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

Tarea 22 Matlab

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

MATLAB:

Ejercicio 6.12
%%%6.12 Linea Corta/Media/Larga
clear all
clc
l=175;%longitud de linea

%Impedancia serie total por fase Z=zl


Z=35+140*i;
mZ=abs(Z);
aZ=angle(Z)*180/pi;
%Admintancia en paralelo total por fase al neutro Y=yl
Y=(930*(10^(-6)))*i;
mY=abs(Y);
aY=angle(Y)*180/pi;

%Impedancia por Unidad de longitud por fase


z=Z/l; % z=A3+i*(A3(Xa)+A4)
mz=abs(z);
az=angle(z)*180/pi;

%Admitancia Paralelo por Unidad por fase al neutro


y=Y/l;% y=j(1/(A3(Xd)+A5))
my=abs(y);
ay=angle(y)*180/pi;

%Constante de Propagacion y
CP=sqrt(y*z);

%Constante de propagacion total yl


yl=CP*l;
myl=abs(yl);
ayl=angle(yl)*180/pi;

%Impedancia Caracteristica de la Linea Zc=sqrt(z/y)


Zc=sqrt(z/y);
mZc=abs(Zc);
aZc=angle(Zc)*180/pi;

%% Linea Corta

%Corriente Ir
Cor=(40000)/(sqrt(3)*220*0.9);
AngCor=acosd(0.9);
Ir=104.94-50.821*i;

%%Linea Corta
Vr=220000/(sqrt(3));
Is1=Ir;
mIs1=abs(Ir);
aIs1=angle(Ir)*180/pi;

Vs1= Vr+Ir*Z; %%Vs=Vr + Ir*Z


mVs1=abs(Vs1);
aVs1=angle(Vs1)*180/pi;

%Voltaje Linea Linea


Vll1=(sqrt(3))*mVs1/1000;%en KV

%% LINEA MEDIA

A=((Z*Y)/2+1);
D=A;
B=Z;
C=Y*(1+((Z*Y)/4));

Vs2=A*Vr+B*Ir;
mVs2=abs(Vs2);
aVs2=angle(Vs2)*180/pi;
%Voltaje Linea Linea
Vll2=(sqrt(3))*mVs2/1000;%en KV

Is2=Vr*C+Ir*D;
mIs2=abs(Is2);
aIs2=angle(Is2)*180/pi;

%% LINEA LARGA

A1=cosh(yl);
B=Zc*sinh(yl);
C=sinh(yl)/Zc;
D=cosh(yl);

Vs3=A*Vr+B*Ir;
mVs3=abs(Vs3)/1000;%en kV
aVs3=angle(Vs3)*180/pi;
%Voltaje Linea Linea
Vll3=(sqrt(3))*mVs3;%en kV

Is3=Vr*C+Ir*D;
mIs3=abs(Is3);
aIs3=angle(Is3)*180/pi;

%% RESPUESTAS

disp('A) Aproximacion de Linea Corta')


disp(Vs1)
disp('Vs Magnitud (en kV):')
t1=mVs1/1000;
disp(t1)
disp('Vs Angulo:')
disp(aVs1)
disp('|Vs|')
disp(Vll1)

disp('B) Aproximacion de Linea Media')


disp(Vs2)
disp('Vs Magnitud (en kV):')
t2=mVs2/1000;
disp(t2)
disp('Vs Angulo:')
disp(aVs2)
disp('|Vs|')
disp(Vll2)

disp('C) Aproximacion de Linea Larga')


disp(Vs3)
disp('Vs Magnitud (en kV):')
disp(mVs3)
disp('Vs Angulo:')
disp(aVs3)
disp('|Vs|')
disp(Vll3)
Ejercicio 6.20
clear all
clc
disp('Ejercicio 6.20')
%%Ejercicio 6_20
l=31.0686;%50km en millas
distancia1=16.4042;%5m en pies

%% A
%ya que sabemos que X = 4*pi*f*10^(-7)*(Dm/Ds) tenemos lo siguiente

Xp=0.0754*log(20.67/0.0217);% X para Partridge 0.752 xq queremos en ohms/km


multiplicamos x*1000
Xo=0.0754*log(20.67/0.0284);% X para Osprey
vs=1;
vr=1;
Prp=((vs*vr)/(Xp))*cos(45);% Potencia partridge
Pro=((vs*vr)/(Xo))*cos(45);% Potencia Osprey
Pi=(Pro-Prp)/Prp*100;%Porcentaje de Incremento= (Valor Final-Valor
Inicial)/Valor Inicial*100

%% B Un segundo conductor Partridge


%se coloca en un arreglo de dos conductores separados a 40cm
distancia2=1.31234;%40cm en pies
Ds=sqrt(distancia2*0.0217);% Ds=RMG=sqrt(rmgviejo*distancia)
Xn=0.0754*log(20.67/Ds);% X nuevo
Prn=((vs*vr)/(Xn))*cos(45);
Pi2=(Prn-Prp)/Prp*100;% porcentaje de incremento 2

%% C Se incrementa el voltaje un espaciamento entre los conductores de 8m


distancia3=26.2467;%8m en pies
Deq2=(26.2467*26.2467*(26.2467+26.2467))^(1/3);
Xn2=0.0754*log(Deq2/0.0217);
Prn2=((vs*vr)/(Xn2))*cos(45);
Pi3=abs((Prn2-Prp)/Prp)*100;% porcentaje de incremento 3
PiV=(230-138)/138*100;% porcentaje de incremento del voltaje
Incremento=Pi3*PiV;

%% RESPUESTAS
disp('A)')
disp('Porcentaje en Incremento de Potencia:')
disp(Pi)
disp('B)')
disp('Porcentaje en Incremento de Potencia:')
disp(Pi2)
disp('C)')
disp('Porcentaje en Incremento de Potencia:')
disp(Incremento)
EJERCICIO 6.24
%% Ejericio 6.24
clear all
clc
disp('Ejericio 6.24')
l=300;% Distancia de la linea
%%La admitancia paralelo de la linea es
y=6.87*10^(-6);
%%para la totalidad de la linea
Bc=y*l;
%%para el 60% de compensacion
Bl=0.6*Bc;
%%Debido a la tabla A.6 tenemos que
A=1;
D=A;
B=0;
C=-j*Bl;

%% Respuesta
disp('A:')
disp(A)
disp('B:')
disp(B)
disp('C:')
disp(C)
disp('D:')
disp(D)

También podría gustarte