%半波振子元阵子的图像
a = [0:0.1:2*pi];
b = pi*cos(a)/2;
c = abs(cos(b)./sin(a));
subplot(1,3,1);
polar(a,c);
  
%阵因子图像
subplot(1,3,2);
d = abs(cos(pi*cos(a))*2);
polar(a,d);
  
%方向图的 E 图像
subplot(1,3,3);
e = c.*d;
polar(a,e);