// Programa da Prática 2 - Cálculo do campo
elétrico de uma LT
clear
clc
close
// definição de constantes para o programa
eps_0 = 8.854D-12
ro_L = 10.13D-6
// coordenadas de cada um dos condutores
xa = -10
ya = 0
xb = 0
yb = 0
xc = 10
yc = 0
// Variáveis auxiliares para as cores dos gráficos e
legendas
cont_graf = 0;
cor1 = 0;
cor2 = 0;
cor3 = 0;
// Mensagem na tela
mprintf('\n0 - Sair \n')
mprintf('1 - Inserir dados \n')
mprintf('2 - Obter perfil de campo elétrico \n')
// Leitura da opção selecionada
opcao = input("Digite a opção desejada: ")
while opcao<>0
clc
select opcao
// Inserir dados
case 1 then
hp = input("Digite a altura (em metros) para a
análise: ")
while (hp<0)
mprintf('A altura deve ser maior do que
zero.')
hp = input("Digite a altura (em metros) para
a análise: ")
end
yp = -12+hp;
// Obter perfil de campo elétrico
case 2 then
cont = 1;
for xp=-20:0.1:20,
// Campo elétrico devido a linha A
r = sqrt((xp-xa)^2+(yp-ya)^2)
Ea_x = ro_L/(2*%pi*eps_0*r^2)*(xp-xa)
Ea_y = ro_L/(2*%pi*eps_0*r^2)*(yp-ya)
Ea = Ea_x + %i*Ea_y
vet_Eax(cont) = Ea_x;
vet_Eay(cont) = Ea_y;
Ea_rot = Ea
vet_Ea(cont) = Ea;
vet_Ea_rot(cont) = Ea_rot;
// Campo elétrico devido a linha B
r = sqrt((xp-xb)^2+(yp-yb)^2)
Eb_x = ro_L/(2*%pi*eps_0*r^2)*(xp-xb)
Eb_y = ro_L/(2*%pi*eps_0*r^2)*(yp-yb)
Eb = Eb_x + %i*Eb_y
vet_Ebx(cont) = Eb_x;
vet_Eby(cont) = Eb_y;
Eb_rot = Eb*exp(%i*(120)/180*%pi) //
Defasagem de 120 graus
vet_Eb(cont) = Eb;
vet_Eb_rot(cont) = Eb_rot;
// Campo elétrico devido a linha C
r = sqrt((xp-xc)^2+(yp-yc)^2)
Ec_x = ro_L/(2*%pi*eps_0*r^2)*(xp-xc)
Ec_y = ro_L/(2*%pi*eps_0*r^2)*(yp-yc)
Ec = Ec_x + %i*Ec_y
vet_Ecx(cont) = Ec_x;
vet_Ecy(cont) = Ec_y;
Ec_rot = Ec*exp(%i*(-120)/180*%pi) //
Defasagem de -120 graus
vet_Ec(cont) = Ec;
vet_Ec_rot(cont) = Ec_rot;
// Cálculo do campo elétrico resultante
E_res = Ea_rot + Eb_rot + Ec_rot
// Magnitude do campo elétrico
Mag_E_res = abs(E_res)
vet_xp(cont) = xp;
vet_mag_E(cont) = Mag_E_res
cont = cont + 1;
end
cor = [cor3 cor2 cor1];
plot(vet_xp,vet_mag_E,'color',cor)
xgrid
cont_graf = cont_graf + 1;
aux_leg(1,cont_graf) = strcat(string(hp)+" m");
legend(aux_leg)
// Código para alterar a cor do próximo gráfico
cor1 = cor1 + 1;
if cor1==2 then
cor1 = 0;
cor2 = cor2 + 1;
if cor2==2 then
cor2 = 0;
cor3 = cor3 + 1;
if cor3==2 then
cor3 = 1;
end
end
end
else
mprintf('Digite uma opção válida!')
end
mprintf('\n0 - Sair \n')
mprintf('1 - Inserir dados \n')
mprintf('2 - Obter perfil de campo elétrico \n')
opcao = input("Digite a opção desejada: ")
end
