Matlab - Grafica en 3D

F=I1*exp(((x1)*sin(TETA).*cos(PHI)))*i

clear all

close all

TETA=pi/4;

x1=rand;

I1=rand;

PHI=pi;

%El codigo de la radiación es la siguiente y si lo muestra correctamente:

E=(sin(TETA).*cos(cos(TETA))./(sin(TETA)^2));

%y al multiplicar por F es decir por:

F=I1*exp(((x1)*sin(TETA).*cos(PHI)))*i;

%me mustra el error

%Al preparar el angulo de elevacion y el azimut tengo

phi=(0:.01:1)*2*pi;%fi 

teta=(0:.01:1)*pi;%theta

[PHI,TETA]=meshgrid(phi,teta);

%y para graficar::

PR=rand(size(TETA));

%cambio de coordenadas esfericas a rectangulares para poder graficar 

X=PR.*sin(TETA).*cos(PHI);

Y=PR.*sin(TETA).*sin(PHI);

Z=PR.*cos(TETA);

% grafica en 3d 

surface(X,Y,Z.*abs(E))

grid on

view(8,18)

>> radiacion

clc;  clear;  close all

TETA=pi/4;

PHI=pi;

phi=(0:.01:1)*2*pi;%fi 

teta=(0:.01:1)*pi;%theta

[PHI,TETA]=meshgrid(phi,teta);

N=10;

I1=0.0116-0.0044i;  x1=0;

I2=0.0040+0.0010i;  x2=0.5;

I3=-0.0011+0.0002i; x3=1;

I4=0.0005-0.0002i;  x4=1.5;

I5=-0.0003+0.0001i; x5=2;

I6=0.0003-0.0001i;  x6=2.5;

I7=-0.0002+0.0001i; x7=3;

I8=0.0002-0.0001i;  x8=3.5;

I9=-0.0002+0.0000i; x9=4;

I10=0.0002-0.0000i; x10=4.5;

F=I1*exp(((x1)*sin(TETA).*cos(PHI)))*i+I2*exp(((x2)*sin(TETA).*cos(PHI)))*i+I3*exp(((x3)*sin(TETA).*cos(PHI)))*i+I4*exp(((x4)*sin(TETA).*cos(PHI)))*i+I5*exp(((x5)*sin(TETA).*cos(PHI)))*i+I6*exp(((x6)*sin(TETA).*cos(PHI)))*i+I7*exp(((x7)*sin(TETA).*cos(PHI)))*i+I8*exp(((x8)*sin(TETA).*cos(PHI)))*i+I9*exp(((x9)*sin(TETA).*cos(PHI)))*i+I10*exp(((x10)*sin(TETA).*cos(PHI)))*i;

%creacion del campo E en el espacio 

E=(sin(TETA).*cos(cos(TETA))./(sin(TETA)^2));

PR=F.*E;

%y para graficar::

PR=rand(size(TETA));

%cambio de coordenadas esfericas a rectangulares para poder graficar 

X=PR.*sin(TETA).*cos(PHI);

Y=PR.*sin(TETA).*sin(PHI);

Z=PR.*cos(TETA);

% grafica en 3d 

surface(X,Y,Z.*abs(PR))

grid on

view(8,18)

xlabel('X')

ylabel('Y')

zlabel('Z')