半波对称振子方向图_Matlab画平行二元半波振子阵立体方向图

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);

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;

2.

平行二元半波振子阵立体方向图：