cell2.m

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

[theta,phi]=meshgrid(0:0.02*pi:2*pi,0:pi/12:1.1*pi);
the=pi/2-theta;
k=2*pi;
a=sin(the).*cos(phi);
b=sin(the).*sin(phi);
%Fac=[1,2,1;2,3,2;1,2,1];
%Fac=[0,0,0;0,1,0;0,0,0];
Fac=[1,2,1;2,-8.9849,2;1,2,1];
PF=exp(i/3*pi*(-a-b))*Fac(1,1)+...
    exp(i/3*pi*(-b))*Fac(1,2)+...
    exp(i/3*pi*(a-b))*Fac(1,3)+...
    exp(i/3*pi*(-a))*Fac(2,1)+...
    ones(size(a))*Fac(2,2)+...
    exp(i/3*pi*(a))*Fac(2,3)+...
    exp(i/3*pi*(-a+b))*Fac(3,1)+...
    exp(i/3*pi*(b))*Fac(3,2)+...
    exp(i/3*pi*(a+b))*Fac(3,3);
E=(1+cos(the)).*(4/k/k).*sin(pi/6*a)./a.*sin(pi/6*b)./b.*PF;
    AE=abs(E)/max(max(abs(E)));
    [x,y,z]=sph2cart(phi,theta,AE);
 mesh(x,y,z,'EdgeColor','black');
 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;
% a=sin(the).*cos(phi);
% b=sin(the).*sin(phi);
% E=(1+cos(the)).*(4/k/k).*sin(pi/6*a)./a.*sin(pi/6*b)./b;
%     AE1=abs(E)/max(max(abs(E)));
figure;polar(theta(10,:),AE(10,:),'-k')