第二次.zip_matlab_平面阵列_方向图_立体方向图_阵列方向图

clear;clc; sita=meshgrid(eps:pi/180:pi); fai=meshgrid(eps:2*pi/180:2*pi)'; l=0.25; d=1; beta=0; m=1; r1=abs(cos(2*pi*l*cos(sita))./sin(sita)); r2=2.*cos(pi./4+0.75.*pi.*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); surf(x2,y2,z2);axis([-1,1,-1,1,-1,1]);title('平行二元半波阵子阵的阵因子');shading interp; subplot(2,2,3); surf(x3,y3,z3);axis([-1,1,-1,1,-1,1]);title('平行二元半波阵子阵立体方向图');shading interp;