Matlab画平行二元半波振子阵立体方向图.docxVIP

1. 平行二元半波振子阵立体方向图的源代码： clear;clc; sita=meshgrid(eps:pi/180:pi); fai=meshgrid(eps:2*pi/180:2*pi); l=0.25;%对称振子的长度为0.25 d=1;%二元阵的间隔距离此处取一倍波长 beta=0;%电流初始相位差为0 m=1;%电流的幅度比为1 r1=abs(cos(2*pi*l*cos(sita))-cos(2*pi*l))./abs(sin(sita)); r2=sqrt(1+m*m+2*m*cos(beta+2*pi*d*sin(sita).*sin(fai))); r3=r1.*r2; r1max=max(max(r1));r2max=max(max(r2));r3max=max(max(r3)); [x1,y1,z1]=sph2cart(fai,pi/2-sita,r1/r1max); [x2,y2,z2]=sph2cart(fai,pi/2-sita,r2/r2max); [x3,y3,z3]=sph2cart(fai,pi/2-sita,r3/r3max); subplot(2,2,1); surf(x1,y1,z1);axis([-1,1,-1,1,-1,1]);title(平行二元半波振子阵的元因子);shading interp; subplot(2,2,2); su