Matlab - Grafica Matlab

Hola

Un problemita con un for y una grafica que resulta de la siguiente codificacion:

function CAMPOS_ E

volt=input('Ingrese la tensión de la Linea de Transmisión (kV): ');
for u=1:3
dfases(u)=input('Ingrese posición x del conductor fase %u (m): ');
h(u)=input('Ingrese altura del conductor fase %u (m): ');
dfasesq(u)=dfases(u);
hq(u)=-h(u);
end
haz=input('Ingrese número de haz: ');
dsc=input('Ingrese diámetro de subconductor (cm): ');
sk=input('Ingrese distancia entre subconductores (cm): ');
xm=input('Ingrese coordenada x del punto de mediciòn M (m): ');
hm=input('Ingrese coordenada y del punto de mediciòn M (m): ');

%Ubicación de los conductores de fases a, b, y c.
%Fase a
fprintf('\n\nLas coordenadas del conductor fase a (xa,ya) son: ');
fprintf('xa = %1.1f ', dfases(1));
fprintf('; ');
fprintf('ya = %1.1f ', h(1));
%Fase b
fprintf('\nLas coordenadas del conductor fase b (xb,yb) son: ');
fprintf('xb = %1.1f ', dfases(2));
fprintf('; ');
fprintf('yb = %1.1f ', h(2));
%Fase c
fprintf('\nLas coordenadas del conductor fase c (xc,yc) son: ');
fprintf('xc = %1.1f ', dfases(3));
fprintf('; ');
fprintf('yc = %1.1f ', h(3));

%Ubicación de las imágenes de los conductores de fases a´, b´, y c´.
%Fase a´
fprintf('\n\nLas coordenadas de la imagen del conductor fase a´ (xa´,ya´) son: ');
fprintf('xa´ = %1.1f ', dfasesq(1));
fprintf('; ');
fprintf('ya´ = %1.1f ', -h(1));
%Fase b´
fprintf('\nLas coordenadas de la imagen del conductor fase b´ (xb´,yb´) son: ');
fprintf('xb´ = %1.1f ', dfasesq(2));
fprintf('; ');
fprintf('yb´ = %1.1f ', -h(1));
%Fase c
fprintf('\nLas coordenadas de la imagen del conductor fase c´ (xc´,yc´) son: ');
fprintf('xc´ = %1.1f ', dfasesq(3));
fprintf('; ');
fprintf('yc´ = %1.1f ', -h(3));

%Diámetro de un solo conductor
db=(sk*10^(-2))/sin(pi/haz);
fprintf('\n');
fprintf('\n\nEl diámetro del haz es (m): %6.2f\n', db);

%Diámetro equivalente del conductor
deq=db*((haz*dsc*10^(-2)/db)^(1/haz));
fprintf('\n');
fprintf('El diámetro equivalente del conductor es (m): %6.2f\n', deq);

e=8.854*10^(-12);

for k=1:3
for j=1:3
if j==k
p1(j,k)=log(4*h(k)/deq)/(2*pi*e);
end
end
end

for k=1:3
for j=1:3
if j~=k
p2(j,k)=log(sqrt(((dfases(j)-dfases(k))^2)+(((h(j)-h(k))^2)))/sqrt(((dfasesq(j)-dfases(k))^2)+(((hq(j)-h(k))^2))))/(2*pi*e);
end
end
end

fprintf('\n\nLa Matriz de Coeficientes Potenciales (P) es: ');
P=p1-p2

fprintf('\nLa Matriz de Capacitancias (C) es: ');
C=inv(P)

fprintf('\nLa Matriz de Tensión V es: ');
Vra=-0.5*volt/sqrt(3); Vrb=1*volt/sqrt(3); Vrc=-0.5*volt/sqrt(3);
Vr=[Vra; Vrb; Vrc];

Via=0.5*volt*i; Vib=0*volt*i; Vic=-0.5*volt*i;
Vi=[Via; Vib; Vic];

V=Vr+Vi

fprintf('\nLa Matriz de Cargas Q es: ');
Qr=1000*C*Vr;
Qi=1000*C*Vi;
Q=1000*C*V

