El tonillo sin fin

% matriz de transformacion

a = (-1/2) + (sqrt(3)/2)*(1i)

A = [1 1 1;

1 a^2 a;

1 a a^2]

% Conductores en haz

Haz = 2;

d = 50e-2;

% Terreno

P = 100;

f = 60;

% permeabilidad

U = 4*pi*1e^-7; %H/m

W = (U/(2*pi))*(2*pi*60); % ohm/m

w = W / 0.001 %Ohm/Km

% permitividad

Eo = 1e-9/(36*pi); %f/m

Eo = (Eo)/0.001; % F/km

K = 1 / (2*pi*Eo) % Km/F

% DME Y RMG

RMG_fabricante = 9.5e-3; %m

RMG = sqrt(RMG_fabricante*d);

RMG_p = RMG*exp(-1/4);

DME = (5*5*10)^(1/3);

%Altura de Torre

y = 10;

h = 2*(y);

%D(p,f)

De = 2160*sqrt(p/f); %ft

% Resistencia de la linea

Ra = 0.3;

Ra = Ra/Haz;

% Resistencia de Carson

Re = 986.9*1e^-6*f;

%%%%%%% MODELO SERIE DE LA LINEA - TRANSPUESTA %%%%%

%% Equivalente por Fase

% Resistencia propia y mutua

Rs = Ra * Re;

Rm = Re

% Reactancia propia y mutua

Xs = W*log(De / RMG_p);

Xm = W*log(De / DME);

Resistencia_Fase = [Rs Rm Rm;

Rm Rs Rm;

Rm Rm Rs];

Reactancia_Fase = [Xs Xm Xm;

Xm Xs Rm;

Xm Xm Xs]*(1i);

%% Modelo por secuencia (0,1,2)

Resistencia_ secuencia=inv(A)*Resistencia_Fase*A;

Reactancia_ secuencia=inv(A)*Resistencia_Fase*A;

%% Modelo paralelo

HMG = (h*h*h)^(1/3); %distancia media equivalente entre fases y sus espejos

lab=sqrt((h)^2 + (5)^2)

lac=sqrt((h)^2 + (10)^2)

lbc=sqrt((h)^2 + (5)^2)

LMG = (lab*lac*lbc)^(1/3);

% -matriz de Potencia

Ps = K * log(HMG/RMG);

Pm = K * log(LMG/DME);

P = [Ps Pm Pm;

Pm Ps Pm;

Pm Pm Ps];

% Capacitancia de Fase

capacitancia _Fase = inv(P);

% Capacitancia de Secuencia

capacitancia _Fase = inv(P);

capacitancia_ secuencia=inv(A)*capacitancia _Fase*A;

%%%%%%%%%%%% LINEA NO TRANSPUESTA %%%%%%%%%%%%%%%

%% Modelo Serie

% Resistencia Propia y Mutua

Rs_nT = Ra + Re;

Rm_nT = Re;

Resistencia_fase_nt = [Rs_nT Rm_nT Rm_nT;

Rm_nT Rs_nT Rm_nT;

Rm_nT Rm_nT Rs_nT;]

% -reactancia inductiva

% Reactancia propia

Xs_nT = W*log(De / RMG_p);

% Reactancia mutua

D(1,2) = 5; D(1,3) = 10; D(2,3) = 5;

Xm_nT(1,2) = W*log (De / D(1,2));

Xm_nT(1,3) =

...