K1=abs(W1);
%----------------------
W=a+(beta.*d.*cos(t));
z3=(N2/2).*(W);
z4=(1/2).*(W);
W2=sin(z3)./(N2.*sin(z4)); % 非归一化的阵因子 K2
K2=abs(W2);
%-------------------------
W=a+(beta.*d.*cos(t));
z5=(N3/2).*(W);
z6=(1/2).*(W);
W3=sin(z5)./(N3.*sin(z6)); % 非归一化的阵因子 K3
K3=abs(W3);
%-------------------- 绘图函数
&gure(2)
subplot(131);polar(t,K1);xlabel('f=3GHz,N=4,d=1/2 波长,a= π/2');
subplot(132);polar(t,K2);xlabel('f=3GHz,N=8,d=1/2 波
长,a= π/2');
subplot(133);polar(t,K3);xlabel('f=3GHz,N=12,d=1/2 波
长,a= π/2'); %----------------------
%只有阵列单元方向角 a 改变的天线方向图
N4=10;d1=lambda/4;a1=0;a2=pi/2;a3=pi/2+pi/10;
A=a1+(beta.*d1.*cos(t));
x1=(N4/2).*(A);
x2=(1/2).*(A);
A1=sin(x1)./(N4.*sin(x2)); % 非归一化的阵因子 K4
K4=abs(A1);
%---------------------------
B=a2+
(beta.*d1.*cos(t));
y_1=(N4/2).*(B);
y_2=(1/2).*(B);
B1=sin(y_1)./(N4.*sin(y_2)); % 非归一化的阵因子 K5
K5=abs(B1);
%----------------------------
C=a3+(beta.*d1.*cos(t));
v1=(N4/2).*(C);
v2=(1/2).*(C);
C1=sin(v1)./(N4.*sin(v2)); % 非归一化的阵因子 K6
K6=abs(C1);
%-------------------------- 绘图函数
&gure(3)
subplot(131);polar(t,K4);xlabel('f=3GHz,N=10,d=1/4 波长,a=0');
subplot(132);polar(t,K5);xlabel('f=3GHz,N=10,d=1/4 波长,a= π/2');
subplot(133);polar(t,K6);xlabel('f=3GHz,N=10,d=1/4 波
长,a= π/2+ π/10');