二元阵方向函数：一、半波振子：程序如下（只修改部分参数）：E 面：a=linspace(0,2*pi); b=linspace(0,pi);f1=abs(cos(pi/2.*cos(a))/sin(a));f2=abs(cos(1/2.*(2*pi*sin(a)+pi)));f=f1.*f2;subplot;polar(a,f);title('半波对称阵子E面,d=λ,ζ=π');H面：a=linspace(0,2*pi); b=linspace(0,pi);f=abs(cos(1/2.*(2*pi*cos(a)+pi)));subplot(221);polar(a,f); title('半波对称阵子H面,d=λ,ζ=π');1、d=λ/4 ：E 面：H 面：2、d=λ/2E 面：H 面：3、d=3λ/4E 面：H 面：4、d=λE 面：H 面：二、全波阵子程序如下：E 面：a=linspace(0,2*pi);b=linspace(0,pi);f1=abs((cos(pi*cos(a))+1)/sin(a));f2=abs(cos(1/2.*(pi/2.*sin(a))));f=f1.*f2subplot(221);polar(a,f); title('全波对称阵子E面,d=λ/4,ζ=0');H面：a=linspace(0,2*pi); b=linspace(0,pi);f=2*abs(cos(1/2.*(pi/2.*cos(a))));subplot(221);polar(a,f); title('全波对称阵子H面,d=λ/4,ζ=0');1、d=λ/4E面：H面： 2、d=λ/2E面：H面：3、d=3λ/4E面：H面：4、d=λE面：H面：。