1.m

计算9格均匀电流面元阵方向图的matlab程序

E=1;
theta=0;
phi=0:0.02*pi:pi/2;
[x,y,z]=sph2cart(theta,phi,E);
plot3(x,y,z,'Color','black');
hold all
theta=pi/2;
[u,v,w]=sph2cart(theta,phi,E);
plot3(u,v,w,'Color','black');
hold all
theta=pi/3;
[u,v,w]=sph2cart(theta,phi,E);
plot3(u,v,w,'Color','black');
hold all
the=0:0.02*pi:pi/2;
ph=0;
[u,v,w]=sph2cart(the,ph,E);
plot3(u,v,0*u,'Color','black');
ph=pi/3;
[u,v,w]=sph2cart(the,ph,E);
plot3(u,v,w*ones(size(v)),'Color','black');
line([0,1],[0,0],[0,0],'Color','black');
line([0,0],[0,1],[0,0],'Color','black');
line([0,0],[0,0],[0,1],'Color','black');
[p1,p2,p3]=sph2cart(pi/3,pi/3,1)
line([p1,p1-0.2*sin(pi/3)],[p2,p2+0.2*cos(pi/3)],[p3,p3],'Color','black');
line([p1,p1+0.2*cos(pi/3)],[p2,p2+0.2*sin(pi/3)],[p3,p3-0.2*cos(pi/3)],'Color','black');
line([0,1.12*p1],[0,1.12*p2],[0,1.12*p3],'Color','black');
line([0,cos(pi/3)],[0,sin(pi/3)],[0,0],'Color','black');
axis off