basic4.m

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

[theta,phi]=meshgrid(0:0.02*pi:2*pi,0:pi/12:1.1*pi);
the=pi/2-theta;
AE=platecomp(the,phi);
[x,y,z]=sph2cart(phi,theta,AE);
axis off
axis equal
[theta,phi]=meshgrid(0:0.01*pi:2*pi,0:pi/12:1.1*pi);
the=pi/2-theta;
phi=ones(size(phi))*pi;
AE1=platecomp(the,phi);
phi=ones(size(phi))*pi/4;
AE2=platecomp(the,phi);
% subplot(1,3,1);
% mesh(x,y,z,'EdgeColor','black'); axis off; axis equal;
% subplot(1,3,2); 
% polar(theta(1,:),AE1(1,:),'-k')
% subplot(1,3,3); 
% polar(theta(1,:),AE2(1,:),'-k')

surf(x,y,z); colormap(winter);axis off; axis equal;
 figure;polar(theta(1,:),AE1(1,:),'-k')
% figure;polar(theta(1,:),AE2(1,:),'-k')
zz=inline('platecomp(the,phi).^2.*sin(the)','the','phi');