MATLAB仿真天线阵代码17页Word下载.docx
《MATLAB仿真天线阵代码17页Word下载.docx》由会员分享，可在线阅读，更多相关《MATLAB仿真天线阵代码17页Word下载.docx（12页珍藏版）》请在冰豆网上搜索。
2*pi;
y1=sin（（N1.*W./2））./（N1.*（sin（W./2）））;
%归一化阵因子
y1=abs（y1）;
r1=max（y1）;
y2=sin（（N2.*W./2））./（N2.*（sin（W./2）））;
y2=abs（y2）;
r2=max（y2）;
y3=sin（（N3.*W./2））./（N3.*（sin（W./2）））;
y3=abs（y3）;
r3=max（y3）;
%归一化阵因子绘图程序,
figure
（1）
subplot（311）;
plot（W,y1）;
gridon;
%绘出N=4等幅等矩阵列的归一化阵因子
xlabel（'
f=3GHz,N=4,d=1/2波长,a=π/2'
）
subplot（312）;
plot（W,y2）;
%绘出N=8等幅等矩阵列的归一化阵因子
f=3GHz,N=8,d=1/2波长,a=π/2'
subplot（313）;
plot（W,y3）;
%绘出N=12等幅等矩阵列的归一化阵因子
f=3GHz,N=12,d=1/2波长,a=π/2'
%只有参数N改变的天线方向图
t=0:
0.01:
W=a+（beta.*d.*cos（t））;
z1=（N1/2）.*（W）;
z2=（1/2）.*（W）;
W1=sin（z1）./（N1.*sin（z2））;
%非归一化的阵因子K1
K1=abs（W1）;
z3=（N2/2）.*（W）;
z4=（1/2）.*（W）;
W2=sin（z3）./（N2.*sin（z4））;
%非归一化的阵因子K2
K2=abs（W2）;
z5=（N3/2）.*（W）;
z6=（1/2）.*（W）;
W3=sin（z5）./（N3.*sin（z6））;
%非归一化的阵因子K3
K3=abs（W3）;
%--------------------绘图函数
figure
（2）
subplot（131）;
polar（t,K1）;
xlabel（'
）;
subplot（132）;
polar（t,K2）;
subplot（133）;
polar（t,K3）;
%只有阵列单元方向角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）;
%--------------------------绘图函数
figure（3）
polar（t,K4）;
f=3GHz,N=10,d=1/4波长,a=0'
polar（t,K5）;
f=3GHz,N=10,d=1/4波长,a=π/2'
polar（t,K6）;
f=3GHz,N=10,d=1/4波长,a=π/2+π/10'
%只有阵列单元间隔d改变的天线方向图
N5=20;
d2=lambda/4;
d3=lambda/2;
d4=0.7*lambda;
a4=pi/2;
D=a4+（beta.*d2.*cos（t））;
p1=（N5/2）.*（D）;
p2=（1/2）.*（D）;
D1=sin（p1）./（N5.*sin（p2））;
%非归一化的阵因子K7
K7=abs（D1）;
E=a4+（beta.*d3.*cos（t））;
q1=（N5/2）.*（E）;
q2=（1/2）.*（E）;
E1=sin（q1）./（N5.*sin（q2））;
%非归一化的阵因子K8
K8=abs（E1）;
F=a4+（beta.*d4.*cos（t））;
r_1=（N5/2）.*（F）;
r_2=（1/2）.*（F）;
F1=sin（r_1）./（N5.*sin（r_2））;
%非归一化的阵因子K9
K9=abs（F1）;
%-----------------------绘图函数
figure（4）
polar（t,K7）;
f=3GHz,N=20,d=1/4波长,a=π/2'
polar（t,K8）;
f=3GHz,N=20,d=1/2波长,a=π/2'
polar（t,K9）;
f=3GHz,N=20,d=0.7波长,a=π/2'
%---------------------------------------3D-天线方向图
n_tehta=130;
%--------------------采样视角点的仰角
n_phi=130;
%--------------------采样点的方向角
[tehta,phi]=meshgrid（eps:
pi./（n_tehta-1）:
pi,...%meshgrid函数为矩形区域的设定范围是epf<
tehta<
π0<
phi<
2π
0:
2*pi./（n_phi-1）:
2*pi）;
t3=tehta;
%-------------只有参数N改变的天线方向3D图
M=a+（beta.*d.*cos（t3））;
%----N1=4;
z_1=（N1/2）.*（M）;
z_2=（1/2）.*（M）;
M1=sin（z_1）./（N1.*sin（z_2））;
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）;
z_3=（N2/2）.*（M）;
z_4=（1/2）.*（M）;
M2=sin（z_3）./（N2.*sin（z_4））;
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）;
z_5=（N3/2）.*（M）;
z_6=（1/2）.*（M）;
M3=sin（z_5）./（N3.*sin（z_6））;
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（'
figure（6）
surf（X2,Y2,Z2）;
figure（7）
surf（X3,Y3,Z3）
%--------------------只有阵列单元方向角a改变的天线方向图
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））;
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））;
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））;
K_6=abs（C_1）;
ra
升级会员 