今天来画3D plot
% Load the data
load latticeExample
xx = x(:);yy = y(:);zz = z(:);
% locate the non-zero points
a = find(T~=0);
% plot the non-zero points using a scatter plot
% use the values of T to represent both color and size of symbols
scatter3(xx(a),yy(a),zz(a),1000*T(a),T(a),'filled');
% set the view
view(3);
接下来STREAM SLICE；
%% Generate the data
% define the grid
x_lim=.8; y_lim=.8; z_lim=.8; step=.01;
x=[-x_lim:step:x_lim]; y=[-y_lim:step:y_lim]; z=[-z_lim:step:z_lim];
[x y z]=meshgrid(x,y,z);r=sqrt(x.^2+y.^2+z.^2);
%generate the data
cos_theta=z./r; sin_theta=sqrt(1-(cos_theta).^2); % These are the angles made by the x,y,z axes
cos_phi=x./sqrt(x.^2+y.^2); sin_phi=y./sqrt(x.^2+y.^2);
V_dip_ang=cos_theta; %only the angular dependence of the radiation pattern is considered
theta_hat_x=cos_theta.*cos_phi; theta_hat_y=cos_theta.*sin_phi; theta_hat_z=-sin_theta; % theta_hat direction resolved in x,y, and z directions
r_hat_x=sin_theta.*cos_phi; r_hat_y=sin_theta.*sin_phi; r_hat_z=cos_theta; %r_hat direction resolved in x,y, and z directions
E_ang_x=2*cos_theta.*r_hat_x + sin_theta.*theta_hat_x; %x component of the electric field (only the angular pattern)
E_ang_y=2*cos_theta.*r_hat_y + sin_theta.*theta_hat_y; %y component of the electric field
E_ang_z=2*cos_theta.*r_hat_z + sin_theta.*theta_hat_z; %z component of the electric field
E_ang=sqrt(E_ang_x.^2+E_ang_y.^2+E_ang_z.^2); % The magnitude of the electric field
%% Plotting the absolute value of the potential and the magnitude of the electric field distribution on a spehere
% plot the data
streamslice(x,y,z,E_ang_x,E_ang_y,E_ang_z,[-.3 .1],[],[]); % This is used for vector data
% set the view
campos([-3,-15,5])
box on
% set the background and tick colors
set(gca,'Color',[.8,.8,.8],'XColor','black', ...
'YColor','black','ZColor','black')
%annotate
title({'Visualizing the electric field pattern of a ',...
'dipole antenna using \color{red}streamslice'},'Fontsize',14)
今天时间有点紧，，就先这些吧。。
