0:
2*pi./（n_phi-1）:
2*pi）;
t3=tehta;
%-------------只有参数N改变的天线方向3D图
M=a+（beta.*d.*cos（t3））;%----N1=4;N2=8;N3=12;
z_1=（N1/2）.*（M）;
z_2=（1/2）.*（M）;
M1=sin（z_1）./（N1.*sin（z_2））;%非归一化的阵因子K1
K_1=abs（M1）;
radio_1=K_1;
X1=radio_1.*sin（tehta）.*cos（phi）;
Y1=radio_1.*sin（tehta）.*sin（phi）;
Z1=radio_1.*cos（tehta）;
%-----------------------------------
M=a+（beta.*d.*cos（t3））;
z_3=（N2/2）.*（M）;
z_4=（1/2）.*（M）;
M2=sin（z_3）./（N2.*sin（z_4））;%非归一化的阵因子K2
K_2=abs（M2）;
radio_2=K_2;
X2=radio_2.*sin（tehta）.*cos（phi）;
Y2=radio_2.*sin（tehta）.*sin（phi）;
Z2=radio_2.*cos（tehta）;
%-------------------------------------
M=a+（beta.*d.*cos（t3））;
z_5=（N3/2）.*（M）;
z_6=（1/2）.*（M）;
M3=sin（z_5）./（N3.*sin（z_6））;%非归一化的阵因子K3
K_3=abs（M3）;
radio_3=K_3;
X3=radio_3.*sin（tehta）.*cos（phi）;
Y3=radio_3.*sin（tehta）.*sin（phi）;
Z3=radio_3.*cos（tehta）;
%------------------------------------3D绘图函数
figure（5）
surf（X1,Y1,Z1）;%三维绘图函数surf，采用伪彩色表示曲面的高度
camlightright
light
colorbar
axisimage
rotate3Don
title（'f=3GHz,N=4,d=1/2波长,a=π/2'）;
figure（6）
surf（X2,Y2,Z2）;
camlightright
light
colorbar
axisimage
rotate3Don
title（'f=3GHz,N=8,d=1/2波长,a=π/2'）;
figure（7）
surf（X3,Y3,Z3）
camlightright
light
colorbar
axisimage
rotate3Don
title（'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_3d=a1+（beta.*d1.*cos（t3））;
x_1=（N4/2）.*（A_3d）;
x_2=（1/2）.*（A_3d）;
A_1=sin（x_1）./（N4.*sin（x_2））;%非归一化的阵因子K4
K_4=abs（A_1）;
radio_4=K_4;
X4=radio_4.*sin（tehta）.*cos（phi）;
Y4=radio_4.*sin（tehta）.*sin（phi）;
Z4=radio_4.*cos（tehta）;
%-----------------------------------
B_3d=a2+（beta.*d1.*cos（t3））;
y_1_3d=（N4/2）.*（B_3d）;
y_2_3d=（1/2）.*（B_3d）;
B_1=sin（y_1_3d）./（N4.*sin（y_2_3d））;%非归一化的阵因子K5
K_5=abs（B_1）;
radio_5=K_5;
X5=radio_5.*sin（tehta）.*cos（phi）;
Y5=radio_5.*sin（tehta）.*sin（phi）;
Z5=radio_5.*cos（tehta）;
%------------------------------------
C_3d=a3+（beta.*d1.*cos（t3））;
v_1=（N4/2）.*（C_3d）;
v_2=（1/2）.*（C_3d）;
C_1=sin（v_1）./（N4.*sin（v_2））;%非归一化的阵因子K6
K_6=abs（C_1）;
radio_6=K_6;
X6=radio_6.*sin（tehta）.*cos（phi）;
Y6=radio_6.*sin（tehta）.*sin（phi）;
Z6=radio_6.*cos（tehta）;
%-----------------------------------
figure（8）
surf（X4,Y4,Z4）;%三维绘图函数surf，采用伪彩色表示曲面的高度
camlightright
light
colorbar
axisimage
rotate3Don
title（'f=3GHz,N=10,d=1/4波长,a=0'）;
figure（9）
surf（X5,Y5,Z5）;
camlightright
light
colorbar
axisimage
rotate3Don
title（'f=3GHz,N=10,d=1/4波长,a=π/2'）;
figure（10）
surf（X6,Y6,Z6）
camlightright
light
colorbar
axisimage
rotate3Don
title（'f=3GHz,N=10,d=1/4波长,a=π/2+π/10'）;
%---------------------------------------------
%-------------------只有阵列单元间隔d改变的天线方向3D图
N5=20;d2=lambda/4;d3=lambda/2;d4=0.7*lambda;a4=pi/2;
D_3d=a4+（beta.*d2.*cos（t3））;
p_1=（N5/2）.*（D_3d）;
p_2=（1/2）.*（D_3d）;
D_1=sin（p_1）./（N5.*sin（p_2））;%非归一化的阵因子K7
K_7=abs（D_1）;
radio_7=K_7;
X7=radio_7.*sin（tehta）.*cos（phi）;
Y7=radio_7.*sin（tehta）.*sin（phi）;
Z7=radio_7.*cos（tehta）;
%---------------------------------------
E_3d=a4+（beta.*d3.*cos（t3））;
q_1=（N5/2）.*（E_3d）;
q_2=（1/2）.*（E_3d）;
E_1=sin（q_1）./（N5.*sin（q_2））;%非归一化的阵因子K8
K_8=abs（E_1）;
radio_8=K_8;
X8=radio_8.*sin（tehta）.*cos（phi）;
Y8=radio_8.*sin（tehta）.*sin（phi）;
Z8=radio_8.*cos（tehta）;
%------------------------------------------
F_3d=a4+（beta.*d4.*cos（t3））;
r_1_3d=（N5/2）.*（F_3d）;
r_2_3d=（1/2）.*（F_3d）;
F_1=sin（r_1_3d）./（N5.*sin（r_2_3d））;%非归一化的阵因子K9
K_9=abs（F_1）;
radio_9=K_9;
X9=radio_9.*sin（tehta）.*cos（phi）;
Y9=radio_9.*sin（tehta）.*sin（phi）;
Z9=radio_9.*cos（tehta）;
%-----------------------------------
figure（11）
surf（X7,Y7,Z7）;%三维绘图函数surf，采用伪彩色表示曲面的高度
camlightright
light
colorbar
axisimage
rotate3Don
title（'f=3GHz,N=20,d=1/4波长,a=π/2'）;
figure（12）
surf（X8,Y8,Z8）;
camlightright
light
colorbar
axisimage
rotate3Don
title（'f=3GHz,N=20,d=1/2波长,a=π/2'）;
figure（13）
surf（X9,Y9,Z9）
camlightright
light
colorbar
axisimage
rotate3Don
title（'f=3GHz,N=20,d=0.7波长,a=π/2'）;
二、
%-----------------均匀直线阵列天线的应用之一：
边射阵
clc
clearall
f=3e10;%30GHz,厘米波
i=1;%天线电流值
lambda=（3e8）/f;%lambda=c/f波长
beta=2.*pi/lambda;
N=15;
t=0:
0.01:
2*pi;
d1=lambda/4;%没有栅瓣效应的边射阵，即间隔d<波长
W1=beta.*d1.*cos（t）;%定义kdcos（方向角）
z1=（N/2）.*W1;
z2=（1/2）.*W1;
F1=sin（z1）./（N.*sin（z2））;
K1=abs（F1）;
d2=lambda*1.5;%有栅瓣现象的边射阵，即间隔d>波长
W2=beta.*d2.*cos（t）;%定义kdcos（方向角）
z3=（N/2）.*W2;
z4=（1/2）.*W2;
F2=sin（z3）./（N.*sin（z4））;
K2=abs（F2）;
figure
（1）
subplot（121）;polar（t,K1）;title（'边射阵f=30GHz,N=15,d=1/4波长'）;
subplot（122）;polar（t,K2）;title（'边射阵（有栅瓣）f=30GHz,N=15,d=1.5倍波长'）;
三、
%-----------------均匀直线阵列天线的应用之二：
普通端射阵
clc
clearall
f=3e10;%30GHz,厘米波
i=1;%天线电流值
lambda=（3e8）/f;%lambda=c/f波长
beta=2.*pi/lambda;
N=15;
t=0:
0.01:
2*pi;
d1=lambda/4;%没有栅瓣效应的普通端射阵，即间隔d<1/2波长
W1=beta.*d1.*cos（t）;%定义kdcos（方向角）
z1=（（N/2）.*W1）+N/2*beta*d1;
z2=（（1/2）.*W1）+1/2*beta*d1;
F1=sin（z1）./（N.*sin（z2））;
K1=abs（F1）;
d2=lambda*0.7;%有栅瓣现象的普通端射阵，即间隔d>1/2波长
W2=beta.*d2.*cos（t）;%定义kdcos（方向角）
z3=（（N/2）.*W2）+N/2*beta*d2;
z4=（（1/2）.*W2）+1/2*beta*d2;
F2=sin（z3）./（N.*sin（z4））;
K2=abs（F2）;
figure
（2）
subplot（121）;polar（t,K1）;title（'普通端射阵f=30GHz,N=15,d=1/4波长'）;
subplot（122）;polar（t,K2）;title（'普通端射阵（有栅瓣）f=30GHz,N=15,d=0.7倍波长'）;
四
%-----------------均匀直线阵列天线的应用之三：
强方向性端射阵
clc
clearall
f=3e10;%30GHz,厘米波
i=1;%天线电流值
lambda=（3e8）/f;%lambda=c/f波长
beta=2.*pi/lambda;
N=15;
t=0:
0.01:
2*pi;
d1=lambda/4;%没有栅瓣效应的强方向性端射阵，即间隔d<（1/2波长）*（1-1/N）
W1=beta.*d1.*cos（t）;%定义kdcos（方向角）
z1=（（N/2）.*W1）+N/2*（beta*d1+pi/N）;
z2=（（1/2）.*W1）+1/2*（beta*d1+pi/N）;
F1=sin（pi/2/N）.*sin（z1）./（sin（z2））;
K1=abs（F1）;
d2=lambda*0.5;%有栅瓣现象的强方向性端射阵，即间隔d>1/2波长*（1-1/N）
W2=beta.*d2.*cos（t）;%定义kdcos（方向角）
z3=（（N/2）.*W2）+N/2*（beta*d2+pi/N）;
z4=（（1/2）.*W2）+1/2*（beta*d2+pi/N）;
F2=sin（pi/2/N）.*sin（z3）./（sin（z4））;
K2=abs（F2）;
figure（3）
subplot（121）;polar（t,K1）;title（'强方向性端射阵f=30GHz,N=15,d=1/4波长'）;
subplot（122）;polar（t,K2）;title（'强方向性端射阵（有栅瓣）f=30GHz,N=15,d=0.5倍波长'）;
升级会员 