%Componente horizontal
for j=1:1
for k=1:3
XH1(k)=((xm-dfases(k))/((dfases(k)-xm)^(2)+(h(k)-hm)^(2)))-((xm-dfases(k))/((dfases(k)-xm)^(2)+(h(k)+hm)^(2)));
f1(j,k)=XH1(k)/(2*pi*e);
end
end
Erx=f1*Qr;
Eix=f1*Qi;

%Componente vertical
for j=1:1
for k=1:3
XH2(k)=((hm-h(k))/((dfases(k)-xm)^(2)+(h(k)-hm)^(2)))-((hm+h(k))/((dfases(k)-xm)^(2)+(h(k)+hm)^(2)));
f2(j,k)=XH2(k)/(2*pi*e);
end
end
Ery=f2*Qr;
Eiy=f2*Qi;

fprintf('\nLas componentes horizontal y vertical del Campo Electrico (V/m) son: ');
Ex=f1*Q
Ey=f2*Q

fprintf('\nEl valor rms del Campo Electrico es (V/m): ');
Erms=sqrt(Erx^(2)+(Eix*(-i))^(2)+Ery^(2)+(Eiy*(-i))^(2))

for j=1:1
for k=1:3
vsp(j,k)=log(sqrt((dfases(k)-xm)^(2)+(h(k)-hm)^(2))/sqrt((dfases(k)-xm)^(2)+(h(k)+hm)^(2)))/(2*pi*e);
end
end
Vspr=vsp*Qr;
Vspi=vsp*Qi;

fprintf('\nEl valor del Potencial Electrico es (V): ');
Vsp=vsp*Q

fprintf('\nEl valor rms del Potencial Electrico es (V): ');
Vrms=sqrt(Vspr^(2)+(Vspi*(-i))^(2))

%GRAFICA DE E
syms x
x=-20:1:20;
Qsi=1000*C*Vi*(i)^(-1);

for j=1:1
for k=1:3
F1(j,k)=(((x-dfases(k))/((dfases(k)-x)^(2)+(h(k)-hm)^(2)))-((x-dfases(k))/((dfases(k)-x)^(2)+(h(k)+hm)^(2))))/(2*pi*e);
end
end

Erx=F1*Qr;A=Erx^(2);
Eix=F1*Qsi;B=Eix^(2);

Eox=(A+B)^(0.5);

for j=1:1
for k=1:3
Fy(j,k)=(((hm-h(k))/((dfases(k)-x)^(2)+(h(k)-hm)^(2)))-((hm+h(k))/((dfases(k)-x)^(2)+(h(k)+hm)^(2))))/(2*pi*e);
end
end

Ery=Fy*Qr;C=Ery^(2);
Eiy=Fy*Qsi;D=Ery^(2);
Eoy=sqrt(C^(2)+D^(2));

E=sqrt(Eox^(2)+Eoy^(2));

plot(x,y)
grid

Les explico en donde esta el problema:

for u=1:3 (al inicio de la codificacion)
deberia salir:
Ingrese posición x del conductor fase 1 (m):
pero sale:
Ingrese posición x del conductor fase %u (m):

Grafica
F1 y F2 son matrices 1x3 y Qr con Qsi son matrices 3x1. la operacion de estos da Erx Eix Eox, Ery Eiy Eoy, finalmente E que depende de x. Luego quiero graficar E vs x... pero nada


Los datos a ingresar son:
Ingrese la tensión de la Linea de Transmisión (kV): 525
Ingrese posición x del conductor fase %u (m): -10
Ingrese altura del conductor fase %u (m): 10.6
Ingrese posición x del conductor fase %u (m): 0
Ingrese altura del conductor fase %u (m): 10.6
Ingrese posición x del conductor fase %u (m): 10
Ingrese altura del conductor fase %u (m): 10.6
Ingrese número de haz: 3
Ingrese diámetro de subconductor (cm): 3.3
Ingrese distancia entre subconductores (cm): 45
Ingrese coordenada x del punto de mediciòn M (m): 20
Ingrese coordenada y del punto de mediciòn M (m): 2

Gracias.

Estimado,

Lograste solucionar el problema con el cálculo campo eléctrico. Ando en busca de un codigo como el tuyo para este tipo de cálculos y para calcular campos magneticos igual. Espero me puedas ayudar.

Saludos.

CLM