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)

...