platecomp.m

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

function AE=platecomp(the,phi)
k=2*pi;
Fac=[0,0,0;0,1,0;0,0,0];
Fac=[1,1,1;1,1,1;1,1,1];
%Fac=[1,2,1;2,3,2;1,2,1];
w=0.8040^2;
Fac=[w/8,w/4,w/8;w/4,-0.8+37/500,w/4;w/8,w/4,w/8];
Fac=[1,2,1;2,-8.9849,2;1,2,1];
%Fac=[1,2,1;2,-1.1231,2;1,2,1];
%Fac=[1,i,1;i,3,i;1,i,1];
E=zeros(size(phi));
a=sin(the).*cos(phi);
b=sin(the).*sin(phi);
for row=1:3
    for col=1:3
        x1=-5/12+1/6*col;
        y1=-5/12+1/6*row;
        x2=-1/4+1/6*col;
        y2=-1/4+1/6*row;
        a1=cos(k*x2.*a)-cos(k*x1.*a); 
        a2=sin(k*x2.*a)-sin(k*x1.*a);
        a3=cos(k*y1.*b)-cos(k*y1.*b);
        a4=sin(k*y2.*b)-sin(k*y1.*b);
        c=(a1+i*a2).*(a3+i*a4);
    E=E+i*c.*(1+cos(the))./(k*a)./(k*b)*Fac(row,col);
    end
end
    AE=abs(E)/max(max(abs(E)));